# 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 4.195095829665661*${_u_distance}
variable latticeconst_converted equal 4.195095829665661*1
lattice       bcc  ${latticeconst_converted}
lattice       bcc  4.19509582966566
Lattice spacing in x,y,z = 4.1950958 4.1950958 4.1950958
region        simbox  block 0 10 0 10 0 10 units lattice
create_box    1 simbox
Created orthogonal box = (0 0 0) to (41.950958 41.950958 41.950958)
  1 by 1 by 1 MPI processor grid
create_atoms  1 box
Created 2000 atoms
  using lattice units in orthogonal box = (0 0 0) to (41.950958 41.950958 41.950958)
  create_atoms CPU = 0.002 seconds

variable mass_converted equal 22.98977*${_u_mass}
variable mass_converted equal 22.98977*1

kim_interactions Na

WARNING: 'kim_<command>' has been renamed to 'kim <command>'. Please update your input.

kim interactions Na
#=== BEGIN kim interactions ==================================
pair_style kim MEAM_LAMMPS_KimKoLee_2020_Na__MO_321355778754_002
pair_coeff * * Na
#=== END kim interactions ====================================


mass          1 ${mass_converted}
mass          1 22.98977

# initial volume
variable      v equal vol        # assign formula
variable      V0 equal ${v}    # evaluate initial value
variable      V0 equal 73828.7742291291    
variable      V0_metal equal ${V0}/(${_u_distance}*${_u_distance}*${_u_distance})
variable      V0_metal equal 73828.7742291291/(${_u_distance}*${_u_distance}*${_u_distance})
variable      V0_metal equal 73828.7742291291/(1*${_u_distance}*${_u_distance})
variable      V0_metal equal 73828.7742291291/(1*1*${_u_distance})
variable      V0_metal equal 73828.7742291291/(1*1*1)
variable      V0_metal_times1000 equal ${V0_metal}*1000
variable      V0_metal_times1000 equal 73828.7742291291*1000

print "Initial system volume: ${V0_metal} Angstroms^3"
Initial system volume: 73828.7742291291 Angstroms^3

# set the time step to 0.001 picoseconds
variable timestep_converted equal 0.001*${_u_time}
variable timestep_converted equal 0.001*1
timestep ${timestep_converted}
timestep 0.001

variable temp_converted equal 253.15*${_u_temperature}
variable temp_converted equal 253.15*1
variable Tdamp_converted equal 0.1*${_u_time}
variable Tdamp_converted equal 0.1*1
variable press_converted equal 0.0*${_u_pressure}
variable press_converted equal 0.0*1
variable Pdamp_converted equal 1*${_u_time}
variable Pdamp_converted equal 1*1

# create initial velocities consistent with the chosen temperature
velocity      all create ${temp_converted} 17 mom yes rot yes
velocity      all create 253.15 17 mom yes rot yes

# set NPT ensemble for all atoms
fix           ensemble all npt temp ${temp_converted} ${temp_converted} ${Tdamp_converted}               iso ${press_converted} ${press_converted} ${Pdamp_converted}
fix           ensemble all npt temp 253.15 ${temp_converted} ${Tdamp_converted}               iso ${press_converted} ${press_converted} ${Pdamp_converted}
fix           ensemble all npt temp 253.15 253.15 ${Tdamp_converted}               iso ${press_converted} ${press_converted} ${Pdamp_converted}
fix           ensemble all npt temp 253.15 253.15 0.1               iso ${press_converted} ${press_converted} ${Pdamp_converted}
fix           ensemble all npt temp 253.15 253.15 0.1               iso 0 ${press_converted} ${Pdamp_converted}
fix           ensemble all npt temp 253.15 253.15 0.1               iso 0 0 ${Pdamp_converted}
fix           ensemble all npt temp 253.15 253.15 0.1               iso 0 0 1

# compute the time averages of pressure, temperature, and volume, respectively
# ignore the first 5000 timesteps
variable      etotal_metal equal etotal/${_u_energy}
variable      etotal_metal equal etotal/1
variable      pe_metal equal pe/${_u_energy}
variable      pe_metal equal pe/1
variable      T_metal equal temp/${_u_temperature}
variable      T_metal equal temp/1
variable      V_metal equal vol/(${_u_distance}*${_u_distance}*${_u_distance})
variable      V_metal equal vol/(1*${_u_distance}*${_u_distance})
variable      V_metal equal vol/(1*1*${_u_distance})
variable      V_metal equal vol/(1*1*1)
variable      P_metal equal press/${_u_pressure}
variable      P_metal equal press/1
fix           avgmyTemp  all ave/time 5 20 100 v_T_metal  ave running start 5000
fix           avgmyPress all ave/time 5 20 100 v_P_metal ave running start 5000
fix           avgmyVol   all ave/time 5 20 100 v_V_metal   ave running start 5000

# extract fix quantities into variables so they can be used in if-else logic later.
variable      T equal f_avgmyTemp
variable      P equal f_avgmyPress
variable      V equal f_avgmyVol

# set error bounds for temperature and pressure in original metal units (K and bar)
variable      T_low equal  "253.15 - 0.2"
variable      T_up  equal  "253.15 + 0.2"
variable      P_low equal  "0.0 - 0.2"
variable      P_up  equal  "0.0 + 0.2"

# print to logfile every 1000 timesteps
thermo_style  custom step etotal v_etotal_metal pe v_pe_metal               temp v_T_metal vol v_V_metal press v_P_metal
thermo        1000

# Run a simulation for at most 2000*1000 timesteps. At each 1000th time step, check
# whether the temperature and pressure have converged. If yes, break.
label top
variable a loop 2000
run 1000

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_321355778754_002#item-citation

CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE

Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule
Neighbor list info ...
  update every 1 steps, delay 10 steps, check yes
  max neighbors/atom: 2000, page size: 100000
  master list distance cutoff = 6.8
  ghost atom cutoff = 6.8
  binsize = 3.4, bins = 13 13 13
  1 neighbor lists, perpetual/occasional/extra = 1 0 0
  (1) pair kim, perpetual
      attributes: full, newton off, cut 6.8
      pair build: full/bin/atomonly
      stencil: full/bin/3d
      bin: standard
Per MPI rank memory allocation (min/avg/max) = 3.326 | 3.326 | 3.326 Mbytes
   Step         TotEng     v_etotal_metal     PotEng       v_pe_metal        Temp        v_T_metal        Volume       v_V_metal        Press        v_P_metal   
         0  -2160.5227     -2160.5227     -2225.9344     -2225.9344      253.15         253.15         73828.774      73828.774      946.34215      946.34215    
      1000  -2085.0463     -2085.0463     -2148.8478     -2148.8478      246.9181       246.9181       77192.139      77192.139     -308.82832     -308.82832    
Loop time of 79.9562 on 1 procs for 1000 steps with 2000 atoms

Performance: 1.081 ns/day, 22.210 hours/ns, 12.507 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    | 79.321     | 79.321     | 79.321     |   0.0 | 99.21
Neigh   | 0.15116    | 0.15116    | 0.15116    |   0.0 |  0.19
Comm    | 0.080459   | 0.080459   | 0.080459   |   0.0 |  0.10
Output  | 0.00021182 | 0.00021182 | 0.00021182 |   0.0 |  0.00
Modify  | 0.34875    | 0.34875    | 0.34875    |   0.0 |  0.44
Other   |            | 0.05486    |            |       |  0.07

