# 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.1465589553117748*${_u_distance} variable latticeconst_converted equal 3.1465589553117748*1 lattice bcc ${latticeconst_converted} lattice bcc 3.14655895531177 Lattice spacing in x,y,z = 3.1465590 3.1465590 3.1465590 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 (31.465590 31.465590 31.465590) 1 by 1 by 1 MPI processor grid create_atoms 1 box Created 2000 atoms using lattice units in orthogonal box = (0.0000000 0.0000000 0.0000000) to (31.465590 31.465590 31.465590) create_atoms CPU = 0.004 seconds variable mass_converted equal 95.94*${_u_mass} variable mass_converted equal 95.94*1 kim_interactions Mo WARNING: 'kim_' has been renamed to 'kim '. Please update your input. kim interactions Mo #=== BEGIN kim interactions ================================== pair_style kim MEAM_LAMMPS_LeeBaskesKim_2001_Mo__MO_805823015127_000 pair_coeff * * Mo #=== END kim interactions ==================================== mass 1 ${mass_converted} mass 1 95.94 # initial volume variable v equal vol # assign formula variable V0 equal ${v} # evaluate initial value variable V0 equal 31153.5555569502 variable V0_metal equal ${V0}/(${_u_distance}*${_u_distance}*${_u_distance}) variable V0_metal equal 31153.5555569502/(${_u_distance}*${_u_distance}*${_u_distance}) variable V0_metal equal 31153.5555569502/(1*${_u_distance}*${_u_distance}) variable V0_metal equal 31153.5555569502/(1*1*${_u_distance}) variable V0_metal equal 31153.5555569502/(1*1*1) variable V0_metal_times1000 equal ${V0_metal}*1000 variable V0_metal_times1000 equal 31153.5555569502*1000 print "Initial system volume: ${V0_metal} Angstroms^3" Initial system volume: 31153.5555569502 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 313.15*${_u_temperature} variable temp_converted equal 313.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 313.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 313.15 ${temp_converted} ${Tdamp_converted} iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 313.15 313.15 ${Tdamp_converted} iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 313.15 313.15 0.1 iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 313.15 313.15 0.1 iso 0 ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 313.15 313.15 0.1 iso 0 0 ${Pdamp_converted} fix ensemble all npt temp 313.15 313.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 "313.15 - 0.2" variable T_up equal "313.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_805823015127_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 = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 0 -13539.085 -13539.085 -13620 -13620 313.15 313.15 31153.556 31153.556 2774.26 2774.26 1000 -13454.118 -13454.118 -13531.665 -13531.665 300.11712 300.11712 31324.333 31324.333 -1328.8309 -1328.8309 Loop time of 97.4138 on 1 procs for 1000 steps with 2000 atoms Performance: 0.887 ns/day, 27.059 hours/ns, 10.265 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 | 96.917 | 96.917 | 96.917 | 0.0 | 99.49 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.077754 | 0.077754 | 0.077754 | 0.0 | 0.08 Output | 0.00020511 | 0.00020511 | 0.00020511 | 0.0 | 0.00 Modify | 0.36088 | 0.36088 | 0.36088 | 0.0 | 0.37 Other | | 0.05752 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 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: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.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 = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 1000 -13454.118 -13454.118 -13531.665 -13531.665 300.11712 300.11712 31324.333 31324.333 -1328.8309 -1328.8309 2000 -13454.652 -13454.652 -13536.457 -13536.457 316.59105 316.59105 31271.775 31271.775 2460.6475 2460.6475 Loop time of 100.307 on 1 procs for 1000 steps with 2000 atoms Performance: 0.861 ns/day, 27.863 hours/ns, 9.969 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 | 99.82 | 99.82 | 99.82 | 0.0 | 99.51 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.074911 | 0.074911 | 0.074911 | 0.0 | 0.07 Output | 0.00022398 | 0.00022398 | 0.00022398 | 0.0 | 0.00 Modify | 0.35427 | 0.35427 | 0.35427 | 0.0 | 0.35 Other | | 0.05723 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 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: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.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 = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 2000 -13454.652 -13454.652 -13536.457 -13536.457 316.59105 316.59105 31271.775 31271.775 2460.6475 2460.6475 3000 -13455.779 -13455.779 -13535.695 -13535.695 309.28408 309.28408 31312.658 31312.658 -843.3457 -843.3457 Loop time of 100.445 on 1 procs for 1000 steps with 2000 atoms Performance: 0.860 ns/day, 27.901 hours/ns, 9.956 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 | 99.958 | 99.958 | 99.958 | 0.0 | 99.51 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.073858 | 0.073858 | 0.073858 | 0.0 | 0.07 Output | 0.00018945 | 0.00018945 | 0.00018945 | 0.0 | 0.00 Modify | 0.35542 | 0.35542 | 0.35542 | 0.0 | 0.35 Other | | 0.05778 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 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: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.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 = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 3000 -13455.779 -13455.779 -13535.695 -13535.695 309.28408 309.28408 31312.658 31312.658 -843.3457 -843.3457 4000 -13454.436 -13454.436 -13526.643 -13526.643 279.44735 279.44735 31313.387 31313.387 -124.83501 -124.83501 Loop time of 98.9743 on 1 procs for 1000 steps with 2000 atoms Performance: 0.873 ns/day, 27.493 hours/ns, 10.104 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 | 98.498 | 98.498 | 98.498 | 0.0 | 99.52 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.073419 | 0.073419 | 0.073419 | 0.0 | 0.07 Output | 0.00022311 | 0.00022311 | 0.00022311 | 0.0 | 0.00 Modify | 0.34616 | 0.34616 | 0.34616 | 0.0 | 0.35 Other | | 0.05695 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 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: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.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 = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 4000 -13454.436 -13454.436 -13526.643 -13526.643 279.44735 279.44735 31313.387 31313.387 -124.83501 -124.83501 5000 -13456.279 -13456.279 -13538.004 -13538.004 316.28131 316.28131 31271.305 31271.305 2393.7871 2393.7871 Loop time of 101.506 on 1 procs for 1000 steps with 2000 atoms Performance: 0.851 ns/day, 28.196 hours/ns, 9.852 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 | 101.01 | 101.01 | 101.01 | 0.0 | 99.51 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.075516 | 0.075516 | 0.075516 | 0.0 | 0.07 Output | 0.00018993 | 0.00018993 | 0.00018993 | 0.0 | 0.00 Modify | 0.36096 | 0.36096 | 0.36096 | 0.0 | 0.36 Other | | 0.05772 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 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: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.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 = 314.227965150513, Press = -543.864168564405 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 = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 5000 -13456.279 -13456.279 -13538.004 -13538.004 316.28131 316.28131 31271.305 31271.305 2393.7871 2393.7871 6000 -13454.544 -13454.544 -13538.158 -13538.158 323.59437 323.59437 31289.674 31289.674 767.64799 767.64799 Loop time of 100.564 on 1 procs for 1000 steps with 2000 atoms Performance: 0.859 ns/day, 27.934 hours/ns, 9.944 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 | 100.06 | 100.06 | 100.06 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.074671 | 0.074671 | 0.074671 | 0.0 | 0.07 Output | 0.00016453 | 0.00016453 | 0.00016453 | 0.0 | 0.00 Modify | 0.37128 | 0.37128 | 0.37128 | 0.0 | 0.37 Other | | 0.05726 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 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: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.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 = 312.287746652785, Press = -15.4216175065583 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 = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 6000 -13454.544 -13454.544 -13538.158 -13538.158 323.59437 323.59437 31289.674 31289.674 767.64799 767.64799 7000 -13452.235 -13452.235 -13534.786 -13534.786 319.47987 319.47987 31343.416 31343.416 -3110.6134 -3110.6134 Loop time of 102.598 on 1 procs for 1000 steps with 2000 atoms Performance: 0.842 ns/day, 28.499 hours/ns, 9.747 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 | 102.08 | 102.08 | 102.08 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.075701 | 0.075701 | 0.075701 | 0.0 | 0.07 Output | 0.00020146 | 0.00020146 | 0.00020146 | 0.0 | 0.00 Modify | 0.38349 | 0.38349 | 0.38349 | 0.0 | 0.37 Other | | 0.05816 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 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: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.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 = 313.060748972671, Press = -20.3965768730545 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 = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 7000 -13452.235 -13452.235 -13534.786 -13534.786 319.47987 319.47987 31343.416 31343.416 -3110.6134 -3110.6134 8000 -13453.521 -13453.521 -13535.739 -13535.739 318.1894 318.1894 31285.475 31285.475 1397.5635 1397.5635 Loop time of 100.683 on 1 procs for 1000 steps with 2000 atoms Performance: 0.858 ns/day, 27.968 hours/ns, 9.932 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 | 100.18 | 100.18 | 100.18 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.0746 | 0.0746 | 0.0746 | 0.0 | 0.07 Output | 0.00015049 | 0.00015049 | 0.00015049 | 0.0 | 0.00 Modify | 0.37262 | 0.37262 | 0.37262 | 0.0 | 0.37 Other | | 0.05769 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 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: 115998.0 ave 115998 max 115998 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 115998 Ave neighs/atom = 57.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 = 312.635748203997, Press = -0.691926476438587 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 = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 8000 -13453.521 -13453.521 -13535.739 -13535.739 318.1894 318.1894 31285.475 31285.475 1397.5635 1397.5635 9000 -13454.863 -13454.863 -13534.967 -13534.967 310.01101 310.01101 31328.753 31328.753 -2032.0606 -2032.0606 Loop time of 103.774 on 1 procs for 1000 steps with 2000 atoms Performance: 0.833 ns/day, 28.826 hours/ns, 9.636 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 | 103.25 | 103.25 | 103.25 | 0.0 | 99.49 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.076362 | 0.076362 | 0.076362 | 0.0 | 0.07 Output | 0.00015624 | 0.00015624 | 0.00015624 | 0.0 | 0.00 Modify | 0.39142 | 0.39142 | 0.39142 | 0.0 | 0.38 Other | | 0.05844 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 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: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.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 = 313.155266242092, Press = -6.33389363357939 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 = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 9000 -13454.863 -13454.863 -13534.967 -13534.967 310.01101 310.01101 31328.753 31328.753 -2032.0606 -2032.0606 10000 -13454.445 -13454.445 -13533.503 -13533.503 305.96346 305.96346 31263.552 31263.552 3255.2636 3255.2636 Loop time of 101.759 on 1 procs for 1000 steps with 2000 atoms Performance: 0.849 ns/day, 28.266 hours/ns, 9.827 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 | 101.25 | 101.25 | 101.25 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.075128 | 0.075128 | 0.075128 | 0.0 | 0.07 Output | 0.0002145 | 0.0002145 | 0.0002145 | 0.0 | 0.00 Modify | 0.37823 | 0.37823 | 0.37823 | 0.0 | 0.37 Other | | 0.05779 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 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: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.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 = 312.967108139006, Press = -4.2392618164284 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 = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 10000 -13454.445 -13454.445 -13533.503 -13533.503 305.96346 305.96346 31263.552 31263.552 3255.2636 3255.2636 11000 -13459.058 -13459.058 -13536.66 -13536.66 300.3267 300.3267 31293.501 31293.501 415.28233 415.28233 Loop time of 100.598 on 1 procs for 1000 steps with 2000 atoms Performance: 0.859 ns/day, 27.944 hours/ns, 9.941 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 | 100.09 | 100.09 | 100.09 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.074539 | 0.074539 | 0.074539 | 0.0 | 0.07 Output | 0.00015619 | 0.00015619 | 0.00015619 | 0.0 | 0.00 Modify | 0.37184 | 0.37184 | 0.37184 | 0.0 | 0.37 Other | | 0.05742 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 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: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.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 = 312.703453124194, Press = 2.74041787666002 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 = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 11000 -13459.058 -13459.058 -13536.66 -13536.66 300.3267 300.3267 31293.501 31293.501 415.28233 415.28233 12000 -13454.635 -13454.635 -13534.157 -13534.157 307.75716 307.75716 31313.967 31313.967 -906.5387 -906.5387 Loop time of 101.363 on 1 procs for 1000 steps with 2000 atoms Performance: 0.852 ns/day, 28.156 hours/ns, 9.866 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 | 100.86 | 100.86 | 100.86 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.074846 | 0.074846 | 0.074846 | 0.0 | 0.07 Output | 0.00015962 | 0.00015962 | 0.00015962 | 0.0 | 0.00 Modify | 0.37551 | 0.37551 | 0.37551 | 0.0 | 0.37 Other | | 0.05697 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 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: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.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 = 312.389184182097, Press = -3.86764195895134 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 = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 12000 -13454.635 -13454.635 -13534.157 -13534.157 307.75716 307.75716 31313.967 31313.967 -906.5387 -906.5387 13000 -13458.02 -13458.02 -13533.762 -13533.762 293.13223 293.13223 31279.537 31279.537 1920.4756 1920.4756 Loop time of 101.04 on 1 procs for 1000 steps with 2000 atoms Performance: 0.855 ns/day, 28.067 hours/ns, 9.897 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 | 100.53 | 100.53 | 100.53 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.07514 | 0.07514 | 0.07514 | 0.0 | 0.07 Output | 0.00015871 | 0.00015871 | 0.00015871 | 0.0 | 0.00 Modify | 0.37582 | 0.37582 | 0.37582 | 0.0 | 0.37 Other | | 0.05735 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 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: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.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 = 312.303202637024, Press = -2.71421337758439 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 = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 13000 -13458.02 -13458.02 -13533.762 -13533.762 293.13223 293.13223 31279.537 31279.537 1920.4756 1920.4756 14000 -13453.486 -13453.486 -13537.315 -13537.315 324.42752 324.42752 31340.029 31340.029 -3194.5698 -3194.5698 Loop time of 101.842 on 1 procs for 1000 steps with 2000 atoms Performance: 0.848 ns/day, 28.290 hours/ns, 9.819 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 | 101.32 | 101.32 | 101.32 | 0.0 | 99.49 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.076855 | 0.076855 | 0.076855 | 0.0 | 0.08 Output | 0.00020034 | 0.00020034 | 0.00020034 | 0.0 | 0.00 Modify | 0.38309 | 0.38309 | 0.38309 | 0.0 | 0.38 Other | | 0.05816 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 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: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.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 = 312.299108739598, Press = -2.80119070608315 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 = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 14000 -13453.486 -13453.486 -13537.315 -13537.315 324.42752 324.42752 31340.029 31340.029 -3194.5698 -3194.5698 15000 -13458.682 -13458.682 -13543.57 -13543.57 328.5231 328.5231 31270.067 31270.067 1990.8261 1990.8261 Loop time of 101.196 on 1 procs for 1000 steps with 2000 atoms Performance: 0.854 ns/day, 28.110 hours/ns, 9.882 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 | 100.69 | 100.69 | 100.69 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.07517 | 0.07517 | 0.07517 | 0.0 | 0.07 Output | 0.00037221 | 0.00037221 | 0.00037221 | 0.0 | 0.00 Modify | 0.37538 | 0.37538 | 0.37538 | 0.0 | 0.37 Other | | 0.05805 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 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: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.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 = 312.075225544419, Press = -3.8026342366017 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 = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 15000 -13458.682 -13458.682 -13543.57 -13543.57 328.5231 328.5231 31270.067 31270.067 1990.8261 1990.8261 16000 -13454.159 -13454.159 -13537.288 -13537.288 321.72012 321.72012 31322.665 31322.665 -1695.3302 -1695.3302 Loop time of 101.278 on 1 procs for 1000 steps with 2000 atoms Performance: 0.853 ns/day, 28.133 hours/ns, 9.874 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 | 100.77 | 100.77 | 100.77 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.075119 | 0.075119 | 0.075119 | 0.0 | 0.07 Output | 0.00016136 | 0.00016136 | 0.00016136 | 0.0 | 0.00 Modify | 0.37791 | 0.37791 | 0.37791 | 0.0 | 0.37 Other | | 0.05773 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 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: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.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 = 312.001738772468, Press = 4.30977673265634 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 = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 16000 -13454.159 -13454.159 -13537.288 -13537.288 321.72012 321.72012 31322.665 31322.665 -1695.3302 -1695.3302 17000 -13458.827 -13458.827 -13539.323 -13539.323 311.52509 311.52509 31290.268 31290.268 649.43848 649.43848 Loop time of 101.567 on 1 procs for 1000 steps with 2000 atoms Performance: 0.851 ns/day, 28.213 hours/ns, 9.846 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 | 101.06 | 101.06 | 101.06 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.0757 | 0.0757 | 0.0757 | 0.0 | 0.07 Output | 0.00015892 | 0.00015892 | 0.00015892 | 0.0 | 0.00 Modify | 0.3761 | 0.3761 | 0.3761 | 0.0 | 0.37 Other | | 0.05765 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 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: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.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 = 311.908586296153, Press = -3.94873498581654 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 = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 17000 -13458.827 -13458.827 -13539.323 -13539.323 311.52509 311.52509 31290.268 31290.268 649.43848 649.43848 18000 -13454.639 -13454.639 -13535.927 -13535.927 314.5922 314.5922 31312.319 31312.319 -924.24875 -924.24875 Loop time of 99.6761 on 1 procs for 1000 steps with 2000 atoms Performance: 0.867 ns/day, 27.688 hours/ns, 10.032 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 | 99.177 | 99.177 | 99.177 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.07398 | 0.07398 | 0.07398 | 0.0 | 0.07 Output | 0.00015682 | 0.00015682 | 0.00015682 | 0.0 | 0.00 Modify | 0.36798 | 0.36798 | 0.36798 | 0.0 | 0.37 Other | | 0.05715 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 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: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.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 = 311.94439516112, Press = -0.474912236656407 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 = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 18000 -13454.639 -13454.639 -13535.927 -13535.927 314.5922 314.5922 31312.319 31312.319 -924.24875 -924.24875 19000 -13455.937 -13455.937 -13535.295 -13535.295 307.12386 307.12386 31289.407 31289.407 1162.7551 1162.7551 Loop time of 99.3686 on 1 procs for 1000 steps with 2000 atoms Performance: 0.869 ns/day, 27.602 hours/ns, 10.064 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 | 98.872 | 98.872 | 98.872 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.074232 | 0.074232 | 0.074232 | 0.0 | 0.07 Output | 0.00015893 | 0.00015893 | 0.00015893 | 0.0 | 0.00 Modify | 0.36571 | 0.36571 | 0.36571 | 0.0 | 0.37 Other | | 0.05642 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 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: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.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 = 312.271969431171, Press = -2.05167206469291 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 = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 19000 -13455.937 -13455.937 -13535.295 -13535.295 307.12386 307.12386 31289.407 31289.407 1162.7551 1162.7551 20000 -13453.241 -13453.241 -13536.225 -13536.225 321.15865 321.15865 31296.28 31296.28 551.59716 551.59716 Loop time of 100.449 on 1 procs for 1000 steps with 2000 atoms Performance: 0.860 ns/day, 27.903 hours/ns, 9.955 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 | 99.945 | 99.945 | 99.945 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.074652 | 0.074652 | 0.074652 | 0.0 | 0.07 Output | 0.00018391 | 0.00018391 | 0.00018391 | 0.0 | 0.00 Modify | 0.37235 | 0.37235 | 0.37235 | 0.0 | 0.37 Other | | 0.05718 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 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: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.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 = 312.367106244393, Press = -1.54273560773782 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 = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 20000 -13453.241 -13453.241 -13536.225 -13536.225 321.15865 321.15865 31296.28 31296.28 551.59716 551.59716 21000 -13453.345 -13453.345 -13533.698 -13533.698 310.97437 310.97437 31299.164 31299.164 393.17568 393.17568 Loop time of 101.055 on 1 procs for 1000 steps with 2000 atoms Performance: 0.855 ns/day, 28.071 hours/ns, 9.896 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 | 100.55 | 100.55 | 100.55 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.075175 | 0.075175 | 0.075175 | 0.0 | 0.07 Output | 0.00015907 | 0.00015907 | 0.00015907 | 0.0 | 0.00 Modify | 0.37567 | 0.37567 | 0.37567 | 0.0 | 0.37 Other | | 0.05729 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 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: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.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 = 312.53168000869, Press = -0.79367489249414 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 = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 21000 -13453.345 -13453.345 -13533.698 -13533.698 310.97437 310.97437 31299.164 31299.164 393.17568 393.17568 22000 -13452.912 -13452.912 -13535.089 -13535.089 318.03412 318.03412 31321.326 31321.326 -1393.7447 -1393.7447 Loop time of 100.567 on 1 procs for 1000 steps with 2000 atoms Performance: 0.859 ns/day, 27.935 hours/ns, 9.944 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 | 100.06 | 100.06 | 100.06 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.075041 | 0.075041 | 0.075041 | 0.0 | 0.07 Output | 0.00015255 | 0.00015255 | 0.00015255 | 0.0 | 0.00 Modify | 0.37399 | 0.37399 | 0.37399 | 0.0 | 0.37 Other | | 0.05704 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 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: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.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 = 312.517246061162, Press = 0.677762645452692 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 = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 22000 -13452.912 -13452.912 -13535.089 -13535.089 318.03412 318.03412 31321.326 31321.326 -1393.7447 -1393.7447 23000 -13454.53 -13454.53 -13536.681 -13536.681 317.93469 317.93469 31293.684 31293.684 687.10598 687.10598 Loop time of 101.117 on 1 procs for 1000 steps with 2000 atoms Performance: 0.854 ns/day, 28.088 hours/ns, 9.889 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 | 100.61 | 100.61 | 100.61 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.075034 | 0.075034 | 0.075034 | 0.0 | 0.07 Output | 0.00017585 | 0.00017585 | 0.00017585 | 0.0 | 0.00 Modify | 0.37654 | 0.37654 | 0.37654 | 0.0 | 0.37 Other | | 0.05796 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 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: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.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 = 312.730618191143, Press = -3.28551318304626 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 = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 23000 -13454.53 -13454.53 -13536.681 -13536.681 317.93469 317.93469 31293.684 31293.684 687.10598 687.10598 24000 -13452.014 -13452.014 -13538.646 -13538.646 335.27304 335.27304 31308.572 31308.572 -660.18872 -660.18872 Loop time of 100.08 on 1 procs for 1000 steps with 2000 atoms Performance: 0.863 ns/day, 27.800 hours/ns, 9.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 | 99.575 | 99.575 | 99.575 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.074824 | 0.074824 | 0.074824 | 0.0 | 0.07 Output | 0.00025963 | 0.00025963 | 0.00025963 | 0.0 | 0.00 Modify | 0.37257 | 0.37257 | 0.37257 | 0.0 | 0.37 Other | | 0.05709 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 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: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.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 = 312.780693457656, Press = 0.261491025578771 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 = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 24000 -13452.014 -13452.014 -13538.646 -13538.646 335.27304 335.27304 31308.572 31308.572 -660.18872 -660.18872 25000 -13456.835 -13456.835 -13535.424 -13535.424 304.14728 304.14728 31329.253 31329.253 -2276.8037 -2276.8037 Loop time of 99.4299 on 1 procs for 1000 steps with 2000 atoms Performance: 0.869 ns/day, 27.619 hours/ns, 10.057 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 | 98.931 | 98.931 | 98.931 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.073913 | 0.073913 | 0.073913 | 0.0 | 0.07 Output | 0.00015939 | 0.00015939 | 0.00015939 | 0.0 | 0.00 Modify | 0.36775 | 0.36775 | 0.36775 | 0.0 | 0.37 Other | | 0.05676 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 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: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.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 = 312.911802858601, Press = -3.68936217252647 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 = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 25000 -13456.835 -13456.835 -13535.424 -13535.424 304.14728 304.14728 31329.253 31329.253 -2276.8037 -2276.8037 26000 -13452.764 -13452.764 -13534.537 -13534.537 316.46961 316.46961 31250.951 31250.951 4368.5678 4368.5678 Loop time of 99.5662 on 1 procs for 1000 steps with 2000 atoms Performance: 0.868 ns/day, 27.657 hours/ns, 10.044 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 | 99.069 | 99.069 | 99.069 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.073495 | 0.073495 | 0.073495 | 0.0 | 0.07 Output | 0.00015697 | 0.00015697 | 0.00015697 | 0.0 | 0.00 Modify | 0.36645 | 0.36645 | 0.36645 | 0.0 | 0.37 Other | | 0.05714 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 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: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.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 = 312.953477048826, Press = 1.72720896233588 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 = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 26000 -13452.764 -13452.764 -13534.537 -13534.537 316.46961 316.46961 31250.951 31250.951 4368.5678 4368.5678 27000 -13456.409 -13456.409 -13533.145 -13533.145 296.97811 296.97811 31320.23 31320.23 -1265.171 -1265.171 Loop time of 99.1875 on 1 procs for 1000 steps with 2000 atoms Performance: 0.871 ns/day, 27.552 hours/ns, 10.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 | 98.69 | 98.69 | 98.69 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.073697 | 0.073697 | 0.073697 | 0.0 | 0.07 Output | 0.00015677 | 0.00015677 | 0.00015677 | 0.0 | 0.00 Modify | 0.36677 | 0.36677 | 0.36677 | 0.0 | 0.37 Other | | 0.05703 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 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: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.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 = 313.121986383466, Press = -1.25624110713638 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 = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 27000 -13456.409 -13456.409 -13533.145 -13533.145 296.97811 296.97811 31320.23 31320.23 -1265.171 -1265.171 28000 -13453.773 -13453.773 -13533.185 -13533.185 307.33368 307.33368 31314.647 31314.647 -856.72653 -856.72653 Loop time of 99.6301 on 1 procs for 1000 steps with 2000 atoms Performance: 0.867 ns/day, 27.675 hours/ns, 10.037 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 | 99.131 | 99.131 | 99.131 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.074067 | 0.074067 | 0.074067 | 0.0 | 0.07 Output | 0.00015856 | 0.00015856 | 0.00015856 | 0.0 | 0.00 Modify | 0.3677 | 0.3677 | 0.3677 | 0.0 | 0.37 Other | | 0.0572 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 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: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.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 = 313.121788986267, Press = -1.32069593284881 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 = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 28000 -13453.773 -13453.773 -13533.185 -13533.185 307.33368 307.33368 31314.647 31314.647 -856.72653 -856.72653 29000 -13457.029 -13457.029 -13540.898 -13540.898 324.57934 324.57934 31299.98 31299.98 -193.83299 -193.83299 Loop time of 99.6253 on 1 procs for 1000 steps with 2000 atoms Performance: 0.867 ns/day, 27.674 hours/ns, 10.038 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 | 99.126 | 99.126 | 99.126 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.074179 | 0.074179 | 0.074179 | 0.0 | 0.07 Output | 0.00015681 | 0.00015681 | 0.00015681 | 0.0 | 0.00 Modify | 0.3681 | 0.3681 | 0.3681 | 0.0 | 0.37 Other | | 0.0571 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 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: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.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 = 313.120567232253, Press = 0.530554439460748 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 = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 29000 -13457.029 -13457.029 -13540.898 -13540.898 324.57934 324.57934 31299.98 31299.98 -193.83299 -193.83299 30000 -13452.953 -13452.953 -13534.007 -13534.007 313.68792 313.68792 31311.111 31311.111 -443.97141 -443.97141 Loop time of 100.129 on 1 procs for 1000 steps with 2000 atoms Performance: 0.863 ns/day, 27.814 hours/ns, 9.987 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 | 99.631 | 99.631 | 99.631 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.073991 | 0.073991 | 0.073991 | 0.0 | 0.07 Output | 0.00015465 | 0.00015465 | 0.00015465 | 0.0 | 0.00 Modify | 0.36731 | 0.36731 | 0.36731 | 0.0 | 0.37 Other | | 0.05738 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 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: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.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 = 313.103678714925, Press = -2.98141258226742 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 = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 30000 -13452.953 -13452.953 -13534.007 -13534.007 313.68792 313.68792 31311.111 31311.111 -443.97141 -443.97141 31000 -13456.484 -13456.484 -13537.924 -13537.924 315.18122 315.18122 31268.127 31268.127 2609.3285 2609.3285 Loop time of 101.037 on 1 procs for 1000 steps with 2000 atoms Performance: 0.855 ns/day, 28.066 hours/ns, 9.897 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 | 100.53 | 100.53 | 100.53 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.07481 | 0.07481 | 0.07481 | 0.0 | 0.07 Output | 0.00015603 | 0.00015603 | 0.00015603 | 0.0 | 0.00 Modify | 0.37773 | 0.37773 | 0.37773 | 0.0 | 0.37 Other | | 0.05742 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 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: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.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 = 313.1807460548, Press = 1.61202570031983 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 = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 31000 -13456.484 -13456.484 -13537.924 -13537.924 315.18122 315.18122 31268.127 31268.127 2609.3285 2609.3285 32000 -13453.874 -13453.874 -13534.236 -13534.236 311.00948 311.00948 31330.701 31330.701 -2235.9428 -2235.9428 Loop time of 99.4561 on 1 procs for 1000 steps with 2000 atoms Performance: 0.869 ns/day, 27.627 hours/ns, 10.055 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 | 98.957 | 98.957 | 98.957 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.074627 | 0.074627 | 0.074627 | 0.0 | 0.08 Output | 0.000157 | 0.000157 | 0.000157 | 0.0 | 0.00 Modify | 0.36745 | 0.36745 | 0.36745 | 0.0 | 0.37 Other | | 0.0574 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 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: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.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 = 313.230558671457, Press = -1.69188052400936 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 = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 32000 -13453.874 -13453.874 -13534.236 -13534.236 311.00948 311.00948 31330.701 31330.701 -2235.9428 -2235.9428 33000 -13455.842 -13455.842 -13535.65 -13535.65 308.86628 308.86628 31278.976 31278.976 1986.4914 1986.4914 Loop time of 99.8079 on 1 procs for 1000 steps with 2000 atoms Performance: 0.866 ns/day, 27.724 hours/ns, 10.019 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 | 99.307 | 99.307 | 99.307 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.073983 | 0.073983 | 0.073983 | 0.0 | 0.07 Output | 0.00015971 | 0.00015971 | 0.00015971 | 0.0 | 0.00 Modify | 0.36932 | 0.36932 | 0.36932 | 0.0 | 0.37 Other | | 0.05722 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 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: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.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 = 313.152824557775, Press = 0.601342933531986 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 = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 33000 -13455.842 -13455.842 -13535.65 -13535.65 308.86628 308.86628 31278.976 31278.976 1986.4914 1986.4914 34000 -13454.014 -13454.014 -13535.67 -13535.67 316.01926 316.01926 31325.37 31325.37 -1845.71 -1845.71 Loop time of 99.4116 on 1 procs for 1000 steps with 2000 atoms Performance: 0.869 ns/day, 27.614 hours/ns, 10.059 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 | 98.913 | 98.913 | 98.913 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.074378 | 0.074378 | 0.074378 | 0.0 | 0.07 Output | 0.00015803 | 0.00015803 | 0.00015803 | 0.0 | 0.00 Modify | 0.36682 | 0.36682 | 0.36682 | 0.0 | 0.37 Other | | 0.05729 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 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: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.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 = 313.149980633851, Press = -1.25880354979782 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 = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 34000 -13454.014 -13454.014 -13535.67 -13535.67 316.01926 316.01926 31325.37 31325.37 -1845.71 -1845.71 35000 -13456.595 -13456.595 -13537.552 -13537.552 313.31342 313.31342 31290.541 31290.541 851.26155 851.26155 Loop time of 100.573 on 1 procs for 1000 steps with 2000 atoms Performance: 0.859 ns/day, 27.937 hours/ns, 9.943 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 | 100.07 | 100.07 | 100.07 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.074781 | 0.074781 | 0.074781 | 0.0 | 0.07 Output | 0.00015306 | 0.00015306 | 0.00015306 | 0.0 | 0.00 Modify | 0.37462 | 0.37462 | 0.37462 | 0.0 | 0.37 Other | | 0.05729 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 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: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.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 = 313.136954449696, Press = -1.20600184137593 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 = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 35000 -13456.595 -13456.595 -13537.552 -13537.552 313.31342 313.31342 31290.541 31290.541 851.26155 851.26155 36000 -13453.637 -13453.637 -13534.863 -13534.863 314.35418 314.35418 31311.724 31311.724 -650.64683 -650.64683 Loop time of 99.8723 on 1 procs for 1000 steps with 2000 atoms Performance: 0.865 ns/day, 27.742 hours/ns, 10.013 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 | 99.367 | 99.367 | 99.367 | 0.0 | 99.49 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.074461 | 0.074461 | 0.074461 | 0.0 | 0.07 Output | 0.0001549 | 0.0001549 | 0.0001549 | 0.0 | 0.00 Modify | 0.37329 | 0.37329 | 0.37329 | 0.0 | 0.37 Other | | 0.05749 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 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: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.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" 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_T313.15.out" else "print 'not_converged' file output/vol_T313.15.out" print '${V}' file output/vol_T313.15.out 31301.866442317 print "LAMMPS calculation completed" LAMMPS calculation completed quit 0