# Variables that can be adjusted by kim-lammps-preprocessor to switch unit sets for
# Simulator Models
variable _u_distance equal 1.0
variable _u_energy equal 1.0
variable _u_mass equal 1.0
variable _u_time equal 1.0
variable _u_pressure equal 1.0
variable _u_temperature equal 1.0

# This line may be swapped out by kim-lammps-preprocessor if running against a Simulator
# Model whose atom_style is not 'atomic'
atom_style atomic

# 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 5.7913857132196505*${_u_distance}
variable latticeconst_converted equal 5.7913857132196505*1
lattice       bcc  ${latticeconst_converted}
lattice       bcc  5.79138571321965
Lattice spacing in x,y,z = 5.79139 5.79139 5.79139
region        simbox  block 0 10 0 10 0 10 units lattice
create_box    1 simbox
Created orthogonal box = (0 0 0) to (57.9139 57.9139 57.9139)
  1 by 1 by 1 MPI processor grid
create_atoms  1 box
Created 2000 atoms
  create_atoms CPU = 0.000323057 secs

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

# specify which KIM Model to use
pair_style    kim Morse_Shifted_GirifalcoWeizer_1959LowCutoff_Rb__MO_754498969542_003
WARNING: KIM Model does not provide `partialParticleVirial'; virial per atom will be zero (src/KIM/pair_kim.cpp:1089)
pair_coeff    * * Rb
mass          1 ${mass_converted}
mass          1 85.4678

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

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

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

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

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

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

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

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

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

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

# Run a simulation for at most 2000*1000 timesteps. At each 1000th time step, check
# whether the temperature and pressure have converged. If yes, break.
label top
variable a loop 2000
run 1000
Neighbor list info ...
  update every 1 steps, delay 10 steps, check yes
  max neighbors/atom: 2000, page size: 100000
  master list distance cutoff = 21.5283
  ghost atom cutoff = 21.5283
  binsize = 10.7642, bins = 6 6 6
  1 neighbor lists, perpetual/occasional/extra = 1 0 0
  (1) pair kim, perpetual
      attributes: full, newton off, cut 21.5283
      pair build: full/bin/atomonly
      stencil: full/bin/3d
      bin: standard
Per MPI rank memory allocation (min/avg/max) = 6.658 | 6.658 | 6.658 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
       0   -1497.9773   -1497.9773    -1563.389    -1563.389       253.15       253.15    194243.94    194243.94    359.68842    359.68842 
    1000   -1414.0099   -1414.0099   -1481.0999   -1481.0999    259.64552    259.64552    206749.12    206749.12   -9.8475135   -9.8475135 
Loop time of 43.3621 on 1 procs for 1000 steps with 2000 atoms

Performance: 1.993 ns/day, 12.045 hours/ns, 23.062 timesteps/s
46.3% CPU use with 1 MPI tasks x 1 OpenMP threads

MPI task timing breakdown:
Section |  min time  |  avg time  |  max time  |%varavg| %total
---------------------------------------------------------------
Pair    | 42.699     | 42.699     | 42.699     |   0.0 | 98.47
Neigh   | 0.31394    | 0.31394    | 0.31394    |   0.0 |  0.72
Comm    | 0.13105    | 0.13105    | 0.13105    |   0.0 |  0.30
Output  | 3.6001e-05 | 3.6001e-05 | 3.6001e-05 |   0.0 |  0.00
Modify  | 0.16429    | 0.16429    | 0.16429    |   0.0 |  0.38
Other   |            | 0.05361    |            |       |  0.12

Nlocal:    2000 ave 2000 max 2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    8645 ave 8645 max 8645 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Neighs:    0 ave 0 max 0 min
Histogram: 1 0 0 0 0 0 0 0 0 0
FullNghs:  797492 ave 797492 max 797492 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 797492
Ave neighs/atom = 398.746
Neighbor list builds = 7
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) = 6.662 | 6.662 | 6.662 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    1000   -1414.0099   -1414.0099   -1481.0999   -1481.0999    259.64552    259.64552    206749.12    206749.12   -9.8475135   -9.8475135 
    2000   -1428.0223   -1428.0223   -1494.8278   -1494.8278    258.54389    258.54389    205809.79    205809.79   -31.617093   -31.617093 
Loop time of 42.8447 on 1 procs for 1000 steps with 2000 atoms

Performance: 2.017 ns/day, 11.901 hours/ns, 23.340 timesteps/s
45.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    | 42.252     | 42.252     | 42.252     |   0.0 | 98.62
Neigh   | 0.29213    | 0.29213    | 0.29213    |   0.0 |  0.68
Comm    | 0.19277    | 0.19277    | 0.19277    |   0.0 |  0.45
Output  | 5.6982e-05 | 5.6982e-05 | 5.6982e-05 |   0.0 |  0.00
Modify  | 0.093654   | 0.093654   | 0.093654   |   0.0 |  0.22
Other   |            | 0.01373    |            |       |  0.03

Nlocal:    2000 ave 2000 max 2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    8670 ave 8670 max 8670 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Neighs:    0 ave 0 max 0 min
Histogram: 1 0 0 0 0 0 0 0 0 0
FullNghs:  800662 ave 800662 max 800662 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 800662
Ave neighs/atom = 400.331
Neighbor list builds = 6
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) = 6.662 | 6.662 | 6.662 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    2000   -1428.0223   -1428.0223   -1494.8278   -1494.8278    258.54389    258.54389    205809.79    205809.79   -31.617093   -31.617093 
    3000   -1428.3587   -1428.3587   -1491.8202   -1491.8202    245.60259    245.60259    206188.37    206188.37   -50.913838   -50.913838 
Loop time of 45.3318 on 1 procs for 1000 steps with 2000 atoms

Performance: 1.906 ns/day, 12.592 hours/ns, 22.060 timesteps/s
43.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    | 44.622     | 44.622     | 44.622     |   0.0 | 98.43
Neigh   | 0.26352    | 0.26352    | 0.26352    |   0.0 |  0.58
Comm    | 0.14773    | 0.14773    | 0.14773    |   0.0 |  0.33
Output  | 3.0994e-05 | 3.0994e-05 | 3.0994e-05 |   0.0 |  0.00
Modify  | 0.24504    | 0.24504    | 0.24504    |   0.0 |  0.54
Other   |            | 0.05374    |            |       |  0.12

Nlocal:    2000 ave 2000 max 2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    8663 ave 8663 max 8663 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Neighs:    0 ave 0 max 0 min
Histogram: 1 0 0 0 0 0 0 0 0 0
FullNghs:  802850 ave 802850 max 802850 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 802850
Ave neighs/atom = 401.425
Neighbor list builds = 6
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) = 6.662 | 6.662 | 6.662 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    3000   -1428.3587   -1428.3587   -1491.8202   -1491.8202    245.60259    245.60259    206188.37    206188.37   -50.913838   -50.913838 
    4000   -1426.4452   -1426.4452   -1490.7312   -1490.7312    248.79318    248.79318    205273.72    205273.72    26.356279    26.356279 
Loop time of 42.1711 on 1 procs for 1000 steps with 2000 atoms

Performance: 2.049 ns/day, 11.714 hours/ns, 23.713 timesteps/s
46.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    | 41.687     | 41.687     | 41.687     |   0.0 | 98.85
Neigh   | 0.2659     | 0.2659     | 0.2659     |   0.0 |  0.63
Comm    | 0.070613   | 0.070613   | 0.070613   |   0.0 |  0.17
Output  | 3.7909e-05 | 3.7909e-05 | 3.7909e-05 |   0.0 |  0.00
Modify  | 0.13404    | 0.13404    | 0.13404    |   0.0 |  0.32
Other   |            | 0.01364    |            |       |  0.03

Nlocal:    2000 ave 2000 max 2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    8648 ave 8648 max 8648 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Neighs:    0 ave 0 max 0 min
Histogram: 1 0 0 0 0 0 0 0 0 0
FullNghs:  802762 ave 802762 max 802762 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 802762
Ave neighs/atom = 401.381
Neighbor list builds = 6
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) = 6.662 | 6.662 | 6.662 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    4000   -1426.4452   -1426.4452   -1490.7312   -1490.7312    248.79318    248.79318    205273.72    205273.72    26.356279    26.356279 
    5000   -1424.4695   -1424.4695   -1491.9114   -1491.9114    261.00736    261.00736    205395.68    205395.68    20.754797    20.754797 
Loop time of 41.2131 on 1 procs for 1000 steps with 2000 atoms

Performance: 2.096 ns/day, 11.448 hours/ns, 24.264 timesteps/s
47.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    | 40.621     | 40.621     | 40.621     |   0.0 | 98.56
Neigh   | 0.28448    | 0.28448    | 0.28448    |   0.0 |  0.69
Comm    | 0.15018    | 0.15018    | 0.15018    |   0.0 |  0.36
Output  | 2.8133e-05 | 2.8133e-05 | 2.8133e-05 |   0.0 |  0.00
Modify  | 0.14405    | 0.14405    | 0.14405    |   0.0 |  0.35
Other   |            | 0.01356    |            |       |  0.03

Nlocal:    2000 ave 2000 max 2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    8677 ave 8677 max 8677 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Neighs:    0 ave 0 max 0 min
Histogram: 1 0 0 0 0 0 0 0 0 0
FullNghs:  800338 ave 800338 max 800338 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 800338
Ave neighs/atom = 400.169
Neighbor list builds = 6
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 = 257.27890751862, Press = 12.6967483439
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 6.662 | 6.662 | 6.662 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    5000   -1424.4695   -1424.4695   -1491.9114   -1491.9114    261.00736    261.00736    205395.68    205395.68    20.754797    20.754797 
    6000   -1425.5884   -1425.5884   -1491.4015   -1491.4015    254.70367    254.70367    204746.38    204746.38    66.566672    66.566672 
Loop time of 42.9852 on 1 procs for 1000 steps with 2000 atoms

Performance: 2.010 ns/day, 11.940 hours/ns, 23.264 timesteps/s
46.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    | 42.404     | 42.404     | 42.404     |   0.0 | 98.65
Neigh   | 0.34942    | 0.34942    | 0.34942    |   0.0 |  0.81
Comm    | 0.070443   | 0.070443   | 0.070443   |   0.0 |  0.16
Output  | 2.2173e-05 | 2.2173e-05 | 2.2173e-05 |   0.0 |  0.00
Modify  | 0.1479     | 0.1479     | 0.1479     |   0.0 |  0.34
Other   |            | 0.01355    |            |       |  0.03

Nlocal:    2000 ave 2000 max 2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    8696 ave 8696 max 8696 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Neighs:    0 ave 0 max 0 min
Histogram: 1 0 0 0 0 0 0 0 0 0
FullNghs:  805182 ave 805182 max 805182 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 805182
Ave neighs/atom = 402.591
Neighbor list builds = 7
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 = 254.044432187214, Press = 1.29059036541614
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 6.662 | 6.662 | 6.662 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    6000   -1425.5884   -1425.5884   -1491.4015   -1491.4015    254.70367    254.70367    204746.38    204746.38    66.566672    66.566672 
    7000   -1428.0675   -1428.0675   -1493.6508   -1493.6508    253.81435    253.81435    204671.63    204671.63    50.261682    50.261682 
Loop time of 41.5357 on 1 procs for 1000 steps with 2000 atoms

Performance: 2.080 ns/day, 11.538 hours/ns, 24.076 timesteps/s
47.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    | 40.799     | 40.799     | 40.799     |   0.0 | 98.23
Neigh   | 0.32627    | 0.32627    | 0.32627    |   0.0 |  0.79
Comm    | 0.11063    | 0.11063    | 0.11063    |   0.0 |  0.27
Output  | 2.6941e-05 | 2.6941e-05 | 2.6941e-05 |   0.0 |  0.00
Modify  | 0.25593    | 0.25593    | 0.25593    |   0.0 |  0.62
Other   |            | 0.04363    |            |       |  0.11

Nlocal:    2000 ave 2000 max 2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    8705 ave 8705 max 8705 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Neighs:    0 ave 0 max 0 min
Histogram: 1 0 0 0 0 0 0 0 0 0
FullNghs:  799802 ave 799802 max 799802 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 799802
Ave neighs/atom = 399.901
Neighbor list builds = 6
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 253.22509696623, Press = -2.31861063347147
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 6.662 | 6.662 | 6.662 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    7000   -1428.0675   -1428.0675   -1493.6508   -1493.6508    253.81435    253.81435    204671.63    204671.63    50.261682    50.261682 
    8000   -1427.1303   -1427.1303   -1491.8306   -1491.8306      250.397      250.397    206056.99    206056.99   -43.681614   -43.681614 
Loop time of 42.2075 on 1 procs for 1000 steps with 2000 atoms

Performance: 2.047 ns/day, 11.724 hours/ns, 23.692 timesteps/s
47.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    | 41.547     | 41.547     | 41.547     |   0.0 | 98.43
Neigh   | 0.31697    | 0.31697    | 0.31697    |   0.0 |  0.75
Comm    | 0.20229    | 0.20229    | 0.20229    |   0.0 |  0.48
Output  | 2.7895e-05 | 2.7895e-05 | 2.7895e-05 |   0.0 |  0.00
Modify  | 0.12774    | 0.12774    | 0.12774    |   0.0 |  0.30
Other   |            | 0.01362    |            |       |  0.03

Nlocal:    2000 ave 2000 max 2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    8683 ave 8683 max 8683 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Neighs:    0 ave 0 max 0 min
Histogram: 1 0 0 0 0 0 0 0 0 0
FullNghs:  798408 ave 798408 max 798408 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 798408
Ave neighs/atom = 399.204
Neighbor list builds = 7
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 253.150628758323, Press = -1.579708962662
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 6.662 | 6.662 | 6.662 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    8000   -1427.1303   -1427.1303   -1491.8306   -1491.8306      250.397      250.397    206056.99    206056.99   -43.681614   -43.681614 
    9000   -1428.1149   -1428.1149   -1493.7342   -1493.7342    253.95345    253.95345    206473.78    206473.78   -82.516418   -82.516418 
Loop time of 44.1947 on 1 procs for 1000 steps with 2000 atoms

Performance: 1.955 ns/day, 12.276 hours/ns, 22.627 timesteps/s
44.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    | 43.589     | 43.589     | 43.589     |   0.0 | 98.63
Neigh   | 0.32401    | 0.32401    | 0.32401    |   0.0 |  0.73
Comm    | 0.05058    | 0.05058    | 0.05058    |   0.0 |  0.11
Output  | 2.9087e-05 | 2.9087e-05 | 2.9087e-05 |   0.0 |  0.00
Modify  | 0.19744    | 0.19744    | 0.19744    |   0.0 |  0.45
Other   |            | 0.03381    |            |       |  0.08

Nlocal:    2000 ave 2000 max 2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    8687 ave 8687 max 8687 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Neighs:    0 ave 0 max 0 min
Histogram: 1 0 0 0 0 0 0 0 0 0
FullNghs:  795998 ave 795998 max 795998 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 795998
Ave neighs/atom = 397.999
Neighbor list builds = 7
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 252.980838216196, Press = -0.207934212188254
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 6.662 | 6.662 | 6.662 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    9000   -1428.1149   -1428.1149   -1493.7342   -1493.7342    253.95345    253.95345    206473.78    206473.78   -82.516418   -82.516418 
   10000   -1428.9244   -1428.9244    -1493.404    -1493.404    249.54249    249.54249    205834.67    205834.67   -40.689464   -40.689464 
Loop time of 43.9961 on 1 procs for 1000 steps with 2000 atoms

Performance: 1.964 ns/day, 12.221 hours/ns, 22.729 timesteps/s
44.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    | 43.378     | 43.378     | 43.378     |   0.0 | 98.60
Neigh   | 0.30583    | 0.30583    | 0.30583    |   0.0 |  0.70
Comm    | 0.16084    | 0.16084    | 0.16084    |   0.0 |  0.37
Output  | 2.7895e-05 | 2.7895e-05 | 2.7895e-05 |   0.0 |  0.00
Modify  | 0.13745    | 0.13745    | 0.13745    |   0.0 |  0.31
Other   |            | 0.01352    |            |       |  0.03

Nlocal:    2000 ave 2000 max 2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    8701 ave 8701 max 8701 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Neighs:    0 ave 0 max 0 min
Histogram: 1 0 0 0 0 0 0 0 0 0
FullNghs:  798876 ave 798876 max 798876 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 798876
Ave neighs/atom = 399.438
Neighbor list builds = 7
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 253.125599406284, Press = 0.260115807877841
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 6.662 | 6.662 | 6.662 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   10000   -1428.9244   -1428.9244    -1493.404    -1493.404    249.54249    249.54249    205834.67    205834.67   -40.689464   -40.689464 
   11000   -1426.1885   -1426.1885   -1492.3043   -1492.3043    255.87485    255.87485    206097.99    206097.99   -46.343623   -46.343623 
Loop time of 41.0314 on 1 procs for 1000 steps with 2000 atoms

Performance: 2.106 ns/day, 11.398 hours/ns, 24.372 timesteps/s
48.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    | 40.409     | 40.409     | 40.409     |   0.0 | 98.48
Neigh   | 0.35128    | 0.35128    | 0.35128    |   0.0 |  0.86
Comm    | 0.091281   | 0.091281   | 0.091281   |   0.0 |  0.22
Output  | 6.6042e-05 | 6.6042e-05 | 6.6042e-05 |   0.0 |  0.00
Modify  | 0.16604    | 0.16604    | 0.16604    |   0.0 |  0.40
Other   |            | 0.01362    |            |       |  0.03

Nlocal:    2000 ave 2000 max 2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    8683 ave 8683 max 8683 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Neighs:    0 ave 0 max 0 min
Histogram: 1 0 0 0 0 0 0 0 0 0
FullNghs:  798024 ave 798024 max 798024 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 798024
Ave neighs/atom = 399.012
Neighbor list builds = 7
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 252.833614547772, Press = 0.148518386314887
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 6.662 | 6.662 | 6.662 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   11000   -1426.1885   -1426.1885   -1492.3043   -1492.3043    255.87485    255.87485    206097.99    206097.99   -46.343623   -46.343623 
   12000   -1424.8488   -1424.8488   -1490.8231   -1490.8231    255.32718    255.32718    205914.74    205914.74   -21.176255   -21.176255 
Loop time of 42.0793 on 1 procs for 1000 steps with 2000 atoms

Performance: 2.053 ns/day, 11.689 hours/ns, 23.765 timesteps/s
46.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    | 41.586     | 41.586     | 41.586     |   0.0 | 98.83
Neigh   | 0.22946    | 0.22946    | 0.22946    |   0.0 |  0.55
Comm    | 0.092106   | 0.092106   | 0.092106   |   0.0 |  0.22
Output  | 2.6941e-05 | 2.6941e-05 | 2.6941e-05 |   0.0 |  0.00
Modify  | 0.13808    | 0.13808    | 0.13808    |   0.0 |  0.33
Other   |            | 0.03382    |            |       |  0.08

Nlocal:    2000 ave 2000 max 2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    8686 ave 8686 max 8686 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Neighs:    0 ave 0 max 0 min
Histogram: 1 0 0 0 0 0 0 0 0 0
FullNghs:  801016 ave 801016 max 801016 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 801016
Ave neighs/atom = 400.508
Neighbor list builds = 6
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 252.844849503159, Press = 0.449241206219857
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 6.662 | 6.662 | 6.662 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   12000   -1424.8488   -1424.8488   -1490.8231   -1490.8231    255.32718    255.32718    205914.74    205914.74   -21.176255   -21.176255 
   13000   -1427.1365   -1427.1365    -1492.432    -1492.432     252.7006     252.7006    205518.04    205518.04   -5.3319786   -5.3319786 
Loop time of 39.4854 on 1 procs for 1000 steps with 2000 atoms

Performance: 2.188 ns/day, 10.968 hours/ns, 25.326 timesteps/s
50.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    | 38.94      | 38.94      | 38.94      |   0.0 | 98.62
Neigh   | 0.25248    | 0.25248    | 0.25248    |   0.0 |  0.64
Comm    | 0.070706   | 0.070706   | 0.070706   |   0.0 |  0.18
Output  | 2.6941e-05 | 2.6941e-05 | 2.6941e-05 |   0.0 |  0.00
Modify  | 0.20809    | 0.20809    | 0.20809    |   0.0 |  0.53
Other   |            | 0.0137     |            |       |  0.03

Nlocal:    2000 ave 2000 max 2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    8695 ave 8695 max 8695 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Neighs:    0 ave 0 max 0 min
Histogram: 1 0 0 0 0 0 0 0 0 0
FullNghs:  802722 ave 802722 max 802722 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 802722
Ave neighs/atom = 401.361
Neighbor list builds = 6
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 252.967584849601, Press = 0.47364010536385
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 6.662 | 6.662 | 6.662 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   13000   -1427.1365   -1427.1365    -1492.432    -1492.432     252.7006     252.7006    205518.04    205518.04   -5.3319786   -5.3319786 
   14000   -1427.0855   -1427.0855   -1491.4603   -1491.4603    249.13709    249.13709    205214.73    205214.73    18.552063    18.552063 
Loop time of 39.5698 on 1 procs for 1000 steps with 2000 atoms

Performance: 2.183 ns/day, 10.992 hours/ns, 25.272 timesteps/s
50.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    | 38.967     | 38.967     | 38.967     |   0.0 | 98.48
Neigh   | 0.27073    | 0.27073    | 0.27073    |   0.0 |  0.68
Comm    | 0.13047    | 0.13047    | 0.13047    |   0.0 |  0.33
Output  | 2.7895e-05 | 2.7895e-05 | 2.7895e-05 |   0.0 |  0.00
Modify  | 0.1881     | 0.1881     | 0.1881     |   0.0 |  0.48
Other   |            | 0.01374    |            |       |  0.03

Nlocal:    2000 ave 2000 max 2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    8690 ave 8690 max 8690 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Neighs:    0 ave 0 max 0 min
Histogram: 1 0 0 0 0 0 0 0 0 0
FullNghs:  803288 ave 803288 max 803288 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 803288
Ave neighs/atom = 401.644
Neighbor list builds = 6
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 252.936202866327, Press = 0.53223021926533
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 6.662 | 6.662 | 6.662 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   14000   -1427.0855   -1427.0855   -1491.4603   -1491.4603    249.13709    249.13709    205214.73    205214.73    18.552063    18.552063 
   15000   -1425.9896   -1425.9896   -1492.4067   -1492.4067    257.04119    257.04119    204712.04    204712.04    62.794274    62.794274 
Loop time of 39.6656 on 1 procs for 1000 steps with 2000 atoms

Performance: 2.178 ns/day, 11.018 hours/ns, 25.211 timesteps/s
49.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    | 39.004     | 39.004     | 39.004     |   0.0 | 98.33
Neigh   | 0.30067    | 0.30067    | 0.30067    |   0.0 |  0.76
Comm    | 0.14845    | 0.14845    | 0.14845    |   0.0 |  0.37
Output  | 2.9087e-05 | 2.9087e-05 | 2.9087e-05 |   0.0 |  0.00
Modify  | 0.15885    | 0.15885    | 0.15885    |   0.0 |  0.40
Other   |            | 0.05372    |            |       |  0.14

Nlocal:    2000 ave 2000 max 2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    8704 ave 8704 max 8704 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Neighs:    0 ave 0 max 0 min
Histogram: 1 0 0 0 0 0 0 0 0 0
FullNghs:  804182 ave 804182 max 804182 min
Histogram: 1 0 0 0 0 0 0 0 0 0

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

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

print "LAMMPS calculation completed"
LAMMPS calculation completed
quit 0