# periodic boundary conditions along all three dimensions
boundary p p p

# Set neighbor skin
variable neigh_skin equal 2.0*${_u_distance}
variable neigh_skin equal 2.0*1
neighbor ${neigh_skin} bin
neighbor 2 bin

# create a supercell with cubic lattice (fcc, bcc, sc, or diamond)
# using 10*10*10 conventional (orthogonal) unit cells
variable latticeconst_converted equal 3.889087237417698*${_u_distance}
variable latticeconst_converted equal 3.889087237417698*1
lattice       fcc  ${latticeconst_converted}
lattice       fcc  3.8890872374177
Lattice spacing in x,y,z = 3.8890872 3.8890872 3.8890872
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 (38.890872 38.890872 38.890872)
  1 by 1 by 1 MPI processor grid
create_atoms  1 box
Created 4000 atoms
  using lattice units in orthogonal box = (0.0000000 0.0000000 0.0000000) to (38.890872 38.890872 38.890872)
  create_atoms CPU = 0.007 seconds

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

kim_interactions Pd

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

kim interactions Pd
#=== BEGIN kim interactions ==================================
pair_style kim MEAM_LAMMPS_KoLee_2013_VPdY__MO_046547823135_001
pair_coeff * * Pd
#=== END kim interactions ====================================


mass          1 ${mass_converted}
mass          1 106.42

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

print "Initial system volume: ${V0_metal} Angstroms^3"
Initial system volume: 58822.4426779163 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_046547823135_001#item-citation

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

Neighbor list info ...
  update every 1 steps, delay 10 steps, check yes
  max neighbors/atom: 2000, page size: 100000
  master list distance cutoff = 8
  ghost atom cutoff = 8
  binsize = 4, bins = 10 10 10
  1 neighbor lists, perpetual/occasional/extra = 1 0 0
  (1) pair kim, perpetual
      attributes: full, newton off, cut 8
      pair build: full/bin/atomonly
      stencil: full/bin/3d
      bin: standard
Per MPI rank memory allocation (min/avg/max) = 5.481 | 5.481 | 5.481 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
       0   -15509.144   -15509.144       -15640       -15640       253.15       253.15    58822.443    58822.443    2376.1662    2376.1662 
    1000   -15378.534   -15378.534   -15507.752   -15507.752    249.98121    249.98121    59512.247    59512.247   -1434.2802   -1434.2802 
Loop time of 436.72 on 1 procs for 1000 steps with 4000 atoms

Performance: 0.198 ns/day, 121.311 hours/ns, 2.290 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    | 435.59     | 435.59     | 435.59     |   0.0 | 99.74
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.22207    | 0.22207    | 0.22207    |   0.0 |  0.05
Output  | 0.00026513 | 0.00026513 | 0.00026513 |   0.0 |  0.00
Modify  | 0.79804    | 0.79804    | 0.79804    |   0.0 |  0.18
Other   |            | 0.1096     |            |       |  0.03

Nlocal:        4000.00 ave        4000 max        4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        8195.00 ave        8195 max        8195 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:      560000.0 ave      560000 max      560000 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 560000
Ave neighs/atom = 140.00000
Neighbor list builds = 0
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 0, Press = 0
next a
jump SELF top
variable a loop 2000
run 1000
Neighbor list info ...
  update every 1 steps, delay 10 steps, check yes
  max neighbors/atom: 2000, page size: 100000
  master list distance cutoff = 8
  ghost atom cutoff = 8
  binsize = 4, bins = 10 10 10
  1 neighbor lists, perpetual/occasional/extra = 1 0 0
  (1) pair kim, perpetual
      attributes: full, newton off, cut 8
      pair build: full/bin/atomonly
      stencil: full/bin/3d
      bin: standard
Per MPI rank memory allocation (min/avg/max) = 5.481 | 5.481 | 5.481 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    1000   -15378.534   -15378.534   -15507.752   -15507.752    249.98121    249.98121    59512.247    59512.247   -1434.2802   -1434.2802 
    2000   -15381.663   -15381.663   -15512.953   -15512.953    253.98887    253.98887    59442.732    59442.732    165.41949    165.41949 
Loop time of 431.828 on 1 procs for 1000 steps with 4000 atoms

Performance: 0.200 ns/day, 119.952 hours/ns, 2.316 timesteps/s
99.9% CPU use with 1 MPI tasks x 1 OpenMP threads

MPI task timing breakdown:
Section |  min time  |  avg time  |  max time  |%varavg| %total
---------------------------------------------------------------
Pair    | 430.76     | 430.76     | 430.76     |   0.0 | 99.75
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.21357    | 0.21357    | 0.21357    |   0.0 |  0.05
Output  | 0.00026863 | 0.00026863 | 0.00026863 |   0.0 |  0.00
Modify  | 0.74465    | 0.74465    | 0.74465    |   0.0 |  0.17
Other   |            | 0.106      |            |       |  0.02

Nlocal:        4000.00 ave        4000 max        4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        8194.00 ave        8194 max        8194 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:      559528.0 ave      559528 max      559528 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 559528
Ave neighs/atom = 139.88200
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 = 8
  ghost atom cutoff = 8
  binsize = 4, bins = 10 10 10
  1 neighbor lists, perpetual/occasional/extra = 1 0 0
  (1) pair kim, perpetual
      attributes: full, newton off, cut 8
      pair build: full/bin/atomonly
      stencil: full/bin/3d
      bin: standard
