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

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

# specify which KIM Model to use
pair_style    kim EMT_Asap_Standard_JacobsenStoltzeNorskov_1996_Ag__MO_303974873468_001
pair_coeff    * * Ag
mass          1 ${mass_converted}
mass          1 107.8682

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

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

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

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

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

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

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

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

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

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

# Run a simulation for at most 2000*1000 timesteps. At each 1000th time step, check
# whether the temperature and pressure have converged. If yes, break.
label top
variable a loop 2000
run 1000
Neighbor list info ...
  update every 1 steps, delay 10 steps, check yes
  max neighbors/atom: 2000, page size: 100000
  master list distance cutoff = 7.61876
  ghost atom cutoff = 7.61876
  binsize = 3.80938, bins = 11 11 11
  1 neighbor lists, perpetual/occasional/extra = 1 0 0
  (1) pair kim, perpetual
      attributes: full, newton off, cut 7.61876
      pair build: full/bin/atomonly
      stencil: full/bin/3d
      bin: standard
Per MPI rank memory allocation (min/avg/max) = 5.369 | 5.369 | 5.369 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
       0   -11689.934   -11689.934   -11841.467   -11841.467       293.15       293.15    67099.313    67099.313     2412.161     2412.161 
    1000   -11517.727   -11517.727   -11673.764   -11673.764    301.86338    301.86338      68663.3      68663.3   -238.23318   -238.23318 
Loop time of 80.511 on 1 procs for 1000 steps with 4000 atoms

Performance: 1.073 ns/day, 22.364 hours/ns, 12.421 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    | 79.829     | 79.829     | 79.829     |   0.0 | 99.15
Neigh   | 0.034553   | 0.034553   | 0.034553   |   0.0 |  0.04
Comm    | 0.11939    | 0.11939    | 0.11939    |   0.0 |  0.15
Output  | 3.7193e-05 | 3.7193e-05 | 3.7193e-05 |   0.0 |  0.00
Modify  | 0.46765    | 0.46765    | 0.46765    |   0.0 |  0.58
Other   |            | 0.06075    |            |       |  0.08

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

Total # of neighbors = 414884
Ave neighs/atom = 103.721
Neighbor list builds = 1
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) = 5.369 | 5.369 | 5.369 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    1000   -11517.727   -11517.727   -11673.764   -11673.764    301.86338    301.86338      68663.3      68663.3   -238.23318   -238.23318 
    2000   -11539.417   -11539.417   -11685.807   -11685.807    283.20243    283.20243    68547.322    68547.322   -155.28727   -155.28727 
Loop time of 84.4567 on 1 procs for 1000 steps with 4000 atoms

Performance: 1.023 ns/day, 23.460 hours/ns, 11.840 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    | 83.919     | 83.919     | 83.919     |   0.0 | 99.36
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.11982    | 0.11982    | 0.11982    |   0.0 |  0.14
Output  | 3.4094e-05 | 3.4094e-05 | 3.4094e-05 |   0.0 |  0.00
Modify  | 0.33755    | 0.33755    | 0.33755    |   0.0 |  0.40
Other   |            | 0.08056    |            |       |  0.10

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

Total # of neighbors = 414842
Ave neighs/atom = 103.71
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) = 5.369 | 5.369 | 5.369 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    2000   -11539.417   -11539.417   -11685.807   -11685.807    283.20243    283.20243    68547.322    68547.322   -155.28727   -155.28727 
    3000   -11527.316   -11527.316   -11681.021   -11681.021    297.35313    297.35313    68506.113    68506.113    1133.1569    1133.1569 
Loop time of 84.45 on 1 procs for 1000 steps with 4000 atoms

Performance: 1.023 ns/day, 23.458 hours/ns, 11.841 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    | 83.844     | 83.844     | 83.844     |   0.0 | 99.28
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.098744   | 0.098744   | 0.098744   |   0.0 |  0.12
Output  | 3.6955e-05 | 3.6955e-05 | 3.6955e-05 |   0.0 |  0.00
Modify  | 0.44663    | 0.44663    | 0.44663    |   0.0 |  0.53
Other   |            | 0.06044    |            |       |  0.07

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

