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

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

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

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

print "Initial system volume: ${V0_metal} Angstroms^3"
Initial system volume: 66983.4414171395 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 313.15*${_u_temperature}
variable temp_converted equal 313.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 313.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 313.15 ${temp_converted} ${Tdamp_converted}               iso ${press_converted} ${press_converted} ${Pdamp_converted}
fix           ensemble all npt temp 313.15 313.15 ${Tdamp_converted}               iso ${press_converted} ${press_converted} ${Pdamp_converted}
fix           ensemble all npt temp 313.15 313.15 0.1               iso ${press_converted} ${press_converted} ${Pdamp_converted}
fix           ensemble all npt temp 313.15 313.15 0.1               iso 0 ${press_converted} ${Pdamp_converted}
fix           ensemble all npt temp 313.15 313.15 0.1               iso 0 0 ${Pdamp_converted}
fix           ensemble all npt temp 313.15 313.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  "313.15 - 0.2"
variable      T_up  equal  "313.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 = 9.80032
  ghost atom cutoff = 9.80032
  binsize = 4.90016, bins = 9 9 9
  1 neighbor lists, perpetual/occasional/extra = 1 0 0
  (1) pair kim, perpetual
      attributes: full, newton off, cut 9.80032
      pair build: full/bin/atomonly
      stencil: full/bin/3d
      bin: standard
Per MPI rank memory allocation (min/avg/max) = 7.004 | 7.004 | 7.004 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
       0   -10639.793   -10639.793   -10801.664   -10801.664       313.15       313.15    66983.441    66983.441    2581.1687    2581.1687 
    1000   -10462.554   -10462.554   -10628.593   -10628.593    321.21365    321.21365    68226.559    68226.559   -1686.3801   -1686.3801 
Loop time of 45.3914 on 1 procs for 1000 steps with 4000 atoms

Performance: 1.903 ns/day, 12.609 hours/ns, 22.031 timesteps/s
37.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    | 44.774     | 44.774     | 44.774     |   0.0 | 98.64
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.15187    | 0.15187    | 0.15187    |   0.0 |  0.33
Output  | 4.1008e-05 | 4.1008e-05 | 4.1008e-05 |   0.0 |  0.00
Modify  | 0.41502    | 0.41502    | 0.41502    |   0.0 |  0.91
Other   |            | 0.05085    |            |       |  0.11

Nlocal:    4000 ave 4000 max 4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    8195 ave 8195 max 8195 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:  896000 ave 896000 max 896000 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 896000
Ave neighs/atom = 224
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) = 7.004 | 7.004 | 7.004 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    1000   -10462.554   -10462.554   -10628.593   -10628.593    321.21365    321.21365    68226.559    68226.559   -1686.3801   -1686.3801 
    2000   -10476.754   -10476.754   -10631.848   -10631.848    300.04024    300.04024    67948.281    67948.281    722.49402    722.49402 
Loop time of 44.3166 on 1 procs for 1000 steps with 4000 atoms

Performance: 1.950 ns/day, 12.310 hours/ns, 22.565 timesteps/s
38.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    | 43.841     | 43.841     | 43.841     |   0.0 | 98.93
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.1113     | 0.1113     | 0.1113     |   0.0 |  0.25
Output  | 4.2915e-05 | 4.2915e-05 | 4.2915e-05 |   0.0 |  0.00
Modify  | 0.34363    | 0.34363    | 0.34363    |   0.0 |  0.78
Other   |            | 0.02053    |            |       |  0.05

Nlocal:    4000 ave 4000 max 4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    8195 ave 8195 max 8195 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:  896418 ave 896418 max 896418 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 896418
Ave neighs/atom = 224.105
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) = 7.004 | 7.004 | 7.004 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    2000   -10476.754   -10476.754   -10631.848   -10631.848    300.04024    300.04024    67948.281    67948.281    722.49402    722.49402 
    3000   -10470.784   -10470.784   -10633.338   -10633.338    314.47148    314.47148    68040.739    68040.739    -191.9935    -191.9935 