Per MPI rank memory allocation (min/avg/max) = 5.481 | 5.481 | 5.481 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    2000   -15381.663   -15381.663   -15512.953   -15512.953    253.98887    253.98887    59442.732    59442.732    165.41949    165.41949 
    3000    -15384.52    -15384.52   -15513.108   -15513.108    248.76219    248.76219    59489.607    59489.607   -1289.9546   -1289.9546 
Loop time of 460.47 on 1 procs for 1000 steps with 4000 atoms

Performance: 0.188 ns/day, 127.908 hours/ns, 2.172 timesteps/s
99.8% CPU use with 1 MPI tasks x 1 OpenMP threads

MPI task timing breakdown:
Section |  min time  |  avg time  |  max time  |%varavg| %total
---------------------------------------------------------------
Pair    | 459.33     | 459.33     | 459.33     |   0.0 | 99.75
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.22124    | 0.22124    | 0.22124    |   0.0 |  0.05
Output  | 0.00021591 | 0.00021591 | 0.00021591 |   0.0 |  0.00
Modify  | 0.80587    | 0.80587    | 0.80587    |   0.0 |  0.18
Other   |            | 0.1112     |            |       |  0.02

Nlocal:        4000.00 ave        4000 max        4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        8193.00 ave        8193 max        8193 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:      559728.0 ave      559728 max      559728 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 559728
Ave neighs/atom = 139.93200
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 = 8
  ghost atom cutoff = 8
  binsize = 4, bins = 10 10 10
  1 neighbor lists, perpetual/occasional/extra = 1 0 0
  (1) pair kim, perpetual
      attributes: full, newton off, cut 8
      pair build: full/bin/atomonly
      stencil: full/bin/3d
      bin: standard
Per MPI rank memory allocation (min/avg/max) = 5.481 | 5.481 | 5.481 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    3000    -15384.52    -15384.52   -15513.108   -15513.108    248.76219    248.76219    59489.607    59489.607   -1289.9546   -1289.9546 
    4000   -15378.007   -15378.007   -15509.872   -15509.872    255.10329    255.10329    59500.739    59500.739    -1085.756    -1085.756 
Loop time of 458.677 on 1 procs for 1000 steps with 4000 atoms

Performance: 0.188 ns/day, 127.410 hours/ns, 2.180 timesteps/s
99.8% CPU use with 1 MPI tasks x 1 OpenMP threads

MPI task timing breakdown:
Section |  min time  |  avg time  |  max time  |%varavg| %total
---------------------------------------------------------------
Pair    | 457.55     | 457.55     | 457.55     |   0.0 | 99.75
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.21962    | 0.21962    | 0.21962    |   0.0 |  0.05
Output  | 0.00022071 | 0.00022071 | 0.00022071 |   0.0 |  0.00
Modify  | 0.79824    | 0.79824    | 0.79824    |   0.0 |  0.17
Other   |            | 0.1107     |            |       |  0.02

Nlocal:        4000.00 ave        4000 max        4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        8182.00 ave        8182 max        8182 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:      559686.0 ave      559686 max      559686 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 559686
Ave neighs/atom = 139.92150
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 = 8
  ghost atom cutoff = 8
  binsize = 4, bins = 10 10 10
  1 neighbor lists, perpetual/occasional/extra = 1 0 0
  (1) pair kim, perpetual
      attributes: full, newton off, cut 8
      pair build: full/bin/atomonly
      stencil: full/bin/3d
      bin: standard
Per MPI rank memory allocation (min/avg/max) = 5.481 | 5.481 | 5.481 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    4000   -15378.007   -15378.007   -15509.872   -15509.872    255.10329    255.10329    59500.739    59500.739    -1085.756    -1085.756 
    5000   -15384.743   -15384.743   -15517.017   -15517.017    255.89225    255.89225    59421.115    59421.115    428.84186    428.84186 
Loop time of 452.816 on 1 procs for 1000 steps with 4000 atoms

Performance: 0.191 ns/day, 125.782 hours/ns, 2.208 timesteps/s
99.8% CPU use with 1 MPI tasks x 1 OpenMP threads

MPI task timing breakdown:
Section |  min time  |  avg time  |  max time  |%varavg| %total
---------------------------------------------------------------
Pair    | 451.69     | 451.69     | 451.69     |   0.0 | 99.75
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.21938    | 0.21938    | 0.21938    |   0.0 |  0.05
Output  | 0.00021427 | 0.00021427 | 0.00021427 |   0.0 |  0.00
Modify  | 0.79656    | 0.79656    | 0.79656    |   0.0 |  0.18
Other   |            | 0.1101     |            |       |  0.02

Nlocal:        4000.00 ave        4000 max        4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        8194.00 ave        8194 max        8194 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:      559742.0 ave      559742 max      559742 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 559742
Ave neighs/atom = 139.93550
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 = 246.820192707338, Press = -722.44096775727
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 = 8
  ghost atom cutoff = 8
  binsize = 4, bins = 10 10 10
  1 neighbor lists, perpetual/occasional/extra = 1 0 0
  (1) pair kim, perpetual
      attributes: full, newton off, cut 8
      pair build: full/bin/atomonly
      stencil: full/bin/3d
      bin: standard