Nlocal:           2000 ave        2000 max        2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:           2838 ave        2838 max        2838 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Neighs:              0 ave           0 max           0 min
Histogram: 1 0 0 0 0 0 0 0 0 0
FullNghs:        62838 ave       62838 max       62838 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 62838
Ave neighs/atom = 31.419
Neighbor list builds = 11
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
Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule
Neighbor list info ...
  update every 1 steps, delay 10 steps, check yes
  max neighbors/atom: 2000, page size: 100000
  master list distance cutoff = 6.8
  ghost atom cutoff = 6.8
  binsize = 3.4, bins = 13 13 13
  1 neighbor lists, perpetual/occasional/extra = 1 0 0
  (1) pair kim, perpetual
      attributes: full, newton off, cut 6.8
      pair build: full/bin/atomonly
      stencil: full/bin/3d
      bin: standard
Per MPI rank memory allocation (min/avg/max) = 3.325 | 3.325 | 3.325 Mbytes
   Step         TotEng     v_etotal_metal     PotEng       v_pe_metal        Temp        v_T_metal        Volume       v_V_metal        Press        v_P_metal   
      1000  -2085.0463     -2085.0463     -2148.8478     -2148.8478      246.9181       246.9181       77192.139      77192.139     -308.82832     -308.82832    
      2000  -2091.894      -2091.894      -2157.4622     -2157.4622      253.75571      253.75571      76726.53       76726.53      -88.441783     -88.441783    
Loop time of 79.5486 on 1 procs for 1000 steps with 2000 atoms

Performance: 1.086 ns/day, 22.097 hours/ns, 12.571 timesteps/s
100.1% CPU use with 1 MPI tasks x 1 OpenMP threads

MPI task timing breakdown:
Section |  min time  |  avg time  |  max time  |%varavg| %total
---------------------------------------------------------------
Pair    | 78.932     | 78.932     | 78.932     |   0.0 | 99.22
Neigh   | 0.15138    | 0.15138    | 0.15138    |   0.0 |  0.19
Comm    | 0.077514   | 0.077514   | 0.077514   |   0.0 |  0.10
Output  | 0.00020003 | 0.00020003 | 0.00020003 |   0.0 |  0.00
Modify  | 0.33338    | 0.33338    | 0.33338    |   0.0 |  0.42
Other   |            | 0.05443    |            |       |  0.07

Nlocal:           2000 ave        2000 max        2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:           2910 ave        2910 max        2910 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Neighs:              0 ave           0 max           0 min
Histogram: 1 0 0 0 0 0 0 0 0 0
FullNghs:        62550 ave       62550 max       62550 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 62550
Ave neighs/atom = 31.275
Neighbor list builds = 11
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
Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule
Neighbor list info ...
  update every 1 steps, delay 10 steps, check yes
  max neighbors/atom: 2000, page size: 100000
  master list distance cutoff = 6.8
  ghost atom cutoff = 6.8
  binsize = 3.4, bins = 13 13 13
  1 neighbor lists, perpetual/occasional/extra = 1 0 0
  (1) pair kim, perpetual
      attributes: full, newton off, cut 6.8
      pair build: full/bin/atomonly
      stencil: full/bin/3d
      bin: standard
Per MPI rank memory allocation (min/avg/max) = 3.325 | 3.325 | 3.325 Mbytes
   Step         TotEng     v_etotal_metal     PotEng       v_pe_metal        Temp        v_T_metal        Volume       v_V_metal        Press        v_P_metal   
      2000  -2091.894      -2091.894      -2157.4622     -2157.4622      253.75571      253.75571      76726.53       76726.53      -88.441783     -88.441783    
      3000  -2091.3569     -2091.3569     -2157.1199     -2157.1199      254.50979      254.50979      76674.969      76674.969      120.9712       120.9712     
Loop time of 79.5847 on 1 procs for 1000 steps with 2000 atoms

Performance: 1.086 ns/day, 22.107 hours/ns, 12.565 timesteps/s
100.1% CPU use with 1 MPI tasks x 1 OpenMP threads

MPI task timing breakdown:
Section |  min time  |  avg time  |  max time  |%varavg| %total
---------------------------------------------------------------
Pair    | 78.967     | 78.967     | 78.967     |   0.0 | 99.22
Neigh   | 0.15166    | 0.15166    | 0.15166    |   0.0 |  0.19
Comm    | 0.077945   | 0.077945   | 0.077945   |   0.0 |  0.10
Output  | 0.00019988 | 0.00019988 | 0.00019988 |   0.0 |  0.00
Modify  | 0.33339    | 0.33339    | 0.33339    |   0.0 |  0.42
Other   |            | 0.05456    |            |       |  0.07

Nlocal:           2000 ave        2000 max        2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:           2881 ave        2881 max        2881 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Neighs:              0 ave           0 max           0 min
Histogram: 1 0 0 0 0 0 0 0 0 0
FullNghs:        62658 ave       62658 max       62658 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 62658
Ave neighs/atom = 31.329
Neighbor list builds = 11
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
Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule
Neighbor list info ...
  update every 1 steps, delay 10 steps, check yes
  max neighbors/atom: 2000, page size: 100000
  master list distance cutoff = 6.8
  ghost atom cutoff = 6.8
  binsize = 3.4, bins = 13 13 13
  1 neighbor lists, perpetual/occasional/extra = 1 0 0
  (1) pair kim, perpetual
      attributes: full, newton off, cut 6.8
      pair build: full/bin/atomonly
      stencil: full/bin/3d
      bin: standard
Per MPI rank memory allocation (min/avg/max) = 3.325 | 3.325 | 3.325 Mbytes
   Step         TotEng     v_etotal_metal     PotEng       v_pe_metal        Temp        v_T_metal        Volume       v_V_metal        Press        v_P_metal   
      3000  -2091.3569     -2091.3569     -2157.1199     -2157.1199      254.50979      254.50979      76674.969      76674.969      120.9712       120.9712     
      4000  -2091.1035     -2091.1035     -2155.71       -2155.71        250.03365      250.03365      76732.593      76732.593      8.621149       8.621149     
Loop time of 79.6373 on 1 procs for 1000 steps with 2000 atoms

Performance: 1.085 ns/day, 22.121 hours/ns, 12.557 timesteps/s
100.1% CPU use with 1 MPI tasks x 1 OpenMP threads

MPI task timing breakdown:
Section |  min time  |  avg time  |  max time  |%varavg| %total
---------------------------------------------------------------
Pair    | 79.02      | 79.02      | 79.02      |   0.0 | 99.22
Neigh   | 0.15161    | 0.15161    | 0.15161    |   0.0 |  0.19
Comm    | 0.077459   | 0.077459   | 0.077459   |   0.0 |  0.10
Output  | 0.00020041 | 0.00020041 | 0.00020041 |   0.0 |  0.00
Modify  | 0.33404    | 0.33404    | 0.33404    |   0.0 |  0.42
Other   |            | 0.0543     |            |       |  0.07

Nlocal:           2000 ave        2000 max        2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:           2890 ave        2890 max        2890 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Neighs:              0 ave           0 max           0 min
Histogram: 1 0 0 0 0 0 0 0 0 0
FullNghs:        62152 ave       62152 max       62152 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 62152
Ave neighs/atom = 31.076
Neighbor list builds = 11
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
Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule
Neighbor list info ...
  update every 1 steps, delay 10 steps, check yes
  max neighbors/atom: 2000, page size: 100000
  master list distance cutoff = 6.8
  ghost atom cutoff = 6.8
  binsize = 3.4, bins = 13 13 13
  1 neighbor lists, perpetual/occasional/extra = 1 0 0
  (1) pair kim, perpetual
      attributes: full, newton off, cut 6.8
      pair build: full/bin/atomonly
      stencil: full/bin/3d
      bin: standard
Per MPI rank memory allocation (min/avg/max) = 3.325 | 3.325 | 3.325 Mbytes
   Step         TotEng     v_etotal_metal     PotEng       v_pe_metal        Temp        v_T_metal        Volume       v_V_metal        Press        v_P_metal   
      4000  -2091.1035     -2091.1035     -2155.71       -2155.71        250.03365      250.03365      76732.593      76732.593      8.621149       8.621149     
      5000  -2090.0365     -2090.0365     -2155.2527     -2155.2527      252.39338      252.39338      76883.778      76883.778     -133.75085     -133.75085    
