# 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 6.083469396829612*${_u_distance}
variable latticeconst_converted equal 6.083469396829612*1
lattice       fcc  ${latticeconst_converted}
lattice       fcc  6.08346939682961
Lattice spacing in x,y,z = 6.08347 6.08347 6.08347
region        simbox  block 0 10 0 10 0 10 units lattice
create_box    1 simbox
Created orthogonal box = (0 0 0) to (60.8347 60.8347 60.8347)
  1 by 1 by 1 MPI processor grid
create_atoms  1 box
Created 4000 atoms
  create_atoms CPU = 0.000452995 secs

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

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

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

print "Initial system volume: ${V0_metal} Angstroms^3"
Initial system volume: 225140.684924316 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
Neighbor list info ...
  update every 1 steps, delay 10 steps, check yes
  max neighbors/atom: 2000, page size: 100000
  master list distance cutoff = 17.4319
  ghost atom cutoff = 17.4319
  binsize = 8.71595, bins = 7 7 7
  1 neighbor lists, perpetual/occasional/extra = 1 0 0
  (1) pair kim, perpetual
      attributes: full, newton off, cut 17.4319
      pair build: full/bin/atomonly
      stencil: full/bin/3d
      bin: standard
Per MPI rank memory allocation (min/avg/max) = 9.433 | 9.433 | 9.433 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
       0    -6697.234    -6697.234   -6838.4284   -6838.4284       273.15       273.15    225140.68    225140.68     669.8581     669.8581 
    1000   -6540.2745   -6540.2745   -6681.9971   -6681.9971    274.17198    274.17198    230227.51    230227.51   -277.18417   -277.18417 
Loop time of 84.9203 on 1 procs for 1000 steps with 4000 atoms

Performance: 1.017 ns/day, 23.589 hours/ns, 11.776 timesteps/s
39.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    | 83.646     | 83.646     | 83.646     |   0.0 | 98.50
Neigh   | 0.56695    | 0.56695    | 0.56695    |   0.0 |  0.67
Comm    | 0.23034    | 0.23034    | 0.23034    |   0.0 |  0.27
Output  | 4.5061e-05 | 4.5061e-05 | 4.5061e-05 |   0.0 |  0.00
Modify  | 0.41526    | 0.41526    | 0.41526    |   0.0 |  0.49
Other   |            | 0.06157    |            |       |  0.07

Nlocal:    4000 ave 4000 max 4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    10895 ave 10895 max 10895 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:  1.50944e+06 ave 1.50944e+06 max 1.50944e+06 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 1509438
Ave neighs/atom = 377.36
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) = 9.436 | 9.436 | 9.436 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    1000   -6540.2745   -6540.2745   -6681.9971   -6681.9971    274.17198    274.17198    230227.51    230227.51   -277.18417   -277.18417 
    2000   -6546.6903   -6546.6903   -6692.1106   -6692.1106    281.32545    281.32545    229403.68    229403.68    22.478051    22.478051 
Loop time of 84.1665 on 1 procs for 1000 steps with 4000 atoms

Performance: 1.027 ns/day, 23.380 hours/ns, 11.881 timesteps/s
39.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    | 83.047     | 83.047     | 83.047     |   0.0 | 98.67
Neigh   | 0.55433    | 0.55433    | 0.55433    |   0.0 |  0.66
Comm    | 0.16816    | 0.16816    | 0.16816    |   0.0 |  0.20
Output  | 4.2915e-05 | 4.2915e-05 | 4.2915e-05 |   0.0 |  0.00
Modify  | 0.29556    | 0.29556    | 0.29556    |   0.0 |  0.35
Other   |            | 0.1015     |            |       |  0.12

Nlocal:    4000 ave 4000 max 4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    10895 ave 10895 max 10895 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:  1.5124e+06 ave 1.5124e+06 max 1.5124e+06 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 1512400
Ave neighs/atom = 378.1
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) = 9.436 | 9.436 | 9.436 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    2000   -6546.6903   -6546.6903   -6692.1106   -6692.1106    281.32545    281.32545    229403.68    229403.68    22.478051    22.478051 
    3000   -6553.7568   -6553.7568   -6698.3244   -6698.3244    279.67571    279.67571    229445.82    229445.82   -48.493505   -48.493505 
Loop time of 83.6556 on 1 procs for 1000 steps with 4000 atoms

