# 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 2.8809780031442638*${_u_distance}
variable latticeconst_converted equal 2.8809780031442638*1
lattice       bcc  ${latticeconst_converted}
lattice       bcc  2.88097800314426
Lattice spacing in x,y,z = 2.8809780 2.8809780 2.8809780
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 (28.809780 28.809780 28.809780)
  1 by 1 by 1 MPI processor grid
create_atoms  1 box
Created 2000 atoms
  create_atoms CPU = 0.000 seconds

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

kim_interactions Cr
#=== BEGIN kim_interactions ==================================
pair_style kim MEAM_LAMMPS_WuLeeSu_2017_NiCrFe__MO_912636107108_000
pair_coeff * * Cr
#=== END kim_interactions ====================================


mass          1 ${mass_converted}
mass          1 51.9961

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

print "Initial system volume: ${V0_metal} Angstroms^3"
Initial system volume: 23912.2161128496 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 293.15*${_u_temperature}
variable temp_converted equal 293.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 293.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 293.15 ${temp_converted} ${Tdamp_converted}               iso ${press_converted} ${press_converted} ${Pdamp_converted}
fix           ensemble all npt temp 293.15 293.15 ${Tdamp_converted}               iso ${press_converted} ${press_converted} ${Pdamp_converted}
fix           ensemble all npt temp 293.15 293.15 0.1               iso ${press_converted} ${press_converted} ${Pdamp_converted}
fix           ensemble all npt temp 293.15 293.15 0.1               iso 0 ${press_converted} ${Pdamp_converted}
fix           ensemble all npt temp 293.15 293.15 0.1               iso 0 0 ${Pdamp_converted}
fix           ensemble all npt temp 293.15 293.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  "293.15 - 0.2"
variable      T_up  equal  "293.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
Neighbor list info ...
  update every 1 steps, delay 10 steps, check yes
  max neighbors/atom: 2000, page size: 100000
  master list distance cutoff = 6
  ghost atom cutoff = 6
  binsize = 3, bins = 10 10 10
  1 neighbor lists, perpetual/occasional/extra = 1 0 0
  (1) pair kim, perpetual
      attributes: full, newton off, cut 6
      pair build: full/bin/atomonly
      stencil: full/bin/3d
      bin: standard
Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
       0   -8124.2527   -8124.2527        -8200        -8200       293.15       293.15    23912.216    23912.216     3383.522     3383.522 
    1000    -8045.708    -8045.708   -8121.8384   -8121.8384    294.63219    294.63219    24200.826    24200.826   -2946.9142   -2946.9142 
Loop time of 28.3758 on 1 procs for 1000 steps with 2000 atoms

Performance: 3.045 ns/day, 7.882 hours/ns, 35.241 timesteps/s
74.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    | 28.234     | 28.234     | 28.234     |   0.0 | 99.50
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.026335   | 0.026335   | 0.026335   |   0.0 |  0.09
Output  | 4.21e-05   | 4.21e-05   | 4.21e-05   |   0.0 |  0.00
Modify  | 0.10471    | 0.10471    | 0.10471    |   0.0 |  0.37
Other   |            | 0.01099    |            |       |  0.04

Nlocal:        2000.00 ave        2000 max        2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        4119.00 ave        4119 max        4119 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:      128000.0 ave      128000 max      128000 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 128000
Ave neighs/atom = 64.000000
Neighbor list builds = 0
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 0, Press = 0
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    1000    -8045.708    -8045.708   -8121.8384   -8121.8384    294.63219    294.63219    24200.826    24200.826   -2946.9142   -2946.9142 
    2000   -8044.8375   -8044.8375   -8124.4146   -8124.4146    307.97166    307.97166    24142.102    24142.102    1344.2956    1344.2956 
Loop time of 31.8949 on 1 procs for 1000 steps with 2000 atoms

Performance: 2.709 ns/day, 8.860 hours/ns, 31.353 timesteps/s
69.4% CPU use with 1 MPI tasks x 1 OpenMP threads

MPI task timing breakdown:
Section |  min time  |  avg time  |  max time  |%varavg| %total
---------------------------------------------------------------
Pair    | 31.716     | 31.716     | 31.716     |   0.0 | 99.44
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.046526   | 0.046526   | 0.046526   |   0.0 |  0.15
Output  | 3.69e-05   | 3.69e-05   | 3.69e-05   |   0.0 |  0.00
Modify  | 0.12081    | 0.12081    | 0.12081    |   0.0 |  0.38
Other   |            | 0.01116    |            |       |  0.03