Loop time of 81.2309 on 1 procs for 1000 steps with 2000 atoms

Performance: 1.064 ns/day, 22.564 hours/ns, 12.311 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    | 80.584     | 80.584     | 80.584     |   0.0 | 99.20
Neigh   | 0.16715    | 0.16715    | 0.16715    |   0.0 |  0.21
Comm    | 0.079121   | 0.079121   | 0.079121   |   0.0 |  0.10
Output  | 0.00021312 | 0.00021312 | 0.00021312 |   0.0 |  0.00
Modify  | 0.3452     | 0.3452     | 0.3452     |   0.0 |  0.42
Other   |            | 0.05503    |            |       |  0.07

Nlocal:           2000 ave        2000 max        2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:           2880 ave        2880 max        2880 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Neighs:              0 ave           0 max           0 min
Histogram: 1 0 0 0 0 0 0 0 0 0
FullNghs:        62518 ave       62518 max       62518 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 62518
Ave neighs/atom = 31.259
Neighbor list builds = 12
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 252.097065193239, Press = -79.9181505919714
next a
jump SELF top
variable a loop 2000
run 1000
Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule
Neighbor list info ...
  update every 1 steps, delay 10 steps, check yes
  max neighbors/atom: 2000, page size: 100000
  master list distance cutoff = 6.8
  ghost atom cutoff = 6.8
  binsize = 3.4, bins = 13 13 13
  1 neighbor lists, perpetual/occasional/extra = 1 0 0
  (1) pair kim, perpetual
      attributes: full, newton off, cut 6.8
      pair build: full/bin/atomonly
      stencil: full/bin/3d
      bin: standard
Per MPI rank memory allocation (min/avg/max) = 3.325 | 3.325 | 3.325 Mbytes
   Step         TotEng     v_etotal_metal     PotEng       v_pe_metal        Temp        v_T_metal        Volume       v_V_metal        Press        v_P_metal   
      5000  -2090.0365     -2090.0365     -2155.2527     -2155.2527      252.39338      252.39338      76883.778      76883.778     -133.75085     -133.75085    
      6000  -2091.4631     -2091.4631     -2157.5091     -2157.5091      255.60476      255.60476      76442.764      76442.764      256.82924      256.82924    
Loop time of 92.9591 on 1 procs for 1000 steps with 2000 atoms

Performance: 0.929 ns/day, 25.822 hours/ns, 10.757 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    | 92.201     | 92.201     | 92.201     |   0.0 | 99.18
Neigh   | 0.17943    | 0.17943    | 0.17943    |   0.0 |  0.19
Comm    | 0.088223   | 0.088223   | 0.088223   |   0.0 |  0.09
Output  | 0.00021689 | 0.00021689 | 0.00021689 |   0.0 |  0.00
Modify  | 0.43048    | 0.43048    | 0.43048    |   0.0 |  0.46
Other   |            | 0.05933    |            |       |  0.06

Nlocal:           2000 ave        2000 max        2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:           2902 ave        2902 max        2902 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Neighs:              0 ave           0 max           0 min
Histogram: 1 0 0 0 0 0 0 0 0 0
FullNghs:        62918 ave       62918 max       62918 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 62918
Ave neighs/atom = 31.459
Neighbor list builds = 12
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 252.606960305817, Press = -8.63613127301441
next a
jump SELF top
variable a loop 2000
run 1000
Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule
Neighbor list info ...
  update every 1 steps, delay 10 steps, check yes
  max neighbors/atom: 2000, page size: 100000
  master list distance cutoff = 6.8
  ghost atom cutoff = 6.8
  binsize = 3.4, bins = 13 13 13
  1 neighbor lists, perpetual/occasional/extra = 1 0 0
  (1) pair kim, perpetual
      attributes: full, newton off, cut 6.8
      pair build: full/bin/atomonly
      stencil: full/bin/3d
      bin: standard
Per MPI rank memory allocation (min/avg/max) = 3.325 | 3.325 | 3.325 Mbytes
   Step         TotEng     v_etotal_metal     PotEng       v_pe_metal        Temp        v_T_metal        Volume       v_V_metal        Press        v_P_metal   
      6000  -2091.4631     -2091.4631     -2157.5091     -2157.5091      255.60476      255.60476      76442.764      76442.764      256.82924      256.82924    
      7000  -2090.2771     -2090.2771     -2156.3728     -2156.3728      255.79724      255.79724      77071.499      77071.499     -410.03172     -410.03172    
Loop time of 92.902 on 1 procs for 1000 steps with 2000 atoms

Performance: 0.930 ns/day, 25.806 hours/ns, 10.764 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    | 92.139     | 92.139     | 92.139     |   0.0 | 99.18
Neigh   | 0.17853    | 0.17853    | 0.17853    |   0.0 |  0.19
Comm    | 0.089142   | 0.089142   | 0.089142   |   0.0 |  0.10
Output  | 0.00026687 | 0.00026687 | 0.00026687 |   0.0 |  0.00
Modify  | 0.43544    | 0.43544    | 0.43544    |   0.0 |  0.47
Other   |            | 0.06008    |            |       |  0.06

Nlocal:           2000 ave        2000 max        2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:           2881 ave        2881 max        2881 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Neighs:              0 ave           0 max           0 min
Histogram: 1 0 0 0 0 0 0 0 0 0
FullNghs:        61952 ave       61952 max       61952 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 61952
Ave neighs/atom = 30.976
Neighbor list builds = 12
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 252.695223580037, Press = -16.1211779390054
next a
jump SELF top
variable a loop 2000
run 1000
Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule
Neighbor list info ...
  update every 1 steps, delay 10 steps, check yes
  max neighbors/atom: 2000, page size: 100000
  master list distance cutoff = 6.8
  ghost atom cutoff = 6.8
  binsize = 3.4, bins = 13 13 13
  1 neighbor lists, perpetual/occasional/extra = 1 0 0
  (1) pair kim, perpetual
      attributes: full, newton off, cut 6.8
      pair build: full/bin/atomonly
      stencil: full/bin/3d
      bin: standard
Per MPI rank memory allocation (min/avg/max) = 3.325 | 3.325 | 3.325 Mbytes
   Step         TotEng     v_etotal_metal     PotEng       v_pe_metal        Temp        v_T_metal        Volume       v_V_metal        Press        v_P_metal   
      7000  -2090.2771     -2090.2771     -2156.3728     -2156.3728      255.79724      255.79724      77071.499      77071.499     -410.03172     -410.03172    
      8000  -2088.4589     -2088.4589     -2153.6659     -2153.6659      252.35803      252.35803      76871.358      76871.358      62.83069       62.83069     
Loop time of 80.4924 on 1 procs for 1000 steps with 2000 atoms

Performance: 1.073 ns/day, 22.359 hours/ns, 12.424 timesteps/s
100.1% CPU use with 1 MPI tasks x 1 OpenMP threads

MPI task timing breakdown:
Section |  min time  |  avg time  |  max time  |%varavg| %total
---------------------------------------------------------------
Pair    | 79.833     | 79.833     | 79.833     |   0.0 | 99.18
Neigh   | 0.16646    | 0.16646    | 0.16646    |   0.0 |  0.21
Comm    | 0.080035   | 0.080035   | 0.080035   |   0.0 |  0.10
Output  | 0.00016746 | 0.00016746 | 0.00016746 |   0.0 |  0.00
Modify  | 0.35868    | 0.35868    | 0.35868    |   0.0 |  0.45
Other   |            | 0.05387    |            |       |  0.07