Performance: 1.033 ns/day, 23.238 hours/ns, 11.954 timesteps/s
39.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    | 82.453     | 82.453     | 82.453     |   0.0 | 98.56
Neigh   | 0.47499    | 0.47499    | 0.47499    |   0.0 |  0.57
Comm    | 0.19833    | 0.19833    | 0.19833    |   0.0 |  0.24
Output  | 4.2915e-05 | 4.2915e-05 | 4.2915e-05 |   0.0 |  0.00
Modify  | 0.39708    | 0.39708    | 0.39708    |   0.0 |  0.47
Other   |            | 0.1318     |            |       |  0.16

Nlocal:    4000 ave 4000 max 4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    10895 ave 10895 max 10895 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:  1.51179e+06 ave 1.51179e+06 max 1.51179e+06 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 1511788
Ave neighs/atom = 377.947
Neighbor list builds = 5
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) = 9.436 | 9.436 | 9.436 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    3000   -6553.7568   -6553.7568   -6698.3244   -6698.3244    279.67571    279.67571    229445.82    229445.82   -48.493505   -48.493505 
    4000   -6560.3339   -6560.3339   -6696.4534   -6696.4534     263.3323     263.3323    229184.89    229184.89    44.326615    44.326615 
Loop time of 83.9714 on 1 procs for 1000 steps with 4000 atoms

Performance: 1.029 ns/day, 23.325 hours/ns, 11.909 timesteps/s
39.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    | 82.951     | 82.951     | 82.951     |   0.0 | 98.78
Neigh   | 0.45668    | 0.45668    | 0.45668    |   0.0 |  0.54
Comm    | 0.14891    | 0.14891    | 0.14891    |   0.0 |  0.18
Output  | 3.5048e-05 | 3.5048e-05 | 3.5048e-05 |   0.0 |  0.00
Modify  | 0.35327    | 0.35327    | 0.35327    |   0.0 |  0.42
Other   |            | 0.06163    |            |       |  0.07

Nlocal:    4000 ave 4000 max 4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    10895 ave 10895 max 10895 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:  1.51107e+06 ave 1.51107e+06 max 1.51107e+06 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 1511068
Ave neighs/atom = 377.767
Neighbor list builds = 5
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) = 9.436 | 9.436 | 9.436 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    4000   -6560.3339   -6560.3339   -6696.4534   -6696.4534     263.3323     263.3323    229184.89    229184.89    44.326615    44.326615 
    5000   -6546.8247   -6546.8247   -6689.8888   -6689.8888    276.76714    276.76714    229305.82    229305.82    100.23471    100.23471 
Loop time of 83.6476 on 1 procs for 1000 steps with 4000 atoms

Performance: 1.033 ns/day, 23.235 hours/ns, 11.955 timesteps/s
40.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    | 82.657     | 82.657     | 82.657     |   0.0 | 98.82
Neigh   | 0.40375    | 0.40375    | 0.40375    |   0.0 |  0.48
Comm    | 0.14809    | 0.14809    | 0.14809    |   0.0 |  0.18
Output  | 4.1962e-05 | 4.1962e-05 | 4.1962e-05 |   0.0 |  0.00
Modify  | 0.37662    | 0.37662    | 0.37662    |   0.0 |  0.45
Other   |            | 0.06169    |            |       |  0.07

Nlocal:    4000 ave 4000 max 4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    10895 ave 10895 max 10895 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:  1.50979e+06 ave 1.50979e+06 max 1.50979e+06 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 1509790
Ave neighs/atom = 377.447
Neighbor list builds = 5
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 = 280.804481082746, Press = -104.44547118837
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 9.436 | 9.436 | 9.436 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    5000   -6546.8247   -6546.8247   -6689.8888   -6689.8888    276.76714    276.76714    229305.82    229305.82    100.23471    100.23471 
    6000   -6558.5103   -6558.5103   -6694.4752   -6694.4752     263.0333     263.0333    229776.23    229776.23   -222.91983   -222.91983 
Loop time of 83.0781 on 1 procs for 1000 steps with 4000 atoms

Performance: 1.040 ns/day, 23.077 hours/ns, 12.037 timesteps/s
40.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    | 81.888     | 81.888     | 81.888     |   0.0 | 98.57
Neigh   | 0.45973    | 0.45973    | 0.45973    |   0.0 |  0.55
Comm    | 0.20884    | 0.20884    | 0.20884    |   0.0 |  0.25
Output  | 2.5988e-05 | 2.5988e-05 | 2.5988e-05 |   0.0 |  0.00
Modify  | 0.50059    | 0.50059    | 0.50059    |   0.0 |  0.60
Other   |            | 0.02132    |            |       |  0.03

