# Variables that can be adjusted by kim-lammps-preprocessor to switch unit sets for # Simulator Models variable _u_distance equal 1.0 variable _u_energy equal 1.0 variable _u_mass equal 1.0 variable _u_time equal 1.0 variable _u_pressure equal 1.0 variable _u_temperature equal 1.0 # This line may be swapped out by kim-lammps-preprocessor if running against a Simulator # Model whose atom_style is not 'atomic' atom_style atomic # 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.139999896287918*${_u_distance} variable latticeconst_converted equal 3.139999896287918*1 lattice bcc ${latticeconst_converted} lattice bcc 3.13999989628792 Lattice spacing in x,y,z = 3.14 3.14 3.14 region simbox block 0 10 0 10 0 10 units lattice create_box 1 simbox Created orthogonal box = (0 0 0) to (31.4 31.4 31.4) 1 by 1 by 1 MPI processor grid create_atoms 1 box Created 2000 atoms create_atoms CPU = 0.000290871 secs variable mass_converted equal 183.84*${_u_mass} variable mass_converted equal 183.84*1 # specify which KIM Model to use pair_style kim EAM_Dynamo_BonnyTerentyev_2014EAM2_W__MO_626183701337_000 pair_coeff * * W mass 1 ${mass_converted} mass 1 183.84 # initial volume variable v equal vol # assign formula variable V0 equal ${v} # evaluate initial value variable V0 equal 30959.1409323212 variable V0_metal equal ${V0}/(${_u_distance}*${_u_distance}*${_u_distance}) variable V0_metal equal 30959.1409323212/(${_u_distance}*${_u_distance}*${_u_distance}) variable V0_metal equal 30959.1409323212/(1*${_u_distance}*${_u_distance}) variable V0_metal equal 30959.1409323212/(1*1*${_u_distance}) variable V0_metal equal 30959.1409323212/(1*1*1) variable V0_metal_times1000 equal ${V0_metal}*1000 variable V0_metal_times1000 equal 30959.1409323212*1000 print "Initial system volume: ${V0_metal} Angstroms^3" Initial system volume: 30959.1409323212 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 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7.46044 ghost atom cutoff = 7.46044 binsize = 3.73022, bins = 9 9 9 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7.46044 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.133 | 4.133 | 4.133 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 0 -17734.548 -17734.548 -17799.96 -17799.96 253.15 253.15 30959.141 30959.141 2256.7181 2256.7181 1000 -17683.349 -17683.349 -17747.52 -17747.52 248.34964 248.34964 31155.321 31155.321 1547.9938 1547.9938 Loop time of 11.7167 on 1 procs for 1000 steps with 2000 atoms Performance: 7.374 ns/day, 3.255 hours/ns, 85.348 timesteps/s 42.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 | 11.475 | 11.475 | 11.475 | 0.0 | 97.94 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.044831 | 0.044831 | 0.044831 | 0.0 | 0.38 Output | 3.2902e-05 | 3.2902e-05 | 3.2902e-05 | 0.0 | 0.00 Modify | 0.18569 | 0.18569 | 0.18569 | 0.0 | 1.58 Other | | 0.01107 | | | 0.09 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 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 Per MPI rank memory allocation (min/avg/max) = 4.133 | 4.133 | 4.133 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 1000 -17683.349 -17683.349 -17747.52 -17747.52 248.34964 248.34964 31155.321 31155.321 1547.9938 1547.9938 2000 -17676.497 -17676.497 -17744.039 -17744.039 261.39647 261.39647 31186.411 31186.411 -488.16491 -488.16491 Loop time of 12.1448 on 1 procs for 1000 steps with 2000 atoms Performance: 7.114 ns/day, 3.374 hours/ns, 82.340 timesteps/s 42.5% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 11.782 | 11.782 | 11.782 | 0.0 | 97.01 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.025082 | 0.025082 | 0.025082 | 0.0 | 0.21 Output | 3.5048e-05 | 3.5048e-05 | 3.5048e-05 | 0.0 | 0.00 Modify | 0.30686 | 0.30686 | 0.30686 | 0.0 | 2.53 Other | | 0.03109 | | | 0.26 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 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 Per MPI rank memory allocation (min/avg/max) = 4.133 | 4.133 | 4.133 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 2000 -17676.497 -17676.497 -17744.039 -17744.039 261.39647 261.39647 31186.411 31186.411 -488.16491 -488.16491 3000 -17680.268 -17680.268 -17746.385 -17746.385 255.87954 255.87954 31203.77 31203.77 -4062.0975 -4062.0975 Loop time of 12.6819 on 1 procs for 1000 steps with 2000 atoms Performance: 6.813 ns/day, 3.523 hours/ns, 78.852 timesteps/s 40.5% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 12.4 | 12.4 | 12.4 | 0.0 | 97.78 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.064869 | 0.064869 | 0.064869 | 0.0 | 0.51 Output | 3.2902e-05 | 3.2902e-05 | 3.2902e-05 | 0.0 | 0.00 Modify | 0.20605 | 0.20605 | 0.20605 | 0.0 | 1.62 Other | | 0.01107 | | | 0.09 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 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 Per MPI rank memory allocation (min/avg/max) = 4.133 | 4.133 | 4.133 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 3000 -17680.268 -17680.268 -17746.385 -17746.385 255.87954 255.87954 31203.77 31203.77 -4062.0975 -4062.0975 4000 -17680.761 -17680.761 -17744.922 -17744.922 248.31138 248.31138 31166.648 31166.648 594.60533 594.60533 Loop time of 12.596 on 1 procs for 1000 steps with 2000 atoms Performance: 6.859 ns/day, 3.499 hours/ns, 79.390 timesteps/s 41.1% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 12.354 | 12.354 | 12.354 | 0.0 | 98.08 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.044764 | 0.044764 | 0.044764 | 0.0 | 0.36 Output | 3.3855e-05 | 3.3855e-05 | 3.3855e-05 | 0.0 | 0.00 Modify | 0.18569 | 0.18569 | 0.18569 | 0.0 | 1.47 Other | | 0.0115 | | | 0.09 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 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 Per MPI rank memory allocation (min/avg/max) = 4.133 | 4.133 | 4.133 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 4000 -17680.761 -17680.761 -17744.922 -17744.922 248.31138 248.31138 31166.648 31166.648 594.60533 594.60533 5000 -17677.667 -17677.667 -17744.173 -17744.173 257.38375 257.38375 31196.896 31196.896 -2588.813 -2588.813 Loop time of 11.8818 on 1 procs for 1000 steps with 2000 atoms Performance: 7.272 ns/day, 3.300 hours/ns, 84.162 timesteps/s 43.3% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 11.517 | 11.517 | 11.517 | 0.0 | 96.93 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.085152 | 0.085152 | 0.085152 | 0.0 | 0.72 Output | 2.9087e-05 | 2.9087e-05 | 2.9087e-05 | 0.0 | 0.00 Modify | 0.26758 | 0.26758 | 0.26758 | 0.0 | 2.25 Other | | 0.01186 | | | 0.10 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 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 = 257.277671149549, Press = -312.078049435732 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.133 | 4.133 | 4.133 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 5000 -17677.667 -17677.667 -17744.173 -17744.173 257.38375 257.38375 31196.896 31196.896 -2588.813 -2588.813 6000 -17680.851 -17680.851 -17745.347 -17745.347 249.60687 249.60687 31183.958 31183.958 -1073.6014 -1073.6014 Loop time of 12.36 on 1 procs for 1000 steps with 2000 atoms Performance: 6.990 ns/day, 3.433 hours/ns, 80.906 timesteps/s 42.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 | 12.066 | 12.066 | 12.066 | 0.0 | 97.62 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.04488 | 0.04488 | 0.04488 | 0.0 | 0.36 Output | 2.8849e-05 | 2.8849e-05 | 2.8849e-05 | 0.0 | 0.00 Modify | 0.23771 | 0.23771 | 0.23771 | 0.0 | 1.92 Other | | 0.01107 | | | 0.09 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 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 = 251.660254793671, Press = -26.8195669945299 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.133 | 4.133 | 4.133 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 6000 -17680.851 -17680.851 -17745.347 -17745.347 249.60687 249.60687 31183.958 31183.958 -1073.6014 -1073.6014 7000 -17679.313 -17679.313 -17743.999 -17743.999 250.34186 250.34186 31156.949 31156.949 2112.4664 2112.4664 Loop time of 11.8743 on 1 procs for 1000 steps with 2000 atoms Performance: 7.276 ns/day, 3.298 hours/ns, 84.215 timesteps/s 44.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 | 11.647 | 11.647 | 11.647 | 0.0 | 98.08 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.025234 | 0.025234 | 0.025234 | 0.0 | 0.21 Output | 2.8849e-05 | 2.8849e-05 | 2.8849e-05 | 0.0 | 0.00 Modify | 0.17085 | 0.17085 | 0.17085 | 0.0 | 1.44 Other | | 0.03135 | | | 0.26 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 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.395161936293, Press = -19.7507960237535 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.133 | 4.133 | 4.133 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 7000 -17679.313 -17679.313 -17743.999 -17743.999 250.34186 250.34186 31156.949 31156.949 2112.4664 2112.4664 8000 -17677.645 -17677.645 -17744.311 -17744.311 258.00737 258.00737 31176.019 31176.019 -103.31808 -103.31808 Loop time of 11.678 on 1 procs for 1000 steps with 2000 atoms Performance: 7.399 ns/day, 3.244 hours/ns, 85.631 timesteps/s 44.5% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 11.43 | 11.43 | 11.43 | 0.0 | 97.88 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.025774 | 0.025774 | 0.025774 | 0.0 | 0.22 Output | 2.3127e-05 | 2.3127e-05 | 2.3127e-05 | 0.0 | 0.00 Modify | 0.21062 | 0.21062 | 0.21062 | 0.0 | 1.80 Other | | 0.01113 | | | 0.10 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 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.149081490064, Press = -8.02108187310298 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.133 | 4.133 | 4.133 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 8000 -17677.645 -17677.645 -17744.311 -17744.311 258.00737 258.00737 31176.019 31176.019 -103.31808 -103.31808 9000 -17679.755 -17679.755 -17744.868 -17744.868 251.9971 251.9971 31195.798 31195.798 -2658.9304 -2658.9304 Loop time of 11.5739 on 1 procs for 1000 steps with 2000 atoms Performance: 7.465 ns/day, 3.215 hours/ns, 86.402 timesteps/s 44.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 | 11.387 | 11.387 | 11.387 | 0.0 | 98.39 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.034929 | 0.034929 | 0.034929 | 0.0 | 0.30 Output | 2.6941e-05 | 2.6941e-05 | 2.6941e-05 | 0.0 | 0.00 Modify | 0.14042 | 0.14042 | 0.14042 | 0.0 | 1.21 Other | | 0.01118 | | | 0.10 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 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.664910330385, Press = -13.8926179920046 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.133 | 4.133 | 4.133 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 9000 -17679.755 -17679.755 -17744.868 -17744.868 251.9971 251.9971 31195.798 31195.798 -2658.9304 -2658.9304 10000 -17680.631 -17680.631 -17747.173 -17747.173 257.52367 257.52367 31138.56 31138.56 3819.7739 3819.7739 Loop time of 12.3751 on 1 procs for 1000 steps with 2000 atoms Performance: 6.982 ns/day, 3.438 hours/ns, 80.807 timesteps/s 41.7% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 12.07 | 12.07 | 12.07 | 0.0 | 97.54 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.044793 | 0.044793 | 0.044793 | 0.0 | 0.36 Output | 3.7909e-05 | 3.7909e-05 | 3.7909e-05 | 0.0 | 0.00 Modify | 0.22888 | 0.22888 | 0.22888 | 0.0 | 1.85 Other | | 0.03119 | | | 0.25 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 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.752911372625, Press = -6.82861694602702 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.133 | 4.133 | 4.133 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 10000 -17680.631 -17680.631 -17747.173 -17747.173 257.52367 257.52367 31138.56 31138.56 3819.7739 3819.7739 11000 -17679.325 -17679.325 -17744.441 -17744.441 252.00535 252.00535 31162.72 31162.72 2225.4319 2225.4319 Loop time of 11.9787 on 1 procs for 1000 steps with 2000 atoms Performance: 7.213 ns/day, 3.327 hours/ns, 83.481 timesteps/s 43.2% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 11.793 | 11.793 | 11.793 | 0.0 | 98.45 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.024683 | 0.024683 | 0.024683 | 0.0 | 0.21 Output | 2.2173e-05 | 2.2173e-05 | 2.2173e-05 | 0.0 | 0.00 Modify | 0.14969 | 0.14969 | 0.14969 | 0.0 | 1.25 Other | | 0.01118 | | | 0.09 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 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.773446575409, Press = -4.21008663493705 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.133 | 4.133 | 4.133 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 11000 -17679.325 -17679.325 -17744.441 -17744.441 252.00535 252.00535 31162.72 31162.72 2225.4319 2225.4319 12000 -17681.077 -17681.077 -17744.685 -17744.685 246.17081 246.17081 31159.779 31159.779 1860.9827 1860.9827 Loop time of 12.31 on 1 procs for 1000 steps with 2000 atoms Performance: 7.019 ns/day, 3.419 hours/ns, 81.235 timesteps/s 42.4% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 12.005 | 12.005 | 12.005 | 0.0 | 97.52 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.044937 | 0.044937 | 0.044937 | 0.0 | 0.37 Output | 2.3127e-05 | 2.3127e-05 | 2.3127e-05 | 0.0 | 0.00 Modify | 0.22882 | 0.22882 | 0.22882 | 0.0 | 1.86 Other | | 0.03114 | | | 0.25 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 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.906816458156, Press = 3.07117006434031 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.133 | 4.133 | 4.133 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 12000 -17681.077 -17681.077 -17744.685 -17744.685 246.17081 246.17081 31159.779 31159.779 1860.9827 1860.9827 13000 -17677.527 -17677.527 -17744.217 -17744.217 258.09821 258.09821 31203.371 31203.371 -3078.9199 -3078.9199 Loop time of 11.7042 on 1 procs for 1000 steps with 2000 atoms Performance: 7.382 ns/day, 3.251 hours/ns, 85.439 timesteps/s 44.4% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 11.518 | 11.518 | 11.518 | 0.0 | 98.41 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.044954 | 0.044954 | 0.044954 | 0.0 | 0.38 Output | 2.9087e-05 | 2.9087e-05 | 2.9087e-05 | 0.0 | 0.00 Modify | 0.12985 | 0.12985 | 0.12985 | 0.0 | 1.11 Other | | 0.01126 | | | 0.10 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 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.278324963716, Press = -4.90276521144102 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.133 | 4.133 | 4.133 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 13000 -17677.527 -17677.527 -17744.217 -17744.217 258.09821 258.09821 31203.371 31203.371 -3078.9199 -3078.9199 14000 -17680.491 -17680.491 -17745.176 -17745.176 250.33794 250.33794 31174.451 31174.451 -130.83037 -130.83037 Loop time of 11.9862 on 1 procs for 1000 steps with 2000 atoms Performance: 7.208 ns/day, 3.330 hours/ns, 83.429 timesteps/s 43.4% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 11.721 | 11.721 | 11.721 | 0.0 | 97.79 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.045039 | 0.045039 | 0.045039 | 0.0 | 0.38 Output | 2.0981e-05 | 2.0981e-05 | 2.0981e-05 | 0.0 | 0.00 Modify | 0.18912 | 0.18912 | 0.18912 | 0.0 | 1.58 Other | | 0.03117 | | | 0.26 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 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.105543714668, Press = -1.16735104507695 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.133 | 4.133 | 4.133 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 14000 -17680.491 -17680.491 -17745.176 -17745.176 250.33794 250.33794 31174.451 31174.451 -130.83037 -130.83037 15000 -17680.694 -17680.694 -17745.293 -17745.293 250.00174 250.00174 31165.87 31165.87 1109.0148 1109.0148 Loop time of 11.8661 on 1 procs for 1000 steps with 2000 atoms Performance: 7.281 ns/day, 3.296 hours/ns, 84.273 timesteps/s 43.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 | 11.653 | 11.653 | 11.653 | 0.0 | 98.21 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.060775 | 0.060775 | 0.060775 | 0.0 | 0.51 Output | 2.9087e-05 | 2.9087e-05 | 2.9087e-05 | 0.0 | 0.00 Modify | 0.12095 | 0.12095 | 0.12095 | 0.0 | 1.02 Other | | 0.03124 | | | 0.26 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 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.012937223243, Press = -2.72536280687294 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.133 | 4.133 | 4.133 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 15000 -17680.694 -17680.694 -17745.293 -17745.293 250.00174 250.00174 31165.87 31165.87 1109.0148 1109.0148 16000 -17677.73 -17677.73 -17743.026 -17743.026 252.70213 252.70213 31176.137 31176.137 743.13051 743.13051 Loop time of 11.3756 on 1 procs for 1000 steps with 2000 atoms Performance: 7.595 ns/day, 3.160 hours/ns, 87.907 timesteps/s 45.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 | 11.079 | 11.079 | 11.079 | 0.0 | 97.39 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.065083 | 0.065083 | 0.065083 | 0.0 | 0.57 Output | 2.2173e-05 | 2.2173e-05 | 2.2173e-05 | 0.0 | 0.00 Modify | 0.2201 | 0.2201 | 0.2201 | 0.0 | 1.93 Other | | 0.01154 | | | 0.10 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 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.824888620414, Press = -6.12918388067518 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.133 | 4.133 | 4.133 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 16000 -17677.73 -17677.73 -17743.026 -17743.026 252.70213 252.70213 31176.137 31176.137 743.13051 743.13051 17000 -17676.644 -17676.644 -17744.827 -17744.827 263.8758 263.8758 31164.318 31164.318 1911.7794 1911.7794 Loop time of 12.0642 on 1 procs for 1000 steps with 2000 atoms Performance: 7.162 ns/day, 3.351 hours/ns, 82.890 timesteps/s 42.7% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 11.78 | 11.78 | 11.78 | 0.0 | 97.64 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.084876 | 0.084876 | 0.084876 | 0.0 | 0.70 Output | 2.1935e-05 | 2.1935e-05 | 2.1935e-05 | 0.0 | 0.00 Modify | 0.16849 | 0.16849 | 0.16849 | 0.0 | 1.40 Other | | 0.03111 | | | 0.26 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 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.879665099743, Press = -7.28790435380808 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.133 | 4.133 | 4.133 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 17000 -17676.644 -17676.644 -17744.827 -17744.827 263.8758 263.8758 31164.318 31164.318 1911.7794 1911.7794 18000 -17679.519 -17679.519 -17746.235 -17746.235 258.1953 258.1953 31183.208 31183.208 -1508.4954 -1508.4954 Loop time of 11.6529 on 1 procs for 1000 steps with 2000 atoms Performance: 7.414 ns/day, 3.237 hours/ns, 85.815 timesteps/s 44.5% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 11.366 | 11.366 | 11.366 | 0.0 | 97.53 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.065129 | 0.065129 | 0.065129 | 0.0 | 0.56 Output | 2.8133e-05 | 2.8133e-05 | 2.8133e-05 | 0.0 | 0.00 Modify | 0.20884 | 0.20884 | 0.20884 | 0.0 | 1.79 Other | | 0.01332 | | | 0.11 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 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.110731197825, Press = -1.03461487274497 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.133 | 4.133 | 4.133 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 18000 -17679.519 -17679.519 -17746.235 -17746.235 258.1953 258.1953 31183.208 31183.208 -1508.4954 -1508.4954 19000 -17678.567 -17678.567 -17744.316 -17744.316 254.4543 254.4543 31177.121 31177.121 404.12812 404.12812 Loop time of 11.5567 on 1 procs for 1000 steps with 2000 atoms Performance: 7.476 ns/day, 3.210 hours/ns, 86.530 timesteps/s 44.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 | 11.279 | 11.279 | 11.279 | 0.0 | 97.60 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.065021 | 0.065021 | 0.065021 | 0.0 | 0.56 Output | 4.3869e-05 | 4.3869e-05 | 4.3869e-05 | 0.0 | 0.00 Modify | 0.20156 | 0.20156 | 0.20156 | 0.0 | 1.74 Other | | 0.01107 | | | 0.10 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 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.168512926006, Press = -0.37619753583247 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.133 | 4.133 | 4.133 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 19000 -17678.567 -17678.567 -17744.316 -17744.316 254.4543 254.4543 31177.121 31177.121 404.12812 404.12812 20000 -17679.401 -17679.401 -17746.426 -17746.426 259.39283 259.39283 31158.476 31158.476 1149.1367 1149.1367 Loop time of 11.7924 on 1 procs for 1000 steps with 2000 atoms Performance: 7.327 ns/day, 3.276 hours/ns, 84.800 timesteps/s 44.2% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 11.527 | 11.527 | 11.527 | 0.0 | 97.75 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.045077 | 0.045077 | 0.045077 | 0.0 | 0.38 Output | 2.9087e-05 | 2.9087e-05 | 2.9087e-05 | 0.0 | 0.00 Modify | 0.20881 | 0.20881 | 0.20881 | 0.0 | 1.77 Other | | 0.01116 | | | 0.09 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 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.202626930052, Press = 0.457608604530193 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.133 | 4.133 | 4.133 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 20000 -17679.401 -17679.401 -17746.426 -17746.426 259.39283 259.39283 31158.476 31158.476 1149.1367 1149.1367 21000 -17673.281 -17673.281 -17741.702 -17741.702 264.79709 264.79709 31178.06 31178.06 629.97624 629.97624 Loop time of 11.3387 on 1 procs for 1000 steps with 2000 atoms Performance: 7.620 ns/day, 3.150 hours/ns, 88.194 timesteps/s 45.5% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 11.134 | 11.134 | 11.134 | 0.0 | 98.19 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.064707 | 0.064707 | 0.064707 | 0.0 | 0.57 Output | 2.9087e-05 | 2.9087e-05 | 2.9087e-05 | 0.0 | 0.00 Modify | 0.12881 | 0.12881 | 0.12881 | 0.0 | 1.14 Other | | 0.01112 | | | 0.10 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 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.28849534421, Press = -3.42831642794751 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.133 | 4.133 | 4.133 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 21000 -17673.281 -17673.281 -17741.702 -17741.702 264.79709 264.79709 31178.06 31178.06 629.97624 629.97624 22000 -17680.73 -17680.73 -17744.769 -17744.769 247.83689 247.83689 31176.889 31176.889 58.758317 58.758317 Loop time of 11.3624 on 1 procs for 1000 steps with 2000 atoms Performance: 7.604 ns/day, 3.156 hours/ns, 88.010 timesteps/s 45.5% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 11.007 | 11.007 | 11.007 | 0.0 | 96.87 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.064868 | 0.064868 | 0.064868 | 0.0 | 0.57 Output | 3.0041e-05 | 3.0041e-05 | 3.0041e-05 | 0.0 | 0.00 Modify | 0.25931 | 0.25931 | 0.25931 | 0.0 | 2.28 Other | | 0.03112 | | | 0.27 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 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.398981518942, Press = -0.416455052922611 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.133 | 4.133 | 4.133 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 22000 -17680.73 -17680.73 -17744.769 -17744.769 247.83689 247.83689 31176.889 31176.889 58.758317 58.758317 23000 -17682.525 -17682.525 -17746.111 -17746.111 246.08518 246.08518 31180.928 31180.928 -1668.5318 -1668.5318 Loop time of 12.0362 on 1 procs for 1000 steps with 2000 atoms Performance: 7.178 ns/day, 3.343 hours/ns, 83.083 timesteps/s 43.5% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 11.75 | 11.75 | 11.75 | 0.0 | 97.62 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.065248 | 0.065248 | 0.065248 | 0.0 | 0.54 Output | 2.8849e-05 | 2.8849e-05 | 2.8849e-05 | 0.0 | 0.00 Modify | 0.18996 | 0.18996 | 0.18996 | 0.0 | 1.58 Other | | 0.03143 | | | 0.26 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 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.357650774909, Press = -1.28412108994626 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.133 | 4.133 | 4.133 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 23000 -17682.525 -17682.525 -17746.111 -17746.111 246.08518 246.08518 31180.928 31180.928 -1668.5318 -1668.5318 24000 -17678.811 -17678.811 -17743.663 -17743.663 250.98528 250.98528 31184.966 31184.966 -1656.4571 -1656.4571 Loop time of 11.7636 on 1 procs for 1000 steps with 2000 atoms Performance: 7.345 ns/day, 3.268 hours/ns, 85.008 timesteps/s 44.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 | 11.599 | 11.599 | 11.599 | 0.0 | 98.60 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.024971 | 0.024971 | 0.024971 | 0.0 | 0.21 Output | 2.1935e-05 | 2.1935e-05 | 2.1935e-05 | 0.0 | 0.00 Modify | 0.088785 | 0.088785 | 0.088785 | 0.0 | 0.75 Other | | 0.05109 | | | 0.43 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 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.257587299843, Press = -4.00709102455134 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.133 | 4.133 | 4.133 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 24000 -17678.811 -17678.811 -17743.663 -17743.663 250.98528 250.98528 31184.966 31184.966 -1656.4571 -1656.4571 25000 -17679.324 -17679.324 -17744.32 -17744.32 251.54197 251.54197 31187.669 31187.669 -1612.9427 -1612.9427 Loop time of 11.0834 on 1 procs for 1000 steps with 2000 atoms Performance: 7.795 ns/day, 3.079 hours/ns, 90.225 timesteps/s 47.3% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 10.798 | 10.798 | 10.798 | 0.0 | 97.42 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.02519 | 0.02519 | 0.02519 | 0.0 | 0.23 Output | 2.7895e-05 | 2.7895e-05 | 2.7895e-05 | 0.0 | 0.00 Modify | 0.20913 | 0.20913 | 0.20913 | 0.0 | 1.89 Other | | 0.05127 | | | 0.46 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 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.15021736731, Press = -0.449037735475256 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.133 | 4.133 | 4.133 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 25000 -17679.324 -17679.324 -17744.32 -17744.32 251.54197 251.54197 31187.669 31187.669 -1612.9427 -1612.9427 26000 -17680.203 -17680.203 -17745.574 -17745.574 252.99259 252.99259 31156.255 31156.255 2020.2823 2020.2823 Loop time of 12.1392 on 1 procs for 1000 steps with 2000 atoms Performance: 7.117 ns/day, 3.372 hours/ns, 82.378 timesteps/s 43.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 | 11.793 | 11.793 | 11.793 | 0.0 | 97.14 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.045456 | 0.045456 | 0.045456 | 0.0 | 0.37 Output | 2.9087e-05 | 2.9087e-05 | 2.9087e-05 | 0.0 | 0.00 Modify | 0.24995 | 0.24995 | 0.24995 | 0.0 | 2.06 Other | | 0.05117 | | | 0.42 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 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.168556522555, Press = -1.72705560748649 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.133 | 4.133 | 4.133 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 26000 -17680.203 -17680.203 -17745.574 -17745.574 252.99259 252.99259 31156.255 31156.255 2020.2823 2020.2823 27000 -17678.949 -17678.949 -17745.129 -17745.129 256.12481 256.12481 31182.68 31182.68 -938.06403 -938.06403 Loop time of 10.6023 on 1 procs for 1000 steps with 2000 atoms Performance: 8.149 ns/day, 2.945 hours/ns, 94.320 timesteps/s 48.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 | 10.377 | 10.377 | 10.377 | 0.0 | 97.88 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.044867 | 0.044867 | 0.044867 | 0.0 | 0.42 Output | 4.6968e-05 | 4.6968e-05 | 4.6968e-05 | 0.0 | 0.00 Modify | 0.169 | 0.169 | 0.169 | 0.0 | 1.59 Other | | 0.0113 | | | 0.11 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 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.232179078592, Press = -0.543763373434322 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.133 | 4.133 | 4.133 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 27000 -17678.949 -17678.949 -17745.129 -17745.129 256.12481 256.12481 31182.68 31182.68 -938.06403 -938.06403 28000 -17678.343 -17678.343 -17745.883 -17745.883 261.38546 261.38546 31165.533 31165.533 1147.3758 1147.3758 Loop time of 11.3533 on 1 procs for 1000 steps with 2000 atoms Performance: 7.610 ns/day, 3.154 hours/ns, 88.080 timesteps/s 45.5% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 11.088 | 11.088 | 11.088 | 0.0 | 97.67 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.085158 | 0.085158 | 0.085158 | 0.0 | 0.75 Output | 2.2888e-05 | 2.2888e-05 | 2.2888e-05 | 0.0 | 0.00 Modify | 0.16869 | 0.16869 | 0.16869 | 0.0 | 1.49 Other | | 0.01108 | | | 0.10 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 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.22217763859, Press = -2.51658506395311 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.133 | 4.133 | 4.133 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 28000 -17678.343 -17678.343 -17745.883 -17745.883 261.38546 261.38546 31165.533 31165.533 1147.3758 1147.3758 29000 -17679.839 -17679.839 -17745.501 -17745.501 254.12086 254.12086 31177.884 31177.884 -1567.176 -1567.176 Loop time of 10.247 on 1 procs for 1000 steps with 2000 atoms Performance: 8.432 ns/day, 2.846 hours/ns, 97.590 timesteps/s 50.6% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 9.9318 | 9.9318 | 9.9318 | 0.0 | 96.92 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.11497 | 0.11497 | 0.11497 | 0.0 | 1.12 Output | 2.6941e-05 | 2.6941e-05 | 2.6941e-05 | 0.0 | 0.00 Modify | 0.14901 | 0.14901 | 0.14901 | 0.0 | 1.45 Other | | 0.05116 | | | 0.50 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 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.228945167802, Press = -2.48753288335844 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.133 | 4.133 | 4.133 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 29000 -17679.839 -17679.839 -17745.501 -17745.501 254.12086 254.12086 31177.884 31177.884 -1567.176 -1567.176 30000 -17682.115 -17682.115 -17746.333 -17746.333 248.53008 248.53008 31147.11 31147.11 3036.0786 3036.0786 Loop time of 10.316 on 1 procs for 1000 steps with 2000 atoms Performance: 8.375 ns/day, 2.866 hours/ns, 96.937 timesteps/s 50.1% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 10.061 | 10.061 | 10.061 | 0.0 | 97.52 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.055267 | 0.055267 | 0.055267 | 0.0 | 0.54 Output | 2.7895e-05 | 2.7895e-05 | 2.7895e-05 | 0.0 | 0.00 Modify | 0.16891 | 0.16891 | 0.16891 | 0.0 | 1.64 Other | | 0.03114 | | | 0.30 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 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 31175.0347462469 print "LAMMPS calculation completed" LAMMPS calculation completed quit 0