Loop time of 44.6212 on 1 procs for 1000 steps with 4000 atoms

Performance: 1.936 ns/day, 12.395 hours/ns, 22.411 timesteps/s
38.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    | 44.184     | 44.184     | 44.184     |   0.0 | 99.02
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.071578   | 0.071578   | 0.071578   |   0.0 |  0.16
Output  | 4.1008e-05 | 4.1008e-05 | 4.1008e-05 |   0.0 |  0.00
Modify  | 0.3445     | 0.3445     | 0.3445     |   0.0 |  0.77
Other   |            | 0.02064    |            |       |  0.05

Nlocal:    4000 ave 4000 max 4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    8195 ave 8195 max 8195 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:  898630 ave 898630 max 898630 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 898630
Ave neighs/atom = 224.657
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) = 7.004 | 7.004 | 7.004 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    3000   -10470.784   -10470.784   -10633.338   -10633.338    314.47148    314.47148    68040.739    68040.739    -191.9935    -191.9935 
    4000     -10470.4     -10470.4    -10634.99    -10634.99    318.40974    318.40974    68002.836    68002.836    121.83817    121.83817 
Loop time of 41.871 on 1 procs for 1000 steps with 4000 atoms

Performance: 2.063 ns/day, 11.631 hours/ns, 23.883 timesteps/s
40.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    | 41.37      | 41.37      | 41.37      |   0.0 | 98.80
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.091863   | 0.091863   | 0.091863   |   0.0 |  0.22
Output  | 4.2915e-05 | 4.2915e-05 | 4.2915e-05 |   0.0 |  0.00
Modify  | 0.33613    | 0.33613    | 0.33613    |   0.0 |  0.80
Other   |            | 0.07339    |            |       |  0.18

Nlocal:    4000 ave 4000 max 4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    8195 ave 8195 max 8195 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:  897892 ave 897892 max 897892 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 897892
Ave neighs/atom = 224.473
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) = 7.004 | 7.004 | 7.004 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    4000     -10470.4     -10470.4    -10634.99    -10634.99    318.40974    318.40974    68002.836    68002.836    121.83817    121.83817 
    5000   -10474.771   -10474.771   -10635.591   -10635.591    311.11718    311.11718    68056.474    68056.474   -519.79716   -519.79716 
Loop time of 42.8225 on 1 procs for 1000 steps with 4000 atoms

Performance: 2.018 ns/day, 11.895 hours/ns, 23.352 timesteps/s
39.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    | 42.235     | 42.235     | 42.235     |   0.0 | 98.63
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.1517     | 0.1517     | 0.1517     |   0.0 |  0.35
Output  | 3.0994e-05 | 3.0994e-05 | 3.0994e-05 |   0.0 |  0.00
Modify  | 0.41518    | 0.41518    | 0.41518    |   0.0 |  0.97
Other   |            | 0.02075    |            |       |  0.05

Nlocal:    4000 ave 4000 max 4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    8195 ave 8195 max 8195 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:  898108 ave 898108 max 898108 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 898108
Ave neighs/atom = 224.527
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 = 313.465886021245, Press = -101.469101779806
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 7.004 | 7.004 | 7.004 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    5000   -10474.771   -10474.771   -10635.591   -10635.591    311.11718    311.11718    68056.474    68056.474   -519.79716   -519.79716 
    6000    -10468.91    -10468.91   -10631.282   -10631.282    314.11832    314.11832    67986.169    67986.169    494.22352    494.22352 
Loop time of 41.6529 on 1 procs for 1000 steps with 4000 atoms

Performance: 2.074 ns/day, 11.570 hours/ns, 24.008 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    | 41.092     | 41.092     | 41.092     |   0.0 | 98.65
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.1819     | 0.1819     | 0.1819     |   0.0 |  0.44
Output  | 2.5988e-05 | 2.5988e-05 | 2.5988e-05 |   0.0 |  0.00
Modify  | 0.32881    | 0.32881    | 0.32881    |   0.0 |  0.79
Other   |            | 0.05069    |            |       |  0.12