Nlocal:    4000 ave 4000 max 4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    10895 ave 10895 max 10895 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:  1.51166e+06 ave 1.51166e+06 max 1.51166e+06 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 1511662
Ave neighs/atom = 377.916
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 = 273.111412943677, Press = 2.33534720079382
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 9.436 | 9.436 | 9.436 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    6000   -6558.5103   -6558.5103   -6694.4752   -6694.4752     263.0333     263.0333    229776.23    229776.23   -222.91983   -222.91983 
    7000   -6544.9427   -6544.9427   -6692.5024   -6692.5024    285.46426    285.46426    229064.58    229064.58    206.75161    206.75161 
Loop time of 81.8746 on 1 procs for 1000 steps with 4000 atoms

Performance: 1.055 ns/day, 22.743 hours/ns, 12.214 timesteps/s
40.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    | 80.903     | 80.903     | 80.903     |   0.0 | 98.81
Neigh   | 0.38451    | 0.38451    | 0.38451    |   0.0 |  0.47
Comm    | 0.088648   | 0.088648   | 0.088648   |   0.0 |  0.11
Output  | 3.1948e-05 | 3.1948e-05 | 3.1948e-05 |   0.0 |  0.00
Modify  | 0.4558     | 0.4558     | 0.4558     |   0.0 |  0.56
Other   |            | 0.04221    |            |       |  0.05

Nlocal:    4000 ave 4000 max 4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    10895 ave 10895 max 10895 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:  1.51144e+06 ave 1.51144e+06 max 1.51144e+06 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 1511438
Ave neighs/atom = 377.86
Neighbor list builds = 5
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.206737785602, Press = 0.232942688858635
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 9.436 | 9.436 | 9.436 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    7000   -6544.9427   -6544.9427   -6692.5024   -6692.5024    285.46426    285.46426    229064.58    229064.58    206.75161    206.75161 
    8000   -6554.3566   -6554.3566   -6694.9954   -6694.9954    272.07524    272.07524    229600.43    229600.43   -125.75079   -125.75079 
Loop time of 80.5281 on 1 procs for 1000 steps with 4000 atoms

Performance: 1.073 ns/day, 22.369 hours/ns, 12.418 timesteps/s
41.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    | 79.541     | 79.541     | 79.541     |   0.0 | 98.77
Neigh   | 0.4004     | 0.4004     | 0.4004     |   0.0 |  0.50
Comm    | 0.14836    | 0.14836    | 0.14836    |   0.0 |  0.18
Output  | 3.2902e-05 | 3.2902e-05 | 3.2902e-05 |   0.0 |  0.00
Modify  | 0.41649    | 0.41649    | 0.41649    |   0.0 |  0.52
Other   |            | 0.02159    |            |       |  0.03

Nlocal:    4000 ave 4000 max 4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    10895 ave 10895 max 10895 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:  1.51298e+06 ave 1.51298e+06 max 1.51298e+06 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 1512976
Ave neighs/atom = 378.244
Neighbor list builds = 5
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.338487039528, Press = -1.68666681766131
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 9.436 | 9.436 | 9.436 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    8000   -6554.3566   -6554.3566   -6694.9954   -6694.9954    272.07524    272.07524    229600.43    229600.43   -125.75079   -125.75079 
    9000   -6546.5439   -6546.5439   -6689.7023   -6689.7023     276.9495     276.9495    229476.21    229476.21    11.407486    11.407486 
Loop time of 77.1281 on 1 procs for 1000 steps with 4000 atoms

Performance: 1.120 ns/day, 21.424 hours/ns, 12.965 timesteps/s
43.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    | 76.241     | 76.241     | 76.241     |   0.0 | 98.85
Neigh   | 0.38865    | 0.38865    | 0.38865    |   0.0 |  0.50
Comm    | 0.17643    | 0.17643    | 0.17643    |   0.0 |  0.23
Output  | 2.6941e-05 | 2.6941e-05 | 2.6941e-05 |   0.0 |  0.00
Modify  | 0.28003    | 0.28003    | 0.28003    |   0.0 |  0.36
Other   |            | 0.04186    |            |       |  0.05