Nlocal:        2000.00 ave        2000 max        2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        4119.00 ave        4119 max        4119 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:      127948.0 ave      127948 max      127948 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 127948
Ave neighs/atom = 63.974000
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
Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    2000   -8044.8375   -8044.8375   -8124.4146   -8124.4146    307.97166    307.97166    24142.102    24142.102    1344.2956    1344.2956 
    3000   -8047.7457   -8047.7457   -8124.0785   -8124.0785    295.41594    295.41594    24178.527    24178.527   -1845.3074   -1845.3074 
Loop time of 29.1782 on 1 procs for 1000 steps with 2000 atoms

Performance: 2.961 ns/day, 8.105 hours/ns, 34.272 timesteps/s
75.6% CPU use with 1 MPI tasks x 1 OpenMP threads

MPI task timing breakdown:
Section |  min time  |  avg time  |  max time  |%varavg| %total
---------------------------------------------------------------
Pair    | 29.016     | 29.016     | 29.016     |   0.0 | 99.44
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.026338   | 0.026338   | 0.026338   |   0.0 |  0.09
Output  | 4.39e-05   | 4.39e-05   | 4.39e-05   |   0.0 |  0.00
Modify  | 0.12484    | 0.12484    | 0.12484    |   0.0 |  0.43
Other   |            | 0.011      |            |       |  0.04

Nlocal:        2000.00 ave        2000 max        2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        4119.00 ave        4119 max        4119 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:      127986.0 ave      127986 max      127986 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 127986
Ave neighs/atom = 63.993000
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
Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    3000   -8047.7457   -8047.7457   -8124.0785   -8124.0785    295.41594    295.41594    24178.527    24178.527   -1845.3074   -1845.3074 
    4000    -8043.845    -8043.845   -8124.4005   -8124.4005    311.75789    311.75789    24152.024    24152.024    354.00348    354.00348 
Loop time of 30.8847 on 1 procs for 1000 steps with 2000 atoms

Performance: 2.798 ns/day, 8.579 hours/ns, 32.378 timesteps/s
72.6% CPU use with 1 MPI tasks x 1 OpenMP threads

MPI task timing breakdown:
Section |  min time  |  avg time  |  max time  |%varavg| %total
---------------------------------------------------------------
Pair    | 30.722     | 30.722     | 30.722     |   0.0 | 99.47
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.046597   | 0.046597   | 0.046597   |   0.0 |  0.15
Output  | 6.77e-05   | 6.77e-05   | 6.77e-05   |   0.0 |  0.00
Modify  | 0.10502    | 0.10502    | 0.10502    |   0.0 |  0.34
Other   |            | 0.01112    |            |       |  0.04

Nlocal:        2000.00 ave        2000 max        2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        4118.00 ave        4118 max        4118 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:      127962.0 ave      127962 max      127962 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 127962
Ave neighs/atom = 63.981000
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
Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    4000    -8043.845    -8043.845   -8124.4005   -8124.4005    311.75789    311.75789    24152.024    24152.024    354.00348    354.00348 
    5000    -8048.463    -8048.463    -8121.569    -8121.569    282.92751    282.92751    24162.982    24162.982   -90.119616   -90.119616 
Loop time of 31.5137 on 1 procs for 1000 steps with 2000 atoms

Performance: 2.742 ns/day, 8.754 hours/ns, 31.732 timesteps/s
71.5% CPU use with 1 MPI tasks x 1 OpenMP threads

MPI task timing breakdown:
Section |  min time  |  avg time  |  max time  |%varavg| %total
---------------------------------------------------------------
Pair    | 31.329     | 31.329     | 31.329     |   0.0 | 99.41
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.026742   | 0.026742   | 0.026742   |   0.0 |  0.08
Output  | 2.91e-05   | 2.91e-05   | 2.91e-05   |   0.0 |  0.00
Modify  | 0.14704    | 0.14704    | 0.14704    |   0.0 |  0.47
Other   |            | 0.01117    |            |       |  0.04

Nlocal:        2000.00 ave        2000 max        2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        4119.00 ave        4119 max        4119 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:      127962.0 ave      127962 max      127962 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 127962
Ave neighs/atom = 63.981000
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 = 290.933034313485, Press = -221.597693055038
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    5000    -8048.463    -8048.463    -8121.569    -8121.569    282.92751    282.92751    24162.982    24162.982   -90.119616   -90.119616 
    6000   -8045.2071   -8045.2071   -8122.8389   -8122.8389     300.4431     300.4431    24145.115    24145.115    918.15207    918.15207 
Loop time of 27.8274 on 1 procs for 1000 steps with 2000 atoms