Nlocal:           2000 ave        2000 max        2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:           2895 ave        2895 max        2895 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Neighs:              0 ave           0 max           0 min
Histogram: 1 0 0 0 0 0 0 0 0 0
FullNghs:        62706 ave       62706 max       62706 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 62706
Ave neighs/atom = 31.353
Neighbor list builds = 12
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 252.586235230256, Press = 1.84008808390517
next a
jump SELF top
variable a loop 2000
run 1000
Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule
Neighbor list info ...
  update every 1 steps, delay 10 steps, check yes
  max neighbors/atom: 2000, page size: 100000
  master list distance cutoff = 6.8
  ghost atom cutoff = 6.8
  binsize = 3.4, bins = 13 13 13
  1 neighbor lists, perpetual/occasional/extra = 1 0 0
  (1) pair kim, perpetual
      attributes: full, newton off, cut 6.8
      pair build: full/bin/atomonly
      stencil: full/bin/3d
      bin: standard
Per MPI rank memory allocation (min/avg/max) = 3.325 | 3.325 | 3.325 Mbytes
   Step         TotEng     v_etotal_metal     PotEng       v_pe_metal        Temp        v_T_metal        Volume       v_V_metal        Press        v_P_metal   
      8000  -2088.4589     -2088.4589     -2153.6659     -2153.6659      252.35803      252.35803      76871.358      76871.358      62.83069       62.83069     
      9000  -2090.9839     -2090.9839     -2156.7846     -2156.7846      254.65576      254.65576      76561.131      76561.131      177.34429      177.34429    
Loop time of 80.1121 on 1 procs for 1000 steps with 2000 atoms

Performance: 1.078 ns/day, 22.253 hours/ns, 12.483 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    | 79.457     | 79.457     | 79.457     |   0.0 | 99.18
Neigh   | 0.16613    | 0.16613    | 0.16613    |   0.0 |  0.21
Comm    | 0.079996   | 0.079996   | 0.079996   |   0.0 |  0.10
Output  | 0.00016917 | 0.00016917 | 0.00016917 |   0.0 |  0.00
Modify  | 0.35469    | 0.35469    | 0.35469    |   0.0 |  0.44
Other   |            | 0.0538     |            |       |  0.07

Nlocal:           2000 ave        2000 max        2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:           2880 ave        2880 max        2880 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Neighs:              0 ave           0 max           0 min
Histogram: 1 0 0 0 0 0 0 0 0 0
FullNghs:        62870 ave       62870 max       62870 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 62870
Ave neighs/atom = 31.435
Neighbor list builds = 12
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 252.474279271227, Press = -1.1448609357623
next a
jump SELF top
variable a loop 2000
run 1000
Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule
Neighbor list info ...
  update every 1 steps, delay 10 steps, check yes
  max neighbors/atom: 2000, page size: 100000
  master list distance cutoff = 6.8
  ghost atom cutoff = 6.8
  binsize = 3.4, bins = 13 13 13
  1 neighbor lists, perpetual/occasional/extra = 1 0 0
  (1) pair kim, perpetual
      attributes: full, newton off, cut 6.8
      pair build: full/bin/atomonly
      stencil: full/bin/3d
      bin: standard
Per MPI rank memory allocation (min/avg/max) = 3.325 | 3.325 | 3.325 Mbytes
   Step         TotEng     v_etotal_metal     PotEng       v_pe_metal        Temp        v_T_metal        Volume       v_V_metal        Press        v_P_metal   
      9000  -2090.9839     -2090.9839     -2156.7846     -2156.7846      254.65576      254.65576      76561.131      76561.131      177.34429      177.34429    
     10000  -2091.1025     -2091.1025     -2155.4673     -2155.4673      249.09813      249.09813      76770.747      76770.747     -81.59001      -81.59001     
Loop time of 80.026 on 1 procs for 1000 steps with 2000 atoms

Performance: 1.080 ns/day, 22.229 hours/ns, 12.496 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    | 79.385     | 79.385     | 79.385     |   0.0 | 99.20
Neigh   | 0.15214    | 0.15214    | 0.15214    |   0.0 |  0.19
Comm    | 0.079243   | 0.079243   | 0.079243   |   0.0 |  0.10
Output  | 0.00017052 | 0.00017052 | 0.00017052 |   0.0 |  0.00
Modify  | 0.35523    | 0.35523    | 0.35523    |   0.0 |  0.44
Other   |            | 0.05425    |            |       |  0.07

Nlocal:           2000 ave        2000 max        2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:           2884 ave        2884 max        2884 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Neighs:              0 ave           0 max           0 min
Histogram: 1 0 0 0 0 0 0 0 0 0
FullNghs:        62572 ave       62572 max       62572 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 62572
Ave neighs/atom = 31.286
Neighbor list builds = 11
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 252.339126408344, Press = -3.27499840632134
next a
jump SELF top
variable a loop 2000
run 1000
Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule
Neighbor list info ...
  update every 1 steps, delay 10 steps, check yes
  max neighbors/atom: 2000, page size: 100000
  master list distance cutoff = 6.8
  ghost atom cutoff = 6.8
  binsize = 3.4, bins = 13 13 13
  1 neighbor lists, perpetual/occasional/extra = 1 0 0
  (1) pair kim, perpetual
      attributes: full, newton off, cut 6.8
      pair build: full/bin/atomonly
      stencil: full/bin/3d
      bin: standard
Per MPI rank memory allocation (min/avg/max) = 3.325 | 3.325 | 3.325 Mbytes
   Step         TotEng     v_etotal_metal     PotEng       v_pe_metal        Temp        v_T_metal        Volume       v_V_metal        Press        v_P_metal   
     10000  -2091.1025     -2091.1025     -2155.4673     -2155.4673      249.09813      249.09813      76770.747      76770.747     -81.59001      -81.59001     
     11000  -2087.3335     -2087.3335     -2153.5781     -2153.5781      256.37349      256.37349      76927.132      76927.132     -178.70415     -178.70415    
Loop time of 80.2096 on 1 procs for 1000 steps with 2000 atoms

Performance: 1.077 ns/day, 22.280 hours/ns, 12.467 timesteps/s
100.1% CPU use with 1 MPI tasks x 1 OpenMP threads

MPI task timing breakdown:
Section |  min time  |  avg time  |  max time  |%varavg| %total
---------------------------------------------------------------
Pair    | 79.554     | 79.554     | 79.554     |   0.0 | 99.18
Neigh   | 0.16596    | 0.16596    | 0.16596    |   0.0 |  0.21
Comm    | 0.080082   | 0.080082   | 0.080082   |   0.0 |  0.10
Output  | 0.00016713 | 0.00016713 | 0.00016713 |   0.0 |  0.00
Modify  | 0.35524    | 0.35524    | 0.35524    |   0.0 |  0.44
Other   |            | 0.05414    |            |       |  0.07

Nlocal:           2000 ave        2000 max        2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:           2873 ave        2873 max        2873 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Neighs:              0 ave           0 max           0 min
Histogram: 1 0 0 0 0 0 0 0 0 0
FullNghs:        62440 ave       62440 max       62440 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 62440
Ave neighs/atom = 31.22
Neighbor list builds = 12
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 252.504047243545, Press = -0.983172080650373
next a
jump SELF top
variable a loop 2000
run 1000
Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule
Neighbor list info ...
  update every 1 steps, delay 10 steps, check yes
  max neighbors/atom: 2000, page size: 100000
  master list distance cutoff = 6.8
  ghost atom cutoff = 6.8
  binsize = 3.4, bins = 13 13 13
  1 neighbor lists, perpetual/occasional/extra = 1 0 0
  (1) pair kim, perpetual
      attributes: full, newton off, cut 6.8
      pair build: full/bin/atomonly
      stencil: full/bin/3d
      bin: standard
