# Variables that can be adjusted by kim-lammps-preprocessor to switch unit sets for # Simulator Models variable _u_distance equal 1.0 variable _u_energy equal 1.0 variable _u_mass equal 1.0 variable _u_time equal 1.0 variable _u_pressure equal 1.0 variable _u_temperature equal 1.0 # This line may be swapped out by kim-lammps-preprocessor if running against a Simulator # Model whose atom_style is not 'atomic' atom_style atomic # periodic boundary conditions along all three dimensions boundary p p p # Set neighbor skin variable neigh_skin equal 2.0*${_u_distance} variable neigh_skin equal 2.0*1 neighbor ${neigh_skin} bin neighbor 2 bin # create a supercell with cubic lattice (fcc, bcc, sc, or diamond) # using 10*10*10 conventional (orthogonal) unit cells variable latticeconst_converted equal 6.126234233379364*${_u_distance} variable latticeconst_converted equal 6.126234233379364*1 lattice fcc ${latticeconst_converted} lattice fcc 6.12623423337936 Lattice spacing in x,y,z = 6.12623 6.12623 6.12623 region simbox block 0 10 0 10 0 10 units lattice create_box 1 simbox Created orthogonal box = (0 0 0) to (61.2623 61.2623 61.2623) 1 by 1 by 1 MPI processor grid create_atoms 1 box Created 4000 atoms create_atoms CPU = 0.000381947 secs variable mass_converted equal 87.62*${_u_mass} variable mass_converted equal 87.62*1 # specify which KIM Model to use pair_style kim Morse_Shifted_GirifalcoWeizer_1959LowCutoff_Sr__MO_801083489225_003 WARNING: KIM Model does not provide `partialParticleVirial'; virial per atom will be zero (src/KIM/pair_kim.cpp:979) pair_coeff * * Sr mass 1 ${mass_converted} mass 1 87.62 # initial volume variable v equal vol # assign formula variable V0 equal ${v} # evaluate initial value variable V0 equal 229922.140227974 variable V0_metal equal ${V0}/(${_u_distance}*${_u_distance}*${_u_distance}) variable V0_metal equal 229922.140227974/(${_u_distance}*${_u_distance}*${_u_distance}) variable V0_metal equal 229922.140227974/(1*${_u_distance}*${_u_distance}) variable V0_metal equal 229922.140227974/(1*1*${_u_distance}) variable V0_metal equal 229922.140227974/(1*1*1) variable V0_metal_times1000 equal ${V0_metal}*1000 variable V0_metal_times1000 equal 229922.140227974*1000 print "Initial system volume: ${V0_metal} Angstroms^3" Initial system volume: 229922.140227974 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 = 14.2649 ghost atom cutoff = 14.2649 binsize = 7.13245, bins = 9 9 9 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 14.2649 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 6.623 | 6.623 | 6.623 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 0 -6250.8786 -6250.8786 -6402.4111 -6402.4111 293.15 293.15 229922.14 229922.14 703.95166 703.95166 1000 -6081.654 -6081.654 -6237.2837 -6237.2837 301.07629 301.07629 236144.47 236144.47 -265.75929 -265.75929 Loop time of 16.1152 on 1 procs for 1000 steps with 4000 atoms Performance: 5.361 ns/day, 4.476 hours/ns, 62.053 timesteps/s 99.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 | 15.761 | 15.761 | 15.761 | 0.0 | 97.80 Neigh | 0.14116 | 0.14116 | 0.14116 | 0.0 | 0.88 Comm | 0.04906 | 0.04906 | 0.04906 | 0.0 | 0.30 Output | 4.3869e-05 | 4.3869e-05 | 4.3869e-05 | 0.0 | 0.00 Modify | 0.14523 | 0.14523 | 0.14523 | 0.0 | 0.90 Other | | 0.01858 | | | 0.12 Nlocal: 4000 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8195 ave 8195 max 8195 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 813474 ave 813474 max 813474 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 813474 Ave neighs/atom = 203.369 Neighbor list builds = 6 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 0, Press = 0 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 6.623 | 6.623 | 6.623 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 1000 -6081.654 -6081.654 -6237.2837 -6237.2837 301.07629 301.07629 236144.47 236144.47 -265.75929 -265.75929 2000 -6089.1657 -6089.1657 -6245.12 -6245.12 301.70429 301.70429 235286.39 235286.39 50.714994 50.714994 Loop time of 18.3306 on 1 procs for 1000 steps with 4000 atoms Performance: 4.713 ns/day, 5.092 hours/ns, 54.554 timesteps/s 92.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 | 17.931 | 17.931 | 17.931 | 0.0 | 97.82 Neigh | 0.13868 | 0.13868 | 0.13868 | 0.0 | 0.76 Comm | 0.064695 | 0.064695 | 0.064695 | 0.0 | 0.35 Output | 0.0026259 | 0.0026259 | 0.0026259 | 0.0 | 0.01 Modify | 0.16284 | 0.16284 | 0.16284 | 0.0 | 0.89 Other | | 0.03101 | | | 0.17 Nlocal: 4000 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8195 ave 8195 max 8195 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 813254 ave 813254 max 813254 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 813254 Ave neighs/atom = 203.314 Neighbor list builds = 6 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 0, Press = 0 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 6.623 | 6.623 | 6.623 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 2000 -6089.1657 -6089.1657 -6245.12 -6245.12 301.70429 301.70429 235286.39 235286.39 50.714994 50.714994 3000 -6100.3068 -6100.3068 -6248.7831 -6248.7831 287.23759 287.23759 235327.05 235327.05 -40.830442 -40.830442 Loop time of 15.8797 on 1 procs for 1000 steps with 4000 atoms Performance: 5.441 ns/day, 4.411 hours/ns, 62.973 timesteps/s 99.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 | 15.543 | 15.543 | 15.543 | 0.0 | 97.88 Neigh | 0.12777 | 0.12777 | 0.12777 | 0.0 | 0.80 Comm | 0.048633 | 0.048633 | 0.048633 | 0.0 | 0.31 Output | 2.9802e-05 | 2.9802e-05 | 2.9802e-05 | 0.0 | 0.00 Modify | 0.14248 | 0.14248 | 0.14248 | 0.0 | 0.90 Other | | 0.0179 | | | 0.11 Nlocal: 4000 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8195 ave 8195 max 8195 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 814498 ave 814498 max 814498 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 814498 Ave neighs/atom = 203.625 Neighbor list builds = 6 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 0, Press = 0 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 6.623 | 6.623 | 6.623 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 3000 -6100.3068 -6100.3068 -6248.7831 -6248.7831 287.23759 287.23759 235327.05 235327.05 -40.830442 -40.830442 4000 -6100.528 -6100.528 -6251.8238 -6251.8238 292.69213 292.69213 235131.57 235131.57 23.126506 23.126506 Loop time of 15.7163 on 1 procs for 1000 steps with 4000 atoms Performance: 5.497 ns/day, 4.366 hours/ns, 63.628 timesteps/s 99.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 | 15.373 | 15.373 | 15.373 | 0.0 | 97.81 Neigh | 0.13844 | 0.13844 | 0.13844 | 0.0 | 0.88 Comm | 0.047204 | 0.047204 | 0.047204 | 0.0 | 0.30 Output | 3.9816e-05 | 3.9816e-05 | 3.9816e-05 | 0.0 | 0.00 Modify | 0.14044 | 0.14044 | 0.14044 | 0.0 | 0.89 Other | | 0.0177 | | | 0.11 Nlocal: 4000 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8195 ave 8195 max 8195 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 813956 ave 813956 max 813956 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 813956 Ave neighs/atom = 203.489 Neighbor list builds = 6 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 0, Press = 0 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 6.623 | 6.623 | 6.623 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 4000 -6100.528 -6100.528 -6251.8238 -6251.8238 292.69213 292.69213 235131.57 235131.57 23.126506 23.126506 5000 -6088.2975 -6088.2975 -6241.1188 -6241.1188 295.64326 295.64326 235587.38 235587.38 -54.75742 -54.75742 Loop time of 15.6675 on 1 procs for 1000 steps with 4000 atoms Performance: 5.515 ns/day, 4.352 hours/ns, 63.826 timesteps/s 100.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 | 15.325 | 15.325 | 15.325 | 0.0 | 97.81 Neigh | 0.13577 | 0.13577 | 0.13577 | 0.0 | 0.87 Comm | 0.047066 | 0.047066 | 0.047066 | 0.0 | 0.30 Output | 2.5988e-05 | 2.5988e-05 | 2.5988e-05 | 0.0 | 0.00 Modify | 0.14227 | 0.14227 | 0.14227 | 0.0 | 0.91 Other | | 0.01765 | | | 0.11 Nlocal: 4000 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8195 ave 8195 max 8195 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 814556 ave 814556 max 814556 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 814556 Ave neighs/atom = 203.639 Neighbor list builds = 6 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 294.454023161355, Press = 75.0162899678082 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 6.623 | 6.623 | 6.623 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 5000 -6088.2975 -6088.2975 -6241.1188 -6241.1188 295.64326 295.64326 235587.38 235587.38 -54.75742 -54.75742 6000 -6100.3253 -6100.3253 -6250.0041 -6250.0041 289.56366 289.56366 234895.57 234895.57 156.42372 156.42372 Loop time of 16.0257 on 1 procs for 1000 steps with 4000 atoms Performance: 5.391 ns/day, 4.452 hours/ns, 62.400 timesteps/s 100.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 | 15.681 | 15.681 | 15.681 | 0.0 | 97.85 Neigh | 0.13214 | 0.13214 | 0.13214 | 0.0 | 0.82 Comm | 0.04767 | 0.04767 | 0.04767 | 0.0 | 0.30 Output | 2.5988e-05 | 2.5988e-05 | 2.5988e-05 | 0.0 | 0.00 Modify | 0.14659 | 0.14659 | 0.14659 | 0.0 | 0.91 Other | | 0.018 | | | 0.11 Nlocal: 4000 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8195 ave 8195 max 8195 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 812692 ave 812692 max 812692 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 812692 Ave neighs/atom = 203.173 Neighbor list builds = 6 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 293.353831083939, Press = -6.6498321093993 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 6.623 | 6.623 | 6.623 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 6000 -6100.3253 -6100.3253 -6250.0041 -6250.0041 289.56366 289.56366 234895.57 234895.57 156.42372 156.42372 7000 -6091.5013 -6091.5013 -6242.245 -6242.245 291.62381 291.62381 235722.37 235722.37 -140.33077 -140.33077 Loop time of 15.5975 on 1 procs for 1000 steps with 4000 atoms Performance: 5.539 ns/day, 4.333 hours/ns, 64.113 timesteps/s 100.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 | 15.268 | 15.268 | 15.268 | 0.0 | 97.89 Neigh | 0.12276 | 0.12276 | 0.12276 | 0.0 | 0.79 Comm | 0.045898 | 0.045898 | 0.045898 | 0.0 | 0.29 Output | 3.9101e-05 | 3.9101e-05 | 3.9101e-05 | 0.0 | 0.00 Modify | 0.14339 | 0.14339 | 0.14339 | 0.0 | 0.92 Other | | 0.01733 | | | 0.11 Nlocal: 4000 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8195 ave 8195 max 8195 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 815530 ave 815530 max 815530 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 815530 Ave neighs/atom = 203.882 Neighbor list builds = 6 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 293.141390356073, Press = 6.00761361125712 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 6.623 | 6.623 | 6.623 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 7000 -6091.5013 -6091.5013 -6242.245 -6242.245 291.62381 291.62381 235722.37 235722.37 -140.33077 -140.33077 8000 -6094.9019 -6094.9019 -6246.3007 -6246.3007 292.89119 292.89119 235145.11 235145.11 82.64261 82.64261 Loop time of 15.5272 on 1 procs for 1000 steps with 4000 atoms Performance: 5.564 ns/day, 4.313 hours/ns, 64.403 timesteps/s 100.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 | 15.186 | 15.186 | 15.186 | 0.0 | 97.80 Neigh | 0.13473 | 0.13473 | 0.13473 | 0.0 | 0.87 Comm | 0.046235 | 0.046235 | 0.046235 | 0.0 | 0.30 Output | 3.0041e-05 | 3.0041e-05 | 3.0041e-05 | 0.0 | 0.00 Modify | 0.14287 | 0.14287 | 0.14287 | 0.0 | 0.92 Other | | 0.01715 | | | 0.11 Nlocal: 4000 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8195 ave 8195 max 8195 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 813232 ave 813232 max 813232 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 813232 Ave neighs/atom = 203.308 Neighbor list builds = 6 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 293.277211530386, Press = -2.12158248748636 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 6.623 | 6.623 | 6.623 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 8000 -6094.9019 -6094.9019 -6246.3007 -6246.3007 292.89119 292.89119 235145.11 235145.11 82.64261 82.64261 9000 -6093.3625 -6093.3625 -6243.9601 -6243.9601 291.34119 291.34119 235462.19 235462.19 -52.445496 -52.445496 Loop time of 15.2268 on 1 procs for 1000 steps with 4000 atoms Performance: 5.674 ns/day, 4.230 hours/ns, 65.674 timesteps/s 100.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.889 | 14.889 | 14.889 | 0.0 | 97.78 Neigh | 0.1356 | 0.1356 | 0.1356 | 0.0 | 0.89 Comm | 0.045621 | 0.045621 | 0.045621 | 0.0 | 0.30 Output | 2.9087e-05 | 2.9087e-05 | 2.9087e-05 | 0.0 | 0.00 Modify | 0.14004 | 0.14004 | 0.14004 | 0.0 | 0.92 Other | | 0.01686 | | | 0.11 Nlocal: 4000 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8195 ave 8195 max 8195 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 814270 ave 814270 max 814270 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 814270 Ave neighs/atom = 203.567 Neighbor list builds = 6 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 293.180128419822, Press = 1.0587135065666 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 6.623 | 6.623 | 6.623 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 9000 -6093.3625 -6093.3625 -6243.9601 -6243.9601 291.34119 291.34119 235462.19 235462.19 -52.445496 -52.445496 10000 -6097.9844 -6097.9844 -6247.6691 -6247.6691 289.57533 289.57533 235145.9 235145.9 57.322192 57.322192 Loop time of 15.2372 on 1 procs for 1000 steps with 4000 atoms Performance: 5.670 ns/day, 4.233 hours/ns, 65.629 timesteps/s 100.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.896 | 14.896 | 14.896 | 0.0 | 97.76 Neigh | 0.13848 | 0.13848 | 0.13848 | 0.0 | 0.91 Comm | 0.045503 | 0.045503 | 0.045503 | 0.0 | 0.30 Output | 2.9087e-05 | 2.9087e-05 | 2.9087e-05 | 0.0 | 0.00 Modify | 0.14037 | 0.14037 | 0.14037 | 0.0 | 0.92 Other | | 0.01703 | | | 0.11 Nlocal: 4000 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8195 ave 8195 max 8195 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 813400 ave 813400 max 813400 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 813400 Ave neighs/atom = 203.35 Neighbor list builds = 6 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 293.379250918381, Press = -0.847830184101612 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 6.623 | 6.623 | 6.623 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 10000 -6097.9844 -6097.9844 -6247.6691 -6247.6691 289.57533 289.57533 235145.9 235145.9 57.322192 57.322192 11000 -6091.7752 -6091.7752 -6245.8431 -6245.8431 298.05474 298.05474 235444.45 235444.45 -42.171736 -42.171736 Loop time of 15.0944 on 1 procs for 1000 steps with 4000 atoms Performance: 5.724 ns/day, 4.193 hours/ns, 66.250 timesteps/s 100.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.763 | 14.763 | 14.763 | 0.0 | 97.80 Neigh | 0.12873 | 0.12873 | 0.12873 | 0.0 | 0.85 Comm | 0.045786 | 0.045786 | 0.045786 | 0.0 | 0.30 Output | 2.9802e-05 | 2.9802e-05 | 2.9802e-05 | 0.0 | 0.00 Modify | 0.14019 | 0.14019 | 0.14019 | 0.0 | 0.93 Other | | 0.01713 | | | 0.11 Nlocal: 4000 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8195 ave 8195 max 8195 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 814922 ave 814922 max 814922 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 814922 Ave neighs/atom = 203.731 Neighbor list builds = 6 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 293.567378819486, Press = 1.0511988757491 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 6.623 | 6.623 | 6.623 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 11000 -6091.7752 -6091.7752 -6245.8431 -6245.8431 298.05474 298.05474 235444.45 235444.45 -42.171736 -42.171736 12000 -6095.1553 -6095.1553 -6246.4525 -6246.4525 292.69479 292.69479 235165.84 235165.84 68.915903 68.915903 Loop time of 17.435 on 1 procs for 1000 steps with 4000 atoms Performance: 4.956 ns/day, 4.843 hours/ns, 57.356 timesteps/s 91.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 | 17.063 | 17.063 | 17.063 | 0.0 | 97.87 Neigh | 0.13859 | 0.13859 | 0.13859 | 0.0 | 0.79 Comm | 0.047768 | 0.047768 | 0.047768 | 0.0 | 0.27 Output | 3.1948e-05 | 3.1948e-05 | 3.1948e-05 | 0.0 | 0.00 Modify | 0.16737 | 0.16737 | 0.16737 | 0.0 | 0.96 Other | | 0.01802 | | | 0.10 Nlocal: 4000 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8195 ave 8195 max 8195 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 813614 ave 813614 max 813614 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 813614 Ave neighs/atom = 203.404 Neighbor list builds = 5 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 293.525701002814, Press = -1.08266955311725 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 6.623 | 6.623 | 6.623 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 12000 -6095.1553 -6095.1553 -6246.4525 -6246.4525 292.69479 292.69479 235165.84 235165.84 68.915903 68.915903 13000 -6097.8322 -6097.8322 -6248.1945 -6248.1945 290.8861 290.8861 235491.42 235491.42 -110.02812 -110.02812 Loop time of 21.8663 on 1 procs for 1000 steps with 4000 atoms Performance: 3.951 ns/day, 6.074 hours/ns, 45.732 timesteps/s 79.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 | 21.39 | 21.39 | 21.39 | 0.0 | 97.82 Neigh | 0.1748 | 0.1748 | 0.1748 | 0.0 | 0.80 Comm | 0.08264 | 0.08264 | 0.08264 | 0.0 | 0.38 Output | 3.0041e-05 | 3.0041e-05 | 3.0041e-05 | 0.0 | 0.00 Modify | 0.17881 | 0.17881 | 0.17881 | 0.0 | 0.82 Other | | 0.03999 | | | 0.18 Nlocal: 4000 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8195 ave 8195 max 8195 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 815794 ave 815794 max 815794 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 815794 Ave neighs/atom = 203.948 Neighbor list builds = 6 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 293.447199518989, Press = 1.63866383444413 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 6.623 | 6.623 | 6.623 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 13000 -6097.8322 -6097.8322 -6248.1945 -6248.1945 290.8861 290.8861 235491.42 235491.42 -110.02812 -110.02812 14000 -6092.0678 -6092.0678 -6247.0458 -6247.0458 299.81538 299.81538 235120.07 235120.07 96.085528 96.085528 Loop time of 20.4464 on 1 procs for 1000 steps with 4000 atoms Performance: 4.226 ns/day, 5.680 hours/ns, 48.908 timesteps/s 84.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 | 20.031 | 20.031 | 20.031 | 0.0 | 97.97 Neigh | 0.14294 | 0.14294 | 0.14294 | 0.0 | 0.70 Comm | 0.05337 | 0.05337 | 0.05337 | 0.0 | 0.26 Output | 2.5034e-05 | 2.5034e-05 | 2.5034e-05 | 0.0 | 0.00 Modify | 0.199 | 0.199 | 0.199 | 0.0 | 0.97 Other | | 0.02012 | | | 0.10 Nlocal: 4000 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8195 ave 8195 max 8195 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 812774 ave 812774 max 812774 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 812774 Ave neighs/atom = 203.194 Neighbor list builds = 6 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 293.392233699345, Press = -0.742551742919957 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 6.623 | 6.623 | 6.623 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 14000 -6092.0678 -6092.0678 -6247.0458 -6247.0458 299.81538 299.81538 235120.07 235120.07 96.085528 96.085528 15000 -6094.5399 -6094.5399 -6244.3062 -6244.3062 289.73322 289.73322 235493.11 235493.11 -67.862413 -67.862413 Loop time of 18.3402 on 1 procs for 1000 steps with 4000 atoms Performance: 4.711 ns/day, 5.095 hours/ns, 54.525 timesteps/s 94.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 | 17.951 | 17.951 | 17.951 | 0.0 | 97.88 Neigh | 0.1342 | 0.1342 | 0.1342 | 0.0 | 0.73 Comm | 0.052332 | 0.052332 | 0.052332 | 0.0 | 0.29 Output | 3.1948e-05 | 3.1948e-05 | 3.1948e-05 | 0.0 | 0.00 Modify | 0.1825 | 0.1825 | 0.1825 | 0.0 | 1.00 Other | | 0.0199 | | | 0.11 Nlocal: 4000 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8195 ave 8195 max 8195 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 814732 ave 814732 max 814732 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 814732 Ave neighs/atom = 203.683 Neighbor list builds = 6 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 293.452469355403, Press = 0.709342673868427 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 6.623 | 6.623 | 6.623 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 15000 -6094.5399 -6094.5399 -6244.3062 -6244.3062 289.73322 289.73322 235493.11 235493.11 -67.862413 -67.862413 16000 -6093.5245 -6093.5245 -6246.0439 -6246.0439 295.05909 295.05909 235283.19 235283.19 27.823402 27.823402 Loop time of 17.3971 on 1 procs for 1000 steps with 4000 atoms Performance: 4.966 ns/day, 4.833 hours/ns, 57.481 timesteps/s 99.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 | 17.022 | 17.022 | 17.022 | 0.0 | 97.85 Neigh | 0.14315 | 0.14315 | 0.14315 | 0.0 | 0.82 Comm | 0.052232 | 0.052232 | 0.052232 | 0.0 | 0.30 Output | 3.0994e-05 | 3.0994e-05 | 3.0994e-05 | 0.0 | 0.00 Modify | 0.15953 | 0.15953 | 0.15953 | 0.0 | 0.92 Other | | 0.01989 | | | 0.11 Nlocal: 4000 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8195 ave 8195 max 8195 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 813742 ave 813742 max 813742 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 813742 Ave neighs/atom = 203.435 Neighbor list builds = 6 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 293.485779282185, Press = -0.0975378747753312 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 6.623 | 6.623 | 6.623 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 16000 -6093.5245 -6093.5245 -6246.0439 -6246.0439 295.05909 295.05909 235283.19 235283.19 27.823402 27.823402 17000 -6095.1407 -6095.1407 -6247.7336 -6247.7336 295.20129 295.20129 235387.58 235387.58 -42.695359 -42.695359 Loop time of 17.3559 on 1 procs for 1000 steps with 4000 atoms Performance: 4.978 ns/day, 4.821 hours/ns, 57.617 timesteps/s 100.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 | 17 | 17 | 17 | 0.0 | 97.95 Neigh | 0.1241 | 0.1241 | 0.1241 | 0.0 | 0.72 Comm | 0.051923 | 0.051923 | 0.051923 | 0.0 | 0.30 Output | 3.0041e-05 | 3.0041e-05 | 3.0041e-05 | 0.0 | 0.00 Modify | 0.15965 | 0.15965 | 0.15965 | 0.0 | 0.92 Other | | 0.01994 | | | 0.11 Nlocal: 4000 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8195 ave 8195 max 8195 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 814284 ave 814284 max 814284 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 814284 Ave neighs/atom = 203.571 Neighbor list builds = 5 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 293.465076556483, Press = 0.246292319964638 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 6.623 | 6.623 | 6.623 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 17000 -6095.1407 -6095.1407 -6247.7336 -6247.7336 295.20129 295.20129 235387.58 235387.58 -42.695359 -42.695359 18000 -6096.8614 -6096.8614 -6246.8903 -6246.8903 290.24102 290.24102 235136.94 235136.94 73.384598 73.384598 Loop time of 17.426 on 1 procs for 1000 steps with 4000 atoms Performance: 4.958 ns/day, 4.841 hours/ns, 57.386 timesteps/s 99.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 | 17.071 | 17.071 | 17.071 | 0.0 | 97.96 Neigh | 0.12146 | 0.12146 | 0.12146 | 0.0 | 0.70 Comm | 0.052125 | 0.052125 | 0.052125 | 0.0 | 0.30 Output | 2.5988e-05 | 2.5988e-05 | 2.5988e-05 | 0.0 | 0.00 Modify | 0.16154 | 0.16154 | 0.16154 | 0.0 | 0.93 Other | | 0.02003 | | | 0.11 Nlocal: 4000 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8195 ave 8195 max 8195 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 813576 ave 813576 max 813576 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 813576 Ave neighs/atom = 203.394 Neighbor list builds = 5 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 293.43351530685, Press = 0.225717501310757 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 6.623 | 6.623 | 6.623 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 18000 -6096.8614 -6096.8614 -6246.8903 -6246.8903 290.24102 290.24102 235136.94 235136.94 73.384598 73.384598 19000 -6093.7798 -6093.7798 -6246.4832 -6246.4832 295.41512 295.41512 235521.82 235521.82 -98.888406 -98.888406 Loop time of 17.4158 on 1 procs for 1000 steps with 4000 atoms Performance: 4.961 ns/day, 4.838 hours/ns, 57.419 timesteps/s 99.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 | 17.037 | 17.037 | 17.037 | 0.0 | 97.82 Neigh | 0.14841 | 0.14841 | 0.14841 | 0.0 | 0.85 Comm | 0.052024 | 0.052024 | 0.052024 | 0.0 | 0.30 Output | 3.1948e-05 | 3.1948e-05 | 3.1948e-05 | 0.0 | 0.00 Modify | 0.15904 | 0.15904 | 0.15904 | 0.0 | 0.91 Other | | 0.01967 | | | 0.11 Nlocal: 4000 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8195 ave 8195 max 8195 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 813740 ave 813740 max 813740 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 813740 Ave neighs/atom = 203.435 Neighbor list builds = 6 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 293.356544757187, Press = -0.13227068455531 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 6.623 | 6.623 | 6.623 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 19000 -6093.7798 -6093.7798 -6246.4832 -6246.4832 295.41512 295.41512 235521.82 235521.82 -98.888406 -98.888406 20000 -6094.0829 -6094.0829 -6244.6595 -6244.6595 291.30047 291.30047 235053.2 235053.2 134.82041 134.82041 Loop time of 17.4776 on 1 procs for 1000 steps with 4000 atoms Performance: 4.943 ns/day, 4.855 hours/ns, 57.216 timesteps/s 99.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 | 17.103 | 17.103 | 17.103 | 0.0 | 97.86 Neigh | 0.14385 | 0.14385 | 0.14385 | 0.0 | 0.82 Comm | 0.05214 | 0.05214 | 0.05214 | 0.0 | 0.30 Output | 3.314e-05 | 3.314e-05 | 3.314e-05 | 0.0 | 0.00 Modify | 0.15895 | 0.15895 | 0.15895 | 0.0 | 0.91 Other | | 0.01975 | | | 0.11 Nlocal: 4000 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8195 ave 8195 max 8195 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 814726 ave 814726 max 814726 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 814726 Ave neighs/atom = 203.681 Neighbor list builds = 6 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 293.350103537609, Press = 0.480120305307059 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 6.623 | 6.623 | 6.623 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 20000 -6094.0829 -6094.0829 -6244.6595 -6244.6595 291.30047 291.30047 235053.2 235053.2 134.82041 134.82041 21000 -6097.6539 -6097.6539 -6248.9295 -6248.9295 292.65289 292.65289 235567.58 235567.58 -145.30596 -145.30596 Loop time of 18.83 on 1 procs for 1000 steps with 4000 atoms Performance: 4.588 ns/day, 5.231 hours/ns, 53.107 timesteps/s 91.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 | 18.446 | 18.446 | 18.446 | 0.0 | 97.96 Neigh | 0.15295 | 0.15295 | 0.15295 | 0.0 | 0.81 Comm | 0.052244 | 0.052244 | 0.052244 | 0.0 | 0.28 Output | 3.314e-05 | 3.314e-05 | 3.314e-05 | 0.0 | 0.00 Modify | 0.15858 | 0.15858 | 0.15858 | 0.0 | 0.84 Other | | 0.02001 | | | 0.11 Nlocal: 4000 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8195 ave 8195 max 8195 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 813018 ave 813018 max 813018 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 813018 Ave neighs/atom = 203.255 Neighbor list builds = 6 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 293.354580561096, Press = -0.414775570934004 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 6.623 | 6.623 | 6.623 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 21000 -6097.6539 -6097.6539 -6248.9295 -6248.9295 292.65289 292.65289 235567.58 235567.58 -145.30596 -145.30596 22000 -6092.8709 -6092.8709 -6243.1752 -6243.1752 290.77394 290.77394 235077.3 235077.3 150.34696 150.34696 Loop time of 28.3282 on 1 procs for 1000 steps with 4000 atoms Performance: 3.050 ns/day, 7.869 hours/ns, 35.300 timesteps/s 61.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 | 27.664 | 27.664 | 27.664 | 0.0 | 97.65 Neigh | 0.26499 | 0.26499 | 0.26499 | 0.0 | 0.94 Comm | 0.080801 | 0.080801 | 0.080801 | 0.0 | 0.29 Output | 4.5061e-05 | 4.5061e-05 | 4.5061e-05 | 0.0 | 0.00 Modify | 0.29836 | 0.29836 | 0.29836 | 0.0 | 1.05 Other | | 0.02018 | | | 0.07 Nlocal: 4000 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8195 ave 8195 max 8195 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 814898 ave 814898 max 814898 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 814898 Ave neighs/atom = 203.725 Neighbor list builds = 6 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 293.30042833005, Press = 0.542620590873831 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 6.623 | 6.623 | 6.623 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 22000 -6092.8709 -6092.8709 -6243.1752 -6243.1752 290.77394 290.77394 235077.3 235077.3 150.34696 150.34696 23000 -6089.9297 -6089.9297 -6244.0698 -6244.0698 298.19445 298.19445 235589.56 235589.56 -92.810626 -92.810626 Loop time of 30.4579 on 1 procs for 1000 steps with 4000 atoms Performance: 2.837 ns/day, 8.461 hours/ns, 32.832 timesteps/s 57.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 | 29.769 | 29.769 | 29.769 | 0.0 | 97.74 Neigh | 0.2295 | 0.2295 | 0.2295 | 0.0 | 0.75 Comm | 0.11311 | 0.11311 | 0.11311 | 0.0 | 0.37 Output | 3.2902e-05 | 3.2902e-05 | 3.2902e-05 | 0.0 | 0.00 Modify | 0.32542 | 0.32542 | 0.32542 | 0.0 | 1.07 Other | | 0.02053 | | | 0.07 Nlocal: 4000 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8195 ave 8195 max 8195 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 813046 ave 813046 max 813046 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 813046 Ave neighs/atom = 203.262 Neighbor list builds = 6 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 293.289441895362, Press = -0.26567506831219 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 6.623 | 6.623 | 6.623 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 23000 -6089.9297 -6089.9297 -6244.0698 -6244.0698 298.19445 298.19445 235589.56 235589.56 -92.810626 -92.810626 24000 -6093.4443 -6093.4443 -6246.3673 -6246.3673 295.83988 295.83988 235152.51 235152.51 89.335836 89.335836 Loop time of 26.1656 on 1 procs for 1000 steps with 4000 atoms Performance: 3.302 ns/day, 7.268 hours/ns, 38.218 timesteps/s 66.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 | 25.597 | 25.597 | 25.597 | 0.0 | 97.83 Neigh | 0.19765 | 0.19765 | 0.19765 | 0.0 | 0.76 Comm | 0.072677 | 0.072677 | 0.072677 | 0.0 | 0.28 Output | 3.2902e-05 | 3.2902e-05 | 3.2902e-05 | 0.0 | 0.00 Modify | 0.2582 | 0.2582 | 0.2582 | 0.0 | 0.99 Other | | 0.04003 | | | 0.15 Nlocal: 4000 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8195 ave 8195 max 8195 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 815176 ave 815176 max 815176 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 815176 Ave neighs/atom = 203.794 Neighbor list builds = 6 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 293.242559398914, Press = 0.400362306936302 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 6.623 | 6.623 | 6.623 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 24000 -6093.4443 -6093.4443 -6246.3673 -6246.3673 295.83988 295.83988 235152.51 235152.51 89.335836 89.335836 25000 -6095.2685 -6095.2685 -6246.3725 -6246.3725 292.32097 292.32097 235588.83 235588.83 -132.9713 -132.9713 Loop time of 26.6305 on 1 procs for 1000 steps with 4000 atoms Performance: 3.244 ns/day, 7.397 hours/ns, 37.551 timesteps/s 65.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 | 26.017 | 26.017 | 26.017 | 0.0 | 97.70 Neigh | 0.20311 | 0.20311 | 0.20311 | 0.0 | 0.76 Comm | 0.092522 | 0.092522 | 0.092522 | 0.0 | 0.35 Output | 3.2902e-05 | 3.2902e-05 | 3.2902e-05 | 0.0 | 0.00 Modify | 0.29765 | 0.29765 | 0.29765 | 0.0 | 1.12 Other | | 0.01983 | | | 0.07 Nlocal: 4000 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8195 ave 8195 max 8195 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 812692 ave 812692 max 812692 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 812692 Ave neighs/atom = 203.173 Neighbor list builds = 6 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 293.212346112079, Press = -0.590305115266853 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 6.623 | 6.623 | 6.623 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 25000 -6095.2685 -6095.2685 -6246.3725 -6246.3725 292.32097 292.32097 235588.83 235588.83 -132.9713 -132.9713 26000 -6094.6893 -6094.6893 -6246.8664 -6246.8664 294.39697 294.39697 235099.61 235099.61 100.57769 100.57769 Loop time of 26.5701 on 1 procs for 1000 steps with 4000 atoms Performance: 3.252 ns/day, 7.381 hours/ns, 37.636 timesteps/s 65.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 | 25.912 | 25.912 | 25.912 | 0.0 | 97.52 Neigh | 0.22629 | 0.22629 | 0.22629 | 0.0 | 0.85 Comm | 0.13293 | 0.13293 | 0.13293 | 0.0 | 0.50 Output | 2.5988e-05 | 2.5988e-05 | 2.5988e-05 | 0.0 | 0.00 Modify | 0.27821 | 0.27821 | 0.27821 | 0.0 | 1.05 Other | | 0.02033 | | | 0.08 Nlocal: 4000 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8195 ave 8195 max 8195 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 815542 ave 815542 max 815542 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 815542 Ave neighs/atom = 203.886 Neighbor list builds = 6 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 293.219965319692, Press = 0.671982586511759 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 6.623 | 6.623 | 6.623 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 26000 -6094.6893 -6094.6893 -6246.8664 -6246.8664 294.39697 294.39697 235099.61 235099.61 100.57769 100.57769 27000 -6094.21 -6094.21 -6248.1551 -6248.1551 297.81734 297.81734 235461.06 235461.06 -76.702783 -76.702783 Loop time of 24.8915 on 1 procs for 1000 steps with 4000 atoms Performance: 3.471 ns/day, 6.914 hours/ns, 40.174 timesteps/s 69.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 | 24.38 | 24.38 | 24.38 | 0.0 | 97.95 Neigh | 0.18044 | 0.18044 | 0.18044 | 0.0 | 0.72 Comm | 0.072618 | 0.072618 | 0.072618 | 0.0 | 0.29 Output | 3.314e-05 | 3.314e-05 | 3.314e-05 | 0.0 | 0.00 Modify | 0.23826 | 0.23826 | 0.23826 | 0.0 | 0.96 Other | | 0.01979 | | | 0.08 Nlocal: 4000 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8195 ave 8195 max 8195 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 813402 ave 813402 max 813402 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 813402 Ave neighs/atom = 203.351 Neighbor list builds = 6 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 293.197243120994, Press = -0.318016037686276 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 6.623 | 6.623 | 6.623 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 27000 -6094.21 -6094.21 -6248.1551 -6248.1551 297.81734 297.81734 235461.06 235461.06 -76.702783 -76.702783 28000 -6095.481 -6095.481 -6245.7854 -6245.7854 290.77399 290.77399 235210.29 235210.29 48.405297 48.405297 Loop time of 22.7261 on 1 procs for 1000 steps with 4000 atoms Performance: 3.802 ns/day, 6.313 hours/ns, 44.002 timesteps/s 76.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 | 22.27 | 22.27 | 22.27 | 0.0 | 97.99 Neigh | 0.13332 | 0.13332 | 0.13332 | 0.0 | 0.59 Comm | 0.092891 | 0.092891 | 0.092891 | 0.0 | 0.41 Output | 2.9087e-05 | 2.9087e-05 | 2.9087e-05 | 0.0 | 0.00 Modify | 0.20904 | 0.20904 | 0.20904 | 0.0 | 0.92 Other | | 0.02052 | | | 0.09 Nlocal: 4000 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8195 ave 8195 max 8195 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 813864 ave 813864 max 813864 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 813864 Ave neighs/atom = 203.466 Neighbor list builds = 5 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 293.108264335758, Press = 0.256698086871175 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 6.623 | 6.623 | 6.623 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 28000 -6095.481 -6095.481 -6245.7854 -6245.7854 290.77399 290.77399 235210.29 235210.29 48.405297 48.405297 29000 -6095.3572 -6095.3572 -6244.5791 -6244.5791 288.67983 288.67983 235403.24 235403.24 -31.027171 -31.027171 Loop time of 25.1115 on 1 procs for 1000 steps with 4000 atoms Performance: 3.441 ns/day, 6.975 hours/ns, 39.822 timesteps/s 69.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 | 22.919 | 22.919 | 22.919 | 0.0 | 91.27 Neigh | 0.25135 | 0.25135 | 0.25135 | 0.0 | 1.00 Comm | 0.076241 | 0.076241 | 0.076241 | 0.0 | 0.30 Output | 0.017597 | 0.017597 | 0.017597 | 0.0 | 0.07 Modify | 1.8162 | 1.8162 | 1.8162 | 0.0 | 7.23 Other | | 0.03092 | | | 0.12 Nlocal: 4000 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8195 ave 8195 max 8195 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 813778 ave 813778 max 813778 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 813778 Ave neighs/atom = 203.445 Neighbor list builds = 6 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" 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 235327.867125249 print "LAMMPS calculation completed" LAMMPS calculation completed quit 0