# 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 2.880978018045426*${_u_distance} variable latticeconst_converted equal 2.880978018045426*1 lattice bcc ${latticeconst_converted} lattice bcc 2.88097801804543 Lattice spacing in x,y,z = 2.8809780 2.8809780 2.8809780 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 (28.809780 28.809780 28.809780) 1 by 1 by 1 MPI processor grid create_atoms 1 box Created 2000 atoms using lattice units in orthogonal box = (0.0000000 0.0000000 0.0000000) to (28.809780 28.809780 28.809780) create_atoms CPU = 0.002 seconds variable mass_converted equal 51.9961*${_u_mass} variable mass_converted equal 51.9961*1 kim_interactions Cr WARNING: 'kim_' has been renamed to 'kim '. Please update your input. kim interactions Cr #=== BEGIN kim interactions ================================== pair_style kim MEAM_LAMMPS_ChoiKimSeol_2017_CoCr__MO_410167849923_001 pair_coeff * * Cr #=== END kim interactions ==================================== mass 1 ${mass_converted} mass 1 51.9961 # initial volume variable v equal vol # assign formula variable V0 equal ${v} # evaluate initial value variable V0 equal 23912.2164838903 variable V0_metal equal ${V0}/(${_u_distance}*${_u_distance}*${_u_distance}) variable V0_metal equal 23912.2164838903/(${_u_distance}*${_u_distance}*${_u_distance}) variable V0_metal equal 23912.2164838903/(1*${_u_distance}*${_u_distance}) variable V0_metal equal 23912.2164838903/(1*1*${_u_distance}) variable V0_metal equal 23912.2164838903/(1*1*1) variable V0_metal_times1000 equal ${V0_metal}*1000 variable V0_metal_times1000 equal 23912.2164838903*1000 print "Initial system volume: ${V0_metal} Angstroms^3" Initial system volume: 23912.2164838903 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 293.15*${_u_temperature} variable temp_converted equal 293.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 293.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 293.15 ${temp_converted} ${Tdamp_converted} iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 293.15 293.15 ${Tdamp_converted} iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 293.15 293.15 0.1 iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 293.15 293.15 0.1 iso 0 ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 293.15 293.15 0.1 iso 0 0 ${Pdamp_converted} fix ensemble all npt temp 293.15 293.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 "293.15 - 0.2" variable T_up equal "293.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_410167849923_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 ghost atom cutoff = 6 binsize = 3, bins = 10 10 10 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 0 -8124.2527 -8124.2527 -8200 -8200 293.15 293.15 23912.216 23912.216 3383.4925 3383.4925 1000 -8045.708 -8045.708 -8121.8384 -8121.8384 294.63219 294.63219 24200.826 24200.826 -2946.9127 -2946.9127 Loop time of 113.737 on 1 procs for 1000 steps with 2000 atoms Performance: 0.760 ns/day, 31.594 hours/ns, 8.792 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 | 113.16 | 113.16 | 113.16 | 0.0 | 99.49 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.10726 | 0.10726 | 0.10726 | 0.0 | 0.09 Output | 0.00019578 | 0.00019578 | 0.00019578 | 0.0 | 0.00 Modify | 0.40425 | 0.40425 | 0.40425 | 0.0 | 0.36 Other | | 0.06402 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119.00 ave 4119 max 4119 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: 128000.0 ave 128000 max 128000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 128000 Ave neighs/atom = 64.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 ghost atom cutoff = 6 binsize = 3, bins = 10 10 10 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 1000 -8045.708 -8045.708 -8121.8384 -8121.8384 294.63219 294.63219 24200.826 24200.826 -2946.9127 -2946.9127 2000 -8044.8375 -8044.8375 -8124.4146 -8124.4146 307.97166 307.97166 24142.102 24142.102 1344.2982 1344.2982 Loop time of 119.056 on 1 procs for 1000 steps with 2000 atoms Performance: 0.726 ns/day, 33.071 hours/ns, 8.399 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 | 118.46 | 118.46 | 118.46 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.10856 | 0.10856 | 0.10856 | 0.0 | 0.09 Output | 0.00023279 | 0.00023279 | 0.00023279 | 0.0 | 0.00 Modify | 0.41859 | 0.41859 | 0.41859 | 0.0 | 0.35 Other | | 0.06557 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119.00 ave 4119 max 4119 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: 127948.0 ave 127948 max 127948 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 127948 Ave neighs/atom = 63.974000 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 ghost atom cutoff = 6 binsize = 3, bins = 10 10 10 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 2000 -8044.8375 -8044.8375 -8124.4146 -8124.4146 307.97166 307.97166 24142.102 24142.102 1344.2982 1344.2982 3000 -8047.7457 -8047.7457 -8124.0785 -8124.0785 295.41594 295.41594 24178.527 24178.527 -1845.3096 -1845.3096 Loop time of 119.115 on 1 procs for 1000 steps with 2000 atoms Performance: 0.725 ns/day, 33.087 hours/ns, 8.395 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 | 118.52 | 118.52 | 118.52 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.10822 | 0.10822 | 0.10822 | 0.0 | 0.09 Output | 0.00018986 | 0.00018986 | 0.00018986 | 0.0 | 0.00 Modify | 0.41661 | 0.41661 | 0.41661 | 0.0 | 0.35 Other | | 0.06641 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119.00 ave 4119 max 4119 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: 127986.0 ave 127986 max 127986 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 127986 Ave neighs/atom = 63.993000 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 ghost atom cutoff = 6 binsize = 3, bins = 10 10 10 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 3000 -8047.7457 -8047.7457 -8124.0785 -8124.0785 295.41594 295.41594 24178.527 24178.527 -1845.3096 -1845.3096 4000 -8043.845 -8043.845 -8124.4005 -8124.4005 311.75789 311.75789 24152.024 24152.024 354.00448 354.00448 Loop time of 120.018 on 1 procs for 1000 steps with 2000 atoms Performance: 0.720 ns/day, 33.338 hours/ns, 8.332 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 | 119.42 | 119.42 | 119.42 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.10976 | 0.10976 | 0.10976 | 0.0 | 0.09 Output | 0.00019128 | 0.00019128 | 0.00019128 | 0.0 | 0.00 Modify | 0.42434 | 0.42434 | 0.42434 | 0.0 | 0.35 Other | | 0.06723 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4118.00 ave 4118 max 4118 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: 127962.0 ave 127962 max 127962 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 127962 Ave neighs/atom = 63.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 = 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 ghost atom cutoff = 6 binsize = 3, bins = 10 10 10 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 4000 -8043.845 -8043.845 -8124.4005 -8124.4005 311.75789 311.75789 24152.024 24152.024 354.00448 354.00448 5000 -8048.463 -8048.463 -8121.569 -8121.569 282.92751 282.92751 24162.982 24162.982 -90.11951 -90.11951 Loop time of 112.485 on 1 procs for 1000 steps with 2000 atoms Performance: 0.768 ns/day, 31.246 hours/ns, 8.890 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 | 111.93 | 111.93 | 111.93 | 0.0 | 99.51 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.10287 | 0.10287 | 0.10287 | 0.0 | 0.09 Output | 0.00015326 | 0.00015326 | 0.00015326 | 0.0 | 0.00 Modify | 0.38784 | 0.38784 | 0.38784 | 0.0 | 0.34 Other | | 0.0633 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119.00 ave 4119 max 4119 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: 127962.0 ave 127962 max 127962 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 127962 Ave neighs/atom = 63.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 = 290.933033202602, Press = -221.597616405018 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 ghost atom cutoff = 6 binsize = 3, bins = 10 10 10 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 5000 -8048.463 -8048.463 -8121.569 -8121.569 282.92751 282.92751 24162.982 24162.982 -90.11951 -90.11951 6000 -8045.2071 -8045.2071 -8122.8389 -8122.8389 300.4431 300.4431 24145.115 24145.115 918.15183 918.15183 Loop time of 102.124 on 1 procs for 1000 steps with 2000 atoms Performance: 0.846 ns/day, 28.368 hours/ns, 9.792 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 | 101.62 | 101.62 | 101.62 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.095373 | 0.095373 | 0.095373 | 0.0 | 0.09 Output | 0.0001581 | 0.0001581 | 0.0001581 | 0.0 | 0.00 Modify | 0.35277 | 0.35277 | 0.35277 | 0.0 | 0.35 Other | | 0.05992 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119.00 ave 4119 max 4119 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: 127942.0 ave 127942 max 127942 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 127942 Ave neighs/atom = 63.971000 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 = 292.504653397178, Press = -7.42973413523656 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 ghost atom cutoff = 6 binsize = 3, bins = 10 10 10 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 6000 -8045.2071 -8045.2071 -8122.8389 -8122.8389 300.4431 300.4431 24145.115 24145.115 918.15183 918.15183 7000 -8046.9158 -8046.9158 -8120.9246 -8120.9246 286.42179 286.42179 24219.553 24219.553 -4696.6896 -4696.6896 Loop time of 118.362 on 1 procs for 1000 steps with 2000 atoms Performance: 0.730 ns/day, 32.878 hours/ns, 8.449 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 | 117.76 | 117.76 | 117.76 | 0.0 | 99.49 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.10744 | 0.10744 | 0.10744 | 0.0 | 0.09 Output | 0.00021717 | 0.00021717 | 0.00021717 | 0.0 | 0.00 Modify | 0.43042 | 0.43042 | 0.43042 | 0.0 | 0.36 Other | | 0.06502 | | | 0.05 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119.00 ave 4119 max 4119 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: 127946.0 ave 127946 max 127946 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 127946 Ave neighs/atom = 63.973000 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 = 291.274344467968, Press = -21.0464262167075 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 ghost atom cutoff = 6 binsize = 3, bins = 10 10 10 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 7000 -8046.9158 -8046.9158 -8120.9246 -8120.9246 286.42179 286.42179 24219.553 24219.553 -4696.6896 -4696.6896 8000 -8045.4272 -8045.4272 -8120.2693 -8120.2693 289.64659 289.64659 24125.177 24125.177 3036.7396 3036.7396 Loop time of 105.222 on 1 procs for 1000 steps with 2000 atoms Performance: 0.821 ns/day, 29.228 hours/ns, 9.504 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 | 104.69 | 104.69 | 104.69 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.097545 | 0.097545 | 0.097545 | 0.0 | 0.09 Output | 0.00015425 | 0.00015425 | 0.00015425 | 0.0 | 0.00 Modify | 0.36853 | 0.36853 | 0.36853 | 0.0 | 0.35 Other | | 0.06109 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119.00 ave 4119 max 4119 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: 127958.0 ave 127958 max 127958 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 127958 Ave neighs/atom = 63.979000 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 = 291.943266976616, Press = 6.72298177126145 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 ghost atom cutoff = 6 binsize = 3, bins = 10 10 10 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 8000 -8045.4272 -8045.4272 -8120.2693 -8120.2693 289.64659 289.64659 24125.177 24125.177 3036.7396 3036.7396 9000 -8044.3937 -8044.3937 -8119.495 -8119.495 290.64973 290.64973 24183.023 24183.023 -1318.1263 -1318.1263 Loop time of 119.377 on 1 procs for 1000 steps with 2000 atoms Performance: 0.724 ns/day, 33.160 hours/ns, 8.377 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 | 118.76 | 118.76 | 118.76 | 0.0 | 99.49 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.10894 | 0.10894 | 0.10894 | 0.0 | 0.09 Output | 0.00015844 | 0.00015844 | 0.00015844 | 0.0 | 0.00 Modify | 0.43738 | 0.43738 | 0.43738 | 0.0 | 0.37 Other | | 0.06568 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119.00 ave 4119 max 4119 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: 127952.0 ave 127952 max 127952 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 127952 Ave neighs/atom = 63.976000 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 = 292.087784065411, Press = -11.2111869487754 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 ghost atom cutoff = 6 binsize = 3, bins = 10 10 10 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 9000 -8044.3937 -8044.3937 -8119.495 -8119.495 290.64973 290.64973 24183.023 24183.023 -1318.1263 -1318.1263 10000 -8046.9057 -8046.9057 -8122.7984 -8122.7984 293.71275 293.71275 24141.53 24141.53 1259.8139 1259.8139 Loop time of 119.571 on 1 procs for 1000 steps with 2000 atoms Performance: 0.723 ns/day, 33.214 hours/ns, 8.363 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 | 118.96 | 118.96 | 118.96 | 0.0 | 99.49 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.10827 | 0.10827 | 0.10827 | 0.0 | 0.09 Output | 0.00015413 | 0.00015413 | 0.00015413 | 0.0 | 0.00 Modify | 0.43877 | 0.43877 | 0.43877 | 0.0 | 0.37 Other | | 0.06562 | | | 0.05 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119.00 ave 4119 max 4119 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: 127956.0 ave 127956 max 127956 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 127956 Ave neighs/atom = 63.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 = 292.577108433244, Press = 3.67625144953318 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 ghost atom cutoff = 6 binsize = 3, bins = 10 10 10 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 10000 -8046.9057 -8046.9057 -8122.7984 -8122.7984 293.71275 293.71275 24141.53 24141.53 1259.8139 1259.8139 11000 -8043.4839 -8043.4839 -8119.1161 -8119.1161 292.70439 292.70439 24186.146 24186.146 -1451.2472 -1451.2472 Loop time of 119.515 on 1 procs for 1000 steps with 2000 atoms Performance: 0.723 ns/day, 33.199 hours/ns, 8.367 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 | 118.9 | 118.9 | 118.9 | 0.0 | 99.49 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.10834 | 0.10834 | 0.10834 | 0.0 | 0.09 Output | 0.00021551 | 0.00021551 | 0.00021551 | 0.0 | 0.00 Modify | 0.4388 | 0.4388 | 0.4388 | 0.0 | 0.37 Other | | 0.06608 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119.00 ave 4119 max 4119 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: 127950.0 ave 127950 max 127950 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 127950 Ave neighs/atom = 63.975000 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 = 293.007095309663, Press = -4.92196738835384 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 ghost atom cutoff = 6 binsize = 3, bins = 10 10 10 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 11000 -8043.4839 -8043.4839 -8119.1161 -8119.1161 292.70439 292.70439 24186.146 24186.146 -1451.2472 -1451.2472 12000 -8048.8788 -8048.8788 -8122.9156 -8122.9156 286.53024 286.53024 24135.174 24135.174 1786.3509 1786.3509 Loop time of 119.767 on 1 procs for 1000 steps with 2000 atoms Performance: 0.721 ns/day, 33.268 hours/ns, 8.350 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 | 119.15 | 119.15 | 119.15 | 0.0 | 99.48 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.10937 | 0.10937 | 0.10937 | 0.0 | 0.09 Output | 0.00019713 | 0.00019713 | 0.00019713 | 0.0 | 0.00 Modify | 0.44301 | 0.44301 | 0.44301 | 0.0 | 0.37 Other | | 0.06615 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119.00 ave 4119 max 4119 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: 127948.0 ave 127948 max 127948 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 127948 Ave neighs/atom = 63.974000 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 = 293.048167172909, Press = 6.5382989982553 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 ghost atom cutoff = 6 binsize = 3, bins = 10 10 10 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 12000 -8048.8788 -8048.8788 -8122.9156 -8122.9156 286.53024 286.53024 24135.174 24135.174 1786.3509 1786.3509 13000 -8045.4818 -8045.4818 -8121.3733 -8121.3733 293.70779 293.70779 24179.609 24179.609 -1489.326 -1489.326 Loop time of 102.335 on 1 procs for 1000 steps with 2000 atoms Performance: 0.844 ns/day, 28.426 hours/ns, 9.772 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 | 101.83 | 101.83 | 101.83 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.095273 | 0.095273 | 0.095273 | 0.0 | 0.09 Output | 0.00015449 | 0.00015449 | 0.00015449 | 0.0 | 0.00 Modify | 0.35389 | 0.35389 | 0.35389 | 0.0 | 0.35 Other | | 0.06071 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119.00 ave 4119 max 4119 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: 127946.0 ave 127946 max 127946 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 127946 Ave neighs/atom = 63.973000 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 = 292.828865740952, Press = -9.38420528882142 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 ghost atom cutoff = 6 binsize = 3, bins = 10 10 10 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 13000 -8045.4818 -8045.4818 -8121.3733 -8121.3733 293.70779 293.70779 24179.609 24179.609 -1489.326 -1489.326 14000 -8046.8846 -8046.8846 -8123.3605 -8123.3605 295.96944 295.96944 24153.636 24153.636 366.08764 366.08764 Loop time of 101.511 on 1 procs for 1000 steps with 2000 atoms Performance: 0.851 ns/day, 28.197 hours/ns, 9.851 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 | 101 | 101 | 101 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.0949 | 0.0949 | 0.0949 | 0.0 | 0.09 Output | 0.00018802 | 0.00018802 | 0.00018802 | 0.0 | 0.00 Modify | 0.35166 | 0.35166 | 0.35166 | 0.0 | 0.35 Other | | 0.06042 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119.00 ave 4119 max 4119 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: 127938.0 ave 127938 max 127938 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 127938 Ave neighs/atom = 63.969000 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 = 292.809514196427, Press = 6.19793066752117 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 ghost atom cutoff = 6 binsize = 3, bins = 10 10 10 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 14000 -8046.8846 -8046.8846 -8123.3605 -8123.3605 295.96944 295.96944 24153.636 24153.636 366.08764 366.08764 15000 -8046.103 -8046.103 -8123.2859 -8123.2859 298.70576 298.70576 24160.267 24160.267 -157.89749 -157.89749 Loop time of 101.47 on 1 procs for 1000 steps with 2000 atoms Performance: 0.851 ns/day, 28.186 hours/ns, 9.855 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 | 100.96 | 100.96 | 100.96 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.094518 | 0.094518 | 0.094518 | 0.0 | 0.09 Output | 0.00015426 | 0.00015426 | 0.00015426 | 0.0 | 0.00 Modify | 0.35175 | 0.35175 | 0.35175 | 0.0 | 0.35 Other | | 0.06015 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119.00 ave 4119 max 4119 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: 127954.0 ave 127954 max 127954 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 127954 Ave neighs/atom = 63.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 = 293.026510842404, Press = -6.2477011435766 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 ghost atom cutoff = 6 binsize = 3, bins = 10 10 10 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 15000 -8046.103 -8046.103 -8123.2859 -8123.2859 298.70576 298.70576 24160.267 24160.267 -157.89749 -157.89749 16000 -8042.0418 -8042.0418 -8118.3716 -8118.3716 295.40432 295.40432 24179.008 24179.008 -599.72756 -599.72756 Loop time of 101.452 on 1 procs for 1000 steps with 2000 atoms Performance: 0.852 ns/day, 28.181 hours/ns, 9.857 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 | 100.95 | 100.95 | 100.95 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.09458 | 0.09458 | 0.09458 | 0.0 | 0.09 Output | 0.00015244 | 0.00015244 | 0.00015244 | 0.0 | 0.00 Modify | 0.35071 | 0.35071 | 0.35071 | 0.0 | 0.35 Other | | 0.06074 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119.00 ave 4119 max 4119 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: 127958.0 ave 127958 max 127958 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 127958 Ave neighs/atom = 63.979000 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 = 293.136742922294, Press = 5.4391811762672 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 ghost atom cutoff = 6 binsize = 3, bins = 10 10 10 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 16000 -8042.0418 -8042.0418 -8118.3716 -8118.3716 295.40432 295.40432 24179.008 24179.008 -599.72756 -599.72756 17000 -8047.3658 -8047.3658 -8123.7206 -8123.7206 295.50089 295.50089 24148.931 24148.931 541.54267 541.54267 Loop time of 108.949 on 1 procs for 1000 steps with 2000 atoms Performance: 0.793 ns/day, 30.264 hours/ns, 9.179 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 | 108.4 | 108.4 | 108.4 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.10075 | 0.10075 | 0.10075 | 0.0 | 0.09 Output | 0.00021384 | 0.00021384 | 0.00021384 | 0.0 | 0.00 Modify | 0.38645 | 0.38645 | 0.38645 | 0.0 | 0.35 Other | | 0.06185 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119.00 ave 4119 max 4119 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: 127912.0 ave 127912 max 127912 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 127912 Ave neighs/atom = 63.956000 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 = 293.285482325559, Press = -4.9338614388353 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 ghost atom cutoff = 6 binsize = 3, bins = 10 10 10 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 17000 -8047.3658 -8047.3658 -8123.7206 -8123.7206 295.50089 295.50089 24148.931 24148.931 541.54267 541.54267 18000 -8043.4316 -8043.4316 -8120.6062 -8120.6062 298.67354 298.67354 24173.822 24173.822 -655.29264 -655.29264 Loop time of 119.82 on 1 procs for 1000 steps with 2000 atoms Performance: 0.721 ns/day, 33.283 hours/ns, 8.346 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 | 119.21 | 119.21 | 119.21 | 0.0 | 99.49 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.10923 | 0.10923 | 0.10923 | 0.0 | 0.09 Output | 0.00015537 | 0.00015537 | 0.00015537 | 0.0 | 0.00 Modify | 0.4381 | 0.4381 | 0.4381 | 0.0 | 0.37 Other | | 0.06596 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119.00 ave 4119 max 4119 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: 127970.0 ave 127970 max 127970 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 127970 Ave neighs/atom = 63.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 = 293.435806303066, Press = 1.44349906209016 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 ghost atom cutoff = 6 binsize = 3, bins = 10 10 10 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 18000 -8043.4316 -8043.4316 -8120.6062 -8120.6062 298.67354 298.67354 24173.822 24173.822 -655.29264 -655.29264 19000 -8046.1048 -8046.1048 -8122.2402 -8122.2402 294.65181 294.65181 24157.083 24157.083 253.70238 253.70238 Loop time of 118.716 on 1 procs for 1000 steps with 2000 atoms Performance: 0.728 ns/day, 32.977 hours/ns, 8.423 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 | 118.11 | 118.11 | 118.11 | 0.0 | 99.49 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.10808 | 0.10808 | 0.10808 | 0.0 | 0.09 Output | 0.00018751 | 0.00018751 | 0.00018751 | 0.0 | 0.00 Modify | 0.43645 | 0.43645 | 0.43645 | 0.0 | 0.37 Other | | 0.06546 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119.00 ave 4119 max 4119 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: 127946.0 ave 127946 max 127946 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 127946 Ave neighs/atom = 63.973000 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 = 293.302298139793, Press = -1.86260677192351 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 ghost atom cutoff = 6 binsize = 3, bins = 10 10 10 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 19000 -8046.1048 -8046.1048 -8122.2402 -8122.2402 294.65181 294.65181 24157.083 24157.083 253.70238 253.70238 20000 -8047.4028 -8047.4028 -8121.0492 -8121.0492 285.01924 285.01924 24162.678 24162.678 -152.98306 -152.98306 Loop time of 119.128 on 1 procs for 1000 steps with 2000 atoms Performance: 0.725 ns/day, 33.091 hours/ns, 8.394 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 | 118.51 | 118.51 | 118.51 | 0.0 | 99.48 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.10827 | 0.10827 | 0.10827 | 0.0 | 0.09 Output | 0.00020021 | 0.00020021 | 0.00020021 | 0.0 | 0.00 Modify | 0.44113 | 0.44113 | 0.44113 | 0.0 | 0.37 Other | | 0.06569 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119.00 ave 4119 max 4119 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: 127958.0 ave 127958 max 127958 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 127958 Ave neighs/atom = 63.979000 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 = 293.191946220261, Press = 0.280627994737257 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 ghost atom cutoff = 6 binsize = 3, bins = 10 10 10 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 20000 -8047.4028 -8047.4028 -8121.0492 -8121.0492 285.01924 285.01924 24162.678 24162.678 -152.98306 -152.98306 21000 -8043.8371 -8043.8371 -8121.7941 -8121.7941 301.70187 301.70187 24134.7 24134.7 2226.0632 2226.0632 Loop time of 118.766 on 1 procs for 1000 steps with 2000 atoms Performance: 0.727 ns/day, 32.990 hours/ns, 8.420 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 | 118.15 | 118.15 | 118.15 | 0.0 | 99.48 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.10882 | 0.10882 | 0.10882 | 0.0 | 0.09 Output | 0.00019933 | 0.00019933 | 0.00019933 | 0.0 | 0.00 Modify | 0.44152 | 0.44152 | 0.44152 | 0.0 | 0.37 Other | | 0.06519 | | | 0.05 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119.00 ave 4119 max 4119 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: 127942.0 ave 127942 max 127942 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 127942 Ave neighs/atom = 63.971000 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 = 293.125577818049, Press = -0.800003858177119 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 ghost atom cutoff = 6 binsize = 3, bins = 10 10 10 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 21000 -8043.8371 -8043.8371 -8121.7941 -8121.7941 301.70187 301.70187 24134.7 24134.7 2226.0632 2226.0632 22000 -8044.4633 -8044.4633 -8119.7731 -8119.7731 291.45694 291.45694 24216.309 24216.309 -3880.9549 -3880.9549 Loop time of 119.447 on 1 procs for 1000 steps with 2000 atoms Performance: 0.723 ns/day, 33.180 hours/ns, 8.372 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 | 118.83 | 118.83 | 118.83 | 0.0 | 99.48 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.10939 | 0.10939 | 0.10939 | 0.0 | 0.09 Output | 0.00015461 | 0.00015461 | 0.00015461 | 0.0 | 0.00 Modify | 0.44543 | 0.44543 | 0.44543 | 0.0 | 0.37 Other | | 0.06619 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119.00 ave 4119 max 4119 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: 127952.0 ave 127952 max 127952 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 127952 Ave neighs/atom = 63.976000 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 = 293.156042906722, Press = -1.79485482417518 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 ghost atom cutoff = 6 binsize = 3, bins = 10 10 10 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 22000 -8044.4633 -8044.4633 -8119.7731 -8119.7731 291.45694 291.45694 24216.309 24216.309 -3880.9549 -3880.9549 23000 -8046.0112 -8046.0112 -8123.1305 -8123.1305 298.45971 298.45971 24130.979 24130.979 2212.7717 2212.7717 Loop time of 102.217 on 1 procs for 1000 steps with 2000 atoms Performance: 0.845 ns/day, 28.393 hours/ns, 9.783 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 | 101.71 | 101.71 | 101.71 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.095707 | 0.095707 | 0.095707 | 0.0 | 0.09 Output | 0.00015407 | 0.00015407 | 0.00015407 | 0.0 | 0.00 Modify | 0.35393 | 0.35393 | 0.35393 | 0.0 | 0.35 Other | | 0.05981 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119.00 ave 4119 max 4119 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: 127926.0 ave 127926 max 127926 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 127926 Ave neighs/atom = 63.963000 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 = 293.339904632838, Press = 0.922519211701235 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 ghost atom cutoff = 6 binsize = 3, bins = 10 10 10 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 23000 -8046.0112 -8046.0112 -8123.1305 -8123.1305 298.45971 298.45971 24130.979 24130.979 2212.7717 2212.7717 24000 -8044.8797 -8044.8797 -8120.8869 -8120.8869 294.15562 294.15562 24175.779 24175.779 -1004.121 -1004.121 Loop time of 107.385 on 1 procs for 1000 steps with 2000 atoms Performance: 0.805 ns/day, 29.829 hours/ns, 9.312 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 | 106.85 | 106.85 | 106.85 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.099383 | 0.099383 | 0.099383 | 0.0 | 0.09 Output | 0.00015623 | 0.00015623 | 0.00015623 | 0.0 | 0.00 Modify | 0.37793 | 0.37793 | 0.37793 | 0.0 | 0.35 Other | | 0.06235 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119.00 ave 4119 max 4119 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: 127940.0 ave 127940 max 127940 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 127940 Ave neighs/atom = 63.970000 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 = 293.327743899551, Press = -1.93816858589161 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 ghost atom cutoff = 6 binsize = 3, bins = 10 10 10 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 24000 -8044.8797 -8044.8797 -8120.8869 -8120.8869 294.15562 294.15562 24175.779 24175.779 -1004.121 -1004.121 25000 -8051.8831 -8051.8831 -8125.1603 -8125.1603 283.59015 283.59015 24144.292 24144.292 435.92983 435.92983 Loop time of 112.789 on 1 procs for 1000 steps with 2000 atoms Performance: 0.766 ns/day, 31.330 hours/ns, 8.866 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 | 112.21 | 112.21 | 112.21 | 0.0 | 99.49 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.10369 | 0.10369 | 0.10369 | 0.0 | 0.09 Output | 0.00015168 | 0.00015168 | 0.00015168 | 0.0 | 0.00 Modify | 0.40662 | 0.40662 | 0.40662 | 0.0 | 0.36 Other | | 0.06382 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119.00 ave 4119 max 4119 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: 127942.0 ave 127942 max 127942 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 127942 Ave neighs/atom = 63.971000 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_T293.15.out" else "print 'not_converged' file output/vol_T293.15.out" print '${V}' file output/vol_T293.15.out 24160.9892469373 print "LAMMPS calculation completed" LAMMPS calculation completed quit 0