Per MPI rank memory allocation (min/avg/max) = 5.481 | 5.481 | 5.481 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    5000   -15384.743   -15384.743   -15517.017   -15517.017    255.89225    255.89225    59421.115    59421.115    428.84186    428.84186 
    6000    -15380.51    -15380.51   -15511.542   -15511.542    253.48905    253.48905    59433.587    59433.587    831.19914    831.19914 
Loop time of 439.305 on 1 procs for 1000 steps with 4000 atoms

Performance: 0.197 ns/day, 122.029 hours/ns, 2.276 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    | 438.2      | 438.2      | 438.2      |   0.0 | 99.75
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.21396    | 0.21396    | 0.21396    |   0.0 |  0.05
Output  | 0.00017746 | 0.00017746 | 0.00017746 |   0.0 |  0.00
Modify  | 0.78712    | 0.78712    | 0.78712    |   0.0 |  0.18
Other   |            | 0.1072     |            |       |  0.02

Nlocal:        4000.00 ave        4000 max        4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        8176.00 ave        8176 max        8176 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:      559764.0 ave      559764 max      559764 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 559764
Ave neighs/atom = 139.94100
Neighbor list builds = 0
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 253.019683424471, Press = -12.2707486507869
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 = 8
  ghost atom cutoff = 8
  binsize = 4, bins = 10 10 10
  1 neighbor lists, perpetual/occasional/extra = 1 0 0
  (1) pair kim, perpetual
      attributes: full, newton off, cut 8
      pair build: full/bin/atomonly
      stencil: full/bin/3d
      bin: standard
Per MPI rank memory allocation (min/avg/max) = 5.481 | 5.481 | 5.481 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    6000    -15380.51    -15380.51   -15511.542   -15511.542    253.48905    253.48905    59433.587    59433.587    831.19914    831.19914 
    7000   -15382.705   -15382.705   -15514.061   -15514.061    254.11861    254.11861    59427.399    59427.399    667.86101    667.86101 
Loop time of 442.92 on 1 procs for 1000 steps with 4000 atoms

Performance: 0.195 ns/day, 123.033 hours/ns, 2.258 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    | 441.8      | 441.8      | 441.8      |   0.0 | 99.75
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.21669    | 0.21669    | 0.21669    |   0.0 |  0.05
Output  | 0.00018037 | 0.00018037 | 0.00018037 |   0.0 |  0.00
Modify  | 0.80032    | 0.80032    | 0.80032    |   0.0 |  0.18
Other   |            | 0.1082     |            |       |  0.02

Nlocal:        4000.00 ave        4000 max        4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        8192.00 ave        8192 max        8192 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:      559692.0 ave      559692 max      559692 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 559692
Ave neighs/atom = 139.92300
Neighbor list builds = 0
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 252.943374184735, Press = -16.0180051297072
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 = 8
  ghost atom cutoff = 8
  binsize = 4, bins = 10 10 10
  1 neighbor lists, perpetual/occasional/extra = 1 0 0
  (1) pair kim, perpetual
      attributes: full, newton off, cut 8
      pair build: full/bin/atomonly
      stencil: full/bin/3d
      bin: standard
Per MPI rank memory allocation (min/avg/max) = 5.481 | 5.481 | 5.481 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    7000   -15382.705   -15382.705   -15514.061   -15514.061    254.11861    254.11861    59427.399    59427.399    667.86101    667.86101 
    8000   -15386.227   -15386.227   -15518.017   -15518.017    254.95651    254.95651    59454.363    59454.363   -680.23233   -680.23233 
Loop time of 397.254 on 1 procs for 1000 steps with 4000 atoms

Performance: 0.217 ns/day, 110.348 hours/ns, 2.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    | 396.26     | 396.26     | 396.26     |   0.0 | 99.75
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.19839    | 0.19839    | 0.19839    |   0.0 |  0.05
Output  | 0.00021227 | 0.00021227 | 0.00021227 |   0.0 |  0.00
Modify  | 0.69181    | 0.69181    | 0.69181    |   0.0 |  0.17
Other   |            | 0.1009     |            |       |  0.03

Nlocal:        4000.00 ave        4000 max        4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        8185.00 ave        8185 max        8185 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:      559776.0 ave      559776 max      559776 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 559776
Ave neighs/atom = 139.94400
Neighbor list builds = 0
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 253.095234546895, Press = -15.6473222516108
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 = 8
  ghost atom cutoff = 8
  binsize = 4, bins = 10 10 10
  1 neighbor lists, perpetual/occasional/extra = 1 0 0
  (1) pair kim, perpetual
      attributes: full, newton off, cut 8
      pair build: full/bin/atomonly
      stencil: full/bin/3d
      bin: standard
Per MPI rank memory allocation (min/avg/max) = 5.481 | 5.481 | 5.481 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    8000   -15386.227   -15386.227   -15518.017   -15518.017    254.95651    254.95651    59454.363    59454.363   -680.23233   -680.23233 
    9000   -15381.042   -15381.042    -15512.62    -15512.62    254.54615    254.54615    59503.143    59503.143   -1604.0377   -1604.0377 
Loop time of 390.55 on 1 procs for 1000 steps with 4000 atoms