Total # of neighbors = 416090
Ave neighs/atom = 104.022
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) = 5.369 | 5.369 | 5.369 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    3000   -11527.316   -11527.316   -11681.021   -11681.021    297.35313    297.35313    68506.113    68506.113    1133.1569    1133.1569 
    4000   -11532.144   -11532.144   -11683.875   -11683.875    293.53387    293.53387    68576.655    68576.655   -209.20047   -209.20047 
Loop time of 84.1712 on 1 procs for 1000 steps with 4000 atoms

Performance: 1.026 ns/day, 23.381 hours/ns, 11.881 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    | 83.491     | 83.491     | 83.491     |   0.0 | 99.19
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.059724   | 0.059724   | 0.059724   |   0.0 |  0.07
Output  | 4.2915e-05 | 4.2915e-05 | 4.2915e-05 |   0.0 |  0.00
Modify  | 0.51931    | 0.51931    | 0.51931    |   0.0 |  0.62
Other   |            | 0.1008     |            |       |  0.12

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

Total # of neighbors = 417198
Ave neighs/atom = 104.299
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) = 5.369 | 5.369 | 5.369 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    4000   -11532.144   -11532.144   -11683.875   -11683.875    293.53387    293.53387    68576.655    68576.655   -209.20047   -209.20047 
    5000   -11532.223   -11532.223   -11683.701   -11683.701    293.04446    293.04446    68633.706    68633.706   -1009.4231   -1009.4231 
Loop time of 85.2823 on 1 procs for 1000 steps with 4000 atoms

Performance: 1.013 ns/day, 23.690 hours/ns, 11.726 timesteps/s
33.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    | 84.634     | 84.634     | 84.634     |   0.0 | 99.24
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.11934    | 0.11934    | 0.11934    |   0.0 |  0.14
Output  | 2.9802e-05 | 2.9802e-05 | 2.9802e-05 |   0.0 |  0.00
Modify  | 0.50827    | 0.50827    | 0.50827    |   0.0 |  0.60
Other   |            | 0.02068    |            |       |  0.02

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

Total # of neighbors = 415260
Ave neighs/atom = 103.815
Neighbor list builds = 0
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 292.330068201002, Press = 474.767170819499
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 5.369 | 5.369 | 5.369 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    5000   -11532.223   -11532.223   -11683.701   -11683.701    293.04446    293.04446    68633.706    68633.706   -1009.4231   -1009.4231 
    6000    -11530.33    -11530.33   -11683.873   -11683.873    297.04063    297.04063    68606.626    68606.626   -636.97084   -636.97084 
Loop time of 80.8522 on 1 procs for 1000 steps with 4000 atoms

Performance: 1.069 ns/day, 22.459 hours/ns, 12.368 timesteps/s
34.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    | 80.127     | 80.127     | 80.127     |   0.0 | 99.10
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.14398    | 0.14398    | 0.14398    |   0.0 |  0.18
Output  | 3.4094e-05 | 3.4094e-05 | 3.4094e-05 |   0.0 |  0.00
Modify  | 0.54114    | 0.54114    | 0.54114    |   0.0 |  0.67
Other   |            | 0.0403     |            |       |  0.05

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

Total # of neighbors = 414372
Ave neighs/atom = 103.593
Neighbor list builds = 0
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 293.215554376203, Press = 28.5491677751056
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 5.369 | 5.369 | 5.369 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    6000    -11530.33    -11530.33   -11683.873   -11683.873    297.04063    297.04063    68606.626    68606.626   -636.97084   -636.97084 
    7000   -11534.289   -11534.289   -11684.792   -11684.792    291.15744    291.15744    68504.613    68504.613    642.26842    642.26842 
Loop time of 80.9337 on 1 procs for 1000 steps with 4000 atoms

