# 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.917371600866318*${_u_distance} variable latticeconst_converted equal 3.917371600866318*1 lattice fcc ${latticeconst_converted} lattice fcc 3.91737160086632 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.006 seconds variable mass_converted equal 63.546*${_u_mass} variable mass_converted equal 63.546*1 kim_interactions Cu WARNING: 'kim_' has been renamed to 'kim '. Please update your input. kim interactions Cu #=== BEGIN kim interactions ================================== pair_style kim MEAM_LAMMPS_KimSeolJi_2017_PtCu__MO_070797404269_001 pair_coeff * * Cu #=== END kim interactions ==================================== mass 1 ${mass_converted} mass 1 63.546 # initial volume variable v equal vol # assign formula variable V0 equal ${v} # evaluate initial value variable V0 equal 60115.2021282468 variable V0_metal equal ${V0}/(${_u_distance}*${_u_distance}*${_u_distance}) variable V0_metal equal 60115.2021282468/(${_u_distance}*${_u_distance}*${_u_distance}) variable V0_metal equal 60115.2021282468/(1*${_u_distance}*${_u_distance}) variable V0_metal equal 60115.2021282468/(1*1*${_u_distance}) variable V0_metal equal 60115.2021282468/(1*1*1) variable V0_metal_times1000 equal ${V0_metal}*1000 variable V0_metal_times1000 equal 60115.2021282468*1000 print "Initial system volume: ${V0_metal} Angstroms^3" Initial system volume: 60115.2021282468 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 253.15*${_u_temperature} variable temp_converted equal 253.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 253.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 253.15 ${temp_converted} ${Tdamp_converted} iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 253.15 253.15 ${Tdamp_converted} iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 253.15 253.15 0.1 iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 253.15 253.15 0.1 iso 0 ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 253.15 253.15 0.1 iso 0 0 ${Pdamp_converted} fix ensemble all npt temp 253.15 253.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 "253.15 - 0.2" variable T_up equal "253.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_070797404269_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 -22949.144 -22949.144 -23080 -23080 253.15 253.15 60115.202 60115.202 2324.976 2324.976 1000 -22814.25 -22814.25 -22943.803 -22943.803 250.62955 250.62955 60712.385 60712.385 -2014.365 -2014.365 Loop time of 184.638 on 1 procs for 1000 steps with 4000 atoms Performance: 0.468 ns/day, 51.288 hours/ns, 5.416 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 | 183.68 | 183.68 | 183.68 | 0.0 | 99.48 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.15712 | 0.15712 | 0.15712 | 0.0 | 0.09 Output | 0.00022137 | 0.00022137 | 0.00022137 | 0.0 | 0.00 Modify | 0.70491 | 0.70491 | 0.70491 | 0.0 | 0.38 Other | | 0.1002 | | | 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 -22814.25 -22814.25 -22943.803 -22943.803 250.62955 250.62955 60712.385 60712.385 -2014.365 -2014.365 2000 -22820.451 -22820.451 -22945.194 -22945.194 241.32461 241.32461 60617.546 60617.546 374.08488 374.08488 Loop time of 193.043 on 1 procs for 1000 steps with 4000 atoms Performance: 0.448 ns/day, 53.623 hours/ns, 5.180 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 | 192.11 | 192.11 | 192.11 | 0.0 | 99.52 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.15158 | 0.15158 | 0.15158 | 0.0 | 0.08 Output | 0.00021786 | 0.00021786 | 0.00021786 | 0.0 | 0.00 Modify | 0.6844 | 0.6844 | 0.6844 | 0.0 | 0.35 Other | | 0.09903 | | | 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: 311958.0 ave 311958 max 311958 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311958 Ave neighs/atom = 77.989500 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 -22820.451 -22820.451 -22945.194 -22945.194 241.32461 241.32461 60617.546 60617.546 374.08488 374.08488 3000 -22820.062 -22820.062 -22947.321 -22947.321 246.19076 246.19076 60623.167 60623.167 3.5702571 3.5702571 Loop time of 192.899 on 1 procs for 1000 steps with 4000 atoms Performance: 0.448 ns/day, 53.583 hours/ns, 5.184 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 | 191.96 | 191.96 | 191.96 | 0.0 | 99.51 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.15213 | 0.15213 | 0.15213 | 0.0 | 0.08 Output | 0.00021204 | 0.00021204 | 0.00021204 | 0.0 | 0.00 Modify | 0.68643 | 0.68643 | 0.68643 | 0.0 | 0.36 Other | | 0.09891 | | | 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: 311966.0 ave 311966 max 311966 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311966 Ave neighs/atom = 77.991500 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 -22820.062 -22820.062 -22947.321 -22947.321 246.19076 246.19076 60623.167 60623.167 3.5702571 3.5702571 4000 -22817.9 -22817.9 -22956.286 -22956.286 267.71714 267.71714 60625.034 60625.034 -704.04546 -704.04546 Loop time of 194.051 on 1 procs for 1000 steps with 4000 atoms Performance: 0.445 ns/day, 53.903 hours/ns, 5.153 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 | 193.12 | 193.12 | 193.12 | 0.0 | 99.52 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.15185 | 0.15185 | 0.15185 | 0.0 | 0.08 Output | 0.0002131 | 0.0002131 | 0.0002131 | 0.0 | 0.00 Modify | 0.68512 | 0.68512 | 0.68512 | 0.0 | 0.35 Other | | 0.09889 | | | 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 = 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 -22817.9 -22817.9 -22956.286 -22956.286 267.71714 267.71714 60625.034 60625.034 -704.04546 -704.04546 5000 -22821.886 -22821.886 -22948.752 -22948.752 245.42999 245.42999 60595.895 60595.895 699.56055 699.56055 Loop time of 192.062 on 1 procs for 1000 steps with 4000 atoms Performance: 0.450 ns/day, 53.351 hours/ns, 5.207 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 | 191.13 | 191.13 | 191.13 | 0.0 | 99.51 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.14954 | 0.14954 | 0.14954 | 0.0 | 0.08 Output | 0.00018299 | 0.00018299 | 0.00018299 | 0.0 | 0.00 Modify | 0.68684 | 0.68684 | 0.68684 | 0.0 | 0.36 Other | | 0.0986 | | | 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: 311980.0 ave 311980 max 311980 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311980 Ave neighs/atom = 77.995000 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 = 250.147235594238, Press = -3.5759416020137 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 -22821.886 -22821.886 -22948.752 -22948.752 245.42999 245.42999 60595.895 60595.895 699.56055 699.56055 6000 -22816.648 -22816.648 -22944.203 -22944.203 246.76385 246.76385 60639.022 60639.022 -153.74279 -153.74279 Loop time of 193.495 on 1 procs for 1000 steps with 4000 atoms Performance: 0.447 ns/day, 53.749 hours/ns, 5.168 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 | 192.53 | 192.53 | 192.53 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.1517 | 0.1517 | 0.1517 | 0.0 | 0.08 Output | 0.00022242 | 0.00022242 | 0.00022242 | 0.0 | 0.00 Modify | 0.71401 | 0.71401 | 0.71401 | 0.0 | 0.37 Other | | 0.0986 | | | 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: 311966.0 ave 311966 max 311966 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311966 Ave neighs/atom = 77.991500 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 = 252.505087212013, Press = -30.432156288179 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 -22816.648 -22816.648 -22944.203 -22944.203 246.76385 246.76385 60639.022 60639.022 -153.74279 -153.74279 7000 -22823.324 -22823.324 -22947.994 -22947.994 241.18167 241.18167 60630.633 60630.633 -296.45169 -296.45169 Loop time of 191.649 on 1 procs for 1000 steps with 4000 atoms Performance: 0.451 ns/day, 53.236 hours/ns, 5.218 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 | 190.7 | 190.7 | 190.7 | 0.0 | 99.51 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.14953 | 0.14953 | 0.14953 | 0.0 | 0.08 Output | 0.00017272 | 0.00017272 | 0.00017272 | 0.0 | 0.00 Modify | 0.69736 | 0.69736 | 0.69736 | 0.0 | 0.36 Other | | 0.09789 | | | 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: 311952.0 ave 311952 max 311952 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311952 Ave neighs/atom = 77.988000 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 = 252.922874874903, Press = 19.5547880930129 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 -22823.324 -22823.324 -22947.994 -22947.994 241.18167 241.18167 60630.633 60630.633 -296.45169 -296.45169 8000 -22818.198 -22818.198 -22952.003 -22952.003 258.85599 258.85599 60573.288 60573.288 1247.1891 1247.1891 Loop time of 188.661 on 1 procs for 1000 steps with 4000 atoms Performance: 0.458 ns/day, 52.406 hours/ns, 5.301 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 | 187.73 | 187.73 | 187.73 | 0.0 | 99.51 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.14825 | 0.14825 | 0.14825 | 0.0 | 0.08 Output | 0.00024487 | 0.00024487 | 0.00024487 | 0.0 | 0.00 Modify | 0.68802 | 0.68802 | 0.68802 | 0.0 | 0.36 Other | | 0.09649 | | | 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: 311968.0 ave 311968 max 311968 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311968 Ave neighs/atom = 77.992000 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 = 252.737490947677, Press = -4.24545764538524 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 -22818.198 -22818.198 -22952.003 -22952.003 258.85599 258.85599 60573.288 60573.288 1247.1891 1247.1891 9000 -22820.982 -22820.982 -22951.443 -22951.443 252.38627 252.38627 60630.546 60630.546 -536.39291 -536.39291 Loop time of 189.173 on 1 procs for 1000 steps with 4000 atoms Performance: 0.457 ns/day, 52.548 hours/ns, 5.286 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 | 188.23 | 188.23 | 188.23 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.14915 | 0.14915 | 0.14915 | 0.0 | 0.08 Output | 0.00022069 | 0.00022069 | 0.00022069 | 0.0 | 0.00 Modify | 0.69361 | 0.69361 | 0.69361 | 0.0 | 0.37 Other | | 0.09724 | | | 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: 311982.0 ave 311982 max 311982 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311982 Ave neighs/atom = 77.995500 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 = 252.631453177882, Press = -2.62030677421006 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 -22820.982 -22820.982 -22951.443 -22951.443 252.38627 252.38627 60630.546 60630.546 -536.39291 -536.39291 10000 -22816.56 -22816.56 -22946.993 -22946.993 252.33187 252.33187 60618.759 60618.759 256.84433 256.84433 Loop time of 190.784 on 1 procs for 1000 steps with 4000 atoms Performance: 0.453 ns/day, 52.996 hours/ns, 5.242 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 | 189.83 | 189.83 | 189.83 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.15002 | 0.15002 | 0.15002 | 0.0 | 0.08 Output | 0.00018234 | 0.00018234 | 0.00018234 | 0.0 | 0.00 Modify | 0.70144 | 0.70144 | 0.70144 | 0.0 | 0.37 Other | | 0.09849 | | | 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: 311974.0 ave 311974 max 311974 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311974 Ave neighs/atom = 77.993500 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 = 252.72769290962, Press = 3.13652439479713 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 -22816.56 -22816.56 -22946.993 -22946.993 252.33187 252.33187 60618.759 60618.759 256.84433 256.84433 11000 -22818.597 -22818.597 -22950.666 -22950.666 255.49609 255.49609 60605.865 60605.865 251.38418 251.38418 Loop time of 191.672 on 1 procs for 1000 steps with 4000 atoms Performance: 0.451 ns/day, 53.242 hours/ns, 5.217 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 | 190.72 | 190.72 | 190.72 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.14983 | 0.14983 | 0.14983 | 0.0 | 0.08 Output | 0.00021797 | 0.00021797 | 0.00021797 | 0.0 | 0.00 Modify | 0.70208 | 0.70208 | 0.70208 | 0.0 | 0.37 Other | | 0.09847 | | | 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: 311966.0 ave 311966 max 311966 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311966 Ave neighs/atom = 77.991500 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 = 253.199934559892, Press = -3.00543679476961 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 -22818.597 -22818.597 -22950.666 -22950.666 255.49609 255.49609 60605.865 60605.865 251.38418 251.38418 12000 -22820.183 -22820.183 -22950.31 -22950.31 251.73842 251.73842 60639.358 60639.358 -641.0932 -641.0932 Loop time of 193.297 on 1 procs for 1000 steps with 4000 atoms Performance: 0.447 ns/day, 53.694 hours/ns, 5.173 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 | 192.34 | 192.34 | 192.34 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.15222 | 0.15222 | 0.15222 | 0.0 | 0.08 Output | 0.00017749 | 0.00017749 | 0.00017749 | 0.0 | 0.00 Modify | 0.70922 | 0.70922 | 0.70922 | 0.0 | 0.37 Other | | 0.09855 | | | 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: 311970.0 ave 311970 max 311970 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311970 Ave neighs/atom = 77.992500 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 = 253.240676405407, Press = 0.697860649989379 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 -22820.183 -22820.183 -22950.31 -22950.31 251.73842 251.73842 60639.358 60639.358 -641.0932 -641.0932 13000 -22812.714 -22812.714 -22945.307 -22945.307 256.5113 256.5113 60556.271 60556.271 2245.154 2245.154 Loop time of 191.052 on 1 procs for 1000 steps with 4000 atoms Performance: 0.452 ns/day, 53.070 hours/ns, 5.234 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 | 190.11 | 190.11 | 190.11 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.14989 | 0.14989 | 0.14989 | 0.0 | 0.08 Output | 0.00017498 | 0.00017498 | 0.00017498 | 0.0 | 0.00 Modify | 0.699 | 0.699 | 0.699 | 0.0 | 0.37 Other | | 0.09733 | | | 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 = 253.401760750742, Press = -0.936786184847796 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 -22812.714 -22812.714 -22945.307 -22945.307 256.5113 256.5113 60556.271 60556.271 2245.154 2245.154 14000 -22819.604 -22819.604 -22952.804 -22952.804 257.68577 257.68577 60653.382 60653.382 -1314.1501 -1314.1501 Loop time of 193.514 on 1 procs for 1000 steps with 4000 atoms Performance: 0.446 ns/day, 53.754 hours/ns, 5.168 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 | 192.56 | 192.56 | 192.56 | 0.0 | 99.51 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.15156 | 0.15156 | 0.15156 | 0.0 | 0.08 Output | 0.00017752 | 0.00017752 | 0.00017752 | 0.0 | 0.00 Modify | 0.70573 | 0.70573 | 0.70573 | 0.0 | 0.36 Other | | 0.09837 | | | 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: 311966.0 ave 311966 max 311966 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311966 Ave neighs/atom = 77.991500 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 = 253.456480695402, Press = -3.29742069656204 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 -22819.604 -22819.604 -22952.804 -22952.804 257.68577 257.68577 60653.382 60653.382 -1314.1501 -1314.1501 15000 -22820.287 -22820.287 -22951.68 -22951.68 254.18907 254.18907 60611.045 60611.045 19.762665 19.762665 Loop time of 193.554 on 1 procs for 1000 steps with 4000 atoms Performance: 0.446 ns/day, 53.765 hours/ns, 5.167 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 | 192.59 | 192.59 | 192.59 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.15193 | 0.15193 | 0.15193 | 0.0 | 0.08 Output | 0.00047611 | 0.00047611 | 0.00047611 | 0.0 | 0.00 Modify | 0.71183 | 0.71183 | 0.71183 | 0.0 | 0.37 Other | | 0.09868 | | | 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 = 253.465517423709, Press = 0.786774440905962 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 -22820.287 -22820.287 -22951.68 -22951.68 254.18907 254.18907 60611.045 60611.045 19.762665 19.762665 16000 -22820.253 -22820.253 -22950.652 -22950.652 252.26561 252.26561 60591.05 60591.05 839.21485 839.21485 Loop time of 190.776 on 1 procs for 1000 steps with 4000 atoms Performance: 0.453 ns/day, 52.993 hours/ns, 5.242 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 | 189.83 | 189.83 | 189.83 | 0.0 | 99.51 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.1501 | 0.1501 | 0.1501 | 0.0 | 0.08 Output | 0.00017791 | 0.00017791 | 0.00017791 | 0.0 | 0.00 Modify | 0.69602 | 0.69602 | 0.69602 | 0.0 | 0.36 Other | | 0.09785 | | | 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: 311978.0 ave 311978 max 311978 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311978 Ave neighs/atom = 77.994500 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 = 253.28020589888, Press = -1.92782722888815 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 -22820.253 -22820.253 -22950.652 -22950.652 252.26561 252.26561 60591.05 60591.05 839.21485 839.21485 17000 -22819.965 -22819.965 -22952.15 -22952.15 255.72242 255.72242 60636.394 60636.394 -668.06133 -668.06133 Loop time of 189.485 on 1 procs for 1000 steps with 4000 atoms Performance: 0.456 ns/day, 52.635 hours/ns, 5.277 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 | 188.54 | 188.54 | 188.54 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.1502 | 0.1502 | 0.1502 | 0.0 | 0.08 Output | 0.00017756 | 0.00017756 | 0.00017756 | 0.0 | 0.00 Modify | 0.69499 | 0.69499 | 0.69499 | 0.0 | 0.37 Other | | 0.09926 | | | 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: 311982.0 ave 311982 max 311982 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311982 Ave neighs/atom = 77.995500 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 = 253.278550605344, Press = -0.440898520572761 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 -22819.965 -22819.965 -22952.15 -22952.15 255.72242 255.72242 60636.394 60636.394 -668.06133 -668.06133 18000 -22817.698 -22817.698 -22949.36 -22949.36 254.70846 254.70846 60604.736 60604.736 454.30608 454.30608 Loop time of 191.468 on 1 procs for 1000 steps with 4000 atoms Performance: 0.451 ns/day, 53.186 hours/ns, 5.223 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 | 190.52 | 190.52 | 190.52 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.14985 | 0.14985 | 0.14985 | 0.0 | 0.08 Output | 0.0001737 | 0.0001737 | 0.0001737 | 0.0 | 0.00 Modify | 0.70171 | 0.70171 | 0.70171 | 0.0 | 0.37 Other | | 0.09768 | | | 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: 311976.0 ave 311976 max 311976 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311976 Ave neighs/atom = 77.994000 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 = 253.258147067539, Press = -0.538834683529463 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 -22817.698 -22817.698 -22949.36 -22949.36 254.70846 254.70846 60604.736 60604.736 454.30608 454.30608 19000 -22822.397 -22822.397 -22951.552 -22951.552 249.8583 249.8583 60627.428 60627.428 -506.79516 -506.79516 Loop time of 194.123 on 1 procs for 1000 steps with 4000 atoms Performance: 0.445 ns/day, 53.923 hours/ns, 5.151 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 | 193.16 | 193.16 | 193.16 | 0.0 | 99.51 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.15127 | 0.15127 | 0.15127 | 0.0 | 0.08 Output | 0.00018313 | 0.00018313 | 0.00018313 | 0.0 | 0.00 Modify | 0.70997 | 0.70997 | 0.70997 | 0.0 | 0.37 Other | | 0.09862 | | | 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 = 253.292328452514, Press = -0.682802559481962 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 -22822.397 -22822.397 -22951.552 -22951.552 249.8583 249.8583 60627.428 60627.428 -506.79516 -506.79516 20000 -22821.554 -22821.554 -22952.499 -22952.499 253.32149 253.32149 60580.527 60580.527 795.4008 795.4008 Loop time of 195.405 on 1 procs for 1000 steps with 4000 atoms Performance: 0.442 ns/day, 54.279 hours/ns, 5.118 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 | 194.44 | 194.44 | 194.44 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.153 | 0.153 | 0.153 | 0.0 | 0.08 Output | 0.00017556 | 0.00017556 | 0.00017556 | 0.0 | 0.00 Modify | 0.7168 | 0.7168 | 0.7168 | 0.0 | 0.37 Other | | 0.09924 | | | 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: 311974.0 ave 311974 max 311974 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311974 Ave neighs/atom = 77.993500 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 = 253.259299452302, Press = 1.170949112167 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 -22821.554 -22821.554 -22952.499 -22952.499 253.32149 253.32149 60580.527 60580.527 795.4008 795.4008 21000 -22819.284 -22819.284 -22950.717 -22950.717 254.26494 254.26494 60626.859 60626.859 -319.17215 -319.17215 Loop time of 190.488 on 1 procs for 1000 steps with 4000 atoms Performance: 0.454 ns/day, 52.913 hours/ns, 5.250 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 | 189.55 | 189.55 | 189.55 | 0.0 | 99.51 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.14875 | 0.14875 | 0.14875 | 0.0 | 0.08 Output | 0.00017878 | 0.00017878 | 0.00017878 | 0.0 | 0.00 Modify | 0.69479 | 0.69479 | 0.69479 | 0.0 | 0.36 Other | | 0.09784 | | | 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: 311984.0 ave 311984 max 311984 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311984 Ave neighs/atom = 77.996000 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 = 253.143882228426, Press = -2.49498575825488 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 -22819.284 -22819.284 -22950.717 -22950.717 254.26494 254.26494 60626.859 60626.859 -319.17215 -319.17215 22000 -22819.988 -22819.988 -22948.446 -22948.446 248.5118 248.5118 60647.337 60647.337 -828.91646 -828.91646 Loop time of 191.9 on 1 procs for 1000 steps with 4000 atoms Performance: 0.450 ns/day, 53.306 hours/ns, 5.211 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 | 190.95 | 190.95 | 190.95 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.15028 | 0.15028 | 0.15028 | 0.0 | 0.08 Output | 0.00017806 | 0.00017806 | 0.00017806 | 0.0 | 0.00 Modify | 0.70228 | 0.70228 | 0.70228 | 0.0 | 0.37 Other | | 0.09783 | | | 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: 311980.0 ave 311980 max 311980 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311980 Ave neighs/atom = 77.995000 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 = 253.16129172821, Press = 1.22308504660717 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 -22819.988 -22819.988 -22948.446 -22948.446 248.5118 248.5118 60647.337 60647.337 -828.91646 -828.91646 23000 -22817.866 -22817.866 -22949.768 -22949.768 255.17271 255.17271 60571.885 60571.885 1396.9821 1396.9821 Loop time of 193.408 on 1 procs for 1000 steps with 4000 atoms Performance: 0.447 ns/day, 53.724 hours/ns, 5.170 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 | 192.44 | 192.44 | 192.44 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.15254 | 0.15254 | 0.15254 | 0.0 | 0.08 Output | 0.00018364 | 0.00018364 | 0.00018364 | 0.0 | 0.00 Modify | 0.71505 | 0.71505 | 0.71505 | 0.0 | 0.37 Other | | 0.09935 | | | 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: 311978.0 ave 311978 max 311978 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311978 Ave neighs/atom = 77.994500 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 = 253.133763617647, Press = -0.374190081245376 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 -22817.866 -22817.866 -22949.768 -22949.768 255.17271 255.17271 60571.885 60571.885 1396.9821 1396.9821 24000 -22820.649 -22820.649 -22951.312 -22951.312 252.77699 252.77699 60635.793 60635.793 -685.64944 -685.64944 Loop time of 194.394 on 1 procs for 1000 steps with 4000 atoms Performance: 0.444 ns/day, 53.998 hours/ns, 5.144 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 | 193.42 | 193.42 | 193.42 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.15396 | 0.15396 | 0.15396 | 0.0 | 0.08 Output | 0.00021794 | 0.00021794 | 0.00021794 | 0.0 | 0.00 Modify | 0.72225 | 0.72225 | 0.72225 | 0.0 | 0.37 Other | | 0.1007 | | | 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: 311986.0 ave 311986 max 311986 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311986 Ave neighs/atom = 77.996500 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 = 253.05870981408, Press = -0.585694984677946 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 -22820.649 -22820.649 -22951.312 -22951.312 252.77699 252.77699 60635.793 60635.793 -685.64944 -685.64944 25000 -22820.563 -22820.563 -22952.284 -22952.284 254.82268 254.82268 60590.401 60590.401 581.02363 581.02363 Loop time of 193.482 on 1 procs for 1000 steps with 4000 atoms Performance: 0.447 ns/day, 53.745 hours/ns, 5.168 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 | 192.52 | 192.52 | 192.52 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.15096 | 0.15096 | 0.15096 | 0.0 | 0.08 Output | 0.00022594 | 0.00022594 | 0.00022594 | 0.0 | 0.00 Modify | 0.71032 | 0.71032 | 0.71032 | 0.0 | 0.37 Other | | 0.09873 | | | 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: 311980.0 ave 311980 max 311980 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311980 Ave neighs/atom = 77.995000 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_T253.15.out" else "print 'not_converged' file output/vol_T253.15.out" print '${V}' file output/vol_T253.15.out 60616.6673748625 print "LAMMPS calculation completed" LAMMPS calculation completed quit 0