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

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

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

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

print "Initial system volume: ${V0_metal} Angstroms^3"
Initial system volume: 32174.6944849008 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 = 8.78363
  ghost atom cutoff = 8.78363
  binsize = 4.39181, bins = 8 8 8
  1 neighbor lists, perpetual/occasional/extra = 1 0 0
  (1) pair kim, perpetual
      attributes: full, newton off, cut 8.78363
      pair build: full/bin/atomonly
      stencil: full/bin/3d
      bin: standard
Per MPI rank memory allocation (min/avg/max) = 4.577 | 4.577 | 4.577 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
       0   -16475.007   -16475.007   -16555.922   -16555.922       313.15       313.15    32174.694    32174.694    2686.1605    2686.1605 
    1000   -16391.539   -16391.539   -16474.747   -16474.747    322.02049    322.02049    32327.856    32327.856   -2654.4841   -2654.4841 
Loop time of 13.061 on 1 procs for 1000 steps with 2000 atoms

Performance: 6.615 ns/day, 3.628 hours/ns, 76.564 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    | 12.773     | 12.773     | 12.773     |   0.0 | 97.79
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.073108   | 0.073108   | 0.073108   |   0.0 |  0.56
Output  | 3.6001e-05 | 3.6001e-05 | 3.6001e-05 |   0.0 |  0.00
Modify  | 0.20313    | 0.20313    | 0.20313    |   0.0 |  1.56
Other   |            | 0.01209    |            |       |  0.09

Nlocal:    2000 ave 2000 max 2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    5471 ave 5471 max 5471 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:  336000 ave 336000 max 336000 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 336000
Ave neighs/atom = 168
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) = 4.577 | 4.577 | 4.577 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    1000   -16391.539   -16391.539   -16474.747   -16474.747    322.02049    322.02049    32327.856    32327.856   -2654.4841   -2654.4841 
    2000   -16391.826   -16391.826   -16469.637   -16469.637    301.13793    301.13793     32263.59     32263.59    4169.8082    4169.8082 
Loop time of 13.9103 on 1 procs for 1000 steps with 2000 atoms

Performance: 6.211 ns/day, 3.864 hours/ns, 71.889 timesteps/s
35.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    | 13.482     | 13.482     | 13.482     |   0.0 | 96.92
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.032953   | 0.032953   | 0.032953   |   0.0 |  0.24
Output  | 3.8147e-05 | 3.8147e-05 | 3.8147e-05 |   0.0 |  0.00
Modify  | 0.30348    | 0.30348    | 0.30348    |   0.0 |  2.18
Other   |            | 0.092      |            |       |  0.66

Nlocal:    2000 ave 2000 max 2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    5471 ave 5471 max 5471 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:  336030 ave 336030 max 336030 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 336030
Ave neighs/atom = 168.015
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) = 4.577 | 4.577 | 4.577 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    2000   -16391.826   -16391.826   -16469.637   -16469.637    301.13793    301.13793     32263.59     32263.59    4169.8082    4169.8082 
    3000   -16394.796   -16394.796   -16474.615   -16474.615    308.90743    308.90743    32321.118    32321.118   -2025.4903   -2025.4903 
Loop time of 13.5209 on 1 procs for 1000 steps with 2000 atoms

Performance: 6.390 ns/day, 3.756 hours/ns, 73.959 timesteps/s
36.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    | 13.092     | 13.092     | 13.092     |   0.0 | 96.83
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.093171   | 0.093171   | 0.093171   |   0.0 |  0.69
Output  | 4.0054e-05 | 4.0054e-05 | 4.0054e-05 |   0.0 |  0.00
Modify  | 0.32326    | 0.32326    | 0.32326    |   0.0 |  2.39
Other   |            | 0.012      |            |       |  0.09

Nlocal:    2000 ave 2000 max 2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    5471 ave 5471 max 5471 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:  336034 ave 336034 max 336034 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 336034
Ave neighs/atom = 168.017
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) = 4.577 | 4.577 | 4.577 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    3000   -16394.796   -16394.796   -16474.615   -16474.615    308.90743    308.90743    32321.118    32321.118   -2025.4903   -2025.4903 
    4000   -16392.841   -16392.841   -16472.597   -16472.597    308.66491    308.66491    32282.398    32282.398    2024.3145    2024.3145 
