# 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.031088918447495*${_u_distance} variable latticeconst_converted equal 3.031088918447495*1 lattice bcc ${latticeconst_converted} lattice bcc 3.03108891844749 Lattice spacing in x,y,z = 3.0310889 3.0310889 3.0310889 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 (30.310889 30.310889 30.310889) 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 (30.310889 30.310889 30.310889) create_atoms CPU = 0.002 seconds variable mass_converted equal 50.9415*${_u_mass} variable mass_converted equal 50.9415*1 kim_interactions V WARNING: 'kim_' has been renamed to 'kim '. Please update your input. kim interactions V #=== BEGIN kim interactions ================================== pair_style kim MEAM_LAMMPS_OhSeolLee_2020_CoV__MO_771146361182_001 pair_coeff * * V #=== END kim interactions ==================================== mass 1 ${mass_converted} mass 1 50.9415 # initial volume variable v equal vol # assign formula variable V0 equal ${v} # evaluate initial value variable V0 equal 27848.1295338222 variable V0_metal equal ${V0}/(${_u_distance}*${_u_distance}*${_u_distance}) variable V0_metal equal 27848.1295338222/(${_u_distance}*${_u_distance}*${_u_distance}) variable V0_metal equal 27848.1295338222/(1*${_u_distance}*${_u_distance}) variable V0_metal equal 27848.1295338222/(1*1*${_u_distance}) variable V0_metal equal 27848.1295338222/(1*1*1) variable V0_metal_times1000 equal ${V0_metal}*1000 variable V0_metal_times1000 equal 27848.1295338222*1000 print "Initial system volume: ${V0_metal} Angstroms^3" Initial system volume: 27848.1295338222 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 313.15*${_u_temperature} variable temp_converted equal 313.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 313.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 313.15 ${temp_converted} ${Tdamp_converted} iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 313.15 313.15 ${Tdamp_converted} iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 313.15 313.15 0.1 iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 313.15 313.15 0.1 iso 0 ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 313.15 313.15 0.1 iso 0 0 ${Pdamp_converted} fix ensemble all npt temp 313.15 313.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 "313.15 - 0.2" variable T_up equal "313.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_771146361182_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.8 ghost atom cutoff = 6.8 binsize = 3.4, bins = 9 9 9 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 0 -10519.085 -10519.085 -10600 -10600 313.15 313.15 27848.13 27848.13 3103.496 3103.496 1000 -10430.487 -10430.487 -10519.103 -10519.103 342.95272 342.95272 28029.245 28029.245 1409.3767 1409.3767 Loop time of 212.19 on 1 procs for 1000 steps with 2000 atoms Performance: 0.407 ns/day, 58.942 hours/ns, 4.713 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 | 211.7 | 211.7 | 211.7 | 0.0 | 99.77 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.096451 | 0.096451 | 0.096451 | 0.0 | 0.05 Output | 0.00019508 | 0.00019508 | 0.00019508 | 0.0 | 0.00 Modify | 0.33355 | 0.33355 | 0.33355 | 0.0 | 0.16 Other | | 0.05784 | | | 0.03 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: 224000.0 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112.00000 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.8 ghost atom cutoff = 6.8 binsize = 3.4, bins = 9 9 9 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 1000 -10430.487 -10430.487 -10519.103 -10519.103 342.95272 342.95272 28029.245 28029.245 1409.3767 1409.3767 2000 -10435.086 -10435.086 -10514.734 -10514.734 308.24802 308.24802 28054.474 28054.474 119.97806 119.97806 Loop time of 202.22 on 1 procs for 1000 steps with 2000 atoms Performance: 0.427 ns/day, 56.172 hours/ns, 4.945 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 | 201.73 | 201.73 | 201.73 | 0.0 | 99.76 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.098199 | 0.098199 | 0.098199 | 0.0 | 0.05 Output | 0.00022279 | 0.00022279 | 0.00022279 | 0.0 | 0.00 Modify | 0.33404 | 0.33404 | 0.33404 | 0.0 | 0.17 Other | | 0.0568 | | | 0.03 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: 198018.0 ave 198018 max 198018 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 198018 Ave neighs/atom = 99.009000 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.8 ghost atom cutoff = 6.8 binsize = 3.4, bins = 9 9 9 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 2000 -10435.086 -10435.086 -10514.734 -10514.734 308.24802 308.24802 28054.474 28054.474 119.97806 119.97806 3000 -10434.002 -10434.002 -10513.77 -10513.77 308.70949 308.70949 28079.402 28079.402 -1327.6428 -1327.6428 Loop time of 201.072 on 1 procs for 1000 steps with 2000 atoms Performance: 0.430 ns/day, 55.853 hours/ns, 4.973 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 | 200.59 | 200.59 | 200.59 | 0.0 | 99.76 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.094804 | 0.094804 | 0.094804 | 0.0 | 0.05 Output | 0.00022307 | 0.00022307 | 0.00022307 | 0.0 | 0.00 Modify | 0.33216 | 0.33216 | 0.33216 | 0.0 | 0.17 Other | | 0.05923 | | | 0.03 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: 197436.0 ave 197436 max 197436 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 197436 Ave neighs/atom = 98.718000 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.8 ghost atom cutoff = 6.8 binsize = 3.4, bins = 9 9 9 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 3000 -10434.002 -10434.002 -10513.77 -10513.77 308.70949 308.70949 28079.402 28079.402 -1327.6428 -1327.6428 4000 -10433.64 -10433.64 -10513.992 -10513.992 310.97131 310.97131 28074.124 28074.124 -686.52154 -686.52154 Loop time of 200.567 on 1 procs for 1000 steps with 2000 atoms Performance: 0.431 ns/day, 55.713 hours/ns, 4.986 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 | 200.08 | 200.08 | 200.08 | 0.0 | 99.76 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.095556 | 0.095556 | 0.095556 | 0.0 | 0.05 Output | 0.00018704 | 0.00018704 | 0.00018704 | 0.0 | 0.00 Modify | 0.33326 | 0.33326 | 0.33326 | 0.0 | 0.17 Other | | 0.05896 | | | 0.03 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: 196584.0 ave 196584 max 196584 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 196584 Ave neighs/atom = 98.292000 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.8 ghost atom cutoff = 6.8 binsize = 3.4, bins = 9 9 9 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 4000 -10433.64 -10433.64 -10513.992 -10513.992 310.97131 310.97131 28074.124 28074.124 -686.52154 -686.52154 5000 -10432.335 -10432.335 -10516.603 -10516.603 326.1258 326.1258 28074.347 28074.347 -781.89775 -781.89775 Loop time of 201.641 on 1 procs for 1000 steps with 2000 atoms Performance: 0.428 ns/day, 56.011 hours/ns, 4.959 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 | 201.15 | 201.15 | 201.15 | 0.0 | 99.76 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.095473 | 0.095473 | 0.095473 | 0.0 | 0.05 Output | 0.000153 | 0.000153 | 0.000153 | 0.0 | 0.00 Modify | 0.33645 | 0.33645 | 0.33645 | 0.0 | 0.17 Other | | 0.05863 | | | 0.03 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: 196744.0 ave 196744 max 196744 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 196744 Ave neighs/atom = 98.372000 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 = 312.26952671485, Press = -211.389961337569 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.8 ghost atom cutoff = 6.8 binsize = 3.4, bins = 9 9 9 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 5000 -10432.335 -10432.335 -10516.603 -10516.603 326.1258 326.1258 28074.347 28074.347 -781.89775 -781.89775 6000 -10432.993 -10432.993 -10514.536 -10514.536 315.57816 315.57816 28065.559 28065.559 -312.45625 -312.45625 Loop time of 201.329 on 1 procs for 1000 steps with 2000 atoms Performance: 0.429 ns/day, 55.925 hours/ns, 4.967 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 | 200.82 | 200.82 | 200.82 | 0.0 | 99.75 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.095037 | 0.095037 | 0.095037 | 0.0 | 0.05 Output | 0.00015332 | 0.00015332 | 0.00015332 | 0.0 | 0.00 Modify | 0.35096 | 0.35096 | 0.35096 | 0.0 | 0.17 Other | | 0.05863 | | | 0.03 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: 197316.0 ave 197316 max 197316 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 197316 Ave neighs/atom = 98.658000 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 = 312.226995363938, Press = -53.6554022871394 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.8 ghost atom cutoff = 6.8 binsize = 3.4, bins = 9 9 9 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 6000 -10432.993 -10432.993 -10514.536 -10514.536 315.57816 315.57816 28065.559 28065.559 -312.45625 -312.45625 7000 -10434.333 -10434.333 -10509.712 -10509.712 291.72295 291.72295 28012.801 28012.801 2520.5813 2520.5813 Loop time of 201.34 on 1 procs for 1000 steps with 2000 atoms Performance: 0.429 ns/day, 55.928 hours/ns, 4.967 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 | 200.83 | 200.83 | 200.83 | 0.0 | 99.75 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.096635 | 0.096635 | 0.096635 | 0.0 | 0.05 Output | 0.00015579 | 0.00015579 | 0.00015579 | 0.0 | 0.00 Modify | 0.35115 | 0.35115 | 0.35115 | 0.0 | 0.17 Other | | 0.0577 | | | 0.03 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: 197090.0 ave 197090 max 197090 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 197090 Ave neighs/atom = 98.545000 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 = 311.824751284512, Press = -31.1038201289018 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.8 ghost atom cutoff = 6.8 binsize = 3.4, bins = 9 9 9 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 7000 -10434.333 -10434.333 -10509.712 -10509.712 291.72295 291.72295 28012.801 28012.801 2520.5813 2520.5813 8000 -10434.762 -10434.762 -10513.019 -10513.019 302.86402 302.86402 28019.053 28019.053 1826.6998 1826.6998 Loop time of 201.314 on 1 procs for 1000 steps with 2000 atoms Performance: 0.429 ns/day, 55.921 hours/ns, 4.967 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 | 200.81 | 200.81 | 200.81 | 0.0 | 99.75 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.096401 | 0.096401 | 0.096401 | 0.0 | 0.05 Output | 0.00019434 | 0.00019434 | 0.00019434 | 0.0 | 0.00 Modify | 0.35087 | 0.35087 | 0.35087 | 0.0 | 0.17 Other | | 0.05805 | | | 0.03 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: 197254.0 ave 197254 max 197254 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 197254 Ave neighs/atom = 98.627000 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 = 312.560214129015, Press = 2.23021368237731 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.8 ghost atom cutoff = 6.8 binsize = 3.4, bins = 9 9 9 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 8000 -10434.762 -10434.762 -10513.019 -10513.019 302.86402 302.86402 28019.053 28019.053 1826.6998 1826.6998 9000 -10432.771 -10432.771 -10512.384 -10512.384 308.11013 308.11013 28041.755 28041.755 1111.7157 1111.7157 Loop time of 201.415 on 1 procs for 1000 steps with 2000 atoms Performance: 0.429 ns/day, 55.949 hours/ns, 4.965 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 | 200.91 | 200.91 | 200.91 | 0.0 | 99.75 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.096913 | 0.096913 | 0.096913 | 0.0 | 0.05 Output | 0.00015365 | 0.00015365 | 0.00015365 | 0.0 | 0.00 Modify | 0.35125 | 0.35125 | 0.35125 | 0.0 | 0.17 Other | | 0.05732 | | | 0.03 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: 197428.0 ave 197428 max 197428 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 197428 Ave neighs/atom = 98.714000 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 = 312.361938616921, Press = 7.26342214638796 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.8 ghost atom cutoff = 6.8 binsize = 3.4, bins = 9 9 9 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 9000 -10432.771 -10432.771 -10512.384 -10512.384 308.11013 308.11013 28041.755 28041.755 1111.7157 1111.7157 10000 -10433.784 -10433.784 -10513.706 -10513.706 309.30696 309.30696 28049.248 28049.248 437.22461 437.22461 Loop time of 201.445 on 1 procs for 1000 steps with 2000 atoms Performance: 0.429 ns/day, 55.957 hours/ns, 4.964 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 | 200.94 | 200.94 | 200.94 | 0.0 | 99.75 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.095993 | 0.095993 | 0.095993 | 0.0 | 0.05 Output | 0.00015304 | 0.00015304 | 0.00015304 | 0.0 | 0.00 Modify | 0.35196 | 0.35196 | 0.35196 | 0.0 | 0.17 Other | | 0.05705 | | | 0.03 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: 197492.0 ave 197492 max 197492 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 197492 Ave neighs/atom = 98.746000 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 = 312.441036872674, Press = 7.43335182600618 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.8 ghost atom cutoff = 6.8 binsize = 3.4, bins = 9 9 9 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 10000 -10433.784 -10433.784 -10513.706 -10513.706 309.30696 309.30696 28049.248 28049.248 437.22461 437.22461 11000 -10435.268 -10435.268 -10514.952 -10514.952 308.38632 308.38632 28066.241 28066.241 -703.28293 -703.28293 Loop time of 201.744 on 1 procs for 1000 steps with 2000 atoms Performance: 0.428 ns/day, 56.040 hours/ns, 4.957 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 | 201.24 | 201.24 | 201.24 | 0.0 | 99.75 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.096084 | 0.096084 | 0.096084 | 0.0 | 0.05 Output | 0.00015533 | 0.00015533 | 0.00015533 | 0.0 | 0.00 Modify | 0.35168 | 0.35168 | 0.35168 | 0.0 | 0.17 Other | | 0.05736 | | | 0.03 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: 197202.0 ave 197202 max 197202 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 197202 Ave neighs/atom = 98.601000 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 = 312.54386592065, Press = 5.95406323500658 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.8 ghost atom cutoff = 6.8 binsize = 3.4, bins = 9 9 9 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 11000 -10435.268 -10435.268 -10514.952 -10514.952 308.38632 308.38632 28066.241 28066.241 -703.28293 -703.28293 12000 -10434.921 -10434.921 -10513.066 -10513.066 302.42823 302.42823 28078.876 28078.876 -1177.3243 -1177.3243 Loop time of 201.562 on 1 procs for 1000 steps with 2000 atoms Performance: 0.429 ns/day, 55.989 hours/ns, 4.961 timesteps/s 100.1% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 201.06 | 201.06 | 201.06 | 0.0 | 99.75 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.095704 | 0.095704 | 0.095704 | 0.0 | 0.05 Output | 0.00015619 | 0.00015619 | 0.00015619 | 0.0 | 0.00 Modify | 0.35146 | 0.35146 | 0.35146 | 0.0 | 0.17 Other | | 0.05759 | | | 0.03 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: 196952.0 ave 196952 max 196952 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 196952 Ave neighs/atom = 98.476000 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 = 312.46529032074, Press = 3.47706682243328 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.8 ghost atom cutoff = 6.8 binsize = 3.4, bins = 9 9 9 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 12000 -10434.921 -10434.921 -10513.066 -10513.066 302.42823 302.42823 28078.876 28078.876 -1177.3243 -1177.3243 13000 -10434.499 -10434.499 -10517.057 -10517.057 319.50564 319.50564 28082.935 28082.935 -1813.8817 -1813.8817 Loop time of 201.491 on 1 procs for 1000 steps with 2000 atoms Performance: 0.429 ns/day, 55.970 hours/ns, 4.963 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 | 200.99 | 200.99 | 200.99 | 0.0 | 99.75 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.096435 | 0.096435 | 0.096435 | 0.0 | 0.05 Output | 0.00015323 | 0.00015323 | 0.00015323 | 0.0 | 0.00 Modify | 0.35086 | 0.35086 | 0.35086 | 0.0 | 0.17 Other | | 0.05804 | | | 0.03 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: 196800.0 ave 196800 max 196800 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 196800 Ave neighs/atom = 98.400000 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 = 312.537025794587, Press = -0.12158930902034 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.8 ghost atom cutoff = 6.8 binsize = 3.4, bins = 9 9 9 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 13000 -10434.499 -10434.499 -10517.057 -10517.057 319.50564 319.50564 28082.935 28082.935 -1813.8817 -1813.8817 14000 -10431.84 -10431.84 -10511.659 -10511.659 308.91094 308.91094 28085.559 28085.559 -1415.8888 -1415.8888 Loop time of 201.662 on 1 procs for 1000 steps with 2000 atoms Performance: 0.428 ns/day, 56.017 hours/ns, 4.959 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 | 201.16 | 201.16 | 201.16 | 0.0 | 99.75 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.096679 | 0.096679 | 0.096679 | 0.0 | 0.05 Output | 0.00015598 | 0.00015598 | 0.00015598 | 0.0 | 0.00 Modify | 0.35185 | 0.35185 | 0.35185 | 0.0 | 0.17 Other | | 0.0575 | | | 0.03 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: 196906.0 ave 196906 max 196906 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 196906 Ave neighs/atom = 98.453000 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 = 312.699876226167, Press = -5.48492568420539 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.8 ghost atom cutoff = 6.8 binsize = 3.4, bins = 9 9 9 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 14000 -10431.84 -10431.84 -10511.659 -10511.659 308.91094 308.91094 28085.559 28085.559 -1415.8888 -1415.8888 15000 -10431.819 -10431.819 -10513.797 -10513.797 317.26662 317.26662 28043.295 28043.295 713.6323 713.6323 Loop time of 201.425 on 1 procs for 1000 steps with 2000 atoms Performance: 0.429 ns/day, 55.951 hours/ns, 4.965 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 | 200.92 | 200.92 | 200.92 | 0.0 | 99.75 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.096387 | 0.096387 | 0.096387 | 0.0 | 0.05 Output | 0.00015379 | 0.00015379 | 0.00015379 | 0.0 | 0.00 Modify | 0.35098 | 0.35098 | 0.35098 | 0.0 | 0.17 Other | | 0.05765 | | | 0.03 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: 196806.0 ave 196806 max 196806 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 196806 Ave neighs/atom = 98.403000 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 = 312.710325834135, Press = -7.12930466406833 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.8 ghost atom cutoff = 6.8 binsize = 3.4, bins = 9 9 9 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 15000 -10431.819 -10431.819 -10513.797 -10513.797 317.26662 317.26662 28043.295 28043.295 713.6323 713.6323 16000 -10432.744 -10432.744 -10515.994 -10515.994 322.18509 322.18509 28021.252 28021.252 1875.4532 1875.4532 Loop time of 201.56 on 1 procs for 1000 steps with 2000 atoms Performance: 0.429 ns/day, 55.989 hours/ns, 4.961 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 | 201.05 | 201.05 | 201.05 | 0.0 | 99.75 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.096382 | 0.096382 | 0.096382 | 0.0 | 0.05 Output | 0.00015479 | 0.00015479 | 0.00015479 | 0.0 | 0.00 Modify | 0.351 | 0.351 | 0.351 | 0.0 | 0.17 Other | | 0.05769 | | | 0.03 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: 197436.0 ave 197436 max 197436 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 197436 Ave neighs/atom = 98.718000 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 = 312.68476510458, Press = -1.92995329718866 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.8 ghost atom cutoff = 6.8 binsize = 3.4, bins = 9 9 9 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 16000 -10432.744 -10432.744 -10515.994 -10515.994 322.18509 322.18509 28021.252 28021.252 1875.4532 1875.4532 17000 -10435.506 -10435.506 -10517.335 -10517.335 316.68658 316.68658 28019.94 28019.94 1706.6938 1706.6938 Loop time of 201.823 on 1 procs for 1000 steps with 2000 atoms Performance: 0.428 ns/day, 56.062 hours/ns, 4.955 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 | 201.32 | 201.32 | 201.32 | 0.0 | 99.75 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.096974 | 0.096974 | 0.096974 | 0.0 | 0.05 Output | 0.00015434 | 0.00015434 | 0.00015434 | 0.0 | 0.00 Modify | 0.35152 | 0.35152 | 0.35152 | 0.0 | 0.17 Other | | 0.05726 | | | 0.03 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: 197992.0 ave 197992 max 197992 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 197992 Ave neighs/atom = 98.996000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 312.783084982376, Press = 1.5848630091128 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.8 ghost atom cutoff = 6.8 binsize = 3.4, bins = 9 9 9 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 17000 -10435.506 -10435.506 -10517.335 -10517.335 316.68658 316.68658 28019.94 28019.94 1706.6938 1706.6938 18000 -10430.961 -10430.961 -10512.933 -10512.933 317.23871 317.23871 28037.197 28037.197 1417.3153 1417.3153 Loop time of 201.368 on 1 procs for 1000 steps with 2000 atoms Performance: 0.429 ns/day, 55.936 hours/ns, 4.966 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 | 200.86 | 200.86 | 200.86 | 0.0 | 99.75 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.095888 | 0.095888 | 0.095888 | 0.0 | 0.05 Output | 0.00015788 | 0.00015788 | 0.00015788 | 0.0 | 0.00 Modify | 0.35013 | 0.35013 | 0.35013 | 0.0 | 0.17 Other | | 0.05927 | | | 0.03 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: 197980.0 ave 197980 max 197980 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 197980 Ave neighs/atom = 98.990000 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 = 312.918971016514, Press = 3.69421932053474 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.8 ghost atom cutoff = 6.8 binsize = 3.4, bins = 9 9 9 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 18000 -10430.961 -10430.961 -10512.933 -10512.933 317.23871 317.23871 28037.197 28037.197 1417.3153 1417.3153 19000 -10433.14 -10433.14 -10515.345 -10515.345 318.14273 318.14273 28059.257 28059.257 -61.873821 -61.873821 Loop time of 201.26 on 1 procs for 1000 steps with 2000 atoms Performance: 0.429 ns/day, 55.906 hours/ns, 4.969 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 | 200.76 | 200.76 | 200.76 | 0.0 | 99.75 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.095363 | 0.095363 | 0.095363 | 0.0 | 0.05 Output | 0.00019356 | 0.00019356 | 0.00019356 | 0.0 | 0.00 Modify | 0.34939 | 0.34939 | 0.34939 | 0.0 | 0.17 Other | | 0.05924 | | | 0.03 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: 197616.0 ave 197616 max 197616 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 197616 Ave neighs/atom = 98.808000 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 = 313.052901856102, Press = 4.0617829955525 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.8 ghost atom cutoff = 6.8 binsize = 3.4, bins = 9 9 9 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 19000 -10433.14 -10433.14 -10515.345 -10515.345 318.14273 318.14273 28059.257 28059.257 -61.873821 -61.873821 20000 -10435.104 -10435.104 -10515.328 -10515.328 310.4757 310.4757 28083.384 28083.384 -1527.5355 -1527.5355 Loop time of 201.39 on 1 procs for 1000 steps with 2000 atoms Performance: 0.429 ns/day, 55.942 hours/ns, 4.965 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 | 200.89 | 200.89 | 200.89 | 0.0 | 99.75 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.094831 | 0.094831 | 0.094831 | 0.0 | 0.05 Output | 0.00019351 | 0.00019351 | 0.00019351 | 0.0 | 0.00 Modify | 0.3496 | 0.3496 | 0.3496 | 0.0 | 0.17 Other | | 0.05955 | | | 0.03 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: 197274.0 ave 197274 max 197274 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 197274 Ave neighs/atom = 98.637000 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 = 313.124507349299, Press = 2.22923381488773 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.8 ghost atom cutoff = 6.8 binsize = 3.4, bins = 9 9 9 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 20000 -10435.104 -10435.104 -10515.328 -10515.328 310.4757 310.4757 28083.384 28083.384 -1527.5355 -1527.5355 21000 -10433.39 -10433.39 -10514.842 -10514.842 315.22563 315.22563 28113.169 28113.169 -3006.9679 -3006.9679 Loop time of 200.354 on 1 procs for 1000 steps with 2000 atoms Performance: 0.431 ns/day, 55.654 hours/ns, 4.991 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 | 199.85 | 199.85 | 199.85 | 0.0 | 99.75 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.094534 | 0.094534 | 0.094534 | 0.0 | 0.05 Output | 0.00015423 | 0.00015423 | 0.00015423 | 0.0 | 0.00 Modify | 0.34931 | 0.34931 | 0.34931 | 0.0 | 0.17 Other | | 0.05912 | | | 0.03 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: 196944.0 ave 196944 max 196944 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 196944 Ave neighs/atom = 98.472000 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 = 313.084320890389, Press = -0.777392814997088 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.8 ghost atom cutoff = 6.8 binsize = 3.4, bins = 9 9 9 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 21000 -10433.39 -10433.39 -10514.842 -10514.842 315.22563 315.22563 28113.169 28113.169 -3006.9679 -3006.9679 22000 -10431.334 -10431.334 -10514.721 -10514.721 322.71785 322.71785 28103.254 28103.254 -2412.3862 -2412.3862 Loop time of 200.306 on 1 procs for 1000 steps with 2000 atoms Performance: 0.431 ns/day, 55.641 hours/ns, 4.992 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 | 199.8 | 199.8 | 199.8 | 0.0 | 99.75 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.097074 | 0.097074 | 0.097074 | 0.0 | 0.05 Output | 0.00015244 | 0.00015244 | 0.00015244 | 0.0 | 0.00 Modify | 0.35298 | 0.35298 | 0.35298 | 0.0 | 0.18 Other | | 0.06016 | | | 0.03 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: 196434.0 ave 196434 max 196434 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 196434 Ave neighs/atom = 98.217000 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 = 313.12648171881, Press = -3.45424863920882 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.8 ghost atom cutoff = 6.8 binsize = 3.4, bins = 9 9 9 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 22000 -10431.334 -10431.334 -10514.721 -10514.721 322.71785 322.71785 28103.254 28103.254 -2412.3862 -2412.3862 23000 -10434.27 -10434.27 -10513.32 -10513.32 305.9302 305.9302 28058.887 28058.887 -3.6688373 -3.6688373 Loop time of 200.432 on 1 procs for 1000 steps with 2000 atoms Performance: 0.431 ns/day, 55.675 hours/ns, 4.989 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 | 199.92 | 199.92 | 199.92 | 0.0 | 99.75 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.096086 | 0.096086 | 0.096086 | 0.0 | 0.05 Output | 0.00015285 | 0.00015285 | 0.00015285 | 0.0 | 0.00 Modify | 0.35156 | 0.35156 | 0.35156 | 0.0 | 0.18 Other | | 0.06045 | | | 0.03 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: 196492.0 ave 196492 max 196492 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 196492 Ave neighs/atom = 98.246000 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 = 313.225874664051, Press = -1.86065479576017 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.8 ghost atom cutoff = 6.8 binsize = 3.4, bins = 9 9 9 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 23000 -10434.27 -10434.27 -10513.32 -10513.32 305.9302 305.9302 28058.887 28058.887 -3.6688373 -3.6688373 24000 -10433.498 -10433.498 -10514.442 -10514.442 313.26374 313.26374 28049.945 28049.945 410.85443 410.85443 Loop time of 200.576 on 1 procs for 1000 steps with 2000 atoms Performance: 0.431 ns/day, 55.715 hours/ns, 4.986 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 | 200.07 | 200.07 | 200.07 | 0.0 | 99.75 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.096721 | 0.096721 | 0.096721 | 0.0 | 0.05 Output | 0.00015454 | 0.00015454 | 0.00015454 | 0.0 | 0.00 Modify | 0.35161 | 0.35161 | 0.35161 | 0.0 | 0.18 Other | | 0.05999 | | | 0.03 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: 197194.0 ave 197194 max 197194 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 197194 Ave neighs/atom = 98.597000 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 = 313.227553952824, Press = -0.355178620565415 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.8 ghost atom cutoff = 6.8 binsize = 3.4, bins = 9 9 9 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 24000 -10433.498 -10433.498 -10514.442 -10514.442 313.26374 313.26374 28049.945 28049.945 410.85443 410.85443 25000 -10436.963 -10436.963 -10517.215 -10517.215 310.58345 310.58345 28038.69 28038.69 928.18218 928.18218 Loop time of 200.965 on 1 procs for 1000 steps with 2000 atoms Performance: 0.430 ns/day, 55.824 hours/ns, 4.976 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 | 200.45 | 200.45 | 200.45 | 0.0 | 99.75 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.096562 | 0.096562 | 0.096562 | 0.0 | 0.05 Output | 0.00016582 | 0.00016582 | 0.00016582 | 0.0 | 0.00 Modify | 0.35405 | 0.35405 | 0.35405 | 0.0 | 0.18 Other | | 0.06 | | | 0.03 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: 197416.0 ave 197416 max 197416 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 197416 Ave neighs/atom = 98.708000 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_T313.15.out" else "print 'not_converged' file output/vol_T313.15.out" print '${V}' file output/vol_T313.15.out 28056.6841807349 print "LAMMPS calculation completed" LAMMPS calculation completed quit 0