Performance: 0.221 ns/day, 108.486 hours/ns, 2.560 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    | 389.58     | 389.58     | 389.58     |   0.0 | 99.75
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.19776    | 0.19776    | 0.19776    |   0.0 |  0.05
Output  | 0.00018179 | 0.00018179 | 0.00018179 |   0.0 |  0.00
Modify  | 0.67498    | 0.67498    | 0.67498    |   0.0 |  0.17
Other   |            | 0.09791    |            |       |  0.03

Nlocal:        4000.00 ave        4000 max        4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        8193.00 ave        8193 max        8193 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:      559732.0 ave      559732 max      559732 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 559732
Ave neighs/atom = 139.93300
Neighbor list builds = 0
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 252.727453444013, Press = -7.69436984496239
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 = 8
  ghost atom cutoff = 8
  binsize = 4, bins = 10 10 10
  1 neighbor lists, perpetual/occasional/extra = 1 0 0
  (1) pair kim, perpetual
      attributes: full, newton off, cut 8
      pair build: full/bin/atomonly
      stencil: full/bin/3d
      bin: standard
Per MPI rank memory allocation (min/avg/max) = 5.481 | 5.481 | 5.481 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    9000   -15381.042   -15381.042    -15512.62    -15512.62    254.54615    254.54615    59503.143    59503.143   -1604.0377   -1604.0377 
   10000   -15385.179   -15385.179   -15513.093   -15513.093    247.45691    247.45691    59476.161    59476.161   -805.37092   -805.37092 
Loop time of 449.624 on 1 procs for 1000 steps with 4000 atoms

Performance: 0.192 ns/day, 124.895 hours/ns, 2.224 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    | 448.48     | 448.48     | 448.48     |   0.0 | 99.75
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.21914    | 0.21914    | 0.21914    |   0.0 |  0.05
Output  | 0.00028579 | 0.00028579 | 0.00028579 |   0.0 |  0.00
Modify  | 0.81337    | 0.81337    | 0.81337    |   0.0 |  0.18
Other   |            | 0.1089     |            |       |  0.02

Nlocal:        4000.00 ave        4000 max        4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        8192.00 ave        8192 max        8192 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:      559648.0 ave      559648 max      559648 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 559648
Ave neighs/atom = 139.91200
Neighbor list builds = 0
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 252.806798068556, Press = -1.41420046177416
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 = 8
  ghost atom cutoff = 8
  binsize = 4, bins = 10 10 10
  1 neighbor lists, perpetual/occasional/extra = 1 0 0
  (1) pair kim, perpetual
      attributes: full, newton off, cut 8
      pair build: full/bin/atomonly
      stencil: full/bin/3d
      bin: standard
Per MPI rank memory allocation (min/avg/max) = 5.481 | 5.481 | 5.481 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   10000   -15385.179   -15385.179   -15513.093   -15513.093    247.45691    247.45691    59476.161    59476.161   -805.37092   -805.37092 
   11000   -15378.093   -15378.093   -15510.986   -15510.986    257.08991    257.08991    59417.965    59417.965    1438.3011    1438.3011 
Loop time of 444.197 on 1 procs for 1000 steps with 4000 atoms

Performance: 0.195 ns/day, 123.388 hours/ns, 2.251 timesteps/s
99.8% CPU use with 1 MPI tasks x 1 OpenMP threads

MPI task timing breakdown:
Section |  min time  |  avg time  |  max time  |%varavg| %total
---------------------------------------------------------------
Pair    | 443.07     | 443.07     | 443.07     |   0.0 | 99.75
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.2186     | 0.2186     | 0.2186     |   0.0 |  0.05
Output  | 0.00029591 | 0.00029591 | 0.00029591 |   0.0 |  0.00
Modify  | 0.79809    | 0.79809    | 0.79809    |   0.0 |  0.18
Other   |            | 0.1087     |            |       |  0.02

Nlocal:        4000.00 ave        4000 max        4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        8194.00 ave        8194 max        8194 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:      559718.0 ave      559718 max      559718 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 559718
Ave neighs/atom = 139.92950
Neighbor list builds = 0
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 252.64993859732, Press = 0.313530999528962
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 = 8
  ghost atom cutoff = 8
  binsize = 4, bins = 10 10 10
  1 neighbor lists, perpetual/occasional/extra = 1 0 0
  (1) pair kim, perpetual
      attributes: full, newton off, cut 8
      pair build: full/bin/atomonly
      stencil: full/bin/3d
      bin: standard
Per MPI rank memory allocation (min/avg/max) = 5.481 | 5.481 | 5.481 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   11000   -15378.093   -15378.093   -15510.986   -15510.986    257.08991    257.08991    59417.965    59417.965    1438.3011    1438.3011 
   12000   -15383.759   -15383.759   -15513.217   -15513.217    250.44451    250.44451    59417.203    59417.203     1063.908     1063.908 
Loop time of 420.89 on 1 procs for 1000 steps with 4000 atoms

Performance: 0.205 ns/day, 116.914 hours/ns, 2.376 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    | 419.83     | 419.83     | 419.83     |   0.0 | 99.75
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.20863    | 0.20863    | 0.20863    |   0.0 |  0.05
Output  | 0.00018294 | 0.00018294 | 0.00018294 |   0.0 |  0.00
Modify  | 0.74454    | 0.74454    | 0.74454    |   0.0 |  0.18
Other   |            | 0.1036     |            |       |  0.02