Loop time of 14.3495 on 1 procs for 1000 steps with 2000 atoms

Performance: 6.021 ns/day, 3.986 hours/ns, 69.689 timesteps/s
35.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    | 14.009     | 14.009     | 14.009     |   0.0 | 97.63
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.053631   | 0.053631   | 0.053631   |   0.0 |  0.37
Output  | 3.8147e-05 | 3.8147e-05 | 3.8147e-05 |   0.0 |  0.00
Modify  | 0.25426    | 0.25426    | 0.25426    |   0.0 |  1.77
Other   |            | 0.03227    |            |       |  0.22

Nlocal:    2000 ave 2000 max 2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    5471 ave 5471 max 5471 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:  336006 ave 336006 max 336006 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 336006
Ave neighs/atom = 168.003
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) = 4.577 | 4.577 | 4.577 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    4000   -16392.841   -16392.841   -16472.597   -16472.597    308.66491    308.66491    32282.398    32282.398    2024.3145    2024.3145 
    5000    -16395.66    -16395.66   -16473.683   -16473.683    301.95663    301.95663    32308.788    32308.788   -811.18834   -811.18834 
Loop time of 13.9624 on 1 procs for 1000 steps with 2000 atoms

Performance: 6.188 ns/day, 3.878 hours/ns, 71.621 timesteps/s
35.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    | 13.633     | 13.633     | 13.633     |   0.0 | 97.64
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.07313    | 0.07313    | 0.07313    |   0.0 |  0.52
Output  | 2.6941e-05 | 2.6941e-05 | 2.6941e-05 |   0.0 |  0.00
Modify  | 0.24424    | 0.24424    | 0.24424    |   0.0 |  1.75
Other   |            | 0.01222    |            |       |  0.09

Nlocal:    2000 ave 2000 max 2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    5471 ave 5471 max 5471 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:  336028 ave 336028 max 336028 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 336028
Ave neighs/atom = 168.014
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 = 306.256135645777, Press = 77.8326110364853
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 4.577 | 4.577 | 4.577 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    5000    -16395.66    -16395.66   -16473.683   -16473.683    301.95663    301.95663    32308.788    32308.788   -811.18834   -811.18834 
    6000   -16392.791   -16392.791   -16474.818   -16474.818    317.45577    317.45577    32291.676    32291.676    868.74472    868.74472 
Loop time of 13.6903 on 1 procs for 1000 steps with 2000 atoms

Performance: 6.311 ns/day, 3.803 hours/ns, 73.044 timesteps/s
36.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    | 13.198     | 13.198     | 13.198     |   0.0 | 96.41
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.093121   | 0.093121   | 0.093121   |   0.0 |  0.68
Output  | 2.8849e-05 | 2.8849e-05 | 2.8849e-05 |   0.0 |  0.00
Modify  | 0.36675    | 0.36675    | 0.36675    |   0.0 |  2.68
Other   |            | 0.0321     |            |       |  0.23

Nlocal:    2000 ave 2000 max 2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    5471 ave 5471 max 5471 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:  336026 ave 336026 max 336026 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 336026
Ave neighs/atom = 168.013
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.587014010914, Press = -29.754133553825
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 4.577 | 4.577 | 4.577 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    6000   -16392.791   -16392.791   -16474.818   -16474.818    317.45577    317.45577    32291.676    32291.676    868.74472    868.74472 
    7000   -16396.233   -16396.233   -16476.799   -16476.799    311.79836    311.79836    32319.955    32319.955   -2173.6107   -2173.6107 
Loop time of 11.2358 on 1 procs for 1000 steps with 2000 atoms

Performance: 7.690 ns/day, 3.121 hours/ns, 89.001 timesteps/s
43.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    | 10.914     | 10.914     | 10.914     |   0.0 | 97.14
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.042992   | 0.042992   | 0.042992   |   0.0 |  0.38
Output  | 2.7895e-05 | 2.7895e-05 | 2.7895e-05 |   0.0 |  0.00
Modify  | 0.24628    | 0.24628    | 0.24628    |   0.0 |  2.19
Other   |            | 0.03217    |            |       |  0.29

