# 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.0499990060925475*${_u_distance}
variable latticeconst_converted equal 4.0499990060925475*1
lattice       fcc  ${latticeconst_converted}
lattice       fcc  4.04999900609255
Lattice spacing in x,y,z = 4.0499990 4.0499990 4.0499990
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 (40.499990 40.499990 40.499990)
  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 (40.499990 40.499990 40.499990)
  create_atoms CPU = 0.003 seconds

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

kim_interactions Al

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

kim interactions Al
#=== BEGIN kim interactions ==================================
pair_style kim MEAM_LAMMPS_DickelBaskesAslam_2018_MgAlZn__MO_093637366498_001
pair_coeff * * Al
#=== END kim interactions ====================================


mass          1 ${mass_converted}
mass          1 26.981538

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

print "Initial system volume: ${V0_metal} Angstroms^3"
Initial system volume: 66430.0760923112 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_093637366498_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 = 11 11 11
  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.370 | 5.370 | 5.370 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
       0   -13309.144   -13309.144       -13440       -13440       253.15       253.15    66430.076    66430.076    2104.0264    2104.0264 
    1000   -13168.698   -13168.698   -13305.677   -13305.677    264.99516    264.99516    67485.005    67485.005   -775.22965   -775.22965 
Loop time of 737.173 on 1 procs for 1000 steps with 4000 atoms

Performance: 0.117 ns/day, 204.770 hours/ns, 1.357 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    | 736.09     | 736.09     | 736.09     |   0.0 | 99.85
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.1655     | 0.1655     | 0.1655     |   0.0 |  0.02
Output  | 0.00027671 | 0.00027671 | 0.00027671 |   0.0 |  0.00
Modify  | 0.81177    | 0.81177    | 0.81177    |   0.0 |  0.11
Other   |            | 0.1083     |            |       |  0.01

Nlocal:        4000.00 ave        4000 max        4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        5841.00 ave        5841 max        5841 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Neighs:         0.00000 ave           0 max           0 min
Histogram: 1 0 0 0 0 0 0 0 0 0
FullNghs:      536000.0 ave      536000 max      536000 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 536000
Ave neighs/atom = 134.00000
Neighbor list builds = 0
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 0, Press = 0
next a
jump SELF top
variable a loop 2000
run 1000
Neighbor list info ...
  update every 1 steps, delay 10 steps, check yes
  max neighbors/atom: 2000, page size: 100000
  master list distance cutoff = 8
  ghost atom cutoff = 8
  binsize = 4, bins = 11 11 11
  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.370 | 5.370 | 5.370 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    1000   -13168.698   -13168.698   -13305.677   -13305.677    264.99516    264.99516    67485.005    67485.005   -775.22965   -775.22965 
    2000   -13179.503   -13179.503   -13311.777   -13311.777    255.89202    255.89202    67400.777    67400.777   -242.38543   -242.38543 
Loop time of 753.254 on 1 procs for 1000 steps with 4000 atoms

Performance: 0.115 ns/day, 209.237 hours/ns, 1.328 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    | 752.19     | 752.19     | 752.19     |   0.0 | 99.86
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.16325    | 0.16325    | 0.16325    |   0.0 |  0.02
Output  | 0.00021358 | 0.00021358 | 0.00021358 |   0.0 |  0.00
Modify  | 0.79622    | 0.79622    | 0.79622    |   0.0 |  0.11
Other   |            | 0.1068     |            |       |  0.01

Nlocal:        4000.00 ave        4000 max        4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        5897.00 ave        5897 max        5897 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:      538504.0 ave      538504 max      538504 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 538504
Ave neighs/atom = 134.62600
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 = 11 11 11
  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.371 | 5.371 | 5.371 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    2000   -13179.503   -13179.503   -13311.777   -13311.777    255.89202    255.89202    67400.777    67400.777   -242.38543   -242.38543 
    3000   -13175.784   -13175.784    -13304.06    -13304.06    248.15923    248.15923      67373.4      67373.4    429.70871    429.70871 
Loop time of 705.236 on 1 procs for 1000 steps with 4000 atoms

Performance: 0.123 ns/day, 195.899 hours/ns, 1.418 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    | 704.24     | 704.24     | 704.24     |   0.0 | 99.86
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.16001    | 0.16001    | 0.16001    |   0.0 |  0.02
Output  | 0.00021443 | 0.00021443 | 0.00021443 |   0.0 |  0.00
Modify  | 0.73145    | 0.73145    | 0.73145    |   0.0 |  0.10
Other   |            | 0.1022     |            |       |  0.01

