# 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 4.044650822877884*${_u_distance} variable latticeconst_converted equal 4.044650822877884*1 lattice fcc ${latticeconst_converted} lattice fcc 4.04465082287788 Lattice spacing in x,y,z = 4.0446508 4.0446508 4.0446508 region simbox block 0 10 0 10 0 10 units lattice create_box 1 simbox Created orthogonal box = (0.0000000 0.0000000 0.0000000) to (40.446508 40.446508 40.446508) 1 by 1 by 1 MPI processor grid create_atoms 1 box Created 4000 atoms using lattice units in orthogonal box = (0.0000000 0.0000000 0.0000000) to (40.446508 40.446508 40.446508) create_atoms CPU = 0.003 seconds variable mass_converted equal 26.981538*${_u_mass} variable mass_converted equal 26.981538*1 kim_interactions Al WARNING: 'kim_' has been renamed to 'kim '. Please update your input. kim interactions Al #=== BEGIN kim interactions ================================== pair_style kim MEAM_LAMMPS_KimKimLee_2009_AlMg__MO_058537087384_001 pair_coeff * * Al #=== END kim interactions ==================================== mass 1 ${mass_converted} mass 1 26.981538 # initial volume variable v equal vol # assign formula variable V0 equal ${v} # evaluate initial value variable V0 equal 66167.2528701086 variable V0_metal equal ${V0}/(${_u_distance}*${_u_distance}*${_u_distance}) variable V0_metal equal 66167.2528701086/(${_u_distance}*${_u_distance}*${_u_distance}) variable V0_metal equal 66167.2528701086/(1*${_u_distance}*${_u_distance}) variable V0_metal equal 66167.2528701086/(1*1*${_u_distance}) variable V0_metal equal 66167.2528701086/(1*1*1) variable V0_metal_times1000 equal ${V0_metal}*1000 variable V0_metal_times1000 equal 66167.2528701086*1000 print "Initial system volume: ${V0_metal} Angstroms^3" Initial system volume: 66167.2528701086 Angstroms^3 # set the time step to 0.001 picoseconds variable timestep_converted equal 0.001*${_u_time} variable timestep_converted equal 0.001*1 timestep ${timestep_converted} timestep 0.001 variable temp_converted equal 273.15*${_u_temperature} variable temp_converted equal 273.15*1 variable Tdamp_converted equal 0.1*${_u_time} variable Tdamp_converted equal 0.1*1 variable press_converted equal 0.0*${_u_pressure} variable press_converted equal 0.0*1 variable Pdamp_converted equal 1*${_u_time} variable Pdamp_converted equal 1*1 # create initial velocities consistent with the chosen temperature velocity all create ${temp_converted} 17 mom yes rot yes velocity all create 273.15 17 mom yes rot yes # set NPT ensemble for all atoms fix ensemble all npt temp ${temp_converted} ${temp_converted} ${Tdamp_converted} iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 273.15 ${temp_converted} ${Tdamp_converted} iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 273.15 273.15 ${Tdamp_converted} iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 273.15 273.15 0.1 iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 273.15 273.15 0.1 iso 0 ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 273.15 273.15 0.1 iso 0 0 ${Pdamp_converted} fix ensemble all npt temp 273.15 273.15 0.1 iso 0 0 1 # compute the time averages of pressure, temperature, and volume, respectively # ignore the first 5000 timesteps variable etotal_metal equal etotal/${_u_energy} variable etotal_metal equal etotal/1 variable pe_metal equal pe/${_u_energy} variable pe_metal equal pe/1 variable T_metal equal temp/${_u_temperature} variable T_metal equal temp/1 variable V_metal equal vol/(${_u_distance}*${_u_distance}*${_u_distance}) variable V_metal equal vol/(1*${_u_distance}*${_u_distance}) variable V_metal equal vol/(1*1*${_u_distance}) variable V_metal equal vol/(1*1*1) variable P_metal equal press/${_u_pressure} variable P_metal equal press/1 fix avgmyTemp all ave/time 5 20 100 v_T_metal ave running start 5000 fix avgmyPress all ave/time 5 20 100 v_P_metal ave running start 5000 fix avgmyVol all ave/time 5 20 100 v_V_metal ave running start 5000 # extract fix quantities into variables so they can be used in if-else logic later. variable T equal f_avgmyTemp variable P equal f_avgmyPress variable V equal f_avgmyVol # set error bounds for temperature and pressure in original metal units (K and bar) variable T_low equal "273.15 - 0.2" variable T_up equal "273.15 + 0.2" variable P_low equal "0.0 - 0.2" variable P_up equal "0.0 + 0.2" # print to logfile every 1000 timesteps thermo_style custom step etotal v_etotal_metal pe v_pe_metal temp v_T_metal vol v_V_metal press v_P_metal thermo 1000 # Run a simulation for at most 2000*1000 timesteps. At each 1000th time step, check # whether the temperature and pressure have converged. If yes, break. label top variable a loop 2000 run 1000 CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE Your simulation uses code contributions which should be cited: - OpenKIM: https://doi.org/10.1007/s11837-011-0102-6 @Article{tadmor:elliott:2011, author = {E. B. Tadmor and R. S. Elliott and J. P. Sethna and R. E. Miller and C. A. Becker}, title = {The potential of atomistic simulations and the {K}nowledgebase of {I}nteratomic {M}odels}, journal = {{JOM}}, year = 2011, volume = 63, number = 17, pages = {17}, doi = {10.1007/s11837-011-0102-6} } - OpenKIM potential: https://openkim.org/cite/MO_058537087384_001#item-citation CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 12 12 12 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) = 4.610 | 4.610 | 4.610 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 0 -13298.806 -13298.806 -13440 -13440 273.15 273.15 66167.253 66167.253 2279.2357 2279.2357 1000 -13142.643 -13142.643 -13293.567 -13293.567 291.97237 291.97237 67223.723 67223.723 796.51666 796.51666 Loop time of 408.197 on 1 procs for 1000 steps with 4000 atoms Performance: 0.212 ns/day, 113.388 hours/ns, 2.450 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 | 407.28 | 407.28 | 407.28 | 0.0 | 99.77 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.14905 | 0.14905 | 0.14905 | 0.0 | 0.04 Output | 0.00027159 | 0.00027159 | 0.00027159 | 0.0 | 0.00 Modify | 0.67581 | 0.67581 | 0.67581 | 0.0 | 0.17 Other | | 0.09638 | | | 0.02 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 312000.0 ave 312000 max 312000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 312000 Ave neighs/atom = 78.000000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 0, Press = 0 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 12 12 12 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) = 4.610 | 4.610 | 4.610 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 1000 -13142.643 -13142.643 -13293.567 -13293.567 291.97237 291.97237 67223.723 67223.723 796.51666 796.51666 2000 -13158.024 -13158.024 -13294.654 -13294.654 264.3213 264.3213 67217.302 67217.302 548.5469 548.5469 Loop time of 428.366 on 1 procs for 1000 steps with 4000 atoms Performance: 0.202 ns/day, 118.990 hours/ns, 2.334 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 | 427.43 | 427.43 | 427.43 | 0.0 | 99.78 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.15416 | 0.15416 | 0.15416 | 0.0 | 0.04 Output | 0.00021685 | 0.00021685 | 0.00021685 | 0.0 | 0.00 Modify | 0.68827 | 0.68827 | 0.68827 | 0.0 | 0.16 Other | | 0.09555 | | | 0.02 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 323510.0 ave 323510 max 323510 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 323510 Ave neighs/atom = 80.877500 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 0, Press = 0 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 12 12 12 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) = 4.610 | 4.610 | 4.610 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 2000 -13158.024 -13158.024 -13294.654 -13294.654 264.3213 264.3213 67217.302 67217.302 548.5469 548.5469 3000 -13148.569 -13148.569 -13286.727 -13286.727 267.27541 267.27541 67264.825 67264.825 485.79871 485.79871 Loop time of 455.972 on 1 procs for 1000 steps with 4000 atoms Performance: 0.189 ns/day, 126.659 hours/ns, 2.193 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 | 454.95 | 454.95 | 454.95 | 0.0 | 99.77 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16465 | 0.16465 | 0.16465 | 0.0 | 0.04 Output | 0.00027534 | 0.00027534 | 0.00027534 | 0.0 | 0.00 Modify | 0.76084 | 0.76084 | 0.76084 | 0.0 | 0.17 Other | | 0.1002 | | | 0.02 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 323620.0 ave 323620 max 323620 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 323620 Ave neighs/atom = 80.905000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 0, Press = 0 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 12 12 12 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) = 4.610 | 4.610 | 4.610 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 3000 -13148.569 -13148.569 -13286.727 -13286.727 267.27541 267.27541 67264.825 67264.825 485.79871 485.79871 4000 -13155.15 -13155.15 -13294.921 -13294.921 270.39709 270.39709 67251.479 67251.479 124.87675 124.87675 Loop time of 469.609 on 1 procs for 1000 steps with 4000 atoms Performance: 0.184 ns/day, 130.447 hours/ns, 2.129 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 | 468.52 | 468.52 | 468.52 | 0.0 | 99.77 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16645 | 0.16645 | 0.16645 | 0.0 | 0.04 Output | 0.00028204 | 0.00028204 | 0.00028204 | 0.0 | 0.00 Modify | 0.81196 | 0.81196 | 0.81196 | 0.0 | 0.17 Other | | 0.1061 | | | 0.02 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 323480.0 ave 323480 max 323480 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 323480 Ave neighs/atom = 80.870000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 0, Press = 0 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 12 12 12 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) = 4.610 | 4.610 | 4.610 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 4000 -13155.15 -13155.15 -13294.921 -13294.921 270.39709 270.39709 67251.479 67251.479 124.87675 124.87675 5000 -13152.5 -13152.5 -13293.216 -13293.216 272.22618 272.22618 67249.256 67249.256 307.60016 307.60016 Loop time of 459.539 on 1 procs for 1000 steps with 4000 atoms Performance: 0.188 ns/day, 127.650 hours/ns, 2.176 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 | 458.5 | 458.5 | 458.5 | 0.0 | 99.77 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16117 | 0.16117 | 0.16117 | 0.0 | 0.04 Output | 0.00022978 | 0.00022978 | 0.00022978 | 0.0 | 0.00 Modify | 0.77331 | 0.77331 | 0.77331 | 0.0 | 0.17 Other | | 0.1031 | | | 0.02 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 323308.0 ave 323308 max 323308 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 323308 Ave neighs/atom = 80.827000 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 = 274.980909931612, Press = 252.621083305634 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 12 12 12 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) = 4.610 | 4.610 | 4.610 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 5000 -13152.5 -13152.5 -13293.216 -13293.216 272.22618 272.22618 67249.256 67249.256 307.60016 307.60016 6000 -13150.846 -13150.846 -13290.927 -13290.927 270.99449 270.99449 67245.185 67245.185 512.89628 512.89628 Loop time of 473.272 on 1 procs for 1000 steps with 4000 atoms Performance: 0.183 ns/day, 131.464 hours/ns, 2.113 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 | 472.15 | 472.15 | 472.15 | 0.0 | 99.76 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16577 | 0.16577 | 0.16577 | 0.0 | 0.04 Output | 0.00017861 | 0.00017861 | 0.00017861 | 0.0 | 0.00 Modify | 0.84557 | 0.84557 | 0.84557 | 0.0 | 0.18 Other | | 0.1087 | | | 0.02 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 323684.0 ave 323684 max 323684 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 323684 Ave neighs/atom = 80.921000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 272.680665166355, Press = 41.1093622121164 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 12 12 12 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) = 4.610 | 4.610 | 4.610 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 6000 -13150.846 -13150.846 -13290.927 -13290.927 270.99449 270.99449 67245.185 67245.185 512.89628 512.89628 7000 -13154.276 -13154.276 -13293.203 -13293.203 268.76433 268.76433 67254.082 67254.082 153.95151 153.95151 Loop time of 473.546 on 1 procs for 1000 steps with 4000 atoms Performance: 0.182 ns/day, 131.541 hours/ns, 2.112 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 | 472.42 | 472.42 | 472.42 | 0.0 | 99.76 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.1663 | 0.1663 | 0.1663 | 0.0 | 0.04 Output | 0.00017846 | 0.00017846 | 0.00017846 | 0.0 | 0.00 Modify | 0.85105 | 0.85105 | 0.85105 | 0.0 | 0.18 Other | | 0.1077 | | | 0.02 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 323410.0 ave 323410 max 323410 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 323410 Ave neighs/atom = 80.852500 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 273.170652872484, Press = 28.6957456771841 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 12 12 12 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) = 4.610 | 4.610 | 4.610 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 7000 -13154.276 -13154.276 -13293.203 -13293.203 268.76433 268.76433 67254.082 67254.082 153.95151 153.95151 8000 -13151.428 -13151.428 -13294.172 -13294.172 276.14942 276.14942 67298.294 67298.294 -317.79722 -317.79722 Loop time of 465.039 on 1 procs for 1000 steps with 4000 atoms Performance: 0.186 ns/day, 129.177 hours/ns, 2.150 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 | 463.94 | 463.94 | 463.94 | 0.0 | 99.76 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16202 | 0.16202 | 0.16202 | 0.0 | 0.03 Output | 0.00027413 | 0.00027413 | 0.00027413 | 0.0 | 0.00 Modify | 0.8261 | 0.8261 | 0.8261 | 0.0 | 0.18 Other | | 0.1058 | | | 0.02 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 323348.0 ave 323348 max 323348 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 323348 Ave neighs/atom = 80.837000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 272.958619260602, Press = 13.7004190636696 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 12 12 12 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) = 4.610 | 4.610 | 4.610 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 8000 -13151.428 -13151.428 -13294.172 -13294.172 276.14942 276.14942 67298.294 67298.294 -317.79722 -317.79722 9000 -13152.828 -13152.828 -13295.295 -13295.295 275.61251 275.61251 67288.272 67288.272 -256.96951 -256.96951 Loop time of 471.898 on 1 procs for 1000 steps with 4000 atoms Performance: 0.183 ns/day, 131.083 hours/ns, 2.119 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 | 470.81 | 470.81 | 470.81 | 0.0 | 99.77 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16307 | 0.16307 | 0.16307 | 0.0 | 0.03 Output | 0.00023322 | 0.00023322 | 0.00023322 | 0.0 | 0.00 Modify | 0.81943 | 0.81943 | 0.81943 | 0.0 | 0.17 Other | | 0.1056 | | | 0.02 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 323256.0 ave 323256 max 323256 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 323256 Ave neighs/atom = 80.814000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 272.779282379509, Press = 7.77809755483617 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 12 12 12 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) = 4.610 | 4.610 | 4.610 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 9000 -13152.828 -13152.828 -13295.295 -13295.295 275.61251 275.61251 67288.272 67288.272 -256.96951 -256.96951 10000 -13149.107 -13149.107 -13290.709 -13290.709 273.93852 273.93852 67298.696 67298.696 -66.669171 -66.669171 Loop time of 473.892 on 1 procs for 1000 steps with 4000 atoms Performance: 0.182 ns/day, 131.637 hours/ns, 2.110 timesteps/s 99.9% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 472.8 | 472.8 | 472.8 | 0.0 | 99.77 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.1655 | 0.1655 | 0.1655 | 0.0 | 0.03 Output | 0.00027541 | 0.00027541 | 0.00027541 | 0.0 | 0.00 Modify | 0.82377 | 0.82377 | 0.82377 | 0.0 | 0.17 Other | | 0.1056 | | | 0.02 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 323260.0 ave 323260 max 323260 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 323260 Ave neighs/atom = 80.815000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 272.629787332656, Press = 3.62836503073575 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 12 12 12 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) = 4.610 | 4.610 | 4.610 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 10000 -13149.107 -13149.107 -13290.709 -13290.709 273.93852 273.93852 67298.696 67298.696 -66.669171 -66.669171 11000 -13154.954 -13154.954 -13295.291 -13295.291 271.49109 271.49109 67277.595 67277.595 -236.46146 -236.46146 Loop time of 455.238 on 1 procs for 1000 steps with 4000 atoms Performance: 0.190 ns/day, 126.455 hours/ns, 2.197 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 | 454.19 | 454.19 | 454.19 | 0.0 | 99.77 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.15863 | 0.15863 | 0.15863 | 0.0 | 0.03 Output | 0.00018009 | 0.00018009 | 0.00018009 | 0.0 | 0.00 Modify | 0.78212 | 0.78212 | 0.78212 | 0.0 | 0.17 Other | | 0.1036 | | | 0.02 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 323178.0 ave 323178 max 323178 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 323178 Ave neighs/atom = 80.794500 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 272.737666727895, Press = 1.99806537507249 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 12 12 12 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) = 4.610 | 4.610 | 4.610 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 11000 -13154.954 -13154.954 -13295.291 -13295.291 271.49109 271.49109 67277.595 67277.595 -236.46146 -236.46146 12000 -13150.595 -13150.595 -13292.457 -13292.457 274.44161 274.44161 67300.348 67300.348 -184.13692 -184.13692 Loop time of 465.749 on 1 procs for 1000 steps with 4000 atoms Performance: 0.186 ns/day, 129.375 hours/ns, 2.147 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 | 464.68 | 464.68 | 464.68 | 0.0 | 99.77 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.162 | 0.162 | 0.162 | 0.0 | 0.03 Output | 0.00018056 | 0.00018056 | 0.00018056 | 0.0 | 0.00 Modify | 0.80285 | 0.80285 | 0.80285 | 0.0 | 0.17 Other | | 0.1036 | | | 0.02 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 323306.0 ave 323306 max 323306 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 323306 Ave neighs/atom = 80.826500 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 272.723426284027, Press = 2.48356255542784 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 12 12 12 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) = 4.610 | 4.610 | 4.610 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 12000 -13150.595 -13150.595 -13292.457 -13292.457 274.44161 274.44161 67300.348 67300.348 -184.13692 -184.13692 13000 -13158.918 -13158.918 -13297.587 -13297.587 268.26483 268.26483 67278.01 67278.01 -387.04183 -387.04183 Loop time of 442.039 on 1 procs for 1000 steps with 4000 atoms Performance: 0.195 ns/day, 122.789 hours/ns, 2.262 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 | 441.04 | 441.04 | 441.04 | 0.0 | 99.77 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.15381 | 0.15381 | 0.15381 | 0.0 | 0.03 Output | 0.00017615 | 0.00017615 | 0.00017615 | 0.0 | 0.00 Modify | 0.74402 | 0.74402 | 0.74402 | 0.0 | 0.17 Other | | 0.0997 | | | 0.02 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 323170.0 ave 323170 max 323170 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 323170 Ave neighs/atom = 80.792500 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 272.609831535214, Press = 2.87494562521845 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 12 12 12 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) = 4.610 | 4.610 | 4.610 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 13000 -13158.918 -13158.918 -13297.587 -13297.587 268.26483 268.26483 67278.01 67278.01 -387.04183 -387.04183 14000 -13150.105 -13150.105 -13291.811 -13291.811 274.14009 274.14009 67341.138 67341.138 -671.65694 -671.65694 Loop time of 472.524 on 1 procs for 1000 steps with 4000 atoms Performance: 0.183 ns/day, 131.257 hours/ns, 2.116 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 | 471.42 | 471.42 | 471.42 | 0.0 | 99.77 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16399 | 0.16399 | 0.16399 | 0.0 | 0.03 Output | 0.00023542 | 0.00023542 | 0.00023542 | 0.0 | 0.00 Modify | 0.83267 | 0.83267 | 0.83267 | 0.0 | 0.18 Other | | 0.1063 | | | 0.02 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 323214.0 ave 323214 max 323214 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 323214 Ave neighs/atom = 80.803500 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 272.603101235258, Press = 1.18673635041062 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 12 12 12 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) = 4.610 | 4.610 | 4.610 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 14000 -13150.105 -13150.105 -13291.811 -13291.811 274.14009 274.14009 67341.138 67341.138 -671.65694 -671.65694 15000 -13156.683 -13156.683 -13294.549 -13294.549 266.71095 266.71095 67347.828 67347.828 -1012.3921 -1012.3921 Loop time of 468.337 on 1 procs for 1000 steps with 4000 atoms Performance: 0.184 ns/day, 130.094 hours/ns, 2.135 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 | 467.24 | 467.24 | 467.24 | 0.0 | 99.77 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16407 | 0.16407 | 0.16407 | 0.0 | 0.04 Output | 0.00017872 | 0.00017872 | 0.00017872 | 0.0 | 0.00 Modify | 0.82431 | 0.82431 | 0.82431 | 0.0 | 0.18 Other | | 0.1055 | | | 0.02 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 323074.0 ave 323074 max 323074 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 323074 Ave neighs/atom = 80.768500 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 272.519500258284, Press = -0.779162999765944 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 12 12 12 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) = 4.610 | 4.610 | 4.610 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 15000 -13156.683 -13156.683 -13294.549 -13294.549 266.71095 266.71095 67347.828 67347.828 -1012.3921 -1012.3921 16000 -13150.456 -13150.456 -13291.126 -13291.126 272.13458 272.13458 67363.825 67363.825 -886.05444 -886.05444 Loop time of 465.278 on 1 procs for 1000 steps with 4000 atoms Performance: 0.186 ns/day, 129.244 hours/ns, 2.149 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 | 464.18 | 464.18 | 464.18 | 0.0 | 99.76 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16533 | 0.16533 | 0.16533 | 0.0 | 0.04 Output | 0.00022669 | 0.00022669 | 0.00022669 | 0.0 | 0.00 Modify | 0.82417 | 0.82417 | 0.82417 | 0.0 | 0.18 Other | | 0.1059 | | | 0.02 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 323000.0 ave 323000 max 323000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 323000 Ave neighs/atom = 80.750000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 272.5159412, Press = -1.50264951383958 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 12 12 12 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) = 4.610 | 4.610 | 4.610 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 16000 -13150.456 -13150.456 -13291.126 -13291.126 272.13458 272.13458 67363.825 67363.825 -886.05444 -886.05444 17000 -13153.263 -13153.263 -13293.88 -13293.88 272.03269 272.03269 67293.779 67293.779 -254.64669 -254.64669 Loop time of 468.306 on 1 procs for 1000 steps with 4000 atoms Performance: 0.184 ns/day, 130.085 hours/ns, 2.135 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 | 467.21 | 467.21 | 467.21 | 0.0 | 99.77 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16246 | 0.16246 | 0.16246 | 0.0 | 0.03 Output | 0.00018156 | 0.00018156 | 0.00018156 | 0.0 | 0.00 Modify | 0.82404 | 0.82404 | 0.82404 | 0.0 | 0.18 Other | | 0.1056 | | | 0.02 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 323098.0 ave 323098 max 323098 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 323098 Ave neighs/atom = 80.774500 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 272.447378583263, Press = -0.367809041782955 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 12 12 12 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) = 4.610 | 4.610 | 4.610 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 17000 -13153.263 -13153.263 -13293.88 -13293.88 272.03269 272.03269 67293.779 67293.779 -254.64669 -254.64669 18000 -13150.27 -13150.27 -13292.642 -13292.642 275.42707 275.42707 67266.749 67266.749 207.74618 207.74618 Loop time of 473.833 on 1 procs for 1000 steps with 4000 atoms Performance: 0.182 ns/day, 131.620 hours/ns, 2.110 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 | 472.72 | 472.72 | 472.72 | 0.0 | 99.77 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16613 | 0.16613 | 0.16613 | 0.0 | 0.04 Output | 0.00017945 | 0.00017945 | 0.00017945 | 0.0 | 0.00 Modify | 0.83847 | 0.83847 | 0.83847 | 0.0 | 0.18 Other | | 0.107 | | | 0.02 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 323320.0 ave 323320 max 323320 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 323320 Ave neighs/atom = 80.830000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 272.435480521222, Press = -0.158563158450825 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 12 12 12 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) = 4.610 | 4.610 | 4.610 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 18000 -13150.27 -13150.27 -13292.642 -13292.642 275.42707 275.42707 67266.749 67266.749 207.74618 207.74618 19000 -13148.043 -13148.043 -13289.747 -13289.747 274.13667 274.13667 67247.069 67247.069 651.51936 651.51936 Loop time of 468.502 on 1 procs for 1000 steps with 4000 atoms Performance: 0.184 ns/day, 130.140 hours/ns, 2.134 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 | 467.41 | 467.41 | 467.41 | 0.0 | 99.77 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16425 | 0.16425 | 0.16425 | 0.0 | 0.04 Output | 0.00022856 | 0.00022856 | 0.00022856 | 0.0 | 0.00 Modify | 0.82696 | 0.82696 | 0.82696 | 0.0 | 0.18 Other | | 0.1057 | | | 0.02 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 323516.0 ave 323516 max 323516 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 323516 Ave neighs/atom = 80.879000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 272.464336320854, Press = 0.526717675388053 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 12 12 12 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) = 4.610 | 4.610 | 4.610 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 19000 -13148.043 -13148.043 -13289.747 -13289.747 274.13667 274.13667 67247.069 67247.069 651.51936 651.51936 20000 -13152.703 -13152.703 -13290.511 -13290.511 266.59843 266.59843 67202.937 67202.937 982.33472 982.33472 Loop time of 451.442 on 1 procs for 1000 steps with 4000 atoms Performance: 0.191 ns/day, 125.401 hours/ns, 2.215 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 | 450.4 | 450.4 | 450.4 | 0.0 | 99.77 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.15912 | 0.15912 | 0.15912 | 0.0 | 0.04 Output | 0.00022898 | 0.00022898 | 0.00022898 | 0.0 | 0.00 Modify | 0.77841 | 0.77841 | 0.77841 | 0.0 | 0.17 Other | | 0.1032 | | | 0.02 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 323408.0 ave 323408 max 323408 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 323408 Ave neighs/atom = 80.852000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 272.625904589655, Press = 0.871054347782962 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 12 12 12 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) = 4.610 | 4.610 | 4.610 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 20000 -13152.703 -13152.703 -13290.511 -13290.511 266.59843 266.59843 67202.937 67202.937 982.33472 982.33472 21000 -13149.421 -13149.421 -13291.078 -13291.078 274.04539 274.04539 67172.928 67172.928 1319.9674 1319.9674 Loop time of 475.145 on 1 procs for 1000 steps with 4000 atoms Performance: 0.182 ns/day, 131.985 hours/ns, 2.105 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 | 474 | 474 | 474 | 0.0 | 99.76 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16645 | 0.16645 | 0.16645 | 0.0 | 0.04 Output | 0.00018124 | 0.00018124 | 0.00018124 | 0.0 | 0.00 Modify | 0.86643 | 0.86643 | 0.86643 | 0.0 | 0.18 Other | | 0.1096 | | | 0.02 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 323616.0 ave 323616 max 323616 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 323616 Ave neighs/atom = 80.904000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 272.667477295433, Press = 1.60798481317329 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 12 12 12 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) = 4.610 | 4.610 | 4.610 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 21000 -13149.421 -13149.421 -13291.078 -13291.078 274.04539 274.04539 67172.928 67172.928 1319.9674 1319.9674 22000 -13153.475 -13153.475 -13293.725 -13293.725 271.32218 271.32218 67153.86 67153.86 1316.0166 1316.0166 Loop time of 459.241 on 1 procs for 1000 steps with 4000 atoms Performance: 0.188 ns/day, 127.567 hours/ns, 2.178 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 | 458.16 | 458.16 | 458.16 | 0.0 | 99.76 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16101 | 0.16101 | 0.16101 | 0.0 | 0.04 Output | 0.00017686 | 0.00017686 | 0.00017686 | 0.0 | 0.00 Modify | 0.81491 | 0.81491 | 0.81491 | 0.0 | 0.18 Other | | 0.1056 | | | 0.02 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 323784.0 ave 323784 max 323784 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 323784 Ave neighs/atom = 80.946000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 272.714777835033, Press = 2.30359229480986 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 12 12 12 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) = 4.610 | 4.610 | 4.610 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 22000 -13153.475 -13153.475 -13293.725 -13293.725 271.32218 271.32218 67153.86 67153.86 1316.0166 1316.0166 23000 -13148.681 -13148.681 -13291.17 -13291.17 275.65483 275.65483 67256.54 67256.54 365.58497 365.58497 Loop time of 442.473 on 1 procs for 1000 steps with 4000 atoms Performance: 0.195 ns/day, 122.909 hours/ns, 2.260 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 | 441.47 | 441.47 | 441.47 | 0.0 | 99.77 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.15481 | 0.15481 | 0.15481 | 0.0 | 0.03 Output | 0.00024286 | 0.00024286 | 0.00024286 | 0.0 | 0.00 Modify | 0.75027 | 0.75027 | 0.75027 | 0.0 | 0.17 Other | | 0.1011 | | | 0.02 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 323598.0 ave 323598 max 323598 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 323598 Ave neighs/atom = 80.899500 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 272.801934646502, Press = 1.53383007976885 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 12 12 12 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) = 4.610 | 4.610 | 4.610 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 23000 -13148.681 -13148.681 -13291.17 -13291.17 275.65483 275.65483 67256.54 67256.54 365.58497 365.58497 24000 -13152.424 -13152.424 -13293.456 -13293.456 272.83679 272.83679 67279.312 67279.312 -53.882632 -53.882632 Loop time of 472.041 on 1 procs for 1000 steps with 4000 atoms Performance: 0.183 ns/day, 131.123 hours/ns, 2.118 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 | 470.93 | 470.93 | 470.93 | 0.0 | 99.76 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16389 | 0.16389 | 0.16389 | 0.0 | 0.03 Output | 0.00017624 | 0.00017624 | 0.00017624 | 0.0 | 0.00 Modify | 0.84313 | 0.84313 | 0.84313 | 0.0 | 0.18 Other | | 0.108 | | | 0.02 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 323556.0 ave 323556 max 323556 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 323556 Ave neighs/atom = 80.889000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 272.850766055833, Press = 1.15163381609551 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 12 12 12 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) = 4.610 | 4.610 | 4.610 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 24000 -13152.424 -13152.424 -13293.456 -13293.456 272.83679 272.83679 67279.312 67279.312 -53.882632 -53.882632 25000 -13148.086 -13148.086 -13290.725 -13290.725 275.94491 275.94491 67303.161 67303.161 -138.33147 -138.33147 Loop time of 470.343 on 1 procs for 1000 steps with 4000 atoms Performance: 0.184 ns/day, 130.651 hours/ns, 2.126 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 | 469.23 | 469.23 | 469.23 | 0.0 | 99.76 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.1645 | 0.1645 | 0.1645 | 0.0 | 0.03 Output | 0.00022534 | 0.00022534 | 0.00022534 | 0.0 | 0.00 Modify | 0.8406 | 0.8406 | 0.8406 | 0.0 | 0.18 Other | | 0.108 | | | 0.02 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 323306.0 ave 323306 max 323306 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 323306 Ave neighs/atom = 80.826500 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 272.957041431475, Press = 0.335861806566941 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 12 12 12 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) = 4.610 | 4.610 | 4.610 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 25000 -13148.086 -13148.086 -13290.725 -13290.725 275.94491 275.94491 67303.161 67303.161 -138.33147 -138.33147 26000 -13153.452 -13153.452 -13294.166 -13294.166 272.22047 272.22047 67258.466 67258.466 130.84384 130.84384 Loop time of 429.946 on 1 procs for 1000 steps with 4000 atoms Performance: 0.201 ns/day, 119.429 hours/ns, 2.326 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 | 428.98 | 428.98 | 428.98 | 0.0 | 99.77 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.14933 | 0.14933 | 0.14933 | 0.0 | 0.03 Output | 0.00017554 | 0.00017554 | 0.00017554 | 0.0 | 0.00 Modify | 0.71972 | 0.71972 | 0.71972 | 0.0 | 0.17 Other | | 0.1 | | | 0.02 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 323252.0 ave 323252 max 323252 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 323252 Ave neighs/atom = 80.813000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 272.967641206552, Press = 0.375162470801439 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 12 12 12 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) = 4.610 | 4.610 | 4.610 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 26000 -13153.452 -13153.452 -13294.166 -13294.166 272.22047 272.22047 67258.466 67258.466 130.84384 130.84384 27000 -13151.058 -13151.058 -13293.093 -13293.093 274.77653 274.77653 67247.35 67247.35 318.8669 318.8669 Loop time of 447.981 on 1 procs for 1000 steps with 4000 atoms Performance: 0.193 ns/day, 124.439 hours/ns, 2.232 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 | 446.95 | 446.95 | 446.95 | 0.0 | 99.77 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.15524 | 0.15524 | 0.15524 | 0.0 | 0.03 Output | 0.00024401 | 0.00024401 | 0.00024401 | 0.0 | 0.00 Modify | 0.76841 | 0.76841 | 0.76841 | 0.0 | 0.17 Other | | 0.1032 | | | 0.02 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 323346.0 ave 323346 max 323346 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 323346 Ave neighs/atom = 80.836500 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 273.055451739071, Press = -0.484370195839726 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 12 12 12 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) = 4.610 | 4.610 | 4.610 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 27000 -13151.058 -13151.058 -13293.093 -13293.093 274.77653 274.77653 67247.35 67247.35 318.8669 318.8669 28000 -13148.355 -13148.355 -13289.703 -13289.703 273.44719 273.44719 67235.613 67235.613 717.24627 717.24627 Loop time of 467.305 on 1 procs for 1000 steps with 4000 atoms Performance: 0.185 ns/day, 129.807 hours/ns, 2.140 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 | 466.2 | 466.2 | 466.2 | 0.0 | 99.76 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.1639 | 0.1639 | 0.1639 | 0.0 | 0.04 Output | 0.00022527 | 0.00022527 | 0.00022527 | 0.0 | 0.00 Modify | 0.83026 | 0.83026 | 0.83026 | 0.0 | 0.18 Other | | 0.1069 | | | 0.02 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 323530.0 ave 323530 max 323530 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 323530 Ave neighs/atom = 80.882500 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 273.083101409147, Press = -0.439531140112757 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 12 12 12 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) = 4.610 | 4.610 | 4.610 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 28000 -13148.355 -13148.355 -13289.703 -13289.703 273.44719 273.44719 67235.613 67235.613 717.24627 717.24627 29000 -13154.572 -13154.572 -13294.057 -13294.057 269.84354 269.84354 67193.069 67193.069 862.72691 862.72691 Loop time of 434.2 on 1 procs for 1000 steps with 4000 atoms Performance: 0.199 ns/day, 120.611 hours/ns, 2.303 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 | 433.2 | 433.2 | 433.2 | 0.0 | 99.77 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.15266 | 0.15266 | 0.15266 | 0.0 | 0.04 Output | 0.00017654 | 0.00017654 | 0.00017654 | 0.0 | 0.00 Modify | 0.74249 | 0.74249 | 0.74249 | 0.0 | 0.17 Other | | 0.1017 | | | 0.02 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 323656.0 ave 323656 max 323656 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 323656 Ave neighs/atom = 80.914000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 273.111574533939, Press = -0.405170579043809 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7 ghost atom cutoff = 7 binsize = 3.5, bins = 12 12 12 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) = 4.610 | 4.610 | 4.610 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 29000 -13154.572 -13154.572 -13294.057 -13294.057 269.84354 269.84354 67193.069 67193.069 862.72691 862.72691 30000 -13150.076 -13150.076 -13292.139 -13292.139 274.83083 274.83083 67191.817 67191.817 1044.0559 1044.0559 Loop time of 473.93 on 1 procs for 1000 steps with 4000 atoms Performance: 0.182 ns/day, 131.647 hours/ns, 2.110 timesteps/s 99.9% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 472.81 | 472.81 | 472.81 | 0.0 | 99.76 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16468 | 0.16468 | 0.16468 | 0.0 | 0.03 Output | 0.00017851 | 0.00017851 | 0.00017851 | 0.0 | 0.00 Modify | 0.84453 | 0.84453 | 0.84453 | 0.0 | 0.18 Other | | 0.1076 | | | 0.02 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 323628.0 ave 323628 max 323628 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 323628 Ave neighs/atom = 80.907000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" jump SELF break # Write final averaged volume to file if temperature and volume have converged; otherwise wirte a # flag to indicate non-convergence. variable myStep equal step if "${myStep} < 2000000" then "print '${V}' file output/vol_T273.15.out" else "print 'not_converged' file output/vol_T273.15.out" print '${V}' file output/vol_T273.15.out 67274.6867865151 print "LAMMPS calculation completed" LAMMPS calculation completed quit 0