# 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.9173715263605122*${_u_distance}
variable latticeconst_converted equal 3.9173715263605122*1
lattice       fcc  ${latticeconst_converted}
lattice       fcc  3.91737152636051
Lattice spacing in x,y,z = 3.9173715 3.9173715 3.9173715
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 (39.173715 39.173715 39.173715)
  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 (39.173715 39.173715 39.173715)
  create_atoms CPU = 0.005 seconds

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

kim_interactions Pt

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

kim interactions Pt
#=== BEGIN kim interactions ==================================
pair_style kim MEAM_LAMMPS_KimSeolJi_2017_PtAl__MO_793141037706_001
pair_coeff * * Pt
#=== END kim interactions ====================================


mass          1 ${mass_converted}
mass          1 195.078

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

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

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

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

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

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

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

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

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

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

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

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

Your simulation uses code contributions which should be cited:

- OpenKIM: https://doi.org/10.1007/s11837-011-0102-6

@Article{tadmor:elliott:2011,
 author = {E. B. Tadmor and R. S. Elliott and J. P. Sethna and R. E. Miller and C. A. Becker},
 title = {The potential of atomistic simulations and the {K}nowledgebase of {I}nteratomic {M}odels},
 journal = {{JOM}},
 year =    2011,
 volume =  63,
 number =  17,
 pages =   {17},
 doi =     {10.1007/s11837-011-0102-6}
}

- OpenKIM potential: https://openkim.org/cite/MO_793141037706_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   -22938.806   -22938.806       -23080       -23080       273.15       273.15    60115.199    60115.199    2508.7109    2508.7109 
    1000   -22789.496   -22789.496   -22940.184   -22940.184    291.51711    291.51711    60599.473    60599.473    489.84391    489.84391 
Loop time of 417.923 on 1 procs for 1000 steps with 4000 atoms

Performance: 0.207 ns/day, 116.090 hours/ns, 2.393 timesteps/s
100.0% CPU use with 1 MPI tasks x 1 OpenMP threads

MPI task timing breakdown:
Section |  min time  |  avg time  |  max time  |%varavg| %total
---------------------------------------------------------------
Pair    | 416.83     | 416.83     | 416.83     |   0.0 | 99.74
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.21545    | 0.21545    | 0.21545    |   0.0 |  0.05
Output  | 0.0002588  | 0.0002588  | 0.0002588  |   0.0 |  0.00
Modify  | 0.76974    | 0.76974    | 0.76974    |   0.0 |  0.18
Other   |            | 0.1065     |            |       |  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   -22789.496   -22789.496   -22940.184   -22940.184    291.51711    291.51711    60599.473    60599.473    489.84391    489.84391 
    2000   -22797.661   -22797.661   -22943.729   -22943.729    282.57906    282.57906    60617.878    60617.878   -1031.1647   -1031.1647 
Loop time of 418.161 on 1 procs for 1000 steps with 4000 atoms

Performance: 0.207 ns/day, 116.156 hours/ns, 2.391 timesteps/s
99.9% CPU use with 1 MPI tasks x 1 OpenMP threads

MPI task timing breakdown:
Section |  min time  |  avg time  |  max time  |%varavg| %total
---------------------------------------------------------------
Pair    | 417.1      | 417.1      | 417.1      |   0.0 | 99.75
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.20955    | 0.20955    | 0.20955    |   0.0 |  0.05
Output  | 0.00031466 | 0.00031466 | 0.00031466 |   0.0 |  0.00
Modify  | 0.74807    | 0.74807    | 0.74807    |   0.0 |  0.18
Other   |            | 0.1072     |            |       |  0.03

Nlocal:        4000.00 ave        4000 max        4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        8179.00 ave        8179 max        8179 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:      559068.0 ave      559068 max      559068 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 559068
Ave neighs/atom = 139.76700
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   -22797.661   -22797.661   -22943.729   -22943.729    282.57906    282.57906    60617.878    60617.878   -1031.1647   -1031.1647 
    3000   -22798.999   -22798.999   -22943.632   -22943.632    279.80238    279.80238    60631.018    60631.018   -1688.4924   -1688.4924 
Loop time of 434.257 on 1 procs for 1000 steps with 4000 atoms

Performance: 0.199 ns/day, 120.627 hours/ns, 2.303 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    | 433.15     | 433.15     | 433.15     |   0.0 | 99.74
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.21663    | 0.21663    | 0.21663    |   0.0 |  0.05
Output  | 0.00032232 | 0.00032232 | 0.00032232 |   0.0 |  0.00
Modify  | 0.78304    | 0.78304    | 0.78304    |   0.0 |  0.18
Other   |            | 0.1089     |            |       |  0.03