Nlocal:        4000.00 ave        4000 max        4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        5960.00 ave        5960 max        5960 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:      538526.0 ave      538526 max      538526 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 538526
Ave neighs/atom = 134.63150
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 = 11 11 11
  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.370 | 5.370 | 5.370 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    3000   -13175.784   -13175.784    -13304.06    -13304.06    248.15923    248.15923      67373.4      67373.4    429.70871    429.70871 
    4000     -13174.8     -13174.8   -13305.887   -13305.887    253.59679    253.59679     67427.33     67427.33   -177.53597   -177.53597 
Loop time of 686.089 on 1 procs for 1000 steps with 4000 atoms

Performance: 0.126 ns/day, 190.580 hours/ns, 1.458 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    | 685.14     | 685.14     | 685.14     |   0.0 | 99.86
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.1533     | 0.1533     | 0.1533     |   0.0 |  0.02
Output  | 0.00025714 | 0.00025714 | 0.00025714 |   0.0 |  0.00
Modify  | 0.69517    | 0.69517    | 0.69517    |   0.0 |  0.10
Other   |            | 0.09873    |            |       |  0.01

Nlocal:        4000.00 ave        4000 max        4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        5928.00 ave        5928 max        5928 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:      538908.0 ave      538908 max      538908 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 538908
Ave neighs/atom = 134.72700
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 = 11 11 11
  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.371 | 5.371 | 5.371 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    4000     -13174.8     -13174.8   -13305.887   -13305.887    253.59679    253.59679     67427.33     67427.33   -177.53597   -177.53597 
    5000   -13178.126   -13178.126   -13309.835   -13309.835    254.79939    254.79939    67368.247    67368.247    189.87564    189.87564 
Loop time of 650.813 on 1 procs for 1000 steps with 4000 atoms

Performance: 0.133 ns/day, 180.781 hours/ns, 1.537 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    | 649.92     | 649.92     | 649.92     |   0.0 | 99.86
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.14839    | 0.14839    | 0.14839    |   0.0 |  0.02
Output  | 0.00018133 | 0.00018133 | 0.00018133 |   0.0 |  0.00
Modify  | 0.64578    | 0.64578    | 0.64578    |   0.0 |  0.10
Other   |            | 0.09491    |            |       |  0.01

Nlocal:        4000.00 ave        4000 max        4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        5968.00 ave        5968 max        5968 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:      538632.0 ave      538632 max      538632 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 538632
Ave neighs/atom = 134.65800
Neighbor list builds = 0
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 253.29487616908, Press = -12.9847315368401
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 = 11 11 11
  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.371 | 5.371 | 5.371 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    5000   -13178.126   -13178.126   -13309.835   -13309.835    254.79939    254.79939    67368.247    67368.247    189.87564    189.87564 
    6000   -13174.304   -13174.304   -13307.289   -13307.289    257.26717    257.26717    67426.416    67426.416    -289.2266    -289.2266 
Loop time of 681.277 on 1 procs for 1000 steps with 4000 atoms

Performance: 0.127 ns/day, 189.244 hours/ns, 1.468 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    | 680.32     | 680.32     | 680.32     |   0.0 | 99.86
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.15262    | 0.15262    | 0.15262    |   0.0 |  0.02
Output  | 0.00017984 | 0.00017984 | 0.00017984 |   0.0 |  0.00
Modify  | 0.70933    | 0.70933    | 0.70933    |   0.0 |  0.10
Other   |            | 0.09897    |            |       |  0.01

Nlocal:        4000.00 ave        4000 max        4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        5960.00 ave        5960 max        5960 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:      538626.0 ave      538626 max      538626 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 538626
Ave neighs/atom = 134.65650
Neighbor list builds = 0
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 252.937270996606, Press = 7.47715836279101
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 = 11 11 11
  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.370 | 5.370 | 5.370 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    6000   -13174.304   -13174.304   -13307.289   -13307.289    257.26717    257.26717    67426.416    67426.416    -289.2266    -289.2266 
    7000   -13178.589   -13178.589   -13307.183   -13307.183    248.77523    248.77523    67422.028    67422.028   -276.92105   -276.92105 