Nlocal:        4000.00 ave        4000 max        4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        8195.00 ave        8195 max        8195 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:      559702.0 ave      559702 max      559702 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 559702
Ave neighs/atom = 139.92550
Neighbor list builds = 0
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 252.82980930975, Press = -4.70503627227098
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 = 8
  ghost atom cutoff = 8
  binsize = 4, bins = 10 10 10
  1 neighbor lists, perpetual/occasional/extra = 1 0 0
  (1) pair kim, perpetual
      attributes: full, newton off, cut 8
      pair build: full/bin/atomonly
      stencil: full/bin/3d
      bin: standard
Per MPI rank memory allocation (min/avg/max) = 5.481 | 5.481 | 5.481 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   12000   -15383.759   -15383.759   -15513.217   -15513.217    250.44451    250.44451    59417.203    59417.203     1063.908     1063.908 
   13000   -15375.625   -15375.625   -15510.792   -15510.792    261.49102    261.49102    59449.047    59449.047     399.9696     399.9696 
Loop time of 441.632 on 1 procs for 1000 steps with 4000 atoms

Performance: 0.196 ns/day, 122.676 hours/ns, 2.264 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    | 440.52     | 440.52     | 440.52     |   0.0 | 99.75
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.21603    | 0.21603    | 0.21603    |   0.0 |  0.05
Output  | 0.00018353 | 0.00018353 | 0.00018353 |   0.0 |  0.00
Modify  | 0.79248    | 0.79248    | 0.79248    |   0.0 |  0.18
Other   |            | 0.1079     |            |       |  0.02

Nlocal:        4000.00 ave        4000 max        4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        8191.00 ave        8191 max        8191 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:      559716.0 ave      559716 max      559716 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 559716
Ave neighs/atom = 139.92900
Neighbor list builds = 0
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 252.925194548823, Press = -4.15384165741438
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 = 8
  ghost atom cutoff = 8
  binsize = 4, bins = 10 10 10
  1 neighbor lists, perpetual/occasional/extra = 1 0 0
  (1) pair kim, perpetual
      attributes: full, newton off, cut 8
      pair build: full/bin/atomonly
      stencil: full/bin/3d
      bin: standard
Per MPI rank memory allocation (min/avg/max) = 5.481 | 5.481 | 5.481 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   13000   -15375.625   -15375.625   -15510.792   -15510.792    261.49102    261.49102    59449.047    59449.047     399.9696     399.9696 
   14000   -15383.066   -15383.066   -15513.967   -15513.967    253.23741    253.23741    59491.805    59491.805    -1335.646    -1335.646 
Loop time of 398.513 on 1 procs for 1000 steps with 4000 atoms

Performance: 0.217 ns/day, 110.698 hours/ns, 2.509 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    | 397.52     | 397.52     | 397.52     |   0.0 | 99.75
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.20172    | 0.20172    | 0.20172    |   0.0 |  0.05
Output  | 0.00018352 | 0.00018352 | 0.00018352 |   0.0 |  0.00
Modify  | 0.69016    | 0.69016    | 0.69016    |   0.0 |  0.17
Other   |            | 0.09978    |            |       |  0.03

Nlocal:        4000.00 ave        4000 max        4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        8187.00 ave        8187 max        8187 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:      559634.0 ave      559634 max      559634 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 559634
Ave neighs/atom = 139.90850
Neighbor list builds = 0
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 252.989950860512, Press = -5.10745913832736
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 = 8
  ghost atom cutoff = 8
  binsize = 4, bins = 10 10 10
  1 neighbor lists, perpetual/occasional/extra = 1 0 0
  (1) pair kim, perpetual
      attributes: full, newton off, cut 8
      pair build: full/bin/atomonly
      stencil: full/bin/3d
      bin: standard
Per MPI rank memory allocation (min/avg/max) = 5.481 | 5.481 | 5.481 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   14000   -15383.066   -15383.066   -15513.967   -15513.967    253.23741    253.23741    59491.805    59491.805    -1335.646    -1335.646 
   15000   -15384.442   -15384.442   -15512.563   -15512.563    247.85788    247.85788    59475.623    59475.623   -674.59294   -674.59294 
Loop time of 446.872 on 1 procs for 1000 steps with 4000 atoms

Performance: 0.193 ns/day, 124.131 hours/ns, 2.238 timesteps/s
99.8% CPU use with 1 MPI tasks x 1 OpenMP threads

MPI task timing breakdown:
Section |  min time  |  avg time  |  max time  |%varavg| %total
---------------------------------------------------------------
Pair    | 445.74     | 445.74     | 445.74     |   0.0 | 99.75
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.21943    | 0.21943    | 0.21943    |   0.0 |  0.05
Output  | 0.00027992 | 0.00027992 | 0.00027992 |   0.0 |  0.00
Modify  | 0.80348    | 0.80348    | 0.80348    |   0.0 |  0.18
Other   |            | 0.1088     |            |       |  0.02

Nlocal:        4000.00 ave        4000 max        4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        8191.00 ave        8191 max        8191 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:      559728.0 ave      559728 max      559728 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 559728
Ave neighs/atom = 139.93200
Neighbor list builds = 0
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 252.964872742179, Press = -1.31624288361837
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 = 8
  ghost atom cutoff = 8
  binsize = 4, bins = 10 10 10
  1 neighbor lists, perpetual/occasional/extra = 1 0 0
  (1) pair kim, perpetual
      attributes: full, newton off, cut 8
      pair build: full/bin/atomonly
      stencil: full/bin/3d
      bin: standard