Per MPI rank memory allocation (min/avg/max) = 3.325 | 3.325 | 3.325 Mbytes
   Step         TotEng     v_etotal_metal     PotEng       v_pe_metal        Temp        v_T_metal        Volume       v_V_metal        Press        v_P_metal   
     11000  -2087.3335     -2087.3335     -2153.5781     -2153.5781      256.37349      256.37349      76927.132      76927.132     -178.70415     -178.70415    
     12000  -2090.738      -2090.738      -2155.5001     -2155.5001      250.63604      250.63604      76739.224      76739.224      29.855674      29.855674    
Loop time of 80.015 on 1 procs for 1000 steps with 2000 atoms

Performance: 1.080 ns/day, 22.226 hours/ns, 12.498 timesteps/s
100.1% CPU use with 1 MPI tasks x 1 OpenMP threads

MPI task timing breakdown:
Section |  min time  |  avg time  |  max time  |%varavg| %total
---------------------------------------------------------------
Pair    | 79.358     | 79.358     | 79.358     |   0.0 | 99.18
Neigh   | 0.16603    | 0.16603    | 0.16603    |   0.0 |  0.21
Comm    | 0.080249   | 0.080249   | 0.080249   |   0.0 |  0.10
Output  | 0.00016865 | 0.00016865 | 0.00016865 |   0.0 |  0.00
Modify  | 0.35637    | 0.35637    | 0.35637    |   0.0 |  0.45
Other   |            | 0.05404    |            |       |  0.07

Nlocal:           2000 ave        2000 max        2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:           2891 ave        2891 max        2891 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Neighs:              0 ave           0 max           0 min
Histogram: 1 0 0 0 0 0 0 0 0 0
FullNghs:        62530 ave       62530 max       62530 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 62530
Ave neighs/atom = 31.265
Neighbor list builds = 12
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 252.70374315876, Press = -2.00466473039047
next a
jump SELF top
variable a loop 2000
run 1000
Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule
Neighbor list info ...
  update every 1 steps, delay 10 steps, check yes
  max neighbors/atom: 2000, page size: 100000
  master list distance cutoff = 6.8
  ghost atom cutoff = 6.8
  binsize = 3.4, bins = 13 13 13
  1 neighbor lists, perpetual/occasional/extra = 1 0 0
  (1) pair kim, perpetual
      attributes: full, newton off, cut 6.8
      pair build: full/bin/atomonly
      stencil: full/bin/3d
      bin: standard
Per MPI rank memory allocation (min/avg/max) = 3.325 | 3.325 | 3.325 Mbytes
   Step         TotEng     v_etotal_metal     PotEng       v_pe_metal        Temp        v_T_metal        Volume       v_V_metal        Press        v_P_metal   
     12000  -2090.738      -2090.738      -2155.5001     -2155.5001      250.63604      250.63604      76739.224      76739.224      29.855674      29.855674    
     13000  -2089.4871     -2089.4871     -2154.9031     -2154.9031      253.16699      253.16699      76781.58       76781.58      -113.52474     -113.52474    
Loop time of 79.9667 on 1 procs for 1000 steps with 2000 atoms

Performance: 1.080 ns/day, 22.213 hours/ns, 12.505 timesteps/s
100.1% CPU use with 1 MPI tasks x 1 OpenMP threads

MPI task timing breakdown:
Section |  min time  |  avg time  |  max time  |%varavg| %total
---------------------------------------------------------------
Pair    | 79.324     | 79.324     | 79.324     |   0.0 | 99.20
Neigh   | 0.153      | 0.153      | 0.153      |   0.0 |  0.19
Comm    | 0.079776   | 0.079776   | 0.079776   |   0.0 |  0.10
Output  | 0.00016971 | 0.00016971 | 0.00016971 |   0.0 |  0.00
Modify  | 0.35595    | 0.35595    | 0.35595    |   0.0 |  0.45
Other   |            | 0.05379    |            |       |  0.07

Nlocal:           2000 ave        2000 max        2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:           2900 ave        2900 max        2900 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Neighs:              0 ave           0 max           0 min
Histogram: 1 0 0 0 0 0 0 0 0 0
FullNghs:        62662 ave       62662 max       62662 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 62662
Ave neighs/atom = 31.331
Neighbor list builds = 11
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 252.688201465621, Press = -1.72541702320006
next a
jump SELF top
variable a loop 2000
run 1000
Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule
Neighbor list info ...
  update every 1 steps, delay 10 steps, check yes
  max neighbors/atom: 2000, page size: 100000
  master list distance cutoff = 6.8
  ghost atom cutoff = 6.8
  binsize = 3.4, bins = 13 13 13
  1 neighbor lists, perpetual/occasional/extra = 1 0 0
  (1) pair kim, perpetual
      attributes: full, newton off, cut 6.8
      pair build: full/bin/atomonly
      stencil: full/bin/3d
      bin: standard
Per MPI rank memory allocation (min/avg/max) = 3.325 | 3.325 | 3.325 Mbytes
   Step         TotEng     v_etotal_metal     PotEng       v_pe_metal        Temp        v_T_metal        Volume       v_V_metal        Press        v_P_metal   
     13000  -2089.4871     -2089.4871     -2154.9031     -2154.9031      253.16699      253.16699      76781.58       76781.58      -113.52474     -113.52474    
     14000  -2088.9672     -2088.9672     -2154.5977     -2154.5977      253.99687      253.99687      76742.915      76742.915      29.864722      29.864722    
Loop time of 88.1564 on 1 procs for 1000 steps with 2000 atoms

Performance: 0.980 ns/day, 24.488 hours/ns, 11.343 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    | 87.452     | 87.452     | 87.452     |   0.0 | 99.20
Neigh   | 0.15979    | 0.15979    | 0.15979    |   0.0 |  0.18
Comm    | 0.085232   | 0.085232   | 0.085232   |   0.0 |  0.10
Output  | 0.00021336 | 0.00021336 | 0.00021336 |   0.0 |  0.00
Modify  | 0.40132    | 0.40132    | 0.40132    |   0.0 |  0.46
Other   |            | 0.05833    |            |       |  0.07

Nlocal:           2000 ave        2000 max        2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:           2892 ave        2892 max        2892 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Neighs:              0 ave           0 max           0 min
Histogram: 1 0 0 0 0 0 0 0 0 0
FullNghs:        62474 ave       62474 max       62474 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 62474
Ave neighs/atom = 31.237
Neighbor list builds = 11
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 252.803999606099, Press = -1.91717721592774
next a
jump SELF top
variable a loop 2000
run 1000
Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule
Neighbor list info ...
  update every 1 steps, delay 10 steps, check yes
  max neighbors/atom: 2000, page size: 100000
  master list distance cutoff = 6.8
  ghost atom cutoff = 6.8
  binsize = 3.4, bins = 13 13 13
  1 neighbor lists, perpetual/occasional/extra = 1 0 0
  (1) pair kim, perpetual
      attributes: full, newton off, cut 6.8
      pair build: full/bin/atomonly
      stencil: full/bin/3d
      bin: standard
Per MPI rank memory allocation (min/avg/max) = 3.325 | 3.325 | 3.325 Mbytes
   Step         TotEng     v_etotal_metal     PotEng       v_pe_metal        Temp        v_T_metal        Volume       v_V_metal        Press        v_P_metal   
     14000  -2088.9672     -2088.9672     -2154.5977     -2154.5977      253.99687      253.99687      76742.915      76742.915      29.864722      29.864722    
     15000  -2090.7599     -2090.7599     -2155.3943     -2155.3943      250.14183      250.14183      76781.476      76781.476     -16.501887     -16.501887    
Loop time of 80.853 on 1 procs for 1000 steps with 2000 atoms

