# 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 3.9173715934157367*${_u_distance} variable latticeconst_converted equal 3.9173715934157367*1 lattice fcc ${latticeconst_converted} lattice fcc 3.91737159341574 Lattice spacing in x,y,z = 3.9173716 3.9173716 3.9173716 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 (39.173716 39.173716 39.173716) 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 (39.173716 39.173716 39.173716) create_atoms CPU = 0.004 seconds variable mass_converted equal 195.078*${_u_mass} variable mass_converted equal 195.078*1 kim_interactions Pt WARNING: 'kim_' has been renamed to 'kim '. Please update your input. kim interactions Pt #=== BEGIN kim interactions ================================== pair_style kim MEAM_LAMMPS_KimLee_2006_PtFe__MO_343168101490_001 pair_coeff * * Pt #=== END kim interactions ==================================== mass 1 ${mass_converted} mass 1 195.078 # initial volume variable v equal vol # assign formula variable V0 equal ${v} # evaluate initial value variable V0 equal 60115.2017852415 variable V0_metal equal ${V0}/(${_u_distance}*${_u_distance}*${_u_distance}) variable V0_metal equal 60115.2017852415/(${_u_distance}*${_u_distance}*${_u_distance}) variable V0_metal equal 60115.2017852415/(1*${_u_distance}*${_u_distance}) variable V0_metal equal 60115.2017852415/(1*1*${_u_distance}) variable V0_metal equal 60115.2017852415/(1*1*1) variable V0_metal_times1000 equal ${V0_metal}*1000 variable V0_metal_times1000 equal 60115.2017852415*1000 print "Initial system volume: ${V0_metal} Angstroms^3" Initial system volume: 60115.2017852415 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 293.15*${_u_temperature} variable temp_converted equal 293.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 293.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 293.15 ${temp_converted} ${Tdamp_converted} iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 293.15 293.15 ${Tdamp_converted} iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 293.15 293.15 0.1 iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 293.15 293.15 0.1 iso 0 ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 293.15 293.15 0.1 iso 0 0 ${Pdamp_converted} fix ensemble all npt temp 293.15 293.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 "293.15 - 0.2" variable T_up equal "293.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_343168101490_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 = 6.4 ghost atom cutoff = 6.4 binsize = 3.2, bins = 13 13 13 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.4 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 0 -22928.467 -22928.467 -23080 -23080 293.15 293.15 60115.202 60115.202 2692.3908 2692.3908 1000 -22768.752 -22768.752 -22929.582 -22929.582 311.13782 311.13782 60647.371 60647.371 -89.208743 -89.208743 Loop time of 188.991 on 1 procs for 1000 steps with 4000 atoms Performance: 0.457 ns/day, 52.497 hours/ns, 5.291 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 | 188.02 | 188.02 | 188.02 | 0.0 | 99.49 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.15654 | 0.15654 | 0.15654 | 0.0 | 0.08 Output | 0.00022407 | 0.00022407 | 0.00022407 | 0.0 | 0.00 Modify | 0.71782 | 0.71782 | 0.71782 | 0.0 | 0.38 Other | | 0.09744 | | | 0.05 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 = 6.4 ghost atom cutoff = 6.4 binsize = 3.2, bins = 13 13 13 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.4 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 1000 -22768.752 -22768.752 -22929.582 -22929.582 311.13782 311.13782 60647.371 60647.371 -89.208743 -89.208743 2000 -22777.151 -22777.151 -22933.308 -22933.308 302.09624 302.09624 60659.733 60659.733 -1364.6647 -1364.6647 Loop time of 186.231 on 1 procs for 1000 steps with 4000 atoms Performance: 0.464 ns/day, 51.731 hours/ns, 5.370 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 | 185.35 | 185.35 | 185.35 | 0.0 | 99.53 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.14573 | 0.14573 | 0.14573 | 0.0 | 0.08 Output | 0.00022926 | 0.00022926 | 0.00022926 | 0.0 | 0.00 Modify | 0.6405 | 0.6405 | 0.6405 | 0.0 | 0.34 Other | | 0.09226 | | | 0.05 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: 310748.0 ave 310748 max 310748 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 310748 Ave neighs/atom = 77.687000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 0, Press = 0 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.4 ghost atom cutoff = 6.4 binsize = 3.2, bins = 13 13 13 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.4 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 2000 -22777.151 -22777.151 -22933.308 -22933.308 302.09624 302.09624 60659.733 60659.733 -1364.6647 -1364.6647 3000 -22778.632 -22778.632 -22933.568 -22933.568 299.73463 299.73463 60671.347 60671.347 -2001.8154 -2001.8154 Loop time of 187.024 on 1 procs for 1000 steps with 4000 atoms Performance: 0.462 ns/day, 51.951 hours/ns, 5.347 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 | 186.14 | 186.14 | 186.14 | 0.0 | 99.53 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.14836 | 0.14836 | 0.14836 | 0.0 | 0.08 Output | 0.0002138 | 0.0002138 | 0.0002138 | 0.0 | 0.00 Modify | 0.64461 | 0.64461 | 0.64461 | 0.0 | 0.34 Other | | 0.09342 | | | 0.05 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: 310948.0 ave 310948 max 310948 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 310948 Ave neighs/atom = 77.737000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 0, Press = 0 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.4 ghost atom cutoff = 6.4 binsize = 3.2, bins = 13 13 13 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.4 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 3000 -22778.632 -22778.632 -22933.568 -22933.568 299.73463 299.73463 60671.347 60671.347 -2001.8154 -2001.8154 4000 -22774.074 -22774.074 -22931.335 -22931.335 304.23199 304.23199 60585.726 60585.726 2642.9195 2642.9195 Loop time of 185.543 on 1 procs for 1000 steps with 4000 atoms Performance: 0.466 ns/day, 51.540 hours/ns, 5.390 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 | 184.67 | 184.67 | 184.67 | 0.0 | 99.53 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.14645 | 0.14645 | 0.14645 | 0.0 | 0.08 Output | 0.0002174 | 0.0002174 | 0.0002174 | 0.0 | 0.00 Modify | 0.63803 | 0.63803 | 0.63803 | 0.0 | 0.34 Other | | 0.09337 | | | 0.05 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: 310662.0 ave 310662 max 310662 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 310662 Ave neighs/atom = 77.665500 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 0, Press = 0 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.4 ghost atom cutoff = 6.4 binsize = 3.2, bins = 13 13 13 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.4 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 4000 -22774.074 -22774.074 -22931.335 -22931.335 304.23199 304.23199 60585.726 60585.726 2642.9195 2642.9195 5000 -22778.839 -22778.839 -22927.907 -22927.907 288.38119 288.38119 60654.154 60654.154 -502.07711 -502.07711 Loop time of 188.007 on 1 procs for 1000 steps with 4000 atoms Performance: 0.460 ns/day, 52.224 hours/ns, 5.319 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 | 187.11 | 187.11 | 187.11 | 0.0 | 99.52 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.14619 | 0.14619 | 0.14619 | 0.0 | 0.08 Output | 0.00022738 | 0.00022738 | 0.00022738 | 0.0 | 0.00 Modify | 0.65265 | 0.65265 | 0.65265 | 0.0 | 0.35 Other | | 0.09438 | | | 0.05 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: 311044.0 ave 311044 max 311044 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311044 Ave neighs/atom = 77.761000 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 = 296.183634582989, Press = 188.057458933668 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.4 ghost atom cutoff = 6.4 binsize = 3.2, bins = 13 13 13 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.4 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 5000 -22778.839 -22778.839 -22927.907 -22927.907 288.38119 288.38119 60654.154 60654.154 -502.07711 -502.07711 6000 -22778.061 -22778.061 -22931.269 -22931.269 296.39023 296.39023 60633.757 60633.757 136.72911 136.72911 Loop time of 217.016 on 1 procs for 1000 steps with 4000 atoms Performance: 0.398 ns/day, 60.282 hours/ns, 4.608 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 | 215.93 | 215.93 | 215.93 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.1646 | 0.1646 | 0.1646 | 0.0 | 0.08 Output | 0.00027946 | 0.00027946 | 0.00027946 | 0.0 | 0.00 Modify | 0.81701 | 0.81701 | 0.81701 | 0.0 | 0.38 Other | | 0.1067 | | | 0.05 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: 310626.0 ave 310626 max 310626 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 310626 Ave neighs/atom = 77.656500 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 = 293.579279189468, Press = 49.0079147018937 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.4 ghost atom cutoff = 6.4 binsize = 3.2, bins = 13 13 13 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.4 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 6000 -22778.061 -22778.061 -22931.269 -22931.269 296.39023 296.39023 60633.757 60633.757 136.72911 136.72911 7000 -22775.928 -22775.928 -22927.187 -22927.187 292.62024 292.62024 60613.362 60613.362 1571.149 1571.149 Loop time of 210.29 on 1 procs for 1000 steps with 4000 atoms Performance: 0.411 ns/day, 58.414 hours/ns, 4.755 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 | 209.24 | 209.24 | 209.24 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16021 | 0.16021 | 0.16021 | 0.0 | 0.08 Output | 0.00022765 | 0.00022765 | 0.00022765 | 0.0 | 0.00 Modify | 0.7858 | 0.7858 | 0.7858 | 0.0 | 0.37 Other | | 0.104 | | | 0.05 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: 310808.0 ave 310808 max 310808 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 310808 Ave neighs/atom = 77.702000 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 = 293.010404725018, Press = 18.8566337240152 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.4 ghost atom cutoff = 6.4 binsize = 3.2, bins = 13 13 13 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.4 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 7000 -22775.928 -22775.928 -22927.187 -22927.187 292.62024 292.62024 60613.362 60613.362 1571.149 1571.149 8000 -22777.617 -22777.617 -22928.377 -22928.377 291.65423 291.65423 60638.935 60638.935 349.26894 349.26894 Loop time of 213.149 on 1 procs for 1000 steps with 4000 atoms Performance: 0.405 ns/day, 59.208 hours/ns, 4.692 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 | 212.08 | 212.08 | 212.08 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16179 | 0.16179 | 0.16179 | 0.0 | 0.08 Output | 0.00023738 | 0.00023738 | 0.00023738 | 0.0 | 0.00 Modify | 0.7998 | 0.7998 | 0.7998 | 0.0 | 0.38 Other | | 0.1041 | | | 0.05 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: 310614.0 ave 310614 max 310614 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 310614 Ave neighs/atom = 77.653500 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 = 293.339274649503, Press = 4.40988760113228 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.4 ghost atom cutoff = 6.4 binsize = 3.2, bins = 13 13 13 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.4 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 8000 -22777.617 -22777.617 -22928.377 -22928.377 291.65423 291.65423 60638.935 60638.935 349.26894 349.26894 9000 -22779.149 -22779.149 -22931.341 -22931.341 294.42597 294.42597 60623.183 60623.183 608.21352 608.21352 Loop time of 200.885 on 1 procs for 1000 steps with 4000 atoms Performance: 0.430 ns/day, 55.801 hours/ns, 4.978 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 | 199.88 | 199.88 | 199.88 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.15568 | 0.15568 | 0.15568 | 0.0 | 0.08 Output | 0.00018004 | 0.00018004 | 0.00018004 | 0.0 | 0.00 Modify | 0.74801 | 0.74801 | 0.74801 | 0.0 | 0.37 Other | | 0.1001 | | | 0.05 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: 310834.0 ave 310834 max 310834 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 310834 Ave neighs/atom = 77.708500 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 = 293.080789606182, Press = 18.4693861142647 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.4 ghost atom cutoff = 6.4 binsize = 3.2, bins = 13 13 13 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.4 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 9000 -22779.149 -22779.149 -22931.341 -22931.341 294.42597 294.42597 60623.183 60623.183 608.21352 608.21352 10000 -22776.631 -22776.631 -22927.23 -22927.23 291.34354 291.34354 60626.95 60626.95 958.33614 958.33614 Loop time of 218.457 on 1 procs for 1000 steps with 4000 atoms Performance: 0.396 ns/day, 60.682 hours/ns, 4.578 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 | 217.36 | 217.36 | 217.36 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16452 | 0.16452 | 0.16452 | 0.0 | 0.08 Output | 0.00023356 | 0.00023356 | 0.00023356 | 0.0 | 0.00 Modify | 0.82261 | 0.82261 | 0.82261 | 0.0 | 0.38 Other | | 0.1065 | | | 0.05 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: 310744.0 ave 310744 max 310744 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 310744 Ave neighs/atom = 77.686000 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 = 293.264855850157, Press = 9.73180655713151 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.4 ghost atom cutoff = 6.4 binsize = 3.2, bins = 13 13 13 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.4 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 10000 -22776.631 -22776.631 -22927.23 -22927.23 291.34354 291.34354 60626.95 60626.95 958.33614 958.33614 11000 -22777.835 -22777.835 -22930.387 -22930.387 295.12389 295.12389 60637.166 60637.166 161.07546 161.07546 Loop time of 218.477 on 1 procs for 1000 steps with 4000 atoms Performance: 0.395 ns/day, 60.688 hours/ns, 4.577 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 | 217.37 | 217.37 | 217.37 | 0.0 | 99.49 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16517 | 0.16517 | 0.16517 | 0.0 | 0.08 Output | 0.00023643 | 0.00023643 | 0.00023643 | 0.0 | 0.00 Modify | 0.83363 | 0.83363 | 0.83363 | 0.0 | 0.38 Other | | 0.1068 | | | 0.05 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: 310770.0 ave 310770 max 310770 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 310770 Ave neighs/atom = 77.692500 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 = 293.201582555162, Press = 0.566835596701931 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.4 ghost atom cutoff = 6.4 binsize = 3.2, bins = 13 13 13 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.4 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 11000 -22777.835 -22777.835 -22930.387 -22930.387 295.12389 295.12389 60637.166 60637.166 161.07546 161.07546 12000 -22779.839 -22779.839 -22928.394 -22928.394 287.38918 287.38918 60655.832 60655.832 -494.50422 -494.50422 Loop time of 215.694 on 1 procs for 1000 steps with 4000 atoms Performance: 0.401 ns/day, 59.915 hours/ns, 4.636 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 | 214.58 | 214.58 | 214.58 | 0.0 | 99.48 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.1672 | 0.1672 | 0.1672 | 0.0 | 0.08 Output | 0.00023159 | 0.00023159 | 0.00023159 | 0.0 | 0.00 Modify | 0.84151 | 0.84151 | 0.84151 | 0.0 | 0.39 Other | | 0.1041 | | | 0.05 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: 310648.0 ave 310648 max 310648 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 310648 Ave neighs/atom = 77.662000 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 = 293.068859461708, Press = 0.599578202415252 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.4 ghost atom cutoff = 6.4 binsize = 3.2, bins = 13 13 13 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.4 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 12000 -22779.839 -22779.839 -22928.394 -22928.394 287.38918 287.38918 60655.832 60655.832 -494.50422 -494.50422 13000 -22774.886 -22774.886 -22929.627 -22929.627 299.35739 299.35739 60666.711 60666.711 -1109.7223 -1109.7223 Loop time of 214.367 on 1 procs for 1000 steps with 4000 atoms Performance: 0.403 ns/day, 59.546 hours/ns, 4.665 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 | 213.29 | 213.29 | 213.29 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16336 | 0.16336 | 0.16336 | 0.0 | 0.08 Output | 0.0001787 | 0.0001787 | 0.0001787 | 0.0 | 0.00 Modify | 0.80547 | 0.80547 | 0.80547 | 0.0 | 0.38 Other | | 0.1047 | | | 0.05 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: 310716.0 ave 310716 max 310716 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 310716 Ave neighs/atom = 77.679000 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 = 292.891309609792, Press = 2.97002258551121 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.4 ghost atom cutoff = 6.4 binsize = 3.2, bins = 13 13 13 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.4 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 13000 -22774.886 -22774.886 -22929.627 -22929.627 299.35739 299.35739 60666.711 60666.711 -1109.7223 -1109.7223 14000 -22776.27 -22776.27 -22931.208 -22931.208 299.73849 299.73849 60645.314 60645.314 -332.79639 -332.79639 Loop time of 217.401 on 1 procs for 1000 steps with 4000 atoms Performance: 0.397 ns/day, 60.389 hours/ns, 4.600 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 | 216.32 | 216.32 | 216.32 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16301 | 0.16301 | 0.16301 | 0.0 | 0.07 Output | 0.00022277 | 0.00022277 | 0.00022277 | 0.0 | 0.00 Modify | 0.81615 | 0.81615 | 0.81615 | 0.0 | 0.38 Other | | 0.1061 | | | 0.05 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: 310690.0 ave 310690 max 310690 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 310690 Ave neighs/atom = 77.672500 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 = 292.921586030265, Press = 3.43257036551497 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.4 ghost atom cutoff = 6.4 binsize = 3.2, bins = 13 13 13 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.4 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 14000 -22776.27 -22776.27 -22931.208 -22931.208 299.73849 299.73849 60645.314 60645.314 -332.79639 -332.79639 15000 -22775.583 -22775.583 -22927.397 -22927.397 293.69599 293.69599 60637.402 60637.402 576.57348 576.57348 Loop time of 211.044 on 1 procs for 1000 steps with 4000 atoms Performance: 0.409 ns/day, 58.623 hours/ns, 4.738 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 | 210 | 210 | 210 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16001 | 0.16001 | 0.16001 | 0.0 | 0.08 Output | 0.00048406 | 0.00048406 | 0.00048406 | 0.0 | 0.00 Modify | 0.78263 | 0.78263 | 0.78263 | 0.0 | 0.37 Other | | 0.1027 | | | 0.05 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: 310798.0 ave 310798 max 310798 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 310798 Ave neighs/atom = 77.699500 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 = 292.939935668181, Press = 3.53908965914714 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.4 ghost atom cutoff = 6.4 binsize = 3.2, bins = 13 13 13 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.4 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 15000 -22775.583 -22775.583 -22927.397 -22927.397 293.69599 293.69599 60637.402 60637.402 576.57348 576.57348 16000 -22779.251 -22779.251 -22928.895 -22928.895 289.49736 289.49736 60595.436 60595.436 2176.3283 2176.3283 Loop time of 205.85 on 1 procs for 1000 steps with 4000 atoms Performance: 0.420 ns/day, 57.180 hours/ns, 4.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 | 204.82 | 204.82 | 204.82 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.15989 | 0.15989 | 0.15989 | 0.0 | 0.08 Output | 0.00023357 | 0.00023357 | 0.00023357 | 0.0 | 0.00 Modify | 0.7704 | 0.7704 | 0.7704 | 0.0 | 0.37 Other | | 0.1027 | | | 0.05 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: 310682.0 ave 310682 max 310682 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 310682 Ave neighs/atom = 77.670500 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 = 292.92658036828, Press = 4.92897282391379 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.4 ghost atom cutoff = 6.4 binsize = 3.2, bins = 13 13 13 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.4 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 16000 -22779.251 -22779.251 -22928.895 -22928.895 289.49736 289.49736 60595.436 60595.436 2176.3283 2176.3283 17000 -22778.58 -22778.58 -22931.132 -22931.132 295.12194 295.12194 60591.148 60591.148 2070.777 2070.777 Loop time of 213.741 on 1 procs for 1000 steps with 4000 atoms Performance: 0.404 ns/day, 59.373 hours/ns, 4.679 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 | 212.64 | 212.64 | 212.64 | 0.0 | 99.48 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16637 | 0.16637 | 0.16637 | 0.0 | 0.08 Output | 0.00028564 | 0.00028564 | 0.00028564 | 0.0 | 0.00 Modify | 0.83168 | 0.83168 | 0.83168 | 0.0 | 0.39 Other | | 0.1051 | | | 0.05 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: 310716.0 ave 310716 max 310716 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 310716 Ave neighs/atom = 77.679000 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 = 292.794582093889, Press = 2.53368678697106 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.4 ghost atom cutoff = 6.4 binsize = 3.2, bins = 13 13 13 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.4 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 17000 -22778.58 -22778.58 -22931.132 -22931.132 295.12194 295.12194 60591.148 60591.148 2070.777 2070.777 18000 -22777.039 -22777.039 -22927.026 -22927.026 290.15929 290.15929 60591.524 60591.524 2654.0412 2654.0412 Loop time of 213.183 on 1 procs for 1000 steps with 4000 atoms Performance: 0.405 ns/day, 59.218 hours/ns, 4.691 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 | 212.12 | 212.12 | 212.12 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16009 | 0.16009 | 0.16009 | 0.0 | 0.08 Output | 0.00017918 | 0.00017918 | 0.00017918 | 0.0 | 0.00 Modify | 0.79922 | 0.79922 | 0.79922 | 0.0 | 0.37 Other | | 0.105 | | | 0.05 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: 310820.0 ave 310820 max 310820 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 310820 Ave neighs/atom = 77.705000 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 = 292.815738265581, Press = -0.158588950700413 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.4 ghost atom cutoff = 6.4 binsize = 3.2, bins = 13 13 13 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.4 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 18000 -22777.039 -22777.039 -22927.026 -22927.026 290.15929 290.15929 60591.524 60591.524 2654.0412 2654.0412 19000 -22781.465 -22781.465 -22930.657 -22930.657 288.62087 288.62087 60636.345 60636.345 85.824909 85.824909 Loop time of 202.849 on 1 procs for 1000 steps with 4000 atoms Performance: 0.426 ns/day, 56.347 hours/ns, 4.930 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 | 201.84 | 201.84 | 201.84 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.15735 | 0.15735 | 0.15735 | 0.0 | 0.08 Output | 0.0001786 | 0.0001786 | 0.0001786 | 0.0 | 0.00 Modify | 0.75572 | 0.75572 | 0.75572 | 0.0 | 0.37 Other | | 0.09914 | | | 0.05 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: 310706.0 ave 310706 max 310706 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 310706 Ave neighs/atom = 77.676500 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 = 292.705432959829, Press = 3.06418396477821 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.4 ghost atom cutoff = 6.4 binsize = 3.2, bins = 13 13 13 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.4 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 19000 -22781.465 -22781.465 -22930.657 -22930.657 288.62087 288.62087 60636.345 60636.345 85.824909 85.824909 20000 -22780.136 -22780.136 -22930.175 -22930.175 290.25997 290.25997 60626.259 60626.259 624.65473 624.65473 Loop time of 186.171 on 1 procs for 1000 steps with 4000 atoms Performance: 0.464 ns/day, 51.714 hours/ns, 5.371 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 | 185.27 | 185.27 | 185.27 | 0.0 | 99.51 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.14599 | 0.14599 | 0.14599 | 0.0 | 0.08 Output | 0.00018234 | 0.00018234 | 0.00018234 | 0.0 | 0.00 Modify | 0.66504 | 0.66504 | 0.66504 | 0.0 | 0.36 Other | | 0.09393 | | | 0.05 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: 310904.0 ave 310904 max 310904 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 310904 Ave neighs/atom = 77.726000 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 = 292.718427135571, Press = 1.16820110132623 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.4 ghost atom cutoff = 6.4 binsize = 3.2, bins = 13 13 13 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.4 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 20000 -22780.136 -22780.136 -22930.175 -22930.175 290.25997 290.25997 60626.259 60626.259 624.65473 624.65473 21000 -22774.513 -22774.513 -22929.226 -22929.226 299.30315 299.30315 60641.072 60641.072 -38.987651 -38.987651 Loop time of 195.999 on 1 procs for 1000 steps with 4000 atoms Performance: 0.441 ns/day, 54.444 hours/ns, 5.102 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 | 195.03 | 195.03 | 195.03 | 0.0 | 99.51 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.15154 | 0.15154 | 0.15154 | 0.0 | 0.08 Output | 0.00024989 | 0.00024989 | 0.00024989 | 0.0 | 0.00 Modify | 0.71492 | 0.71492 | 0.71492 | 0.0 | 0.36 Other | | 0.09804 | | | 0.05 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: 310766.0 ave 310766 max 310766 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 310766 Ave neighs/atom = 77.691500 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 = 292.779136735021, Press = 4.32670311309112 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.4 ghost atom cutoff = 6.4 binsize = 3.2, bins = 13 13 13 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.4 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 21000 -22774.513 -22774.513 -22929.226 -22929.226 299.30315 299.30315 60641.072 60641.072 -38.987651 -38.987651 22000 -22779.406 -22779.406 -22928.613 -22928.613 288.65062 288.65062 60636.019 60636.019 266.83043 266.83043 Loop time of 214.328 on 1 procs for 1000 steps with 4000 atoms Performance: 0.403 ns/day, 59.536 hours/ns, 4.666 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 | 213.24 | 213.24 | 213.24 | 0.0 | 99.49 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16388 | 0.16388 | 0.16388 | 0.0 | 0.08 Output | 0.00024909 | 0.00024909 | 0.00024909 | 0.0 | 0.00 Modify | 0.81506 | 0.81506 | 0.81506 | 0.0 | 0.38 Other | | 0.1059 | | | 0.05 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: 310764.0 ave 310764 max 310764 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 310764 Ave neighs/atom = 77.691000 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 = 292.817680475795, Press = 3.14644076316482 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.4 ghost atom cutoff = 6.4 binsize = 3.2, bins = 13 13 13 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.4 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 22000 -22779.406 -22779.406 -22928.613 -22928.613 288.65062 288.65062 60636.019 60636.019 266.83043 266.83043 23000 -22776.738 -22776.738 -22929.817 -22929.817 296.14263 296.14263 60663.176 60663.176 -981.85477 -981.85477 Loop time of 206.578 on 1 procs for 1000 steps with 4000 atoms Performance: 0.418 ns/day, 57.383 hours/ns, 4.841 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 | 205.55 | 205.55 | 205.55 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.15814 | 0.15814 | 0.15814 | 0.0 | 0.08 Output | 0.00021874 | 0.00021874 | 0.00021874 | 0.0 | 0.00 Modify | 0.766 | 0.766 | 0.766 | 0.0 | 0.37 Other | | 0.1016 | | | 0.05 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: 310746.0 ave 310746 max 310746 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 310746 Ave neighs/atom = 77.686500 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 = 292.851527651228, Press = 2.26759447914331 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.4 ghost atom cutoff = 6.4 binsize = 3.2, bins = 13 13 13 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.4 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 23000 -22776.738 -22776.738 -22929.817 -22929.817 296.14263 296.14263 60663.176 60663.176 -981.85477 -981.85477 24000 -22778.925 -22778.925 -22928.242 -22928.242 288.86382 288.86382 60654.396 60654.396 -536.30291 -536.30291 Loop time of 207.904 on 1 procs for 1000 steps with 4000 atoms Performance: 0.416 ns/day, 57.751 hours/ns, 4.810 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 | 206.87 | 206.87 | 206.87 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.15902 | 0.15902 | 0.15902 | 0.0 | 0.08 Output | 0.00017973 | 0.00017973 | 0.00017973 | 0.0 | 0.00 Modify | 0.77146 | 0.77146 | 0.77146 | 0.0 | 0.37 Other | | 0.1027 | | | 0.05 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: 310788.0 ave 310788 max 310788 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 310788 Ave neighs/atom = 77.697000 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 = 292.93789422705, Press = 0.597212290774878 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.4 ghost atom cutoff = 6.4 binsize = 3.2, bins = 13 13 13 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.4 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 24000 -22778.925 -22778.925 -22928.242 -22928.242 288.86382 288.86382 60654.396 60654.396 -536.30291 -536.30291 25000 -22778.433 -22778.433 -22929.108 -22929.108 291.49018 291.49018 60607.256 60607.256 1670.4793 1670.4793 Loop time of 219.007 on 1 procs for 1000 steps with 4000 atoms Performance: 0.395 ns/day, 60.835 hours/ns, 4.566 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 | 217.91 | 217.91 | 217.91 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16575 | 0.16575 | 0.16575 | 0.0 | 0.08 Output | 0.00027871 | 0.00027871 | 0.00027871 | 0.0 | 0.00 Modify | 0.82702 | 0.82702 | 0.82702 | 0.0 | 0.38 Other | | 0.1077 | | | 0.05 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: 310648.0 ave 310648 max 310648 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 310648 Ave neighs/atom = 77.662000 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 = 292.944723875393, Press = 1.64073838678936 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.4 ghost atom cutoff = 6.4 binsize = 3.2, bins = 13 13 13 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.4 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 25000 -22778.433 -22778.433 -22929.108 -22929.108 291.49018 291.49018 60607.256 60607.256 1670.4793 1670.4793 26000 -22777.905 -22777.905 -22930.446 -22930.446 295.10146 295.10146 60638.568 60638.568 10.63213 10.63213 Loop time of 212.318 on 1 procs for 1000 steps with 4000 atoms Performance: 0.407 ns/day, 58.977 hours/ns, 4.710 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 | 211.25 | 211.25 | 211.25 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16307 | 0.16307 | 0.16307 | 0.0 | 0.08 Output | 0.00017889 | 0.00017889 | 0.00017889 | 0.0 | 0.00 Modify | 0.80023 | 0.80023 | 0.80023 | 0.0 | 0.38 Other | | 0.1045 | | | 0.05 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: 310784.0 ave 310784 max 310784 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 310784 Ave neighs/atom = 77.696000 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 = 292.971474313986, Press = -0.252816732900898 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.4 ghost atom cutoff = 6.4 binsize = 3.2, bins = 13 13 13 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.4 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 26000 -22777.905 -22777.905 -22930.446 -22930.446 295.10146 295.10146 60638.568 60638.568 10.63213 10.63213 27000 -22774.72 -22774.72 -22926.964 -22926.964 294.52586 294.52586 60614.475 60614.475 1638.5934 1638.5934 Loop time of 196.777 on 1 procs for 1000 steps with 4000 atoms Performance: 0.439 ns/day, 54.660 hours/ns, 5.082 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 | 195.8 | 195.8 | 195.8 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.15533 | 0.15533 | 0.15533 | 0.0 | 0.08 Output | 0.00022068 | 0.00022068 | 0.00022068 | 0.0 | 0.00 Modify | 0.72189 | 0.72189 | 0.72189 | 0.0 | 0.37 Other | | 0.09685 | | | 0.05 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: 310848.0 ave 310848 max 310848 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 310848 Ave neighs/atom = 77.712000 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 = 293.066028396067, Press = -0.218742281970894 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.4 ghost atom cutoff = 6.4 binsize = 3.2, bins = 13 13 13 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.4 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 27000 -22774.72 -22774.72 -22926.964 -22926.964 294.52586 294.52586 60614.475 60614.475 1638.5934 1638.5934 28000 -22781.798 -22781.798 -22929.163 -22929.163 285.08879 285.08879 60631.835 60631.835 364.27991 364.27991 Loop time of 190.722 on 1 procs for 1000 steps with 4000 atoms Performance: 0.453 ns/day, 52.978 hours/ns, 5.243 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 | 189.78 | 189.78 | 189.78 | 0.0 | 99.51 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.15124 | 0.15124 | 0.15124 | 0.0 | 0.08 Output | 0.00018356 | 0.00018356 | 0.00018356 | 0.0 | 0.00 Modify | 0.69088 | 0.69088 | 0.69088 | 0.0 | 0.36 Other | | 0.09548 | | | 0.05 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: 310732.0 ave 310732 max 310732 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 310732 Ave neighs/atom = 77.683000 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 = 293.03982519451, Press = 0.589428663828067 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.4 ghost atom cutoff = 6.4 binsize = 3.2, bins = 13 13 13 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.4 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 28000 -22781.798 -22781.798 -22929.163 -22929.163 285.08879 285.08879 60631.835 60631.835 364.27991 364.27991 29000 -22777.449 -22777.449 -22930.35 -22930.35 295.79882 295.79882 60650.714 60650.714 -510.70522 -510.70522 Loop time of 185.535 on 1 procs for 1000 steps with 4000 atoms Performance: 0.466 ns/day, 51.537 hours/ns, 5.390 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 | 184.63 | 184.63 | 184.63 | 0.0 | 99.51 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.14921 | 0.14921 | 0.14921 | 0.0 | 0.08 Output | 0.00017901 | 0.00017901 | 0.00017901 | 0.0 | 0.00 Modify | 0.66644 | 0.66644 | 0.66644 | 0.0 | 0.36 Other | | 0.0919 | | | 0.05 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: 310792.0 ave 310792 max 310792 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 310792 Ave neighs/atom = 77.698000 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 = 293.032579916058, Press = -0.885503670114265 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.4 ghost atom cutoff = 6.4 binsize = 3.2, bins = 13 13 13 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.4 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 29000 -22777.449 -22777.449 -22930.35 -22930.35 295.79882 295.79882 60650.714 60650.714 -510.70522 -510.70522 30000 -22780.047 -22780.047 -22930.177 -22930.177 290.43628 290.43628 60659.398 60659.398 -994.06043 -994.06043 Loop time of 186.291 on 1 procs for 1000 steps with 4000 atoms Performance: 0.464 ns/day, 51.747 hours/ns, 5.368 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 | 185.38 | 185.38 | 185.38 | 0.0 | 99.51 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.14685 | 0.14685 | 0.14685 | 0.0 | 0.08 Output | 0.00022258 | 0.00022258 | 0.00022258 | 0.0 | 0.00 Modify | 0.66691 | 0.66691 | 0.66691 | 0.0 | 0.36 Other | | 0.09441 | | | 0.05 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: 310768.0 ave 310768 max 310768 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 310768 Ave neighs/atom = 77.692000 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 = 293.100064140074, Press = -0.511801029041529 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.4 ghost atom cutoff = 6.4 binsize = 3.2, bins = 13 13 13 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.4 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 30000 -22780.047 -22780.047 -22930.177 -22930.177 290.43628 290.43628 60659.398 60659.398 -994.06043 -994.06043 31000 -22778.362 -22778.362 -22930.102 -22930.102 293.55084 293.55084 60649.405 60649.405 -573.08207 -573.08207 Loop time of 197.971 on 1 procs for 1000 steps with 4000 atoms Performance: 0.436 ns/day, 54.992 hours/ns, 5.051 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 | 196.99 | 196.99 | 196.99 | 0.0 | 99.51 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.15068 | 0.15068 | 0.15068 | 0.0 | 0.08 Output | 0.00017821 | 0.00017821 | 0.00017821 | 0.0 | 0.00 Modify | 0.72627 | 0.72627 | 0.72627 | 0.0 | 0.37 Other | | 0.1001 | | | 0.05 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: 310702.0 ave 310702 max 310702 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 310702 Ave neighs/atom = 77.675500 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_T293.15.out" else "print 'not_converged' file output/vol_T293.15.out" print '${V}' file output/vol_T293.15.out 60641.8339459364 print "LAMMPS calculation completed" LAMMPS calculation completed quit 0