Performance: 3.105 ns/day, 7.730 hours/ns, 35.936 timesteps/s
80.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    | 27.674     | 27.674     | 27.674     |   0.0 | 99.45
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.032866   | 0.032866   | 0.032866   |   0.0 |  0.12
Output  | 3.27e-05   | 3.27e-05   | 3.27e-05   |   0.0 |  0.00
Modify  | 0.10921    | 0.10921    | 0.10921    |   0.0 |  0.39
Other   |            | 0.01094    |            |       |  0.04

Nlocal:        2000.00 ave        2000 max        2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        4119.00 ave        4119 max        4119 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:      127942.0 ave      127942 max      127942 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 127942
Ave neighs/atom = 63.971000
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 = 292.504653275098, Press = -7.42972221210893
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    6000   -8045.2071   -8045.2071   -8122.8389   -8122.8389     300.4431     300.4431    24145.115    24145.115    918.15207    918.15207 
    7000   -8046.9158   -8046.9158   -8120.9246   -8120.9246    286.42179    286.42179    24219.553    24219.553   -4696.6846   -4696.6846 
Loop time of 29.5391 on 1 procs for 1000 steps with 2000 atoms

Performance: 2.925 ns/day, 8.205 hours/ns, 33.853 timesteps/s
76.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    | 29.312     | 29.312     | 29.312     |   0.0 | 99.23
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.066292   | 0.066292   | 0.066292   |   0.0 |  0.22
Output  | 3.14e-05   | 3.14e-05   | 3.14e-05   |   0.0 |  0.00
Modify  | 0.12972    | 0.12972    | 0.12972    |   0.0 |  0.44
Other   |            | 0.03109    |            |       |  0.11

Nlocal:        2000.00 ave        2000 max        2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        4119.00 ave        4119 max        4119 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:      127946.0 ave      127946 max      127946 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 127946
Ave neighs/atom = 63.973000
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 = 291.274344859836, Press = -21.0464029274886
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    7000   -8046.9158   -8046.9158   -8120.9246   -8120.9246    286.42179    286.42179    24219.553    24219.553   -4696.6846   -4696.6846 
    8000   -8045.4272   -8045.4272   -8120.2693   -8120.2693    289.64659    289.64659    24125.177    24125.177    3036.7385    3036.7385 
Loop time of 30.3448 on 1 procs for 1000 steps with 2000 atoms

Performance: 2.847 ns/day, 8.429 hours/ns, 32.955 timesteps/s
73.2% CPU use with 1 MPI tasks x 1 OpenMP threads

MPI task timing breakdown:
Section |  min time  |  avg time  |  max time  |%varavg| %total
---------------------------------------------------------------
Pair    | 30.151     | 30.151     | 30.151     |   0.0 | 99.36
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.027426   | 0.027426   | 0.027426   |   0.0 |  0.09
Output  | 3.12e-05   | 3.12e-05   | 3.12e-05   |   0.0 |  0.00
Modify  | 0.1554     | 0.1554     | 0.1554     |   0.0 |  0.51
Other   |            | 0.01116    |            |       |  0.04

Nlocal:        2000.00 ave        2000 max        2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        4119.00 ave        4119 max        4119 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:      127958.0 ave      127958 max      127958 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 127958
Ave neighs/atom = 63.979000
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 = 291.943267328621, Press = 6.7229657832848
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    8000   -8045.4272   -8045.4272   -8120.2693   -8120.2693    289.64659    289.64659    24125.177    24125.177    3036.7385    3036.7385 
    9000   -8044.3937   -8044.3937    -8119.495    -8119.495    290.64971    290.64971    24183.023    24183.023   -1318.1243   -1318.1243 
Loop time of 27.5505 on 1 procs for 1000 steps with 2000 atoms

Performance: 3.136 ns/day, 7.653 hours/ns, 36.297 timesteps/s
80.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    | 27.394     | 27.394     | 27.394     |   0.0 | 99.43
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.026401   | 0.026401   | 0.026401   |   0.0 |  0.10
Output  | 2.96e-05   | 2.96e-05   | 2.96e-05   |   0.0 |  0.00
Modify  | 0.11912    | 0.11912    | 0.11912    |   0.0 |  0.43
Other   |            | 0.01105    |            |       |  0.04

Nlocal:        2000.00 ave        2000 max        2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        4119.00 ave        4119 max        4119 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:      127952.0 ave      127952 max      127952 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 127952
Ave neighs/atom = 63.976000
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 = 292.087782066433, Press = -11.2112060339322
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    9000   -8044.3937   -8044.3937    -8119.495    -8119.495    290.64971    290.64971    24183.023    24183.023   -1318.1243   -1318.1243 
   10000   -8046.9057   -8046.9057   -8122.7984   -8122.7984    293.71268    293.71268     24141.53     24141.53    1259.8072    1259.8072 