Nlocal:    4000 ave 4000 max 4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    8195 ave 8195 max 8195 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:  897656 ave 897656 max 897656 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 897656
Ave neighs/atom = 224.414
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 = 312.790775625203, Press = -14.7075881408817
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 7.004 | 7.004 | 7.004 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    6000    -10468.91    -10468.91   -10631.282   -10631.282    314.11832    314.11832    67986.169    67986.169    494.22352    494.22352 
    7000   -10475.032   -10475.032   -10635.039   -10635.039     309.5439     309.5439    68032.657    68032.657   -237.63358   -237.63358 
Loop time of 49.5598 on 1 procs for 1000 steps with 4000 atoms

Performance: 1.743 ns/day, 13.767 hours/ns, 20.178 timesteps/s
34.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    | 48.905     | 48.905     | 48.905     |   0.0 | 98.68
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.16447    | 0.16447    | 0.16447    |   0.0 |  0.33
Output  | 3.1948e-05 | 3.1948e-05 | 3.1948e-05 |   0.0 |  0.00
Modify  | 0.4091     | 0.4091     | 0.4091     |   0.0 |  0.83
Other   |            | 0.08078    |            |       |  0.16

Nlocal:    4000 ave 4000 max 4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    8195 ave 8195 max 8195 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:  898096 ave 898096 max 898096 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 898096
Ave neighs/atom = 224.524
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 = 313.054687358545, Press = 7.26619172245757
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 7.004 | 7.004 | 7.004 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    7000   -10475.032   -10475.032   -10635.039   -10635.039     309.5439     309.5439    68032.657    68032.657   -237.63358   -237.63358 
    8000   -10471.396   -10471.396   -10633.922   -10633.922    314.41776    314.41776    68024.812    68024.812   -64.491446   -64.491446 
Loop time of 51.2964 on 1 procs for 1000 steps with 4000 atoms

Performance: 1.684 ns/day, 14.249 hours/ns, 19.495 timesteps/s
33.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    | 50.723     | 50.723     | 50.723     |   0.0 | 98.88
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.13236    | 0.13236    | 0.13236    |   0.0 |  0.26
Output  | 5.0068e-05 | 5.0068e-05 | 5.0068e-05 |   0.0 |  0.00
Modify  | 0.40026    | 0.40026    | 0.40026    |   0.0 |  0.78
Other   |            | 0.041      |            |       |  0.08

Nlocal:    4000 ave 4000 max 4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    8195 ave 8195 max 8195 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:  897802 ave 897802 max 897802 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 897802
Ave neighs/atom = 224.451
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 = 313.099958236645, Press = -20.0046569116553
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 7.004 | 7.004 | 7.004 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    8000   -10471.396   -10471.396   -10633.922   -10633.922    314.41776    314.41776    68024.812    68024.812   -64.491446   -64.491446 
    9000   -10476.301   -10476.301   -10637.953   -10637.953    312.72736    312.72736    67951.036    67951.036    387.06924    387.06924 
Loop time of 51.4481 on 1 procs for 1000 steps with 4000 atoms

Performance: 1.679 ns/day, 14.291 hours/ns, 19.437 timesteps/s
33.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    | 50.7       | 50.7       | 50.7       |   0.0 | 98.55
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.14209    | 0.14209    | 0.14209    |   0.0 |  0.28
Output  | 3.0994e-05 | 3.0994e-05 | 3.0994e-05 |   0.0 |  0.00
Modify  | 0.58521    | 0.58521    | 0.58521    |   0.0 |  1.14
Other   |            | 0.02084    |            |       |  0.04