Performance: 1.069 ns/day, 22.459 hours/ns, 12.368 timesteps/s
100.1% CPU use with 1 MPI tasks x 1 OpenMP threads

MPI task timing breakdown:
Section |  min time  |  avg time  |  max time  |%varavg| %total
---------------------------------------------------------------
Pair    | 80.191     | 80.191     | 80.191     |   0.0 | 99.18
Neigh   | 0.16643    | 0.16643    | 0.16643    |   0.0 |  0.21
Comm    | 0.079433   | 0.079433   | 0.079433   |   0.0 |  0.10
Output  | 0.00016849 | 0.00016849 | 0.00016849 |   0.0 |  0.00
Modify  | 0.3606     | 0.3606     | 0.3606     |   0.0 |  0.45
Other   |            | 0.05506    |            |       |  0.07

Nlocal:           2000 ave        2000 max        2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:           2878 ave        2878 max        2878 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Neighs:              0 ave           0 max           0 min
Histogram: 1 0 0 0 0 0 0 0 0 0
FullNghs:        62510 ave       62510 max       62510 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 62510
Ave neighs/atom = 31.255
Neighbor list builds = 12
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 252.925311386175, Press = -1.04379985516776
next a
jump SELF top
variable a loop 2000
run 1000
Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule
Neighbor list info ...
  update every 1 steps, delay 10 steps, check yes
  max neighbors/atom: 2000, page size: 100000
  master list distance cutoff = 6.8
  ghost atom cutoff = 6.8
  binsize = 3.4, bins = 13 13 13
  1 neighbor lists, perpetual/occasional/extra = 1 0 0
  (1) pair kim, perpetual
      attributes: full, newton off, cut 6.8
      pair build: full/bin/atomonly
      stencil: full/bin/3d
      bin: standard
Per MPI rank memory allocation (min/avg/max) = 3.325 | 3.325 | 3.325 Mbytes
   Step         TotEng     v_etotal_metal     PotEng       v_pe_metal        Temp        v_T_metal        Volume       v_V_metal        Press        v_P_metal   
     15000  -2090.7599     -2090.7599     -2155.3943     -2155.3943      250.14183      250.14183      76781.476      76781.476     -16.501887     -16.501887    
     16000  -2088.4408     -2088.4408     -2154.1148     -2154.1148      254.16513      254.16513      76732.822      76732.822      129.72355      129.72355    
Loop time of 79.5851 on 1 procs for 1000 steps with 2000 atoms

Performance: 1.086 ns/day, 22.107 hours/ns, 12.565 timesteps/s
100.1% CPU use with 1 MPI tasks x 1 OpenMP threads

MPI task timing breakdown:
Section |  min time  |  avg time  |  max time  |%varavg| %total
---------------------------------------------------------------
Pair    | 78.945     | 78.945     | 78.945     |   0.0 | 99.20
Neigh   | 0.15129    | 0.15129    | 0.15129    |   0.0 |  0.19
Comm    | 0.078113   | 0.078113   | 0.078113   |   0.0 |  0.10
Output  | 0.00016753 | 0.00016753 | 0.00016753 |   0.0 |  0.00
Modify  | 0.35518    | 0.35518    | 0.35518    |   0.0 |  0.45
Other   |            | 0.05529    |            |       |  0.07

Nlocal:           2000 ave        2000 max        2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:           2873 ave        2873 max        2873 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Neighs:              0 ave           0 max           0 min
Histogram: 1 0 0 0 0 0 0 0 0 0
FullNghs:        62680 ave       62680 max       62680 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 62680
Ave neighs/atom = 31.34
Neighbor list builds = 11
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 252.826583000049, Press = -2.44531685267877
next a
jump SELF top
variable a loop 2000
run 1000
Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule
Neighbor list info ...
  update every 1 steps, delay 10 steps, check yes
  max neighbors/atom: 2000, page size: 100000
  master list distance cutoff = 6.8
  ghost atom cutoff = 6.8
  binsize = 3.4, bins = 13 13 13
  1 neighbor lists, perpetual/occasional/extra = 1 0 0
  (1) pair kim, perpetual
      attributes: full, newton off, cut 6.8
      pair build: full/bin/atomonly
      stencil: full/bin/3d
      bin: standard
Per MPI rank memory allocation (min/avg/max) = 3.325 | 3.325 | 3.325 Mbytes
   Step         TotEng     v_etotal_metal     PotEng       v_pe_metal        Temp        v_T_metal        Volume       v_V_metal        Press        v_P_metal   
     16000  -2088.4408     -2088.4408     -2154.1148     -2154.1148      254.16513      254.16513      76732.822      76732.822      129.72355      129.72355    
     17000  -2089.8906     -2089.8906     -2155.3622     -2155.3622      253.38179      253.38179      76961.899      76961.899     -158.10733     -158.10733    
Loop time of 79.6101 on 1 procs for 1000 steps with 2000 atoms

Performance: 1.085 ns/day, 22.114 hours/ns, 12.561 timesteps/s
100.1% CPU use with 1 MPI tasks x 1 OpenMP threads

MPI task timing breakdown:
Section |  min time  |  avg time  |  max time  |%varavg| %total
---------------------------------------------------------------
Pair    | 78.957     | 78.957     | 78.957     |   0.0 | 99.18
Neigh   | 0.16515    | 0.16515    | 0.16515    |   0.0 |  0.21
Comm    | 0.078161   | 0.078161   | 0.078161   |   0.0 |  0.10
Output  | 0.00019829 | 0.00019829 | 0.00019829 |   0.0 |  0.00
Modify  | 0.35481    | 0.35481    | 0.35481    |   0.0 |  0.45
Other   |            | 0.05495    |            |       |  0.07

Nlocal:           2000 ave        2000 max        2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:           2873 ave        2873 max        2873 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Neighs:              0 ave           0 max           0 min
Histogram: 1 0 0 0 0 0 0 0 0 0
FullNghs:        62590 ave       62590 max       62590 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 62590
Ave neighs/atom = 31.295
Neighbor list builds = 12
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 252.857899150285, Press = -0.61532259674229
next a
jump SELF top
variable a loop 2000
run 1000
Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule
Neighbor list info ...
  update every 1 steps, delay 10 steps, check yes
  max neighbors/atom: 2000, page size: 100000
  master list distance cutoff = 6.8
  ghost atom cutoff = 6.8
  binsize = 3.4, bins = 13 13 13
  1 neighbor lists, perpetual/occasional/extra = 1 0 0
  (1) pair kim, perpetual
      attributes: full, newton off, cut 6.8
      pair build: full/bin/atomonly
      stencil: full/bin/3d
      bin: standard
Per MPI rank memory allocation (min/avg/max) = 3.325 | 3.325 | 3.325 Mbytes
   Step         TotEng     v_etotal_metal     PotEng       v_pe_metal        Temp        v_T_metal        Volume       v_V_metal        Press        v_P_metal   
     17000  -2089.8906     -2089.8906     -2155.3622     -2155.3622      253.38179      253.38179      76961.899      76961.899     -158.10733     -158.10733    
     18000  -2086.2864     -2086.2864     -2153.1769     -2153.1769      258.8732       258.8732       76662.999      76662.999      287.00819      287.00819    
Loop time of 79.9098 on 1 procs for 1000 steps with 2000 atoms

Performance: 1.081 ns/day, 22.197 hours/ns, 12.514 timesteps/s
100.1% CPU use with 1 MPI tasks x 1 OpenMP threads