Loop time of 29.0496 on 1 procs for 1000 steps with 2000 atoms

Performance: 2.974 ns/day, 8.069 hours/ns, 34.424 timesteps/s
76.6% CPU use with 1 MPI tasks x 1 OpenMP threads

MPI task timing breakdown:
Section |  min time  |  avg time  |  max time  |%varavg| %total
---------------------------------------------------------------
Pair    | 28.892     | 28.892     | 28.892     |   0.0 | 99.46
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.026337   | 0.026337   | 0.026337   |   0.0 |  0.09
Output  | 2.72e-05   | 2.72e-05   | 2.72e-05   |   0.0 |  0.00
Modify  | 0.12057    | 0.12057    | 0.12057    |   0.0 |  0.42
Other   |            | 0.01097    |            |       |  0.04

Nlocal:        2000.00 ave        2000 max        2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        4119.00 ave        4119 max        4119 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:      127956.0 ave      127956 max      127956 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 127956
Ave neighs/atom = 63.978000
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 = 292.57709586908, Press = 3.67624876794802
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   10000   -8046.9057   -8046.9057   -8122.7984   -8122.7984    293.71268    293.71268     24141.53     24141.53    1259.8072    1259.8072 
   11000   -8043.4838   -8043.4838   -8119.1161   -8119.1161    292.70459    292.70459    24186.146    24186.146   -1451.2318   -1451.2318 
Loop time of 29.7528 on 1 procs for 1000 steps with 2000 atoms

Performance: 2.904 ns/day, 8.265 hours/ns, 33.610 timesteps/s
74.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    | 29.581     | 29.581     | 29.581     |   0.0 | 99.42
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.047331   | 0.047331   | 0.047331   |   0.0 |  0.16
Output  | 4.66e-05   | 4.66e-05   | 4.66e-05   |   0.0 |  0.00
Modify  | 0.11327    | 0.11327    | 0.11327    |   0.0 |  0.38
Other   |            | 0.01092    |            |       |  0.04

Nlocal:        2000.00 ave        2000 max        2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        4119.00 ave        4119 max        4119 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:      127950.0 ave      127950 max      127950 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 127950
Ave neighs/atom = 63.975000
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 = 293.007093505434, Press = -4.92190832287155
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   11000   -8043.4838   -8043.4838   -8119.1161   -8119.1161    292.70459    292.70459    24186.146    24186.146   -1451.2318   -1451.2318 
   12000   -8048.8787   -8048.8787   -8122.9156   -8122.9156     286.5304     286.5304    24135.174    24135.174    1786.3189    1786.3189 
Loop time of 29.1387 on 1 procs for 1000 steps with 2000 atoms

Performance: 2.965 ns/day, 8.094 hours/ns, 34.319 timesteps/s
77.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    | 28.932     | 28.932     | 28.932     |   0.0 | 99.29
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.046372   | 0.046372   | 0.046372   |   0.0 |  0.16
Output  | 3.1e-05    | 3.1e-05    | 3.1e-05    |   0.0 |  0.00
Modify  | 0.14966    | 0.14966    | 0.14966    |   0.0 |  0.51
Other   |            | 0.01111    |            |       |  0.04

Nlocal:        2000.00 ave        2000 max        2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        4119.00 ave        4119 max        4119 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:      127948.0 ave      127948 max      127948 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 127948
Ave neighs/atom = 63.974000
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 = 293.048176719083, Press = 6.53830070638005
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   12000   -8048.8787   -8048.8787   -8122.9156   -8122.9156     286.5304     286.5304    24135.174    24135.174    1786.3189    1786.3189 
   13000   -8045.4818   -8045.4818   -8121.3733   -8121.3733    293.70784    293.70784    24179.609    24179.609   -1489.2901   -1489.2901 
Loop time of 29.4819 on 1 procs for 1000 steps with 2000 atoms

Performance: 2.931 ns/day, 8.189 hours/ns, 33.919 timesteps/s
75.4% CPU use with 1 MPI tasks x 1 OpenMP threads

MPI task timing breakdown:
Section |  min time  |  avg time  |  max time  |%varavg| %total
---------------------------------------------------------------
Pair    | 29.315     | 29.315     | 29.315     |   0.0 | 99.43
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.046218   | 0.046218   | 0.046218   |   0.0 |  0.16
Output  | 2.71e-05   | 2.71e-05   | 2.71e-05   |   0.0 |  0.00
Modify  | 0.10932    | 0.10932    | 0.10932    |   0.0 |  0.37
Other   |            | 0.01114    |            |       |  0.04