Nlocal:    4000 ave 4000 max 4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    8195 ave 8195 max 8195 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:  898008 ave 898008 max 898008 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 898008
Ave neighs/atom = 224.502
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 = 312.697086165299, Press = 8.866125093486
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 7.004 | 7.004 | 7.004 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    9000   -10476.301   -10476.301   -10637.953   -10637.953    312.72736    312.72736    67951.036    67951.036    387.06924    387.06924 
   10000   -10469.173   -10469.173   -10630.413   -10630.413    311.92948    311.92948    68061.846    68061.846   -271.02096   -271.02096 
Loop time of 51.3858 on 1 procs for 1000 steps with 4000 atoms

Performance: 1.681 ns/day, 14.274 hours/ns, 19.461 timesteps/s
33.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    | 50.576     | 50.576     | 50.576     |   0.0 | 98.42
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.23218    | 0.23218    | 0.23218    |   0.0 |  0.45
Output  | 2.5988e-05 | 2.5988e-05 | 2.5988e-05 |   0.0 |  0.00
Modify  | 0.51644    | 0.51644    | 0.51644    |   0.0 |  1.01
Other   |            | 0.06095    |            |       |  0.12

Nlocal:    4000 ave 4000 max 4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    8195 ave 8195 max 8195 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:  898332 ave 898332 max 898332 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 898332
Ave neighs/atom = 224.583
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 = 312.691889794815, Press = -4.40486716651125
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 7.004 | 7.004 | 7.004 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   10000   -10469.173   -10469.173   -10630.413   -10630.413    311.92948    311.92948    68061.846    68061.846   -271.02096   -271.02096 
   11000   -10474.761   -10474.761   -10634.761   -10634.761    309.52924    309.52924    67984.022    67984.022     236.7747     236.7747 
Loop time of 51.3174 on 1 procs for 1000 steps with 4000 atoms

Performance: 1.684 ns/day, 14.255 hours/ns, 19.487 timesteps/s
33.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    | 50.725     | 50.725     | 50.725     |   0.0 | 98.85
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.11187    | 0.11187    | 0.11187    |   0.0 |  0.22
Output  | 3.3855e-05 | 3.3855e-05 | 3.3855e-05 |   0.0 |  0.00
Modify  | 0.35919    | 0.35919    | 0.35919    |   0.0 |  0.70
Other   |            | 0.1211     |            |       |  0.24

Nlocal:    4000 ave 4000 max 4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    8195 ave 8195 max 8195 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:  897638 ave 897638 max 897638 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 897638
Ave neighs/atom = 224.41
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 = 312.781542625052, Press = -2.23068433210558
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 7.004 | 7.004 | 7.004 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   11000   -10474.761   -10474.761   -10634.761   -10634.761    309.52924    309.52924    67984.022    67984.022     236.7747     236.7747 
   12000   -10470.868   -10470.868   -10631.337   -10631.337    310.43813    310.43813    68032.817    68032.817   -34.901893   -34.901893 
Loop time of 51.4532 on 1 procs for 1000 steps with 4000 atoms

Performance: 1.679 ns/day, 14.293 hours/ns, 19.435 timesteps/s
33.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    | 50.81      | 50.81      | 50.81      |   0.0 | 98.75
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.10176    | 0.10176    | 0.10176    |   0.0 |  0.20
Output  | 2.6941e-05 | 2.6941e-05 | 2.6941e-05 |   0.0 |  0.00
Modify  | 0.46038    | 0.46038    | 0.46038    |   0.0 |  0.89
Other   |            | 0.08087    |            |       |  0.16

Nlocal:    4000 ave 4000 max 4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    8195 ave 8195 max 8195 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:  898262 ave 898262 max 898262 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 898262
Ave neighs/atom = 224.565
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 = 312.595437073053, Press = -0.295134691036733
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 7.004 | 7.004 | 7.004 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   12000   -10470.868   -10470.868   -10631.337   -10631.337    310.43813    310.43813    68032.817    68032.817   -34.901893   -34.901893 
   13000   -10475.517   -10475.517   -10636.025   -10636.025    310.51392    310.51392    67976.789    67976.789    251.90072    251.90072 
