# 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.429999992251397*${_u_distance} variable latticeconst_converted equal 5.429999992251397*1 lattice diamond ${latticeconst_converted} lattice diamond 5.4299999922514 Lattice spacing in x,y,z = 5.43 5.43 5.43 region simbox block 0 10 0 10 0 10 units lattice create_box 1 simbox Created orthogonal box = (0 0 0) to (54.3 54.3 54.3) 1 by 1 by 1 MPI processor grid create_atoms 1 box Created 8000 atoms using lattice units in orthogonal box = (0 0 0) to (54.3 54.3 54.3) 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 MEAM_LAMMPS_AslamBaskesDickel_2019_FeMnSiC__MO_427873955970_002 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 160103.0063146 variable V0_metal equal ${V0}/(${_u_distance}*${_u_distance}*${_u_distance}) variable V0_metal equal 160103.0063146/(${_u_distance}*${_u_distance}*${_u_distance}) variable V0_metal equal 160103.0063146/(1*${_u_distance}*${_u_distance}) variable V0_metal equal 160103.0063146/(1*1*${_u_distance}) variable V0_metal equal 160103.0063146/(1*1*1) variable V0_metal_times1000 equal ${V0_metal}*1000 variable V0_metal_times1000 equal 160103.0063146*1000 print "Initial system volume: ${V0_metal} Angstroms^3" Initial system volume: 160103.0063146 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 253.15*${_u_temperature} variable temp_converted equal 253.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 253.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 253.15 ${temp_converted} ${Tdamp_converted} iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 253.15 253.15 ${Tdamp_converted} iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 253.15 253.15 0.1 iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 253.15 253.15 0.1 iso 0 ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 253.15 253.15 0.1 iso 0 0 ${Pdamp_converted} fix ensemble all npt temp 253.15 253.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 "253.15 - 0.2" variable T_up equal "253.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_427873955970_002#item-citation CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 16 16 16 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 7.894 | 7.894 | 7.894 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 0 -36778.255 -36778.255 -37039.999 -37039.999 253.15 253.15 160103.01 160103.01 1746.2118 1746.2118 1000 -36482.93 -36482.93 -36756.284 -36756.284 264.37743 264.37743 162562.01 162562.01 -622.8048 -622.8048 Loop time of 30.6408 on 1 procs for 1000 steps with 8000 atoms Performance: 2.820 ns/day, 8.511 hours/ns, 32.636 timesteps/s 99.5% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 30.341 | 30.341 | 30.341 | 0.0 | 99.02 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.03497 | 0.03497 | 0.03497 | 0.0 | 0.11 Output | 0.00024289 | 0.00024289 | 0.00024289 | 0.0 | 0.00 Modify | 0.24199 | 0.24199 | 0.24199 | 0.0 | 0.79 Other | | 0.02291 | | | 0.07 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8601 ave 8601 max 8601 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 560000 ave 560000 max 560000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 560000 Ave neighs/atom = 70 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 0, Press = 0 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 16 16 16 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 8.146 | 8.146 | 8.146 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 1000 -36482.93 -36482.93 -36756.284 -36756.284 264.37743 264.37743 162562.01 162562.01 -622.8048 -622.8048 2000 -36511.647 -36511.647 -36779.387 -36779.387 258.94772 258.94772 162361.59 162361.59 -318.66855 -318.66855 Loop time of 41.1551 on 1 procs for 1000 steps with 8000 atoms Performance: 2.099 ns/day, 11.432 hours/ns, 24.298 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 | 40.863 | 40.863 | 40.863 | 0.0 | 99.29 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.033064 | 0.033064 | 0.033064 | 0.0 | 0.08 Output | 0.00011244 | 0.00011244 | 0.00011244 | 0.0 | 0.00 Modify | 0.23772 | 0.23772 | 0.23772 | 0.0 | 0.58 Other | | 0.02156 | | | 0.05 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8601 ave 8601 max 8601 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 572516 ave 572516 max 572516 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 572516 Ave neighs/atom = 71.5645 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 0, Press = 0 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 16 16 16 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 8.146 | 8.146 | 8.146 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 2000 -36511.647 -36511.647 -36779.387 -36779.387 258.94772 258.94772 162361.59 162361.59 -318.66855 -318.66855 3000 -36493.651 -36493.651 -36760.3 -36760.3 257.8925 257.8925 162462.96 162462.96 -214.78674 -214.78674 Loop time of 41.2873 on 1 procs for 1000 steps with 8000 atoms Performance: 2.093 ns/day, 11.469 hours/ns, 24.221 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 | 40.995 | 40.995 | 40.995 | 0.0 | 99.29 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.033151 | 0.033151 | 0.033151 | 0.0 | 0.08 Output | 0.00013311 | 0.00013311 | 0.00013311 | 0.0 | 0.00 Modify | 0.23758 | 0.23758 | 0.23758 | 0.0 | 0.58 Other | | 0.02165 | | | 0.05 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8601 ave 8601 max 8601 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 571894 ave 571894 max 571894 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 571894 Ave neighs/atom = 71.48675 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 0, Press = 0 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 16 16 16 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 8.146 | 8.146 | 8.146 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 3000 -36493.651 -36493.651 -36760.3 -36760.3 257.8925 257.8925 162462.96 162462.96 -214.78674 -214.78674 4000 -36507.956 -36507.956 -36767.65 -36767.65 251.16557 251.16557 162299.41 162299.41 370.82187 370.82187 Loop time of 41.3973 on 1 procs for 1000 steps with 8000 atoms Performance: 2.087 ns/day, 11.499 hours/ns, 24.156 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 | 41.104 | 41.104 | 41.104 | 0.0 | 99.29 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.033275 | 0.033275 | 0.033275 | 0.0 | 0.08 Output | 0.00011143 | 0.00011143 | 0.00011143 | 0.0 | 0.00 Modify | 0.23809 | 0.23809 | 0.23809 | 0.0 | 0.58 Other | | 0.02161 | | | 0.05 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8599 ave 8599 max 8599 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 572632 ave 572632 max 572632 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 572632 Ave neighs/atom = 71.579 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 0, Press = 0 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 16 16 16 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 8.146 | 8.146 | 8.146 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 4000 -36507.956 -36507.956 -36767.65 -36767.65 251.16557 251.16557 162299.41 162299.41 370.82187 370.82187 5000 -36497.963 -36497.963 -36766.025 -36766.025 259.25953 259.25953 162408.85 162408.85 -115.86022 -115.86022 Loop time of 41.3873 on 1 procs for 1000 steps with 8000 atoms Performance: 2.088 ns/day, 11.496 hours/ns, 24.162 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 | 41.094 | 41.094 | 41.094 | 0.0 | 99.29 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.033249 | 0.033249 | 0.033249 | 0.0 | 0.08 Output | 5.3801e-05 | 5.3801e-05 | 5.3801e-05 | 0.0 | 0.00 Modify | 0.23867 | 0.23867 | 0.23867 | 0.0 | 0.58 Other | | 0.02171 | | | 0.05 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8599 ave 8599 max 8599 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 573462 ave 573462 max 573462 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 573462 Ave neighs/atom = 71.68275 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 = 255.747141570106, Press = 116.854831322096 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 16 16 16 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 8.146 | 8.146 | 8.146 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 5000 -36497.963 -36497.963 -36766.025 -36766.025 259.25953 259.25953 162408.85 162408.85 -115.86022 -115.86022 6000 -36504.859 -36504.859 -36773.321 -36773.321 259.64631 259.64631 162431.86 162431.86 -544.80955 -544.80955 Loop time of 41.4463 on 1 procs for 1000 steps with 8000 atoms Performance: 2.085 ns/day, 11.513 hours/ns, 24.128 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 | 41.145 | 41.145 | 41.145 | 0.0 | 99.27 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.033483 | 0.033483 | 0.033483 | 0.0 | 0.08 Output | 8.1814e-05 | 8.1814e-05 | 8.1814e-05 | 0.0 | 0.00 Modify | 0.24579 | 0.24579 | 0.24579 | 0.0 | 0.59 Other | | 0.02171 | | | 0.05 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8601 ave 8601 max 8601 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 573046 ave 573046 max 573046 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 573046 Ave neighs/atom = 71.63075 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 = 253.049200774923, Press = -8.16446784346203 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 16 16 16 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 8.146 | 8.146 | 8.146 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 6000 -36504.859 -36504.859 -36773.321 -36773.321 259.64631 259.64631 162431.86 162431.86 -544.80955 -544.80955 7000 -36501.306 -36501.306 -36762.474 -36762.474 252.59226 252.59226 162387.68 162387.68 97.92379 97.92379 Loop time of 41.3512 on 1 procs for 1000 steps with 8000 atoms Performance: 2.089 ns/day, 11.486 hours/ns, 24.183 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 | 41.05 | 41.05 | 41.05 | 0.0 | 99.27 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.033543 | 0.033543 | 0.033543 | 0.0 | 0.08 Output | 8.2204e-05 | 8.2204e-05 | 8.2204e-05 | 0.0 | 0.00 Modify | 0.24564 | 0.24564 | 0.24564 | 0.0 | 0.59 Other | | 0.02173 | | | 0.05 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8600 ave 8600 max 8600 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 572636 ave 572636 max 572636 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 572636 Ave neighs/atom = 71.5795 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 = 253.299678722262, Press = 16.2028130003484 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 16 16 16 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 8.146 | 8.146 | 8.146 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 7000 -36501.306 -36501.306 -36762.474 -36762.474 252.59226 252.59226 162387.68 162387.68 97.92379 97.92379 8000 -36502.982 -36502.982 -36763.125 -36763.125 251.60021 251.60021 162436.41 162436.41 -234.96384 -234.96384 Loop time of 41.4474 on 1 procs for 1000 steps with 8000 atoms Performance: 2.085 ns/day, 11.513 hours/ns, 24.127 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 | 41.147 | 41.147 | 41.147 | 0.0 | 99.27 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.033409 | 0.033409 | 0.033409 | 0.0 | 0.08 Output | 5.2879e-05 | 5.2879e-05 | 5.2879e-05 | 0.0 | 0.00 Modify | 0.24537 | 0.24537 | 0.24537 | 0.0 | 0.59 Other | | 0.02167 | | | 0.05 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8598 ave 8598 max 8598 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 573772 ave 573772 max 573772 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 573772 Ave neighs/atom = 71.7215 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 = 253.110661586299, Press = -6.9345598378301 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 16 16 16 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 8.146 | 8.146 | 8.146 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 8000 -36502.982 -36502.982 -36763.125 -36763.125 251.60021 251.60021 162436.41 162436.41 -234.96384 -234.96384 9000 -36503.163 -36503.163 -36767.164 -36767.164 255.33236 255.33236 162364.4 162364.4 115.81839 115.81839 Loop time of 41.4334 on 1 procs for 1000 steps with 8000 atoms Performance: 2.085 ns/day, 11.509 hours/ns, 24.135 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 | 41.133 | 41.133 | 41.133 | 0.0 | 99.27 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.033446 | 0.033446 | 0.033446 | 0.0 | 0.08 Output | 8.3927e-05 | 8.3927e-05 | 8.3927e-05 | 0.0 | 0.00 Modify | 0.24529 | 0.24529 | 0.24529 | 0.0 | 0.59 Other | | 0.02168 | | | 0.05 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8601 ave 8601 max 8601 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 573118 ave 573118 max 573118 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 573118 Ave neighs/atom = 71.63975 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 = 253.209882075816, Press = 2.18093746210203 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 16 16 16 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 8.146 | 8.146 | 8.146 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 9000 -36503.163 -36503.163 -36767.164 -36767.164 255.33236 255.33236 162364.4 162364.4 115.81839 115.81839 10000 -36500.725 -36500.725 -36762.009 -36762.009 252.70482 252.70482 162408.54 162408.54 5.1915805 5.1915805 Loop time of 41.3295 on 1 procs for 1000 steps with 8000 atoms Performance: 2.091 ns/day, 11.480 hours/ns, 24.196 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 | 41.029 | 41.029 | 41.029 | 0.0 | 99.27 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.033352 | 0.033352 | 0.033352 | 0.0 | 0.08 Output | 8.2975e-05 | 8.2975e-05 | 8.2975e-05 | 0.0 | 0.00 Modify | 0.24539 | 0.24539 | 0.24539 | 0.0 | 0.59 Other | | 0.02169 | | | 0.05 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8601 ave 8601 max 8601 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 572374 ave 572374 max 572374 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 572374 Ave neighs/atom = 71.54675 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 = 253.134660444942, Press = 0.623050905660905 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 16 16 16 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 8.146 | 8.146 | 8.146 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 10000 -36500.725 -36500.725 -36762.009 -36762.009 252.70482 252.70482 162408.54 162408.54 5.1915805 5.1915805 11000 -36504.597 -36504.597 -36761.65 -36761.65 248.61268 248.61268 162374.11 162374.11 143.33757 143.33757 Loop time of 41.4516 on 1 procs for 1000 steps with 8000 atoms Performance: 2.084 ns/day, 11.514 hours/ns, 24.124 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 | 41.151 | 41.151 | 41.151 | 0.0 | 99.27 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.033501 | 0.033501 | 0.033501 | 0.0 | 0.08 Output | 5.5935e-05 | 5.5935e-05 | 5.5935e-05 | 0.0 | 0.00 Modify | 0.24577 | 0.24577 | 0.24577 | 0.0 | 0.59 Other | | 0.02165 | | | 0.05 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8600 ave 8600 max 8600 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 573386 ave 573386 max 573386 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 573386 Ave neighs/atom = 71.67325 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 = 253.18970931906, Press = 2.0447222812454 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 16 16 16 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 8.146 | 8.146 | 8.146 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 11000 -36504.597 -36504.597 -36761.65 -36761.65 248.61268 248.61268 162374.11 162374.11 143.33757 143.33757 12000 -36496.503 -36496.503 -36767.388 -36767.388 261.98937 261.98937 162536.74 162536.74 -884.67157 -884.67157 Loop time of 41.4947 on 1 procs for 1000 steps with 8000 atoms Performance: 2.082 ns/day, 11.526 hours/ns, 24.099 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 | 41.194 | 41.194 | 41.194 | 0.0 | 99.28 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.033333 | 0.033333 | 0.033333 | 0.0 | 0.08 Output | 9.035e-05 | 9.035e-05 | 9.035e-05 | 0.0 | 0.00 Modify | 0.24547 | 0.24547 | 0.24547 | 0.0 | 0.59 Other | | 0.02174 | | | 0.05 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8595 ave 8595 max 8595 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 573810 ave 573810 max 573810 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 573810 Ave neighs/atom = 71.72625 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 = 253.222655592517, Press = -0.28876762289915 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 16 16 16 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 8.146 | 8.146 | 8.146 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 12000 -36496.503 -36496.503 -36767.388 -36767.388 261.98937 261.98937 162536.74 162536.74 -884.67157 -884.67157 13000 -36503.446 -36503.446 -36761.109 -36761.109 249.20186 249.20186 162225.33 162225.33 1068.2608 1068.2608 Loop time of 41.3521 on 1 procs for 1000 steps with 8000 atoms Performance: 2.089 ns/day, 11.487 hours/ns, 24.183 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 | 41.051 | 41.051 | 41.051 | 0.0 | 99.27 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.033447 | 0.033447 | 0.033447 | 0.0 | 0.08 Output | 8.4439e-05 | 8.4439e-05 | 8.4439e-05 | 0.0 | 0.00 Modify | 0.24578 | 0.24578 | 0.24578 | 0.0 | 0.59 Other | | 0.02178 | | | 0.05 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8601 ave 8601 max 8601 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 571768 ave 571768 max 571768 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 571768 Ave neighs/atom = 71.471 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 = 253.312136758983, Press = 1.2493426712424 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 16 16 16 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 8.146 | 8.146 | 8.146 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 13000 -36503.446 -36503.446 -36761.109 -36761.109 249.20186 249.20186 162225.33 162225.33 1068.2608 1068.2608 14000 -36499.07 -36499.07 -36764.481 -36764.481 256.69648 256.69648 162565.99 162565.99 -1028.0107 -1028.0107 Loop time of 41.5168 on 1 procs for 1000 steps with 8000 atoms Performance: 2.081 ns/day, 11.532 hours/ns, 24.087 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 | 41.216 | 41.216 | 41.216 | 0.0 | 99.28 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.033477 | 0.033477 | 0.033477 | 0.0 | 0.08 Output | 0.00012482 | 0.00012482 | 0.00012482 | 0.0 | 0.00 Modify | 0.24512 | 0.24512 | 0.24512 | 0.0 | 0.59 Other | | 0.02171 | | | 0.05 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8599 ave 8599 max 8599 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 574908 ave 574908 max 574908 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 574908 Ave neighs/atom = 71.8635 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 = 253.374322934501, Press = -0.240226081602244 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 16 16 16 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 8.146 | 8.146 | 8.146 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 14000 -36499.07 -36499.07 -36764.481 -36764.481 256.69648 256.69648 162565.99 162565.99 -1028.0107 -1028.0107 15000 -36500.511 -36500.511 -36766.039 -36766.039 256.80879 256.80879 162248.8 162248.8 819.40362 819.40362 Loop time of 41.3954 on 1 procs for 1000 steps with 8000 atoms Performance: 2.087 ns/day, 11.499 hours/ns, 24.157 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 | 41.095 | 41.095 | 41.095 | 0.0 | 99.27 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.0334 | 0.0334 | 0.0334 | 0.0 | 0.08 Output | 5.6717e-05 | 5.6717e-05 | 5.6717e-05 | 0.0 | 0.00 Modify | 0.24549 | 0.24549 | 0.24549 | 0.0 | 0.59 Other | | 0.02168 | | | 0.05 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8601 ave 8601 max 8601 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 572758 ave 572758 max 572758 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 572758 Ave neighs/atom = 71.59475 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 = 253.32601311811, Press = 1.52644851256645 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 16 16 16 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 8.146 | 8.146 | 8.146 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 15000 -36500.511 -36500.511 -36766.039 -36766.039 256.80879 256.80879 162248.8 162248.8 819.40362 819.40362 16000 -36498.7 -36498.7 -36760.74 -36760.74 253.43576 253.43576 162489.23 162489.23 -415.32558 -415.32558 Loop time of 41.5591 on 1 procs for 1000 steps with 8000 atoms Performance: 2.079 ns/day, 11.544 hours/ns, 24.062 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 | 41.258 | 41.258 | 41.258 | 0.0 | 99.28 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.033448 | 0.033448 | 0.033448 | 0.0 | 0.08 Output | 8.4348e-05 | 8.4348e-05 | 8.4348e-05 | 0.0 | 0.00 Modify | 0.24603 | 0.24603 | 0.24603 | 0.0 | 0.59 Other | | 0.02172 | | | 0.05 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8601 ave 8601 max 8601 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 574048 ave 574048 max 574048 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 574048 Ave neighs/atom = 71.756 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 = 253.378825065805, Press = 0.294609815324774 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 16 16 16 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 8.146 | 8.146 | 8.146 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 16000 -36498.7 -36498.7 -36760.74 -36760.74 253.43576 253.43576 162489.23 162489.23 -415.32558 -415.32558 17000 -36502.182 -36502.182 -36762.3 -36762.3 251.57661 251.57661 162366.13 162366.13 243.51117 243.51117 Loop time of 41.5066 on 1 procs for 1000 steps with 8000 atoms Performance: 2.082 ns/day, 11.530 hours/ns, 24.093 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 | 41.205 | 41.205 | 41.205 | 0.0 | 99.27 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.033405 | 0.033405 | 0.033405 | 0.0 | 0.08 Output | 7.996e-05 | 7.996e-05 | 7.996e-05 | 0.0 | 0.00 Modify | 0.24601 | 0.24601 | 0.24601 | 0.0 | 0.59 Other | | 0.02172 | | | 0.05 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8601 ave 8601 max 8601 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 572510 ave 572510 max 572510 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 572510 Ave neighs/atom = 71.56375 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 = 253.360327964887, Press = -0.41034123585271 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 16 16 16 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 8.146 | 8.146 | 8.146 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 17000 -36502.182 -36502.182 -36762.3 -36762.3 251.57661 251.57661 162366.13 162366.13 243.51117 243.51117 18000 -36500.126 -36500.126 -36760.298 -36760.298 251.62891 251.62891 162462.28 162462.28 -250.01967 -250.01967 Loop time of 41.5163 on 1 procs for 1000 steps with 8000 atoms Performance: 2.081 ns/day, 11.532 hours/ns, 24.087 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 | 41.215 | 41.215 | 41.215 | 0.0 | 99.27 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.03344 | 0.03344 | 0.03344 | 0.0 | 0.08 Output | 5.2178e-05 | 5.2178e-05 | 5.2178e-05 | 0.0 | 0.00 Modify | 0.24604 | 0.24604 | 0.24604 | 0.0 | 0.59 Other | | 0.02164 | | | 0.05 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8600 ave 8600 max 8600 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 573670 ave 573670 max 573670 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 573670 Ave neighs/atom = 71.70875 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 = 253.392444841875, Press = 0.701374332975029 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 16 16 16 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 8.146 | 8.146 | 8.146 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 18000 -36500.126 -36500.126 -36760.298 -36760.298 251.62891 251.62891 162462.28 162462.28 -250.01967 -250.01967 19000 -36499.908 -36499.908 -36760.651 -36760.651 252.18099 252.18099 162251.33 162251.33 998.07787 998.07787 Loop time of 41.5078 on 1 procs for 1000 steps with 8000 atoms Performance: 2.082 ns/day, 11.530 hours/ns, 24.092 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 | 41.206 | 41.206 | 41.206 | 0.0 | 99.27 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.033657 | 0.033657 | 0.033657 | 0.0 | 0.08 Output | 8.5892e-05 | 8.5892e-05 | 8.5892e-05 | 0.0 | 0.00 Modify | 0.24608 | 0.24608 | 0.24608 | 0.0 | 0.59 Other | | 0.02169 | | | 0.05 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8600 ave 8600 max 8600 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 573200 ave 573200 max 573200 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 573200 Ave neighs/atom = 71.65 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 = 253.35456551276, Press = -0.89484174779178 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 16 16 16 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 8.146 | 8.146 | 8.146 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 19000 -36499.908 -36499.908 -36760.651 -36760.651 252.18099 252.18099 162251.33 162251.33 998.07787 998.07787 20000 -36504.817 -36504.817 -36767.443 -36767.443 254.00239 254.00239 162506.4 162506.4 -815.81669 -815.81669 Loop time of 41.6173 on 1 procs for 1000 steps with 8000 atoms Performance: 2.076 ns/day, 11.560 hours/ns, 24.028 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 | 41.316 | 41.316 | 41.316 | 0.0 | 99.28 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.033471 | 0.033471 | 0.033471 | 0.0 | 0.08 Output | 8.4779e-05 | 8.4779e-05 | 8.4779e-05 | 0.0 | 0.00 Modify | 0.24617 | 0.24617 | 0.24617 | 0.0 | 0.59 Other | | 0.02178 | | | 0.05 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8601 ave 8601 max 8601 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 574386 ave 574386 max 574386 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 574386 Ave neighs/atom = 71.79825 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 = 253.384960180788, Press = 1.86792340198623 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 16 16 16 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 8.146 | 8.146 | 8.146 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 20000 -36504.817 -36504.817 -36767.443 -36767.443 254.00239 254.00239 162506.4 162506.4 -815.81669 -815.81669 21000 -36501.039 -36501.039 -36759.269 -36759.269 249.75062 249.75062 162336.86 162336.86 508.8412 508.8412 Loop time of 41.6304 on 1 procs for 1000 steps with 8000 atoms Performance: 2.075 ns/day, 11.564 hours/ns, 24.021 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 | 41.329 | 41.329 | 41.329 | 0.0 | 99.27 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.033714 | 0.033714 | 0.033714 | 0.0 | 0.08 Output | 5.2229e-05 | 5.2229e-05 | 5.2229e-05 | 0.0 | 0.00 Modify | 0.24644 | 0.24644 | 0.24644 | 0.0 | 0.59 Other | | 0.02171 | | | 0.05 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8601 ave 8601 max 8601 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 572668 ave 572668 max 572668 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 572668 Ave neighs/atom = 71.5835 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 = 253.365047737815, Press = -0.109132121505246 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 16 16 16 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 8.146 | 8.146 | 8.146 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 21000 -36501.039 -36501.039 -36759.269 -36759.269 249.75062 249.75062 162336.86 162336.86 508.8412 508.8412 22000 -36504.929 -36504.929 -36761.216 -36761.216 247.87087 247.87087 162449.26 162449.26 -278.75559 -278.75559 Loop time of 41.7328 on 1 procs for 1000 steps with 8000 atoms Performance: 2.070 ns/day, 11.592 hours/ns, 23.962 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 | 41.431 | 41.431 | 41.431 | 0.0 | 99.28 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.033526 | 0.033526 | 0.033526 | 0.0 | 0.08 Output | 8.2505e-05 | 8.2505e-05 | 8.2505e-05 | 0.0 | 0.00 Modify | 0.24677 | 0.24677 | 0.24677 | 0.0 | 0.59 Other | | 0.02178 | | | 0.05 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8601 ave 8601 max 8601 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 573834 ave 573834 max 573834 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 573834 Ave neighs/atom = 71.72925 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 = 253.283150432458, Press = 0.956051475212002 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 16 16 16 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 8.146 | 8.146 | 8.146 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 22000 -36504.929 -36504.929 -36761.216 -36761.216 247.87087 247.87087 162449.26 162449.26 -278.75559 -278.75559 23000 -36501.16 -36501.16 -36761.949 -36761.949 252.22494 252.22494 162348.36 162348.36 393.15769 393.15769 Loop time of 41.8221 on 1 procs for 1000 steps with 8000 atoms Performance: 2.066 ns/day, 11.617 hours/ns, 23.911 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 | 41.518 | 41.518 | 41.518 | 0.0 | 99.27 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.033691 | 0.033691 | 0.033691 | 0.0 | 0.08 Output | 8.9889e-05 | 8.9889e-05 | 8.9889e-05 | 0.0 | 0.00 Modify | 0.24812 | 0.24812 | 0.24812 | 0.0 | 0.59 Other | | 0.02189 | | | 0.05 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8601 ave 8601 max 8601 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 573820 ave 573820 max 573820 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 573820 Ave neighs/atom = 71.7275 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 = 253.259798939192, Press = -0.396965587350203 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 16 16 16 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 8.146 | 8.146 | 8.146 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 23000 -36501.16 -36501.16 -36761.949 -36761.949 252.22494 252.22494 162348.36 162348.36 393.15769 393.15769 24000 -36502.716 -36502.716 -36762.321 -36762.321 251.08051 251.08051 162389.66 162389.66 107.65595 107.65595 Loop time of 41.8078 on 1 procs for 1000 steps with 8000 atoms Performance: 2.067 ns/day, 11.613 hours/ns, 23.919 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 | 41.505 | 41.505 | 41.505 | 0.0 | 99.27 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.033574 | 0.033574 | 0.033574 | 0.0 | 0.08 Output | 5.312e-05 | 5.312e-05 | 5.312e-05 | 0.0 | 0.00 Modify | 0.24774 | 0.24774 | 0.24774 | 0.0 | 0.59 Other | | 0.02186 | | | 0.05 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8601 ave 8601 max 8601 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 573364 ave 573364 max 573364 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 573364 Ave neighs/atom = 71.6705 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 = 253.247785841404, Press = 0.732546006061106 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 16 16 16 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 8.146 | 8.146 | 8.146 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 24000 -36502.716 -36502.716 -36762.321 -36762.321 251.08051 251.08051 162389.66 162389.66 107.65595 107.65595 25000 -36501.739 -36501.739 -36760.926 -36760.926 250.67573 250.67573 162415.9 162415.9 -32.593509 -32.593509 Loop time of 41.8601 on 1 procs for 1000 steps with 8000 atoms Performance: 2.064 ns/day, 11.628 hours/ns, 23.889 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 | 41.556 | 41.556 | 41.556 | 0.0 | 99.27 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.033839 | 0.033839 | 0.033839 | 0.0 | 0.08 Output | 8.5671e-05 | 8.5671e-05 | 8.5671e-05 | 0.0 | 0.00 Modify | 0.2481 | 0.2481 | 0.2481 | 0.0 | 0.59 Other | | 0.02186 | | | 0.05 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 8599 ave 8599 max 8599 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 573548 ave 573548 max 573548 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 573548 Ave neighs/atom = 71.6935 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_T253.15.out" else "print 'not_converged' file output/vol_T253.15.out" print '${V}' file output/vol_T253.15.out 162397.176129139 print "LAMMPS calculation completed" LAMMPS calculation completed quit 0