# periodic boundary conditions along all three dimensions boundary p p p # Set neighbor skin variable neigh_skin equal 2.0*${_u_distance} variable neigh_skin equal 2.0*1 neighbor ${neigh_skin} bin neighbor 2 bin # create a supercell with cubic lattice (fcc, bcc, sc, or diamond) # using 10*10*10 conventional (orthogonal) unit cells variable latticeconst_converted equal 3.190829664468765*${_u_distance} variable latticeconst_converted equal 3.190829664468765*1 lattice bcc ${latticeconst_converted} lattice bcc 3.19082966446876 Lattice spacing in x,y,z = 3.1908297 3.1908297 3.1908297 region simbox block 0 10 0 10 0 10 units lattice create_box 1 simbox Created orthogonal box = (0 0 0) to (31.908297 31.908297 31.908297) 1 by 1 by 1 MPI processor grid create_atoms 1 box Created 2000 atoms using lattice units in orthogonal box = (0 0 0) to (31.908297 31.908297 31.908297) create_atoms CPU = 0.000 seconds variable mass_converted equal 183.84*${_u_mass} variable mass_converted equal 183.84*1 kim_interactions W WARNING: 'kim_' has been renamed to 'kim '. Please update your input. kim interactions W #=== BEGIN kim interactions ================================== pair_style kim MEAM_LAMMPS_Lenosky_2017_W__MO_999198119251_002 pair_coeff * * W #=== END kim interactions ==================================== mass 1 ${mass_converted} mass 1 183.84 # initial volume variable v equal vol # assign formula variable V0 equal ${v} # evaluate initial value variable V0 equal 32487.0938338165 variable V0_metal equal ${V0}/(${_u_distance}*${_u_distance}*${_u_distance}) variable V0_metal equal 32487.0938338165/(${_u_distance}*${_u_distance}*${_u_distance}) variable V0_metal equal 32487.0938338165/(1*${_u_distance}*${_u_distance}) variable V0_metal equal 32487.0938338165/(1*1*${_u_distance}) variable V0_metal equal 32487.0938338165/(1*1*1) variable V0_metal_times1000 equal ${V0_metal}*1000 variable V0_metal_times1000 equal 32487.0938338165*1000 print "Initial system volume: ${V0_metal} Angstroms^3" Initial system volume: 32487.0938338165 Angstroms^3 # set the time step to 0.001 picoseconds variable timestep_converted equal 0.001*${_u_time} variable timestep_converted equal 0.001*1 timestep ${timestep_converted} timestep 0.001 variable temp_converted equal 313.15*${_u_temperature} variable temp_converted equal 313.15*1 variable Tdamp_converted equal 0.1*${_u_time} variable Tdamp_converted equal 0.1*1 variable press_converted equal 0.0*${_u_pressure} variable press_converted equal 0.0*1 variable Pdamp_converted equal 1*${_u_time} variable Pdamp_converted equal 1*1 # create initial velocities consistent with the chosen temperature velocity all create ${temp_converted} 17 mom yes rot yes velocity all create 313.15 17 mom yes rot yes # set NPT ensemble for all atoms fix ensemble all npt temp ${temp_converted} ${temp_converted} ${Tdamp_converted} iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 313.15 ${temp_converted} ${Tdamp_converted} iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 313.15 313.15 ${Tdamp_converted} iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 313.15 313.15 0.1 iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 313.15 313.15 0.1 iso 0 ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 313.15 313.15 0.1 iso 0 0 ${Pdamp_converted} fix ensemble all npt temp 313.15 313.15 0.1 iso 0 0 1 # compute the time averages of pressure, temperature, and volume, respectively # ignore the first 5000 timesteps variable etotal_metal equal etotal/${_u_energy} variable etotal_metal equal etotal/1 variable pe_metal equal pe/${_u_energy} variable pe_metal equal pe/1 variable T_metal equal temp/${_u_temperature} variable T_metal equal temp/1 variable V_metal equal vol/(${_u_distance}*${_u_distance}*${_u_distance}) variable V_metal equal vol/(1*${_u_distance}*${_u_distance}) variable V_metal equal vol/(1*1*${_u_distance}) variable V_metal equal vol/(1*1*1) variable P_metal equal press/${_u_pressure} variable P_metal equal press/1 fix avgmyTemp all ave/time 5 20 100 v_T_metal ave running start 5000 fix avgmyPress all ave/time 5 20 100 v_P_metal ave running start 5000 fix avgmyVol all ave/time 5 20 100 v_V_metal ave running start 5000 # extract fix quantities into variables so they can be used in if-else logic later. variable T equal f_avgmyTemp variable P equal f_avgmyPress variable V equal f_avgmyVol # set error bounds for temperature and pressure in original metal units (K and bar) variable T_low equal "313.15 - 0.2" variable T_up equal "313.15 + 0.2" variable P_low equal "0.0 - 0.2" variable P_up equal "0.0 + 0.2" # print to logfile every 1000 timesteps thermo_style custom step etotal v_etotal_metal pe v_pe_metal temp v_T_metal vol v_V_metal press v_P_metal thermo 1000 # Run a simulation for at most 2000*1000 timesteps. At each 1000th time step, check # whether the temperature and pressure have converged. If yes, break. label top variable a loop 2000 run 1000 CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE Your simulation uses code contributions which should be cited: - OpenKIM: https://doi.org/10.1007/s11837-011-0102-6 @Article{tadmor:elliott:2011, author = {E. B. Tadmor and R. S. Elliott and J. P. Sethna and R. E. Miller and C. A. Becker}, title = {The potential of atomistic simulations and the {K}nowledgebase of {I}nteratomic {M}odels}, journal = {{JOM}}, year = 2011, volume = 63, number = 17, pages = {17}, doi = {10.1007/s11837-011-0102-6} } - OpenKIM potential: https://openkim.org/cite/MO_999198119251_002#item-citation CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7.5 ghost atom cutoff = 7.5 binsize = 3.75, bins = 9 9 9 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 0 -17613.702 -17613.702 -17694.618 -17694.618 313.15 313.15 32487.094 32487.094 2660.2826 2660.2826 1000 -17531.605 -17531.605 -17614.013 -17614.013 318.92441 318.92441 32527.302 32527.302 2670.1189 2670.1189 Loop time of 98.5009 on 1 procs for 1000 steps with 2000 atoms Performance: 0.877 ns/day, 27.361 hours/ns, 10.152 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 | 98.42 | 98.42 | 98.42 | 0.0 | 99.92 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.015144 | 0.015144 | 0.015144 | 0.0 | 0.02 Output | 0.00011778 | 0.00011778 | 0.00011778 | 0.0 | 0.00 Modify | 0.059959 | 0.059959 | 0.059959 | 0.0 | 0.06 Other | | 0.005536 | | | 0.01 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 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: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 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 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7.5 ghost atom cutoff = 7.5 binsize = 3.75, bins = 9 9 9 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 1000 -17531.605 -17531.605 -17614.013 -17614.013 318.92441 318.92441 32527.302 32527.302 2670.1189 2670.1189 2000 -17531.516 -17531.516 -17621.685 -17621.685 348.95988 348.95988 32539.514 32539.514 884.47214 884.47214 Loop time of 100.033 on 1 procs for 1000 steps with 2000 atoms Performance: 0.864 ns/day, 27.787 hours/ns, 9.997 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 | 99.953 | 99.953 | 99.953 | 0.0 | 99.92 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.01499 | 0.01499 | 0.01499 | 0.0 | 0.01 Output | 3.2832e-05 | 3.2832e-05 | 3.2832e-05 | 0.0 | 0.00 Modify | 0.059836 | 0.059836 | 0.059836 | 0.0 | 0.06 Other | | 0.005358 | | | 0.01 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 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: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 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 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7.5 ghost atom cutoff = 7.5 binsize = 3.75, bins = 9 9 9 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 2000 -17531.516 -17531.516 -17621.685 -17621.685 348.95988 348.95988 32539.514 32539.514 884.47214 884.47214 3000 -17531.986 -17531.986 -17616.482 -17616.482 327.00939 327.00939 32515.294 32515.294 3373.4529 3373.4529 Loop time of 100.104 on 1 procs for 1000 steps with 2000 atoms Performance: 0.863 ns/day, 27.807 hours/ns, 9.990 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 | 100.02 | 100.02 | 100.02 | 0.0 | 99.92 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.014963 | 0.014963 | 0.014963 | 0.0 | 0.01 Output | 3.4074e-05 | 3.4074e-05 | 3.4074e-05 | 0.0 | 0.00 Modify | 0.059616 | 0.059616 | 0.059616 | 0.0 | 0.06 Other | | 0.005365 | | | 0.01 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 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: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 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 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7.5 ghost atom cutoff = 7.5 binsize = 3.75, bins = 9 9 9 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 3000 -17531.986 -17531.986 -17616.482 -17616.482 327.00939 327.00939 32515.294 32515.294 3373.4529 3373.4529 4000 -17530.964 -17530.964 -17611.311 -17611.311 310.95208 310.95208 32557.352 32557.352 -188.67959 -188.67959 Loop time of 100.081 on 1 procs for 1000 steps with 2000 atoms Performance: 0.863 ns/day, 27.800 hours/ns, 9.992 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 | 100 | 100 | 100 | 0.0 | 99.92 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.015027 | 0.015027 | 0.015027 | 0.0 | 0.02 Output | 3.4175e-05 | 3.4175e-05 | 3.4175e-05 | 0.0 | 0.00 Modify | 0.060062 | 0.060062 | 0.060062 | 0.0 | 0.06 Other | | 0.005353 | | | 0.01 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 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: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 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 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7.5 ghost atom cutoff = 7.5 binsize = 3.75, bins = 9 9 9 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 4000 -17530.964 -17530.964 -17611.311 -17611.311 310.95208 310.95208 32557.352 32557.352 -188.67959 -188.67959 5000 -17533.16 -17533.16 -17617.055 -17617.055 324.67994 324.67994 32557.071 32557.071 -418.37137 -418.37137 Loop time of 100.129 on 1 procs for 1000 steps with 2000 atoms Performance: 0.863 ns/day, 27.814 hours/ns, 9.987 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 | 100.05 | 100.05 | 100.05 | 0.0 | 99.92 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.014983 | 0.014983 | 0.014983 | 0.0 | 0.01 Output | 2.6299e-05 | 2.6299e-05 | 2.6299e-05 | 0.0 | 0.00 Modify | 0.060652 | 0.060652 | 0.060652 | 0.0 | 0.06 Other | | 0.005332 | | | 0.01 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 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: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 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 = 309.257443187525, Press = -136.809404674847 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7.5 ghost atom cutoff = 7.5 binsize = 3.75, bins = 9 9 9 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 5000 -17533.16 -17533.16 -17617.055 -17617.055 324.67994 324.67994 32557.071 32557.071 -418.37137 -418.37137 6000 -17532.845 -17532.845 -17615.043 -17615.043 318.11458 318.11458 32578.702 32578.702 -2501.0911 -2501.0911 Loop time of 100.113 on 1 procs for 1000 steps with 2000 atoms Performance: 0.863 ns/day, 27.809 hours/ns, 9.989 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 | 100.03 | 100.03 | 100.03 | 0.0 | 99.92 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.015 | 0.015 | 0.015 | 0.0 | 0.01 Output | 2.9214e-05 | 2.9214e-05 | 2.9214e-05 | 0.0 | 0.00 Modify | 0.061866 | 0.061866 | 0.061866 | 0.0 | 0.06 Other | | 0.005374 | | | 0.01 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 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: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 312.638482568981, Press = -51.4385908916566 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7.5 ghost atom cutoff = 7.5 binsize = 3.75, bins = 9 9 9 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 6000 -17532.845 -17532.845 -17615.043 -17615.043 318.11458 318.11458 32578.702 32578.702 -2501.0911 -2501.0911 7000 -17531.447 -17531.447 -17607.689 -17607.689 295.06497 295.06497 32536.882 32536.882 1808.77 1808.77 Loop time of 100.109 on 1 procs for 1000 steps with 2000 atoms Performance: 0.863 ns/day, 27.808 hours/ns, 9.989 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 | 100.03 | 100.03 | 100.03 | 0.0 | 99.92 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.015056 | 0.015056 | 0.015056 | 0.0 | 0.02 Output | 2.8684e-05 | 2.8684e-05 | 2.8684e-05 | 0.0 | 0.00 Modify | 0.061758 | 0.061758 | 0.061758 | 0.0 | 0.06 Other | | 0.005365 | | | 0.01 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 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: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 312.674688847653, Press = -27.5743483177687 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7.5 ghost atom cutoff = 7.5 binsize = 3.75, bins = 9 9 9 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 7000 -17531.447 -17531.447 -17607.689 -17607.689 295.06497 295.06497 32536.882 32536.882 1808.77 1808.77 8000 -17533.889 -17533.889 -17613.394 -17613.394 307.69495 307.69495 32569.645 32569.645 -1811.5491 -1811.5491 Loop time of 100.076 on 1 procs for 1000 steps with 2000 atoms Performance: 0.863 ns/day, 27.799 hours/ns, 9.992 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 | 99.994 | 99.994 | 99.994 | 0.0 | 99.92 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.014986 | 0.014986 | 0.014986 | 0.0 | 0.01 Output | 2.7021e-05 | 2.7021e-05 | 2.7021e-05 | 0.0 | 0.00 Modify | 0.061843 | 0.061843 | 0.061843 | 0.0 | 0.06 Other | | 0.005376 | | | 0.01 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 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: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 312.243314744399, Press = -7.49384442539178 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7.5 ghost atom cutoff = 7.5 binsize = 3.75, bins = 9 9 9 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 8000 -17533.889 -17533.889 -17613.394 -17613.394 307.69495 307.69495 32569.645 32569.645 -1811.5491 -1811.5491 9000 -17531.778 -17531.778 -17611.606 -17611.606 308.94363 308.94363 32549.015 32549.015 647.55122 647.55122 Loop time of 100.081 on 1 procs for 1000 steps with 2000 atoms Performance: 0.863 ns/day, 27.800 hours/ns, 9.992 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 | 99.999 | 99.999 | 99.999 | 0.0 | 99.92 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.01504 | 0.01504 | 0.01504 | 0.0 | 0.02 Output | 2.9676e-05 | 2.9676e-05 | 2.9676e-05 | 0.0 | 0.00 Modify | 0.061842 | 0.061842 | 0.061842 | 0.0 | 0.06 Other | | 0.005368 | | | 0.01 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 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: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 312.347211979902, Press = 14.7245670770408 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7.5 ghost atom cutoff = 7.5 binsize = 3.75, bins = 9 9 9 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 9000 -17531.778 -17531.778 -17611.606 -17611.606 308.94363 308.94363 32549.015 32549.015 647.55122 647.55122 10000 -17529.171 -17529.171 -17612.795 -17612.795 323.63172 323.63172 32598.128 32598.128 -3983.223 -3983.223 Loop time of 100.068 on 1 procs for 1000 steps with 2000 atoms Performance: 0.863 ns/day, 27.797 hours/ns, 9.993 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 | 99.986 | 99.986 | 99.986 | 0.0 | 99.92 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.01497 | 0.01497 | 0.01497 | 0.0 | 0.01 Output | 2.9115e-05 | 2.9115e-05 | 2.9115e-05 | 0.0 | 0.00 Modify | 0.061758 | 0.061758 | 0.061758 | 0.0 | 0.06 Other | | 0.005386 | | | 0.01 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 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: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 312.489070203691, Press = -7.62646196753232 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7.5 ghost atom cutoff = 7.5 binsize = 3.75, bins = 9 9 9 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 10000 -17529.171 -17529.171 -17612.795 -17612.795 323.63172 323.63172 32598.128 32598.128 -3983.223 -3983.223 11000 -17532.204 -17532.204 -17614.262 -17614.262 317.57378 317.57378 32501.903 32501.903 4412.8234 4412.8234 Loop time of 100.109 on 1 procs for 1000 steps with 2000 atoms Performance: 0.863 ns/day, 27.808 hours/ns, 9.989 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 | 100.03 | 100.03 | 100.03 | 0.0 | 99.92 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.015015 | 0.015015 | 0.015015 | 0.0 | 0.01 Output | 2.697e-05 | 2.697e-05 | 2.697e-05 | 0.0 | 0.00 Modify | 0.061714 | 0.061714 | 0.061714 | 0.0 | 0.06 Other | | 0.005378 | | | 0.01 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 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: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 312.543603104067, Press = 7.58609669347483 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7.5 ghost atom cutoff = 7.5 binsize = 3.75, bins = 9 9 9 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 11000 -17532.204 -17532.204 -17614.262 -17614.262 317.57378 317.57378 32501.903 32501.903 4412.8234 4412.8234 12000 -17530.772 -17530.772 -17616.786 -17616.786 332.88423 332.88423 32555.465 32555.465 -181.3268 -181.3268 Loop time of 100.06 on 1 procs for 1000 steps with 2000 atoms Performance: 0.863 ns/day, 27.795 hours/ns, 9.994 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 | 99.978 | 99.978 | 99.978 | 0.0 | 99.92 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.014996 | 0.014996 | 0.014996 | 0.0 | 0.01 Output | 3.5216e-05 | 3.5216e-05 | 3.5216e-05 | 0.0 | 0.00 Modify | 0.061961 | 0.061961 | 0.061961 | 0.0 | 0.06 Other | | 0.005381 | | | 0.01 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 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: 224002 ave 224002 max 224002 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224002 Ave neighs/atom = 112.001 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 312.769427055901, Press = 1.2045146368954 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7.5 ghost atom cutoff = 7.5 binsize = 3.75, bins = 9 9 9 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 12000 -17530.772 -17530.772 -17616.786 -17616.786 332.88423 332.88423 32555.465 32555.465 -181.3268 -181.3268 13000 -17528.932 -17528.932 -17616.154 -17616.154 337.55774 337.55774 32541.395 32541.395 1198.1134 1198.1134 Loop time of 100.105 on 1 procs for 1000 steps with 2000 atoms Performance: 0.863 ns/day, 27.807 hours/ns, 9.990 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 | 100.02 | 100.02 | 100.02 | 0.0 | 99.92 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.014977 | 0.014977 | 0.014977 | 0.0 | 0.01 Output | 2.8132e-05 | 2.8132e-05 | 2.8132e-05 | 0.0 | 0.00 Modify | 0.061751 | 0.061751 | 0.061751 | 0.0 | 0.06 Other | | 0.005364 | | | 0.01 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 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: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 312.878917178364, Press = -5.87779849521185 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7.5 ghost atom cutoff = 7.5 binsize = 3.75, bins = 9 9 9 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 13000 -17528.932 -17528.932 -17616.154 -17616.154 337.55774 337.55774 32541.395 32541.395 1198.1134 1198.1134 14000 -17530.825 -17530.825 -17618.721 -17618.721 340.16333 340.16333 32510.817 32510.817 3699.3176 3699.3176 Loop time of 100.096 on 1 procs for 1000 steps with 2000 atoms Performance: 0.863 ns/day, 27.805 hours/ns, 9.990 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 | 100.01 | 100.01 | 100.01 | 0.0 | 99.92 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.015005 | 0.015005 | 0.015005 | 0.0 | 0.01 Output | 5.0856e-05 | 5.0856e-05 | 5.0856e-05 | 0.0 | 0.00 Modify | 0.061696 | 0.061696 | 0.061696 | 0.0 | 0.06 Other | | 0.005396 | | | 0.01 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 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: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 313.045429912709, Press = 0.473742344868467 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7.5 ghost atom cutoff = 7.5 binsize = 3.75, bins = 9 9 9 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 14000 -17530.825 -17530.825 -17618.721 -17618.721 340.16333 340.16333 32510.817 32510.817 3699.3176 3699.3176 15000 -17528.864 -17528.864 -17611.895 -17611.895 321.33877 321.33877 32550.082 32550.082 723.99477 723.99477 Loop time of 100.065 on 1 procs for 1000 steps with 2000 atoms Performance: 0.863 ns/day, 27.796 hours/ns, 9.993 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 | 99.983 | 99.983 | 99.983 | 0.0 | 99.92 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.015006 | 0.015006 | 0.015006 | 0.0 | 0.01 Output | 2.9505e-05 | 2.9505e-05 | 2.9505e-05 | 0.0 | 0.00 Modify | 0.06182 | 0.06182 | 0.06182 | 0.0 | 0.06 Other | | 0.005368 | | | 0.01 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 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: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 313.040648441087, Press = 3.73625707521291 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7.5 ghost atom cutoff = 7.5 binsize = 3.75, bins = 9 9 9 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 15000 -17528.864 -17528.864 -17611.895 -17611.895 321.33877 321.33877 32550.082 32550.082 723.99477 723.99477 16000 -17534.786 -17534.786 -17612.675 -17612.675 301.43991 301.43991 32557.168 32557.168 -593.46603 -593.46603 Loop time of 100.094 on 1 procs for 1000 steps with 2000 atoms Performance: 0.863 ns/day, 27.804 hours/ns, 9.991 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 | 100.01 | 100.01 | 100.01 | 0.0 | 99.92 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.014939 | 0.014939 | 0.014939 | 0.0 | 0.01 Output | 3.4555e-05 | 3.4555e-05 | 3.4555e-05 | 0.0 | 0.00 Modify | 0.061786 | 0.061786 | 0.061786 | 0.0 | 0.06 Other | | 0.005364 | | | 0.01 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 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: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 312.830776554284, Press = -2.1477308549846 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7.5 ghost atom cutoff = 7.5 binsize = 3.75, bins = 9 9 9 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 16000 -17534.786 -17534.786 -17612.675 -17612.675 301.43991 301.43991 32557.168 32557.168 -593.46603 -593.46603 17000 -17531.801 -17531.801 -17609.866 -17609.866 302.12097 302.12097 32534.918 32534.918 1928.4555 1928.4555 Loop time of 100.069 on 1 procs for 1000 steps with 2000 atoms Performance: 0.863 ns/day, 27.797 hours/ns, 9.993 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 | 99.987 | 99.987 | 99.987 | 0.0 | 99.92 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.014996 | 0.014996 | 0.014996 | 0.0 | 0.01 Output | 3.206e-05 | 3.206e-05 | 3.206e-05 | 0.0 | 0.00 Modify | 0.062247 | 0.062247 | 0.062247 | 0.0 | 0.06 Other | | 0.005424 | | | 0.01 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 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: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 312.64546378619, Press = -4.98432195737954 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7.5 ghost atom cutoff = 7.5 binsize = 3.75, bins = 9 9 9 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 17000 -17531.801 -17531.801 -17609.866 -17609.866 302.12097 302.12097 32534.918 32534.918 1928.4555 1928.4555 18000 -17531.1 -17531.1 -17611.437 -17611.437 310.91087 310.91087 32562.959 32562.959 -405.33183 -405.33183 Loop time of 100.077 on 1 procs for 1000 steps with 2000 atoms Performance: 0.863 ns/day, 27.799 hours/ns, 9.992 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 | 99.995 | 99.995 | 99.995 | 0.0 | 99.92 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.015055 | 0.015055 | 0.015055 | 0.0 | 0.02 Output | 4.9092e-05 | 4.9092e-05 | 4.9092e-05 | 0.0 | 0.00 Modify | 0.062144 | 0.062144 | 0.062144 | 0.0 | 0.06 Other | | 0.005429 | | | 0.01 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 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: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 312.620198856344, Press = 0.101458833202277 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7.5 ghost atom cutoff = 7.5 binsize = 3.75, bins = 9 9 9 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 18000 -17531.1 -17531.1 -17611.437 -17611.437 310.91087 310.91087 32562.959 32562.959 -405.33183 -405.33183 19000 -17529.16 -17529.16 -17610.076 -17610.076 313.15283 313.15283 32594.185 32594.185 -3646.058 -3646.058 Loop time of 100.07 on 1 procs for 1000 steps with 2000 atoms Performance: 0.863 ns/day, 27.797 hours/ns, 9.993 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 | 99.987 | 99.987 | 99.987 | 0.0 | 99.92 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.015072 | 0.015072 | 0.015072 | 0.0 | 0.02 Output | 4.4484e-05 | 4.4484e-05 | 4.4484e-05 | 0.0 | 0.00 Modify | 0.06229 | 0.06229 | 0.06229 | 0.0 | 0.06 Other | | 0.00545 | | | 0.01 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 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: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 312.849542416129, Press = -4.19283006708606 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7.5 ghost atom cutoff = 7.5 binsize = 3.75, bins = 9 9 9 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 19000 -17529.16 -17529.16 -17610.076 -17610.076 313.15283 313.15283 32594.185 32594.185 -3646.058 -3646.058 20000 -17531.574 -17531.574 -17610.582 -17610.582 305.77012 305.77012 32540.527 32540.527 1212.9317 1212.9317 Loop time of 100.097 on 1 procs for 1000 steps with 2000 atoms Performance: 0.863 ns/day, 27.805 hours/ns, 9.990 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 | 100.01 | 100.01 | 100.01 | 0.0 | 99.92 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.01501 | 0.01501 | 0.01501 | 0.0 | 0.01 Output | 3.6288e-05 | 3.6288e-05 | 3.6288e-05 | 0.0 | 0.00 Modify | 0.062132 | 0.062132 | 0.062132 | 0.0 | 0.06 Other | | 0.005422 | | | 0.01 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 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: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 313.080173427974, Press = -5.19586469814698 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7.5 ghost atom cutoff = 7.5 binsize = 3.75, bins = 9 9 9 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 20000 -17531.574 -17531.574 -17610.582 -17610.582 305.77012 305.77012 32540.527 32540.527 1212.9317 1212.9317 21000 -17532.193 -17532.193 -17612.123 -17612.123 309.33498 309.33498 32532.051 32532.051 2289.1905 2289.1905 Loop time of 100.057 on 1 procs for 1000 steps with 2000 atoms Performance: 0.864 ns/day, 27.794 hours/ns, 9.994 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 | 99.974 | 99.974 | 99.974 | 0.0 | 99.92 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.014977 | 0.014977 | 0.014977 | 0.0 | 0.01 Output | 4.1848e-05 | 4.1848e-05 | 4.1848e-05 | 0.0 | 0.00 Modify | 0.062198 | 0.062198 | 0.062198 | 0.0 | 0.06 Other | | 0.005466 | | | 0.01 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 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: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 313.082102314641, Press = 2.91222912410675 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7.5 ghost atom cutoff = 7.5 binsize = 3.75, bins = 9 9 9 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 21000 -17532.193 -17532.193 -17612.123 -17612.123 309.33498 309.33498 32532.051 32532.051 2289.1905 2289.1905 22000 -17527.69 -17527.69 -17612.433 -17612.433 327.96565 327.96565 32582.269 32582.269 -2463.9286 -2463.9286 Loop time of 100.081 on 1 procs for 1000 steps with 2000 atoms Performance: 0.863 ns/day, 27.800 hours/ns, 9.992 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 | 99.998 | 99.998 | 99.998 | 0.0 | 99.92 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.015067 | 0.015067 | 0.015067 | 0.0 | 0.02 Output | 4.4623e-05 | 4.4623e-05 | 4.4623e-05 | 0.0 | 0.00 Modify | 0.062413 | 0.062413 | 0.062413 | 0.0 | 0.06 Other | | 0.005425 | | | 0.01 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 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: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 313.282065282248, Press = -0.576801172469292 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7.5 ghost atom cutoff = 7.5 binsize = 3.75, bins = 9 9 9 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 22000 -17527.69 -17527.69 -17612.433 -17612.433 327.96565 327.96565 32582.269 32582.269 -2463.9286 -2463.9286 23000 -17534.294 -17534.294 -17611.284 -17611.284 297.95833 297.95833 32546.955 32546.955 756.21192 756.21192 Loop time of 100.041 on 1 procs for 1000 steps with 2000 atoms Performance: 0.864 ns/day, 27.789 hours/ns, 9.996 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 | 99.958 | 99.958 | 99.958 | 0.0 | 99.92 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.015018 | 0.015018 | 0.015018 | 0.0 | 0.02 Output | 3.2741e-05 | 3.2741e-05 | 3.2741e-05 | 0.0 | 0.00 Modify | 0.062501 | 0.062501 | 0.062501 | 0.0 | 0.06 Other | | 0.005428 | | | 0.01 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 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: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 313.34962567132, Press = 1.65519911623719 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7.5 ghost atom cutoff = 7.5 binsize = 3.75, bins = 9 9 9 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 23000 -17534.294 -17534.294 -17611.284 -17611.284 297.95833 297.95833 32546.955 32546.955 756.21192 756.21192 24000 -17531.014 -17531.014 -17612.721 -17612.721 316.21676 316.21676 32592.897 32592.897 -3755.2175 -3755.2175 Loop time of 100.038 on 1 procs for 1000 steps with 2000 atoms Performance: 0.864 ns/day, 27.788 hours/ns, 9.996 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 | 99.955 | 99.955 | 99.955 | 0.0 | 99.92 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.01497 | 0.01497 | 0.01497 | 0.0 | 0.01 Output | 4.292e-05 | 4.292e-05 | 4.292e-05 | 0.0 | 0.00 Modify | 0.062478 | 0.062478 | 0.062478 | 0.0 | 0.06 Other | | 0.005418 | | | 0.01 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 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: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 313.235002426884, Press = -1.37736778046816 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7.5 ghost atom cutoff = 7.5 binsize = 3.75, bins = 9 9 9 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 24000 -17531.014 -17531.014 -17612.721 -17612.721 316.21676 316.21676 32592.897 32592.897 -3755.2175 -3755.2175 25000 -17533.694 -17533.694 -17612.34 -17612.34 304.36666 304.36666 32532.221 32532.221 1948.5023 1948.5023 Loop time of 100.189 on 1 procs for 1000 steps with 2000 atoms Performance: 0.862 ns/day, 27.830 hours/ns, 9.981 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 | 100.11 | 100.11 | 100.11 | 0.0 | 99.92 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.01524 | 0.01524 | 0.01524 | 0.0 | 0.02 Output | 4.6167e-05 | 4.6167e-05 | 4.6167e-05 | 0.0 | 0.00 Modify | 0.062944 | 0.062944 | 0.062944 | 0.0 | 0.06 Other | | 0.00554 | | | 0.01 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 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: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 313.142994289089, Press = -1.30064212095865 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7.5 ghost atom cutoff = 7.5 binsize = 3.75, bins = 9 9 9 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 25000 -17533.694 -17533.694 -17612.34 -17612.34 304.36666 304.36666 32532.221 32532.221 1948.5023 1948.5023 26000 -17530.219 -17530.219 -17612.475 -17612.475 318.34189 318.34189 32558.797 32558.797 -392.11137 -392.11137 Loop time of 100.33 on 1 procs for 1000 steps with 2000 atoms Performance: 0.861 ns/day, 27.869 hours/ns, 9.967 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 | 100.25 | 100.25 | 100.25 | 0.0 | 99.92 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.015173 | 0.015173 | 0.015173 | 0.0 | 0.02 Output | 4.4654e-05 | 4.4654e-05 | 4.4654e-05 | 0.0 | 0.00 Modify | 0.062988 | 0.062988 | 0.062988 | 0.0 | 0.06 Other | | 0.005498 | | | 0.01 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 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: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 313.115213613704, Press = 1.77736603482914 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7.5 ghost atom cutoff = 7.5 binsize = 3.75, bins = 9 9 9 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 26000 -17530.219 -17530.219 -17612.475 -17612.475 318.34189 318.34189 32558.797 32558.797 -392.11137 -392.11137 27000 -17535.87 -17535.87 -17619.311 -17619.311 322.92419 322.92419 32563.143 32563.143 -1128.7633 -1128.7633 Loop time of 100.298 on 1 procs for 1000 steps with 2000 atoms Performance: 0.861 ns/day, 27.861 hours/ns, 9.970 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 | 100.22 | 100.22 | 100.22 | 0.0 | 99.92 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.015057 | 0.015057 | 0.015057 | 0.0 | 0.02 Output | 4.4905e-05 | 4.4905e-05 | 4.4905e-05 | 0.0 | 0.00 Modify | 0.062685 | 0.062685 | 0.062685 | 0.0 | 0.06 Other | | 0.005387 | | | 0.01 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 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: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 313.058977920846, Press = -1.59481788308375 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7.5 ghost atom cutoff = 7.5 binsize = 3.75, bins = 9 9 9 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 27000 -17535.87 -17535.87 -17619.311 -17619.311 322.92419 322.92419 32563.143 32563.143 -1128.7633 -1128.7633 28000 -17531.247 -17531.247 -17614.474 -17614.474 322.09477 322.09477 32524.518 32524.518 2763.2813 2763.2813 Loop time of 100.399 on 1 procs for 1000 steps with 2000 atoms Performance: 0.861 ns/day, 27.888 hours/ns, 9.960 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 | 100.32 | 100.32 | 100.32 | 0.0 | 99.92 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.01504 | 0.01504 | 0.01504 | 0.0 | 0.01 Output | 3.8261e-05 | 3.8261e-05 | 3.8261e-05 | 0.0 | 0.00 Modify | 0.062513 | 0.062513 | 0.062513 | 0.0 | 0.06 Other | | 0.005402 | | | 0.01 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 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: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 312.998631069208, Press = -0.337034962009776 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7.5 ghost atom cutoff = 7.5 binsize = 3.75, bins = 9 9 9 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 28000 -17531.247 -17531.247 -17614.474 -17614.474 322.09477 322.09477 32524.518 32524.518 2763.2813 2763.2813 29000 -17533.93 -17533.93 -17610.921 -17610.921 297.95994 297.95994 32550.396 32550.396 175.58747 175.58747 Loop time of 100.28 on 1 procs for 1000 steps with 2000 atoms Performance: 0.862 ns/day, 27.856 hours/ns, 9.972 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 | 100.2 | 100.2 | 100.2 | 0.0 | 99.92 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.015021 | 0.015021 | 0.015021 | 0.0 | 0.01 Output | 4.1829e-05 | 4.1829e-05 | 4.1829e-05 | 0.0 | 0.00 Modify | 0.062314 | 0.062314 | 0.062314 | 0.0 | 0.06 Other | | 0.005353 | | | 0.01 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 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: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 313.0101519046, Press = 1.06143255682375 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7.5 ghost atom cutoff = 7.5 binsize = 3.75, bins = 9 9 9 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 29000 -17533.93 -17533.93 -17610.921 -17610.921 297.95994 297.95994 32550.396 32550.396 175.58747 175.58747 30000 -17530.886 -17530.886 -17614.226 -17614.226 322.535 322.535 32621.302 32621.302 -6002.9801 -6002.9801 Loop time of 100.265 on 1 procs for 1000 steps with 2000 atoms Performance: 0.862 ns/day, 27.851 hours/ns, 9.974 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 | 100.18 | 100.18 | 100.18 | 0.0 | 99.92 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.015023 | 0.015023 | 0.015023 | 0.0 | 0.01 Output | 4.3232e-05 | 4.3232e-05 | 4.3232e-05 | 0.0 | 0.00 Modify | 0.062377 | 0.062377 | 0.062377 | 0.0 | 0.06 Other | | 0.005383 | | | 0.01 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 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: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 313.019636017737, Press = -1.0565387281549 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7.5 ghost atom cutoff = 7.5 binsize = 3.75, bins = 9 9 9 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 30000 -17530.886 -17530.886 -17614.226 -17614.226 322.535 322.535 32621.302 32621.302 -6002.9801 -6002.9801 31000 -17532.023 -17532.023 -17613.457 -17613.457 315.15952 315.15952 32546.453 32546.453 675.72143 675.72143 Loop time of 100.399 on 1 procs for 1000 steps with 2000 atoms Performance: 0.861 ns/day, 27.889 hours/ns, 9.960 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 | 100.32 | 100.32 | 100.32 | 0.0 | 99.92 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.01508 | 0.01508 | 0.01508 | 0.0 | 0.02 Output | 3.7741e-05 | 3.7741e-05 | 3.7741e-05 | 0.0 | 0.00 Modify | 0.062393 | 0.062393 | 0.062393 | 0.0 | 0.06 Other | | 0.005386 | | | 0.01 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 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: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 313.053962768657, Press = -2.80716905144018 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7.5 ghost atom cutoff = 7.5 binsize = 3.75, bins = 9 9 9 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 31000 -17532.023 -17532.023 -17613.457 -17613.457 315.15952 315.15952 32546.453 32546.453 675.72143 675.72143 32000 -17530.345 -17530.345 -17614.027 -17614.027 323.85839 323.85839 32558.371 32558.371 -392.0452 -392.0452 Loop time of 100.375 on 1 procs for 1000 steps with 2000 atoms Performance: 0.861 ns/day, 27.882 hours/ns, 9.963 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 | 100.29 | 100.29 | 100.29 | 0.0 | 99.92 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.015012 | 0.015012 | 0.015012 | 0.0 | 0.01 Output | 4.252e-05 | 4.252e-05 | 4.252e-05 | 0.0 | 0.00 Modify | 0.062492 | 0.062492 | 0.062492 | 0.0 | 0.06 Other | | 0.005386 | | | 0.01 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 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: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 313.114742597123, Press = 0.763194379368641 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7.5 ghost atom cutoff = 7.5 binsize = 3.75, bins = 9 9 9 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 32000 -17530.345 -17530.345 -17614.027 -17614.027 323.85839 323.85839 32558.371 32558.371 -392.0452 -392.0452 33000 -17530.265 -17530.265 -17610.653 -17610.653 311.10957 311.10957 32565.806 32565.806 -1095.4298 -1095.4298 Loop time of 100.199 on 1 procs for 1000 steps with 2000 atoms Performance: 0.862 ns/day, 27.833 hours/ns, 9.980 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 | 100.12 | 100.12 | 100.12 | 0.0 | 99.92 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.014991 | 0.014991 | 0.014991 | 0.0 | 0.01 Output | 3.9734e-05 | 3.9734e-05 | 3.9734e-05 | 0.0 | 0.00 Modify | 0.062314 | 0.062314 | 0.062314 | 0.0 | 0.06 Other | | 0.005392 | | | 0.01 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 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: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" jump SELF break # Write final averaged volume to file if temperature and volume have converged; otherwise wirte a # flag to indicate non-convergence. variable myStep equal step if "${myStep} < 2000000" then "print '${V}' file output/vol_T313.15.out" else "print 'not_converged' file output/vol_T313.15.out" print '${V}' file output/vol_T313.15.out 32553.8968181098 print "LAMMPS calculation completed" LAMMPS calculation completed quit 0