Per MPI rank memory allocation (min/avg/max) = 5.481 | 5.481 | 5.481 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   15000   -15384.442   -15384.442   -15512.563   -15512.563    247.85788    247.85788    59475.623    59475.623   -674.59294   -674.59294 
   16000   -15381.131   -15381.131   -15512.953   -15512.953    255.01801    255.01801    59455.509    59455.509   -37.624199   -37.624199 
Loop time of 459.788 on 1 procs for 1000 steps with 4000 atoms

Performance: 0.188 ns/day, 127.719 hours/ns, 2.175 timesteps/s
99.8% CPU use with 1 MPI tasks x 1 OpenMP threads

MPI task timing breakdown:
Section |  min time  |  avg time  |  max time  |%varavg| %total
---------------------------------------------------------------
Pair    | 458.62     | 458.62     | 458.62     |   0.0 | 99.75
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.22257    | 0.22257    | 0.22257    |   0.0 |  0.05
Output  | 0.00022803 | 0.00022803 | 0.00022803 |   0.0 |  0.00
Modify  | 0.83254    | 0.83254    | 0.83254    |   0.0 |  0.18
Other   |            | 0.1104     |            |       |  0.02

Nlocal:        4000.00 ave        4000 max        4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        8186.00 ave        8186 max        8186 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:      559656.0 ave      559656 max      559656 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 559656
Ave neighs/atom = 139.91400
Neighbor list builds = 0
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 252.914718906213, Press = -1.61613673250122
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 = 8
  ghost atom cutoff = 8
  binsize = 4, bins = 10 10 10
  1 neighbor lists, perpetual/occasional/extra = 1 0 0
  (1) pair kim, perpetual
      attributes: full, newton off, cut 8
      pair build: full/bin/atomonly
      stencil: full/bin/3d
      bin: standard
Per MPI rank memory allocation (min/avg/max) = 5.481 | 5.481 | 5.481 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   16000   -15381.131   -15381.131   -15512.953   -15512.953    255.01801    255.01801    59455.509    59455.509   -37.624199   -37.624199 
   17000   -15382.409   -15382.409   -15512.254   -15512.254    251.19275    251.19275    59453.167    59453.167   -98.702452   -98.702452 
Loop time of 453.708 on 1 procs for 1000 steps with 4000 atoms

Performance: 0.190 ns/day, 126.030 hours/ns, 2.204 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    | 452.55     | 452.55     | 452.55     |   0.0 | 99.75
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.22179    | 0.22179    | 0.22179    |   0.0 |  0.05
Output  | 0.00017927 | 0.00017927 | 0.00017927 |   0.0 |  0.00
Modify  | 0.82276    | 0.82276    | 0.82276    |   0.0 |  0.18
Other   |            | 0.109      |            |       |  0.02

Nlocal:        4000.00 ave        4000 max        4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        8188.00 ave        8188 max        8188 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:      559734.0 ave      559734 max      559734 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 559734
Ave neighs/atom = 139.93350
Neighbor list builds = 0
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 252.955746756689, Press = -1.1237843254233
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 = 8
  ghost atom cutoff = 8
  binsize = 4, bins = 10 10 10
  1 neighbor lists, perpetual/occasional/extra = 1 0 0
  (1) pair kim, perpetual
      attributes: full, newton off, cut 8
      pair build: full/bin/atomonly
      stencil: full/bin/3d
      bin: standard
Per MPI rank memory allocation (min/avg/max) = 5.481 | 5.481 | 5.481 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   17000   -15382.409   -15382.409   -15512.254   -15512.254    251.19275    251.19275    59453.167    59453.167   -98.702452   -98.702452 
   18000   -15381.051   -15381.051   -15512.511   -15512.511    254.31828    254.31828     59450.24     59450.24    59.320092    59.320092 
Loop time of 453.881 on 1 procs for 1000 steps with 4000 atoms

Performance: 0.190 ns/day, 126.078 hours/ns, 2.203 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    | 452.73     | 452.73     | 452.73     |   0.0 | 99.75
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.22121    | 0.22121    | 0.22121    |   0.0 |  0.05
Output  | 0.00027046 | 0.00027046 | 0.00027046 |   0.0 |  0.00
Modify  | 0.8239     | 0.8239     | 0.8239     |   0.0 |  0.18
Other   |            | 0.1099     |            |       |  0.02

Nlocal:        4000.00 ave        4000 max        4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        8191.00 ave        8191 max        8191 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:      559732.0 ave      559732 max      559732 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 559732
Ave neighs/atom = 139.93300
Neighbor list builds = 0
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 252.994437596598, Press = -2.88285828165938
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 = 8
  ghost atom cutoff = 8
  binsize = 4, bins = 10 10 10
  1 neighbor lists, perpetual/occasional/extra = 1 0 0
  (1) pair kim, perpetual
      attributes: full, newton off, cut 8
      pair build: full/bin/atomonly
      stencil: full/bin/3d
      bin: standard
Per MPI rank memory allocation (min/avg/max) = 5.481 | 5.481 | 5.481 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   18000   -15381.051   -15381.051   -15512.511   -15512.511    254.31828    254.31828     59450.24     59450.24    59.320092    59.320092 
   19000   -15383.873   -15383.873   -15513.566   -15513.566    250.89879    250.89879    59484.087    59484.087   -1191.9873   -1191.9873 