Nlocal:        4000.00 ave        4000 max        4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        8184.00 ave        8184 max        8184 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:      559094.0 ave      559094 max      559094 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 559094
Ave neighs/atom = 139.77350
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   -22798.999   -22798.999   -22943.632   -22943.632    279.80238    279.80238    60631.018    60631.018   -1688.4924   -1688.4924 
    4000   -22794.726   -22794.726   -22942.034   -22942.034    284.97901    284.97901    60545.098    60545.098    2878.1302    2878.1302 
Loop time of 436.522 on 1 procs for 1000 steps with 4000 atoms

Performance: 0.198 ns/day, 121.256 hours/ns, 2.291 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    | 435.42     | 435.42     | 435.42     |   0.0 | 99.75
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.21635    | 0.21635    | 0.21635    |   0.0 |  0.05
Output  | 0.00021824 | 0.00021824 | 0.00021824 |   0.0 |  0.00
Modify  | 0.77914    | 0.77914    | 0.77914    |   0.0 |  0.18
Other   |            | 0.1093     |            |       |  0.03

Nlocal:        4000.00 ave        4000 max        4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        8171.00 ave        8171 max        8171 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:      558994.0 ave      558994 max      558994 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 558994
Ave neighs/atom = 139.74850
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   -22794.726   -22794.726   -22942.034   -22942.034    284.97901    284.97901    60545.098    60545.098    2878.1302    2878.1302 
    5000   -22799.179   -22799.179   -22938.983   -22938.983      270.461      270.461    60605.748    60605.748    69.248754    69.248754 
Loop time of 406.051 on 1 procs for 1000 steps with 4000 atoms

Performance: 0.213 ns/day, 112.792 hours/ns, 2.463 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    | 405.03     | 405.03     | 405.03     |   0.0 | 99.75
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.20562    | 0.20562    | 0.20562    |   0.0 |  0.05
Output  | 0.00017871 | 0.00017871 | 0.00017871 |   0.0 |  0.00
Modify  | 0.71091    | 0.71091    | 0.71091    |   0.0 |  0.18
Other   |            | 0.102      |            |       |  0.03

Nlocal:        4000.00 ave        4000 max        4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        8183.00 ave        8183 max        8183 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:      559202.0 ave      559202 max      559202 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 559202
Ave neighs/atom = 139.80050
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 = 275.94552130403, Press = 131.896667209406
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   -22799.179   -22799.179   -22938.983   -22938.983      270.461      270.461    60605.748    60605.748    69.248754    69.248754 
    6000   -22798.385   -22798.385   -22940.749   -22940.749    275.41298    275.41298    60617.017    60617.017   -588.31345   -588.31345 
Loop time of 416.505 on 1 procs for 1000 steps with 4000 atoms

Performance: 0.207 ns/day, 115.696 hours/ns, 2.401 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    | 415.42     | 415.42     | 415.42     |   0.0 | 99.74
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.20803    | 0.20803    | 0.20803    |   0.0 |  0.05
Output  | 0.00017829 | 0.00017829 | 0.00017829 |   0.0 |  0.00
Modify  | 0.76501    | 0.76501    | 0.76501    |   0.0 |  0.18
Other   |            | 0.1066     |            |       |  0.03

Nlocal:        4000.00 ave        4000 max        4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        8163.00 ave        8163 max        8163 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:      558924.0 ave      558924 max      558924 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 558924
Ave neighs/atom = 139.73100
Neighbor list builds = 0
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 273.532538134711, Press = 11.9033586715957
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   -22798.385   -22798.385   -22940.749   -22940.749    275.41298    275.41298    60617.017    60617.017   -588.31345   -588.31345 
    7000     -22796.9     -22796.9   -22936.969   -22936.969    270.97172    270.97172    60582.704    60582.704    1460.3447    1460.3447 
Loop time of 425.028 on 1 procs for 1000 steps with 4000 atoms

Performance: 0.203 ns/day, 118.063 hours/ns, 2.353 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    | 423.91     | 423.91     | 423.91     |   0.0 | 99.74
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.21473    | 0.21473    | 0.21473    |   0.0 |  0.05
Output  | 0.00027975 | 0.00027975 | 0.00027975 |   0.0 |  0.00
Modify  | 0.79135    | 0.79135    | 0.79135    |   0.0 |  0.19
Other   |            | 0.1079     |            |       |  0.03

