# 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 2.8809780403971668*${_u_distance} variable latticeconst_converted equal 2.8809780403971668*1 lattice bcc ${latticeconst_converted} lattice bcc 2.88097804039717 Lattice spacing in x,y,z = 2.880978 2.880978 2.880978 region simbox block 0 10 0 10 0 10 units lattice create_box 1 simbox Created orthogonal box = (0 0 0) to (28.80978 28.80978 28.80978) 1 by 1 by 1 MPI processor grid create_atoms 1 box Created 2000 atoms using lattice units in orthogonal box = (0 0 0) to (28.80978 28.80978 28.80978) create_atoms CPU = 0.002 seconds variable mass_converted equal 51.9961*${_u_mass} variable mass_converted equal 51.9961*1 kim_interactions Cr WARNING: 'kim_' has been renamed to 'kim '. Please update your input. kim interactions Cr #=== BEGIN kim interactions ================================== pair_style kim MEAM_LAMMPS_ChoiKimSeol_2017_CrMn__MO_671124822359_002 pair_coeff * * Cr #=== END kim interactions ==================================== mass 1 ${mass_converted} mass 1 51.9961 # initial volume variable v equal vol # assign formula variable V0 equal ${v} # evaluate initial value variable V0 equal 23912.2170404509 variable V0_metal equal ${V0}/(${_u_distance}*${_u_distance}*${_u_distance}) variable V0_metal equal 23912.2170404509/(${_u_distance}*${_u_distance}*${_u_distance}) variable V0_metal equal 23912.2170404509/(1*${_u_distance}*${_u_distance}) variable V0_metal equal 23912.2170404509/(1*1*${_u_distance}) variable V0_metal equal 23912.2170404509/(1*1*1) variable V0_metal_times1000 equal ${V0_metal}*1000 variable V0_metal_times1000 equal 23912.2170404509*1000 print "Initial system volume: ${V0_metal} Angstroms^3" Initial system volume: 23912.2170404509 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_671124822359_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 = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 9 9 9 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.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 0 -8134.5883 -8134.5883 -8200 -8200 253.15 253.15 23912.217 23912.217 2921.8042 2921.8042 1000 -8067.3727 -8067.3727 -8132.8821 -8132.8821 253.52832 253.52832 24084.253 24084.253 3378.7616 3378.7616 Loop time of 163.417 on 1 procs for 1000 steps with 2000 atoms Performance: 0.529 ns/day, 45.394 hours/ns, 6.119 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.86 | 162.86 | 162.86 | 0.0 | 99.66 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.10374 | 0.10374 | 0.10374 | 0.0 | 0.06 Output | 0.00028215 | 0.00028215 | 0.00028215 | 0.0 | 0.00 Modify | 0.38708 | 0.38708 | 0.38708 | 0.0 | 0.24 Other | | 0.06485 | | | 0.04 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 0, Press = 0 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 9 9 9 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.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 1000 -8067.3727 -8067.3727 -8132.8821 -8132.8821 253.52832 253.52832 24084.253 24084.253 3378.7616 3378.7616 2000 -8066.2642 -8066.2642 -8132.8578 -8132.8578 257.72427 257.72427 24108.423 24108.423 1400.3185 1400.3185 Loop time of 146.177 on 1 procs for 1000 steps with 2000 atoms Performance: 0.591 ns/day, 40.605 hours/ns, 6.841 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.69 | 145.69 | 145.69 | 0.0 | 99.66 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.095417 | 0.095417 | 0.095417 | 0.0 | 0.07 Output | 0.00020632 | 0.00020632 | 0.00020632 | 0.0 | 0.00 Modify | 0.3364 | 0.3364 | 0.3364 | 0.0 | 0.23 Other | | 0.05883 | | | 0.04 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 209794 ave 209794 max 209794 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 209794 Ave neighs/atom = 104.897 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 = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 9 9 9 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.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 2000 -8066.2642 -8066.2642 -8132.8578 -8132.8578 257.72427 257.72427 24108.423 24108.423 1400.3185 1400.3185 3000 -8069.2541 -8069.2541 -8130.7811 -8130.7811 238.11595 238.11595 24109.289 24109.289 1120.7114 1120.7114 Loop time of 161.454 on 1 procs for 1000 steps with 2000 atoms Performance: 0.535 ns/day, 44.848 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 | 160.9 | 160.9 | 160.9 | 0.0 | 99.65 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.10438 | 0.10438 | 0.10438 | 0.0 | 0.06 Output | 0.00020933 | 0.00020933 | 0.00020933 | 0.0 | 0.00 Modify | 0.38866 | 0.38866 | 0.38866 | 0.0 | 0.24 Other | | 0.06425 | | | 0.04 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 209432 ave 209432 max 209432 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 209432 Ave neighs/atom = 104.716 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 = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 9 9 9 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.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 3000 -8069.2541 -8069.2541 -8130.7811 -8130.7811 238.11595 238.11595 24109.289 24109.289 1120.7114 1120.7114 4000 -8065.397 -8065.397 -8134.1992 -8134.1992 266.2719 266.2719 24116.548 24116.548 612.18948 612.18948 Loop time of 146.279 on 1 procs for 1000 steps with 2000 atoms Performance: 0.591 ns/day, 40.633 hours/ns, 6.836 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.79 | 145.79 | 145.79 | 0.0 | 99.66 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.095964 | 0.095964 | 0.095964 | 0.0 | 0.07 Output | 0.00024738 | 0.00024738 | 0.00024738 | 0.0 | 0.00 Modify | 0.33707 | 0.33707 | 0.33707 | 0.0 | 0.23 Other | | 0.05942 | | | 0.04 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 209130 ave 209130 max 209130 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 209130 Ave neighs/atom = 104.565 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 = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 9 9 9 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.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 4000 -8065.397 -8065.397 -8134.1992 -8134.1992 266.2719 266.2719 24116.548 24116.548 612.18948 612.18948 5000 -8069.3265 -8069.3265 -8134.6988 -8134.6988 252.99741 252.99741 24107.699 24107.699 990.07035 990.07035 Loop time of 145.86 on 1 procs for 1000 steps with 2000 atoms Performance: 0.592 ns/day, 40.517 hours/ns, 6.856 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.37 | 145.37 | 145.37 | 0.0 | 99.66 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.096194 | 0.096194 | 0.096194 | 0.0 | 0.07 Output | 0.00017141 | 0.00017141 | 0.00017141 | 0.0 | 0.00 Modify | 0.33791 | 0.33791 | 0.33791 | 0.0 | 0.23 Other | | 0.0592 | | | 0.04 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 209376 ave 209376 max 209376 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 209376 Ave neighs/atom = 104.688 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 = 247.834812003668, Press = -394.758070578651 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 = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 9 9 9 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.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 5000 -8069.3265 -8069.3265 -8134.6988 -8134.6988 252.99741 252.99741 24107.699 24107.699 990.07035 990.07035 6000 -8066.372 -8066.372 -8133.542 -8133.542 259.95476 259.95476 24102.169 24102.169 1805.0495 1805.0495 Loop time of 156.546 on 1 procs for 1000 steps with 2000 atoms Performance: 0.552 ns/day, 43.485 hours/ns, 6.388 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 | 155.99 | 155.99 | 155.99 | 0.0 | 99.64 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.10338 | 0.10338 | 0.10338 | 0.0 | 0.07 Output | 0.00017101 | 0.00017101 | 0.00017101 | 0.0 | 0.00 Modify | 0.39371 | 0.39371 | 0.39371 | 0.0 | 0.25 Other | | 0.06201 | | | 0.04 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 209454 ave 209454 max 209454 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 209454 Ave neighs/atom = 104.727 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 = 252.338331136479, Press = -25.3225643954218 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 = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 9 9 9 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.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 6000 -8066.372 -8066.372 -8133.542 -8133.542 259.95476 259.95476 24102.169 24102.169 1805.0495 1805.0495 7000 -8066.866 -8066.866 -8135.6686 -8135.6686 266.27345 266.27345 24085.53 24085.53 2743.1175 2743.1175 Loop time of 147.608 on 1 procs for 1000 steps with 2000 atoms Performance: 0.585 ns/day, 41.002 hours/ns, 6.775 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 | 147.09 | 147.09 | 147.09 | 0.0 | 99.65 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.096719 | 0.096719 | 0.096719 | 0.0 | 0.07 Output | 0.00020505 | 0.00020505 | 0.00020505 | 0.0 | 0.00 Modify | 0.36393 | 0.36393 | 0.36393 | 0.0 | 0.25 Other | | 0.05905 | | | 0.04 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 209618 ave 209618 max 209618 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 209618 Ave neighs/atom = 104.809 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 = 251.663858629339, Press = -12.0244432424192 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 = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 9 9 9 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.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 7000 -8066.866 -8066.866 -8135.6686 -8135.6686 266.27345 266.27345 24085.53 24085.53 2743.1175 2743.1175 8000 -8068.4392 -8068.4392 -8135.0434 -8135.0434 257.76538 257.76538 24086.31 24086.31 2694.6508 2694.6508 Loop time of 146.002 on 1 procs for 1000 steps with 2000 atoms Performance: 0.592 ns/day, 40.556 hours/ns, 6.849 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.49 | 145.49 | 145.49 | 0.0 | 99.65 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.096543 | 0.096543 | 0.096543 | 0.0 | 0.07 Output | 0.00016842 | 0.00016842 | 0.00016842 | 0.0 | 0.00 Modify | 0.35617 | 0.35617 | 0.35617 | 0.0 | 0.24 Other | | 0.05881 | | | 0.04 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 210044 ave 210044 max 210044 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 210044 Ave neighs/atom = 105.022 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 = 252.236332577871, Press = -10.9305808759779 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 = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 9 9 9 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.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 8000 -8068.4392 -8068.4392 -8135.0434 -8135.0434 257.76538 257.76538 24086.31 24086.31 2694.6508 2694.6508 9000 -8066.891 -8066.891 -8131.4145 -8131.4145 249.71248 249.71248 24105.459 24105.459 1737.7051 1737.7051 Loop time of 145.821 on 1 procs for 1000 steps with 2000 atoms Performance: 0.593 ns/day, 40.506 hours/ns, 6.858 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.31 | 145.31 | 145.31 | 0.0 | 99.65 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.096066 | 0.096066 | 0.096066 | 0.0 | 0.07 Output | 0.00020016 | 0.00020016 | 0.00020016 | 0.0 | 0.00 Modify | 0.35667 | 0.35667 | 0.35667 | 0.0 | 0.24 Other | | 0.05894 | | | 0.04 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 210008 ave 210008 max 210008 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 210008 Ave neighs/atom = 105.004 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 = 252.447866773865, Press = -6.58353421345232 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 = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 9 9 9 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.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 9000 -8066.891 -8066.891 -8131.4145 -8131.4145 249.71248 249.71248 24105.459 24105.459 1737.7051 1737.7051 10000 -8068.8352 -8068.8352 -8135.6668 -8135.6668 258.64522 258.64522 24118.506 24118.506 -81.864599 -81.864599 Loop time of 145.676 on 1 procs for 1000 steps with 2000 atoms Performance: 0.593 ns/day, 40.466 hours/ns, 6.865 timesteps/s 100.1% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 145.16 | 145.16 | 145.16 | 0.0 | 99.65 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.096703 | 0.096703 | 0.096703 | 0.0 | 0.07 Output | 0.00016879 | 0.00016879 | 0.00016879 | 0.0 | 0.00 Modify | 0.35657 | 0.35657 | 0.35657 | 0.0 | 0.24 Other | | 0.05896 | | | 0.04 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 209142 ave 209142 max 209142 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 209142 Ave neighs/atom = 104.571 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 = 252.597996474633, Press = -4.93577695029161 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 = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 9 9 9 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.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 10000 -8068.8352 -8068.8352 -8135.6668 -8135.6668 258.64522 258.64522 24118.506 24118.506 -81.864599 -81.864599 11000 -8065.2532 -8065.2532 -8131.2369 -8131.2369 255.36386 255.36386 24120.665 24120.665 846.15671 846.15671 Loop time of 145.941 on 1 procs for 1000 steps with 2000 atoms Performance: 0.592 ns/day, 40.539 hours/ns, 6.852 timesteps/s 100.1% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 145.43 | 145.43 | 145.43 | 0.0 | 99.65 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.096217 | 0.096217 | 0.096217 | 0.0 | 0.07 Output | 0.0001701 | 0.0001701 | 0.0001701 | 0.0 | 0.00 Modify | 0.35718 | 0.35718 | 0.35718 | 0.0 | 0.24 Other | | 0.0583 | | | 0.04 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 209240 ave 209240 max 209240 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 209240 Ave neighs/atom = 104.62 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.026601781319, Press = -5.49155643568945 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 = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 9 9 9 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.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 11000 -8065.2532 -8065.2532 -8131.2369 -8131.2369 255.36386 255.36386 24120.665 24120.665 846.15671 846.15671 12000 -8066.9466 -8066.9466 -8132.8698 -8132.8698 255.12967 255.12967 24112.937 24112.937 977.05805 977.05805 Loop time of 145.602 on 1 procs for 1000 steps with 2000 atoms Performance: 0.593 ns/day, 40.445 hours/ns, 6.868 timesteps/s 100.1% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 145.09 | 145.09 | 145.09 | 0.0 | 99.65 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.096686 | 0.096686 | 0.096686 | 0.0 | 0.07 Output | 0.00016992 | 0.00016992 | 0.00016992 | 0.0 | 0.00 Modify | 0.35653 | 0.35653 | 0.35653 | 0.0 | 0.24 Other | | 0.05899 | | | 0.04 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 208824 ave 208824 max 208824 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 208824 Ave neighs/atom = 104.412 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.008120223404, Press = -4.69923655100586 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 = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 9 9 9 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.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 12000 -8066.9466 -8066.9466 -8132.8698 -8132.8698 255.12967 255.12967 24112.937 24112.937 977.05805 977.05805 13000 -8066.5979 -8066.5979 -8133.1287 -8133.1287 257.48083 257.48083 24113.95 24113.95 1039.3301 1039.3301 Loop time of 150.874 on 1 procs for 1000 steps with 2000 atoms Performance: 0.573 ns/day, 41.910 hours/ns, 6.628 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 | 150.34 | 150.34 | 150.34 | 0.0 | 99.64 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.10001 | 0.10001 | 0.10001 | 0.0 | 0.07 Output | 0.0002305 | 0.0002305 | 0.0002305 | 0.0 | 0.00 Modify | 0.37527 | 0.37527 | 0.37527 | 0.0 | 0.25 Other | | 0.06076 | | | 0.04 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 209344 ave 209344 max 209344 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 209344 Ave neighs/atom = 104.672 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 = 252.956198994266, Press = -4.22741217595175 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 = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 9 9 9 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.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 13000 -8066.5979 -8066.5979 -8133.1287 -8133.1287 257.48083 257.48083 24113.95 24113.95 1039.3301 1039.3301 14000 -8067.5677 -8067.5677 -8133.907 -8133.907 256.74013 256.74013 24097.329 24097.329 2164.8314 2164.8314 Loop time of 171.528 on 1 procs for 1000 steps with 2000 atoms Performance: 0.504 ns/day, 47.647 hours/ns, 5.830 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 | 170.91 | 170.91 | 170.91 | 0.0 | 99.64 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.10964 | 0.10964 | 0.10964 | 0.0 | 0.06 Output | 0.0002157 | 0.0002157 | 0.0002157 | 0.0 | 0.00 Modify | 0.44336 | 0.44336 | 0.44336 | 0.0 | 0.26 Other | | 0.06737 | | | 0.04 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 209348 ave 209348 max 209348 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 209348 Ave neighs/atom = 104.674 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 = 252.779391297363, Press = -2.23223458639935 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 = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 9 9 9 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.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 14000 -8067.5677 -8067.5677 -8133.907 -8133.907 256.74013 256.74013 24097.329 24097.329 2164.8314 2164.8314 15000 -8063.3899 -8063.3899 -8131.5788 -8131.5788 263.898 263.898 24096.967 24096.967 2751.4104 2751.4104 Loop time of 152.008 on 1 procs for 1000 steps with 2000 atoms Performance: 0.568 ns/day, 42.224 hours/ns, 6.579 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 | 151.47 | 151.47 | 151.47 | 0.0 | 99.65 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.098965 | 0.098965 | 0.098965 | 0.0 | 0.07 Output | 0.00021393 | 0.00021393 | 0.00021393 | 0.0 | 0.00 Modify | 0.37619 | 0.37619 | 0.37619 | 0.0 | 0.25 Other | | 0.06037 | | | 0.04 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 209422 ave 209422 max 209422 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 209422 Ave neighs/atom = 104.711 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 = 252.921553546056, Press = -0.195345343880037 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 = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 9 9 9 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.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 15000 -8063.3899 -8063.3899 -8131.5788 -8131.5788 263.898 263.898 24096.967 24096.967 2751.4104 2751.4104 16000 -8067.0421 -8067.0421 -8133.736 -8133.736 258.11204 258.11204 24104.144 24104.144 1619.9012 1619.9012 Loop time of 145.536 on 1 procs for 1000 steps with 2000 atoms Performance: 0.594 ns/day, 40.427 hours/ns, 6.871 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.02 | 145.02 | 145.02 | 0.0 | 99.65 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.09601 | 0.09601 | 0.09601 | 0.0 | 0.07 Output | 0.00016836 | 0.00016836 | 0.00016836 | 0.0 | 0.00 Modify | 0.35654 | 0.35654 | 0.35654 | 0.0 | 0.24 Other | | 0.05847 | | | 0.04 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 209432 ave 209432 max 209432 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 209432 Ave neighs/atom = 104.716 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.02044934094, Press = -0.896968463327436 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 = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 9 9 9 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.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 16000 -8067.0421 -8067.0421 -8133.736 -8133.736 258.11204 258.11204 24104.144 24104.144 1619.9012 1619.9012 17000 -8067.6246 -8067.6246 -8132.9055 -8132.9055 252.64361 252.64361 24108.508 24108.508 1243.7541 1243.7541 Loop time of 145.922 on 1 procs for 1000 steps with 2000 atoms Performance: 0.592 ns/day, 40.534 hours/ns, 6.853 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.41 | 145.41 | 145.41 | 0.0 | 99.65 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.09651 | 0.09651 | 0.09651 | 0.0 | 0.07 Output | 0.00020245 | 0.00020245 | 0.00020245 | 0.0 | 0.00 Modify | 0.35605 | 0.35605 | 0.35605 | 0.0 | 0.24 Other | | 0.05805 | | | 0.04 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 209434 ave 209434 max 209434 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 209434 Ave neighs/atom = 104.717 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.066391206483, Press = -0.780711061853579 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 = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 9 9 9 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.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 17000 -8067.6246 -8067.6246 -8132.9055 -8132.9055 252.64361 252.64361 24108.508 24108.508 1243.7541 1243.7541 18000 -8068.744 -8068.744 -8134.7203 -8134.7203 255.3351 255.3351 24097.044 24097.044 1832.5091 1832.5091 Loop time of 145.998 on 1 procs for 1000 steps with 2000 atoms Performance: 0.592 ns/day, 40.555 hours/ns, 6.849 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.49 | 145.49 | 145.49 | 0.0 | 99.65 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.097005 | 0.097005 | 0.097005 | 0.0 | 0.07 Output | 0.00016726 | 0.00016726 | 0.00016726 | 0.0 | 0.00 Modify | 0.35469 | 0.35469 | 0.35469 | 0.0 | 0.24 Other | | 0.0594 | | | 0.04 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 209412 ave 209412 max 209412 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 209412 Ave neighs/atom = 104.706 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.156839182774, Press = -0.733265380225339 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 = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 9 9 9 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.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 18000 -8068.744 -8068.744 -8134.7203 -8134.7203 255.3351 255.3351 24097.044 24097.044 1832.5091 1832.5091 19000 -8066.557 -8066.557 -8132.9283 -8132.9283 256.86374 256.86374 24105.733 24105.733 1699.7362 1699.7362 Loop time of 145.824 on 1 procs for 1000 steps with 2000 atoms Performance: 0.592 ns/day, 40.507 hours/ns, 6.858 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.31 | 145.31 | 145.31 | 0.0 | 99.65 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.095646 | 0.095646 | 0.095646 | 0.0 | 0.07 Output | 0.00020401 | 0.00020401 | 0.00020401 | 0.0 | 0.00 Modify | 0.35676 | 0.35676 | 0.35676 | 0.0 | 0.24 Other | | 0.05814 | | | 0.04 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 209588 ave 209588 max 209588 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 209588 Ave neighs/atom = 104.794 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.256738061567, Press = 1.16498623903938 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 = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 9 9 9 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.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 19000 -8066.557 -8066.557 -8132.9283 -8132.9283 256.86374 256.86374 24105.733 24105.733 1699.7362 1699.7362 20000 -8070.95 -8070.95 -8136.2247 -8136.2247 252.62019 252.62019 24102.749 24102.749 1038.1733 1038.1733 Loop time of 161.756 on 1 procs for 1000 steps with 2000 atoms Performance: 0.534 ns/day, 44.932 hours/ns, 6.182 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.18 | 161.18 | 161.18 | 0.0 | 99.64 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.10437 | 0.10437 | 0.10437 | 0.0 | 0.06 Output | 0.00019795 | 0.00019795 | 0.00019795 | 0.0 | 0.00 Modify | 0.40821 | 0.40821 | 0.40821 | 0.0 | 0.25 Other | | 0.0641 | | | 0.04 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 209426 ave 209426 max 209426 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 209426 Ave neighs/atom = 104.713 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.152296025879, Press = 2.61781141754734 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 = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 9 9 9 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.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 20000 -8070.95 -8070.95 -8136.2247 -8136.2247 252.62019 252.62019 24102.749 24102.749 1038.1733 1038.1733 21000 -8067.9203 -8067.9203 -8135.5671 -8135.5671 261.80006 261.80006 24114.053 24114.053 406.18015 406.18015 Loop time of 173.037 on 1 procs for 1000 steps with 2000 atoms Performance: 0.499 ns/day, 48.066 hours/ns, 5.779 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 | 172.42 | 172.42 | 172.42 | 0.0 | 99.64 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.11013 | 0.11013 | 0.11013 | 0.0 | 0.06 Output | 0.00016984 | 0.00016984 | 0.00016984 | 0.0 | 0.00 Modify | 0.44305 | 0.44305 | 0.44305 | 0.0 | 0.26 Other | | 0.06765 | | | 0.04 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 209586 ave 209586 max 209586 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 209586 Ave neighs/atom = 104.793 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.083578385764, Press = 2.74413628929256 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 = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 9 9 9 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.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 21000 -8067.9203 -8067.9203 -8135.5671 -8135.5671 261.80006 261.80006 24114.053 24114.053 406.18015 406.18015 22000 -8066.9428 -8066.9428 -8133.0873 -8133.0873 255.98605 255.98605 24119.421 24119.421 438.03044 438.03044 Loop time of 168.582 on 1 procs for 1000 steps with 2000 atoms Performance: 0.513 ns/day, 46.828 hours/ns, 5.932 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 | 167.98 | 167.98 | 167.98 | 0.0 | 99.64 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.10859 | 0.10859 | 0.10859 | 0.0 | 0.06 Output | 0.0002514 | 0.0002514 | 0.0002514 | 0.0 | 0.00 Modify | 0.42958 | 0.42958 | 0.42958 | 0.0 | 0.25 Other | | 0.06598 | | | 0.04 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 209282 ave 209282 max 209282 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 209282 Ave neighs/atom = 104.641 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.031142031947, Press = 0.464264742694407 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 = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 9 9 9 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.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 22000 -8066.9428 -8066.9428 -8133.0873 -8133.0873 255.98605 255.98605 24119.421 24119.421 438.03044 438.03044 23000 -8068.5611 -8068.5611 -8133.0677 -8133.0677 249.64729 249.64729 24121.267 24121.267 336.56417 336.56417 Loop time of 163.9 on 1 procs for 1000 steps with 2000 atoms Performance: 0.527 ns/day, 45.528 hours/ns, 6.101 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 | 163.31 | 163.31 | 163.31 | 0.0 | 99.64 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.1063 | 0.1063 | 0.1063 | 0.0 | 0.06 Output | 0.00020836 | 0.00020836 | 0.00020836 | 0.0 | 0.00 Modify | 0.4198 | 0.4198 | 0.4198 | 0.0 | 0.26 Other | | 0.06456 | | | 0.04 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 209220 ave 209220 max 209220 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 209220 Ave neighs/atom = 104.61 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.041861738228, Press = -0.378403139021122 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 = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 9 9 9 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.134 | 4.134 | 4.134 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 23000 -8068.5611 -8068.5611 -8133.0677 -8133.0677 249.64729 249.64729 24121.267 24121.267 336.56417 336.56417 24000 -8065.7577 -8065.7577 -8132.2111 -8132.2111 257.18132 257.18132 24122.2 24122.2 510.17006 510.17006 Loop time of 162.834 on 1 procs for 1000 steps with 2000 atoms Performance: 0.531 ns/day, 45.232 hours/ns, 6.141 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.25 | 162.25 | 162.25 | 0.0 | 99.64 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.10458 | 0.10458 | 0.10458 | 0.0 | 0.06 Output | 0.00020708 | 0.00020708 | 0.00020708 | 0.0 | 0.00 Modify | 0.41538 | 0.41538 | 0.41538 | 0.0 | 0.26 Other | | 0.06436 | | | 0.04 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 208980 ave 208980 max 208980 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 208980 Ave neighs/atom = 104.49 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 24125.0368946093 print "LAMMPS calculation completed" LAMMPS calculation completed quit 0