Nlocal:        2000.00 ave        2000 max        2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        4119.00 ave        4119 max        4119 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:      127946.0 ave      127946 max      127946 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 127946
Ave neighs/atom = 63.973000
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 = 292.828898293059, Press = -9.38430785333928
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   13000   -8045.4818   -8045.4818   -8121.3733   -8121.3733    293.70784    293.70784    24179.609    24179.609   -1489.2901   -1489.2901 
   14000   -8046.8863   -8046.8863   -8123.3613   -8123.3613    295.96584    295.96584    24153.633    24153.633    366.10328    366.10328 
Loop time of 29.9442 on 1 procs for 1000 steps with 2000 atoms

Performance: 2.885 ns/day, 8.318 hours/ns, 33.395 timesteps/s
74.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    | 29.798     | 29.798     | 29.798     |   0.0 | 99.51
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.026185   | 0.026185   | 0.026185   |   0.0 |  0.09
Output  | 3.1e-05    | 3.1e-05    | 3.1e-05    |   0.0 |  0.00
Modify  | 0.10904    | 0.10904    | 0.10904    |   0.0 |  0.36
Other   |            | 0.01086    |            |       |  0.04

Nlocal:        2000.00 ave        2000 max        2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        4119.00 ave        4119 max        4119 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:      127938.0 ave      127938 max      127938 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 127938
Ave neighs/atom = 63.969000
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 = 292.809427374262, Press = 6.1975491886284
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   14000   -8046.8863   -8046.8863   -8123.3613   -8123.3613    295.96584    295.96584    24153.633    24153.633    366.10328    366.10328 
   15000    -8046.114    -8046.114   -8123.2918   -8123.2918    298.68603    298.68603    24160.254    24160.254   -158.31002   -158.31002 
Loop time of 30.1244 on 1 procs for 1000 steps with 2000 atoms

Performance: 2.868 ns/day, 8.368 hours/ns, 33.196 timesteps/s
74.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    | 29.914     | 29.914     | 29.914     |   0.0 | 99.30
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.066265   | 0.066265   | 0.066265   |   0.0 |  0.22
Output  | 3.18e-05   | 3.18e-05   | 3.18e-05   |   0.0 |  0.00
Modify  | 0.13321    | 0.13321    | 0.13321    |   0.0 |  0.44
Other   |            | 0.01103    |            |       |  0.04

Nlocal:        2000.00 ave        2000 max        2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        4119.00 ave        4119 max        4119 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:      127954.0 ave      127954 max      127954 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 127954
Ave neighs/atom = 63.977000
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 = 293.025428976144, Press = -6.24726784344809
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   15000    -8046.114    -8046.114   -8123.2918   -8123.2918    298.68603    298.68603    24160.254    24160.254   -158.31002   -158.31002 
   16000   -8042.0269   -8042.0269   -8118.3651   -8118.3651    295.43667    295.43667    24179.039    24179.039   -600.32936   -600.32936 
Loop time of 27.5989 on 1 procs for 1000 steps with 2000 atoms

Performance: 3.131 ns/day, 7.666 hours/ns, 36.233 timesteps/s
81.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    | 27.452     | 27.452     | 27.452     |   0.0 | 99.47
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.046349   | 0.046349   | 0.046349   |   0.0 |  0.17
Output  | 3.07e-05   | 3.07e-05   | 3.07e-05   |   0.0 |  0.00
Modify  | 0.089394   | 0.089394   | 0.089394   |   0.0 |  0.32
Other   |            | 0.01098    |            |       |  0.04

Nlocal:        2000.00 ave        2000 max        2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        4119.00 ave        4119 max        4119 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:      127958.0 ave      127958 max      127958 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 127958
Ave neighs/atom = 63.979000
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 = 293.136334246919, Press = 5.42029423612746
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   16000   -8042.0269   -8042.0269   -8118.3651   -8118.3651    295.43667    295.43667    24179.039    24179.039   -600.32936   -600.32936 
   17000   -8047.3662   -8047.3662   -8123.7217   -8123.7217    295.50352    295.50352    24148.841    24148.841    548.47971    548.47971 
Loop time of 26.4477 on 1 procs for 1000 steps with 2000 atoms

Performance: 3.267 ns/day, 7.347 hours/ns, 37.810 timesteps/s
83.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    | 26.321     | 26.321     | 26.321     |   0.0 | 99.52
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.026246   | 0.026246   | 0.026246   |   0.0 |  0.10
Output  | 3.11e-05   | 3.11e-05   | 3.11e-05   |   0.0 |  0.00
Modify  | 0.089221   | 0.089221   | 0.089221   |   0.0 |  0.34
Other   |            | 0.01098    |            |       |  0.04