Nlocal:        4000.00 ave        4000 max        4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        8166.00 ave        8166 max        8166 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:      559048.0 ave      559048 max      559048 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 559048
Ave neighs/atom = 139.76200
Neighbor list builds = 0
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 273.027415005748, Press = 16.8218088061216
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     -22796.9     -22796.9   -22936.969   -22936.969    270.97172    270.97172    60582.704    60582.704    1460.3447    1460.3447 
    8000   -22797.826   -22797.826   -22937.808   -22937.808    270.80583    270.80583    60621.234    60621.234   -312.39972   -312.39972 
Loop time of 406.489 on 1 procs for 1000 steps with 4000 atoms

Performance: 0.213 ns/day, 112.914 hours/ns, 2.460 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    | 405.43     | 405.43     | 405.43     |   0.0 | 99.74
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.20807    | 0.20807    | 0.20807    |   0.0 |  0.05
Output  | 0.00023012 | 0.00023012 | 0.00023012 |   0.0 |  0.00
Modify  | 0.74233    | 0.74233    | 0.74233    |   0.0 |  0.18
Other   |            | 0.1037     |            |       |  0.03

Nlocal:        4000.00 ave        4000 max        4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        8169.00 ave        8169 max        8169 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:      558964.0 ave      558964 max      558964 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 558964
Ave neighs/atom = 139.74100
Neighbor list builds = 0
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 273.35230291343, Press = -0.0165300350611712
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   -22797.826   -22797.826   -22937.808   -22937.808    270.80583    270.80583    60621.234    60621.234   -312.39972   -312.39972 
    9000   -22801.061   -22801.061   -22941.336   -22941.336    271.37238    271.37238    60609.338    60609.338    -316.3414    -316.3414 
Loop time of 446.821 on 1 procs for 1000 steps with 4000 atoms

Performance: 0.193 ns/day, 124.117 hours/ns, 2.238 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    | 445.67     | 445.67     | 445.67     |   0.0 | 99.74
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.21917    | 0.21917    | 0.21917    |   0.0 |  0.05
Output  | 0.00027846 | 0.00027846 | 0.00027846 |   0.0 |  0.00
Modify  | 0.82495    | 0.82495    | 0.82495    |   0.0 |  0.18
Other   |            | 0.1095     |            |       |  0.02

Nlocal:        4000.00 ave        4000 max        4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        8178.00 ave        8178 max        8178 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:      559004.0 ave      559004 max      559004 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 559004
Ave neighs/atom = 139.75100
Neighbor list builds = 0
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 272.98959018625, Press = 13.8872429556421
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   -22801.061   -22801.061   -22941.336   -22941.336    271.37238    271.37238    60609.338    60609.338    -316.3414    -316.3414 
   10000   -22795.712   -22795.712   -22935.439   -22935.439    270.30977    270.30977    60664.137    60664.137    -2085.703    -2085.703 
Loop time of 437.498 on 1 procs for 1000 steps with 4000 atoms

Performance: 0.197 ns/day, 121.527 hours/ns, 2.286 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    | 436.35     | 436.35     | 436.35     |   0.0 | 99.74
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.21966    | 0.21966    | 0.21966    |   0.0 |  0.05
Output  | 0.00017937 | 0.00017937 | 0.00017937 |   0.0 |  0.00
Modify  | 0.81468    | 0.81468    | 0.81468    |   0.0 |  0.19
Other   |            | 0.1093     |            |       |  0.02

Nlocal:        4000.00 ave        4000 max        4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        8184.00 ave        8184 max        8184 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:      559044.0 ave      559044 max      559044 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 559044
Ave neighs/atom = 139.76100
Neighbor list builds = 0
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 272.998511031945, Press = 2.08846264700363
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   -22795.712   -22795.712   -22935.439   -22935.439    270.30977    270.30977    60664.137    60664.137    -2085.703    -2085.703 
   11000   -22804.195   -22804.195   -22943.267   -22943.267     269.0434     269.0434    60605.469    60605.469   -380.72293   -380.72293 
Loop time of 442.793 on 1 procs for 1000 steps with 4000 atoms

Performance: 0.195 ns/day, 122.998 hours/ns, 2.258 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    | 441.64     | 441.64     | 441.64     |   0.0 | 99.74
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.2197     | 0.2197     | 0.2197     |   0.0 |  0.05
Output  | 0.0002588  | 0.0002588  | 0.0002588  |   0.0 |  0.00
Modify  | 0.82365    | 0.82365    | 0.82365    |   0.0 |  0.19
Other   |            | 0.1096     |            |       |  0.02