Nlocal:    2000 ave 2000 max 2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    5471 ave 5471 max 5471 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:  336026 ave 336026 max 336026 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 336026
Ave neighs/atom = 168.013
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.705698151986, Press = 15.5736858652599
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 4.577 | 4.577 | 4.577 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    7000   -16396.233   -16396.233   -16476.799   -16476.799    311.79836    311.79836    32319.955    32319.955   -2173.6107   -2173.6107 
    8000   -16394.785   -16394.785    -16476.34    -16476.34    315.62557    315.62557    32277.891    32277.891    2042.1069    2042.1069 
Loop time of 13.4934 on 1 procs for 1000 steps with 2000 atoms

Performance: 6.403 ns/day, 3.748 hours/ns, 74.110 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    | 13.103     | 13.103     | 13.103     |   0.0 | 97.10
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.11317    | 0.11317    | 0.11317    |   0.0 |  0.84
Output  | 2.8849e-05 | 2.8849e-05 | 2.8849e-05 |   0.0 |  0.00
Modify  | 0.26547    | 0.26547    | 0.26547    |   0.0 |  1.97
Other   |            | 0.01202    |            |       |  0.09

Nlocal:    2000 ave 2000 max 2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    5471 ave 5471 max 5471 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:  336030 ave 336030 max 336030 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 336030
Ave neighs/atom = 168.015
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.403647751196, Press = -23.1774162815115
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 4.577 | 4.577 | 4.577 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    8000   -16394.785   -16394.785    -16476.34    -16476.34    315.62557    315.62557    32277.891    32277.891    2042.1069    2042.1069 
    9000   -16394.247   -16394.247   -16473.948   -16473.948    308.44918    308.44918    32303.161    32303.161   -166.33007   -166.33007 
Loop time of 13.3721 on 1 procs for 1000 steps with 2000 atoms

Performance: 6.461 ns/day, 3.714 hours/ns, 74.783 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    | 12.99      | 12.99      | 12.99      |   0.0 | 97.14
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.1129     | 0.1129     | 0.1129     |   0.0 |  0.84
Output  | 2.3127e-05 | 2.3127e-05 | 2.3127e-05 |   0.0 |  0.00
Modify  | 0.23705    | 0.23705    | 0.23705    |   0.0 |  1.77
Other   |            | 0.03217    |            |       |  0.24

Nlocal:    2000 ave 2000 max 2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    5471 ave 5471 max 5471 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:  336028 ave 336028 max 336028 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 336028
Ave neighs/atom = 168.014
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.015683615058, Press = 9.0205033482872
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 4.577 | 4.577 | 4.577 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    9000   -16394.247   -16394.247   -16473.948   -16473.948    308.44918    308.44918    32303.161    32303.161   -166.33007   -166.33007 
   10000   -16392.743   -16392.743   -16471.665   -16471.665    305.43461    305.43461    32305.378    32305.378   -273.97104   -273.97104 
Loop time of 12.5842 on 1 procs for 1000 steps with 2000 atoms

Performance: 6.866 ns/day, 3.496 hours/ns, 79.465 timesteps/s
39.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    | 12.313     | 12.313     | 12.313     |   0.0 | 97.85
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.032518   | 0.032518   | 0.032518   |   0.0 |  0.26
Output  | 2.9087e-05 | 2.9087e-05 | 2.9087e-05 |   0.0 |  0.00
Modify  | 0.20626    | 0.20626    | 0.20626    |   0.0 |  1.64
Other   |            | 0.03192    |            |       |  0.25

Nlocal:    2000 ave 2000 max 2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    5471 ave 5471 max 5471 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:  336032 ave 336032 max 336032 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 336032
Ave neighs/atom = 168.016
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.976469805571, Press = -5.48963271881269
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 4.577 | 4.577 | 4.577 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   10000   -16392.743   -16392.743   -16471.665   -16471.665    305.43461    305.43461    32305.378    32305.378   -273.97104   -273.97104 
   11000   -16392.076   -16392.076   -16474.583   -16474.583    319.30865    319.30865    32295.631    32295.631    547.09344    547.09344 