Nlocal:        2000.00 ave        2000 max        2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        4119.00 ave        4119 max        4119 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:      127912.0 ave      127912 max      127912 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 127912
Ave neighs/atom = 63.956000
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 = 293.284214103881, Press = -4.94236332418159
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   17000   -8047.3662   -8047.3662   -8123.7217   -8123.7217    295.50352    295.50352    24148.841    24148.841    548.47971    548.47971 
   18000   -8043.4466   -8043.4466   -8120.6137   -8120.6137     298.6445     298.6445    24173.959    24173.959   -667.93738   -667.93738 
Loop time of 28.6126 on 1 procs for 1000 steps with 2000 atoms

Performance: 3.020 ns/day, 7.948 hours/ns, 34.950 timesteps/s
78.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    | 28.466     | 28.466     | 28.466     |   0.0 | 99.49
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.026243   | 0.026243   | 0.026243   |   0.0 |  0.09
Output  | 3.09e-05   | 3.09e-05   | 3.09e-05   |   0.0 |  0.00
Modify  | 0.10964    | 0.10964    | 0.10964    |   0.0 |  0.38
Other   |            | 0.01102    |            |       |  0.04

Nlocal:        2000.00 ave        2000 max        2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        4119.00 ave        4119 max        4119 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:      127970.0 ave      127970 max      127970 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 127970
Ave neighs/atom = 63.985000
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 = 293.434899677863, Press = 1.46897163689878
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   18000   -8043.4466   -8043.4466   -8120.6137   -8120.6137     298.6445     298.6445    24173.959    24173.959   -667.93738   -667.93738 
   19000   -8046.1152   -8046.1152   -8122.2464   -8122.2464    294.63544    294.63544     24157.01     24157.01    258.00467    258.00467 
Loop time of 28.1394 on 1 procs for 1000 steps with 2000 atoms

Performance: 3.070 ns/day, 7.817 hours/ns, 35.537 timesteps/s
78.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    | 27.951     | 27.951     | 27.951     |   0.0 | 99.33
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.046353   | 0.046353   | 0.046353   |   0.0 |  0.16
Output  | 3.12e-05   | 3.12e-05   | 3.12e-05   |   0.0 |  0.00
Modify  | 0.13099    | 0.13099    | 0.13099    |   0.0 |  0.47
Other   |            | 0.01092    |            |       |  0.04

Nlocal:        2000.00 ave        2000 max        2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        4119.00 ave        4119 max        4119 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:      127946.0 ave      127946 max      127946 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 127946
Ave neighs/atom = 63.973000
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 = 293.299672698174, Press = -1.87441716698394
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   19000   -8046.1152   -8046.1152   -8122.2464   -8122.2464    294.63544    294.63544     24157.01     24157.01    258.00467    258.00467 
   20000   -8047.4297   -8047.4297    -8121.065    -8121.065    284.97644    284.97644     24162.52     24162.52   -144.24283   -144.24283 
Loop time of 27.6927 on 1 procs for 1000 steps with 2000 atoms

Performance: 3.120 ns/day, 7.692 hours/ns, 36.111 timesteps/s
80.7% CPU use with 1 MPI tasks x 1 OpenMP threads

MPI task timing breakdown:
Section |  min time  |  avg time  |  max time  |%varavg| %total
---------------------------------------------------------------
Pair    | 27.526     | 27.526     | 27.526     |   0.0 | 99.40
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.046376   | 0.046376   | 0.046376   |   0.0 |  0.17
Output  | 2.46e-05   | 2.46e-05   | 2.46e-05   |   0.0 |  0.00
Modify  | 0.10958    | 0.10958    | 0.10958    |   0.0 |  0.40
Other   |            | 0.01092    |            |       |  0.04

Nlocal:        2000.00 ave        2000 max        2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        4119.00 ave        4119 max        4119 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:      127958.0 ave      127958 max      127958 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 127958
Ave neighs/atom = 63.979000
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 = 293.18925959642, Press = 0.290237536148372
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   20000   -8047.4297   -8047.4297    -8121.065    -8121.065    284.97644    284.97644     24162.52     24162.52   -144.24283   -144.24283 
   21000   -8043.8321   -8043.8321   -8121.7866   -8121.7866    301.69195    301.69195    24135.536    24135.536    2161.4163    2161.4163 
Loop time of 26.0877 on 1 procs for 1000 steps with 2000 atoms