Loop time of 49.3428 on 1 procs for 1000 steps with 4000 atoms

Performance: 1.751 ns/day, 13.706 hours/ns, 20.266 timesteps/s
34.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    | 48.694     | 48.694     | 48.694     |   0.0 | 98.68
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.20856    | 0.20856    | 0.20856    |   0.0 |  0.42
Output  | 2.5988e-05 | 2.5988e-05 | 2.5988e-05 |   0.0 |  0.00
Modify  | 0.37936    | 0.37936    | 0.37936    |   0.0 |  0.77
Other   |            | 0.06091    |            |       |  0.12

Nlocal:    4000 ave 4000 max 4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    8195 ave 8195 max 8195 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:  897856 ave 897856 max 897856 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 897856
Ave neighs/atom = 224.464
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 = 312.486107118472, Press = -3.47548289055458
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 7.004 | 7.004 | 7.004 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   13000   -10475.517   -10475.517   -10636.025   -10636.025    310.51392    310.51392    67976.789    67976.789    251.90072    251.90072 
   14000   -10469.118   -10469.118   -10630.867   -10630.867    312.91499    312.91499    68032.023    68032.023    36.809482    36.809482 
Loop time of 49.4778 on 1 procs for 1000 steps with 4000 atoms

Performance: 1.746 ns/day, 13.744 hours/ns, 20.211 timesteps/s
34.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    | 48.765     | 48.765     | 48.765     |   0.0 | 98.56
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.27226    | 0.27226    | 0.27226    |   0.0 |  0.55
Output  | 3.1948e-05 | 3.1948e-05 | 3.1948e-05 |   0.0 |  0.00
Modify  | 0.39952    | 0.39952    | 0.39952    |   0.0 |  0.81
Other   |            | 0.04089    |            |       |  0.08

Nlocal:    4000 ave 4000 max 4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    8195 ave 8195 max 8195 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:  898080 ave 898080 max 898080 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 898080
Ave neighs/atom = 224.52
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 = 312.760292079805, Press = 2.43142771415145
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 7.004 | 7.004 | 7.004 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   14000   -10469.118   -10469.118   -10630.867   -10630.867    312.91499    312.91499    68032.023    68032.023    36.809482    36.809482 
   15000   -10470.336   -10470.336    -10632.78    -10632.78    314.25929    314.25929    68026.351    68026.351   -67.536649   -67.536649 
Loop time of 46.6122 on 1 procs for 1000 steps with 4000 atoms

Performance: 1.854 ns/day, 12.948 hours/ns, 21.454 timesteps/s
36.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    | 45.998     | 45.998     | 45.998     |   0.0 | 98.68
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.095056   | 0.095056   | 0.095056   |   0.0 |  0.20
Output  | 3.2902e-05 | 3.2902e-05 | 3.2902e-05 |   0.0 |  0.00
Modify  | 0.4989     | 0.4989     | 0.4989     |   0.0 |  1.07
Other   |            | 0.02045    |            |       |  0.04

Nlocal:    4000 ave 4000 max 4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    8195 ave 8195 max 8195 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:  898024 ave 898024 max 898024 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 898024
Ave neighs/atom = 224.506
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 = 312.845794273778, Press = -4.20161990284275
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 7.004 | 7.004 | 7.004 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   15000   -10470.336   -10470.336    -10632.78    -10632.78    314.25929    314.25929    68026.351    68026.351   -67.536649   -67.536649 
   16000   -10473.352   -10473.352   -10634.171   -10634.171    311.11518    311.11518    67959.698    67959.698    534.03174    534.03174 
Loop time of 49.6238 on 1 procs for 1000 steps with 4000 atoms