MPI task timing breakdown:
Section |  min time  |  avg time  |  max time  |%varavg| %total
---------------------------------------------------------------
Pair    | 79.254     | 79.254     | 79.254     |   0.0 | 99.18
Neigh   | 0.16628    | 0.16628    | 0.16628    |   0.0 |  0.21
Comm    | 0.079465   | 0.079465   | 0.079465   |   0.0 |  0.10
Output  | 0.00016851 | 0.00016851 | 0.00016851 |   0.0 |  0.00
Modify  | 0.35573    | 0.35573    | 0.35573    |   0.0 |  0.45
Other   |            | 0.05458    |            |       |  0.07

Nlocal:           2000 ave        2000 max        2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:           2891 ave        2891 max        2891 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Neighs:              0 ave           0 max           0 min
Histogram: 1 0 0 0 0 0 0 0 0 0
FullNghs:        62640 ave       62640 max       62640 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 62640
Ave neighs/atom = 31.32
Neighbor list builds = 12
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 253.127658180367, Press = -0.438412092896876
next a
jump SELF top
variable a loop 2000
run 1000
Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule
Neighbor list info ...
  update every 1 steps, delay 10 steps, check yes
  max neighbors/atom: 2000, page size: 100000
  master list distance cutoff = 6.8
  ghost atom cutoff = 6.8
  binsize = 3.4, bins = 13 13 13
  1 neighbor lists, perpetual/occasional/extra = 1 0 0
  (1) pair kim, perpetual
      attributes: full, newton off, cut 6.8
      pair build: full/bin/atomonly
      stencil: full/bin/3d
      bin: standard
Per MPI rank memory allocation (min/avg/max) = 3.325 | 3.325 | 3.325 Mbytes
   Step         TotEng     v_etotal_metal     PotEng       v_pe_metal        Temp        v_T_metal        Volume       v_V_metal        Press        v_P_metal   
     18000  -2086.2864     -2086.2864     -2153.1769     -2153.1769      258.8732       258.8732       76662.999      76662.999      287.00819      287.00819    
     19000  -2090.7312     -2090.7312     -2158.369      -2158.369       261.76526      261.76526      76577.841      76577.841      0.35699285     0.35699285   
Loop time of 80.2628 on 1 procs for 1000 steps with 2000 atoms

Performance: 1.076 ns/day, 22.295 hours/ns, 12.459 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    | 79.621     | 79.621     | 79.621     |   0.0 | 99.20
Neigh   | 0.1518     | 0.1518     | 0.1518     |   0.0 |  0.19
Comm    | 0.078817   | 0.078817   | 0.078817   |   0.0 |  0.10
Output  | 0.00016613 | 0.00016613 | 0.00016613 |   0.0 |  0.00
Modify  | 0.35699    | 0.35699    | 0.35699    |   0.0 |  0.44
Other   |            | 0.05415    |            |       |  0.07

Nlocal:           2000 ave        2000 max        2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:           2913 ave        2913 max        2913 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Neighs:              0 ave           0 max           0 min
Histogram: 1 0 0 0 0 0 0 0 0 0
FullNghs:        62448 ave       62448 max       62448 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 62448
Ave neighs/atom = 31.224
Neighbor list builds = 11
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 253.206868064198, Press = -1.22725290025692
next a
jump SELF top
variable a loop 2000
run 1000
Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule
Neighbor list info ...
  update every 1 steps, delay 10 steps, check yes
  max neighbors/atom: 2000, page size: 100000
  master list distance cutoff = 6.8
  ghost atom cutoff = 6.8
  binsize = 3.4, bins = 13 13 13
  1 neighbor lists, perpetual/occasional/extra = 1 0 0
  (1) pair kim, perpetual
      attributes: full, newton off, cut 6.8
      pair build: full/bin/atomonly
      stencil: full/bin/3d
      bin: standard
Per MPI rank memory allocation (min/avg/max) = 3.325 | 3.325 | 3.325 Mbytes
   Step         TotEng     v_etotal_metal     PotEng       v_pe_metal        Temp        v_T_metal        Volume       v_V_metal        Press        v_P_metal   
     19000  -2090.7312     -2090.7312     -2158.369      -2158.369       261.76526      261.76526      76577.841      76577.841      0.35699285     0.35699285   
     20000  -2089.1025     -2089.1025     -2154.7161     -2154.7161      253.93135      253.93135      77116.522      77116.522     -285.36832     -285.36832    
Loop time of 79.6163 on 1 procs for 1000 steps with 2000 atoms

Performance: 1.085 ns/day, 22.116 hours/ns, 12.560 timesteps/s
100.1% CPU use with 1 MPI tasks x 1 OpenMP threads

MPI task timing breakdown:
Section |  min time  |  avg time  |  max time  |%varavg| %total
---------------------------------------------------------------
Pair    | 78.966     | 78.966     | 78.966     |   0.0 | 99.18
Neigh   | 0.1652     | 0.1652     | 0.1652     |   0.0 |  0.21
Comm    | 0.07775    | 0.07775    | 0.07775    |   0.0 |  0.10
Output  | 0.00016686 | 0.00016686 | 0.00016686 |   0.0 |  0.00
Modify  | 0.35272    | 0.35272    | 0.35272    |   0.0 |  0.44
Other   |            | 0.05417    |            |       |  0.07

Nlocal:           2000 ave        2000 max        2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:           2879 ave        2879 max        2879 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Neighs:              0 ave           0 max           0 min
Histogram: 1 0 0 0 0 0 0 0 0 0
FullNghs:        62256 ave       62256 max       62256 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 62256
Ave neighs/atom = 31.128
Neighbor list builds = 12
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 253.300130376512, Press = -1.10974448666456
next a
jump SELF top
variable a loop 2000
run 1000
Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule
Neighbor list info ...
  update every 1 steps, delay 10 steps, check yes
  max neighbors/atom: 2000, page size: 100000
  master list distance cutoff = 6.8
  ghost atom cutoff = 6.8
  binsize = 3.4, bins = 13 13 13
  1 neighbor lists, perpetual/occasional/extra = 1 0 0
  (1) pair kim, perpetual
      attributes: full, newton off, cut 6.8
      pair build: full/bin/atomonly
      stencil: full/bin/3d
      bin: standard
Per MPI rank memory allocation (min/avg/max) = 3.325 | 3.325 | 3.325 Mbytes
   Step         TotEng     v_etotal_metal     PotEng       v_pe_metal        Temp        v_T_metal        Volume       v_V_metal        Press        v_P_metal   
     20000  -2089.1025     -2089.1025     -2154.7161     -2154.7161      253.93135      253.93135      77116.522      77116.522     -285.36832     -285.36832    
     21000  -2090.5746     -2090.5746     -2154.9961     -2154.9961      249.31777      249.31777      76563.093      76563.093      242.5626       242.5626     
Loop time of 79.5565 on 1 procs for 1000 steps with 2000 atoms

Performance: 1.086 ns/day, 22.099 hours/ns, 12.570 timesteps/s
100.1% CPU use with 1 MPI tasks x 1 OpenMP threads

MPI task timing breakdown:
Section |  min time  |  avg time  |  max time  |%varavg| %total
---------------------------------------------------------------
Pair    | 78.905     | 78.905     | 78.905     |   0.0 | 99.18
Neigh   | 0.16524    | 0.16524    | 0.16524    |   0.0 |  0.21
Comm    | 0.078241   | 0.078241   | 0.078241   |   0.0 |  0.10
Output  | 0.00016757 | 0.00016757 | 0.00016757 |   0.0 |  0.00
Modify  | 0.35365    | 0.35365    | 0.35365    |   0.0 |  0.44
Other   |            | 0.05449    |            |       |  0.07