Performance: 1.068 ns/day, 22.482 hours/ns, 12.356 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    | 80.171     | 80.171     | 80.171     |   0.0 | 99.06
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.11919    | 0.11919    | 0.11919    |   0.0 |  0.15
Output  | 2.6941e-05 | 2.6941e-05 | 2.6941e-05 |   0.0 |  0.00
Modify  | 0.61366    | 0.61366    | 0.61366    |   0.0 |  0.76
Other   |            | 0.03012    |            |       |  0.04

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

Total # of neighbors = 414972
Ave neighs/atom = 103.743
Neighbor list builds = 0
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 293.150771129663, Press = 25.1081391534555
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 5.369 | 5.369 | 5.369 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    7000   -11534.289   -11534.289   -11684.792   -11684.792    291.15744    291.15744    68504.613    68504.613    642.26842    642.26842 
    8000   -11530.035   -11530.035   -11679.102   -11679.102    288.38087    288.38087    68551.532    68551.532    673.37454    673.37454 
Loop time of 78.9686 on 1 procs for 1000 steps with 4000 atoms

Performance: 1.094 ns/day, 21.936 hours/ns, 12.663 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    | 78.178     | 78.178     | 78.178     |   0.0 | 99.00
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.15911    | 0.15911    | 0.15911    |   0.0 |  0.20
Output  | 3.3855e-05 | 3.3855e-05 | 3.3855e-05 |   0.0 |  0.00
Modify  | 0.5309     | 0.5309     | 0.5309     |   0.0 |  0.67
Other   |            | 0.1005     |            |       |  0.13

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

Total # of neighbors = 416770
Ave neighs/atom = 104.192
Neighbor list builds = 0
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 293.028946116623, Press = 4.73152389623824
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 5.369 | 5.369 | 5.369 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    8000   -11530.035   -11530.035   -11679.102   -11679.102    288.38087    288.38087    68551.532    68551.532    673.37454    673.37454 
    9000    -11535.99    -11535.99   -11685.799   -11685.799    289.81508    289.81508    68565.078    68565.078   -286.76081   -286.76081 
Loop time of 78.6121 on 1 procs for 1000 steps with 4000 atoms

Performance: 1.099 ns/day, 21.837 hours/ns, 12.721 timesteps/s
35.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    | 77.839     | 77.839     | 77.839     |   0.0 | 99.02
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.099366   | 0.099366   | 0.099366   |   0.0 |  0.13
Output  | 4.1008e-05 | 4.1008e-05 | 4.1008e-05 |   0.0 |  0.00
Modify  | 0.55299    | 0.55299    | 0.55299    |   0.0 |  0.70
Other   |            | 0.1206     |            |       |  0.15

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

Total # of neighbors = 416376
Ave neighs/atom = 104.094
Neighbor list builds = 0
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 292.711236676015, Press = 2.93739540474553
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 5.369 | 5.369 | 5.369 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
    9000    -11535.99    -11535.99   -11685.799   -11685.799    289.81508    289.81508    68565.078    68565.078   -286.76081   -286.76081 
   10000   -11526.608   -11526.608   -11680.132   -11680.132     297.0032     297.0032    68527.353    68527.353    993.86966    993.86966 
Loop time of 77.3478 on 1 procs for 1000 steps with 4000 atoms

Performance: 1.117 ns/day, 21.485 hours/ns, 12.929 timesteps/s
36.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    | 76.818     | 76.818     | 76.818     |   0.0 | 99.32
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.069115   | 0.069115   | 0.069115   |   0.0 |  0.09
Output  | 6.0081e-05 | 6.0081e-05 | 6.0081e-05 |   0.0 |  0.00
Modify  | 0.40018    | 0.40018    | 0.40018    |   0.0 |  0.52
Other   |            | 0.06015    |            |       |  0.08

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

Total # of neighbors = 415322
Ave neighs/atom = 103.831
Neighbor list builds = 0
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 292.847015115707, Press = 6.11407681994124
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 5.369 | 5.369 | 5.369 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   10000   -11526.608   -11526.608   -11680.132   -11680.132     297.0032     297.0032    68527.353    68527.353    993.86966    993.86966 
   11000   -11530.977   -11530.977   -11682.104   -11682.104    292.36605    292.36605    68668.248    68668.248   -1290.5946   -1290.5946 