Performance: 1.741 ns/day, 13.784 hours/ns, 20.152 timesteps/s
34.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    | 49.012     | 49.012     | 49.012     |   0.0 | 98.77
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.17188    | 0.17188    | 0.17188    |   0.0 |  0.35
Output  | 2.5988e-05 | 2.5988e-05 | 2.5988e-05 |   0.0 |  0.00
Modify  | 0.37958    | 0.37958    | 0.37958    |   0.0 |  0.76
Other   |            | 0.06073    |            |       |  0.12

Nlocal:    4000 ave 4000 max 4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    8195 ave 8195 max 8195 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:  897946 ave 897946 max 897946 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 897946
Ave neighs/atom = 224.487
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 = 312.919401643558, Press = 2.18016688885896
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 7.004 | 7.004 | 7.004 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   16000   -10473.352   -10473.352   -10634.171   -10634.171    311.11518    311.11518    67959.698    67959.698    534.03174    534.03174 
   17000   -10469.186   -10469.186   -10633.947   -10633.947    318.74064    318.74064    68084.314    68084.314   -644.11693   -644.11693 
Loop time of 49.5384 on 1 procs for 1000 steps with 4000 atoms

Performance: 1.744 ns/day, 13.761 hours/ns, 20.186 timesteps/s
34.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    | 48.796     | 48.796     | 48.796     |   0.0 | 98.50
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.21194    | 0.21194    | 0.21194    |   0.0 |  0.43
Output  | 3.2187e-05 | 3.2187e-05 | 3.2187e-05 |   0.0 |  0.00
Modify  | 0.44976    | 0.44976    | 0.44976    |   0.0 |  0.91
Other   |            | 0.08071    |            |       |  0.16

Nlocal:    4000 ave 4000 max 4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    8195 ave 8195 max 8195 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:  898528 ave 898528 max 898528 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 898528
Ave neighs/atom = 224.632
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 = 312.7724980217, Press = -2.2990892042956
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 7.004 | 7.004 | 7.004 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   17000   -10469.186   -10469.186   -10633.947   -10633.947    318.74064    318.74064    68084.314    68084.314   -644.11693   -644.11693 
   18000   -10471.836   -10471.836   -10635.021   -10635.021    315.69364    315.69364    67954.574    67954.574    589.48712    589.48712 
Loop time of 47.393 on 1 procs for 1000 steps with 4000 atoms

Performance: 1.823 ns/day, 13.165 hours/ns, 21.100 timesteps/s
36.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    | 46.801     | 46.801     | 46.801     |   0.0 | 98.75
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.11175    | 0.11175    | 0.11175    |   0.0 |  0.24
Output  | 3.3855e-05 | 3.3855e-05 | 3.3855e-05 |   0.0 |  0.00
Modify  | 0.43944    | 0.43944    | 0.43944    |   0.0 |  0.93
Other   |            | 0.04082    |            |       |  0.09

Nlocal:    4000 ave 4000 max 4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    8195 ave 8195 max 8195 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:  897624 ave 897624 max 897624 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 897624
Ave neighs/atom = 224.406
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 = 312.847941827615, Press = 1.07450159134199
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 7.004 | 7.004 | 7.004 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   18000   -10471.836   -10471.836   -10635.021   -10635.021    315.69364    315.69364    67954.574    67954.574    589.48712    589.48712 
   19000   -10472.507   -10472.507    -10634.35    -10634.35    313.09629    313.09629     68092.19     68092.19   -785.50976   -785.50976 
Loop time of 46.7106 on 1 procs for 1000 steps with 4000 atoms

Performance: 1.850 ns/day, 12.975 hours/ns, 21.408 timesteps/s
36.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    | 46.051     | 46.051     | 46.051     |   0.0 | 98.59
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.1316     | 0.1316     | 0.1316     |   0.0 |  0.28
Output  | 3.314e-05  | 3.314e-05  | 3.314e-05  |   0.0 |  0.00
Modify  | 0.44691    | 0.44691    | 0.44691    |   0.0 |  0.96
Other   |            | 0.08097    |            |       |  0.17

