# 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.6133156493306164*${_u_distance} variable latticeconst_converted equal 3.6133156493306164*1 lattice fcc ${latticeconst_converted} lattice fcc 3.61331564933062 Lattice spacing in x,y,z = 3.6133156 3.6133156 3.6133156 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 (36.133156 36.133156 36.133156) 1 by 1 by 1 MPI processor grid create_atoms 1 box Created 4000 atoms using lattice units in orthogonal box = (0.0000000 0.0000000 0.0000000) to (36.133156 36.133156 36.133156) create_atoms CPU = 0.005 seconds variable mass_converted equal 63.546*${_u_mass} variable mass_converted equal 63.546*1 kim_interactions Cu WARNING: 'kim_' has been renamed to 'kim '. Please update your input. kim interactions Cu #=== BEGIN kim interactions ================================== pair_style kim MEAM_LAMMPS_KimLee_2008_CuZr__MO_407917731909_000 pair_coeff * * Cu #=== END kim interactions ==================================== mass 1 ${mass_converted} mass 1 63.546 # initial volume variable v equal vol # assign formula variable V0 equal ${v} # evaluate initial value variable V0 equal 47175.6297173106 variable V0_metal equal ${V0}/(${_u_distance}*${_u_distance}*${_u_distance}) variable V0_metal equal 47175.6297173106/(${_u_distance}*${_u_distance}*${_u_distance}) variable V0_metal equal 47175.6297173106/(1*${_u_distance}*${_u_distance}) variable V0_metal equal 47175.6297173106/(1*1*${_u_distance}) variable V0_metal equal 47175.6297173106/(1*1*1) variable V0_metal_times1000 equal ${V0_metal}*1000 variable V0_metal_times1000 equal 47175.6297173106*1000 print "Initial system volume: ${V0_metal} Angstroms^3" Initial system volume: 47175.6297173106 Angstroms^3 # set the time step to 0.001 picoseconds variable timestep_converted equal 0.001*${_u_time} variable timestep_converted equal 0.001*1 timestep ${timestep_converted} timestep 0.001 variable temp_converted equal 293.15*${_u_temperature} variable temp_converted equal 293.15*1 variable Tdamp_converted equal 0.1*${_u_time} variable Tdamp_converted equal 0.1*1 variable press_converted equal 0.0*${_u_pressure} variable press_converted equal 0.0*1 variable Pdamp_converted equal 1*${_u_time} variable Pdamp_converted equal 1*1 # create initial velocities consistent with the chosen temperature velocity all create ${temp_converted} 17 mom yes rot yes velocity all create 293.15 17 mom yes rot yes # set NPT ensemble for all atoms fix ensemble all npt temp ${temp_converted} ${temp_converted} ${Tdamp_converted} iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 293.15 ${temp_converted} ${Tdamp_converted} iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 293.15 293.15 ${Tdamp_converted} iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 293.15 293.15 0.1 iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 293.15 293.15 0.1 iso 0 ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 293.15 293.15 0.1 iso 0 0 ${Pdamp_converted} fix ensemble all npt temp 293.15 293.15 0.1 iso 0 0 1 # compute the time averages of pressure, temperature, and volume, respectively # ignore the first 5000 timesteps variable etotal_metal equal etotal/${_u_energy} variable etotal_metal equal etotal/1 variable pe_metal equal pe/${_u_energy} variable pe_metal equal pe/1 variable T_metal equal temp/${_u_temperature} variable T_metal equal temp/1 variable V_metal equal vol/(${_u_distance}*${_u_distance}*${_u_distance}) variable V_metal equal vol/(1*${_u_distance}*${_u_distance}) variable V_metal equal vol/(1*1*${_u_distance}) variable V_metal equal vol/(1*1*1) variable P_metal equal press/${_u_pressure} variable P_metal equal press/1 fix avgmyTemp all ave/time 5 20 100 v_T_metal ave running start 5000 fix avgmyPress all ave/time 5 20 100 v_P_metal ave running start 5000 fix avgmyVol all ave/time 5 20 100 v_V_metal ave running start 5000 # extract fix quantities into variables so they can be used in if-else logic later. variable T equal f_avgmyTemp variable P equal f_avgmyPress variable V equal f_avgmyVol # set error bounds for temperature and pressure in original metal units (K and bar) variable T_low equal "293.15 - 0.2" variable T_up equal "293.15 + 0.2" variable P_low equal "0.0 - 0.2" variable P_up equal "0.0 + 0.2" # print to logfile every 1000 timesteps thermo_style custom step etotal v_etotal_metal pe v_pe_metal temp v_T_metal vol v_V_metal press v_P_metal thermo 1000 # Run a simulation for at most 2000*1000 timesteps. At each 1000th time step, check # whether the temperature and pressure have converged. If yes, break. label top variable a loop 2000 run 1000 CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE Your simulation uses code contributions which should be cited: - OpenKIM: https://doi.org/10.1007/s11837-011-0102-6 @Article{tadmor:elliott:2011, author = {E. B. Tadmor and R. S. Elliott and J. P. Sethna and R. E. Miller and C. A. Becker}, title = {The potential of atomistic simulations and the {K}nowledgebase of {I}nteratomic {M}odels}, journal = {{JOM}}, year = 2011, volume = 63, number = 17, pages = {17}, doi = {10.1007/s11837-011-0102-6} } - OpenKIM potential: https://openkim.org/cite/MO_407917731909_000#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 = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 5.370 | 5.370 | 5.370 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 0 -14008.467 -14008.467 -14160 -14160 293.15 293.15 47175.63 47175.63 3430.8919 3430.8919 1000 -13849.002 -13849.002 -14003.684 -14003.684 299.24183 299.24183 48019.328 48019.328 -2367.7404 -2367.7404 Loop time of 392.935 on 1 procs for 1000 steps with 4000 atoms Performance: 0.220 ns/day, 109.149 hours/ns, 2.545 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 | 391.94 | 391.94 | 391.94 | 0.0 | 99.75 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.15914 | 0.15914 | 0.15914 | 0.0 | 0.04 Output | 0.00022793 | 0.00022793 | 0.00022793 | 0.0 | 0.00 Modify | 0.73739 | 0.73739 | 0.73739 | 0.0 | 0.19 Other | | 0.1025 | | | 0.03 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 536000.0 ave 536000 max 536000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 536000 Ave neighs/atom = 134.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 = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 5.370 | 5.370 | 5.370 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 1000 -13849.002 -13849.002 -14003.684 -14003.684 299.24183 299.24183 48019.328 48019.328 -2367.7404 -2367.7404 2000 -13859.305 -13859.305 -14012.236 -14012.236 295.85458 295.85458 47901.143 47901.143 142.62974 142.62974 Loop time of 418.268 on 1 procs for 1000 steps with 4000 atoms Performance: 0.207 ns/day, 116.186 hours/ns, 2.391 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 | 417.2 | 417.2 | 417.2 | 0.0 | 99.75 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16888 | 0.16888 | 0.16888 | 0.0 | 0.04 Output | 0.00027248 | 0.00027248 | 0.00027248 | 0.0 | 0.00 Modify | 0.78933 | 0.78933 | 0.78933 | 0.0 | 0.19 Other | | 0.1063 | | | 0.03 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5843.00 ave 5843 max 5843 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: 529250.0 ave 529250 max 529250 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 529250 Ave neighs/atom = 132.31250 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 = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 5.370 | 5.370 | 5.370 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 2000 -13859.305 -13859.305 -14012.236 -14012.236 295.85458 295.85458 47901.143 47901.143 142.62974 142.62974 3000 -13856.435 -13856.435 -14011.373 -14011.373 299.73786 299.73786 47942.165 47942.165 -746.16536 -746.16536 Loop time of 419.508 on 1 procs for 1000 steps with 4000 atoms Performance: 0.206 ns/day, 116.530 hours/ns, 2.384 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 | 418.44 | 418.44 | 418.44 | 0.0 | 99.74 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16846 | 0.16846 | 0.16846 | 0.0 | 0.04 Output | 0.00021278 | 0.00021278 | 0.00021278 | 0.0 | 0.00 Modify | 0.79538 | 0.79538 | 0.79538 | 0.0 | 0.19 Other | | 0.106 | | | 0.03 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 530996.0 ave 530996 max 530996 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 530996 Ave neighs/atom = 132.74900 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 = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 5.370 | 5.370 | 5.370 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 3000 -13856.435 -13856.435 -14011.373 -14011.373 299.73786 299.73786 47942.165 47942.165 -746.16536 -746.16536 4000 -13856.343 -13856.343 -14007.995 -14007.995 293.38063 293.38063 47939.837 47939.837 -455.49599 -455.49599 Loop time of 369.556 on 1 procs for 1000 steps with 4000 atoms Performance: 0.234 ns/day, 102.655 hours/ns, 2.706 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 | 368.66 | 368.66 | 368.66 | 0.0 | 99.76 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.14905 | 0.14905 | 0.14905 | 0.0 | 0.04 Output | 0.00025891 | 0.00025891 | 0.00025891 | 0.0 | 0.00 Modify | 0.65288 | 0.65288 | 0.65288 | 0.0 | 0.18 Other | | 0.09773 | | | 0.03 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5843.00 ave 5843 max 5843 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: 530960.0 ave 530960 max 530960 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 530960 Ave neighs/atom = 132.74000 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 = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 5.370 | 5.370 | 5.370 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 4000 -13856.343 -13856.343 -14007.995 -14007.995 293.38063 293.38063 47939.837 47939.837 -455.49599 -455.49599 5000 -13859.978 -13859.978 -14007.154 -14007.154 284.72245 284.72245 47898.374 47898.374 826.36639 826.36639 Loop time of 381.602 on 1 procs for 1000 steps with 4000 atoms Performance: 0.226 ns/day, 106.001 hours/ns, 2.621 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 | 380.66 | 380.66 | 380.66 | 0.0 | 99.75 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.15359 | 0.15359 | 0.15359 | 0.0 | 0.04 Output | 0.00024826 | 0.00024826 | 0.00024826 | 0.0 | 0.00 Modify | 0.68974 | 0.68974 | 0.68974 | 0.0 | 0.18 Other | | 0.09868 | | | 0.03 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5843.00 ave 5843 max 5843 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: 530340.0 ave 530340 max 530340 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 530340 Ave neighs/atom = 132.58500 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 293.810508600348, Press = 334.949140754789 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 = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 5.370 | 5.370 | 5.370 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 5000 -13859.978 -13859.978 -14007.154 -14007.154 284.72245 284.72245 47898.374 47898.374 826.36639 826.36639 6000 -13853.348 -13853.348 -14008.924 -14008.924 300.97202 300.97202 47960.439 47960.439 -939.9283 -939.9283 Loop time of 431.317 on 1 procs for 1000 steps with 4000 atoms Performance: 0.200 ns/day, 119.810 hours/ns, 2.318 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 | 430.2 | 430.2 | 430.2 | 0.0 | 99.74 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16688 | 0.16688 | 0.16688 | 0.0 | 0.04 Output | 0.00029623 | 0.00029623 | 0.00029623 | 0.0 | 0.00 Modify | 0.84032 | 0.84032 | 0.84032 | 0.0 | 0.19 Other | | 0.1094 | | | 0.03 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5847.00 ave 5847 max 5847 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: 530800.0 ave 530800 max 530800 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 530800 Ave neighs/atom = 132.70000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 292.542087612406, Press = -15.9123680764866 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 = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 5.370 | 5.370 | 5.370 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 6000 -13853.348 -13853.348 -14008.924 -14008.924 300.97202 300.97202 47960.439 47960.439 -939.9283 -939.9283 7000 -13859.988 -13859.988 -14009.647 -14009.647 289.52613 289.52613 47858.292 47858.292 1789.3869 1789.3869 Loop time of 431.859 on 1 procs for 1000 steps with 4000 atoms Performance: 0.200 ns/day, 119.961 hours/ns, 2.316 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 | 430.74 | 430.74 | 430.74 | 0.0 | 99.74 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16931 | 0.16931 | 0.16931 | 0.0 | 0.04 Output | 0.00027468 | 0.00027468 | 0.00027468 | 0.0 | 0.00 Modify | 0.8441 | 0.8441 | 0.8441 | 0.0 | 0.20 Other | | 0.1101 | | | 0.03 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 530778.0 ave 530778 max 530778 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 530778 Ave neighs/atom = 132.69450 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 292.973491953525, Press = 30.9873148035644 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 = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 5.370 | 5.370 | 5.370 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 7000 -13859.988 -13859.988 -14009.647 -14009.647 289.52613 289.52613 47858.292 47858.292 1789.3869 1789.3869 8000 -13858.608 -13858.608 -14007.995 -14007.995 289.0005 289.0005 47971.786 47971.786 -1397.5379 -1397.5379 Loop time of 421.062 on 1 procs for 1000 steps with 4000 atoms Performance: 0.205 ns/day, 116.962 hours/ns, 2.375 timesteps/s 99.8% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 419.95 | 419.95 | 419.95 | 0.0 | 99.73 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16859 | 0.16859 | 0.16859 | 0.0 | 0.04 Output | 0.00028488 | 0.00028488 | 0.00028488 | 0.0 | 0.00 Modify | 0.83909 | 0.83909 | 0.83909 | 0.0 | 0.20 Other | | 0.1086 | | | 0.03 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 531310.0 ave 531310 max 531310 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 531310 Ave neighs/atom = 132.82750 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 293.062142221265, Press = -12.2121180798796 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 = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 5.370 | 5.370 | 5.370 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 8000 -13858.608 -13858.608 -14007.995 -14007.995 289.0005 289.0005 47971.786 47971.786 -1397.5379 -1397.5379 9000 -13855.399 -13855.399 -14005.368 -14005.368 290.12677 290.12677 47901.193 47901.193 1019.3209 1019.3209 Loop time of 429.651 on 1 procs for 1000 steps with 4000 atoms Performance: 0.201 ns/day, 119.348 hours/ns, 2.327 timesteps/s 99.8% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 428.51 | 428.51 | 428.51 | 0.0 | 99.74 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16872 | 0.16872 | 0.16872 | 0.0 | 0.04 Output | 0.00017718 | 0.00017718 | 0.00017718 | 0.0 | 0.00 Modify | 0.85893 | 0.85893 | 0.85893 | 0.0 | 0.20 Other | | 0.1105 | | | 0.03 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5842.00 ave 5842 max 5842 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: 530150.0 ave 530150 max 530150 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 530150 Ave neighs/atom = 132.53750 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 292.743708486337, Press = 11.5909950340859 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 = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 5.370 | 5.370 | 5.370 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 9000 -13855.399 -13855.399 -14005.368 -14005.368 290.12677 290.12677 47901.193 47901.193 1019.3209 1019.3209 10000 -13858.393 -13858.393 -14007.834 -14007.834 289.10437 289.10437 47922.369 47922.369 94.839529 94.839529 Loop time of 426.585 on 1 procs for 1000 steps with 4000 atoms Performance: 0.203 ns/day, 118.496 hours/ns, 2.344 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 | 425.48 | 425.48 | 425.48 | 0.0 | 99.74 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16331 | 0.16331 | 0.16331 | 0.0 | 0.04 Output | 0.00028087 | 0.00028087 | 0.00028087 | 0.0 | 0.00 Modify | 0.82876 | 0.82876 | 0.82876 | 0.0 | 0.19 Other | | 0.1086 | | | 0.03 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 530654.0 ave 530654 max 530654 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 530654 Ave neighs/atom = 132.66350 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 292.911950113028, Press = -0.976297577983023 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 = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 5.370 | 5.370 | 5.370 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 10000 -13858.393 -13858.393 -14007.834 -14007.834 289.10437 289.10437 47922.369 47922.369 94.839529 94.839529 11000 -13858.433 -13858.433 -14007.174 -14007.174 287.74927 287.74927 47929.495 47929.495 -88.938413 -88.938413 Loop time of 423.838 on 1 procs for 1000 steps with 4000 atoms Performance: 0.204 ns/day, 117.733 hours/ns, 2.359 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 | 422.75 | 422.75 | 422.75 | 0.0 | 99.74 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16254 | 0.16254 | 0.16254 | 0.0 | 0.04 Output | 0.0002726 | 0.0002726 | 0.0002726 | 0.0 | 0.00 Modify | 0.81977 | 0.81977 | 0.81977 | 0.0 | 0.19 Other | | 0.1089 | | | 0.03 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 530632.0 ave 530632 max 530632 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 530632 Ave neighs/atom = 132.65800 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 293.236925224355, Press = 3.67590997262094 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 = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 5.370 | 5.370 | 5.370 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 11000 -13858.433 -13858.433 -14007.174 -14007.174 287.74927 287.74927 47929.495 47929.495 -88.938413 -88.938413 12000 -13854.763 -13854.763 -14006.924 -14006.924 294.36568 294.36568 47909.753 47909.753 638.60896 638.60896 Loop time of 433.769 on 1 procs for 1000 steps with 4000 atoms Performance: 0.199 ns/day, 120.491 hours/ns, 2.305 timesteps/s 99.8% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 432.66 | 432.66 | 432.66 | 0.0 | 99.74 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16597 | 0.16597 | 0.16597 | 0.0 | 0.04 Output | 0.0001764 | 0.0001764 | 0.0001764 | 0.0 | 0.00 Modify | 0.83527 | 0.83527 | 0.83527 | 0.0 | 0.19 Other | | 0.108 | | | 0.02 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5843.00 ave 5843 max 5843 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: 530496.0 ave 530496 max 530496 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 530496 Ave neighs/atom = 132.62400 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 293.033828980731, Press = -0.69768942250814 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 = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 5.370 | 5.370 | 5.370 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 12000 -13854.763 -13854.763 -14006.924 -14006.924 294.36568 294.36568 47909.753 47909.753 638.60896 638.60896 13000 -13859.424 -13859.424 -14009.66 -14009.66 290.64131 290.64131 47968.386 47968.386 -1362.3192 -1362.3192 Loop time of 420.027 on 1 procs for 1000 steps with 4000 atoms Performance: 0.206 ns/day, 116.674 hours/ns, 2.381 timesteps/s 99.8% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 418.97 | 418.97 | 418.97 | 0.0 | 99.75 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16064 | 0.16064 | 0.16064 | 0.0 | 0.04 Output | 0.00017912 | 0.00017912 | 0.00017912 | 0.0 | 0.00 Modify | 0.79581 | 0.79581 | 0.79581 | 0.0 | 0.19 Other | | 0.1056 | | | 0.03 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 530658.0 ave 530658 max 530658 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 530658 Ave neighs/atom = 132.66450 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 292.9816082929, Press = 0.913730484803009 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 = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 5.370 | 5.370 | 5.370 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 13000 -13859.424 -13859.424 -14009.66 -14009.66 290.64131 290.64131 47968.386 47968.386 -1362.3192 -1362.3192 14000 -13856.538 -13856.538 -14006.477 -14006.477 290.06785 290.06785 47874.4 47874.4 1695.5172 1695.5172 Loop time of 365.301 on 1 procs for 1000 steps with 4000 atoms Performance: 0.237 ns/day, 101.472 hours/ns, 2.737 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 | 364.39 | 364.39 | 364.39 | 0.0 | 99.75 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.14309 | 0.14309 | 0.14309 | 0.0 | 0.04 Output | 0.00017754 | 0.00017754 | 0.00017754 | 0.0 | 0.00 Modify | 0.66606 | 0.66606 | 0.66606 | 0.0 | 0.18 Other | | 0.09879 | | | 0.03 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 530848.0 ave 530848 max 530848 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 530848 Ave neighs/atom = 132.71200 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 292.933902053877, Press = 3.07433764529737 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 = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 5.370 | 5.370 | 5.370 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 14000 -13856.538 -13856.538 -14006.477 -14006.477 290.06785 290.06785 47874.4 47874.4 1695.5172 1695.5172 15000 -13861.27 -13861.27 -14010.731 -14010.731 289.14301 289.14301 47983.18 47983.18 -1953.9293 -1953.9293 Loop time of 365.12 on 1 procs for 1000 steps with 4000 atoms Performance: 0.237 ns/day, 101.422 hours/ns, 2.739 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 | 364.21 | 364.21 | 364.21 | 0.0 | 99.75 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.14394 | 0.14394 | 0.14394 | 0.0 | 0.04 Output | 0.0001848 | 0.0001848 | 0.0001848 | 0.0 | 0.00 Modify | 0.66534 | 0.66534 | 0.66534 | 0.0 | 0.18 Other | | 0.09807 | | | 0.03 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 530922.0 ave 530922 max 530922 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 530922 Ave neighs/atom = 132.73050 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 292.769164257451, Press = -1.42594151862399 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 = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 5.370 | 5.370 | 5.370 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 15000 -13861.27 -13861.27 -14010.731 -14010.731 289.14301 289.14301 47983.18 47983.18 -1953.9293 -1953.9293 16000 -13851.703 -13851.703 -14004.743 -14004.743 296.06741 296.06741 47878.063 47878.063 1868.5115 1868.5115 Loop time of 376.546 on 1 procs for 1000 steps with 4000 atoms Performance: 0.229 ns/day, 104.596 hours/ns, 2.656 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 | 375.61 | 375.61 | 375.61 | 0.0 | 99.75 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.14784 | 0.14784 | 0.14784 | 0.0 | 0.04 Output | 0.00018014 | 0.00018014 | 0.00018014 | 0.0 | 0.00 Modify | 0.689 | 0.689 | 0.689 | 0.0 | 0.18 Other | | 0.0973 | | | 0.03 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 530492.0 ave 530492 max 530492 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 530492 Ave neighs/atom = 132.62300 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 292.820731648212, Press = 3.36343376498558 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 = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 5.370 | 5.370 | 5.370 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 16000 -13851.703 -13851.703 -14004.743 -14004.743 296.06741 296.06741 47878.063 47878.063 1868.5115 1868.5115 17000 -13858.205 -13858.205 -14010.154 -14010.154 293.95634 293.95634 47936.937 47936.937 -458.09631 -458.09631 Loop time of 367.279 on 1 procs for 1000 steps with 4000 atoms Performance: 0.235 ns/day, 102.022 hours/ns, 2.723 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 | 366.37 | 366.37 | 366.37 | 0.0 | 99.75 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.14697 | 0.14697 | 0.14697 | 0.0 | 0.04 Output | 0.00018234 | 0.00018234 | 0.00018234 | 0.0 | 0.00 Modify | 0.66937 | 0.66937 | 0.66937 | 0.0 | 0.18 Other | | 0.09436 | | | 0.03 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5843.00 ave 5843 max 5843 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: 530974.0 ave 530974 max 530974 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 530974 Ave neighs/atom = 132.74350 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 292.992690527145, Press = -0.528802322431409 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 = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 5.370 | 5.370 | 5.370 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 17000 -13858.205 -13858.205 -14010.154 -14010.154 293.95634 293.95634 47936.937 47936.937 -458.09631 -458.09631 18000 -13850.533 -13850.533 -14004.786 -14004.786 298.41362 298.41362 47931.842 47931.842 284.30468 284.30468 Loop time of 366.943 on 1 procs for 1000 steps with 4000 atoms Performance: 0.235 ns/day, 101.929 hours/ns, 2.725 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 | 366.03 | 366.03 | 366.03 | 0.0 | 99.75 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.14704 | 0.14704 | 0.14704 | 0.0 | 0.04 Output | 0.00022042 | 0.00022042 | 0.00022042 | 0.0 | 0.00 Modify | 0.66856 | 0.66856 | 0.66856 | 0.0 | 0.18 Other | | 0.09523 | | | 0.03 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 530788.0 ave 530788 max 530788 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 530788 Ave neighs/atom = 132.69700 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 293.073928682047, Press = 1.64585055050039 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 = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 5.370 | 5.370 | 5.370 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 18000 -13850.533 -13850.533 -14004.786 -14004.786 298.41362 298.41362 47931.842 47931.842 284.30468 284.30468 19000 -13859.881 -13859.881 -14009.38 -14009.38 289.21496 289.21496 47945.717 47945.717 -656.00392 -656.00392 Loop time of 367.559 on 1 procs for 1000 steps with 4000 atoms Performance: 0.235 ns/day, 102.100 hours/ns, 2.721 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 | 366.65 | 366.65 | 366.65 | 0.0 | 99.75 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.14658 | 0.14658 | 0.14658 | 0.0 | 0.04 Output | 0.00017635 | 0.00017635 | 0.00017635 | 0.0 | 0.00 Modify | 0.66798 | 0.66798 | 0.66798 | 0.0 | 0.18 Other | | 0.09508 | | | 0.03 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 530566.0 ave 530566 max 530566 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 530566 Ave neighs/atom = 132.64150 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 293.156511285101, Press = 0.728744423407984 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 = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 5.370 | 5.370 | 5.370 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 19000 -13859.881 -13859.881 -14009.38 -14009.38 289.21496 289.21496 47945.717 47945.717 -656.00392 -656.00392 20000 -13853.333 -13853.333 -14006.209 -14006.209 295.7492 295.7492 47894.677 47894.677 1177.0955 1177.0955 Loop time of 386.818 on 1 procs for 1000 steps with 4000 atoms Performance: 0.223 ns/day, 107.450 hours/ns, 2.585 timesteps/s 100.0% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 385.85 | 385.85 | 385.85 | 0.0 | 99.75 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.15236 | 0.15236 | 0.15236 | 0.0 | 0.04 Output | 0.00017728 | 0.00017728 | 0.00017728 | 0.0 | 0.00 Modify | 0.72083 | 0.72083 | 0.72083 | 0.0 | 0.19 Other | | 0.09915 | | | 0.03 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 530868.0 ave 530868 max 530868 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 530868 Ave neighs/atom = 132.71700 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 293.194681291064, Press = 2.69135698154323 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 = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 5.370 | 5.370 | 5.370 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 20000 -13853.333 -13853.333 -14006.209 -14006.209 295.7492 295.7492 47894.677 47894.677 1177.0955 1177.0955 21000 -13856.147 -13856.147 -14007.401 -14007.401 292.61186 292.61186 47970.336 47970.336 -1200.9952 -1200.9952 Loop time of 395.957 on 1 procs for 1000 steps with 4000 atoms Performance: 0.218 ns/day, 109.988 hours/ns, 2.526 timesteps/s 99.8% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 394.95 | 394.95 | 394.95 | 0.0 | 99.75 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.15643 | 0.15643 | 0.15643 | 0.0 | 0.04 Output | 0.00017721 | 0.00017721 | 0.00017721 | 0.0 | 0.00 Modify | 0.74964 | 0.74964 | 0.74964 | 0.0 | 0.19 Other | | 0.1027 | | | 0.03 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5842.00 ave 5842 max 5842 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: 530722.0 ave 530722 max 530722 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 530722 Ave neighs/atom = 132.68050 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 293.220660117673, Press = -1.44125948960271 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 = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 5.370 | 5.370 | 5.370 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 21000 -13856.147 -13856.147 -14007.401 -14007.401 292.61186 292.61186 47970.336 47970.336 -1200.9952 -1200.9952 22000 -13860.484 -13860.484 -14010.874 -14010.874 290.93858 290.93858 47869.188 47869.188 1267.4103 1267.4103 Loop time of 434.382 on 1 procs for 1000 steps with 4000 atoms Performance: 0.199 ns/day, 120.662 hours/ns, 2.302 timesteps/s 99.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 | 433.27 | 433.27 | 433.27 | 0.0 | 99.74 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16333 | 0.16333 | 0.16333 | 0.0 | 0.04 Output | 0.00017859 | 0.00017859 | 0.00017859 | 0.0 | 0.00 Modify | 0.84042 | 0.84042 | 0.84042 | 0.0 | 0.19 Other | | 0.1074 | | | 0.02 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5845.00 ave 5845 max 5845 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: 530510.0 ave 530510 max 530510 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 530510 Ave neighs/atom = 132.62750 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 293.188462880543, Press = 3.33858837139715 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 = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 5.370 | 5.370 | 5.370 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 22000 -13860.484 -13860.484 -14010.874 -14010.874 290.93858 290.93858 47869.188 47869.188 1267.4103 1267.4103 23000 -13860.903 -13860.903 -14009.121 -14009.121 286.73755 286.73755 47940.999 47940.999 -620.81721 -620.81721 Loop time of 390.978 on 1 procs for 1000 steps with 4000 atoms Performance: 0.221 ns/day, 108.605 hours/ns, 2.558 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 | 390 | 390 | 390 | 0.0 | 99.75 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.15027 | 0.15027 | 0.15027 | 0.0 | 0.04 Output | 0.00018142 | 0.00018142 | 0.00018142 | 0.0 | 0.00 Modify | 0.72777 | 0.72777 | 0.72777 | 0.0 | 0.19 Other | | 0.0987 | | | 0.03 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 531162.0 ave 531162 max 531162 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 531162 Ave neighs/atom = 132.79050 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 293.153302594588, Press = -1.27926697902654 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 = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 5.370 | 5.370 | 5.370 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 23000 -13860.903 -13860.903 -14009.121 -14009.121 286.73755 286.73755 47940.999 47940.999 -620.81721 -620.81721 24000 -13855.19 -13855.19 -14006.349 -14006.349 292.42598 292.42598 47907.584 47907.584 718.41939 718.41939 Loop time of 416.758 on 1 procs for 1000 steps with 4000 atoms Performance: 0.207 ns/day, 115.766 hours/ns, 2.399 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 | 415.69 | 415.69 | 415.69 | 0.0 | 99.74 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.15868 | 0.15868 | 0.15868 | 0.0 | 0.04 Output | 0.000179 | 0.000179 | 0.000179 | 0.0 | 0.00 Modify | 0.80872 | 0.80872 | 0.80872 | 0.0 | 0.19 Other | | 0.1058 | | | 0.03 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5845.00 ave 5845 max 5845 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: 530626.0 ave 530626 max 530626 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 530626 Ave neighs/atom = 132.65650 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 293.112629309257, Press = 1.39996097100085 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 = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 5.370 | 5.370 | 5.370 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 24000 -13855.19 -13855.19 -14006.349 -14006.349 292.42598 292.42598 47907.584 47907.584 718.41939 718.41939 25000 -13857.713 -13857.713 -14009.743 -14009.743 294.1141 294.1141 47927.224 47927.224 -197.46388 -197.46388 Loop time of 389.231 on 1 procs for 1000 steps with 4000 atoms Performance: 0.222 ns/day, 108.120 hours/ns, 2.569 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 | 388.25 | 388.25 | 388.25 | 0.0 | 99.75 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.15012 | 0.15012 | 0.15012 | 0.0 | 0.04 Output | 0.00022649 | 0.00022649 | 0.00022649 | 0.0 | 0.00 Modify | 0.73064 | 0.73064 | 0.73064 | 0.0 | 0.19 Other | | 0.1027 | | | 0.03 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 530372.0 ave 530372 max 530372 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 530372 Ave neighs/atom = 132.59300 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" jump SELF break # Write final averaged volume to file if temperature and volume have converged; otherwise wirte a # flag to indicate non-convergence. variable myStep equal step if "${myStep} < 2000000" then "print '${V}' file output/vol_T293.15.out" else "print 'not_converged' file output/vol_T293.15.out" print '${V}' file output/vol_T293.15.out 47924.2714369161 print "LAMMPS calculation completed" LAMMPS calculation completed quit 0