Loop time of 648.401 on 1 procs for 1000 steps with 4000 atoms

Performance: 0.133 ns/day, 180.111 hours/ns, 1.542 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    | 647.49     | 647.49     | 647.49     |   0.0 | 99.86
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.14727    | 0.14727    | 0.14727    |   0.0 |  0.02
Output  | 0.00017779 | 0.00017779 | 0.00017779 |   0.0 |  0.00
Modify  | 0.66635    | 0.66635    | 0.66635    |   0.0 |  0.10
Other   |            | 0.09413    |            |       |  0.01

Nlocal:        4000.00 ave        4000 max        4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        5935.00 ave        5935 max        5935 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:      538606.0 ave      538606 max      538606 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 538606
Ave neighs/atom = 134.65150
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.088397440851, Press = -12.7749218396271
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 = 11 11 11
  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.371 | 5.371 | 5.371 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    7000   -13178.589   -13178.589   -13307.183   -13307.183    248.77523    248.77523    67422.028    67422.028   -276.92105   -276.92105 
    8000   -13173.598   -13173.598   -13304.332   -13304.332    252.91396    252.91396    67300.205    67300.205    1415.0297    1415.0297 
Loop time of 650.194 on 1 procs for 1000 steps with 4000 atoms

Performance: 0.133 ns/day, 180.609 hours/ns, 1.538 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    | 649.28     | 649.28     | 649.28     |   0.0 | 99.86
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.1499     | 0.1499     | 0.1499     |   0.0 |  0.02
Output  | 0.0001788  | 0.0001788  | 0.0001788  |   0.0 |  0.00
Modify  | 0.6689     | 0.6689     | 0.6689     |   0.0 |  0.10
Other   |            | 0.0941     |            |       |  0.01

Nlocal:        4000.00 ave        4000 max        4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        5987.00 ave        5987 max        5987 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:      538582.0 ave      538582 max      538582 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 538582
Ave neighs/atom = 134.64550
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.186361265465, Press = 4.11090360091798
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 = 11 11 11
  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.371 | 5.371 | 5.371 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    8000   -13173.598   -13173.598   -13304.332   -13304.332    252.91396    252.91396    67300.205    67300.205    1415.0297    1415.0297 
    9000   -13177.273   -13177.273   -13307.737   -13307.737    252.39281    252.39281    67485.641    67485.641   -1046.2144   -1046.2144 
Loop time of 650.475 on 1 procs for 1000 steps with 4000 atoms

Performance: 0.133 ns/day, 180.688 hours/ns, 1.537 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    | 649.57     | 649.57     | 649.57     |   0.0 | 99.86
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.14903    | 0.14903    | 0.14903    |   0.0 |  0.02
Output  | 0.00018138 | 0.00018138 | 0.00018138 |   0.0 |  0.00
Modify  | 0.66573    | 0.66573    | 0.66573    |   0.0 |  0.10
Other   |            | 0.09407    |            |       |  0.01

Nlocal:        4000.00 ave        4000 max        4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        5954.00 ave        5954 max        5954 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:      538916.0 ave      538916 max      538916 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 538916
Ave neighs/atom = 134.72900
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.198988276243, Press = 2.66977610570949
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 = 11 11 11
  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.370 | 5.370 | 5.370 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    9000   -13177.273   -13177.273   -13307.737   -13307.737    252.39281    252.39281    67485.641    67485.641   -1046.2144   -1046.2144 
   10000   -13175.042   -13175.042    -13309.23    -13309.23    259.59539    259.59539    67403.278    67403.278   -134.60224   -134.60224 
Loop time of 649.941 on 1 procs for 1000 steps with 4000 atoms

Performance: 0.133 ns/day, 180.539 hours/ns, 1.539 timesteps/s
100.0% CPU use with 1 MPI tasks x 1 OpenMP threads

MPI task timing breakdown:
Section |  min time  |  avg time  |  max time  |%varavg| %total
---------------------------------------------------------------
Pair    | 649.03     | 649.03     | 649.03     |   0.0 | 99.86
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.14772    | 0.14772    | 0.14772    |   0.0 |  0.02
Output  | 0.00018067 | 0.00018067 | 0.00018067 |   0.0 |  0.00
Modify  | 0.66455    | 0.66455    | 0.66455    |   0.0 |  0.10
Other   |            | 0.09435    |            |       |  0.01

