# 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 3.917371571063995*${_u_distance} variable latticeconst_converted equal 3.917371571063995*1 lattice fcc ${latticeconst_converted} lattice fcc 3.91737157106399 Lattice spacing in x,y,z = 3.9173716 3.9173716 3.9173716 region simbox block 0 10 0 10 0 10 units lattice create_box 1 simbox Created orthogonal box = (0.0000000 0.0000000 0.0000000) to (39.173716 39.173716 39.173716) 1 by 1 by 1 MPI processor grid create_atoms 1 box Created 4000 atoms using lattice units in orthogonal box = (0.0000000 0.0000000 0.0000000) to (39.173716 39.173716 39.173716) create_atoms CPU = 0.005 seconds variable mass_converted equal 195.078*${_u_mass} variable mass_converted equal 195.078*1 kim_interactions Pt WARNING: 'kim_' has been renamed to 'kim '. Please update your input. kim interactions Pt #=== BEGIN kim interactions ================================== pair_style kim MEAM_LAMMPS_KimSeolJi_2017_PtCo__MO_545073984441_001 pair_coeff * * Pt #=== END kim interactions ==================================== mass 1 ${mass_converted} mass 1 195.078 # initial volume variable v equal vol # assign formula variable V0 equal ${v} # evaluate initial value variable V0 equal 60115.200756225 variable V0_metal equal ${V0}/(${_u_distance}*${_u_distance}*${_u_distance}) variable V0_metal equal 60115.200756225/(${_u_distance}*${_u_distance}*${_u_distance}) variable V0_metal equal 60115.200756225/(1*${_u_distance}*${_u_distance}) variable V0_metal equal 60115.200756225/(1*1*${_u_distance}) variable V0_metal equal 60115.200756225/(1*1*1) variable V0_metal_times1000 equal ${V0_metal}*1000 variable V0_metal_times1000 equal 60115.200756225*1000 print "Initial system volume: ${V0_metal} Angstroms^3" Initial system volume: 60115.200756225 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_545073984441_001#item-citation CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 13 13 13 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 0 -22907.791 -22907.791 -23080 -23080 333.15 333.15 60115.201 60115.201 3059.7538 3059.7538 1000 -22727.252 -22727.252 -22909.134 -22909.134 351.86297 351.86297 60793.1 60793.1 -3709.838 -3709.838 Loop time of 196.292 on 1 procs for 1000 steps with 4000 atoms Performance: 0.440 ns/day, 54.526 hours/ns, 5.094 timesteps/s 100.0% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 195.3 | 195.3 | 195.3 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.15594 | 0.15594 | 0.15594 | 0.0 | 0.08 Output | 0.00022094 | 0.00022094 | 0.00022094 | 0.0 | 0.00 Modify | 0.72991 | 0.72991 | 0.72991 | 0.0 | 0.37 Other | | 0.1012 | | | 0.05 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 312000.0 ave 312000 max 312000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 312000 Ave neighs/atom = 78.000000 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 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 13 13 13 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 1000 -22727.252 -22727.252 -22909.134 -22909.134 351.86297 351.86297 60793.1 60793.1 -3709.838 -3709.838 2000 -22736.054 -22736.054 -22913.154 -22913.154 342.61214 342.61214 60708.229 60708.229 -525.31801 -525.31801 Loop time of 222.336 on 1 procs for 1000 steps with 4000 atoms Performance: 0.389 ns/day, 61.760 hours/ns, 4.498 timesteps/s 99.9% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 221.27 | 221.27 | 221.27 | 0.0 | 99.52 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16566 | 0.16566 | 0.16566 | 0.0 | 0.07 Output | 0.00021806 | 0.00021806 | 0.00021806 | 0.0 | 0.00 Modify | 0.79651 | 0.79651 | 0.79651 | 0.0 | 0.36 Other | | 0.1062 | | | 0.05 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 311912.0 ave 311912 max 311912 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311912 Ave neighs/atom = 77.978000 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 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 13 13 13 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 2000 -22736.054 -22736.054 -22913.154 -22913.154 342.61214 342.61214 60708.229 60708.229 -525.31801 -525.31801 3000 -22738.201 -22738.201 -22913.914 -22913.914 339.92749 339.92749 60709.687 60709.687 -759.35402 -759.35402 Loop time of 218.888 on 1 procs for 1000 steps with 4000 atoms Performance: 0.395 ns/day, 60.802 hours/ns, 4.569 timesteps/s 99.8% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 217.82 | 217.82 | 217.82 | 0.0 | 99.51 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16447 | 0.16447 | 0.16447 | 0.0 | 0.08 Output | 0.00021392 | 0.00021392 | 0.00021392 | 0.0 | 0.00 Modify | 0.79407 | 0.79407 | 0.79407 | 0.0 | 0.36 Other | | 0.1053 | | | 0.05 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 311962.0 ave 311962 max 311962 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311962 Ave neighs/atom = 77.990500 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 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 13 13 13 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 3000 -22738.201 -22738.201 -22913.914 -22913.914 339.92749 339.92749 60709.687 60709.687 -759.35402 -759.35402 4000 -22732.743 -22732.743 -22909.708 -22909.708 342.35041 342.35041 60694.946 60694.946 835.29619 835.29619 Loop time of 210.081 on 1 procs for 1000 steps with 4000 atoms Performance: 0.411 ns/day, 58.356 hours/ns, 4.760 timesteps/s 99.9% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 209.08 | 209.08 | 209.08 | 0.0 | 99.52 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.15929 | 0.15929 | 0.15929 | 0.0 | 0.08 Output | 0.00021399 | 0.00021399 | 0.00021399 | 0.0 | 0.00 Modify | 0.74221 | 0.74221 | 0.74221 | 0.0 | 0.35 Other | | 0.1012 | | | 0.05 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 311934.0 ave 311934 max 311934 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311934 Ave neighs/atom = 77.983500 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 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 13 13 13 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 4000 -22732.743 -22732.743 -22909.708 -22909.708 342.35041 342.35041 60694.946 60694.946 835.29619 835.29619 5000 -22738.273 -22738.273 -22906.813 -22906.813 326.05207 326.05207 60703.787 60703.787 378.62302 378.62302 Loop time of 221.854 on 1 procs for 1000 steps with 4000 atoms Performance: 0.389 ns/day, 61.626 hours/ns, 4.507 timesteps/s 99.8% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 220.78 | 220.78 | 220.78 | 0.0 | 99.52 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16561 | 0.16561 | 0.16561 | 0.0 | 0.07 Output | 0.00030023 | 0.00030023 | 0.00030023 | 0.0 | 0.00 Modify | 0.80305 | 0.80305 | 0.80305 | 0.0 | 0.36 Other | | 0.1066 | | | 0.05 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 311956.0 ave 311956 max 311956 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311956 Ave neighs/atom = 77.989000 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 = 335.64961627685, Press = 323.638384166388 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 13 13 13 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 5000 -22738.273 -22738.273 -22906.813 -22906.813 326.05207 326.05207 60703.787 60703.787 378.62302 378.62302 6000 -22737.637 -22737.637 -22911.02 -22911.02 335.41966 335.41966 60683.166 60683.166 957.09315 957.09315 Loop time of 216.61 on 1 procs for 1000 steps with 4000 atoms Performance: 0.399 ns/day, 60.169 hours/ns, 4.617 timesteps/s 99.9% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 215.54 | 215.54 | 215.54 | 0.0 | 99.51 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16394 | 0.16394 | 0.16394 | 0.0 | 0.08 Output | 0.00021815 | 0.00021815 | 0.00021815 | 0.0 | 0.00 Modify | 0.79688 | 0.79688 | 0.79688 | 0.0 | 0.37 Other | | 0.1052 | | | 0.05 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 311908.0 ave 311908 max 311908 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311908 Ave neighs/atom = 77.977000 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.595166242606, Press = -13.4420850242124 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 13 13 13 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 6000 -22737.637 -22737.637 -22911.02 -22911.02 335.41966 335.41966 60683.166 60683.166 957.09315 957.09315 7000 -22733.996 -22733.996 -22907.92 -22907.92 336.46831 336.46831 60721.483 60721.483 -450.13297 -450.13297 Loop time of 215.236 on 1 procs for 1000 steps with 4000 atoms Performance: 0.401 ns/day, 59.788 hours/ns, 4.646 timesteps/s 99.9% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 214.19 | 214.19 | 214.19 | 0.0 | 99.51 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16034 | 0.16034 | 0.16034 | 0.0 | 0.07 Output | 0.00017715 | 0.00017715 | 0.00017715 | 0.0 | 0.00 Modify | 0.78584 | 0.78584 | 0.78584 | 0.0 | 0.37 Other | | 0.1039 | | | 0.05 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 311922.0 ave 311922 max 311922 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311922 Ave neighs/atom = 77.980500 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.972777973307, Press = 0.691538958367015 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 13 13 13 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 7000 -22733.996 -22733.996 -22907.92 -22907.92 336.46831 336.46831 60721.483 60721.483 -450.13297 -450.13297 8000 -22738.779 -22738.779 -22910.184 -22910.184 331.59347 331.59347 60701.369 60701.369 246.57876 246.57876 Loop time of 216.573 on 1 procs for 1000 steps with 4000 atoms Performance: 0.399 ns/day, 60.159 hours/ns, 4.617 timesteps/s 99.9% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 215.5 | 215.5 | 215.5 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16307 | 0.16307 | 0.16307 | 0.0 | 0.08 Output | 0.00023552 | 0.00023552 | 0.00023552 | 0.0 | 0.00 Modify | 0.80417 | 0.80417 | 0.80417 | 0.0 | 0.37 Other | | 0.1048 | | | 0.05 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 311888.0 ave 311888 max 311888 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311888 Ave neighs/atom = 77.972000 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.230956108026, Press = -15.1089946157047 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 13 13 13 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 8000 -22738.779 -22738.779 -22910.184 -22910.184 331.59347 331.59347 60701.369 60701.369 246.57876 246.57876 9000 -22732.156 -22732.156 -22908.229 -22908.229 340.6235 340.6235 60665.572 60665.572 2131.3408 2131.3408 Loop time of 223.721 on 1 procs for 1000 steps with 4000 atoms Performance: 0.386 ns/day, 62.145 hours/ns, 4.470 timesteps/s 99.8% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 222.61 | 222.61 | 222.61 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16652 | 0.16652 | 0.16652 | 0.0 | 0.07 Output | 0.00017446 | 0.00017446 | 0.00017446 | 0.0 | 0.00 Modify | 0.836 | 0.836 | 0.836 | 0.0 | 0.37 Other | | 0.1083 | | | 0.05 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 311938.0 ave 311938 max 311938 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311938 Ave neighs/atom = 77.984500 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.421333294312, Press = 5.1929138269364 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 13 13 13 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 9000 -22732.156 -22732.156 -22908.229 -22908.229 340.6235 340.6235 60665.572 60665.572 2131.3408 2131.3408 10000 -22738.882 -22738.882 -22909.913 -22909.913 330.87121 330.87121 60748.361 60748.361 -1982.6347 -1982.6347 Loop time of 219.243 on 1 procs for 1000 steps with 4000 atoms Performance: 0.394 ns/day, 60.901 hours/ns, 4.561 timesteps/s 99.9% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 218.15 | 218.15 | 218.15 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16342 | 0.16342 | 0.16342 | 0.0 | 0.07 Output | 0.00017968 | 0.00017968 | 0.00017968 | 0.0 | 0.00 Modify | 0.81818 | 0.81818 | 0.81818 | 0.0 | 0.37 Other | | 0.1064 | | | 0.05 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 311914.0 ave 311914 max 311914 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311914 Ave neighs/atom = 77.978500 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.762933673706, Press = -4.45655819274515 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 13 13 13 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 10000 -22738.882 -22738.882 -22909.913 -22909.913 330.87121 330.87121 60748.361 60748.361 -1982.6347 -1982.6347 11000 -22732.409 -22732.409 -22906.031 -22906.031 335.88498 335.88498 60700.942 60700.942 836.7391 836.7391 Loop time of 189.875 on 1 procs for 1000 steps with 4000 atoms Performance: 0.455 ns/day, 52.743 hours/ns, 5.267 timesteps/s 100.1% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 188.97 | 188.97 | 188.97 | 0.0 | 99.52 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.14558 | 0.14558 | 0.14558 | 0.0 | 0.08 Output | 0.00018002 | 0.00018002 | 0.00018002 | 0.0 | 0.00 Modify | 0.66543 | 0.66543 | 0.66543 | 0.0 | 0.35 Other | | 0.09278 | | | 0.05 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 311940.0 ave 311940 max 311940 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311940 Ave neighs/atom = 77.985000 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.640701121525, Press = -11.2866537630502 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 13 13 13 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 11000 -22732.409 -22732.409 -22906.031 -22906.031 335.88498 335.88498 60700.942 60700.942 836.7391 836.7391 12000 -22738.943 -22738.943 -22911.471 -22911.471 333.76746 333.76746 60711.355 60711.355 -379.70822 -379.70822 Loop time of 208.276 on 1 procs for 1000 steps with 4000 atoms Performance: 0.415 ns/day, 57.854 hours/ns, 4.801 timesteps/s 99.9% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 207.26 | 207.26 | 207.26 | 0.0 | 99.51 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.15655 | 0.15655 | 0.15655 | 0.0 | 0.08 Output | 0.00018111 | 0.00018111 | 0.00018111 | 0.0 | 0.00 Modify | 0.75411 | 0.75411 | 0.75411 | 0.0 | 0.36 Other | | 0.1016 | | | 0.05 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 311924.0 ave 311924 max 311924 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311924 Ave neighs/atom = 77.981000 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.658765420142, Press = 4.76174885331589 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 13 13 13 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 12000 -22738.943 -22738.943 -22911.471 -22911.471 333.76746 333.76746 60711.355 60711.355 -379.70822 -379.70822 13000 -22732.353 -22732.353 -22906.996 -22906.996 337.85861 337.85861 60758.575 60758.575 -1946.8057 -1946.8057 Loop time of 217.283 on 1 procs for 1000 steps with 4000 atoms Performance: 0.398 ns/day, 60.356 hours/ns, 4.602 timesteps/s 99.9% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 216.2 | 216.2 | 216.2 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16498 | 0.16498 | 0.16498 | 0.0 | 0.08 Output | 0.00017354 | 0.00017354 | 0.00017354 | 0.0 | 0.00 Modify | 0.81366 | 0.81366 | 0.81366 | 0.0 | 0.37 Other | | 0.1064 | | | 0.05 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 311926.0 ave 311926 max 311926 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311926 Ave neighs/atom = 77.981500 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.781780993164, Press = -1.52152038791671 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 13 13 13 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 13000 -22732.353 -22732.353 -22906.996 -22906.996 337.85861 337.85861 60758.575 60758.575 -1946.8057 -1946.8057 14000 -22737.154 -22737.154 -22909.555 -22909.555 333.5217 333.5217 60678.346 60678.346 1272.3329 1272.3329 Loop time of 217.503 on 1 procs for 1000 steps with 4000 atoms Performance: 0.397 ns/day, 60.417 hours/ns, 4.598 timesteps/s 100.0% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 216.43 | 216.43 | 216.43 | 0.0 | 99.51 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16138 | 0.16138 | 0.16138 | 0.0 | 0.07 Output | 0.00017808 | 0.00017808 | 0.00017808 | 0.0 | 0.00 Modify | 0.80633 | 0.80633 | 0.80633 | 0.0 | 0.37 Other | | 0.1064 | | | 0.05 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 311908.0 ave 311908 max 311908 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311908 Ave neighs/atom = 77.977000 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.948269403119, Press = -2.60423226193707 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 13 13 13 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 14000 -22737.154 -22737.154 -22909.555 -22909.555 333.5217 333.5217 60678.346 60678.346 1272.3329 1272.3329 15000 -22739.506 -22739.506 -22909.448 -22909.448 328.76452 328.76452 60728.555 60728.555 -1006.425 -1006.425 Loop time of 220.663 on 1 procs for 1000 steps with 4000 atoms Performance: 0.392 ns/day, 61.295 hours/ns, 4.532 timesteps/s 100.0% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 219.57 | 219.57 | 219.57 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16471 | 0.16471 | 0.16471 | 0.0 | 0.07 Output | 0.00053315 | 0.00053315 | 0.00053315 | 0.0 | 0.00 Modify | 0.82608 | 0.82608 | 0.82608 | 0.0 | 0.37 Other | | 0.1062 | | | 0.05 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 311944.0 ave 311944 max 311944 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311944 Ave neighs/atom = 77.986000 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.887449241976, Press = 2.5915244312164 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 13 13 13 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 15000 -22739.506 -22739.506 -22909.448 -22909.448 328.76452 328.76452 60728.555 60728.555 -1006.425 -1006.425 16000 -22741.67 -22741.67 -22910.835 -22910.835 327.26212 327.26212 60723.259 60723.259 -991.71066 -991.71066 Loop time of 213.634 on 1 procs for 1000 steps with 4000 atoms Performance: 0.404 ns/day, 59.343 hours/ns, 4.681 timesteps/s 100.0% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 212.57 | 212.57 | 212.57 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16594 | 0.16594 | 0.16594 | 0.0 | 0.08 Output | 0.00018235 | 0.00018235 | 0.00018235 | 0.0 | 0.00 Modify | 0.79641 | 0.79641 | 0.79641 | 0.0 | 0.37 Other | | 0.1029 | | | 0.05 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 311954.0 ave 311954 max 311954 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311954 Ave neighs/atom = 77.988500 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.81868832886, Press = -4.70767176502436 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 13 13 13 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 16000 -22741.67 -22741.67 -22910.835 -22910.835 327.26212 327.26212 60723.259 60723.259 -991.71066 -991.71066 17000 -22744.101 -22744.101 -22911.572 -22911.572 323.98412 323.98412 60674.397 60674.397 1228.1113 1228.1113 Loop time of 189.133 on 1 procs for 1000 steps with 4000 atoms Performance: 0.457 ns/day, 52.537 hours/ns, 5.287 timesteps/s 100.1% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 188.23 | 188.23 | 188.23 | 0.0 | 99.52 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.14641 | 0.14641 | 0.14641 | 0.0 | 0.08 Output | 0.00017702 | 0.00017702 | 0.00017702 | 0.0 | 0.00 Modify | 0.66396 | 0.66396 | 0.66396 | 0.0 | 0.35 Other | | 0.0925 | | | 0.05 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 311944.0 ave 311944 max 311944 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311944 Ave neighs/atom = 77.986000 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.573957341161, Press = -0.270162696091106 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 13 13 13 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 17000 -22744.101 -22744.101 -22911.572 -22911.572 323.98412 323.98412 60674.397 60674.397 1228.1113 1228.1113 18000 -22735.221 -22735.221 -22907.882 -22907.882 334.02498 334.02498 60693.449 60693.449 946.72585 946.72585 Loop time of 213.201 on 1 procs for 1000 steps with 4000 atoms Performance: 0.405 ns/day, 59.222 hours/ns, 4.690 timesteps/s 99.9% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 212.15 | 212.15 | 212.15 | 0.0 | 99.51 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16106 | 0.16106 | 0.16106 | 0.0 | 0.08 Output | 0.00017531 | 0.00017531 | 0.00017531 | 0.0 | 0.00 Modify | 0.78318 | 0.78318 | 0.78318 | 0.0 | 0.37 Other | | 0.1033 | | | 0.05 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 311938.0 ave 311938 max 311938 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311938 Ave neighs/atom = 77.984500 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.47926609492, Press = 0.759055773994423 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 13 13 13 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 18000 -22735.221 -22735.221 -22907.882 -22907.882 334.02498 334.02498 60693.449 60693.449 946.72585 946.72585 19000 -22742.187 -22742.187 -22913.129 -22913.129 330.69819 330.69819 60751.704 60751.704 -2490.0247 -2490.0247 Loop time of 211.924 on 1 procs for 1000 steps with 4000 atoms Performance: 0.408 ns/day, 58.868 hours/ns, 4.719 timesteps/s 99.9% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 210.89 | 210.89 | 210.89 | 0.0 | 99.51 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.15912 | 0.15912 | 0.15912 | 0.0 | 0.08 Output | 0.00017729 | 0.00017729 | 0.00017729 | 0.0 | 0.00 Modify | 0.77448 | 0.77448 | 0.77448 | 0.0 | 0.37 Other | | 0.1027 | | | 0.05 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 311930.0 ave 311930 max 311930 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311930 Ave neighs/atom = 77.982500 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.395738246372, Press = -3.1063880121466 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 13 13 13 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 19000 -22742.187 -22742.187 -22913.129 -22913.129 330.69819 330.69819 60751.704 60751.704 -2490.0247 -2490.0247 20000 -22739.011 -22739.011 -22911.915 -22911.915 334.49375 334.49375 60625.227 60625.227 3644.542 3644.542 Loop time of 206.646 on 1 procs for 1000 steps with 4000 atoms Performance: 0.418 ns/day, 57.402 hours/ns, 4.839 timesteps/s 99.9% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 205.64 | 205.64 | 205.64 | 0.0 | 99.51 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.15532 | 0.15532 | 0.15532 | 0.0 | 0.08 Output | 0.00018528 | 0.00018528 | 0.00018528 | 0.0 | 0.00 Modify | 0.74762 | 0.74762 | 0.74762 | 0.0 | 0.36 Other | | 0.101 | | | 0.05 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 311932.0 ave 311932 max 311932 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311932 Ave neighs/atom = 77.983000 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.38868539348, Press = -3.49654124498807 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 13 13 13 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 20000 -22739.011 -22739.011 -22911.915 -22911.915 334.49375 334.49375 60625.227 60625.227 3644.542 3644.542 21000 -22741.815 -22741.815 -22911.753 -22911.753 328.75571 328.75571 60635.357 60635.357 3107.8522 3107.8522 Loop time of 210.852 on 1 procs for 1000 steps with 4000 atoms Performance: 0.410 ns/day, 58.570 hours/ns, 4.743 timesteps/s 99.9% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 209.82 | 209.82 | 209.82 | 0.0 | 99.51 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.15911 | 0.15911 | 0.15911 | 0.0 | 0.08 Output | 0.00027942 | 0.00027942 | 0.00027942 | 0.0 | 0.00 Modify | 0.76927 | 0.76927 | 0.76927 | 0.0 | 0.36 Other | | 0.1027 | | | 0.05 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 311956.0 ave 311956 max 311956 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311956 Ave neighs/atom = 77.989000 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.398922522323, Press = 2.78340288864404 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 13 13 13 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 21000 -22741.815 -22741.815 -22911.753 -22911.753 328.75571 328.75571 60635.357 60635.357 3107.8522 3107.8522 22000 -22735.801 -22735.801 -22909.298 -22909.298 335.64153 335.64153 60714.166 60714.166 -145.386 -145.386 Loop time of 218.557 on 1 procs for 1000 steps with 4000 atoms Performance: 0.395 ns/day, 60.710 hours/ns, 4.575 timesteps/s 100.0% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 217.47 | 217.47 | 217.47 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16504 | 0.16504 | 0.16504 | 0.0 | 0.08 Output | 0.00017811 | 0.00017811 | 0.00017811 | 0.0 | 0.00 Modify | 0.81606 | 0.81606 | 0.81606 | 0.0 | 0.37 Other | | 0.1061 | | | 0.05 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 311942.0 ave 311942 max 311942 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311942 Ave neighs/atom = 77.985500 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.232382578581, Press = 0.917162499230255 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 13 13 13 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 22000 -22735.801 -22735.801 -22909.298 -22909.298 335.64153 335.64153 60714.166 60714.166 -145.386 -145.386 23000 -22739.781 -22739.781 -22911.6 -22911.6 332.39473 332.39473 60705.305 60705.305 -181.25003 -181.25003 Loop time of 213.814 on 1 procs for 1000 steps with 4000 atoms Performance: 0.404 ns/day, 59.393 hours/ns, 4.677 timesteps/s 100.0% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 212.76 | 212.76 | 212.76 | 0.0 | 99.51 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.15963 | 0.15963 | 0.15963 | 0.0 | 0.07 Output | 0.0001837 | 0.0001837 | 0.0001837 | 0.0 | 0.00 Modify | 0.78893 | 0.78893 | 0.78893 | 0.0 | 0.37 Other | | 0.1036 | | | 0.05 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 311964.0 ave 311964 max 311964 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311964 Ave neighs/atom = 77.991000 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.20669505598, Press = -0.740674818560954 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 13 13 13 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 23000 -22739.781 -22739.781 -22911.6 -22911.6 332.39473 332.39473 60705.305 60705.305 -181.25003 -181.25003 24000 -22735.783 -22735.783 -22908.524 -22908.524 334.17907 334.17907 60736.736 60736.736 -1256.3348 -1256.3348 Loop time of 218.182 on 1 procs for 1000 steps with 4000 atoms Performance: 0.396 ns/day, 60.606 hours/ns, 4.583 timesteps/s 99.9% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 217.1 | 217.1 | 217.1 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16702 | 0.16702 | 0.16702 | 0.0 | 0.08 Output | 0.00017931 | 0.00017931 | 0.00017931 | 0.0 | 0.00 Modify | 0.81136 | 0.81136 | 0.81136 | 0.0 | 0.37 Other | | 0.1048 | | | 0.05 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 311944.0 ave 311944 max 311944 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311944 Ave neighs/atom = 77.986000 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.232709360541, Press = -2.38446501235527 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 13 13 13 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 24000 -22735.783 -22735.783 -22908.524 -22908.524 334.17907 334.17907 60736.736 60736.736 -1256.3348 -1256.3348 25000 -22741.28 -22741.28 -22910.521 -22910.521 327.40821 327.40821 60658.172 60658.172 2013.0745 2013.0745 Loop time of 197.548 on 1 procs for 1000 steps with 4000 atoms Performance: 0.437 ns/day, 54.874 hours/ns, 5.062 timesteps/s 100.0% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 196.59 | 196.59 | 196.59 | 0.0 | 99.51 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.15371 | 0.15371 | 0.15371 | 0.0 | 0.08 Output | 0.00021952 | 0.00021952 | 0.00021952 | 0.0 | 0.00 Modify | 0.71132 | 0.71132 | 0.71132 | 0.0 | 0.36 Other | | 0.09747 | | | 0.05 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 311898.0 ave 311898 max 311898 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311898 Ave neighs/atom = 77.974500 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.205068113843, Press = -2.2410762254973 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 13 13 13 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 25000 -22741.28 -22741.28 -22910.521 -22910.521 327.40821 327.40821 60658.172 60658.172 2013.0745 2013.0745 26000 -22734.952 -22734.952 -22908.083 -22908.083 334.93292 334.93292 60704.179 60704.179 283.26846 283.26846 Loop time of 218.013 on 1 procs for 1000 steps with 4000 atoms Performance: 0.396 ns/day, 60.559 hours/ns, 4.587 timesteps/s 100.0% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 216.93 | 216.93 | 216.93 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16564 | 0.16564 | 0.16564 | 0.0 | 0.08 Output | 0.00017747 | 0.00017747 | 0.00017747 | 0.0 | 0.00 Modify | 0.8151 | 0.8151 | 0.8151 | 0.0 | 0.37 Other | | 0.1057 | | | 0.05 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 311940.0 ave 311940 max 311940 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311940 Ave neighs/atom = 77.985000 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.294866783582, Press = 1.45842573903793 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 13 13 13 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 26000 -22734.952 -22734.952 -22908.083 -22908.083 334.93292 334.93292 60704.179 60704.179 283.26846 283.26846 27000 -22737.728 -22737.728 -22907.552 -22907.552 328.53521 328.53521 60752.23 60752.23 -1896.0732 -1896.0732 Loop time of 218.095 on 1 procs for 1000 steps with 4000 atoms Performance: 0.396 ns/day, 60.582 hours/ns, 4.585 timesteps/s 100.0% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 217 | 217 | 217 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16789 | 0.16789 | 0.16789 | 0.0 | 0.08 Output | 0.000178 | 0.000178 | 0.000178 | 0.0 | 0.00 Modify | 0.82176 | 0.82176 | 0.82176 | 0.0 | 0.38 Other | | 0.1052 | | | 0.05 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 311924.0 ave 311924 max 311924 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311924 Ave neighs/atom = 77.981000 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 60709.1742579911 print "LAMMPS calculation completed" LAMMPS calculation completed quit 0