# periodic boundary conditions along all three dimensions boundary p p p # Set neighbor skin variable neigh_skin equal 2.0*${_u_distance} variable neigh_skin equal 2.0*1 neighbor ${neigh_skin} bin neighbor 2 bin # create a supercell with cubic lattice (fcc, bcc, sc, or diamond) # using 10*10*10 conventional (orthogonal) unit cells variable latticeconst_converted equal 5.658032611012459*${_u_distance} variable latticeconst_converted equal 5.658032611012459*1 lattice diamond ${latticeconst_converted} lattice diamond 5.65803261101246 Lattice spacing in x,y,z = 5.6580326 5.6580326 5.6580326 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 (56.580326 56.580326 56.580326) 1 by 1 by 1 MPI processor grid create_atoms 1 box Created 8000 atoms using lattice units in orthogonal box = (0.0000000 0.0000000 0.0000000) to (56.580326 56.580326 56.580326) create_atoms CPU = 0.006 seconds variable mass_converted equal 72.64*${_u_mass} variable mass_converted equal 72.64*1 kim_interactions Ge WARNING: 'kim_' has been renamed to 'kim '. Please update your input. kim interactions Ge #=== BEGIN kim interactions ================================== pair_style kim MEAM_LAMMPS_KimShinLee_2008_Ge__MO_657096500078_000 pair_coeff * * Ge #=== END kim interactions ==================================== mass 1 ${mass_converted} mass 1 72.64 # initial volume variable v equal vol # assign formula variable V0 equal ${v} # evaluate initial value variable V0 equal 181132.482255555 variable V0_metal equal ${V0}/(${_u_distance}*${_u_distance}*${_u_distance}) variable V0_metal equal 181132.482255555/(${_u_distance}*${_u_distance}*${_u_distance}) variable V0_metal equal 181132.482255555/(1*${_u_distance}*${_u_distance}) variable V0_metal equal 181132.482255555/(1*1*${_u_distance}) variable V0_metal equal 181132.482255555/(1*1*1) variable V0_metal_times1000 equal ${V0_metal}*1000 variable V0_metal_times1000 equal 181132.482255555*1000 print "Initial system volume: ${V0_metal} Angstroms^3" Initial system volume: 181132.482255555 Angstroms^3 # set the time step to 0.001 picoseconds variable timestep_converted equal 0.001*${_u_time} variable timestep_converted equal 0.001*1 timestep ${timestep_converted} timestep 0.001 variable temp_converted equal 253.15*${_u_temperature} variable temp_converted equal 253.15*1 variable Tdamp_converted equal 0.1*${_u_time} variable Tdamp_converted equal 0.1*1 variable press_converted equal 0.0*${_u_pressure} variable press_converted equal 0.0*1 variable Pdamp_converted equal 1*${_u_time} variable Pdamp_converted equal 1*1 # create initial velocities consistent with the chosen temperature velocity all create ${temp_converted} 17 mom yes rot yes velocity all create 253.15 17 mom yes rot yes # set NPT ensemble for all atoms fix ensemble all npt temp ${temp_converted} ${temp_converted} ${Tdamp_converted} iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 253.15 ${temp_converted} ${Tdamp_converted} iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 253.15 253.15 ${Tdamp_converted} iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 253.15 253.15 0.1 iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 253.15 253.15 0.1 iso 0 ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 253.15 253.15 0.1 iso 0 0 ${Pdamp_converted} fix ensemble all npt temp 253.15 253.15 0.1 iso 0 0 1 # compute the time averages of pressure, temperature, and volume, respectively # ignore the first 5000 timesteps variable etotal_metal equal etotal/${_u_energy} variable etotal_metal equal etotal/1 variable pe_metal equal pe/${_u_energy} variable pe_metal equal pe/1 variable T_metal equal temp/${_u_temperature} variable T_metal equal temp/1 variable V_metal equal vol/(${_u_distance}*${_u_distance}*${_u_distance}) variable V_metal equal vol/(1*${_u_distance}*${_u_distance}) variable V_metal equal vol/(1*1*${_u_distance}) variable V_metal equal vol/(1*1*1) variable P_metal equal press/${_u_pressure} variable P_metal equal press/1 fix avgmyTemp all ave/time 5 20 100 v_T_metal ave running start 5000 fix avgmyPress all ave/time 5 20 100 v_P_metal ave running start 5000 fix avgmyVol all ave/time 5 20 100 v_V_metal ave running start 5000 # extract fix quantities into variables so they can be used in if-else logic later. variable T equal f_avgmyTemp variable P equal f_avgmyPress variable V equal f_avgmyVol # set error bounds for temperature and pressure in original metal units (K and bar) variable T_low equal "253.15 - 0.2" variable T_up equal "253.15 + 0.2" variable P_low equal "0.0 - 0.2" variable P_up equal "0.0 + 0.2" # print to logfile every 1000 timesteps thermo_style custom step etotal v_etotal_metal pe v_pe_metal temp v_T_metal vol v_V_metal press v_P_metal thermo 1000 # Run a simulation for at most 2000*1000 timesteps. At each 1000th time step, check # whether the temperature and pressure have converged. If yes, break. label top variable a loop 2000 run 1000 CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE Your simulation uses code contributions which should be cited: - OpenKIM: https://doi.org/10.1007/s11837-011-0102-6 @Article{tadmor:elliott:2011, author = {E. B. Tadmor and R. S. Elliott and J. P. Sethna and R. E. Miller and C. A. Becker}, title = {The potential of atomistic simulations and the {K}nowledgebase of {I}nteratomic {M}odels}, journal = {{JOM}}, year = 2011, volume = 63, number = 17, pages = {17}, doi = {10.1007/s11837-011-0102-6} } - OpenKIM potential: https://openkim.org/cite/MO_657096500078_000#item-citation CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 18 18 18 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.683 | 4.683 | 4.683 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 0 -30538.255 -30538.255 -30800 -30800 253.15 253.15 181132.48 181132.48 1543.4804 1543.4804 1000 -30254.631 -30254.631 -30530.752 -30530.752 267.0544 267.0544 182129.85 182129.85 -49.723285 -49.723285 Loop time of 125.576 on 1 procs for 1000 steps with 8000 atoms Performance: 0.688 ns/day, 34.882 hours/ns, 7.963 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 | 123.82 | 123.82 | 123.82 | 0.0 | 98.61 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.21491 | 0.21491 | 0.21491 | 0.0 | 0.17 Output | 0.00027756 | 0.00027756 | 0.00027756 | 0.0 | 0.00 Modify | 1.3727 | 1.3727 | 1.3727 | 0.0 | 1.09 Other | | 0.1635 | | | 0.13 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 368000.0 ave 368000 max 368000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 368000 Ave neighs/atom = 46.000000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 0, Press = 0 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 18 18 18 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.684 | 4.684 | 4.684 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 1000 -30254.631 -30254.631 -30530.752 -30530.752 267.0544 267.0544 182129.85 182129.85 -49.723285 -49.723285 2000 -30281.629 -30281.629 -30550.472 -30550.472 260.01454 260.01454 181990.84 181990.84 450.54102 450.54102 Loop time of 122.088 on 1 procs for 1000 steps with 8000 atoms Performance: 0.708 ns/day, 33.913 hours/ns, 8.191 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 | 120.42 | 120.42 | 120.42 | 0.0 | 98.63 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.2117 | 0.2117 | 0.2117 | 0.0 | 0.17 Output | 0.00027859 | 0.00027859 | 0.00027859 | 0.0 | 0.00 Modify | 1.3023 | 1.3023 | 1.3023 | 0.0 | 1.07 Other | | 0.1585 | | | 0.13 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 367998.0 ave 367998 max 367998 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 367998 Ave neighs/atom = 45.999750 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 0, Press = 0 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 18 18 18 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.684 | 4.684 | 4.684 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 2000 -30281.629 -30281.629 -30550.472 -30550.472 260.01454 260.01454 181990.84 181990.84 450.54102 450.54102 3000 -30267.285 -30267.285 -30536.275 -30536.275 260.15684 260.15684 182299.41 182299.41 -687.66585 -687.66585 Loop time of 122.082 on 1 procs for 1000 steps with 8000 atoms Performance: 0.708 ns/day, 33.912 hours/ns, 8.191 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 | 120.41 | 120.41 | 120.41 | 0.0 | 98.63 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.21204 | 0.21204 | 0.21204 | 0.0 | 0.17 Output | 0.00026994 | 0.00026994 | 0.00026994 | 0.0 | 0.00 Modify | 1.3026 | 1.3026 | 1.3026 | 0.0 | 1.07 Other | | 0.1596 | | | 0.13 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 367992.0 ave 367992 max 367992 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 367992 Ave neighs/atom = 45.999000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 0, Press = 0 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 18 18 18 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.684 | 4.684 | 4.684 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 3000 -30267.285 -30267.285 -30536.275 -30536.275 260.15684 260.15684 182299.41 182299.41 -687.66585 -687.66585 4000 -30276.589 -30276.589 -30539.626 -30539.626 254.39907 254.39907 182241.69 182241.69 -448.75433 -448.75433 Loop time of 122.242 on 1 procs for 1000 steps with 8000 atoms Performance: 0.707 ns/day, 33.956 hours/ns, 8.181 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 | 120.57 | 120.57 | 120.57 | 0.0 | 98.63 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.21292 | 0.21292 | 0.21292 | 0.0 | 0.17 Output | 0.00031153 | 0.00031153 | 0.00031153 | 0.0 | 0.00 Modify | 1.2999 | 1.2999 | 1.2999 | 0.0 | 1.06 Other | | 0.159 | | | 0.13 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 367996.0 ave 367996 max 367996 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 367996 Ave neighs/atom = 45.999500 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 0, Press = 0 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 18 18 18 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.684 | 4.684 | 4.684 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 4000 -30276.589 -30276.589 -30539.626 -30539.626 254.39907 254.39907 182241.69 182241.69 -448.75433 -448.75433 5000 -30272.334 -30272.334 -30527.338 -30527.338 246.63058 246.63058 182305.94 182305.94 -785.03168 -785.03168 Loop time of 122.303 on 1 procs for 1000 steps with 8000 atoms Performance: 0.706 ns/day, 33.973 hours/ns, 8.176 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 | 120.62 | 120.62 | 120.62 | 0.0 | 98.63 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.21174 | 0.21174 | 0.21174 | 0.0 | 0.17 Output | 0.00027861 | 0.00027861 | 0.00027861 | 0.0 | 0.00 Modify | 1.3077 | 1.3077 | 1.3077 | 0.0 | 1.07 Other | | 0.1594 | | | 0.13 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 368000.0 ave 368000 max 368000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 368000 Ave neighs/atom = 46.000000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 257.047983002139, Press = 196.114143078798 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 18 18 18 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.684 | 4.684 | 4.684 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 5000 -30272.334 -30272.334 -30527.338 -30527.338 246.63058 246.63058 182305.94 182305.94 -785.03168 -785.03168 6000 -30272.334 -30272.334 -30531.885 -30531.885 251.02828 251.02828 182317.95 182317.95 -747.18488 -747.18488 Loop time of 122.124 on 1 procs for 1000 steps with 8000 atoms Performance: 0.707 ns/day, 33.923 hours/ns, 8.188 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 | 120.43 | 120.43 | 120.43 | 0.0 | 98.61 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.20927 | 0.20927 | 0.20927 | 0.0 | 0.17 Output | 0.00024149 | 0.00024149 | 0.00024149 | 0.0 | 0.00 Modify | 1.3288 | 1.3288 | 1.3288 | 0.0 | 1.09 Other | | 0.1593 | | | 0.13 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 367996.0 ave 367996 max 367996 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 367996 Ave neighs/atom = 45.999500 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 253.093125472314, Press = 5.8207204466941 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 18 18 18 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.684 | 4.684 | 4.684 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 6000 -30272.334 -30272.334 -30531.885 -30531.885 251.02828 251.02828 182317.95 182317.95 -747.18488 -747.18488 7000 -30274.106 -30274.106 -30532.249 -30532.249 249.66632 249.66632 182269.04 182269.04 -562.99484 -562.99484 Loop time of 122.134 on 1 procs for 1000 steps with 8000 atoms Performance: 0.707 ns/day, 33.926 hours/ns, 8.188 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 | 120.44 | 120.44 | 120.44 | 0.0 | 98.61 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.21116 | 0.21116 | 0.21116 | 0.0 | 0.17 Output | 0.00023466 | 0.00023466 | 0.00023466 | 0.0 | 0.00 Modify | 1.3267 | 1.3267 | 1.3267 | 0.0 | 1.09 Other | | 0.1602 | | | 0.13 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 368000.0 ave 368000 max 368000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 368000 Ave neighs/atom = 46.000000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 253.333856154108, Press = 1.75685205890609 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 18 18 18 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.684 | 4.684 | 4.684 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 7000 -30274.106 -30274.106 -30532.249 -30532.249 249.66632 249.66632 182269.04 182269.04 -562.99484 -562.99484 8000 -30273.018 -30273.018 -30531.522 -30531.522 250.01541 250.01541 182221.25 182221.25 -335.82058 -335.82058 Loop time of 122.148 on 1 procs for 1000 steps with 8000 atoms Performance: 0.707 ns/day, 33.930 hours/ns, 8.187 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 | 120.45 | 120.45 | 120.45 | 0.0 | 98.61 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.21177 | 0.21177 | 0.21177 | 0.0 | 0.17 Output | 0.00022997 | 0.00022997 | 0.00022997 | 0.0 | 0.00 Modify | 1.3281 | 1.3281 | 1.3281 | 0.0 | 1.09 Other | | 0.1584 | | | 0.13 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 367996.0 ave 367996 max 367996 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 367996 Ave neighs/atom = 45.999500 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 253.206367491911, Press = -0.622765691278274 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 18 18 18 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.684 | 4.684 | 4.684 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 8000 -30273.018 -30273.018 -30531.522 -30531.522 250.01541 250.01541 182221.25 182221.25 -335.82058 -335.82058 9000 -30274.735 -30274.735 -30532.454 -30532.454 249.25649 249.25649 182243.58 182243.58 -392.61023 -392.61023 Loop time of 122.146 on 1 procs for 1000 steps with 8000 atoms Performance: 0.707 ns/day, 33.929 hours/ns, 8.187 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 | 120.45 | 120.45 | 120.45 | 0.0 | 98.61 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.21 | 0.21 | 0.21 | 0.0 | 0.17 Output | 0.00024517 | 0.00024517 | 0.00024517 | 0.0 | 0.00 Modify | 1.3274 | 1.3274 | 1.3274 | 0.0 | 1.09 Other | | 0.1594 | | | 0.13 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 368000.0 ave 368000 max 368000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 368000 Ave neighs/atom = 46.000000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 253.153812752312, Press = 1.23497450046102 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 18 18 18 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.684 | 4.684 | 4.684 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 9000 -30274.735 -30274.735 -30532.454 -30532.454 249.25649 249.25649 182243.58 182243.58 -392.61023 -392.61023 10000 -30272.944 -30272.944 -30531.357 -30531.357 249.92701 249.92701 182049.64 182049.64 369.13079 369.13079 Loop time of 122.311 on 1 procs for 1000 steps with 8000 atoms Performance: 0.706 ns/day, 33.975 hours/ns, 8.176 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 | 120.61 | 120.61 | 120.61 | 0.0 | 98.61 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.2088 | 0.2088 | 0.2088 | 0.0 | 0.17 Output | 0.00023341 | 0.00023341 | 0.00023341 | 0.0 | 0.00 Modify | 1.3306 | 1.3306 | 1.3306 | 0.0 | 1.09 Other | | 0.1583 | | | 0.13 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 367998.0 ave 367998 max 367998 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 367998 Ave neighs/atom = 45.999750 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 253.139114147614, Press = 0.215895440711275 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 18 18 18 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.684 | 4.684 | 4.684 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 10000 -30272.944 -30272.944 -30531.357 -30531.357 249.92701 249.92701 182049.64 182049.64 369.13079 369.13079 11000 -30273.596 -30273.596 -30533.609 -30533.609 251.47583 251.47583 182077.21 182077.21 230.14956 230.14956 Loop time of 122.282 on 1 procs for 1000 steps with 8000 atoms Performance: 0.707 ns/day, 33.967 hours/ns, 8.178 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 | 120.58 | 120.58 | 120.58 | 0.0 | 98.61 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.21073 | 0.21073 | 0.21073 | 0.0 | 0.17 Output | 0.00023217 | 0.00023217 | 0.00023217 | 0.0 | 0.00 Modify | 1.3284 | 1.3284 | 1.3284 | 0.0 | 1.09 Other | | 0.1594 | | | 0.13 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 368000.0 ave 368000 max 368000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 368000 Ave neighs/atom = 46.000000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 253.230839453383, Press = 1.75365987817893 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 18 18 18 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.684 | 4.684 | 4.684 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 11000 -30273.596 -30273.596 -30533.609 -30533.609 251.47583 251.47583 182077.21 182077.21 230.14956 230.14956 12000 -30276.325 -30276.325 -30535.604 -30535.604 250.76531 250.76531 181974.29 181974.29 605.3878 605.3878 Loop time of 122.299 on 1 procs for 1000 steps with 8000 atoms Performance: 0.706 ns/day, 33.972 hours/ns, 8.177 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 | 120.6 | 120.6 | 120.6 | 0.0 | 98.61 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.21126 | 0.21126 | 0.21126 | 0.0 | 0.17 Output | 0.00024156 | 0.00024156 | 0.00024156 | 0.0 | 0.00 Modify | 1.3296 | 1.3296 | 1.3296 | 0.0 | 1.09 Other | | 0.1589 | | | 0.13 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 367998.0 ave 367998 max 367998 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 367998 Ave neighs/atom = 45.999750 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 253.1993000444, Press = -0.775437123953161 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 18 18 18 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.684 | 4.684 | 4.684 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 12000 -30276.325 -30276.325 -30535.604 -30535.604 250.76531 250.76531 181974.29 181974.29 605.3878 605.3878 13000 -30275.905 -30275.905 -30530.684 -30530.684 246.41262 246.41262 181966.86 181966.86 601.10305 601.10305 Loop time of 122.281 on 1 procs for 1000 steps with 8000 atoms Performance: 0.707 ns/day, 33.967 hours/ns, 8.178 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 | 120.58 | 120.58 | 120.58 | 0.0 | 98.61 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.21157 | 0.21157 | 0.21157 | 0.0 | 0.17 Output | 0.00023416 | 0.00023416 | 0.00023416 | 0.0 | 0.00 Modify | 1.3288 | 1.3288 | 1.3288 | 0.0 | 1.09 Other | | 0.1587 | | | 0.13 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 367998.0 ave 367998 max 367998 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 367998 Ave neighs/atom = 45.999750 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 253.115783086846, Press = 0.348315327520334 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 18 18 18 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.684 | 4.684 | 4.684 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 13000 -30275.905 -30275.905 -30530.684 -30530.684 246.41262 246.41262 181966.86 181966.86 601.10305 601.10305 14000 -30274.282 -30274.282 -30536.371 -30536.371 253.48208 253.48208 181983.31 181983.31 612.65246 612.65246 Loop time of 120.989 on 1 procs for 1000 steps with 8000 atoms Performance: 0.714 ns/day, 33.608 hours/ns, 8.265 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 | 119.32 | 119.32 | 119.32 | 0.0 | 98.62 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.19965 | 0.19965 | 0.19965 | 0.0 | 0.17 Output | 0.00022554 | 0.00022554 | 0.00022554 | 0.0 | 0.00 Modify | 1.3114 | 1.3114 | 1.3114 | 0.0 | 1.08 Other | | 0.1585 | | | 0.13 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 368000.0 ave 368000 max 368000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 368000 Ave neighs/atom = 46.000000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 253.085859244995, Press = 0.88950811145657 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 18 18 18 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.684 | 4.684 | 4.684 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 14000 -30274.282 -30274.282 -30536.371 -30536.371 253.48208 253.48208 181983.31 181983.31 612.65246 612.65246 15000 -30277.849 -30277.849 -30533.863 -30533.863 247.60764 247.60764 181998.33 181998.33 513.97346 513.97346 Loop time of 120.828 on 1 procs for 1000 steps with 8000 atoms Performance: 0.715 ns/day, 33.563 hours/ns, 8.276 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 | 119.16 | 119.16 | 119.16 | 0.0 | 98.62 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.19991 | 0.19991 | 0.19991 | 0.0 | 0.17 Output | 0.00052542 | 0.00052542 | 0.00052542 | 0.0 | 0.00 Modify | 1.3074 | 1.3074 | 1.3074 | 0.0 | 1.08 Other | | 0.1583 | | | 0.13 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 367998.0 ave 367998 max 367998 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 367998 Ave neighs/atom = 45.999750 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 253.06580176201, Press = 1.31346677679157 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 18 18 18 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.684 | 4.684 | 4.684 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 15000 -30277.849 -30277.849 -30533.863 -30533.863 247.60764 247.60764 181998.33 181998.33 513.97346 513.97346 16000 -30272.103 -30272.103 -30533.935 -30533.935 253.23459 253.23459 181915.01 181915.01 877.78575 877.78575 Loop time of 128.494 on 1 procs for 1000 steps with 8000 atoms Performance: 0.672 ns/day, 35.693 hours/ns, 7.782 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 | 126.69 | 126.69 | 126.69 | 0.0 | 98.60 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.21665 | 0.21665 | 0.21665 | 0.0 | 0.17 Output | 0.00030725 | 0.00030725 | 0.00030725 | 0.0 | 0.00 Modify | 1.4148 | 1.4148 | 1.4148 | 0.0 | 1.10 Other | | 0.1698 | | | 0.13 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 367994.0 ave 367994 max 367994 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 367994 Ave neighs/atom = 45.999250 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 253.045287575866, Press = 1.30227932663944 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 18 18 18 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.684 | 4.684 | 4.684 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 16000 -30272.103 -30272.103 -30533.935 -30533.935 253.23459 253.23459 181915.01 181915.01 877.78575 877.78575 17000 -30278.66 -30278.66 -30542.603 -30542.603 255.27579 255.27579 182037.01 182037.01 394.9627 394.9627 Loop time of 143.442 on 1 procs for 1000 steps with 8000 atoms Performance: 0.602 ns/day, 39.845 hours/ns, 6.971 timesteps/s 99.7% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 141.42 | 141.42 | 141.42 | 0.0 | 98.59 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.23279 | 0.23279 | 0.23279 | 0.0 | 0.16 Output | 0.0002368 | 0.0002368 | 0.0002368 | 0.0 | 0.00 Modify | 1.6055 | 1.6055 | 1.6055 | 0.0 | 1.12 Other | | 0.1843 | | | 0.13 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 367996.0 ave 367996 max 367996 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 367996 Ave neighs/atom = 45.999500 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 253.036766291242, Press = 1.77345583380033 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 18 18 18 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.684 | 4.684 | 4.684 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 17000 -30278.66 -30278.66 -30542.603 -30542.603 255.27579 255.27579 182037.01 182037.01 394.9627 394.9627 18000 -30272.098 -30272.098 -30537.428 -30537.428 256.61743 256.61743 181988.01 181988.01 716.88744 716.88744 Loop time of 143.284 on 1 procs for 1000 steps with 8000 atoms Performance: 0.603 ns/day, 39.801 hours/ns, 6.979 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 | 141.22 | 141.22 | 141.22 | 0.0 | 98.56 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.23559 | 0.23559 | 0.23559 | 0.0 | 0.16 Output | 0.00029406 | 0.00029406 | 0.00029406 | 0.0 | 0.00 Modify | 1.6475 | 1.6475 | 1.6475 | 0.0 | 1.15 Other | | 0.1851 | | | 0.13 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 367998.0 ave 367998 max 367998 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 367998 Ave neighs/atom = 45.999750 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 252.955344168593, Press = 2.57846903081654 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 18 18 18 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.684 | 4.684 | 4.684 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 18000 -30272.098 -30272.098 -30537.428 -30537.428 256.61743 256.61743 181988.01 181988.01 716.88744 716.88744 19000 -30272.583 -30272.583 -30534.807 -30534.807 253.61293 253.61293 182063.46 182063.46 305.51221 305.51221 Loop time of 143.018 on 1 procs for 1000 steps with 8000 atoms Performance: 0.604 ns/day, 39.727 hours/ns, 6.992 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 | 140.95 | 140.95 | 140.95 | 0.0 | 98.55 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.23394 | 0.23394 | 0.23394 | 0.0 | 0.16 Output | 0.00023307 | 0.00023307 | 0.00023307 | 0.0 | 0.00 Modify | 1.6505 | 1.6505 | 1.6505 | 0.0 | 1.15 Other | | 0.1847 | | | 0.13 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 367998.0 ave 367998 max 367998 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 367998 Ave neighs/atom = 45.999750 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 252.985771501089, Press = 2.71392574352289 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 18 18 18 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.684 | 4.684 | 4.684 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 19000 -30272.583 -30272.583 -30534.807 -30534.807 253.61293 253.61293 182063.46 182063.46 305.51221 305.51221 20000 -30270.341 -30270.341 -30534.29 -30534.29 255.28207 255.28207 182063.19 182063.19 373.05417 373.05417 Loop time of 140.204 on 1 procs for 1000 steps with 8000 atoms Performance: 0.616 ns/day, 38.946 hours/ns, 7.132 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 | 138.19 | 138.19 | 138.19 | 0.0 | 98.57 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.22915 | 0.22915 | 0.22915 | 0.0 | 0.16 Output | 0.00031625 | 0.00031625 | 0.00031625 | 0.0 | 0.00 Modify | 1.6024 | 1.6024 | 1.6024 | 0.0 | 1.14 Other | | 0.1799 | | | 0.13 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 367992.0 ave 367992 max 367992 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 367992 Ave neighs/atom = 45.999000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 253.103681705364, Press = 2.11577760734045 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 18 18 18 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.684 | 4.684 | 4.684 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 20000 -30270.341 -30270.341 -30534.29 -30534.29 255.28207 255.28207 182063.19 182063.19 373.05417 373.05417 21000 -30279.584 -30279.584 -30540.828 -30540.828 252.66547 252.66547 182067.88 182067.88 292.19818 292.19818 Loop time of 142.488 on 1 procs for 1000 steps with 8000 atoms Performance: 0.606 ns/day, 39.580 hours/ns, 7.018 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 | 140.43 | 140.43 | 140.43 | 0.0 | 98.55 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.23343 | 0.23343 | 0.23343 | 0.0 | 0.16 Output | 0.00023202 | 0.00023202 | 0.00023202 | 0.0 | 0.00 Modify | 1.6455 | 1.6455 | 1.6455 | 0.0 | 1.15 Other | | 0.1833 | | | 0.13 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 367998.0 ave 367998 max 367998 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 367998 Ave neighs/atom = 45.999750 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 253.086508325137, Press = 2.16199963774457 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 18 18 18 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.684 | 4.684 | 4.684 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 21000 -30279.584 -30279.584 -30540.828 -30540.828 252.66547 252.66547 182067.88 182067.88 292.19818 292.19818 22000 -30272.895 -30272.895 -30534.965 -30534.965 253.46488 253.46488 182162.72 182162.72 -110.04571 -110.04571 Loop time of 142.622 on 1 procs for 1000 steps with 8000 atoms Performance: 0.606 ns/day, 39.617 hours/ns, 7.012 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 | 140.55 | 140.55 | 140.55 | 0.0 | 98.55 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.2329 | 0.2329 | 0.2329 | 0.0 | 0.16 Output | 0.00032206 | 0.00032206 | 0.00032206 | 0.0 | 0.00 Modify | 1.6526 | 1.6526 | 1.6526 | 0.0 | 1.16 Other | | 0.1837 | | | 0.13 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 368000.0 ave 368000 max 368000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 368000 Ave neighs/atom = 46.000000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 253.079494579412, Press = 1.74107000007484 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 18 18 18 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.684 | 4.684 | 4.684 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 22000 -30272.895 -30272.895 -30534.965 -30534.965 253.46488 253.46488 182162.72 182162.72 -110.04571 -110.04571 23000 -30277.365 -30277.365 -30534.967 -30534.967 249.14266 249.14266 182209.05 182209.05 -346.93847 -346.93847 Loop time of 130.894 on 1 procs for 1000 steps with 8000 atoms Performance: 0.660 ns/day, 36.359 hours/ns, 7.640 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 | 129.05 | 129.05 | 129.05 | 0.0 | 98.59 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.21449 | 0.21449 | 0.21449 | 0.0 | 0.16 Output | 0.00027669 | 0.00027669 | 0.00027669 | 0.0 | 0.00 Modify | 1.4579 | 1.4579 | 1.4579 | 0.0 | 1.11 Other | | 0.1689 | | | 0.13 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 368000.0 ave 368000 max 368000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 368000 Ave neighs/atom = 46.000000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 253.019397780179, Press = 1.72542016698381 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 18 18 18 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.684 | 4.684 | 4.684 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 23000 -30277.365 -30277.365 -30534.967 -30534.967 249.14266 249.14266 182209.05 182209.05 -346.93847 -346.93847 24000 -30271.835 -30271.835 -30539.194 -30539.194 258.58002 258.58002 182250.07 182250.07 -362.72116 -362.72116 Loop time of 120.749 on 1 procs for 1000 steps with 8000 atoms Performance: 0.716 ns/day, 33.541 hours/ns, 8.282 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 | 119.08 | 119.08 | 119.08 | 0.0 | 98.62 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.20094 | 0.20094 | 0.20094 | 0.0 | 0.17 Output | 0.00022767 | 0.00022767 | 0.00022767 | 0.0 | 0.00 Modify | 1.3102 | 1.3102 | 1.3102 | 0.0 | 1.09 Other | | 0.1583 | | | 0.13 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 367996.0 ave 367996 max 367996 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 367996 Ave neighs/atom = 45.999500 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 253.050017960325, Press = 1.83556008011189 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 18 18 18 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.684 | 4.684 | 4.684 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 24000 -30271.835 -30271.835 -30539.194 -30539.194 258.58002 258.58002 182250.07 182250.07 -362.72116 -362.72116 25000 -30279.086 -30279.086 -30535.415 -30535.415 247.91221 247.91221 182270.9 182270.9 -624.98143 -624.98143 Loop time of 120.769 on 1 procs for 1000 steps with 8000 atoms Performance: 0.715 ns/day, 33.547 hours/ns, 8.280 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 | 119.1 | 119.1 | 119.1 | 0.0 | 98.62 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.19934 | 0.19934 | 0.19934 | 0.0 | 0.17 Output | 0.00022817 | 0.00022817 | 0.00022817 | 0.0 | 0.00 Modify | 1.3086 | 1.3086 | 1.3086 | 0.0 | 1.08 Other | | 0.1587 | | | 0.13 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 368000.0 ave 368000 max 368000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 368000 Ave neighs/atom = 46.000000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 253.008125317691, Press = 0.973748009966577 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 18 18 18 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.684 | 4.684 | 4.684 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 25000 -30279.086 -30279.086 -30535.415 -30535.415 247.91221 247.91221 182270.9 182270.9 -624.98143 -624.98143 26000 -30275.585 -30275.585 -30540.849 -30540.849 256.55268 256.55268 182263.12 182263.12 -450.38255 -450.38255 Loop time of 128.306 on 1 procs for 1000 steps with 8000 atoms Performance: 0.673 ns/day, 35.641 hours/ns, 7.794 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 | 126.51 | 126.51 | 126.51 | 0.0 | 98.60 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.21009 | 0.21009 | 0.21009 | 0.0 | 0.16 Output | 0.00023778 | 0.00023778 | 0.00023778 | 0.0 | 0.00 Modify | 1.4209 | 1.4209 | 1.4209 | 0.0 | 1.11 Other | | 0.1668 | | | 0.13 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 367998.0 ave 367998 max 367998 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 367998 Ave neighs/atom = 45.999750 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 252.975307868795, Press = 0.557954127957422 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 18 18 18 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.684 | 4.684 | 4.684 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 26000 -30275.585 -30275.585 -30540.849 -30540.849 256.55268 256.55268 182263.12 182263.12 -450.38255 -450.38255 27000 -30271.615 -30271.615 -30532.682 -30532.682 252.49458 252.49458 182287.34 182287.34 -624.31855 -624.31855 Loop time of 143.379 on 1 procs for 1000 steps with 8000 atoms Performance: 0.603 ns/day, 39.827 hours/ns, 6.975 timesteps/s 99.7% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 141.33 | 141.33 | 141.33 | 0.0 | 98.57 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.23391 | 0.23391 | 0.23391 | 0.0 | 0.16 Output | 0.00032404 | 0.00032404 | 0.00032404 | 0.0 | 0.00 Modify | 1.6278 | 1.6278 | 1.6278 | 0.0 | 1.14 Other | | 0.1825 | | | 0.13 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 367996.0 ave 367996 max 367996 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 367996 Ave neighs/atom = 45.999500 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 252.950796106052, Press = 0.659920439429283 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 18 18 18 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.684 | 4.684 | 4.684 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 27000 -30271.615 -30271.615 -30532.682 -30532.682 252.49458 252.49458 182287.34 182287.34 -624.31855 -624.31855 28000 -30271.904 -30271.904 -30537.115 -30537.115 256.50278 256.50278 182209.68 182209.68 -267.3497 -267.3497 Loop time of 143.048 on 1 procs for 1000 steps with 8000 atoms Performance: 0.604 ns/day, 39.735 hours/ns, 6.991 timesteps/s 99.6% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 140.96 | 140.96 | 140.96 | 0.0 | 98.54 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.2389 | 0.2389 | 0.2389 | 0.0 | 0.17 Output | 0.00031446 | 0.00031446 | 0.00031446 | 0.0 | 0.00 Modify | 1.6585 | 1.6585 | 1.6585 | 0.0 | 1.16 Other | | 0.186 | | | 0.13 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 368000.0 ave 368000 max 368000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 368000 Ave neighs/atom = 46.000000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 252.964358995194, Press = 0.744966462670184 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 18 18 18 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.684 | 4.684 | 4.684 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 28000 -30271.904 -30271.904 -30537.115 -30537.115 256.50278 256.50278 182209.68 182209.68 -267.3497 -267.3497 29000 -30276.401 -30276.401 -30534.736 -30534.736 249.85192 249.85192 182262.36 182262.36 -507.6469 -507.6469 Loop time of 143.376 on 1 procs for 1000 steps with 8000 atoms Performance: 0.603 ns/day, 39.827 hours/ns, 6.975 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 | 141.3 | 141.3 | 141.3 | 0.0 | 98.55 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.23724 | 0.23724 | 0.23724 | 0.0 | 0.17 Output | 0.00031098 | 0.00031098 | 0.00031098 | 0.0 | 0.00 Modify | 1.6512 | 1.6512 | 1.6512 | 0.0 | 1.15 Other | | 0.1847 | | | 0.13 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 368000.0 ave 368000 max 368000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 368000 Ave neighs/atom = 46.000000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 252.975970648373, Press = 0.33543536057986 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 18 18 18 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.684 | 4.684 | 4.684 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 29000 -30276.401 -30276.401 -30534.736 -30534.736 249.85192 249.85192 182262.36 182262.36 -507.6469 -507.6469 30000 -30275.523 -30275.523 -30534.273 -30534.273 250.25293 250.25293 182183.62 182183.62 -130.38988 -130.38988 Loop time of 125.12 on 1 procs for 1000 steps with 8000 atoms Performance: 0.691 ns/day, 34.756 hours/ns, 7.992 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 | 123.36 | 123.36 | 123.36 | 0.0 | 98.59 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.21796 | 0.21796 | 0.21796 | 0.0 | 0.17 Output | 0.00023278 | 0.00023278 | 0.00023278 | 0.0 | 0.00 Modify | 1.3823 | 1.3823 | 1.3823 | 0.0 | 1.10 Other | | 0.1624 | | | 0.13 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 367998.0 ave 367998 max 367998 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 367998 Ave neighs/atom = 45.999750 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 252.970238226362, Press = 0.404178736107357 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 18 18 18 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.684 | 4.684 | 4.684 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 30000 -30275.523 -30275.523 -30534.273 -30534.273 250.25293 250.25293 182183.62 182183.62 -130.38988 -130.38988 31000 -30272.288 -30272.288 -30532.667 -30532.667 251.82868 251.82868 182187.33 182187.33 -197.366 -197.366 Loop time of 121.958 on 1 procs for 1000 steps with 8000 atoms Performance: 0.708 ns/day, 33.877 hours/ns, 8.200 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 | 120.25 | 120.25 | 120.25 | 0.0 | 98.60 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.21351 | 0.21351 | 0.21351 | 0.0 | 0.18 Output | 0.00023095 | 0.00023095 | 0.00023095 | 0.0 | 0.00 Modify | 1.3345 | 1.3345 | 1.3345 | 0.0 | 1.09 Other | | 0.1592 | | | 0.13 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725.00 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 367998.0 ave 367998 max 367998 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 367998 Ave neighs/atom = 45.999750 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" jump SELF break # Write final averaged volume to file if temperature and volume have converged; otherwise wirte a # flag to indicate non-convergence. variable myStep equal step if "${myStep} < 2000000" then "print '${V}' file output/vol_T253.15.out" else "print 'not_converged' file output/vol_T253.15.out" print '${V}' file output/vol_T253.15.out 182138.391804633 print "LAMMPS calculation completed" LAMMPS calculation completed quit 0