Loop time of 12.5289 on 1 procs for 1000 steps with 2000 atoms

Performance: 6.896 ns/day, 3.480 hours/ns, 79.816 timesteps/s
39.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    | 12.178     | 12.178     | 12.178     |   0.0 | 97.20
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.072513   | 0.072513   | 0.072513   |   0.0 |  0.58
Output  | 2.7895e-05 | 2.7895e-05 | 2.7895e-05 |   0.0 |  0.00
Modify  | 0.26686    | 0.26686    | 0.26686    |   0.0 |  2.13
Other   |            | 0.01192    |            |       |  0.10

Nlocal:    2000 ave 2000 max 2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    5471 ave 5471 max 5471 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:  336046 ave 336046 max 336046 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 336046
Ave neighs/atom = 168.023
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.859409078175, Press = 0.334354571335405
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 4.577 | 4.577 | 4.577 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   11000   -16392.076   -16392.076   -16474.583   -16474.583    319.30865    319.30865    32295.631    32295.631    547.09344    547.09344 
   12000    -16392.99    -16392.99   -16473.374   -16473.374    311.09265    311.09265    32309.247    32309.247   -805.98727   -805.98727 
Loop time of 12.7576 on 1 procs for 1000 steps with 2000 atoms

Performance: 6.772 ns/day, 3.544 hours/ns, 78.384 timesteps/s
38.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    | 12.269     | 12.269     | 12.269     |   0.0 | 96.17
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.17287    | 0.17287    | 0.17287    |   0.0 |  1.36
Output  | 2.2173e-05 | 2.2173e-05 | 2.2173e-05 |   0.0 |  0.00
Modify  | 0.30392    | 0.30392    | 0.30392    |   0.0 |  2.38
Other   |            | 0.01188    |            |       |  0.09

Nlocal:    2000 ave 2000 max 2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    5471 ave 5471 max 5471 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:  336030 ave 336030 max 336030 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 336030
Ave neighs/atom = 168.015
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.861223669778, Press = -2.98225907463284
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 4.577 | 4.577 | 4.577 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   12000    -16392.99    -16392.99   -16473.374   -16473.374    311.09265    311.09265    32309.247    32309.247   -805.98727   -805.98727 
   13000   -16395.581   -16395.581   -16474.826   -16474.826    306.68862    306.68862    32276.259    32276.259     2229.944     2229.944 
Loop time of 13.273 on 1 procs for 1000 steps with 2000 atoms

Performance: 6.509 ns/day, 3.687 hours/ns, 75.341 timesteps/s
37.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    | 13.022     | 13.022     | 13.022     |   0.0 | 98.11
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.052827   | 0.052827   | 0.052827   |   0.0 |  0.40
Output  | 2.8133e-05 | 2.8133e-05 | 2.8133e-05 |   0.0 |  0.00
Modify  | 0.16655    | 0.16655    | 0.16655    |   0.0 |  1.25
Other   |            | 0.0319     |            |       |  0.24

Nlocal:    2000 ave 2000 max 2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    5471 ave 5471 max 5471 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:  336030 ave 336030 max 336030 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 336030
Ave neighs/atom = 168.015
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.687683094222, Press = -0.654833804872145
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 4.577 | 4.577 | 4.577 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   13000   -16395.581   -16395.581   -16474.826   -16474.826    306.68862    306.68862    32276.259    32276.259     2229.944     2229.944 
   14000   -16389.839   -16389.839   -16473.145   -16473.145    322.40149    322.40149    32338.322    32338.322   -3465.5314   -3465.5314 
Loop time of 12.8084 on 1 procs for 1000 steps with 2000 atoms

