# periodic boundary conditions along all three dimensions boundary p p p # Set neighbor skin variable neigh_skin equal 2.0*${_u_distance} variable neigh_skin equal 2.0*1 neighbor ${neigh_skin} bin neighbor 2 bin # create a supercell with cubic lattice (fcc, bcc, sc, or diamond) # using 10*10*10 conventional (orthogonal) unit cells variable latticeconst_converted equal 5.430999964475632*${_u_distance} variable latticeconst_converted equal 5.430999964475632*1 lattice diamond ${latticeconst_converted} lattice diamond 5.43099996447563 Lattice spacing in x,y,z = 5.431 5.431 5.431 region simbox block 0 10 0 10 0 10 units lattice create_box 1 simbox Created orthogonal box = (0 0 0) to (54.31 54.31 54.31) 1 by 1 by 1 MPI processor grid create_atoms 1 box Created 8000 atoms using lattice units in orthogonal box = (0 0 0) to (54.31 54.31 54.31) create_atoms CPU = 0.001 seconds variable mass_converted equal 28.0855*${_u_mass} variable mass_converted equal 28.0855*1 kim_interactions Si WARNING: 'kim_' has been renamed to 'kim '. Please update your input. kim interactions Si #=== BEGIN kim interactions ================================== pair_style kim MEAM_LAMMPS_JelinekGrohHorstemeyer_2012_AlSiMgCuFe__MO_262519520678_002 pair_coeff * * Si #=== END kim interactions ==================================== mass 1 ${mass_converted} mass 1 28.0855 # initial volume variable v equal vol # assign formula variable V0 equal ${v} # evaluate initial value variable V0 equal 160191.474847545 variable V0_metal equal ${V0}/(${_u_distance}*${_u_distance}*${_u_distance}) variable V0_metal equal 160191.474847545/(${_u_distance}*${_u_distance}*${_u_distance}) variable V0_metal equal 160191.474847545/(1*${_u_distance}*${_u_distance}) variable V0_metal equal 160191.474847545/(1*1*${_u_distance}) variable V0_metal equal 160191.474847545/(1*1*1) variable V0_metal_times1000 equal ${V0_metal}*1000 variable V0_metal_times1000 equal 160191.474847545*1000 print "Initial system volume: ${V0_metal} Angstroms^3" Initial system volume: 160191.474847545 Angstroms^3 # set the time step to 0.001 picoseconds variable timestep_converted equal 0.001*${_u_time} variable timestep_converted equal 0.001*1 timestep ${timestep_converted} timestep 0.001 variable temp_converted equal 333.15*${_u_temperature} variable temp_converted equal 333.15*1 variable Tdamp_converted equal 0.1*${_u_time} variable Tdamp_converted equal 0.1*1 variable press_converted equal 0.0*${_u_pressure} variable press_converted equal 0.0*1 variable Pdamp_converted equal 1*${_u_time} variable Pdamp_converted equal 1*1 # create initial velocities consistent with the chosen temperature velocity all create ${temp_converted} 17 mom yes rot yes velocity all create 333.15 17 mom yes rot yes # set NPT ensemble for all atoms fix ensemble all npt temp ${temp_converted} ${temp_converted} ${Tdamp_converted} iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 333.15 ${temp_converted} ${Tdamp_converted} iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 333.15 333.15 ${Tdamp_converted} iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 333.15 333.15 0.1 iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 333.15 333.15 0.1 iso 0 ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 333.15 333.15 0.1 iso 0 0 ${Pdamp_converted} fix ensemble all npt temp 333.15 333.15 0.1 iso 0 0 1 # compute the time averages of pressure, temperature, and volume, respectively # ignore the first 5000 timesteps variable etotal_metal equal etotal/${_u_energy} variable etotal_metal equal etotal/1 variable pe_metal equal pe/${_u_energy} variable pe_metal equal pe/1 variable T_metal equal temp/${_u_temperature} variable T_metal equal temp/1 variable V_metal equal vol/(${_u_distance}*${_u_distance}*${_u_distance}) variable V_metal equal vol/(1*${_u_distance}*${_u_distance}) variable V_metal equal vol/(1*1*${_u_distance}) variable V_metal equal vol/(1*1*1) variable P_metal equal press/${_u_pressure} variable P_metal equal press/1 fix avgmyTemp all ave/time 5 20 100 v_T_metal ave running start 5000 fix avgmyPress all ave/time 5 20 100 v_P_metal ave running start 5000 fix avgmyVol all ave/time 5 20 100 v_V_metal ave running start 5000 # extract fix quantities into variables so they can be used in if-else logic later. variable T equal f_avgmyTemp variable P equal f_avgmyPress variable V equal f_avgmyVol # set error bounds for temperature and pressure in original metal units (K and bar) variable T_low equal "333.15 - 0.2" variable T_up equal "333.15 + 0.2" variable P_low equal "0.0 - 0.2" variable P_up equal "0.0 + 0.2" # print to logfile every 1000 timesteps thermo_style custom step etotal v_etotal_metal pe v_pe_metal temp v_T_metal vol v_V_metal press v_P_metal thermo 1000 # Run a simulation for at most 2000*1000 timesteps. At each 1000th time step, check # whether the temperature and pressure have converged. If yes, break. label top variable a loop 2000 run 1000 CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE Your simulation uses code contributions which should be cited: - OpenKIM: https://doi.org/10.1007/s11837-011-0102-6 @Article{tadmor:elliott:2011, author = {E. B. Tadmor and R. S. Elliott and J. P. Sethna and R. E. Miller and C. A. Becker}, title = {The potential of atomistic simulations and the {K}nowledgebase of {I}nteratomic {M}odels}, journal = {{JOM}}, year = 2011, volume = 63, number = 17, pages = {17}, doi = {10.1007/s11837-011-0102-6} } - OpenKIM potential: https://openkim.org/cite/MO_262519520678_002#item-citation CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 16 16 16 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 7.894 | 7.894 | 7.894 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 0 -36695.539 -36695.539 -37040 -37040 333.15 333.15 160191.47 160191.47 2296.7926 2296.7926 1000 -36312.946 -36312.946 -36671.928 -36671.928 347.19406 347.19406 162513.64 162513.64 -650.79412 -650.79412 Loop time of 33.633 on 1 procs for 1000 steps with 8000 atoms Performance: 2.569 ns/day, 9.342 hours/ns, 29.733 timesteps/s 99.6% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 33.319 | 33.319 | 33.319 | 0.0 | 99.07 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.037409 | 0.037409 | 0.037409 | 0.0 | 0.11 Output | 0.00017191 | 0.00017191 | 0.00017191 | 0.0 | 0.00 Modify | 0.24891 | 0.24891 | 0.24891 | 0.0 | 0.74 Other | | 0.02769 | | | 0.08 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8601 ave 8601 max 8601 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 560000 ave 560000 max 560000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 560000 Ave neighs/atom = 70 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 0, Press = 0 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 16 16 16 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 8.146 | 8.146 | 8.146 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 1000 -36312.946 -36312.946 -36671.928 -36671.928 347.19406 347.19406 162513.64 162513.64 -650.79412 -650.79412 2000 -36350.283 -36350.283 -36699.787 -36699.787 338.02772 338.02772 162312.55 162312.55 -270.66131 -270.66131 Loop time of 41.6744 on 1 procs for 1000 steps with 8000 atoms Performance: 2.073 ns/day, 11.576 hours/ns, 23.996 timesteps/s 99.7% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 41.374 | 41.374 | 41.374 | 0.0 | 99.28 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.037426 | 0.037426 | 0.037426 | 0.0 | 0.09 Output | 0.00010965 | 0.00010965 | 0.00010965 | 0.0 | 0.00 Modify | 0.24068 | 0.24068 | 0.24068 | 0.0 | 0.58 Other | | 0.02181 | | | 0.05 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8594 ave 8594 max 8594 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 577792 ave 577792 max 577792 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 577792 Ave neighs/atom = 72.224 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 0, Press = 0 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 16 16 16 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 8.146 | 8.146 | 8.146 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 2000 -36350.283 -36350.283 -36699.787 -36699.787 338.02772 338.02772 162312.55 162312.55 -270.66131 -270.66131 3000 -36328.58 -36328.58 -36683.431 -36683.431 343.19894 343.19894 162320.51 162320.51 139.61341 139.61341 Loop time of 41.7581 on 1 procs for 1000 steps with 8000 atoms Performance: 2.069 ns/day, 11.599 hours/ns, 23.947 timesteps/s 99.7% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 41.457 | 41.457 | 41.457 | 0.0 | 99.28 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.037284 | 0.037284 | 0.037284 | 0.0 | 0.09 Output | 9.6552e-05 | 9.6552e-05 | 9.6552e-05 | 0.0 | 0.00 Modify | 0.24154 | 0.24154 | 0.24154 | 0.0 | 0.58 Other | | 0.02188 | | | 0.05 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8601 ave 8601 max 8601 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 577934 ave 577934 max 577934 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 577934 Ave neighs/atom = 72.24175 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 0, Press = 0 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 16 16 16 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 8.146 | 8.146 | 8.146 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 3000 -36328.58 -36328.58 -36683.431 -36683.431 343.19894 343.19894 162320.51 162320.51 139.61341 139.61341 4000 -36343.224 -36343.224 -36679.476 -36679.476 325.21024 325.21024 162330.85 162330.85 68.347937 68.347937 Loop time of 41.7615 on 1 procs for 1000 steps with 8000 atoms Performance: 2.069 ns/day, 11.600 hours/ns, 23.945 timesteps/s 99.7% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 41.461 | 41.461 | 41.461 | 0.0 | 99.28 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.037197 | 0.037197 | 0.037197 | 0.0 | 0.09 Output | 9.7424e-05 | 9.7424e-05 | 9.7424e-05 | 0.0 | 0.00 Modify | 0.2413 | 0.2413 | 0.2413 | 0.0 | 0.58 Other | | 0.02205 | | | 0.05 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8599 ave 8599 max 8599 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 578680 ave 578680 max 578680 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 578680 Ave neighs/atom = 72.335 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 0, Press = 0 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 16 16 16 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 8.146 | 8.146 | 8.146 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 4000 -36343.224 -36343.224 -36679.476 -36679.476 325.21024 325.21024 162330.85 162330.85 68.347937 68.347937 5000 -36335.929 -36335.929 -36680.507 -36680.507 333.26287 333.26287 162329.19 162329.19 48.342612 48.342612 Loop time of 41.7632 on 1 procs for 1000 steps with 8000 atoms Performance: 2.069 ns/day, 11.601 hours/ns, 23.945 timesteps/s 99.7% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 41.462 | 41.462 | 41.462 | 0.0 | 99.28 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.036805 | 0.036805 | 0.036805 | 0.0 | 0.09 Output | 6.1737e-05 | 6.1737e-05 | 6.1737e-05 | 0.0 | 0.00 Modify | 0.24224 | 0.24224 | 0.24224 | 0.0 | 0.58 Other | | 0.02189 | | | 0.05 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8594 ave 8594 max 8594 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 579166 ave 579166 max 579166 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 579166 Ave neighs/atom = 72.39575 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 336.594847766683, Press = 55.9708269811527 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 16 16 16 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 8.146 | 8.146 | 8.146 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 5000 -36335.929 -36335.929 -36680.507 -36680.507 333.26287 333.26287 162329.19 162329.19 48.342612 48.342612 6000 -36338.08 -36338.08 -36690.319 -36690.319 340.6727 340.6727 162284.63 162284.63 129.1502 129.1502 Loop time of 41.7299 on 1 procs for 1000 steps with 8000 atoms Performance: 2.070 ns/day, 11.592 hours/ns, 23.964 timesteps/s 99.7% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 41.422 | 41.422 | 41.422 | 0.0 | 99.26 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.036881 | 0.036881 | 0.036881 | 0.0 | 0.09 Output | 7.1064e-05 | 7.1064e-05 | 7.1064e-05 | 0.0 | 0.00 Modify | 0.24923 | 0.24923 | 0.24923 | 0.0 | 0.60 Other | | 0.02211 | | | 0.05 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8601 ave 8601 max 8601 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 578990 ave 578990 max 578990 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 578990 Ave neighs/atom = 72.37375 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 332.784561104135, Press = -8.09869586551353 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 16 16 16 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 8.145 | 8.145 | 8.145 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 6000 -36338.08 -36338.08 -36690.319 -36690.319 340.6727 340.6727 162284.63 162284.63 129.1502 129.1502 7000 -36340.036 -36340.036 -36688.073 -36688.073 336.60814 336.60814 162535.79 162535.79 -1265.714 -1265.714 Loop time of 41.7168 on 1 procs for 1000 steps with 8000 atoms Performance: 2.071 ns/day, 11.588 hours/ns, 23.971 timesteps/s 99.7% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 41.409 | 41.409 | 41.409 | 0.0 | 99.26 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.037148 | 0.037148 | 0.037148 | 0.0 | 0.09 Output | 6.902e-05 | 6.902e-05 | 6.902e-05 | 0.0 | 0.00 Modify | 0.24847 | 0.24847 | 0.24847 | 0.0 | 0.60 Other | | 0.02206 | | | 0.05 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8592 ave 8592 max 8592 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 579078 ave 579078 max 579078 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 579078 Ave neighs/atom = 72.38475 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 333.232859019606, Press = 1.44453051512047 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 16 16 16 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 8.146 | 8.146 | 8.146 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 7000 -36340.036 -36340.036 -36688.073 -36688.073 336.60814 336.60814 162535.79 162535.79 -1265.714 -1265.714 8000 -36334.987 -36334.987 -36678.375 -36678.375 332.11234 332.11234 162245.35 162245.35 607.57464 607.57464 Loop time of 41.709 on 1 procs for 1000 steps with 8000 atoms Performance: 2.071 ns/day, 11.586 hours/ns, 23.976 timesteps/s 99.7% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 41.401 | 41.401 | 41.401 | 0.0 | 99.26 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.037229 | 0.037229 | 0.037229 | 0.0 | 0.09 Output | 5.0686e-05 | 5.0686e-05 | 5.0686e-05 | 0.0 | 0.00 Modify | 0.24875 | 0.24875 | 0.24875 | 0.0 | 0.60 Other | | 0.02197 | | | 0.05 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8598 ave 8598 max 8598 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 576688 ave 576688 max 576688 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 576688 Ave neighs/atom = 72.086 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 333.074675203346, Press = 10.694846049397 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 16 16 16 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 8.146 | 8.146 | 8.146 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 8000 -36334.987 -36334.987 -36678.375 -36678.375 332.11234 332.11234 162245.35 162245.35 607.57464 607.57464 9000 -36343.145 -36343.145 -36680.655 -36680.655 326.42718 326.42718 162269.2 162269.2 363.89078 363.89078 Loop time of 41.7884 on 1 procs for 1000 steps with 8000 atoms Performance: 2.068 ns/day, 11.608 hours/ns, 23.930 timesteps/s 99.7% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 41.481 | 41.481 | 41.481 | 0.0 | 99.27 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.037185 | 0.037185 | 0.037185 | 0.0 | 0.09 Output | 7.8639e-05 | 7.8639e-05 | 7.8639e-05 | 0.0 | 0.00 Modify | 0.24791 | 0.24791 | 0.24791 | 0.0 | 0.59 Other | | 0.02197 | | | 0.05 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8598 ave 8598 max 8598 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 579746 ave 579746 max 579746 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 579746 Ave neighs/atom = 72.46825 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 333.153233775847, Press = -1.08093345767875 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 16 16 16 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 8.146 | 8.146 | 8.146 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 9000 -36343.145 -36343.145 -36680.655 -36680.655 326.42718 326.42718 162269.2 162269.2 363.89078 363.89078 10000 -36333.682 -36333.682 -36674.925 -36674.925 330.03751 330.03751 162403.34 162403.34 -245.23281 -245.23281 Loop time of 41.7976 on 1 procs for 1000 steps with 8000 atoms Performance: 2.067 ns/day, 11.610 hours/ns, 23.925 timesteps/s 99.7% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 41.49 | 41.49 | 41.49 | 0.0 | 99.26 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.037341 | 0.037341 | 0.037341 | 0.0 | 0.09 Output | 7.2036e-05 | 7.2036e-05 | 7.2036e-05 | 0.0 | 0.00 Modify | 0.24811 | 0.24811 | 0.24811 | 0.0 | 0.59 Other | | 0.02201 | | | 0.05 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8600 ave 8600 max 8600 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 579698 ave 579698 max 579698 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 579698 Ave neighs/atom = 72.46225 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 333.143608199849, Press = 0.424165059889542 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 16 16 16 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 8.145 | 8.145 | 8.145 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 10000 -36333.682 -36333.682 -36674.925 -36674.925 330.03751 330.03751 162403.34 162403.34 -245.23281 -245.23281 11000 -36341.501 -36341.501 -36684.395 -36684.395 331.63364 331.63364 162279.92 162279.92 270.18299 270.18299 Loop time of 41.7359 on 1 procs for 1000 steps with 8000 atoms Performance: 2.070 ns/day, 11.593 hours/ns, 23.960 timesteps/s 99.7% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 41.428 | 41.428 | 41.428 | 0.0 | 99.26 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.037209 | 0.037209 | 0.037209 | 0.0 | 0.09 Output | 6.7838e-05 | 6.7838e-05 | 6.7838e-05 | 0.0 | 0.00 Modify | 0.24853 | 0.24853 | 0.24853 | 0.0 | 0.60 Other | | 0.0222 | | | 0.05 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8589 ave 8589 max 8589 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 578896 ave 578896 max 578896 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 578896 Ave neighs/atom = 72.362 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 333.096842727363, Press = 2.04365231786787 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 16 16 16 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 8.146 | 8.146 | 8.146 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 11000 -36341.501 -36341.501 -36684.395 -36684.395 331.63364 331.63364 162279.92 162279.92 270.18299 270.18299 12000 -36338.434 -36338.434 -36680.751 -36680.751 331.07663 331.07663 162296.91 162296.91 247.05328 247.05328 Loop time of 41.7687 on 1 procs for 1000 steps with 8000 atoms Performance: 2.069 ns/day, 11.602 hours/ns, 23.941 timesteps/s 99.7% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 41.461 | 41.461 | 41.461 | 0.0 | 99.26 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.037081 | 0.037081 | 0.037081 | 0.0 | 0.09 Output | 4.8742e-05 | 4.8742e-05 | 4.8742e-05 | 0.0 | 0.00 Modify | 0.24822 | 0.24822 | 0.24822 | 0.0 | 0.59 Other | | 0.02194 | | | 0.05 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8596 ave 8596 max 8596 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 578806 ave 578806 max 578806 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 578806 Ave neighs/atom = 72.35075 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 333.231338326829, Press = -1.93580315756039 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 16 16 16 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 8.146 | 8.146 | 8.146 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 12000 -36338.434 -36338.434 -36680.751 -36680.751 331.07663 331.07663 162296.91 162296.91 247.05328 247.05328 13000 -36343.944 -36343.944 -36684.087 -36684.087 328.9745 328.9745 162514.72 162514.72 -1130.4959 -1130.4959 Loop time of 41.7883 on 1 procs for 1000 steps with 8000 atoms Performance: 2.068 ns/day, 11.608 hours/ns, 23.930 timesteps/s 99.7% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 41.48 | 41.48 | 41.48 | 0.0 | 99.26 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.037277 | 0.037277 | 0.037277 | 0.0 | 0.09 Output | 7.2026e-05 | 7.2026e-05 | 7.2026e-05 | 0.0 | 0.00 Modify | 0.24869 | 0.24869 | 0.24869 | 0.0 | 0.60 Other | | 0.02191 | | | 0.05 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8599 ave 8599 max 8599 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 579028 ave 579028 max 579028 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 579028 Ave neighs/atom = 72.3785 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 333.296827629858, Press = 0.651714899652072 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 16 16 16 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 8.146 | 8.146 | 8.146 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 13000 -36343.944 -36343.944 -36684.087 -36684.087 328.9745 328.9745 162514.72 162514.72 -1130.4959 -1130.4959 14000 -36346.256 -36346.256 -36689.572 -36689.572 332.04237 332.04237 162211.26 162211.26 509.97424 509.97424 Loop time of 41.7396 on 1 procs for 1000 steps with 8000 atoms Performance: 2.070 ns/day, 11.594 hours/ns, 23.958 timesteps/s 99.7% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 41.432 | 41.432 | 41.432 | 0.0 | 99.26 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.036852 | 0.036852 | 0.036852 | 0.0 | 0.09 Output | 0.00013558 | 0.00013558 | 0.00013558 | 0.0 | 0.00 Modify | 0.24871 | 0.24871 | 0.24871 | 0.0 | 0.60 Other | | 0.02194 | | | 0.05 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8595 ave 8595 max 8595 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 577912 ave 577912 max 577912 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 577912 Ave neighs/atom = 72.239 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 333.264656869482, Press = 2.81798560704122 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 16 16 16 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 8.146 | 8.146 | 8.146 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 14000 -36346.256 -36346.256 -36689.572 -36689.572 332.04237 332.04237 162211.26 162211.26 509.97424 509.97424 15000 -36335.629 -36335.629 -36673.142 -36673.142 326.42956 326.42956 162330.74 162330.74 242.5178 242.5178 Loop time of 41.7422 on 1 procs for 1000 steps with 8000 atoms Performance: 2.070 ns/day, 11.595 hours/ns, 23.957 timesteps/s 99.7% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 41.435 | 41.435 | 41.435 | 0.0 | 99.26 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.037137 | 0.037137 | 0.037137 | 0.0 | 0.09 Output | 4.8201e-05 | 4.8201e-05 | 4.8201e-05 | 0.0 | 0.00 Modify | 0.24805 | 0.24805 | 0.24805 | 0.0 | 0.59 Other | | 0.02187 | | | 0.05 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8595 ave 8595 max 8595 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 579242 ave 579242 max 579242 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 579242 Ave neighs/atom = 72.40525 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 333.110106777433, Press = -0.594001281162314 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 16 16 16 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 8.146 | 8.146 | 8.146 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 15000 -36335.629 -36335.629 -36673.142 -36673.142 326.42956 326.42956 162330.74 162330.74 242.5178 242.5178 16000 -36342.437 -36342.437 -36681.092 -36681.092 327.53504 327.53504 162366.45 162366.45 -209.96119 -209.96119 Loop time of 41.7782 on 1 procs for 1000 steps with 8000 atoms Performance: 2.068 ns/day, 11.605 hours/ns, 23.936 timesteps/s 99.7% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 41.471 | 41.471 | 41.471 | 0.0 | 99.26 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.037111 | 0.037111 | 0.037111 | 0.0 | 0.09 Output | 7.6013e-05 | 7.6013e-05 | 7.6013e-05 | 0.0 | 0.00 Modify | 0.24834 | 0.24834 | 0.24834 | 0.0 | 0.59 Other | | 0.02189 | | | 0.05 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8593 ave 8593 max 8593 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 579484 ave 579484 max 579484 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 579484 Ave neighs/atom = 72.4355 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 333.021993294054, Press = 0.622350037061967 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 16 16 16 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 8.146 | 8.146 | 8.146 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 16000 -36342.437 -36342.437 -36681.092 -36681.092 327.53504 327.53504 162366.45 162366.45 -209.96119 -209.96119 17000 -36337.91 -36337.91 -36681.545 -36681.545 332.35136 332.35136 162276.33 162276.33 358.78149 358.78149 Loop time of 41.7407 on 1 procs for 1000 steps with 8000 atoms Performance: 2.070 ns/day, 11.595 hours/ns, 23.957 timesteps/s 99.7% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 41.432 | 41.432 | 41.432 | 0.0 | 99.26 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.037104 | 0.037104 | 0.037104 | 0.0 | 0.09 Output | 6.5844e-05 | 6.5844e-05 | 6.5844e-05 | 0.0 | 0.00 Modify | 0.24902 | 0.24902 | 0.24902 | 0.0 | 0.60 Other | | 0.02207 | | | 0.05 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8598 ave 8598 max 8598 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 579338 ave 579338 max 579338 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 579338 Ave neighs/atom = 72.41725 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 332.941812962821, Press = 1.2413095129551 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 16 16 16 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 8.145 | 8.145 | 8.145 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 17000 -36337.91 -36337.91 -36681.545 -36681.545 332.35136 332.35136 162276.33 162276.33 358.78149 358.78149 18000 -36345.305 -36345.305 -36680.859 -36680.859 324.53539 324.53539 162272.13 162272.13 330.19406 330.19406 Loop time of 41.715 on 1 procs for 1000 steps with 8000 atoms Performance: 2.071 ns/day, 11.588 hours/ns, 23.972 timesteps/s 99.7% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 41.407 | 41.407 | 41.407 | 0.0 | 99.26 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.037184 | 0.037184 | 0.037184 | 0.0 | 0.09 Output | 4.8131e-05 | 4.8131e-05 | 4.8131e-05 | 0.0 | 0.00 Modify | 0.24828 | 0.24828 | 0.24828 | 0.0 | 0.60 Other | | 0.02201 | | | 0.05 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8592 ave 8592 max 8592 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 579928 ave 579928 max 579928 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 579928 Ave neighs/atom = 72.491 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 332.928185824702, Press = -0.9683850278694 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 16 16 16 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 8.145 | 8.145 | 8.145 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 18000 -36345.305 -36345.305 -36680.859 -36680.859 324.53539 324.53539 162272.13 162272.13 330.19406 330.19406 19000 -36328.078 -36328.078 -36682.058 -36682.058 342.35634 342.35634 162532.99 162532.99 -1055.1282 -1055.1282 Loop time of 41.7273 on 1 procs for 1000 steps with 8000 atoms Performance: 2.071 ns/day, 11.591 hours/ns, 23.965 timesteps/s 99.7% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 41.42 | 41.42 | 41.42 | 0.0 | 99.26 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.037346 | 0.037346 | 0.037346 | 0.0 | 0.09 Output | 6.6466e-05 | 6.6466e-05 | 6.6466e-05 | 0.0 | 0.00 Modify | 0.24776 | 0.24776 | 0.24776 | 0.0 | 0.59 Other | | 0.02192 | | | 0.05 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8586 ave 8586 max 8586 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 580052 ave 580052 max 580052 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 580052 Ave neighs/atom = 72.5065 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 332.997354612118, Press = 0.563955946000652 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 16 16 16 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 8.146 | 8.146 | 8.146 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 19000 -36328.078 -36328.078 -36682.058 -36682.058 342.35634 342.35634 162532.99 162532.99 -1055.1282 -1055.1282 20000 -36340.789 -36340.789 -36681.107 -36681.107 329.14248 329.14248 162177.22 162177.22 933.87877 933.87877 Loop time of 41.5964 on 1 procs for 1000 steps with 8000 atoms Performance: 2.077 ns/day, 11.555 hours/ns, 24.041 timesteps/s 99.7% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 41.288 | 41.288 | 41.288 | 0.0 | 99.26 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.03732 | 0.03732 | 0.03732 | 0.0 | 0.09 Output | 6.7097e-05 | 6.7097e-05 | 6.7097e-05 | 0.0 | 0.00 Modify | 0.24854 | 0.24854 | 0.24854 | 0.0 | 0.60 Other | | 0.02199 | | | 0.05 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8599 ave 8599 max 8599 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 577026 ave 577026 max 577026 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 577026 Ave neighs/atom = 72.12825 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 333.066641760677, Press = 2.21067812216785 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 16 16 16 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 8.146 | 8.146 | 8.146 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 20000 -36340.789 -36340.789 -36681.107 -36681.107 329.14248 329.14248 162177.22 162177.22 933.87877 933.87877 21000 -36333.003 -36333.003 -36682.111 -36682.111 337.64387 337.64387 162328.56 162328.56 66.256855 66.256855 Loop time of 41.7262 on 1 procs for 1000 steps with 8000 atoms Performance: 2.071 ns/day, 11.591 hours/ns, 23.966 timesteps/s 99.7% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 41.418 | 41.418 | 41.418 | 0.0 | 99.26 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.037363 | 0.037363 | 0.037363 | 0.0 | 0.09 Output | 6.3119e-05 | 6.3119e-05 | 6.3119e-05 | 0.0 | 0.00 Modify | 0.24841 | 0.24841 | 0.24841 | 0.0 | 0.60 Other | | 0.02211 | | | 0.05 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8596 ave 8596 max 8596 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 580218 ave 580218 max 580218 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 580218 Ave neighs/atom = 72.52725 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 333.074239752253, Press = -0.535839389169172 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 16 16 16 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 8.146 | 8.146 | 8.146 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 21000 -36333.003 -36333.003 -36682.111 -36682.111 337.64387 337.64387 162328.56 162328.56 66.256855 66.256855 22000 -36340.985 -36340.985 -36684.135 -36684.135 331.88114 331.88114 162354.52 162354.52 -177.16771 -177.16771 Loop time of 41.7733 on 1 procs for 1000 steps with 8000 atoms Performance: 2.068 ns/day, 11.604 hours/ns, 23.939 timesteps/s 99.7% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 41.466 | 41.466 | 41.466 | 0.0 | 99.26 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.037314 | 0.037314 | 0.037314 | 0.0 | 0.09 Output | 4.776e-05 | 4.776e-05 | 4.776e-05 | 0.0 | 0.00 Modify | 0.24834 | 0.24834 | 0.24834 | 0.0 | 0.59 Other | | 0.02199 | | | 0.05 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8596 ave 8596 max 8596 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 579242 ave 579242 max 579242 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 579242 Ave neighs/atom = 72.40525 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 333.185180724963, Press = 0.23146758902934 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 16 16 16 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 8.145 | 8.145 | 8.145 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 22000 -36340.985 -36340.985 -36684.135 -36684.135 331.88114 331.88114 162354.52 162354.52 -177.16771 -177.16771 23000 -36333.311 -36333.311 -36673.783 -36673.783 329.29164 329.29164 162305.1 162305.1 419.19946 419.19946 Loop time of 41.7156 on 1 procs for 1000 steps with 8000 atoms Performance: 2.071 ns/day, 11.588 hours/ns, 23.972 timesteps/s 99.7% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 41.408 | 41.408 | 41.408 | 0.0 | 99.26 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.037276 | 0.037276 | 0.037276 | 0.0 | 0.09 Output | 6.6125e-05 | 6.6125e-05 | 6.6125e-05 | 0.0 | 0.00 Modify | 0.24827 | 0.24827 | 0.24827 | 0.0 | 0.60 Other | | 0.02198 | | | 0.05 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8585 ave 8585 max 8585 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 578928 ave 578928 max 578928 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 578928 Ave neighs/atom = 72.366 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 333.176470509516, Press = 0.489042785169561 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 16 16 16 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 8.146 | 8.146 | 8.146 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 23000 -36333.311 -36333.311 -36673.783 -36673.783 329.29164 329.29164 162305.1 162305.1 419.19946 419.19946 24000 -36336.225 -36336.225 -36682.923 -36682.923 335.31363 335.31363 162344.12 162344.12 -34.322445 -34.322445 Loop time of 41.7999 on 1 procs for 1000 steps with 8000 atoms Performance: 2.067 ns/day, 11.611 hours/ns, 23.924 timesteps/s 99.7% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 41.492 | 41.492 | 41.492 | 0.0 | 99.26 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.03747 | 0.03747 | 0.03747 | 0.0 | 0.09 Output | 6.5273e-05 | 6.5273e-05 | 6.5273e-05 | 0.0 | 0.00 Modify | 0.24838 | 0.24838 | 0.24838 | 0.0 | 0.59 Other | | 0.02189 | | | 0.05 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8599 ave 8599 max 8599 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 579370 ave 579370 max 579370 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 579370 Ave neighs/atom = 72.42125 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" jump SELF break # Write final averaged volume to file if temperature and volume have converged; otherwise wirte a # flag to indicate non-convergence. variable myStep equal step if "${myStep} < 2000000" then "print '${V}' file output/vol_T333.15.out" else "print 'not_converged' file output/vol_T333.15.out" print '${V}' file output/vol_T333.15.out 162333.197037151 print "LAMMPS calculation completed" LAMMPS calculation completed quit 0