Nlocal:        4000.00 ave        4000 max        4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        8165.00 ave        8165 max        8165 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:      558912.0 ave      558912 max      558912 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 558912
Ave neighs/atom = 139.72800
Neighbor list builds = 0
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 272.802919438393, Press = -1.61744260840005
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   -22804.195   -22804.195   -22943.267   -22943.267     269.0434     269.0434    60605.469    60605.469   -380.72293   -380.72293 
   12000   -22799.912   -22799.912   -22938.616   -22938.616    268.33234    268.33234    60580.534    60580.534    1426.9363    1426.9363 
Loop time of 437.714 on 1 procs for 1000 steps with 4000 atoms

Performance: 0.197 ns/day, 121.587 hours/ns, 2.285 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    | 436.56     | 436.56     | 436.56     |   0.0 | 99.74
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.21903    | 0.21903    | 0.21903    |   0.0 |  0.05
Output  | 0.00023966 | 0.00023966 | 0.00023966 |   0.0 |  0.00
Modify  | 0.82144    | 0.82144    | 0.82144    |   0.0 |  0.19
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:        8172.00 ave        8172 max        8172 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:      558998.0 ave      558998 max      558998 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 558998
Ave neighs/atom = 139.74950
Neighbor list builds = 0
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 272.812899944052, Press = -0.386427048349793
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   -22799.912   -22799.912   -22938.616   -22938.616    268.33234    268.33234    60580.534    60580.534    1426.9363    1426.9363 
   13000   -22799.088   -22799.088   -22941.514   -22941.514    275.53394    275.53394     60590.72     60590.72    567.64017    567.64017 
Loop time of 415.187 on 1 procs for 1000 steps with 4000 atoms

Performance: 0.208 ns/day, 115.330 hours/ns, 2.409 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    | 414.12     | 414.12     | 414.12     |   0.0 | 99.74
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.20841    | 0.20841    | 0.20841    |   0.0 |  0.05
Output  | 0.00022258 | 0.00022258 | 0.00022258 |   0.0 |  0.00
Modify  | 0.75016    | 0.75016    | 0.75016    |   0.0 |  0.18
Other   |            | 0.1041     |            |       |  0.03

Nlocal:        4000.00 ave        4000 max        4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        8177.00 ave        8177 max        8177 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:      559030.0 ave      559030 max      559030 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 559030
Ave neighs/atom = 139.75750
Neighbor list builds = 0
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 272.862498086691, Press = -0.266659205252523
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   -22799.088   -22799.088   -22941.514   -22941.514    275.53394    275.53394     60590.72     60590.72    567.64017    567.64017 
   14000   -22794.902   -22794.902   -22939.351   -22939.351    279.44728    279.44728     60567.58     60567.58    2050.5218    2050.5218 
Loop time of 380.149 on 1 procs for 1000 steps with 4000 atoms

Performance: 0.227 ns/day, 105.597 hours/ns, 2.631 timesteps/s
100.1% CPU use with 1 MPI tasks x 1 OpenMP threads

MPI task timing breakdown:
Section |  min time  |  avg time  |  max time  |%varavg| %total
---------------------------------------------------------------
Pair    | 379.19     | 379.19     | 379.19     |   0.0 | 99.75
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.19696    | 0.19696    | 0.19696    |   0.0 |  0.05
Output  | 0.00017887 | 0.00017887 | 0.00017887 |   0.0 |  0.00
Modify  | 0.66932    | 0.66932    | 0.66932    |   0.0 |  0.18
Other   |            | 0.09631    |            |       |  0.03

Nlocal:        4000.00 ave        4000 max        4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        8167.00 ave        8167 max        8167 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:      559012.0 ave      559012 max      559012 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 559012
Ave neighs/atom = 139.75300
Neighbor list builds = 0
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 273.109028405836, Press = 2.64151392606304
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   -22794.902   -22794.902   -22939.351   -22939.351    279.44728    279.44728     60567.58     60567.58    2050.5218    2050.5218 
   15000   -22798.592   -22798.592   -22939.599   -22939.599    272.78683    272.78683     60612.56     60612.56   -155.61351   -155.61351 
Loop time of 379.344 on 1 procs for 1000 steps with 4000 atoms

Performance: 0.228 ns/day, 105.373 hours/ns, 2.636 timesteps/s
100.1% CPU use with 1 MPI tasks x 1 OpenMP threads

