# 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 5.430949792265892*${_u_distance} variable latticeconst_converted equal 5.430949792265892*1 lattice diamond ${latticeconst_converted} lattice diamond 5.43094979226589 Lattice spacing in x,y,z = 5.4309498 5.4309498 5.4309498 region simbox block 0 10 0 10 0 10 units lattice create_box 1 simbox Created orthogonal box = (0.0000000 0.0000000 0.0000000) to (54.309498 54.309498 54.309498) 1 by 1 by 1 MPI processor grid create_atoms 1 box Created 8000 atoms create_atoms CPU = 0.001 seconds variable mass_converted equal 28.0855*${_u_mass} variable mass_converted equal 28.0855*1 kim_interactions Si WARNING: 'kim_' has been renamed to 'kim '. Please update your input. kim interactions Si #=== BEGIN kim interactions ================================== pair_style kim SW_BalamaneHauchShi_2017Brittle_Si__MO_381114941873_003 pair_coeff * * Si #=== END kim interactions ==================================== mass 1 ${mass_converted} mass 1 28.0855 # initial volume variable v equal vol # assign formula variable V0 equal ${v} # evaluate initial value variable V0 equal 160187.035286095 variable V0_metal equal ${V0}/(${_u_distance}*${_u_distance}*${_u_distance}) variable V0_metal equal 160187.035286095/(${_u_distance}*${_u_distance}*${_u_distance}) variable V0_metal equal 160187.035286095/(1*${_u_distance}*${_u_distance}) variable V0_metal equal 160187.035286095/(1*1*${_u_distance}) variable V0_metal equal 160187.035286095/(1*1*1) variable V0_metal_times1000 equal ${V0_metal}*1000 variable V0_metal_times1000 equal 160187.035286095*1000 print "Initial system volume: ${V0_metal} Angstroms^3" Initial system volume: 160187.035286095 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 273.15*${_u_temperature} variable temp_converted equal 273.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 273.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 273.15 ${temp_converted} ${Tdamp_converted} iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 273.15 273.15 ${Tdamp_converted} iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 273.15 273.15 0.1 iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 273.15 273.15 0.1 iso 0 ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 273.15 273.15 0.1 iso 0 0 ${Pdamp_converted} fix ensemble all npt temp 273.15 273.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 "273.15 - 0.2" variable T_up equal "273.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_381114941873_003#item-citation CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 5.77118 ghost atom cutoff = 5.77118 binsize = 2.88559, bins = 19 19 19 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.77118 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.307 | 4.307 | 4.307 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 0 -36757.574 -36757.574 -37039.998 -37039.998 273.15 273.15 160187.04 160187.04 1883.178 1883.178 1000 -36452.28 -36452.28 -36751.274 -36751.274 289.17614 289.17614 161016.9 161016.9 -718.55543 -718.55543 Loop time of 8.47834 on 1 procs for 1000 steps with 8000 atoms Performance: 10.191 ns/day, 2.355 hours/ns, 117.948 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 | 8.1199 | 8.1199 | 8.1199 | 0.0 | 95.77 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.042162 | 0.042162 | 0.042162 | 0.0 | 0.50 Output | 4.4e-05 | 4.4e-05 | 4.4e-05 | 0.0 | 0.00 Modify | 0.28564 | 0.28564 | 0.28564 | 0.0 | 3.37 Other | | 0.03061 | | | 0.36 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 272000.0 ave 272000 max 272000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 272000 Ave neighs/atom = 34.000000 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 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 5.77118 ghost atom cutoff = 5.77118 binsize = 2.88559, bins = 19 19 19 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.77118 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.309 | 4.309 | 4.309 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 1000 -36452.28 -36452.28 -36751.274 -36751.274 289.17614 289.17614 161016.9 161016.9 -718.55543 -718.55543 2000 -36481.545 -36481.545 -36753.311 -36753.311 262.84206 262.84206 160853.63 160853.63 155.38355 155.38355 Loop time of 9.19477 on 1 procs for 1000 steps with 8000 atoms Performance: 9.397 ns/day, 2.554 hours/ns, 108.758 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 | 8.8421 | 8.8421 | 8.8421 | 0.0 | 96.16 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.041698 | 0.041698 | 0.041698 | 0.0 | 0.45 Output | 7.06e-05 | 7.06e-05 | 7.06e-05 | 0.0 | 0.00 Modify | 0.2806 | 0.2806 | 0.2806 | 0.0 | 3.05 Other | | 0.03033 | | | 0.33 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 272384.0 ave 272384 max 272384 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 272384 Ave neighs/atom = 34.048000 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 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 5.77118 ghost atom cutoff = 5.77118 binsize = 2.88559, bins = 19 19 19 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.77118 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.309 | 4.309 | 4.309 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 2000 -36481.545 -36481.545 -36753.311 -36753.311 262.84206 262.84206 160853.63 160853.63 155.38355 155.38355 3000 -36467.958 -36467.958 -36758.707 -36758.707 281.20119 281.20119 160862.07 160862.07 178.21887 178.21887 Loop time of 9.20673 on 1 procs for 1000 steps with 8000 atoms Performance: 9.384 ns/day, 2.557 hours/ns, 108.616 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 | 8.85 | 8.85 | 8.85 | 0.0 | 96.13 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.041656 | 0.041656 | 0.041656 | 0.0 | 0.45 Output | 4.57e-05 | 4.57e-05 | 4.57e-05 | 0.0 | 0.00 Modify | 0.28457 | 0.28457 | 0.28457 | 0.0 | 3.09 Other | | 0.03048 | | | 0.33 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 272492.0 ave 272492 max 272492 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 272492 Ave neighs/atom = 34.061500 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 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 5.77118 ghost atom cutoff = 5.77118 binsize = 2.88559, bins = 19 19 19 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.77118 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.309 | 4.309 | 4.309 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 3000 -36467.958 -36467.958 -36758.707 -36758.707 281.20119 281.20119 160862.07 160862.07 178.21887 178.21887 4000 -36472.917 -36472.917 -36748.415 -36748.415 266.45108 266.45108 160948.34 160948.34 -433.44314 -433.44314 Loop time of 9.09351 on 1 procs for 1000 steps with 8000 atoms Performance: 9.501 ns/day, 2.526 hours/ns, 109.969 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 | 8.7443 | 8.7443 | 8.7443 | 0.0 | 96.16 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.041433 | 0.041433 | 0.041433 | 0.0 | 0.46 Output | 4.44e-05 | 4.44e-05 | 4.44e-05 | 0.0 | 0.00 Modify | 0.27767 | 0.27767 | 0.27767 | 0.0 | 3.05 Other | | 0.0301 | | | 0.33 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 272410.0 ave 272410 max 272410 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 272410 Ave neighs/atom = 34.051250 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 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 5.77118 ghost atom cutoff = 5.77118 binsize = 2.88559, bins = 19 19 19 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.77118 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.309 | 4.309 | 4.309 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 4000 -36472.917 -36472.917 -36748.415 -36748.415 266.45108 266.45108 160948.34 160948.34 -433.44314 -433.44314 5000 -36475.08 -36475.08 -36762.806 -36762.806 278.27784 278.27784 160780.04 160780.04 654.29799 654.29799 Loop time of 8.95597 on 1 procs for 1000 steps with 8000 atoms Performance: 9.647 ns/day, 2.488 hours/ns, 111.657 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 | 8.6107 | 8.6107 | 8.6107 | 0.0 | 96.14 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.040072 | 0.040072 | 0.040072 | 0.0 | 0.45 Output | 3.56e-05 | 3.56e-05 | 3.56e-05 | 0.0 | 0.00 Modify | 0.27576 | 0.27576 | 0.27576 | 0.0 | 3.08 Other | | 0.02944 | | | 0.33 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 272528.0 ave 272528 max 272528 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 272528 Ave neighs/atom = 34.066000 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 = 272.207889000892, Press = -42.2379701174314 next a jump SELF 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 = 5.77118 ghost atom cutoff = 5.77118 binsize = 2.88559, bins = 19 19 19 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.77118 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.309 | 4.309 | 4.309 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 5000 -36475.08 -36475.08 -36762.806 -36762.806 278.27784 278.27784 160780.04 160780.04 654.29799 654.29799 6000 -36467.827 -36467.827 -36750.139 -36750.139 273.04223 273.04223 161014.86 161014.86 -822.49261 -822.49261 Loop time of 8.99137 on 1 procs for 1000 steps with 8000 atoms Performance: 9.609 ns/day, 2.498 hours/ns, 111.218 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 | 8.6402 | 8.6402 | 8.6402 | 0.0 | 96.09 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.040262 | 0.040262 | 0.040262 | 0.0 | 0.45 Output | 3.53e-05 | 3.53e-05 | 3.53e-05 | 0.0 | 0.00 Modify | 0.28132 | 0.28132 | 0.28132 | 0.0 | 3.13 Other | | 0.02952 | | | 0.33 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 272420.0 ave 272420 max 272420 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 272420 Ave neighs/atom = 34.052500 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 = 272.856273122449, Press = -31.6949699105577 next a jump SELF 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 = 5.77118 ghost atom cutoff = 5.77118 binsize = 2.88559, bins = 19 19 19 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.77118 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.309 | 4.309 | 4.309 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 6000 -36467.827 -36467.827 -36750.139 -36750.139 273.04223 273.04223 161014.86 161014.86 -822.49261 -822.49261 7000 -36477.36 -36477.36 -36770.053 -36770.053 283.08259 283.08259 160873.5 160873.5 4.3497769 4.3497769 Loop time of 8.93883 on 1 procs for 1000 steps with 8000 atoms Performance: 9.666 ns/day, 2.483 hours/ns, 111.872 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 | 8.589 | 8.589 | 8.589 | 0.0 | 96.09 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.040351 | 0.040351 | 0.040351 | 0.0 | 0.45 Output | 4.39e-05 | 4.39e-05 | 4.39e-05 | 0.0 | 0.00 Modify | 0.27998 | 0.27998 | 0.27998 | 0.0 | 3.13 Other | | 0.02941 | | | 0.33 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 272378.0 ave 272378 max 272378 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 272378 Ave neighs/atom = 34.047250 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 = 273.016867185141, Press = 6.67470040831431 next a jump SELF 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 = 5.77118 ghost atom cutoff = 5.77118 binsize = 2.88559, bins = 19 19 19 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.77118 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.309 | 4.309 | 4.309 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 7000 -36477.36 -36477.36 -36770.053 -36770.053 283.08259 283.08259 160873.5 160873.5 4.3497769 4.3497769 8000 -36468.322 -36468.322 -36741.849 -36741.849 264.54509 264.54509 160852.36 160852.36 278.3772 278.3772 Loop time of 8.85061 on 1 procs for 1000 steps with 8000 atoms Performance: 9.762 ns/day, 2.459 hours/ns, 112.987 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 | 8.5043 | 8.5043 | 8.5043 | 0.0 | 96.09 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.039768 | 0.039768 | 0.039768 | 0.0 | 0.45 Output | 4.98e-05 | 4.98e-05 | 4.98e-05 | 0.0 | 0.00 Modify | 0.27731 | 0.27731 | 0.27731 | 0.0 | 3.13 Other | | 0.02921 | | | 0.33 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 272342.0 ave 272342 max 272342 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 272342 Ave neighs/atom = 34.042750 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 = 273.103994086594, Press = -5.62537044060859 next a jump SELF 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 = 5.77118 ghost atom cutoff = 5.77118 binsize = 2.88559, bins = 19 19 19 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.77118 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.309 | 4.309 | 4.309 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 8000 -36468.322 -36468.322 -36741.849 -36741.849 264.54509 264.54509 160852.36 160852.36 278.3772 278.3772 9000 -36474.557 -36474.557 -36755.965 -36755.965 272.16768 272.16768 160945.07 160945.07 -418.21758 -418.21758 Loop time of 8.71765 on 1 procs for 1000 steps with 8000 atoms Performance: 9.911 ns/day, 2.422 hours/ns, 114.710 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 | 8.375 | 8.375 | 8.375 | 0.0 | 96.07 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.039236 | 0.039236 | 0.039236 | 0.0 | 0.45 Output | 3.56e-05 | 3.56e-05 | 3.56e-05 | 0.0 | 0.00 Modify | 0.27457 | 0.27457 | 0.27457 | 0.0 | 3.15 Other | | 0.02879 | | | 0.33 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 272594.0 ave 272594 max 272594 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 272594 Ave neighs/atom = 34.074250 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 = 273.06225388833, Press = -1.59717663526793 next a jump SELF 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 = 5.77118 ghost atom cutoff = 5.77118 binsize = 2.88559, bins = 19 19 19 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.77118 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.309 | 4.309 | 4.309 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 9000 -36474.557 -36474.557 -36755.965 -36755.965 272.16768 272.16768 160945.07 160945.07 -418.21758 -418.21758 10000 -36473.957 -36473.957 -36759.12 -36759.12 275.79981 275.79981 160806.55 160806.55 492.88391 492.88391 Loop time of 8.8301 on 1 procs for 1000 steps with 8000 atoms Performance: 9.785 ns/day, 2.453 hours/ns, 113.249 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 | 8.4832 | 8.4832 | 8.4832 | 0.0 | 96.07 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.040632 | 0.040632 | 0.040632 | 0.0 | 0.46 Output | 4.5e-05 | 4.5e-05 | 4.5e-05 | 0.0 | 0.00 Modify | 0.27717 | 0.27717 | 0.27717 | 0.0 | 3.14 Other | | 0.02905 | | | 0.33 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 272434.0 ave 272434 max 272434 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 272434 Ave neighs/atom = 34.054250 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 = 273.143520238484, Press = -0.773238291136324 next a jump SELF 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 = 5.77118 ghost atom cutoff = 5.77118 binsize = 2.88559, bins = 19 19 19 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.77118 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.309 | 4.309 | 4.309 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 10000 -36473.957 -36473.957 -36759.12 -36759.12 275.79981 275.79981 160806.55 160806.55 492.88391 492.88391 11000 -36465.726 -36465.726 -36750.229 -36750.229 275.16049 275.16049 160966.1 160966.1 -495.08388 -495.08388 Loop time of 9.03853 on 1 procs for 1000 steps with 8000 atoms Performance: 9.559 ns/day, 2.511 hours/ns, 110.637 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 | 8.6867 | 8.6867 | 8.6867 | 0.0 | 96.11 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.040531 | 0.040531 | 0.040531 | 0.0 | 0.45 Output | 3.54e-05 | 3.54e-05 | 3.54e-05 | 0.0 | 0.00 Modify | 0.28151 | 0.28151 | 0.28151 | 0.0 | 3.11 Other | | 0.02979 | | | 0.33 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 272484.0 ave 272484 max 272484 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 272484 Ave neighs/atom = 34.060500 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 = 273.226833518575, Press = -4.80345898994712 next a jump SELF 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 = 5.77118 ghost atom cutoff = 5.77118 binsize = 2.88559, bins = 19 19 19 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.77118 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.309 | 4.309 | 4.309 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 11000 -36465.726 -36465.726 -36750.229 -36750.229 275.16049 275.16049 160966.1 160966.1 -495.08388 -495.08388 12000 -36472.695 -36472.695 -36762.733 -36762.733 280.51383 280.51383 160912.45 160912.45 -191.76471 -191.76471 Loop time of 8.78566 on 1 procs for 1000 steps with 8000 atoms Performance: 9.834 ns/day, 2.440 hours/ns, 113.822 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 | 8.4414 | 8.4414 | 8.4414 | 0.0 | 96.08 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.039374 | 0.039374 | 0.039374 | 0.0 | 0.45 Output | 3.56e-05 | 3.56e-05 | 3.56e-05 | 0.0 | 0.00 Modify | 0.27594 | 0.27594 | 0.27594 | 0.0 | 3.14 Other | | 0.0289 | | | 0.33 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 272434.0 ave 272434 max 272434 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 272434 Ave neighs/atom = 34.054250 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 = 273.330249733871, Press = 2.84646536962901 next a jump SELF 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 = 5.77118 ghost atom cutoff = 5.77118 binsize = 2.88559, bins = 19 19 19 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.77118 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.309 | 4.309 | 4.309 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 12000 -36472.695 -36472.695 -36762.733 -36762.733 280.51383 280.51383 160912.45 160912.45 -191.76471 -191.76471 13000 -36473.585 -36473.585 -36752.698 -36752.698 269.94761 269.94761 160832.1 160832.1 368.94382 368.94382 Loop time of 8.62318 on 1 procs for 1000 steps with 8000 atoms Performance: 10.019 ns/day, 2.395 hours/ns, 115.966 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 | 8.2855 | 8.2855 | 8.2855 | 0.0 | 96.08 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.038596 | 0.038596 | 0.038596 | 0.0 | 0.45 Output | 3.48e-05 | 3.48e-05 | 3.48e-05 | 0.0 | 0.00 Modify | 0.27097 | 0.27097 | 0.27097 | 0.0 | 3.14 Other | | 0.02811 | | | 0.33 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 272312.0 ave 272312 max 272312 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 272312 Ave neighs/atom = 34.039000 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 = 273.420366953341, Press = -2.13358956541811 next a jump SELF 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 = 5.77118 ghost atom cutoff = 5.77118 binsize = 2.88559, bins = 19 19 19 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.77118 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.309 | 4.309 | 4.309 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 13000 -36473.585 -36473.585 -36752.698 -36752.698 269.94761 269.94761 160832.1 160832.1 368.94382 368.94382 14000 -36466.301 -36466.301 -36753.084 -36753.084 277.36638 277.36638 160942.44 160942.44 -318.64371 -318.64371 Loop time of 8.59319 on 1 procs for 1000 steps with 8000 atoms Performance: 10.054 ns/day, 2.387 hours/ns, 116.371 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 | 8.2545 | 8.2545 | 8.2545 | 0.0 | 96.06 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.038794 | 0.038794 | 0.038794 | 0.0 | 0.45 Output | 3.53e-05 | 3.53e-05 | 3.53e-05 | 0.0 | 0.00 Modify | 0.27167 | 0.27167 | 0.27167 | 0.0 | 3.16 Other | | 0.0282 | | | 0.33 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 272466.0 ave 272466 max 272466 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 272466 Ave neighs/atom = 34.058250 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 = 273.411956231393, Press = -0.874991000853215 next a jump SELF 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 = 5.77118 ghost atom cutoff = 5.77118 binsize = 2.88559, bins = 19 19 19 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.77118 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.309 | 4.309 | 4.309 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 14000 -36466.301 -36466.301 -36753.084 -36753.084 277.36638 277.36638 160942.44 160942.44 -318.64371 -318.64371 15000 -36472.026 -36472.026 -36753.015 -36753.015 271.76286 271.76286 160838.32 160838.32 350.92651 350.92651 Loop time of 8.68982 on 1 procs for 1000 steps with 8000 atoms Performance: 9.943 ns/day, 2.414 hours/ns, 115.077 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 | 8.3496 | 8.3496 | 8.3496 | 0.0 | 96.09 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.038961 | 0.038961 | 0.038961 | 0.0 | 0.45 Output | 3.51e-05 | 3.51e-05 | 3.51e-05 | 0.0 | 0.00 Modify | 0.27288 | 0.27288 | 0.27288 | 0.0 | 3.14 Other | | 0.0283 | | | 0.33 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 272388.0 ave 272388 max 272388 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 272388 Ave neighs/atom = 34.048500 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 = 273.530342498272, Press = -0.318728115931118 next a jump SELF 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 = 5.77118 ghost atom cutoff = 5.77118 binsize = 2.88559, bins = 19 19 19 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.77118 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.309 | 4.309 | 4.309 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 15000 -36472.026 -36472.026 -36753.015 -36753.015 271.76286 271.76286 160838.32 160838.32 350.92651 350.92651 16000 -36468.064 -36468.064 -36746.142 -36746.142 268.94739 268.94739 160925.67 160925.67 -236.72273 -236.72273 Loop time of 8.59796 on 1 procs for 1000 steps with 8000 atoms Performance: 10.049 ns/day, 2.388 hours/ns, 116.307 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 | 8.2584 | 8.2584 | 8.2584 | 0.0 | 96.05 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.038899 | 0.038899 | 0.038899 | 0.0 | 0.45 Output | 4.51e-05 | 4.51e-05 | 4.51e-05 | 0.0 | 0.00 Modify | 0.27238 | 0.27238 | 0.27238 | 0.0 | 3.17 Other | | 0.02823 | | | 0.33 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 272422.0 ave 272422 max 272422 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 272422 Ave neighs/atom = 34.052750 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 = 273.51835871211, Press = -2.28356062097253 next a jump SELF 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 = 5.77118 ghost atom cutoff = 5.77118 binsize = 2.88559, bins = 19 19 19 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.77118 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.309 | 4.309 | 4.309 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 16000 -36468.064 -36468.064 -36746.142 -36746.142 268.94739 268.94739 160925.67 160925.67 -236.72273 -236.72273 17000 -36472.5 -36472.5 -36757.532 -36757.532 275.67283 275.67283 160946.46 160946.46 -415.27098 -415.27098 Loop time of 8.6839 on 1 procs for 1000 steps with 8000 atoms Performance: 9.949 ns/day, 2.412 hours/ns, 115.156 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 | 8.3435 | 8.3435 | 8.3435 | 0.0 | 96.08 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.039027 | 0.039027 | 0.039027 | 0.0 | 0.45 Output | 3.54e-05 | 3.54e-05 | 3.54e-05 | 0.0 | 0.00 Modify | 0.27303 | 0.27303 | 0.27303 | 0.0 | 3.14 Other | | 0.02831 | | | 0.33 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 272484.0 ave 272484 max 272484 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 272484 Ave neighs/atom = 34.060500 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 = 273.484137424363, Press = 1.37530134398808 next a jump SELF 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 = 5.77118 ghost atom cutoff = 5.77118 binsize = 2.88559, bins = 19 19 19 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.77118 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.309 | 4.309 | 4.309 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 17000 -36472.5 -36472.5 -36757.532 -36757.532 275.67283 275.67283 160946.46 160946.46 -415.27098 -415.27098 18000 -36464.072 -36464.072 -36743.656 -36743.656 270.40331 270.40331 160785.79 160785.79 756.11859 756.11859 Loop time of 8.72866 on 1 procs for 1000 steps with 8000 atoms Performance: 9.898 ns/day, 2.425 hours/ns, 114.565 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 | 8.3867 | 8.3867 | 8.3867 | 0.0 | 96.08 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.039223 | 0.039223 | 0.039223 | 0.0 | 0.45 Output | 3.59e-05 | 3.59e-05 | 3.59e-05 | 0.0 | 0.00 Modify | 0.27388 | 0.27388 | 0.27388 | 0.0 | 3.14 Other | | 0.02884 | | | 0.33 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 272400.0 ave 272400 max 272400 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 272400 Ave neighs/atom = 34.050000 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 = 273.525295364415, Press = -1.44112319754964 next a jump SELF 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 = 5.77118 ghost atom cutoff = 5.77118 binsize = 2.88559, bins = 19 19 19 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.77118 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.309 | 4.309 | 4.309 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 18000 -36464.072 -36464.072 -36743.656 -36743.656 270.40331 270.40331 160785.79 160785.79 756.11859 756.11859 19000 -36471.921 -36471.921 -36759.148 -36759.148 277.79515 277.79515 160995.38 160995.38 -741.24322 -741.24322 Loop time of 8.37507 on 1 procs for 1000 steps with 8000 atoms Performance: 10.316 ns/day, 2.326 hours/ns, 119.402 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 | 8.0417 | 8.0417 | 8.0417 | 0.0 | 96.02 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.038194 | 0.038194 | 0.038194 | 0.0 | 0.46 Output | 3.57e-05 | 3.57e-05 | 3.57e-05 | 0.0 | 0.00 Modify | 0.26732 | 0.26732 | 0.26732 | 0.0 | 3.19 Other | | 0.02786 | | | 0.33 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 272510.0 ave 272510 max 272510 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 272510 Ave neighs/atom = 34.063750 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 = 273.516846142743, Press = -0.783922522740365 next a jump SELF 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 = 5.77118 ghost atom cutoff = 5.77118 binsize = 2.88559, bins = 19 19 19 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.77118 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.309 | 4.309 | 4.309 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 19000 -36471.921 -36471.921 -36759.148 -36759.148 277.79515 277.79515 160995.38 160995.38 -741.24322 -741.24322 20000 -36470.231 -36470.231 -36756.523 -36756.523 276.89074 276.89074 160815.97 160815.97 479.28451 479.28451 Loop time of 8.45605 on 1 procs for 1000 steps with 8000 atoms Performance: 10.218 ns/day, 2.349 hours/ns, 118.258 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 | 8.121 | 8.121 | 8.121 | 0.0 | 96.04 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.038222 | 0.038222 | 0.038222 | 0.0 | 0.45 Output | 4.41e-05 | 4.41e-05 | 4.41e-05 | 0.0 | 0.00 Modify | 0.26899 | 0.26899 | 0.26899 | 0.0 | 3.18 Other | | 0.02779 | | | 0.33 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 272438.0 ave 272438 max 272438 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 272438 Ave neighs/atom = 34.054750 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 = 273.519881084122, Press = 0.229389555872102 next a jump SELF 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 = 5.77118 ghost atom cutoff = 5.77118 binsize = 2.88559, bins = 19 19 19 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.77118 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.309 | 4.309 | 4.309 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 20000 -36470.231 -36470.231 -36756.523 -36756.523 276.89074 276.89074 160815.97 160815.97 479.28451 479.28451 21000 -36473.774 -36473.774 -36755.467 -36755.467 272.44271 272.44271 160900.24 160900.24 -101.07594 -101.07594 Loop time of 8.487 on 1 procs for 1000 steps with 8000 atoms Performance: 10.180 ns/day, 2.358 hours/ns, 117.827 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 | 8.1524 | 8.1524 | 8.1524 | 0.0 | 96.06 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.038603 | 0.038603 | 0.038603 | 0.0 | 0.45 Output | 3.6e-05 | 3.6e-05 | 3.6e-05 | 0.0 | 0.00 Modify | 0.26816 | 0.26816 | 0.26816 | 0.0 | 3.16 Other | | 0.02775 | | | 0.33 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 272470.0 ave 272470 max 272470 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 272470 Ave neighs/atom = 34.058750 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 = 273.444712984643, Press = -1.47895415818069 next a jump SELF 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 = 5.77118 ghost atom cutoff = 5.77118 binsize = 2.88559, bins = 19 19 19 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.77118 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.309 | 4.309 | 4.309 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 21000 -36473.774 -36473.774 -36755.467 -36755.467 272.44271 272.44271 160900.24 160900.24 -101.07594 -101.07594 22000 -36471.588 -36471.588 -36759.098 -36759.098 278.06894 278.06894 160938.31 160938.31 -359.40737 -359.40737 Loop time of 8.52032 on 1 procs for 1000 steps with 8000 atoms Performance: 10.140 ns/day, 2.367 hours/ns, 117.366 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 | 8.1851 | 8.1851 | 8.1851 | 0.0 | 96.07 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.038334 | 0.038334 | 0.038334 | 0.0 | 0.45 Output | 3.55e-05 | 3.55e-05 | 3.55e-05 | 0.0 | 0.00 Modify | 0.26899 | 0.26899 | 0.26899 | 0.0 | 3.16 Other | | 0.02784 | | | 0.33 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 272438.0 ave 272438 max 272438 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 272438 Ave neighs/atom = 34.054750 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 = 273.448452543565, Press = 0.390186693534871 next a jump SELF 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 = 5.77118 ghost atom cutoff = 5.77118 binsize = 2.88559, bins = 19 19 19 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.77118 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.309 | 4.309 | 4.309 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 22000 -36471.588 -36471.588 -36759.098 -36759.098 278.06894 278.06894 160938.31 160938.31 -359.40737 -359.40737 23000 -36466.057 -36466.057 -36751.085 -36751.085 275.66837 275.66837 160728.8 160728.8 1117.1744 1117.1744 Loop time of 8.62049 on 1 procs for 1000 steps with 8000 atoms Performance: 10.023 ns/day, 2.395 hours/ns, 116.003 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 | 8.2818 | 8.2818 | 8.2818 | 0.0 | 96.07 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.038718 | 0.038718 | 0.038718 | 0.0 | 0.45 Output | 3.51e-05 | 3.51e-05 | 3.51e-05 | 0.0 | 0.00 Modify | 0.27184 | 0.27184 | 0.27184 | 0.0 | 3.15 Other | | 0.02811 | | | 0.33 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 272354.0 ave 272354 max 272354 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 272354 Ave neighs/atom = 34.044250 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 = 273.430338829487, Press = -0.964121803760383 next a jump SELF 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 = 5.77118 ghost atom cutoff = 5.77118 binsize = 2.88559, bins = 19 19 19 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.77118 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.309 | 4.309 | 4.309 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 23000 -36466.057 -36466.057 -36751.085 -36751.085 275.66837 275.66837 160728.8 160728.8 1117.1744 1117.1744 24000 -36476.073 -36476.073 -36758.563 -36758.563 273.21438 273.21438 161034.51 161034.51 -1031.6762 -1031.6762 Loop time of 8.58667 on 1 procs for 1000 steps with 8000 atoms Performance: 10.062 ns/day, 2.385 hours/ns, 116.460 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 | 8.2484 | 8.2484 | 8.2484 | 0.0 | 96.06 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.0386 | 0.0386 | 0.0386 | 0.0 | 0.45 Output | 3.51e-05 | 3.51e-05 | 3.51e-05 | 0.0 | 0.00 Modify | 0.27163 | 0.27163 | 0.27163 | 0.0 | 3.16 Other | | 0.02797 | | | 0.33 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 272448.0 ave 272448 max 272448 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 272448 Ave neighs/atom = 34.056000 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 = 273.406883497563, Press = -0.818517851793147 next a jump SELF 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 = 5.77118 ghost atom cutoff = 5.77118 binsize = 2.88559, bins = 19 19 19 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.77118 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.309 | 4.309 | 4.309 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 24000 -36476.073 -36476.073 -36758.563 -36758.563 273.21438 273.21438 161034.51 161034.51 -1031.6762 -1031.6762 25000 -36469.604 -36469.604 -36753.157 -36753.157 274.24248 274.24248 160844.2 160844.2 302.59548 302.59548 Loop time of 8.51792 on 1 procs for 1000 steps with 8000 atoms Performance: 10.143 ns/day, 2.366 hours/ns, 117.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 | 8.1819 | 8.1819 | 8.1819 | 0.0 | 96.06 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.038457 | 0.038457 | 0.038457 | 0.0 | 0.45 Output | 4.4e-05 | 4.4e-05 | 4.4e-05 | 0.0 | 0.00 Modify | 0.2697 | 0.2697 | 0.2697 | 0.0 | 3.17 Other | | 0.0278 | | | 0.33 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 272350.0 ave 272350 max 272350 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 272350 Ave neighs/atom = 34.043750 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 = 273.363479608757, Press = 0.0890886720576291 next a jump SELF 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 = 5.77118 ghost atom cutoff = 5.77118 binsize = 2.88559, bins = 19 19 19 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.77118 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.309 | 4.309 | 4.309 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 25000 -36469.604 -36469.604 -36753.157 -36753.157 274.24248 274.24248 160844.2 160844.2 302.59548 302.59548 26000 -36476.257 -36476.257 -36757.425 -36757.425 271.93532 271.93532 160878.43 160878.43 28.431292 28.431292 Loop time of 8.90567 on 1 procs for 1000 steps with 8000 atoms Performance: 9.702 ns/day, 2.474 hours/ns, 112.288 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 | 8.5596 | 8.5596 | 8.5596 | 0.0 | 96.11 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.03934 | 0.03934 | 0.03934 | 0.0 | 0.44 Output | 7.48e-05 | 7.48e-05 | 7.48e-05 | 0.0 | 0.00 Modify | 0.27777 | 0.27777 | 0.27777 | 0.0 | 3.12 Other | | 0.02888 | | | 0.32 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 272480.0 ave 272480 max 272480 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 272480 Ave neighs/atom = 34.060000 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 = 273.291997507393, Press = -0.80150090293515 next a jump SELF 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 = 5.77118 ghost atom cutoff = 5.77118 binsize = 2.88559, bins = 19 19 19 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.77118 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.309 | 4.309 | 4.309 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 26000 -36476.257 -36476.257 -36757.425 -36757.425 271.93532 271.93532 160878.43 160878.43 28.431292 28.431292 27000 -36474.929 -36474.929 -36759.539 -36759.539 275.26412 275.26412 160936.86 160936.86 -369.3931 -369.3931 Loop time of 8.6034 on 1 procs for 1000 steps with 8000 atoms Performance: 10.043 ns/day, 2.390 hours/ns, 116.233 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 | 8.2648 | 8.2648 | 8.2648 | 0.0 | 96.06 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.038745 | 0.038745 | 0.038745 | 0.0 | 0.45 Output | 3.56e-05 | 3.56e-05 | 3.56e-05 | 0.0 | 0.00 Modify | 0.27163 | 0.27163 | 0.27163 | 0.0 | 3.16 Other | | 0.02817 | | | 0.33 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 272448.0 ave 272448 max 272448 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 272448 Ave neighs/atom = 34.056000 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_T273.15.out" else "print 'not_converged' file output/vol_T273.15.out" print '${V}' file output/vol_T273.15.out 160886.546386417 print "LAMMPS calculation completed" LAMMPS calculation completed quit 0