# 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.427092507481575*${_u_distance} variable latticeconst_converted equal 5.427092507481575*1 lattice diamond ${latticeconst_converted} lattice diamond 5.42709250748158 Lattice spacing in x,y,z = 5.4270925 5.4270925 5.4270925 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.270925 54.270925 54.270925) 1 by 1 by 1 MPI processor grid create_atoms 1 box Created 8000 atoms using lattice units in orthogonal box = (0.0000000 0.0000000 0.0000000) to (54.270925 54.270925 54.270925) create_atoms CPU = 0.008 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 MEAM_LAMMPS_Lee_2007_Si__MO_774917820956_000 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 159845.963304675 variable V0_metal equal ${V0}/(${_u_distance}*${_u_distance}*${_u_distance}) variable V0_metal equal 159845.963304675/(${_u_distance}*${_u_distance}*${_u_distance}) variable V0_metal equal 159845.963304675/(1*${_u_distance}*${_u_distance}) variable V0_metal equal 159845.963304675/(1*1*${_u_distance}) variable V0_metal equal 159845.963304675/(1*1*1) variable V0_metal_times1000 equal ${V0_metal}*1000 variable V0_metal_times1000 equal 159845.963304675*1000 print "Initial system volume: ${V0_metal} Angstroms^3" Initial system volume: 159845.963304675 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_774917820956_000#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 = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 17 17 17 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.677 | 4.677 | 4.677 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 0 -36757.576 -36757.576 -37040 -37040 273.15 273.15 159845.96 159845.96 1887.2003 1887.2003 1000 -36450.913 -36450.913 -36755.507 -36755.507 294.5924 294.5924 160516.04 160516.04 -337.04923 -337.04923 Loop time of 146.813 on 1 procs for 1000 steps with 8000 atoms Performance: 0.589 ns/day, 40.781 hours/ns, 6.811 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 | 145.09 | 145.09 | 145.09 | 0.0 | 98.82 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.21816 | 0.21816 | 0.21816 | 0.0 | 0.15 Output | 0.00027922 | 0.00027922 | 0.00027922 | 0.0 | 0.00 Modify | 1.3475 | 1.3475 | 1.3475 | 0.0 | 0.92 Other | | 0.1621 | | | 0.11 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: 368000.0 ave 368000 max 368000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 368000 Ave neighs/atom = 46.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 = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 17 17 17 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.679 | 4.679 | 4.679 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 1000 -36450.913 -36450.913 -36755.507 -36755.507 294.5924 294.5924 160516.04 160516.04 -337.04923 -337.04923 2000 -36481.026 -36481.026 -36749.148 -36749.148 259.31753 259.31753 160503.89 160503.89 -264.04478 -264.04478 Loop time of 144.517 on 1 procs for 1000 steps with 8000 atoms Performance: 0.598 ns/day, 40.143 hours/ns, 6.920 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 | 142.84 | 142.84 | 142.84 | 0.0 | 98.84 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.21411 | 0.21411 | 0.21411 | 0.0 | 0.15 Output | 0.00031446 | 0.00031446 | 0.00031446 | 0.0 | 0.00 Modify | 1.3071 | 1.3071 | 1.3071 | 0.0 | 0.90 Other | | 0.1583 | | | 0.11 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: 372558.0 ave 372558 max 372558 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 372558 Ave neighs/atom = 46.569750 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 = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 17 17 17 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.679 | 4.679 | 4.679 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 2000 -36481.026 -36481.026 -36749.148 -36749.148 259.31753 259.31753 160503.89 160503.89 -264.04478 -264.04478 3000 -36466.413 -36466.413 -36748.454 -36748.454 272.78005 272.78005 160427.92 160427.92 247.57569 247.57569 Loop time of 144.373 on 1 procs for 1000 steps with 8000 atoms Performance: 0.598 ns/day, 40.104 hours/ns, 6.926 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 | 142.7 | 142.7 | 142.7 | 0.0 | 98.84 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.21371 | 0.21371 | 0.21371 | 0.0 | 0.15 Output | 0.00027356 | 0.00027356 | 0.00027356 | 0.0 | 0.00 Modify | 1.3045 | 1.3045 | 1.3045 | 0.0 | 0.90 Other | | 0.159 | | | 0.11 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: 372346.0 ave 372346 max 372346 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 372346 Ave neighs/atom = 46.543250 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 = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 17 17 17 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.679 | 4.679 | 4.679 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 3000 -36466.413 -36466.413 -36748.454 -36748.454 272.78005 272.78005 160427.92 160427.92 247.57569 247.57569 4000 -36473.134 -36473.134 -36755.064 -36755.064 272.67276 272.67276 160509.62 160509.62 -337.35843 -337.35843 Loop time of 144.334 on 1 procs for 1000 steps with 8000 atoms Performance: 0.599 ns/day, 40.093 hours/ns, 6.928 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 | 142.66 | 142.66 | 142.66 | 0.0 | 98.84 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.21291 | 0.21291 | 0.21291 | 0.0 | 0.15 Output | 0.00031012 | 0.00031012 | 0.00031012 | 0.0 | 0.00 Modify | 1.3027 | 1.3027 | 1.3027 | 0.0 | 0.90 Other | | 0.1585 | | | 0.11 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: 372752.0 ave 372752 max 372752 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 372752 Ave neighs/atom = 46.594000 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 = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 17 17 17 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.679 | 4.679 | 4.679 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 4000 -36473.134 -36473.134 -36755.064 -36755.064 272.67276 272.67276 160509.62 160509.62 -337.35843 -337.35843 5000 -36472.808 -36472.808 -36756.473 -36756.473 274.34988 274.34988 160545.5 160545.5 -626.37868 -626.37868 Loop time of 161.434 on 1 procs for 1000 steps with 8000 atoms Performance: 0.535 ns/day, 44.843 hours/ns, 6.194 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 | 159.48 | 159.48 | 159.48 | 0.0 | 98.79 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.23045 | 0.23045 | 0.23045 | 0.0 | 0.14 Output | 0.00035201 | 0.00035201 | 0.00035201 | 0.0 | 0.00 Modify | 1.5475 | 1.5475 | 1.5475 | 0.0 | 0.96 Other | | 0.1771 | | | 0.11 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: 372294.0 ave 372294 max 372294 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 372294 Ave neighs/atom = 46.536750 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.956110788359, Press = -428.192328327475 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 = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 17 17 17 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.679 | 4.679 | 4.679 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 5000 -36472.808 -36472.808 -36756.473 -36756.473 274.34988 274.34988 160545.5 160545.5 -626.37868 -626.37868 6000 -36467.975 -36467.975 -36748.778 -36748.778 271.58253 271.58253 160333.94 160333.94 672.62122 672.62122 Loop time of 166.676 on 1 procs for 1000 steps with 8000 atoms Performance: 0.518 ns/day, 46.299 hours/ns, 6.000 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 | 164.62 | 164.62 | 164.62 | 0.0 | 98.76 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.2352 | 0.2352 | 0.2352 | 0.0 | 0.14 Output | 0.00022913 | 0.00022913 | 0.00022913 | 0.0 | 0.00 Modify | 1.6429 | 1.6429 | 1.6429 | 0.0 | 0.99 Other | | 0.1819 | | | 0.11 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: 372394.0 ave 372394 max 372394 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 372394 Ave neighs/atom = 46.549250 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.749140727159, Press = -19.5563768519375 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 = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 17 17 17 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.679 | 4.679 | 4.679 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 6000 -36467.975 -36467.975 -36748.778 -36748.778 271.58253 271.58253 160333.94 160333.94 672.62122 672.62122 7000 -36475.261 -36475.261 -36753.501 -36753.501 269.10395 269.10395 160492.63 160492.63 -266.24608 -266.24608 Loop time of 166.523 on 1 procs for 1000 steps with 8000 atoms Performance: 0.519 ns/day, 46.257 hours/ns, 6.005 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 | 164.47 | 164.47 | 164.47 | 0.0 | 98.77 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.23521 | 0.23521 | 0.23521 | 0.0 | 0.14 Output | 0.00031385 | 0.00031385 | 0.00031385 | 0.0 | 0.00 Modify | 1.6367 | 1.6367 | 1.6367 | 0.0 | 0.98 Other | | 0.1812 | | | 0.11 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: 373132.0 ave 373132 max 373132 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 373132 Ave neighs/atom = 46.641500 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.095899627594, Press = -2.43717130931009 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 = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 17 17 17 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.679 | 4.679 | 4.679 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 7000 -36475.261 -36475.261 -36753.501 -36753.501 269.10395 269.10395 160492.63 160492.63 -266.24608 -266.24608 8000 -36465.889 -36465.889 -36746.175 -36746.175 271.08278 271.08278 160521.51 160521.51 -314.63887 -314.63887 Loop time of 150.541 on 1 procs for 1000 steps with 8000 atoms Performance: 0.574 ns/day, 41.817 hours/ns, 6.643 timesteps/s 99.9% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 148.74 | 148.74 | 148.74 | 0.0 | 98.80 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.21302 | 0.21302 | 0.21302 | 0.0 | 0.14 Output | 0.00023075 | 0.00023075 | 0.00023075 | 0.0 | 0.00 Modify | 1.4216 | 1.4216 | 1.4216 | 0.0 | 0.94 Other | | 0.1667 | | | 0.11 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: 372540.0 ave 372540 max 372540 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 372540 Ave neighs/atom = 46.567500 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.084686132902, Press = -5.64353353283266 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 = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 17 17 17 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.679 | 4.679 | 4.679 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 8000 -36465.889 -36465.889 -36746.175 -36746.175 271.08278 271.08278 160521.51 160521.51 -314.63887 -314.63887 9000 -36474.261 -36474.261 -36763.603 -36763.603 279.84112 279.84112 160408.42 160408.42 359.98864 359.98864 Loop time of 165.836 on 1 procs for 1000 steps with 8000 atoms Performance: 0.521 ns/day, 46.065 hours/ns, 6.030 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 | 163.77 | 163.77 | 163.77 | 0.0 | 98.76 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.23778 | 0.23778 | 0.23778 | 0.0 | 0.14 Output | 0.00034862 | 0.00034862 | 0.00034862 | 0.0 | 0.00 Modify | 1.6384 | 1.6384 | 1.6384 | 0.0 | 0.99 Other | | 0.1842 | | | 0.11 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: 372374.0 ave 372374 max 372374 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 372374 Ave neighs/atom = 46.546750 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.059006112314, Press = -5.51265193821556 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 = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 17 17 17 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.679 | 4.679 | 4.679 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 9000 -36474.261 -36474.261 -36763.603 -36763.603 279.84112 279.84112 160408.42 160408.42 359.98864 359.98864 10000 -36469.104 -36469.104 -36755.368 -36755.368 276.86466 276.86466 160402.68 160402.68 425.51235 425.51235 Loop time of 167.438 on 1 procs for 1000 steps with 8000 atoms Performance: 0.516 ns/day, 46.511 hours/ns, 5.972 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 | 165.35 | 165.35 | 165.35 | 0.0 | 98.75 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.24131 | 0.24131 | 0.24131 | 0.0 | 0.14 Output | 0.00023015 | 0.00023015 | 0.00023015 | 0.0 | 0.00 Modify | 1.6603 | 1.6603 | 1.6603 | 0.0 | 0.99 Other | | 0.1869 | | | 0.11 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: 371798.0 ave 371798 max 371798 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 371798 Ave neighs/atom = 46.474750 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.168157458971, Press = 1.26499848531382 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 = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 17 17 17 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.679 | 4.679 | 4.679 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 10000 -36469.104 -36469.104 -36755.368 -36755.368 276.86466 276.86466 160402.68 160402.68 425.51235 425.51235 11000 -36470.02 -36470.02 -36759.176 -36759.176 279.66117 279.66117 160583.38 160583.38 -625.70201 -625.70201 Loop time of 167.291 on 1 procs for 1000 steps with 8000 atoms Performance: 0.516 ns/day, 46.470 hours/ns, 5.978 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 | 165.22 | 165.22 | 165.22 | 0.0 | 98.76 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.23831 | 0.23831 | 0.23831 | 0.0 | 0.14 Output | 0.00028289 | 0.00028289 | 0.00028289 | 0.0 | 0.00 Modify | 1.651 | 1.651 | 1.651 | 0.0 | 0.99 Other | | 0.1864 | | | 0.11 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: 372326.0 ave 372326 max 372326 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 372326 Ave neighs/atom = 46.540750 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.203575428609, Press = -1.62999906949034 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 = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 17 17 17 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.679 | 4.679 | 4.679 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 11000 -36470.02 -36470.02 -36759.176 -36759.176 279.66117 279.66117 160583.38 160583.38 -625.70201 -625.70201 12000 -36469.572 -36469.572 -36760.808 -36760.808 281.6723 281.6723 160492.91 160492.91 -81.349965 -81.349965 Loop time of 167.048 on 1 procs for 1000 steps with 8000 atoms Performance: 0.517 ns/day, 46.402 hours/ns, 5.986 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 | 164.98 | 164.98 | 164.98 | 0.0 | 98.76 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.23774 | 0.23774 | 0.23774 | 0.0 | 0.14 Output | 0.00030341 | 0.00030341 | 0.00030341 | 0.0 | 0.00 Modify | 1.6477 | 1.6477 | 1.6477 | 0.0 | 0.99 Other | | 0.1833 | | | 0.11 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: 371840.0 ave 371840 max 371840 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 371840 Ave neighs/atom = 46.480000 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.416103791931, Press = -3.40477260166568 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 = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 17 17 17 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.679 | 4.679 | 4.679 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 12000 -36469.572 -36469.572 -36760.808 -36760.808 281.6723 281.6723 160492.91 160492.91 -81.349965 -81.349965 13000 -36473.073 -36473.073 -36758.405 -36758.405 275.96273 275.96273 160439.4 160439.4 210.37172 210.37172 Loop time of 146.558 on 1 procs for 1000 steps with 8000 atoms Performance: 0.590 ns/day, 40.711 hours/ns, 6.823 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 | 144.82 | 144.82 | 144.82 | 0.0 | 98.81 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.21057 | 0.21057 | 0.21057 | 0.0 | 0.14 Output | 0.00023469 | 0.00023469 | 0.00023469 | 0.0 | 0.00 Modify | 1.365 | 1.365 | 1.365 | 0.0 | 0.93 Other | | 0.1631 | | | 0.11 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: 371776.0 ave 371776 max 371776 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 371776 Ave neighs/atom = 46.472000 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.359287231743, Press = -2.3164254075614 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 = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 17 17 17 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.679 | 4.679 | 4.679 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 13000 -36473.073 -36473.073 -36758.405 -36758.405 275.96273 275.96273 160439.4 160439.4 210.37172 210.37172 14000 -36473.225 -36473.225 -36759.128 -36759.128 276.51473 276.51473 160417.1 160417.1 287.94942 287.94942 Loop time of 167.755 on 1 procs for 1000 steps with 8000 atoms Performance: 0.515 ns/day, 46.599 hours/ns, 5.961 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 | 165.67 | 165.67 | 165.67 | 0.0 | 98.76 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.23677 | 0.23677 | 0.23677 | 0.0 | 0.14 Output | 0.0002329 | 0.0002329 | 0.0002329 | 0.0 | 0.00 Modify | 1.6623 | 1.6623 | 1.6623 | 0.0 | 0.99 Other | | 0.1829 | | | 0.11 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: 372090.0 ave 372090 max 372090 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 372090 Ave neighs/atom = 46.511250 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.294818219128, Press = -1.70493104067534 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 = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 17 17 17 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.679 | 4.679 | 4.679 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 14000 -36473.225 -36473.225 -36759.128 -36759.128 276.51473 276.51473 160417.1 160417.1 287.94942 287.94942 15000 -36465.71 -36465.71 -36753.375 -36753.375 278.21891 278.21891 160457.77 160457.77 43.901863 43.901863 Loop time of 159.471 on 1 procs for 1000 steps with 8000 atoms Performance: 0.542 ns/day, 44.298 hours/ns, 6.271 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 | 157.53 | 157.53 | 157.53 | 0.0 | 98.78 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.22768 | 0.22768 | 0.22768 | 0.0 | 0.14 Output | 0.00053786 | 0.00053786 | 0.00053786 | 0.0 | 0.00 Modify | 1.5427 | 1.5427 | 1.5427 | 0.0 | 0.97 Other | | 0.175 | | | 0.11 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: 372134.0 ave 372134 max 372134 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 372134 Ave neighs/atom = 46.516750 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.27156250036, Press = 0.311690625286921 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 = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 17 17 17 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.679 | 4.679 | 4.679 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 15000 -36465.71 -36465.71 -36753.375 -36753.375 278.21891 278.21891 160457.77 160457.77 43.901863 43.901863 16000 -36468.047 -36468.047 -36749.591 -36749.591 272.29891 272.29891 160646.18 160646.18 -1122.19 -1122.19 Loop time of 165.782 on 1 procs for 1000 steps with 8000 atoms Performance: 0.521 ns/day, 46.051 hours/ns, 6.032 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 | 163.74 | 163.74 | 163.74 | 0.0 | 98.77 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.23335 | 0.23335 | 0.23335 | 0.0 | 0.14 Output | 0.00023022 | 0.00023022 | 0.00023022 | 0.0 | 0.00 Modify | 1.6268 | 1.6268 | 1.6268 | 0.0 | 0.98 Other | | 0.1828 | | | 0.11 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: 372366.0 ave 372366 max 372366 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 372366 Ave neighs/atom = 46.545750 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.266388913338, Press = -0.477427501375859 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 = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 17 17 17 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.679 | 4.679 | 4.679 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 16000 -36468.047 -36468.047 -36749.591 -36749.591 272.29891 272.29891 160646.18 160646.18 -1122.19 -1122.19 17000 -36465.08 -36465.08 -36750.417 -36750.417 275.96762 275.96762 160497.82 160497.82 -167.50071 -167.50071 Loop time of 167.195 on 1 procs for 1000 steps with 8000 atoms Performance: 0.517 ns/day, 46.443 hours/ns, 5.981 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 | 165.13 | 165.13 | 165.13 | 0.0 | 98.76 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.23475 | 0.23475 | 0.23475 | 0.0 | 0.14 Output | 0.00022858 | 0.00022858 | 0.00022858 | 0.0 | 0.00 Modify | 1.6499 | 1.6499 | 1.6499 | 0.0 | 0.99 Other | | 0.1849 | | | 0.11 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: 372392.0 ave 372392 max 372392 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 372392 Ave neighs/atom = 46.549000 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.273753863815, Press = -2.68429248490274 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 = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 17 17 17 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.679 | 4.679 | 4.679 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 17000 -36465.08 -36465.08 -36750.417 -36750.417 275.96762 275.96762 160497.82 160497.82 -167.50071 -167.50071 18000 -36473.748 -36473.748 -36752.702 -36752.702 269.79393 269.79393 160359.86 160359.86 563.96341 563.96341 Loop time of 167.221 on 1 procs for 1000 steps with 8000 atoms Performance: 0.517 ns/day, 46.450 hours/ns, 5.980 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 | 165.15 | 165.15 | 165.15 | 0.0 | 98.76 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.23693 | 0.23693 | 0.23693 | 0.0 | 0.14 Output | 0.00022923 | 0.00022923 | 0.00022923 | 0.0 | 0.00 Modify | 1.6528 | 1.6528 | 1.6528 | 0.0 | 0.99 Other | | 0.1853 | | | 0.11 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: 372498.0 ave 372498 max 372498 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 372498 Ave neighs/atom = 46.562250 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.333032367848, Press = -0.902515829123254 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 = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 17 17 17 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.679 | 4.679 | 4.679 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 18000 -36473.748 -36473.748 -36752.702 -36752.702 269.79393 269.79393 160359.86 160359.86 563.96341 563.96341 19000 -36469.567 -36469.567 -36754.017 -36754.017 275.10985 275.10985 160489.35 160489.35 -149.02099 -149.02099 Loop time of 167.262 on 1 procs for 1000 steps with 8000 atoms Performance: 0.517 ns/day, 46.462 hours/ns, 5.979 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 | 165.2 | 165.2 | 165.2 | 0.0 | 98.76 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.23611 | 0.23611 | 0.23611 | 0.0 | 0.14 Output | 0.00031244 | 0.00031244 | 0.00031244 | 0.0 | 0.00 Modify | 1.6455 | 1.6455 | 1.6455 | 0.0 | 0.98 Other | | 0.1839 | | | 0.11 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: 372916.0 ave 372916 max 372916 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 372916 Ave neighs/atom = 46.614500 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.351483373069, Press = 0.447216557526911 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 = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 17 17 17 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.679 | 4.679 | 4.679 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 19000 -36469.567 -36469.567 -36754.017 -36754.017 275.10985 275.10985 160489.35 160489.35 -149.02099 -149.02099 20000 -36477.957 -36477.957 -36757.366 -36757.366 270.23394 270.23394 160543.06 160543.06 -529.79864 -529.79864 Loop time of 166.606 on 1 procs for 1000 steps with 8000 atoms Performance: 0.519 ns/day, 46.279 hours/ns, 6.002 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 | 164.54 | 164.54 | 164.54 | 0.0 | 98.76 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.23644 | 0.23644 | 0.23644 | 0.0 | 0.14 Output | 0.00033558 | 0.00033558 | 0.00033558 | 0.0 | 0.00 Modify | 1.6453 | 1.6453 | 1.6453 | 0.0 | 0.99 Other | | 0.1842 | | | 0.11 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: 372464.0 ave 372464 max 372464 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 372464 Ave neighs/atom = 46.558000 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.263225201601, Press = -1.44259928355111 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 = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 17 17 17 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.679 | 4.679 | 4.679 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 20000 -36477.957 -36477.957 -36757.366 -36757.366 270.23394 270.23394 160543.06 160543.06 -529.79864 -529.79864 21000 -36472.452 -36472.452 -36752.474 -36752.474 270.82703 270.82703 160416.78 160416.78 264.75203 264.75203 Loop time of 150.011 on 1 procs for 1000 steps with 8000 atoms Performance: 0.576 ns/day, 41.670 hours/ns, 6.666 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 | 148.19 | 148.19 | 148.19 | 0.0 | 98.79 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.22204 | 0.22204 | 0.22204 | 0.0 | 0.15 Output | 0.00023158 | 0.00023158 | 0.00023158 | 0.0 | 0.00 Modify | 1.4268 | 1.4268 | 1.4268 | 0.0 | 0.95 Other | | 0.167 | | | 0.11 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: 372062.0 ave 372062 max 372062 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 372062 Ave neighs/atom = 46.507750 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.254113063504, Press = -1.24080364172884 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 = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 17 17 17 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.679 | 4.679 | 4.679 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 21000 -36472.452 -36472.452 -36752.474 -36752.474 270.82703 270.82703 160416.78 160416.78 264.75203 264.75203 22000 -36470.927 -36470.927 -36757.028 -36757.028 276.70596 276.70596 160439.58 160439.58 191.19636 191.19636 Loop time of 164.952 on 1 procs for 1000 steps with 8000 atoms Performance: 0.524 ns/day, 45.820 hours/ns, 6.062 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 | 162.9 | 162.9 | 162.9 | 0.0 | 98.76 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.23686 | 0.23686 | 0.23686 | 0.0 | 0.14 Output | 0.00023104 | 0.00023104 | 0.00023104 | 0.0 | 0.00 Modify | 1.6281 | 1.6281 | 1.6281 | 0.0 | 0.99 Other | | 0.1836 | | | 0.11 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: 372486.0 ave 372486 max 372486 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 372486 Ave neighs/atom = 46.560750 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.238769999839, Press = -0.425017478671718 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 = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 17 17 17 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.679 | 4.679 | 4.679 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 22000 -36470.927 -36470.927 -36757.028 -36757.028 276.70596 276.70596 160439.58 160439.58 191.19636 191.19636 23000 -36471.162 -36471.162 -36753.526 -36753.526 273.09139 273.09139 160496.7 160496.7 -220.95926 -220.95926 Loop time of 157.548 on 1 procs for 1000 steps with 8000 atoms Performance: 0.548 ns/day, 43.763 hours/ns, 6.347 timesteps/s 99.9% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 155.63 | 155.63 | 155.63 | 0.0 | 98.78 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.22511 | 0.22511 | 0.22511 | 0.0 | 0.14 Output | 0.0002321 | 0.0002321 | 0.0002321 | 0.0 | 0.00 Modify | 1.5205 | 1.5205 | 1.5205 | 0.0 | 0.97 Other | | 0.174 | | | 0.11 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: 372156.0 ave 372156 max 372156 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 372156 Ave neighs/atom = 46.519500 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.24114252559, Press = -0.717333736951136 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 = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 17 17 17 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.679 | 4.679 | 4.679 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 23000 -36471.162 -36471.162 -36753.526 -36753.526 273.09139 273.09139 160496.7 160496.7 -220.95926 -220.95926 24000 -36464.136 -36464.136 -36754.027 -36754.027 280.37179 280.37179 160437.98 160437.98 213.7946 213.7946 Loop time of 163.656 on 1 procs for 1000 steps with 8000 atoms Performance: 0.528 ns/day, 45.460 hours/ns, 6.110 timesteps/s 99.9% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 161.63 | 161.63 | 161.63 | 0.0 | 98.76 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.23337 | 0.23337 | 0.23337 | 0.0 | 0.14 Output | 0.00023069 | 0.00023069 | 0.00023069 | 0.0 | 0.00 Modify | 1.6147 | 1.6147 | 1.6147 | 0.0 | 0.99 Other | | 0.1818 | | | 0.11 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: 372344.0 ave 372344 max 372344 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 372344 Ave neighs/atom = 46.543000 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.226444397419, Press = -0.799025395360884 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 = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 17 17 17 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.679 | 4.679 | 4.679 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 24000 -36464.136 -36464.136 -36754.027 -36754.027 280.37179 280.37179 160437.98 160437.98 213.7946 213.7946 25000 -36472.004 -36472.004 -36755.345 -36755.345 274.03641 274.03641 160502 160502 -223.29722 -223.29722 Loop time of 165.322 on 1 procs for 1000 steps with 8000 atoms Performance: 0.523 ns/day, 45.923 hours/ns, 6.049 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 | 163.28 | 163.28 | 163.28 | 0.0 | 98.76 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.23678 | 0.23678 | 0.23678 | 0.0 | 0.14 Output | 0.00023183 | 0.00023183 | 0.00023183 | 0.0 | 0.00 Modify | 1.6239 | 1.6239 | 1.6239 | 0.0 | 0.98 Other | | 0.1841 | | | 0.11 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: 372148.0 ave 372148 max 372148 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 372148 Ave neighs/atom = 46.518500 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.260379618664, Press = -0.752546793898843 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 = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 17 17 17 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.679 | 4.679 | 4.679 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 25000 -36472.004 -36472.004 -36755.345 -36755.345 274.03641 274.03641 160502 160502 -223.29722 -223.29722 26000 -36468.853 -36468.853 -36752.072 -36752.072 273.91863 273.91863 160512.75 160512.75 -340.69718 -340.69718 Loop time of 144.3 on 1 procs for 1000 steps with 8000 atoms Performance: 0.599 ns/day, 40.083 hours/ns, 6.930 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 | 142.59 | 142.59 | 142.59 | 0.0 | 98.82 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.21061 | 0.21061 | 0.21061 | 0.0 | 0.15 Output | 0.00023454 | 0.00023454 | 0.00023454 | 0.0 | 0.00 Modify | 1.3362 | 1.3362 | 1.3362 | 0.0 | 0.93 Other | | 0.1627 | | | 0.11 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: 372110.0 ave 372110 max 372110 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 372110 Ave neighs/atom = 46.513750 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.247446717758, Press = -1.55829953004814 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 = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 17 17 17 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.679 | 4.679 | 4.679 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 26000 -36468.853 -36468.853 -36752.072 -36752.072 273.91863 273.91863 160512.75 160512.75 -340.69718 -340.69718 27000 -36475.522 -36475.522 -36750.956 -36750.956 266.38897 266.38897 160312.49 160312.49 837.28417 837.28417 Loop time of 142.694 on 1 procs for 1000 steps with 8000 atoms Performance: 0.605 ns/day, 39.637 hours/ns, 7.008 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 | 141.01 | 141.01 | 141.01 | 0.0 | 98.82 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.20568 | 0.20568 | 0.20568 | 0.0 | 0.14 Output | 0.0002271 | 0.0002271 | 0.0002271 | 0.0 | 0.00 Modify | 1.3202 | 1.3202 | 1.3202 | 0.0 | 0.93 Other | | 0.1614 | | | 0.11 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: 372578.0 ave 372578 max 372578 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 372578 Ave neighs/atom = 46.572250 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.186022359776, Press = -0.76671317242318 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 = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 17 17 17 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.679 | 4.679 | 4.679 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 27000 -36475.522 -36475.522 -36750.956 -36750.956 266.38897 266.38897 160312.49 160312.49 837.28417 837.28417 28000 -36472.431 -36472.431 -36754.238 -36754.238 272.5535 272.5535 160442.94 160442.94 37.368708 37.368708 Loop time of 143.493 on 1 procs for 1000 steps with 8000 atoms Performance: 0.602 ns/day, 39.859 hours/ns, 6.969 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 | 141.79 | 141.79 | 141.79 | 0.0 | 98.81 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.21028 | 0.21028 | 0.21028 | 0.0 | 0.15 Output | 0.00023351 | 0.00023351 | 0.00023351 | 0.0 | 0.00 Modify | 1.3353 | 1.3353 | 1.3353 | 0.0 | 0.93 Other | | 0.1611 | | | 0.11 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: 373046.0 ave 373046 max 373046 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 373046 Ave neighs/atom = 46.630750 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 160459.003580683 print "LAMMPS calculation completed" LAMMPS calculation completed quit 0