MPI task timing breakdown:
Section |  min time  |  avg time  |  max time  |%varavg| %total
---------------------------------------------------------------
Pair    | 378.38     | 378.38     | 378.38     |   0.0 | 99.75
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.19786    | 0.19786    | 0.19786    |   0.0 |  0.05
Output  | 0.00018261 | 0.00018261 | 0.00018261 |   0.0 |  0.00
Modify  | 0.66923    | 0.66923    | 0.66923    |   0.0 |  0.18
Other   |            | 0.09639    |            |       |  0.03

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:      559158.0 ave      559158 max      559158 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 559158
Ave neighs/atom = 139.78950
Neighbor list builds = 0
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 273.08606044451, Press = 4.19630467559243
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   -22798.592   -22798.592   -22939.599   -22939.599    272.78683    272.78683     60612.56     60612.56   -155.61351   -155.61351 
   16000   -22798.422   -22798.422   -22939.047   -22939.047    272.04826    272.04826     60632.07     60632.07   -1114.7537   -1114.7537 
Loop time of 378.621 on 1 procs for 1000 steps with 4000 atoms

Performance: 0.228 ns/day, 105.172 hours/ns, 2.641 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    | 377.66     | 377.66     | 377.66     |   0.0 | 99.75
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.19852    | 0.19852    | 0.19852    |   0.0 |  0.05
Output  | 0.00022455 | 0.00022455 | 0.00022455 |   0.0 |  0.00
Modify  | 0.66775    | 0.66775    | 0.66775    |   0.0 |  0.18
Other   |            | 0.09589    |            |       |  0.03

Nlocal:        4000.00 ave        4000 max        4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        8168.00 ave        8168 max        8168 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:      559016.0 ave      559016 max      559016 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 559016
Ave neighs/atom = 139.75400
Neighbor list builds = 0
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 273.058571631352, Press = 1.4783357940203
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   -22798.422   -22798.422   -22939.047   -22939.047    272.04826    272.04826     60632.07     60632.07   -1114.7537   -1114.7537 
   17000   -22795.591   -22795.591   -22939.715   -22939.715    278.81722    278.81722    60656.376    60656.376   -2271.1768   -2271.1768 
Loop time of 379.78 on 1 procs for 1000 steps with 4000 atoms

Performance: 0.228 ns/day, 105.494 hours/ns, 2.633 timesteps/s
100.1% CPU use with 1 MPI tasks x 1 OpenMP threads

MPI task timing breakdown:
Section |  min time  |  avg time  |  max time  |%varavg| %total
---------------------------------------------------------------
Pair    | 378.82     | 378.82     | 378.82     |   0.0 | 99.75
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.19727    | 0.19727    | 0.19727    |   0.0 |  0.05
Output  | 0.00018402 | 0.00018402 | 0.00018402 |   0.0 |  0.00
Modify  | 0.66906    | 0.66906    | 0.66906    |   0.0 |  0.18
Other   |            | 0.09717    |            |       |  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:      558976.0 ave      558976 max      558976 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 558976
Ave neighs/atom = 139.74400
Neighbor list builds = 0
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 273.042816660225, Press = -0.452864199954119
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   -22795.591   -22795.591   -22939.715   -22939.715    278.81722    278.81722    60656.376    60656.376   -2271.1768   -2271.1768 
   18000   -22800.898   -22800.898   -22939.275   -22939.275    267.69923    267.69923    60589.245    60589.245    827.18764    827.18764 
Loop time of 379.557 on 1 procs for 1000 steps with 4000 atoms

Performance: 0.228 ns/day, 105.432 hours/ns, 2.635 timesteps/s
100.1% CPU use with 1 MPI tasks x 1 OpenMP threads

MPI task timing breakdown:
Section |  min time  |  avg time  |  max time  |%varavg| %total
---------------------------------------------------------------
Pair    | 378.6      | 378.6      | 378.6      |   0.0 | 99.75
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.19658    | 0.19658    | 0.19658    |   0.0 |  0.05
Output  | 0.00017963 | 0.00017963 | 0.00017963 |   0.0 |  0.00
Modify  | 0.66769    | 0.66769    | 0.66769    |   0.0 |  0.18
Other   |            | 0.09663    |            |       |  0.03

Nlocal:        4000.00 ave        4000 max        4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        8180.00 ave        8180 max        8180 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:      558928.0 ave      558928 max      558928 min
Histogram: 1 0 0 0 0 0 0 0 0 0

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

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

print "LAMMPS calculation completed"
LAMMPS calculation completed
quit 0