Loop time of 71.9556 on 1 procs for 1000 steps with 4000 atoms

Performance: 1.201 ns/day, 19.988 hours/ns, 13.897 timesteps/s
38.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    | 71.466     | 71.466     | 71.466     |   0.0 | 99.32
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.078867   | 0.078867   | 0.078867   |   0.0 |  0.11
Output  | 3.1948e-05 | 3.1948e-05 | 3.1948e-05 |   0.0 |  0.00
Modify  | 0.34026    | 0.34026    | 0.34026    |   0.0 |  0.47
Other   |            | 0.07039    |            |       |  0.10

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

Total # of neighbors = 417002
Ave neighs/atom = 104.251
Neighbor list builds = 0
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 293.130675451846, Press = 7.442260420176
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 5.369 | 5.369 | 5.369 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   11000   -11530.977   -11530.977   -11682.104   -11682.104    292.36605    292.36605    68668.248    68668.248   -1290.5946   -1290.5946 
   12000   -11532.436   -11532.436   -11685.303   -11685.303    295.73174    295.73174    68512.614    68512.614    515.22273    515.22273 
Loop time of 74.3268 on 1 procs for 1000 steps with 4000 atoms

Performance: 1.162 ns/day, 20.646 hours/ns, 13.454 timesteps/s
37.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    | 73.727     | 73.727     | 73.727     |   0.0 | 99.19
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.11915    | 0.11915    | 0.11915    |   0.0 |  0.16
Output  | 3.314e-05  | 3.314e-05  | 3.314e-05  |   0.0 |  0.00
Modify  | 0.40053    | 0.40053    | 0.40053    |   0.0 |  0.54
Other   |            | 0.08037    |            |       |  0.11

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

Total # of neighbors = 413560
Ave neighs/atom = 103.39
Neighbor list builds = 0
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 292.889214006453, Press = 2.26356959352529
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 5.369 | 5.369 | 5.369 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   12000   -11532.436   -11532.436   -11685.303   -11685.303    295.73174    295.73174    68512.614    68512.614    515.22273    515.22273 
   13000   -11530.324   -11530.324   -11682.359   -11682.359    294.12086    294.12086    68635.714    68635.714   -893.70726   -893.70726 
Loop time of 71.8148 on 1 procs for 1000 steps with 4000 atoms

Performance: 1.203 ns/day, 19.949 hours/ns, 13.925 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    | 71.116     | 71.116     | 71.116     |   0.0 | 99.03
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.15852    | 0.15852    | 0.15852    |   0.0 |  0.22
Output  | 3.2902e-05 | 3.2902e-05 | 3.2902e-05 |   0.0 |  0.00
Modify  | 0.52       | 0.52       | 0.52       |   0.0 |  0.72
Other   |            | 0.02031    |            |       |  0.03

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

Total # of neighbors = 416982
Ave neighs/atom = 104.246
Neighbor list builds = 0
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 292.863176682892, Press = 4.56082765010097
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 5.369 | 5.369 | 5.369 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   13000   -11530.324   -11530.324   -11682.359   -11682.359    294.12086    294.12086    68635.714    68635.714   -893.70726   -893.70726 
   14000   -11530.071   -11530.071   -11681.584   -11681.584     293.1115     293.1115    68500.137    68500.137    1051.2866    1051.2866 
Loop time of 74.315 on 1 procs for 1000 steps with 4000 atoms

Performance: 1.163 ns/day, 20.643 hours/ns, 13.456 timesteps/s
37.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    | 73.764     | 73.764     | 73.764     |   0.0 | 99.26
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.09931    | 0.09931    | 0.09931    |   0.0 |  0.13
Output  | 3.2902e-05 | 3.2902e-05 | 3.2902e-05 |   0.0 |  0.00
Modify  | 0.39147    | 0.39147    | 0.39147    |   0.0 |  0.53
Other   |            | 0.06059    |            |       |  0.08

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