Performance: 6.746 ns/day, 3.558 hours/ns, 78.074 timesteps/s
38.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    | 12.397     | 12.397     | 12.397     |   0.0 | 96.79
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.13275    | 0.13275    | 0.13275    |   0.0 |  1.04
Output  | 2.8849e-05 | 2.8849e-05 | 2.8849e-05 |   0.0 |  0.00
Modify  | 0.22629    | 0.22629    | 0.22629    |   0.0 |  1.77
Other   |            | 0.05205    |            |       |  0.41

Nlocal:    2000 ave 2000 max 2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    5471 ave 5471 max 5471 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:  336034 ave 336034 max 336034 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 336034
Ave neighs/atom = 168.017
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.947471794546, Press = -1.45388455043757
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 4.577 | 4.577 | 4.577 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   14000   -16389.839   -16389.839   -16473.145   -16473.145    322.40149    322.40149    32338.322    32338.322   -3465.5314   -3465.5314 
   15000   -16392.818   -16392.818   -16472.281   -16472.281    307.53292    307.53292    32268.834    32268.834    3345.4193    3345.4193 
Loop time of 12.4777 on 1 procs for 1000 steps with 2000 atoms

Performance: 6.924 ns/day, 3.466 hours/ns, 80.143 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    | 12.247     | 12.247     | 12.247     |   0.0 | 98.15
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.052796   | 0.052796   | 0.052796   |   0.0 |  0.42
Output  | 2.2888e-05 | 2.2888e-05 | 2.2888e-05 |   0.0 |  0.00
Modify  | 0.14642    | 0.14642    | 0.14642    |   0.0 |  1.17
Other   |            | 0.03193    |            |       |  0.26

Nlocal:    2000 ave 2000 max 2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    5471 ave 5471 max 5471 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:  336028 ave 336028 max 336028 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 336028
Ave neighs/atom = 168.014
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.062165910895, Press = -2.1144747802141
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 4.577 | 4.577 | 4.577 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   15000   -16392.818   -16392.818   -16472.281   -16472.281    307.53292    307.53292    32268.834    32268.834    3345.4193    3345.4193 
   16000   -16396.954   -16396.954   -16474.209   -16474.209     298.9841     298.9841    32325.996    32325.996   -2618.6232   -2618.6232 
Loop time of 12.7314 on 1 procs for 1000 steps with 2000 atoms

Performance: 6.786 ns/day, 3.536 hours/ns, 78.546 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    | 12.481     | 12.481     | 12.481     |   0.0 | 98.03
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.072532   | 0.072532   | 0.072532   |   0.0 |  0.57
Output  | 2.8849e-05 | 2.8849e-05 | 2.8849e-05 |   0.0 |  0.00
Modify  | 0.16613    | 0.16613    | 0.16613    |   0.0 |  1.30
Other   |            | 0.01195    |            |       |  0.09

Nlocal:    2000 ave 2000 max 2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    5471 ave 5471 max 5471 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:  336038 ave 336038 max 336038 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 336038
Ave neighs/atom = 168.019
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.046248202554, Press = -0.545750735596884
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 4.577 | 4.577 | 4.577 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   16000   -16396.954   -16396.954   -16474.209   -16474.209     298.9841     298.9841    32325.996    32325.996   -2618.6232   -2618.6232 
   17000   -16392.092   -16392.092   -16473.066   -16473.066    313.37805    313.37805    32276.025    32276.025    2624.0797    2624.0797 
Loop time of 11.8284 on 1 procs for 1000 steps with 2000 atoms

Performance: 7.304 ns/day, 3.286 hours/ns, 84.543 timesteps/s
41.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    | 11.513     | 11.513     | 11.513     |   0.0 | 97.33
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.11245    | 0.11245    | 0.11245    |   0.0 |  0.95
Output  | 2.7895e-05 | 2.7895e-05 | 2.7895e-05 |   0.0 |  0.00
Modify  | 0.19137    | 0.19137    | 0.19137    |   0.0 |  1.62
Other   |            | 0.01181    |            |       |  0.10

Nlocal:    2000 ave 2000 max 2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    5471 ave 5471 max 5471 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:  336036 ave 336036 max 336036 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 336036
Ave neighs/atom = 168.018
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.842721691712, Press = -1.54267969247286
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 4.577 | 4.577 | 4.577 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   17000   -16392.092   -16392.092   -16473.066   -16473.066    313.37805    313.37805    32276.025    32276.025    2624.0797    2624.0797 
   18000   -16397.901   -16397.901   -16477.797   -16477.797    309.20395    309.20395    32326.859    32326.859   -2963.9763   -2963.9763 
