# 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 5.430949784815312*${_u_distance} variable latticeconst_converted equal 5.430949784815312*1 lattice diamond ${latticeconst_converted} lattice diamond 5.43094978481531 Lattice spacing in x,y,z = 5.4309498 5.4309498 5.4309498 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 (54.309498 54.309498 54.309498) 1 by 1 by 1 MPI processor grid create_atoms 1 box Created 8000 atoms create_atoms CPU = 0.001 seconds variable mass_converted equal 28.0855*${_u_mass} variable mass_converted equal 28.0855*1 kim_interactions Si WARNING: 'kim_' has been renamed to 'kim '. Please update your input. kim interactions Si #=== BEGIN kim interactions ================================== pair_style kim SW_HauchHollandMarder_1999Brittle_Si__MO_119167353542_005 pair_coeff * * Si #=== END kim interactions ==================================== mass 1 ${mass_converted} mass 1 28.0855 # initial volume variable v equal vol # assign formula variable V0 equal ${v} # evaluate initial value variable V0 equal 160187.034626825 variable V0_metal equal ${V0}/(${_u_distance}*${_u_distance}*${_u_distance}) variable V0_metal equal 160187.034626825/(${_u_distance}*${_u_distance}*${_u_distance}) variable V0_metal equal 160187.034626825/(1*${_u_distance}*${_u_distance}) variable V0_metal equal 160187.034626825/(1*1*${_u_distance}) variable V0_metal equal 160187.034626825/(1*1*1) variable V0_metal_times1000 equal ${V0_metal}*1000 variable V0_metal_times1000 equal 160187.034626825*1000 print "Initial system volume: ${V0_metal} Angstroms^3" Initial system volume: 160187.034626825 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 273.15*${_u_temperature} variable temp_converted equal 273.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 273.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 273.15 ${temp_converted} ${Tdamp_converted} iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 273.15 273.15 ${Tdamp_converted} iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 273.15 273.15 0.1 iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 273.15 273.15 0.1 iso 0 ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 273.15 273.15 0.1 iso 0 0 ${Pdamp_converted} fix ensemble all npt temp 273.15 273.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 "273.15 - 0.2" variable T_up equal "273.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_119167353542_005#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 = 5.77118 ghost atom cutoff = 5.77118 binsize = 2.88559, bins = 19 19 19 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.77118 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.307 | 4.307 | 4.307 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 0 -34408.776 -34408.776 -34691.2 -34691.2 273.15 273.15 160187.03 160187.03 1883.1827 1883.1827 1000 -34103.837 -34103.837 -34401.615 -34401.615 287.99984 287.99984 161070.74 161070.74 -771.86567 -771.86567 Loop time of 8.54887 on 1 procs for 1000 steps with 8000 atoms Performance: 10.107 ns/day, 2.375 hours/ns, 116.974 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 | 8.1973 | 8.1973 | 8.1973 | 0.0 | 95.89 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.041852 | 0.041852 | 0.041852 | 0.0 | 0.49 Output | 5.19e-05 | 5.19e-05 | 5.19e-05 | 0.0 | 0.00 Modify | 0.28 | 0.28 | 0.28 | 0.0 | 3.28 Other | | 0.02964 | | | 0.35 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 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: 272000.0 ave 272000 max 272000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 272000 Ave neighs/atom = 34.000000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 0, Press = 0 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 5.77118 ghost atom cutoff = 5.77118 binsize = 2.88559, bins = 19 19 19 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.77118 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.309 | 4.309 | 4.309 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 1000 -34103.837 -34103.837 -34401.615 -34401.615 287.99984 287.99984 161070.74 161070.74 -771.86567 -771.86567 2000 -34133.334 -34133.334 -34396.631 -34396.631 254.65075 254.65075 160896.65 160896.65 177.56081 177.56081 Loop time of 9.23311 on 1 procs for 1000 steps with 8000 atoms Performance: 9.358 ns/day, 2.565 hours/ns, 108.306 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 | 8.887 | 8.887 | 8.887 | 0.0 | 96.25 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.040599 | 0.040599 | 0.040599 | 0.0 | 0.44 Output | 6.49e-05 | 6.49e-05 | 6.49e-05 | 0.0 | 0.00 Modify | 0.27663 | 0.27663 | 0.27663 | 0.0 | 3.00 Other | | 0.02882 | | | 0.31 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 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: 272754.0 ave 272754 max 272754 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 272754 Ave neighs/atom = 34.094250 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 = 5.77118 ghost atom cutoff = 5.77118 binsize = 2.88559, bins = 19 19 19 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.77118 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.309 | 4.309 | 4.309 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 2000 -34133.334 -34133.334 -34396.631 -34396.631 254.65075 254.65075 160896.65 160896.65 177.56081 177.56081 3000 -34118.945 -34118.945 -34417.73 -34417.73 288.97432 288.97432 160907.66 160907.66 139.778 139.778 Loop time of 9.20483 on 1 procs for 1000 steps with 8000 atoms Performance: 9.386 ns/day, 2.557 hours/ns, 108.639 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 | 8.8606 | 8.8606 | 8.8606 | 0.0 | 96.26 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.040119 | 0.040119 | 0.040119 | 0.0 | 0.44 Output | 5.38e-05 | 5.38e-05 | 5.38e-05 | 0.0 | 0.00 Modify | 0.27566 | 0.27566 | 0.27566 | 0.0 | 2.99 Other | | 0.02839 | | | 0.31 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 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: 272672.0 ave 272672 max 272672 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 272672 Ave neighs/atom = 34.084000 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 = 5.77118 ghost atom cutoff = 5.77118 binsize = 2.88559, bins = 19 19 19 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.77118 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.309 | 4.309 | 4.309 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 3000 -34118.945 -34118.945 -34417.73 -34417.73 288.97432 288.97432 160907.66 160907.66 139.778 139.778 4000 -34123.744 -34123.744 -34406.423 -34406.423 273.39663 273.39663 160998.04 160998.04 -438.18477 -438.18477 Loop time of 9.29064 on 1 procs for 1000 steps with 8000 atoms Performance: 9.300 ns/day, 2.581 hours/ns, 107.635 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 | 8.9423 | 8.9423 | 8.9423 | 0.0 | 96.25 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.041042 | 0.041042 | 0.041042 | 0.0 | 0.44 Output | 4.37e-05 | 4.37e-05 | 4.37e-05 | 0.0 | 0.00 Modify | 0.27835 | 0.27835 | 0.27835 | 0.0 | 3.00 Other | | 0.02892 | | | 0.31 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 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: 272462.0 ave 272462 max 272462 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 272462 Ave neighs/atom = 34.057750 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 = 5.77118 ghost atom cutoff = 5.77118 binsize = 2.88559, bins = 19 19 19 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.77118 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.309 | 4.309 | 4.309 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 4000 -34123.744 -34123.744 -34406.423 -34406.423 273.39663 273.39663 160998.04 160998.04 -438.18477 -438.18477 5000 -34126.339 -34126.339 -34401.753 -34401.753 266.37077 266.37077 160828.41 160828.41 638.81915 638.81915 Loop time of 9.37027 on 1 procs for 1000 steps with 8000 atoms Performance: 9.221 ns/day, 2.603 hours/ns, 106.721 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 | 9.0191 | 9.0191 | 9.0191 | 0.0 | 96.25 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.041216 | 0.041216 | 0.041216 | 0.0 | 0.44 Output | 3.54e-05 | 3.54e-05 | 3.54e-05 | 0.0 | 0.00 Modify | 0.28055 | 0.28055 | 0.28055 | 0.0 | 2.99 Other | | 0.02932 | | | 0.31 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 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: 272584.0 ave 272584 max 272584 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 272584 Ave neighs/atom = 34.073000 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 = 273.350634849694, Press = -136.60805669215 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 = 5.77118 ghost atom cutoff = 5.77118 binsize = 2.88559, bins = 19 19 19 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.77118 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.309 | 4.309 | 4.309 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 5000 -34126.339 -34126.339 -34401.753 -34401.753 266.37077 266.37077 160828.41 160828.41 638.81915 638.81915 6000 -34119.462 -34119.462 -34402.4 -34402.4 273.647 273.647 161060.66 161060.66 -764.27 -764.27 Loop time of 9.08942 on 1 procs for 1000 steps with 8000 atoms Performance: 9.506 ns/day, 2.525 hours/ns, 110.018 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 | 8.7424 | 8.7424 | 8.7424 | 0.0 | 96.18 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.039961 | 0.039961 | 0.039961 | 0.0 | 0.44 Output | 4.54e-05 | 4.54e-05 | 4.54e-05 | 0.0 | 0.00 Modify | 0.2786 | 0.2786 | 0.2786 | 0.0 | 3.07 Other | | 0.02839 | | | 0.31 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 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: 272606.0 ave 272606 max 272606 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 272606 Ave neighs/atom = 34.075750 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 = 272.69857177986, Press = 8.36948829059871 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 = 5.77118 ghost atom cutoff = 5.77118 binsize = 2.88559, bins = 19 19 19 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.77118 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.309 | 4.309 | 4.309 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 6000 -34119.462 -34119.462 -34402.4 -34402.4 273.647 273.647 161060.66 161060.66 -764.27 -764.27 7000 -34127.853 -34127.853 -34422.603 -34422.603 285.07186 285.07186 160925.03 160925.03 -40.99568 -40.99568 Loop time of 9.4346 on 1 procs for 1000 steps with 8000 atoms Performance: 9.158 ns/day, 2.621 hours/ns, 105.993 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 | 9.0757 | 9.0757 | 9.0757 | 0.0 | 96.20 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.042022 | 0.042022 | 0.042022 | 0.0 | 0.45 Output | 3.59e-05 | 3.59e-05 | 3.59e-05 | 0.0 | 0.00 Modify | 0.28739 | 0.28739 | 0.28739 | 0.0 | 3.05 Other | | 0.02944 | | | 0.31 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 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: 272594.0 ave 272594 max 272594 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 272594 Ave neighs/atom = 34.074250 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 = 273.045520114659, Press = -16.6154287144471 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 = 5.77118 ghost atom cutoff = 5.77118 binsize = 2.88559, bins = 19 19 19 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.77118 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.309 | 4.309 | 4.309 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 7000 -34127.853 -34127.853 -34422.603 -34422.603 285.07186 285.07186 160925.03 160925.03 -40.99568 -40.99568 8000 -34119.399 -34119.399 -34402.505 -34402.505 273.80936 273.80936 160895.16 160895.16 271.31078 271.31078 Loop time of 9.19785 on 1 procs for 1000 steps with 8000 atoms Performance: 9.393 ns/day, 2.555 hours/ns, 108.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 | 8.8483 | 8.8483 | 8.8483 | 0.0 | 96.20 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.040184 | 0.040184 | 0.040184 | 0.0 | 0.44 Output | 4.45e-05 | 4.45e-05 | 4.45e-05 | 0.0 | 0.00 Modify | 0.28083 | 0.28083 | 0.28083 | 0.0 | 3.05 Other | | 0.02851 | | | 0.31 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 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: 272392.0 ave 272392 max 272392 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 272392 Ave neighs/atom = 34.049000 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 = 273.100345588096, Press = -0.540493901142272 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 = 5.77118 ghost atom cutoff = 5.77118 binsize = 2.88559, bins = 19 19 19 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.77118 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.309 | 4.309 | 4.309 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 8000 -34119.399 -34119.399 -34402.505 -34402.505 273.80936 273.80936 160895.16 160895.16 271.31078 271.31078 9000 -34126.352 -34126.352 -34405.494 -34405.494 269.97558 269.97558 160995.58 160995.58 -412.01011 -412.01011 Loop time of 9.45053 on 1 procs for 1000 steps with 8000 atoms Performance: 9.142 ns/day, 2.625 hours/ns, 105.814 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 | 9.094 | 9.094 | 9.094 | 0.0 | 96.23 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.041018 | 0.041018 | 0.041018 | 0.0 | 0.43 Output | 3.57e-05 | 3.57e-05 | 3.57e-05 | 0.0 | 0.00 Modify | 0.28627 | 0.28627 | 0.28627 | 0.0 | 3.03 Other | | 0.02922 | | | 0.31 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 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: 272600.0 ave 272600 max 272600 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 272600 Ave neighs/atom = 34.075000 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 = 273.037400221908, Press = -3.40823833065307 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 = 5.77118 ghost atom cutoff = 5.77118 binsize = 2.88559, bins = 19 19 19 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.77118 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.309 | 4.309 | 4.309 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 9000 -34126.352 -34126.352 -34405.494 -34405.494 269.97558 269.97558 160995.58 160995.58 -412.01011 -412.01011 10000 -34123.475 -34123.475 -34411.241 -34411.241 278.31629 278.31629 160859.45 160859.45 438.60209 438.60209 Loop time of 9.24242 on 1 procs for 1000 steps with 8000 atoms Performance: 9.348 ns/day, 2.567 hours/ns, 108.197 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 | 8.891 | 8.891 | 8.891 | 0.0 | 96.20 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.04061 | 0.04061 | 0.04061 | 0.0 | 0.44 Output | 3.57e-05 | 3.57e-05 | 3.57e-05 | 0.0 | 0.00 Modify | 0.28196 | 0.28196 | 0.28196 | 0.0 | 3.05 Other | | 0.02882 | | | 0.31 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 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: 272532.0 ave 272532 max 272532 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 272532 Ave neighs/atom = 34.066500 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 = 273.13577381866, Press = -3.22580242139122 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 = 5.77118 ghost atom cutoff = 5.77118 binsize = 2.88559, bins = 19 19 19 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.77118 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.309 | 4.309 | 4.309 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 10000 -34123.475 -34123.475 -34411.241 -34411.241 278.31629 278.31629 160859.45 160859.45 438.60209 438.60209 11000 -34115.03 -34115.03 -34402.981 -34402.981 278.49601 278.49601 160996.76 160996.76 -335.31777 -335.31777 Loop time of 9.22983 on 1 procs for 1000 steps with 8000 atoms Performance: 9.361 ns/day, 2.564 hours/ns, 108.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 | 8.8791 | 8.8791 | 8.8791 | 0.0 | 96.20 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.040442 | 0.040442 | 0.040442 | 0.0 | 0.44 Output | 3.51e-05 | 3.51e-05 | 3.51e-05 | 0.0 | 0.00 Modify | 0.28163 | 0.28163 | 0.28163 | 0.0 | 3.05 Other | | 0.02867 | | | 0.31 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 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: 272558.0 ave 272558 max 272558 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 272558 Ave neighs/atom = 34.069750 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 = 273.147714357316, Press = 0.951399964267665 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 = 5.77118 ghost atom cutoff = 5.77118 binsize = 2.88559, bins = 19 19 19 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.77118 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.309 | 4.309 | 4.309 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 11000 -34115.03 -34115.03 -34402.981 -34402.981 278.49601 278.49601 160996.76 160996.76 -335.31777 -335.31777 12000 -34124.556 -34124.556 -34413.734 -34413.734 279.68192 279.68192 160978.53 160978.53 -331.72468 -331.72468 Loop time of 9.17476 on 1 procs for 1000 steps with 8000 atoms Performance: 9.417 ns/day, 2.549 hours/ns, 108.995 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 | 8.8255 | 8.8255 | 8.8255 | 0.0 | 96.19 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.040063 | 0.040063 | 0.040063 | 0.0 | 0.44 Output | 3.53e-05 | 3.53e-05 | 3.53e-05 | 0.0 | 0.00 Modify | 0.28081 | 0.28081 | 0.28081 | 0.0 | 3.06 Other | | 0.02839 | | | 0.31 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 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: 272496.0 ave 272496 max 272496 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 272496 Ave neighs/atom = 34.062000 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 = 273.250399711478, Press = -5.83147750966053 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 = 5.77118 ghost atom cutoff = 5.77118 binsize = 2.88559, bins = 19 19 19 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.77118 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.309 | 4.309 | 4.309 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 12000 -34124.556 -34124.556 -34413.734 -34413.734 279.68192 279.68192 160978.53 160978.53 -331.72468 -331.72468 13000 -34116.793 -34116.793 -34409.083 -34409.083 282.69182 282.69182 160869.15 160869.15 428.39546 428.39546 Loop time of 9.11081 on 1 procs for 1000 steps with 8000 atoms Performance: 9.483 ns/day, 2.531 hours/ns, 109.760 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 | 8.7647 | 8.7647 | 8.7647 | 0.0 | 96.20 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.03977 | 0.03977 | 0.03977 | 0.0 | 0.44 Output | 3.57e-05 | 3.57e-05 | 3.57e-05 | 0.0 | 0.00 Modify | 0.27823 | 0.27823 | 0.27823 | 0.0 | 3.05 Other | | 0.02812 | | | 0.31 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 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: 272542.0 ave 272542 max 272542 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 272542 Ave neighs/atom = 34.067750 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 = 273.372011217816, Press = -0.181037633112522 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 = 5.77118 ghost atom cutoff = 5.77118 binsize = 2.88559, bins = 19 19 19 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.77118 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.309 | 4.309 | 4.309 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 13000 -34116.793 -34116.793 -34409.083 -34409.083 282.69182 282.69182 160869.15 160869.15 428.39546 428.39546 14000 -34120.916 -34120.916 -34406.552 -34406.552 276.2569 276.2569 160989.53 160989.53 -331.14711 -331.14711 Loop time of 9.17391 on 1 procs for 1000 steps with 8000 atoms Performance: 9.418 ns/day, 2.548 hours/ns, 109.005 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 | 8.8272 | 8.8272 | 8.8272 | 0.0 | 96.22 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.039519 | 0.039519 | 0.039519 | 0.0 | 0.43 Output | 3.57e-05 | 3.57e-05 | 3.57e-05 | 0.0 | 0.00 Modify | 0.27914 | 0.27914 | 0.27914 | 0.0 | 3.04 Other | | 0.02802 | | | 0.31 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 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: 272592.0 ave 272592 max 272592 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 272592 Ave neighs/atom = 34.074000 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 = 273.358284494616, Press = -1.22890366204803 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 = 5.77118 ghost atom cutoff = 5.77118 binsize = 2.88559, bins = 19 19 19 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.77118 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.309 | 4.309 | 4.309 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 14000 -34120.916 -34120.916 -34406.552 -34406.552 276.2569 276.2569 160989.53 160989.53 -331.14711 -331.14711 15000 -34123.288 -34123.288 -34401.261 -34401.261 268.8453 268.8453 160894.53 160894.53 264.18755 264.18755 Loop time of 8.52191 on 1 procs for 1000 steps with 8000 atoms Performance: 10.139 ns/day, 2.367 hours/ns, 117.345 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 | 8.1942 | 8.1942 | 8.1942 | 0.0 | 96.16 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.037535 | 0.037535 | 0.037535 | 0.0 | 0.44 Output | 4.48e-05 | 4.48e-05 | 4.48e-05 | 0.0 | 0.00 Modify | 0.26424 | 0.26424 | 0.26424 | 0.0 | 3.10 Other | | 0.02584 | | | 0.30 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 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: 272530.0 ave 272530 max 272530 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 272530 Ave neighs/atom = 34.066250 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 = 273.385405276125, Press = -1.71947128216179 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 = 5.77118 ghost atom cutoff = 5.77118 binsize = 2.88559, bins = 19 19 19 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.77118 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.309 | 4.309 | 4.309 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 15000 -34123.288 -34123.288 -34401.261 -34401.261 268.8453 268.8453 160894.53 160894.53 264.18755 264.18755 16000 -34122.697 -34122.697 -34405.918 -34405.918 273.92027 273.92027 160954.55 160954.55 -130.40394 -130.40394 Loop time of 8.45391 on 1 procs for 1000 steps with 8000 atoms Performance: 10.220 ns/day, 2.348 hours/ns, 118.288 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 | 8.1295 | 8.1295 | 8.1295 | 0.0 | 96.16 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.03706 | 0.03706 | 0.03706 | 0.0 | 0.44 Output | 3.51e-05 | 3.51e-05 | 3.51e-05 | 0.0 | 0.00 Modify | 0.26213 | 0.26213 | 0.26213 | 0.0 | 3.10 Other | | 0.02515 | | | 0.30 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 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: 272596.0 ave 272596 max 272596 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 272596 Ave neighs/atom = 34.074500 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 = 273.369468433761, Press = 0.337920636122811 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 = 5.77118 ghost atom cutoff = 5.77118 binsize = 2.88559, bins = 19 19 19 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.77118 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.309 | 4.309 | 4.309 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 16000 -34122.697 -34122.697 -34405.918 -34405.918 273.92027 273.92027 160954.55 160954.55 -130.40394 -130.40394 17000 -34126.603 -34126.603 -34409.25 -34409.25 273.36576 273.36576 161009.34 161009.34 -496.27315 -496.27315 Loop time of 8.59161 on 1 procs for 1000 steps with 8000 atoms Performance: 10.056 ns/day, 2.387 hours/ns, 116.393 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 | 8.2632 | 8.2632 | 8.2632 | 0.0 | 96.18 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.037247 | 0.037247 | 0.037247 | 0.0 | 0.43 Output | 3.52e-05 | 3.52e-05 | 3.52e-05 | 0.0 | 0.00 Modify | 0.26549 | 0.26549 | 0.26549 | 0.0 | 3.09 Other | | 0.02558 | | | 0.30 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 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: 272556.0 ave 272556 max 272556 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 272556 Ave neighs/atom = 34.069500 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 = 273.432768221064, Press = -2.92932137209773 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 = 5.77118 ghost atom cutoff = 5.77118 binsize = 2.88559, bins = 19 19 19 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.77118 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.309 | 4.309 | 4.309 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 17000 -34126.603 -34126.603 -34409.25 -34409.25 273.36576 273.36576 161009.34 161009.34 -496.27315 -496.27315 18000 -34122.467 -34122.467 -34401.265 -34401.265 269.64278 269.64278 160799.2 160799.2 861.16376 861.16376 Loop time of 8.43145 on 1 procs for 1000 steps with 8000 atoms Performance: 10.247 ns/day, 2.342 hours/ns, 118.604 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 | 8.1073 | 8.1073 | 8.1073 | 0.0 | 96.15 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.036955 | 0.036955 | 0.036955 | 0.0 | 0.44 Output | 3.59e-05 | 3.59e-05 | 3.59e-05 | 0.0 | 0.00 Modify | 0.262 | 0.262 | 0.262 | 0.0 | 3.11 Other | | 0.02521 | | | 0.30 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 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: 272412.0 ave 272412 max 272412 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 272412 Ave neighs/atom = 34.051500 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 = 273.4330166084, Press = -0.0210951696102594 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 = 5.77118 ghost atom cutoff = 5.77118 binsize = 2.88559, bins = 19 19 19 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.77118 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.309 | 4.309 | 4.309 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 18000 -34122.467 -34122.467 -34401.265 -34401.265 269.64278 269.64278 160799.2 160799.2 861.16376 861.16376 19000 -34121.473 -34121.473 -34407.967 -34407.967 277.08661 277.08661 161043.48 161043.48 -666.28754 -666.28754 Loop time of 8.72941 on 1 procs for 1000 steps with 8000 atoms Performance: 9.898 ns/day, 2.425 hours/ns, 114.555 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 | 8.3974 | 8.3974 | 8.3974 | 0.0 | 96.20 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.037531 | 0.037531 | 0.037531 | 0.0 | 0.43 Output | 3.55e-05 | 3.55e-05 | 3.55e-05 | 0.0 | 0.00 Modify | 0.26833 | 0.26833 | 0.26833 | 0.0 | 3.07 Other | | 0.02613 | | | 0.30 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 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: 272680.0 ave 272680 max 272680 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 272680 Ave neighs/atom = 34.085000 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 = 273.390267389278, Press = -0.54105433054321 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 = 5.77118 ghost atom cutoff = 5.77118 binsize = 2.88559, bins = 19 19 19 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.77118 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.309 | 4.309 | 4.309 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 19000 -34121.473 -34121.473 -34407.967 -34407.967 277.08661 277.08661 161043.48 161043.48 -666.28754 -666.28754 20000 -34123.572 -34123.572 -34406.968 -34406.968 274.0898 274.0898 160874.73 160874.73 366.95313 366.95313 Loop time of 8.79356 on 1 procs for 1000 steps with 8000 atoms Performance: 9.825 ns/day, 2.443 hours/ns, 113.720 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 | 8.4593 | 8.4593 | 8.4593 | 0.0 | 96.20 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.03782 | 0.03782 | 0.03782 | 0.0 | 0.43 Output | 3.59e-05 | 3.59e-05 | 3.59e-05 | 0.0 | 0.00 Modify | 0.27024 | 0.27024 | 0.27024 | 0.0 | 3.07 Other | | 0.02614 | | | 0.30 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 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: 272482.0 ave 272482 max 272482 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 272482 Ave neighs/atom = 34.060250 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 = 273.404660444563, Press = -1.49918835736918 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 = 5.77118 ghost atom cutoff = 5.77118 binsize = 2.88559, bins = 19 19 19 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.77118 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.309 | 4.309 | 4.309 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 20000 -34123.572 -34123.572 -34406.968 -34406.968 274.0898 274.0898 160874.73 160874.73 366.95313 366.95313 21000 -34120.31 -34120.31 -34401.751 -34401.751 272.19927 272.19927 160937.37 160937.37 32.343023 32.343023 Loop time of 8.55876 on 1 procs for 1000 steps with 8000 atoms Performance: 10.095 ns/day, 2.377 hours/ns, 116.839 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 | 8.2307 | 8.2307 | 8.2307 | 0.0 | 96.17 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.03723 | 0.03723 | 0.03723 | 0.0 | 0.43 Output | 3.57e-05 | 3.57e-05 | 3.57e-05 | 0.0 | 0.00 Modify | 0.26518 | 0.26518 | 0.26518 | 0.0 | 3.10 Other | | 0.02559 | | | 0.30 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 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: 272574.0 ave 272574 max 272574 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 272574 Ave neighs/atom = 34.071750 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 = 273.39734757244, Press = 0.364529825054857 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 = 5.77118 ghost atom cutoff = 5.77118 binsize = 2.88559, bins = 19 19 19 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.77118 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.309 | 4.309 | 4.309 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 21000 -34120.31 -34120.31 -34401.751 -34401.751 272.19927 272.19927 160937.37 160937.37 32.343023 32.343023 22000 -34116.902 -34116.902 -34401.898 -34401.898 275.63811 275.63811 161012.17 161012.17 -435.36417 -435.36417 Loop time of 8.58078 on 1 procs for 1000 steps with 8000 atoms Performance: 10.069 ns/day, 2.384 hours/ns, 116.539 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 | 8.2519 | 8.2519 | 8.2519 | 0.0 | 96.17 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.037233 | 0.037233 | 0.037233 | 0.0 | 0.43 Output | 4.54e-05 | 4.54e-05 | 4.54e-05 | 0.0 | 0.00 Modify | 0.26617 | 0.26617 | 0.26617 | 0.0 | 3.10 Other | | 0.02548 | | | 0.30 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 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: 272526.0 ave 272526 max 272526 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 272526 Ave neighs/atom = 34.065750 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 = 273.406101076387, Press = -1.45641396807837 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 = 5.77118 ghost atom cutoff = 5.77118 binsize = 2.88559, bins = 19 19 19 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.77118 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.309 | 4.309 | 4.309 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 22000 -34116.902 -34116.902 -34401.898 -34401.898 275.63811 275.63811 161012.17 161012.17 -435.36417 -435.36417 23000 -34126.726 -34126.726 -34405.202 -34405.202 269.3314 269.3314 160773.51 160773.51 986.06176 986.06176 Loop time of 7.97649 on 1 procs for 1000 steps with 8000 atoms Performance: 10.832 ns/day, 2.216 hours/ns, 125.368 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 | 7.6679 | 7.6679 | 7.6679 | 0.0 | 96.13 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.035377 | 0.035377 | 0.035377 | 0.0 | 0.44 Output | 3.53e-05 | 3.53e-05 | 3.53e-05 | 0.0 | 0.00 Modify | 0.25063 | 0.25063 | 0.25063 | 0.0 | 3.14 Other | | 0.02258 | | | 0.28 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 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: 272488.0 ave 272488 max 272488 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 272488 Ave neighs/atom = 34.061000 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 = 273.439054439421, Press = -0.169580910701783 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 = 5.77118 ghost atom cutoff = 5.77118 binsize = 2.88559, bins = 19 19 19 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.77118 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.309 | 4.309 | 4.309 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 23000 -34126.726 -34126.726 -34405.202 -34405.202 269.3314 269.3314 160773.51 160773.51 986.06176 986.06176 24000 -34129.582 -34129.582 -34406.891 -34406.891 268.20271 268.20271 161092.12 161092.12 -1042.5867 -1042.5867 Loop time of 8.26277 on 1 procs for 1000 steps with 8000 atoms Performance: 10.457 ns/day, 2.295 hours/ns, 121.025 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 | 7.946 | 7.946 | 7.946 | 0.0 | 96.17 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.03607 | 0.03607 | 0.03607 | 0.0 | 0.44 Output | 3.53e-05 | 3.53e-05 | 3.53e-05 | 0.0 | 0.00 Modify | 0.25704 | 0.25704 | 0.25704 | 0.0 | 3.11 Other | | 0.02359 | | | 0.29 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 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: 272572.0 ave 272572 max 272572 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 272572 Ave neighs/atom = 34.071500 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 = 273.3961962429, Press = -0.0134127296585117 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 = 5.77118 ghost atom cutoff = 5.77118 binsize = 2.88559, bins = 19 19 19 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.77118 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.309 | 4.309 | 4.309 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 24000 -34129.582 -34129.582 -34406.891 -34406.891 268.20271 268.20271 161092.12 161092.12 -1042.5867 -1042.5867 25000 -34120.268 -34120.268 -34404.435 -34404.435 274.8358 274.8358 160904 160904 210.66526 210.66526 Loop time of 8.09604 on 1 procs for 1000 steps with 8000 atoms Performance: 10.672 ns/day, 2.249 hours/ns, 123.517 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 | 7.7837 | 7.7837 | 7.7837 | 0.0 | 96.14 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.035844 | 0.035844 | 0.035844 | 0.0 | 0.44 Output | 3.49e-05 | 3.49e-05 | 3.49e-05 | 0.0 | 0.00 Modify | 0.2533 | 0.2533 | 0.2533 | 0.0 | 3.13 Other | | 0.02318 | | | 0.29 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 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: 272506.0 ave 272506 max 272506 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 272506 Ave neighs/atom = 34.063250 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 = 273.33683878591, Press = -1.11235519947002 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 = 5.77118 ghost atom cutoff = 5.77118 binsize = 2.88559, bins = 19 19 19 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.77118 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.309 | 4.309 | 4.309 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 25000 -34120.268 -34120.268 -34404.435 -34404.435 274.8358 274.8358 160904 160904 210.66526 210.66526 26000 -34128.215 -34128.215 -34410.864 -34410.864 273.36812 273.36812 160909.77 160909.77 99.047423 99.047423 Loop time of 7.95586 on 1 procs for 1000 steps with 8000 atoms Performance: 10.860 ns/day, 2.210 hours/ns, 125.693 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 | 7.6473 | 7.6473 | 7.6473 | 0.0 | 96.12 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.035332 | 0.035332 | 0.035332 | 0.0 | 0.44 Output | 5.89e-05 | 5.89e-05 | 5.89e-05 | 0.0 | 0.00 Modify | 0.25039 | 0.25039 | 0.25039 | 0.0 | 3.15 Other | | 0.02282 | | | 0.29 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 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: 272624.0 ave 272624 max 272624 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 272624 Ave neighs/atom = 34.078000 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_T273.15.out" else "print 'not_converged' file output/vol_T273.15.out" print '${V}' file output/vol_T273.15.out 160934.404080676 print "LAMMPS calculation completed" LAMMPS calculation completed quit 0