Loop time of 457.447 on 1 procs for 1000 steps with 4000 atoms

Performance: 0.189 ns/day, 127.068 hours/ns, 2.186 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    | 456.29     | 456.29     | 456.29     |   0.0 | 99.75
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.22076    | 0.22076    | 0.22076    |   0.0 |  0.05
Output  | 0.00017962 | 0.00017962 | 0.00017962 |   0.0 |  0.00
Modify  | 0.83013    | 0.83013    | 0.83013    |   0.0 |  0.18
Other   |            | 0.1094     |            |       |  0.02

Nlocal:        4000.00 ave        4000 max        4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        8193.00 ave        8193 max        8193 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:      559698.0 ave      559698 max      559698 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 559698
Ave neighs/atom = 139.92450
Neighbor list builds = 0
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 252.927889074551, Press = -4.22368875256588
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 = 8
  ghost atom cutoff = 8
  binsize = 4, bins = 10 10 10
  1 neighbor lists, perpetual/occasional/extra = 1 0 0
  (1) pair kim, perpetual
      attributes: full, newton off, cut 8
      pair build: full/bin/atomonly
      stencil: full/bin/3d
      bin: standard
Per MPI rank memory allocation (min/avg/max) = 5.481 | 5.481 | 5.481 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   19000   -15383.873   -15383.873   -15513.566   -15513.566    250.89879    250.89879    59484.087    59484.087   -1191.9873   -1191.9873 
   20000   -15380.722   -15380.722   -15510.313   -15510.313     250.7025     250.7025    59482.597    59482.597   -775.67936   -775.67936 
Loop time of 455.387 on 1 procs for 1000 steps with 4000 atoms

Performance: 0.190 ns/day, 126.496 hours/ns, 2.196 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    | 454.22     | 454.22     | 454.22     |   0.0 | 99.74
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.22187    | 0.22187    | 0.22187    |   0.0 |  0.05
Output  | 0.00018856 | 0.00018856 | 0.00018856 |   0.0 |  0.00
Modify  | 0.83184    | 0.83184    | 0.83184    |   0.0 |  0.18
Other   |            | 0.1118     |            |       |  0.02

Nlocal:        4000.00 ave        4000 max        4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        8191.00 ave        8191 max        8191 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:      559688.0 ave      559688 max      559688 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 559688
Ave neighs/atom = 139.92200
Neighbor list builds = 0
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 252.968899960874, Press = 0.5717218041154
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 = 8
  ghost atom cutoff = 8
  binsize = 4, bins = 10 10 10
  1 neighbor lists, perpetual/occasional/extra = 1 0 0
  (1) pair kim, perpetual
      attributes: full, newton off, cut 8
      pair build: full/bin/atomonly
      stencil: full/bin/3d
      bin: standard
Per MPI rank memory allocation (min/avg/max) = 5.481 | 5.481 | 5.481 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   20000   -15380.722   -15380.722   -15510.313   -15510.313     250.7025     250.7025    59482.597    59482.597   -775.67936   -775.67936 
   21000   -15377.575   -15377.575   -15511.606   -15511.606    259.29189    259.29189    59412.572    59412.572    1500.6141    1500.6141 
Loop time of 456.705 on 1 procs for 1000 steps with 4000 atoms

Performance: 0.189 ns/day, 126.862 hours/ns, 2.190 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    | 455.53     | 455.53     | 455.53     |   0.0 | 99.74
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.22306    | 0.22306    | 0.22306    |   0.0 |  0.05
Output  | 0.00023558 | 0.00023558 | 0.00023558 |   0.0 |  0.00
Modify  | 0.8352     | 0.8352     | 0.8352     |   0.0 |  0.18
Other   |            | 0.1121     |            |       |  0.02

Nlocal:        4000.00 ave        4000 max        4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        8182.00 ave        8182 max        8182 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:      559654.0 ave      559654 max      559654 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 559654
Ave neighs/atom = 139.91350
Neighbor list builds = 0
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 253.032557174423, Press = 0.351450839177104
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 = 8
  ghost atom cutoff = 8
  binsize = 4, bins = 10 10 10
  1 neighbor lists, perpetual/occasional/extra = 1 0 0
  (1) pair kim, perpetual
      attributes: full, newton off, cut 8
      pair build: full/bin/atomonly
      stencil: full/bin/3d
      bin: standard
Per MPI rank memory allocation (min/avg/max) = 5.481 | 5.481 | 5.481 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   21000   -15377.575   -15377.575   -15511.606   -15511.606    259.29189    259.29189    59412.572    59412.572    1500.6141    1500.6141 
   22000   -15382.134   -15382.134   -15512.272   -15512.272    251.76009    251.76009    59441.199    59441.199    404.61743    404.61743 
Loop time of 458.218 on 1 procs for 1000 steps with 4000 atoms

Performance: 0.189 ns/day, 127.283 hours/ns, 2.182 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    | 457.05     | 457.05     | 457.05     |   0.0 | 99.75
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.22165    | 0.22165    | 0.22165    |   0.0 |  0.05
Output  | 0.00023307 | 0.00023307 | 0.00023307 |   0.0 |  0.00
Modify  | 0.83016    | 0.83016    | 0.83016    |   0.0 |  0.18
Other   |            | 0.1108     |            |       |  0.02