Nlocal:        4000.00 ave        4000 max        4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        5941.00 ave        5941 max        5941 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:      538534.0 ave      538534 max      538534 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 538534
Ave neighs/atom = 134.63350
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.560875826757, Press = -4.98538839784176
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 = 11 11 11
  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.371 | 5.371 | 5.371 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   10000   -13175.042   -13175.042    -13309.23    -13309.23    259.59539    259.59539    67403.278    67403.278   -134.60224   -134.60224 
   11000   -13175.697   -13175.697   -13308.366   -13308.366    256.65771    256.65771    67373.649    67373.649    287.02954    287.02954 
Loop time of 649.708 on 1 procs for 1000 steps with 4000 atoms

Performance: 0.133 ns/day, 180.474 hours/ns, 1.539 timesteps/s
100.0% CPU use with 1 MPI tasks x 1 OpenMP threads

MPI task timing breakdown:
Section |  min time  |  avg time  |  max time  |%varavg| %total
---------------------------------------------------------------
Pair    | 648.79     | 648.79     | 648.79     |   0.0 | 99.86
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.15092    | 0.15092    | 0.15092    |   0.0 |  0.02
Output  | 0.00018042 | 0.00018042 | 0.00018042 |   0.0 |  0.00
Modify  | 0.66818    | 0.66818    | 0.66818    |   0.0 |  0.10
Other   |            | 0.09548    |            |       |  0.01

Nlocal:        4000.00 ave        4000 max        4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        5967.00 ave        5967 max        5967 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:      538688.0 ave      538688 max      538688 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 538688
Ave neighs/atom = 134.67200
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.510795109725, Press = 0.906537415570443
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 = 11 11 11
  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.371 | 5.371 | 5.371 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   11000   -13175.697   -13175.697   -13308.366   -13308.366    256.65771    256.65771    67373.649    67373.649    287.02954    287.02954 
   12000   -13179.194   -13179.194   -13305.909   -13305.909     245.1378     245.1378    67409.625    67409.625   -103.52225   -103.52225 
Loop time of 741.7 on 1 procs for 1000 steps with 4000 atoms

Performance: 0.116 ns/day, 206.028 hours/ns, 1.348 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    | 740.62     | 740.62     | 740.62     |   0.0 | 99.85
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.16696    | 0.16696    | 0.16696    |   0.0 |  0.02
Output  | 0.00018341 | 0.00018341 | 0.00018341 |   0.0 |  0.00
Modify  | 0.80642    | 0.80642    | 0.80642    |   0.0 |  0.11
Other   |            | 0.1065     |            |       |  0.01

Nlocal:        4000.00 ave        4000 max        4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        5960.00 ave        5960 max        5960 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:      538680.0 ave      538680 max      538680 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 538680
Ave neighs/atom = 134.67000
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.50953287894, Press = 0.45237622538783
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 = 11 11 11
  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.371 | 5.371 | 5.371 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   12000   -13179.194   -13179.194   -13305.909   -13305.909     245.1378     245.1378    67409.625    67409.625   -103.52225   -103.52225 
   13000   -13175.011   -13175.011   -13308.373   -13308.373    257.99788    257.99788    67399.691    67399.691   -46.848653   -46.848653 
Loop time of 710.654 on 1 procs for 1000 steps with 4000 atoms

Performance: 0.122 ns/day, 197.404 hours/ns, 1.407 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    | 709.61     | 709.61     | 709.61     |   0.0 | 99.85
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.16375    | 0.16375    | 0.16375    |   0.0 |  0.02
Output  | 0.00024789 | 0.00024789 | 0.00024789 |   0.0 |  0.00
Modify  | 0.77275    | 0.77275    | 0.77275    |   0.0 |  0.11
Other   |            | 0.1048     |            |       |  0.01

Nlocal:        4000.00 ave        4000 max        4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        5989.00 ave        5989 max        5989 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:      538792.0 ave      538792 max      538792 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 538792
Ave neighs/atom = 134.69800
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.329310513354, Press = -1.58550514990831
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 = 11 11 11
  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.370 | 5.370 | 5.370 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   13000   -13175.011   -13175.011   -13308.373   -13308.373    257.99788    257.99788    67399.691    67399.691   -46.848653   -46.848653 
   14000   -13178.674   -13178.674   -13309.332   -13309.332    252.76561    252.76561    67375.566    67375.566     124.0065     124.0065 