Total # of neighbors = 414302
Ave neighs/atom = 103.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 = 292.914194335286, Press = 1.96409601119542
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 5.369 | 5.369 | 5.369 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   14000   -11530.071   -11530.071   -11681.584   -11681.584     293.1115     293.1115    68500.137    68500.137    1051.2866    1051.2866 
   15000    -11533.36    -11533.36   -11683.857   -11683.857    291.14848    291.14848    68668.306    68668.306   -1585.7416   -1585.7416 
Loop time of 71.8317 on 1 procs for 1000 steps with 4000 atoms

Performance: 1.203 ns/day, 19.953 hours/ns, 13.921 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    | 71.312     | 71.312     | 71.312     |   0.0 | 99.28
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.11849    | 0.11849    | 0.11849    |   0.0 |  0.16
Output  | 3.2902e-05 | 3.2902e-05 | 3.2902e-05 |   0.0 |  0.00
Modify  | 0.38054    | 0.38054    | 0.38054    |   0.0 |  0.53
Other   |            | 0.02025    |            |       |  0.03

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

Total # of neighbors = 417598
Ave neighs/atom = 104.4
Neighbor list builds = 0
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 292.924652481563, Press = 2.37995515412684
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 5.369 | 5.369 | 5.369 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   15000    -11533.36    -11533.36   -11683.857   -11683.857    291.14848    291.14848    68668.306    68668.306   -1585.7416   -1585.7416 
   16000   -11530.775   -11530.775   -11682.751   -11682.751    294.00905    294.00905    68524.343    68524.343    625.85901    625.85901 
Loop time of 69.2891 on 1 procs for 1000 steps with 4000 atoms

Performance: 1.247 ns/day, 19.247 hours/ns, 14.432 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    | 68.86      | 68.86      | 68.86      |   0.0 | 99.38
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.11865    | 0.11865    | 0.11865    |   0.0 |  0.17
Output  | 3.2902e-05 | 3.2902e-05 | 3.2902e-05 |   0.0 |  0.00
Modify  | 0.29019    | 0.29019    | 0.29019    |   0.0 |  0.42
Other   |            | 0.02009    |            |       |  0.03

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

Total # of neighbors = 413570
Ave neighs/atom = 103.392
Neighbor list builds = 0
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 292.98992637404, Press = 2.91072977121248
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 5.369 | 5.369 | 5.369 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   16000   -11530.775   -11530.775   -11682.751   -11682.751    294.00905    294.00905    68524.343    68524.343    625.85901    625.85901 
   17000   -11532.322   -11532.322    -11683.42    -11683.42    292.30989    292.30989    68584.977    68584.977   -273.30407   -273.30407 
Loop time of 66.5892 on 1 procs for 1000 steps with 4000 atoms

Performance: 1.298 ns/day, 18.497 hours/ns, 15.017 timesteps/s
41.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    | 66.15      | 66.15      | 66.15      |   0.0 | 99.34
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.038713   | 0.038713   | 0.038713   |   0.0 |  0.06
Output  | 6.8188e-05 | 6.8188e-05 | 6.8188e-05 |   0.0 |  0.00
Modify  | 0.35986    | 0.35986    | 0.35986    |   0.0 |  0.54
Other   |            | 0.04023    |            |       |  0.06

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

Total # of neighbors = 416980
Ave neighs/atom = 104.245
Neighbor list builds = 0
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 293.075926121596, Press = 0.399906504641818
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 5.369 | 5.369 | 5.369 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   17000   -11532.322   -11532.322    -11683.42    -11683.42    292.30989    292.30989    68584.977    68584.977   -273.30407   -273.30407 
   18000   -11532.397   -11532.397   -11684.118   -11684.118    293.51591    293.51591    68486.455    68486.455    1016.4915    1016.4915 
Loop time of 78.3197 on 1 procs for 1000 steps with 4000 atoms