Loop time of 12.9812 on 1 procs for 1000 steps with 2000 atoms

Performance: 6.656 ns/day, 3.606 hours/ns, 77.034 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    | 12.506     | 12.506     | 12.506     |   0.0 | 96.34
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.13297    | 0.13297    | 0.13297    |   0.0 |  1.02
Output  | 2.718e-05  | 2.718e-05  | 2.718e-05  |   0.0 |  0.00
Modify  | 0.31004    | 0.31004    | 0.31004    |   0.0 |  2.39
Other   |            | 0.03208    |            |       |  0.25

Nlocal:    2000 ave 2000 max 2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    5471 ave 5471 max 5471 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:  336036 ave 336036 max 336036 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 336036
Ave neighs/atom = 168.018
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.752948400739, Press = -0.867801490827813
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 4.577 | 4.577 | 4.577 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   18000   -16397.901   -16397.901   -16477.797   -16477.797    309.20395    309.20395    32326.859    32326.859   -2963.9763   -2963.9763 
   19000   -16392.751   -16392.751   -16474.441   -16474.441    316.14882    316.14882      32280.6      32280.6    1976.7987    1976.7987 
Loop time of 12.8681 on 1 procs for 1000 steps with 2000 atoms

Performance: 6.714 ns/day, 3.574 hours/ns, 77.712 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    | 12.448     | 12.448     | 12.448     |   0.0 | 96.73
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.14267    | 0.14267    | 0.14267    |   0.0 |  1.11
Output  | 2.3127e-05 | 2.3127e-05 | 2.3127e-05 |   0.0 |  0.00
Modify  | 0.24594    | 0.24594    | 0.24594    |   0.0 |  1.91
Other   |            | 0.03193    |            |       |  0.25

Nlocal:    2000 ave 2000 max 2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    5471 ave 5471 max 5471 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:  336016 ave 336016 max 336016 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 336016
Ave neighs/atom = 168.008
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.635908669573, Press = -1.2757138014784
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 4.577 | 4.577 | 4.577 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   19000   -16392.751   -16392.751   -16474.441   -16474.441    316.14882    316.14882      32280.6      32280.6    1976.7987    1976.7987 
   20000   -16394.407   -16394.407   -16475.859   -16475.859    315.22809    315.22809    32313.182    32313.182   -1321.1817   -1321.1817 
Loop time of 12.2388 on 1 procs for 1000 steps with 2000 atoms

Performance: 7.060 ns/day, 3.400 hours/ns, 81.707 timesteps/s
40.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    | 11.838     | 11.838     | 11.838     |   0.0 | 96.73
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.052449   | 0.052449   | 0.052449   |   0.0 |  0.43
Output  | 2.2888e-05 | 2.2888e-05 | 2.2888e-05 |   0.0 |  0.00
Modify  | 0.31594    | 0.31594    | 0.31594    |   0.0 |  2.58
Other   |            | 0.03221    |            |       |  0.26

Nlocal:    2000 ave 2000 max 2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    5471 ave 5471 max 5471 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:  336026 ave 336026 max 336026 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 336026
Ave neighs/atom = 168.013
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.721468576486, Press = -1.22380065960516
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 4.577 | 4.577 | 4.577 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   20000   -16394.407   -16394.407   -16475.859   -16475.859    315.22809    315.22809    32313.182    32313.182   -1321.1817   -1321.1817 
   21000   -16387.784   -16387.784   -16470.277   -16470.277    319.25378    319.25378    32296.046    32296.046    990.53371    990.53371 
Loop time of 13.0036 on 1 procs for 1000 steps with 2000 atoms