Nlocal:    4000 ave 4000 max 4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    10895 ave 10895 max 10895 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:  1.51179e+06 ave 1.51179e+06 max 1.51179e+06 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 1511790
Ave neighs/atom = 377.947
Neighbor list builds = 5
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.575585917423, Press = 0.531261578479066
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 9.436 | 9.436 | 9.436 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    9000   -6546.5439   -6546.5439   -6689.7023   -6689.7023     276.9495     276.9495    229476.21    229476.21    11.407486    11.407486 
   10000   -6552.5635   -6552.5635   -6693.3736   -6693.3736    272.40672    272.40672    229208.31    229208.31    89.649688    89.649688 
Loop time of 76.8596 on 1 procs for 1000 steps with 4000 atoms

Performance: 1.124 ns/day, 21.350 hours/ns, 13.011 timesteps/s
43.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    | 76.027     | 76.027     | 76.027     |   0.0 | 98.92
Neigh   | 0.38259    | 0.38259    | 0.38259    |   0.0 |  0.50
Comm    | 0.1083     | 0.1083     | 0.1083     |   0.0 |  0.14
Output  | 3.2902e-05 | 3.2902e-05 | 3.2902e-05 |   0.0 |  0.00
Modify  | 0.30015    | 0.30015    | 0.30015    |   0.0 |  0.39
Other   |            | 0.04162    |            |       |  0.05

Nlocal:    4000 ave 4000 max 4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    10895 ave 10895 max 10895 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:  1.51038e+06 ave 1.51038e+06 max 1.51038e+06 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 1510380
Ave neighs/atom = 377.595
Neighbor list builds = 5
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.575524066205, Press = -0.756524949946992
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 9.436 | 9.436 | 9.436 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   10000   -6552.5635   -6552.5635   -6693.3736   -6693.3736    272.40672    272.40672    229208.31    229208.31    89.649688    89.649688 
   11000   -6548.3434   -6548.3434   -6693.3797   -6693.3797     280.5825     280.5825    229511.77    229511.77   -34.094202   -34.094202 
Loop time of 84.9545 on 1 procs for 1000 steps with 4000 atoms

Performance: 1.017 ns/day, 23.598 hours/ns, 11.771 timesteps/s
39.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    | 83.773     | 83.773     | 83.773     |   0.0 | 98.61
Neigh   | 0.46895    | 0.46895    | 0.46895    |   0.0 |  0.55
Comm    | 0.12894    | 0.12894    | 0.12894    |   0.0 |  0.15
Output  | 3.2902e-05 | 3.2902e-05 | 3.2902e-05 |   0.0 |  0.00
Modify  | 0.52216    | 0.52216    | 0.52216    |   0.0 |  0.61
Other   |            | 0.06168    |            |       |  0.07

Nlocal:    4000 ave 4000 max 4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    10895 ave 10895 max 10895 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:  1.51196e+06 ave 1.51196e+06 max 1.51196e+06 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 1511964
Ave neighs/atom = 377.991
Neighbor list builds = 5
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.497399025841, Press = -0.189111682110184
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 9.436 | 9.436 | 9.436 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   11000   -6548.3434   -6548.3434   -6693.3797   -6693.3797     280.5825     280.5825    229511.77    229511.77   -34.094202   -34.094202 
   12000   -6547.4777   -6547.4777    -6689.795    -6689.795    275.32246    275.32246    229489.48    229489.48   -1.7791188   -1.7791188 
Loop time of 86.4524 on 1 procs for 1000 steps with 4000 atoms

Performance: 0.999 ns/day, 24.015 hours/ns, 11.567 timesteps/s
38.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    | 85.288     | 85.288     | 85.288     |   0.0 | 98.65
Neigh   | 0.51116    | 0.51116    | 0.51116    |   0.0 |  0.59
Comm    | 0.22931    | 0.22931    | 0.22931    |   0.0 |  0.27
Output  | 3.1948e-05 | 3.1948e-05 | 3.1948e-05 |   0.0 |  0.00
Modify  | 0.38196    | 0.38196    | 0.38196    |   0.0 |  0.44
Other   |            | 0.0416     |            |       |  0.05