Nlocal:           2000 ave        2000 max        2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:           2908 ave        2908 max        2908 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Neighs:              0 ave           0 max           0 min
Histogram: 1 0 0 0 0 0 0 0 0 0
FullNghs:        62812 ave       62812 max       62812 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 62812
Ave neighs/atom = 31.406
Neighbor list builds = 12
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 253.257668126324, Press = -0.253742862526286
next a
jump SELF top
variable a loop 2000
run 1000
Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule
Neighbor list info ...
  update every 1 steps, delay 10 steps, check yes
  max neighbors/atom: 2000, page size: 100000
  master list distance cutoff = 6.8
  ghost atom cutoff = 6.8
  binsize = 3.4, bins = 13 13 13
  1 neighbor lists, perpetual/occasional/extra = 1 0 0
  (1) pair kim, perpetual
      attributes: full, newton off, cut 6.8
      pair build: full/bin/atomonly
      stencil: full/bin/3d
      bin: standard
Per MPI rank memory allocation (min/avg/max) = 3.326 | 3.326 | 3.326 Mbytes
   Step         TotEng     v_etotal_metal     PotEng       v_pe_metal        Temp        v_T_metal        Volume       v_V_metal        Press        v_P_metal   
     21000  -2090.5746     -2090.5746     -2154.9961     -2154.9961      249.31777      249.31777      76563.093      76563.093      242.5626       242.5626     
     22000  -2092.0455     -2092.0455     -2157.3256     -2157.3256      252.64071      252.64071      76520.911      76520.911      104.05708      104.05708    
Loop time of 79.6255 on 1 procs for 1000 steps with 2000 atoms

Performance: 1.085 ns/day, 22.118 hours/ns, 12.559 timesteps/s
100.1% CPU use with 1 MPI tasks x 1 OpenMP threads

MPI task timing breakdown:
Section |  min time  |  avg time  |  max time  |%varavg| %total
---------------------------------------------------------------
Pair    | 78.975     | 78.975     | 78.975     |   0.0 | 99.18
Neigh   | 0.16491    | 0.16491    | 0.16491    |   0.0 |  0.21
Comm    | 0.078073   | 0.078073   | 0.078073   |   0.0 |  0.10
Output  | 0.00016745 | 0.00016745 | 0.00016745 |   0.0 |  0.00
Modify  | 0.35317    | 0.35317    | 0.35317    |   0.0 |  0.44
Other   |            | 0.05392    |            |       |  0.07

Nlocal:           2000 ave        2000 max        2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:           2908 ave        2908 max        2908 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Neighs:              0 ave           0 max           0 min
Histogram: 1 0 0 0 0 0 0 0 0 0
FullNghs:        62348 ave       62348 max       62348 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 62348
Ave neighs/atom = 31.174
Neighbor list builds = 12
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 253.272035599286, Press = -2.07196655088794
next a
jump SELF top
variable a loop 2000
run 1000
Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule
Neighbor list info ...
  update every 1 steps, delay 10 steps, check yes
  max neighbors/atom: 2000, page size: 100000
  master list distance cutoff = 6.8
  ghost atom cutoff = 6.8
  binsize = 3.4, bins = 13 13 13
  1 neighbor lists, perpetual/occasional/extra = 1 0 0
  (1) pair kim, perpetual
      attributes: full, newton off, cut 6.8
      pair build: full/bin/atomonly
      stencil: full/bin/3d
      bin: standard
Per MPI rank memory allocation (min/avg/max) = 3.326 | 3.326 | 3.326 Mbytes
   Step         TotEng     v_etotal_metal     PotEng       v_pe_metal        Temp        v_T_metal        Volume       v_V_metal        Press        v_P_metal   
     22000  -2092.0455     -2092.0455     -2157.3256     -2157.3256      252.64071      252.64071      76520.911      76520.911      104.05708      104.05708    
     23000  -2088.1753     -2088.1753     -2154.0077     -2154.0077      254.77799      254.77799      77035.29       77035.29      -370.17191     -370.17191    
Loop time of 79.9931 on 1 procs for 1000 steps with 2000 atoms

Performance: 1.080 ns/day, 22.220 hours/ns, 12.501 timesteps/s
100.1% CPU use with 1 MPI tasks x 1 OpenMP threads

MPI task timing breakdown:
Section |  min time  |  avg time  |  max time  |%varavg| %total
---------------------------------------------------------------
Pair    | 79.354     | 79.354     | 79.354     |   0.0 | 99.20
Neigh   | 0.1519     | 0.1519     | 0.1519     |   0.0 |  0.19
Comm    | 0.078661   | 0.078661   | 0.078661   |   0.0 |  0.10
Output  | 0.00016731 | 0.00016731 | 0.00016731 |   0.0 |  0.00
Modify  | 0.35466    | 0.35466    | 0.35466    |   0.0 |  0.44
Other   |            | 0.0541     |            |       |  0.07

Nlocal:           2000 ave        2000 max        2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:           2906 ave        2906 max        2906 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Neighs:              0 ave           0 max           0 min
Histogram: 1 0 0 0 0 0 0 0 0 0
FullNghs:        62574 ave       62574 max       62574 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 62574
Ave neighs/atom = 31.287
Neighbor list builds = 11
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 253.203576634184, Press = -0.59651840942918
next a
jump SELF top
variable a loop 2000
run 1000
Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule
Neighbor list info ...
  update every 1 steps, delay 10 steps, check yes
  max neighbors/atom: 2000, page size: 100000
  master list distance cutoff = 6.8
  ghost atom cutoff = 6.8
  binsize = 3.4, bins = 13 13 13
  1 neighbor lists, perpetual/occasional/extra = 1 0 0
  (1) pair kim, perpetual
      attributes: full, newton off, cut 6.8
      pair build: full/bin/atomonly
      stencil: full/bin/3d
      bin: standard
Per MPI rank memory allocation (min/avg/max) = 3.326 | 3.326 | 3.326 Mbytes
   Step         TotEng     v_etotal_metal     PotEng       v_pe_metal        Temp        v_T_metal        Volume       v_V_metal        Press        v_P_metal   
     23000  -2088.1753     -2088.1753     -2154.0077     -2154.0077      254.77799      254.77799      77035.29       77035.29      -370.17191     -370.17191    
     24000  -2087.8864     -2087.8864     -2154.7218     -2154.7218      258.66006      258.66006      76722.651      76722.651      125.01665      125.01665    
Loop time of 79.4378 on 1 procs for 1000 steps with 2000 atoms

Performance: 1.088 ns/day, 22.066 hours/ns, 12.588 timesteps/s
100.1% CPU use with 1 MPI tasks x 1 OpenMP threads

MPI task timing breakdown:
Section |  min time  |  avg time  |  max time  |%varavg| %total
---------------------------------------------------------------
Pair    | 78.802     | 78.802     | 78.802     |   0.0 | 99.20
Neigh   | 0.15154    | 0.15154    | 0.15154    |   0.0 |  0.19
Comm    | 0.07748    | 0.07748    | 0.07748    |   0.0 |  0.10
Output  | 0.0001665  | 0.0001665  | 0.0001665  |   0.0 |  0.00
Modify  | 0.35282    | 0.35282    | 0.35282    |   0.0 |  0.44
Other   |            | 0.05411    |            |       |  0.07

Nlocal:           2000 ave        2000 max        2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:           2879 ave        2879 max        2879 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Neighs:              0 ave           0 max           0 min
Histogram: 1 0 0 0 0 0 0 0 0 0
FullNghs:        62920 ave       62920 max       62920 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 62920
Ave neighs/atom = 31.46
Neighbor list builds = 11
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
jump SELF break

# Write final averaged volume to file if temperature and volume have converged; otherwise wirte a
# flag to indicate non-convergence.
variable      myStep equal step
if "${myStep} < 2000000" then "print '${V}'  file output/vol_T253.15.out"  else  "print 'not_converged'  file output/vol_T253.15.out"
print '${V}'  file output/vol_T253.15.out
76747.7627742931

print "LAMMPS calculation completed"
LAMMPS calculation completed
quit 0