Nlocal:    4000 ave 4000 max 4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    8195 ave 8195 max 8195 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:  898302 ave 898302 max 898302 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 898302
Ave neighs/atom = 224.576
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 = 312.90842541386, Press = -1.83449346001708
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 7.004 | 7.004 | 7.004 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   19000   -10472.507   -10472.507    -10634.35    -10634.35    313.09629    313.09629     68092.19     68092.19   -785.50976   -785.50976 
   20000   -10463.854   -10463.854    -10631.76    -10631.76     324.8269     324.8269    67959.893    67959.893    736.47676    736.47676 
Loop time of 46.457 on 1 procs for 1000 steps with 4000 atoms

Performance: 1.860 ns/day, 12.905 hours/ns, 21.525 timesteps/s
36.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    | 45.835     | 45.835     | 45.835     |   0.0 | 98.66
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.11128    | 0.11128    | 0.11128    |   0.0 |  0.24
Output  | 3.1948e-05 | 3.1948e-05 | 3.1948e-05 |   0.0 |  0.00
Modify  | 0.43       | 0.43       | 0.43       |   0.0 |  0.93
Other   |            | 0.08073    |            |       |  0.17

Nlocal:    4000 ave 4000 max 4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    8195 ave 8195 max 8195 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:  897562 ave 897562 max 897562 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 897562
Ave neighs/atom = 224.391
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 = 312.966201795809, Press = 0.305820636892832
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 7.004 | 7.004 | 7.004 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   20000   -10463.854   -10463.854    -10631.76    -10631.76     324.8269     324.8269    67959.893    67959.893    736.47676    736.47676 
   21000   -10474.004   -10474.004   -10636.642   -10636.642    314.63334    314.63334     68112.75     68112.75   -1119.4313   -1119.4313 
Loop time of 42.9071 on 1 procs for 1000 steps with 4000 atoms

Performance: 2.014 ns/day, 11.919 hours/ns, 23.306 timesteps/s
39.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    | 42.406     | 42.406     | 42.406     |   0.0 | 98.83
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.11156    | 0.11156    | 0.11156    |   0.0 |  0.26
Output  | 3.1948e-05 | 3.1948e-05 | 3.1948e-05 |   0.0 |  0.00
Modify  | 0.34912    | 0.34912    | 0.34912    |   0.0 |  0.81
Other   |            | 0.04065    |            |       |  0.09

Nlocal:    4000 ave 4000 max 4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    8195 ave 8195 max 8195 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:  898560 ave 898560 max 898560 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 898560
Ave neighs/atom = 224.64
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 = 313.064584818345, Press = -1.09319854245415
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 7.004 | 7.004 | 7.004 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   21000   -10474.004   -10474.004   -10636.642   -10636.642    314.63334    314.63334     68112.75     68112.75   -1119.4313   -1119.4313 
   22000   -10469.723   -10469.723   -10633.461   -10633.461    316.76205    316.76205    67944.192    67944.192    760.69142    760.69142 
Loop time of 40.9055 on 1 procs for 1000 steps with 4000 atoms

Performance: 2.112 ns/day, 11.363 hours/ns, 24.447 timesteps/s
42.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    | 40.28      | 40.28      | 40.28      |   0.0 | 98.47
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.14688    | 0.14688    | 0.14688    |   0.0 |  0.36
Output  | 5.1022e-05 | 5.1022e-05 | 5.1022e-05 |   0.0 |  0.00
Modify  | 0.41878    | 0.41878    | 0.41878    |   0.0 |  1.02
Other   |            | 0.05933    |            |       |  0.15

Nlocal:    4000 ave 4000 max 4000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    8195 ave 8195 max 8195 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:  897654 ave 897654 max 897654 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 897654
Ave neighs/atom = 224.413
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_T313.15.out"  else  "print 'not_converged'  file output/vol_T313.15.out"
print '${V}'  file output/vol_T313.15.out
68017.2216462439

print "LAMMPS calculation completed"
LAMMPS calculation completed
quit 0