Performance: 1.103 ns/day, 21.755 hours/ns, 12.768 timesteps/s
35.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    | 77.645     | 77.645     | 77.645     |   0.0 | 99.14
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.10584    | 0.10584    | 0.10584    |   0.0 |  0.14
Output  | 3.4094e-05 | 3.4094e-05 | 3.4094e-05 |   0.0 |  0.00
Modify  | 0.52908    | 0.52908    | 0.52908    |   0.0 |  0.68
Other   |            | 0.04002    |            |       |  0.05

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

Total # of neighbors = 415668
Ave neighs/atom = 103.917
Neighbor list builds = 0
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 293.090937306953, Press = 4.24566105364353
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 5.369 | 5.369 | 5.369 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   18000   -11532.397   -11532.397   -11684.118   -11684.118    293.51591    293.51591    68486.455    68486.455    1016.4915    1016.4915 
   19000    -11528.86    -11528.86   -11680.657   -11680.657    293.66155    293.66155    68565.711    68565.711    276.27935    276.27935 
Loop time of 78.8355 on 1 procs for 1000 steps with 4000 atoms

Performance: 1.096 ns/day, 21.899 hours/ns, 12.685 timesteps/s
35.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    | 78.386     | 78.386     | 78.386     |   0.0 | 99.43
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.089029   | 0.089029   | 0.089029   |   0.0 |  0.11
Output  | 3.4094e-05 | 3.4094e-05 | 3.4094e-05 |   0.0 |  0.00
Modify  | 0.34042    | 0.34042    | 0.34042    |   0.0 |  0.43
Other   |            | 0.02003    |            |       |  0.03

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

Total # of neighbors = 417318
Ave neighs/atom = 104.329
Neighbor list builds = 0
Dangerous builds = 0
if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable"
if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break"
print "flag: Temp = ${T}, Press = ${P}"
flag: Temp = 293.16156295653, Press = 2.08026552976277
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 5.369 | 5.369 | 5.369 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   19000    -11528.86    -11528.86   -11680.657   -11680.657    293.66155    293.66155    68565.711    68565.711    276.27935    276.27935 
   20000   -11532.424   -11532.424   -11684.566   -11684.566    294.32787    294.32787    68536.182    68536.182    256.55018    256.55018 
Loop time of 79.0384 on 1 procs for 1000 steps with 4000 atoms

Performance: 1.093 ns/day, 21.955 hours/ns, 12.652 timesteps/s
35.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    | 78.339     | 78.339     | 78.339     |   0.0 | 99.12
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.13885    | 0.13885    | 0.13885    |   0.0 |  0.18
Output  | 2.5988e-05 | 2.5988e-05 | 2.5988e-05 |   0.0 |  0.00
Modify  | 0.50016    | 0.50016    | 0.50016    |   0.0 |  0.63
Other   |            | 0.06036    |            |       |  0.08

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

Total # of neighbors = 416060
Ave neighs/atom = 104.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 = 293.230490588838, Press = 3.9834501927005
next a
jump SELF top
variable a loop 2000
run 1000
Per MPI rank memory allocation (min/avg/max) = 5.369 | 5.369 | 5.369 Mbytes
Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 
   20000   -11532.424   -11532.424   -11684.566   -11684.566    294.32787    294.32787    68536.182    68536.182    256.55018    256.55018 
   21000   -11528.895   -11528.895   -11678.934   -11678.934    290.25999    290.25999    68594.016    68594.016    107.87126    107.87126 
Loop time of 78.8461 on 1 procs for 1000 steps with 4000 atoms

Performance: 1.096 ns/day, 21.902 hours/ns, 12.683 timesteps/s
35.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    | 78.251     | 78.251     | 78.251     |   0.0 | 99.25
Neigh   | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.11892    | 0.11892    | 0.11892    |   0.0 |  0.15
Output  | 3.1948e-05 | 3.1948e-05 | 3.1948e-05 |   0.0 |  0.00
Modify  | 0.41522    | 0.41522    | 0.41522    |   0.0 |  0.53
Other   |            | 0.06067    |            |       |  0.08

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

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

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

print "LAMMPS calculation completed"
LAMMPS calculation completed
quit 0