Nlocal:    4000 ave 4000 max 4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    10895 ave 10895 max 10895 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:  1.51178e+06 ave 1.51178e+06 max 1.51178e+06 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 1511780
Ave neighs/atom = 377.945
Neighbor list builds = 5
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.453493961567, Press = -0.0267982230289081
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 9.436 | 9.436 | 9.436 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   12000   -6547.4777   -6547.4777    -6689.795    -6689.795    275.32246    275.32246    229489.48    229489.48   -1.7791188   -1.7791188 
   13000   -6552.8236   -6552.8236   -6691.8834   -6691.8834    269.02071    269.02071    229392.41    229392.41    9.0215615    9.0215615 
Loop time of 84.7076 on 1 procs for 1000 steps with 4000 atoms

Performance: 1.020 ns/day, 23.530 hours/ns, 11.805 timesteps/s
39.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    | 83.62      | 83.62      | 83.62      |   0.0 | 98.72
Neigh   | 0.45291    | 0.45291    | 0.45291    |   0.0 |  0.53
Comm    | 0.20931    | 0.20931    | 0.20931    |   0.0 |  0.25
Output  | 2.5034e-05 | 2.5034e-05 | 2.5034e-05 |   0.0 |  0.00
Modify  | 0.36284    | 0.36284    | 0.36284    |   0.0 |  0.43
Other   |            | 0.06204    |            |       |  0.07

Nlocal:    4000 ave 4000 max 4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    10895 ave 10895 max 10895 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:  1.51155e+06 ave 1.51155e+06 max 1.51155e+06 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 1511554
Ave neighs/atom = 377.889
Neighbor list builds = 5
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.303163284339, Press = -0.305957429292585
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 9.436 | 9.436 | 9.436 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   13000   -6552.8236   -6552.8236   -6691.8834   -6691.8834    269.02071    269.02071    229392.41    229392.41    9.0215615    9.0215615 
   14000    -6554.265    -6554.265   -6695.4606   -6695.4606    273.15253    273.15253     229371.1     229371.1   -16.818621   -16.818621 
Loop time of 83.8551 on 1 procs for 1000 steps with 4000 atoms

Performance: 1.030 ns/day, 23.293 hours/ns, 11.925 timesteps/s
40.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    | 82.797     | 82.797     | 82.797     |   0.0 | 98.74
Neigh   | 0.38804    | 0.38804    | 0.38804    |   0.0 |  0.46
Comm    | 0.16895    | 0.16895    | 0.16895    |   0.0 |  0.20
Output  | 2.5988e-05 | 2.5988e-05 | 2.5988e-05 |   0.0 |  0.00
Modify  | 0.47908    | 0.47908    | 0.47908    |   0.0 |  0.57
Other   |            | 0.02189    |            |       |  0.03

Nlocal:    4000 ave 4000 max 4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    10895 ave 10895 max 10895 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:  1.51157e+06 ave 1.51157e+06 max 1.51157e+06 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 1511574
Ave neighs/atom = 377.894
Neighbor list builds = 5
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.183303519187, Press = -0.360570498885071
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 9.436 | 9.436 | 9.436 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   14000    -6554.265    -6554.265   -6695.4606   -6695.4606    273.15253    273.15253     229371.1     229371.1   -16.818621   -16.818621 
   15000   -6553.0163   -6553.0163   -6691.4024   -6691.4024    267.71725    267.71725    229384.14    229384.14    7.7253587    7.7253587 
Loop time of 82.0334 on 1 procs for 1000 steps with 4000 atoms

Performance: 1.053 ns/day, 22.787 hours/ns, 12.190 timesteps/s
40.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    | 80.853     | 80.853     | 80.853     |   0.0 | 98.56
Neigh   | 0.4068     | 0.4068     | 0.4068     |   0.0 |  0.50
Comm    | 0.1773     | 0.1773     | 0.1773     |   0.0 |  0.22
Output  | 3.3855e-05 | 3.3855e-05 | 3.3855e-05 |   0.0 |  0.00
Modify  | 0.4943     | 0.4943     | 0.4943     |   0.0 |  0.60
Other   |            | 0.1018     |            |       |  0.12

Nlocal:    4000 ave 4000 max 4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    10895 ave 10895 max 10895 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:  1.51198e+06 ave 1.51198e+06 max 1.51198e+06 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 1511978
Ave neighs/atom = 377.995
Neighbor list builds = 5
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
229407.280760426

print "LAMMPS calculation completed"
LAMMPS calculation completed
quit 0