Performance: 6.644 ns/day, 3.612 hours/ns, 76.902 timesteps/s
37.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    | 12.675     | 12.675     | 12.675     |   0.0 | 97.47
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.07296    | 0.07296    | 0.07296    |   0.0 |  0.56
Output  | 2.3127e-05 | 2.3127e-05 | 2.3127e-05 |   0.0 |  0.00
Modify  | 0.1835     | 0.1835     | 0.1835     |   0.0 |  1.41
Other   |            | 0.07203    |            |       |  0.55

Nlocal:    2000 ave 2000 max 2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    5471 ave 5471 max 5471 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:  336030 ave 336030 max 336030 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 336030
Ave neighs/atom = 168.015
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.895663149598, Press = -0.774514929058308
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 4.577 | 4.577 | 4.577 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   21000   -16387.784   -16387.784   -16470.277   -16470.277    319.25378    319.25378    32296.046    32296.046    990.53371    990.53371 
   22000   -16395.935   -16395.935   -16476.646   -16476.646    312.35809    312.35809    32312.286    32312.286   -1429.2015   -1429.2015 
Loop time of 12.0383 on 1 procs for 1000 steps with 2000 atoms

Performance: 7.177 ns/day, 3.344 hours/ns, 83.068 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    | 11.758     | 11.758     | 11.758     |   0.0 | 97.67
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.042418   | 0.042418   | 0.042418   |   0.0 |  0.35
Output  | 4.4823e-05 | 4.4823e-05 | 4.4823e-05 |   0.0 |  0.00
Modify  | 0.20567    | 0.20567    | 0.20567    |   0.0 |  1.71
Other   |            | 0.0319     |            |       |  0.26

Nlocal:    2000 ave 2000 max 2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    5471 ave 5471 max 5471 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:  336026 ave 336026 max 336026 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 336026
Ave neighs/atom = 168.013
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.980393223368, Press = -0.874531791084502
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 4.577 | 4.577 | 4.577 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   22000   -16395.935   -16395.935   -16476.646   -16476.646    312.35809    312.35809    32312.286    32312.286   -1429.2015   -1429.2015 
   23000   -16394.739   -16394.739   -16475.783   -16475.783    313.64731    313.64731     32273.06     32273.06    2587.3976    2587.3976 
Loop time of 12.38 on 1 procs for 1000 steps with 2000 atoms

Performance: 6.979 ns/day, 3.439 hours/ns, 80.776 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    | 12.14      | 12.14      | 12.14      |   0.0 | 98.06
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.052267   | 0.052267   | 0.052267   |   0.0 |  0.42
Output  | 3.5048e-05 | 3.5048e-05 | 3.5048e-05 |   0.0 |  0.00
Modify  | 0.15563    | 0.15563    | 0.15563    |   0.0 |  1.26
Other   |            | 0.0319     |            |       |  0.26

Nlocal:    2000 ave 2000 max 2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    5471 ave 5471 max 5471 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:  336026 ave 336026 max 336026 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 336026
Ave neighs/atom = 168.013
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.035988266541, Press = -1.43841538679115
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 4.577 | 4.577 | 4.577 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   23000   -16394.739   -16394.739   -16475.783   -16475.783    313.64731    313.64731     32273.06     32273.06    2587.3976    2587.3976 
   24000    -16393.02    -16393.02   -16474.831   -16474.831    316.61954    316.61954    32334.919    32334.919    -3365.857    -3365.857 
Loop time of 12.347 on 1 procs for 1000 steps with 2000 atoms

Performance: 6.998 ns/day, 3.430 hours/ns, 80.991 timesteps/s
40.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    | 12.016     | 12.016     | 12.016     |   0.0 | 97.32
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.072955   | 0.072955   | 0.072955   |   0.0 |  0.59
Output  | 2.1935e-05 | 2.1935e-05 | 2.1935e-05 |   0.0 |  0.00
Modify  | 0.22617    | 0.22617    | 0.22617    |   0.0 |  1.83
Other   |            | 0.03193    |            |       |  0.26

Nlocal:    2000 ave 2000 max 2000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    5471 ave 5471 max 5471 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:  336026 ave 336026 max 336026 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 336026
Ave neighs/atom = 168.013
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
32300.6457513447

print "LAMMPS calculation completed"
LAMMPS calculation completed
quit 0