Performance: 3.312 ns/day, 7.247 hours/ns, 38.332 timesteps/s
85.2% CPU use with 1 MPI tasks x 1 OpenMP threads

MPI task timing breakdown:
Section |  min time  |  avg time  |  max time  |%varavg| %total
---------------------------------------------------------------
Pair    | 25.935     | 25.935     | 25.935     |   0.0 | 99.42
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.026263   | 0.026263   | 0.026263   |   0.0 |  0.10
Output  | 6.31e-05   | 6.31e-05   | 6.31e-05   |   0.0 |  0.00
Modify  | 0.11007    | 0.11007    | 0.11007    |   0.0 |  0.42
Other   |            | 0.01585    |            |       |  0.06

Nlocal:        2000.00 ave        2000 max        2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        4119.00 ave        4119 max        4119 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:      127942.0 ave      127942 max      127942 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 127942
Ave neighs/atom = 63.971000
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 = 293.12246526595, Press = -0.774692915025372
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   21000   -8043.8321   -8043.8321   -8121.7866   -8121.7866    301.69195    301.69195    24135.536    24135.536    2161.4163    2161.4163 
   22000   -8044.5097   -8044.5097   -8119.7973   -8119.7973    291.37115    291.37115    24216.217    24216.217   -3879.9087   -3879.9087 
Loop time of 26.2606 on 1 procs for 1000 steps with 2000 atoms

Performance: 3.290 ns/day, 7.295 hours/ns, 38.080 timesteps/s
85.2% CPU use with 1 MPI tasks x 1 OpenMP threads

MPI task timing breakdown:
Section |  min time  |  avg time  |  max time  |%varavg| %total
---------------------------------------------------------------
Pair    | 26.105     | 26.105     | 26.105     |   0.0 | 99.41
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.026468   | 0.026468   | 0.026468   |   0.0 |  0.10
Output  | 3.07e-05   | 3.07e-05   | 3.07e-05   |   0.0 |  0.00
Modify  | 0.11776    | 0.11776    | 0.11776    |   0.0 |  0.45
Other   |            | 0.01113    |            |       |  0.04

Nlocal:        2000.00 ave        2000 max        2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        4119.00 ave        4119 max        4119 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:      127952.0 ave      127952 max      127952 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 127952
Ave neighs/atom = 63.976000
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 = 293.152263064806, Press = -1.8575580917873
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   22000   -8044.5097   -8044.5097   -8119.7973   -8119.7973    291.37115    291.37115    24216.217    24216.217   -3879.9087   -3879.9087 
   23000   -8045.9791   -8045.9791   -8123.1171   -8123.1171    298.53196    298.53196    24130.608    24130.608    2245.8801    2245.8801 
Loop time of 25.9689 on 1 procs for 1000 steps with 2000 atoms

Performance: 3.327 ns/day, 7.214 hours/ns, 38.508 timesteps/s
85.5% CPU use with 1 MPI tasks x 1 OpenMP threads

MPI task timing breakdown:
Section |  min time  |  avg time  |  max time  |%varavg| %total
---------------------------------------------------------------
Pair    | 25.842     | 25.842     | 25.842     |   0.0 | 99.51
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.026307   | 0.026307   | 0.026307   |   0.0 |  0.10
Output  | 3.25e-05   | 3.25e-05   | 3.25e-05   |   0.0 |  0.00
Modify  | 0.089554   | 0.089554   | 0.089554   |   0.0 |  0.34
Other   |            | 0.01107    |            |       |  0.04

Nlocal:        2000.00 ave        2000 max        2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        4119.00 ave        4119 max        4119 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:      127926.0 ave      127926 max      127926 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 127926
Ave neighs/atom = 63.963000
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 = 293.337134135036, Press = 0.978892357478476
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   23000   -8045.9791   -8045.9791   -8123.1171   -8123.1171    298.53196    298.53196    24130.608    24130.608    2245.8801    2245.8801 
   24000   -8044.9632   -8044.9632   -8120.9292   -8120.9292    293.99631    293.99631    24175.857    24175.857   -1020.7412   -1020.7412 
Loop time of 25.944 on 1 procs for 1000 steps with 2000 atoms

Performance: 3.330 ns/day, 7.207 hours/ns, 38.545 timesteps/s
85.7% CPU use with 1 MPI tasks x 1 OpenMP threads

MPI task timing breakdown:
Section |  min time  |  avg time  |  max time  |%varavg| %total
---------------------------------------------------------------
Pair    | 25.796     | 25.796     | 25.796     |   0.0 | 99.43
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.046242   | 0.046242   | 0.046242   |   0.0 |  0.18
Output  | 2.55e-05   | 2.55e-05   | 2.55e-05   |   0.0 |  0.00
Modify  | 0.090246   | 0.090246   | 0.090246   |   0.0 |  0.35
Other   |            | 0.01102    |            |       |  0.04