Loop time of 730.798 on 1 procs for 1000 steps with 4000 atoms

Performance: 0.118 ns/day, 202.999 hours/ns, 1.368 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    | 729.75     | 729.75     | 729.75     |   0.0 | 99.86
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.16206    | 0.16206    | 0.16206    |   0.0 |  0.02
Output  | 0.00017702 | 0.00017702 | 0.00017702 |   0.0 |  0.00
Modify  | 0.77982    | 0.77982    | 0.77982    |   0.0 |  0.11
Other   |            | 0.105      |            |       |  0.01

Nlocal:        4000.00 ave        4000 max        4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        5939.00 ave        5939 max        5939 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:      538726.0 ave      538726 max      538726 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 538726
Ave neighs/atom = 134.68150
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.157047172824, Press = 0.619283128501492
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 = 11 11 11
  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.371 | 5.371 | 5.371 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   14000   -13178.674   -13178.674   -13309.332   -13309.332    252.76561    252.76561    67375.566    67375.566     124.0065     124.0065 
   15000   -13176.391   -13176.391   -13306.028   -13306.028    250.79245    250.79245    67419.292    67419.292   -154.29939   -154.29939 
Loop time of 749.277 on 1 procs for 1000 steps with 4000 atoms

Performance: 0.115 ns/day, 208.132 hours/ns, 1.335 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    | 748.19     | 748.19     | 748.19     |   0.0 | 99.85
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.16735    | 0.16735    | 0.16735    |   0.0 |  0.02
Output  | 0.00017895 | 0.00017895 | 0.00017895 |   0.0 |  0.00
Modify  | 0.81499    | 0.81499    | 0.81499    |   0.0 |  0.11
Other   |            | 0.1082     |            |       |  0.01

Nlocal:        4000.00 ave        4000 max        4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        5957.00 ave        5957 max        5957 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:      538608.0 ave      538608 max      538608 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 538608
Ave neighs/atom = 134.65200
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.135406080572, Press = -1.06852926512365
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 = 11 11 11
  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.371 | 5.371 | 5.371 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   15000   -13176.391   -13176.391   -13306.028   -13306.028    250.79245    250.79245    67419.292    67419.292   -154.29939   -154.29939 
   16000   -13178.624   -13178.624   -13308.601   -13308.601    251.44974    251.44974     67379.15     67379.15    192.31406    192.31406 
Loop time of 757.343 on 1 procs for 1000 steps with 4000 atoms

Performance: 0.114 ns/day, 210.373 hours/ns, 1.320 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    | 756.23     | 756.23     | 756.23     |   0.0 | 99.85
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.16918    | 0.16918    | 0.16918    |   0.0 |  0.02
Output  | 0.00018196 | 0.00018196 | 0.00018196 |   0.0 |  0.00
Modify  | 0.834      | 0.834      | 0.834      |   0.0 |  0.11
Other   |            | 0.1083     |            |       |  0.01

Nlocal:        4000.00 ave        4000 max        4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        5985.00 ave        5985 max        5985 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:      538688.0 ave      538688 max      538688 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 538688
Ave neighs/atom = 134.67200
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.166608133292, Press = -0.916188296221891
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 = 11 11 11
  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.371 | 5.371 | 5.371 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   16000   -13178.624   -13178.624   -13308.601   -13308.601    251.44974    251.44974     67379.15     67379.15    192.31406    192.31406 
   17000   -13172.251   -13172.251   -13305.635   -13305.635    258.04123    258.04123    67416.775    67416.775   -24.460322   -24.460322 
Loop time of 708.401 on 1 procs for 1000 steps with 4000 atoms

Performance: 0.122 ns/day, 196.778 hours/ns, 1.412 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    | 707.39     | 707.39     | 707.39     |   0.0 | 99.86
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.1592     | 0.1592     | 0.1592     |   0.0 |  0.02
Output  | 0.00017747 | 0.00017747 | 0.00017747 |   0.0 |  0.00
Modify  | 0.75327    | 0.75327    | 0.75327    |   0.0 |  0.11
Other   |            | 0.1021     |            |       |  0.01

Nlocal:        4000.00 ave        4000 max        4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        5959.00 ave        5959 max        5959 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:      538618.0 ave      538618 max      538618 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 538618
Ave neighs/atom = 134.65450
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
67401.0751348438

print "LAMMPS calculation completed"
LAMMPS calculation completed
quit 0