Nlocal:        4000.00 ave        4000 max        4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        8190.00 ave        8190 max        8190 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:      559694.0 ave      559694 max      559694 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 559694
Ave neighs/atom = 139.92350
Neighbor list builds = 0
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 253.096673193837, Press = -0.929874000263978
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 = 8
  ghost atom cutoff = 8
  binsize = 4, bins = 10 10 10
  1 neighbor lists, perpetual/occasional/extra = 1 0 0
  (1) pair kim, perpetual
      attributes: full, newton off, cut 8
      pair build: full/bin/atomonly
      stencil: full/bin/3d
      bin: standard
Per MPI rank memory allocation (min/avg/max) = 5.481 | 5.481 | 5.481 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   22000   -15382.134   -15382.134   -15512.272   -15512.272    251.76009    251.76009    59441.199    59441.199    404.61743    404.61743 
   23000   -15381.857   -15381.857   -15513.162   -15513.162    254.01732    254.01732    59439.583    59439.583    391.15951    391.15951 
Loop time of 457.724 on 1 procs for 1000 steps with 4000 atoms

Performance: 0.189 ns/day, 127.146 hours/ns, 2.185 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    | 456.57     | 456.57     | 456.57     |   0.0 | 99.75
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.22038    | 0.22038    | 0.22038    |   0.0 |  0.05
Output  | 0.00023391 | 0.00023391 | 0.00023391 |   0.0 |  0.00
Modify  | 0.82679    | 0.82679    | 0.82679    |   0.0 |  0.18
Other   |            | 0.1105     |            |       |  0.02

Nlocal:        4000.00 ave        4000 max        4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        8190.00 ave        8190 max        8190 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:      559668.0 ave      559668 max      559668 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 559668
Ave neighs/atom = 139.91700
Neighbor list builds = 0
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 253.139353348729, Press = -1.69693420151027
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 = 8
  ghost atom cutoff = 8
  binsize = 4, bins = 10 10 10
  1 neighbor lists, perpetual/occasional/extra = 1 0 0
  (1) pair kim, perpetual
      attributes: full, newton off, cut 8
      pair build: full/bin/atomonly
      stencil: full/bin/3d
      bin: standard
Per MPI rank memory allocation (min/avg/max) = 5.481 | 5.481 | 5.481 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   23000   -15381.857   -15381.857   -15513.162   -15513.162    254.01732    254.01732    59439.583    59439.583    391.15951    391.15951 
   24000   -15384.826   -15384.826   -15514.099   -15514.099    250.08799    250.08799    59462.601    59462.601   -518.55749   -518.55749 
Loop time of 459.12 on 1 procs for 1000 steps with 4000 atoms

Performance: 0.188 ns/day, 127.533 hours/ns, 2.178 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    | 457.95     | 457.95     | 457.95     |   0.0 | 99.75
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.2218     | 0.2218     | 0.2218     |   0.0 |  0.05
Output  | 0.00024689 | 0.00024689 | 0.00024689 |   0.0 |  0.00
Modify  | 0.83336    | 0.83336    | 0.83336    |   0.0 |  0.18
Other   |            | 0.1116     |            |       |  0.02

Nlocal:        4000.00 ave        4000 max        4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        8194.00 ave        8194 max        8194 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:      559658.0 ave      559658 max      559658 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 559658
Ave neighs/atom = 139.91450
Neighbor list builds = 0
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 253.158178994251, Press = -1.97351662310079
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 = 8
  ghost atom cutoff = 8
  binsize = 4, bins = 10 10 10
  1 neighbor lists, perpetual/occasional/extra = 1 0 0
  (1) pair kim, perpetual
      attributes: full, newton off, cut 8
      pair build: full/bin/atomonly
      stencil: full/bin/3d
      bin: standard
Per MPI rank memory allocation (min/avg/max) = 5.481 | 5.481 | 5.481 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   24000   -15384.826   -15384.826   -15514.099   -15514.099    250.08799    250.08799    59462.601    59462.601   -518.55749   -518.55749 
   25000    -15380.45    -15380.45   -15511.141   -15511.141    252.83033    252.83033    59480.712    59480.712   -666.32646   -666.32646 
Loop time of 451.466 on 1 procs for 1000 steps with 4000 atoms

Performance: 0.191 ns/day, 125.407 hours/ns, 2.215 timesteps/s
99.8% CPU use with 1 MPI tasks x 1 OpenMP threads

MPI task timing breakdown:
Section |  min time  |  avg time  |  max time  |%varavg| %total
---------------------------------------------------------------
Pair    | 450.31     | 450.31     | 450.31     |   0.0 | 99.74
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.22088    | 0.22088    | 0.22088    |   0.0 |  0.05
Output  | 0.00022894 | 0.00022894 | 0.00022894 |   0.0 |  0.00
Modify  | 0.81933    | 0.81933    | 0.81933    |   0.0 |  0.18
Other   |            | 0.1114     |            |       |  0.02

Nlocal:        4000.00 ave        4000 max        4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        8191.00 ave        8191 max        8191 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:      559694.0 ave      559694 max      559694 min
Histogram: 1 0 0 0 0 0 0 0 0 0

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

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

print "LAMMPS calculation completed"
LAMMPS calculation completed
quit 0