Nlocal:        2000.00 ave        2000 max        2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        4119.00 ave        4119 max        4119 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:      127940.0 ave      127940 max      127940 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 127940
Ave neighs/atom = 63.970000
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 = 293.321914027374, Press = -1.97383432233629
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   24000   -8044.9632   -8044.9632   -8120.9292   -8120.9292    293.99631    293.99631    24175.857    24175.857   -1020.7412   -1020.7412 
   25000   -8050.0806   -8050.0806   -8124.0114   -8124.0114    286.11977    286.11977    24149.806    24149.806    265.48425    265.48425 
Loop time of 24.2613 on 1 procs for 1000 steps with 2000 atoms

Performance: 3.561 ns/day, 6.739 hours/ns, 41.218 timesteps/s
90.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    | 24.134     | 24.134     | 24.134     |   0.0 | 99.48
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.026116   | 0.026116   | 0.026116   |   0.0 |  0.11
Output  | 2.52e-05   | 2.52e-05   | 2.52e-05   |   0.0 |  0.00
Modify  | 0.090042   | 0.090042   | 0.090042   |   0.0 |  0.37
Other   |            | 0.01086    |            |       |  0.04

Nlocal:        2000.00 ave        2000 max        2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        4119.00 ave        4119 max        4119 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:      127942.0 ave      127942 max      127942 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 127942
Ave neighs/atom = 63.971000
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 = 293.315057126884, Press = 0.222349815794548
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   25000   -8050.0806   -8050.0806   -8124.0114   -8124.0114    286.11977    286.11977    24149.806    24149.806    265.48425    265.48425 
   26000    -8045.322    -8045.322   -8121.9769   -8121.9769    296.66235    296.66235    24169.736    24169.736   -632.96031   -632.96031 
Loop time of 24.1705 on 1 procs for 1000 steps with 2000 atoms

Performance: 3.575 ns/day, 6.714 hours/ns, 41.373 timesteps/s
92.5% CPU use with 1 MPI tasks x 1 OpenMP threads

MPI task timing breakdown:
Section |  min time  |  avg time  |  max time  |%varavg| %total
---------------------------------------------------------------
Pair    | 24.044     | 24.044     | 24.044     |   0.0 | 99.47
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.026323   | 0.026323   | 0.026323   |   0.0 |  0.11
Output  | 3.1e-05    | 3.1e-05    | 3.1e-05    |   0.0 |  0.00
Modify  | 0.089512   | 0.089512   | 0.089512   |   0.0 |  0.37
Other   |            | 0.01108    |            |       |  0.05

Nlocal:        2000.00 ave        2000 max        2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        4119.00 ave        4119 max        4119 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:      127956.0 ave      127956 max      127956 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 127956
Ave neighs/atom = 63.978000
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 = 293.174968968998, Press = -1.53544862856665
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   26000    -8045.322    -8045.322   -8121.9769   -8121.9769    296.66235    296.66235    24169.736    24169.736   -632.96031   -632.96031 
   27000   -8047.4322   -8047.4322   -8124.1262   -8124.1262    296.81376    296.81376     24162.84     24162.84   -425.12978   -425.12978 
Loop time of 24.9854 on 1 procs for 1000 steps with 2000 atoms

Performance: 3.458 ns/day, 6.940 hours/ns, 40.023 timesteps/s
89.4% CPU use with 1 MPI tasks x 1 OpenMP threads

MPI task timing breakdown:
Section |  min time  |  avg time  |  max time  |%varavg| %total
---------------------------------------------------------------
Pair    | 24.798     | 24.798     | 24.798     |   0.0 | 99.25
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.026377   | 0.026377   | 0.026377   |   0.0 |  0.11
Output  | 3.44e-05   | 3.44e-05   | 3.44e-05   |   0.0 |  0.00
Modify  | 0.13007    | 0.13007    | 0.13007    |   0.0 |  0.52
Other   |            | 0.03101    |            |       |  0.12

Nlocal:        2000.00 ave        2000 max        2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:        4119.00 ave        4119 max        4119 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:      127930.0 ave      127930 max      127930 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 127930
Ave neighs/atom = 63.965000
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_T293.15.out"  else  "print 'not_converged'  file output/vol_T293.15.out"
print '${V}'  file output/vol_T293.15.out
24160.8828960695

print "LAMMPS calculation completed"
LAMMPS calculation completed
quit 0