# periodic boundary conditions along all three dimensions boundary p p p # Set neighbor skin variable neigh_skin equal 2.0*${_u_distance} variable neigh_skin equal 2.0*1 neighbor ${neigh_skin} bin neighbor 2 bin # create a supercell with cubic lattice (fcc, bcc, sc, or diamond) # using 10*10*10 conventional (orthogonal) unit cells variable latticeconst_converted equal 2.8886456787586208*${_u_distance} variable latticeconst_converted equal 2.8886456787586208*1 lattice bcc ${latticeconst_converted} lattice bcc 2.88864567875862 Lattice spacing in x,y,z = 2.8886457 2.8886457 2.8886457 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 (28.886457 28.886457 28.886457) 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 (28.886457 28.886457 28.886457) create_atoms CPU = 0.003 seconds variable mass_converted equal 55.845*${_u_mass} variable mass_converted equal 55.845*1 kim_interactions Fe WARNING: 'kim_' has been renamed to 'kim '. Please update your input. kim interactions Fe #=== BEGIN kim interactions ================================== pair_style kim Tersoff_LAMMPS_ByggmastarNagelAlbe_2019_FeO__MO_608695023236_000 pair_coeff * * Fe #=== END kim interactions ==================================== mass 1 ${mass_converted} mass 1 55.845 # initial volume variable v equal vol # assign formula variable V0 equal ${v} # evaluate initial value variable V0 equal 24103.6506205884 variable V0_metal equal ${V0}/(${_u_distance}*${_u_distance}*${_u_distance}) variable V0_metal equal 24103.6506205884/(${_u_distance}*${_u_distance}*${_u_distance}) variable V0_metal equal 24103.6506205884/(1*${_u_distance}*${_u_distance}) variable V0_metal equal 24103.6506205884/(1*1*${_u_distance}) variable V0_metal equal 24103.6506205884/(1*1*1) variable V0_metal_times1000 equal ${V0_metal}*1000 variable V0_metal_times1000 equal 24103.6506205884*1000 print "Initial system volume: ${V0_metal} Angstroms^3" Initial system volume: 24103.6506205884 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 333.15*${_u_temperature} variable temp_converted equal 333.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 333.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 333.15 ${temp_converted} ${Tdamp_converted} iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 333.15 333.15 ${Tdamp_converted} iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 333.15 333.15 0.1 iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 333.15 333.15 0.1 iso 0 ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 333.15 333.15 0.1 iso 0 0 ${Pdamp_converted} fix ensemble all npt temp 333.15 333.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 "333.15 - 0.2" variable T_up equal "333.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_608695023236_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.35 ghost atom cutoff = 5.35 binsize = 2.675, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, ghost, cut 5.35 pair build: full/bin/ghost stencil: full/ghost/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.082 | 4.082 | 4.082 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 0 -8271.2334 -8271.2334 -8357.3164 -8357.3164 333.15 333.15 24103.651 24103.651 3814.6516 3814.6516 1000 -8192.7855 -8192.7855 -8271.3453 -8271.3453 304.03446 304.03446 24112.884 24112.884 -1212.9414 -1212.9414 Loop time of 379.5 on 1 procs for 1000 steps with 2000 atoms Performance: 0.228 ns/day, 105.417 hours/ns, 2.635 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 | 379.01 | 379.01 | 379.01 | 0.0 | 99.87 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.077059 | 0.077059 | 0.077059 | 0.0 | 0.02 Output | 0.00019386 | 0.00019386 | 0.00019386 | 0.0 | 0.00 Modify | 0.35743 | 0.35743 | 0.35743 | 0.0 | 0.09 Other | | 0.05719 | | | 0.02 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.35 ghost atom cutoff = 5.35 binsize = 2.675, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, ghost, cut 5.35 pair build: full/bin/ghost stencil: full/ghost/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.082 | 4.082 | 4.082 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 1000 -8192.7855 -8192.7855 -8271.3453 -8271.3453 304.03446 304.03446 24112.884 24112.884 -1212.9414 -1212.9414 2000 -8186.4405 -8186.4405 -8274.9791 -8274.9791 342.65386 342.65386 24097.872 24097.872 823.90599 823.90599 Loop time of 389.865 on 1 procs for 1000 steps with 2000 atoms Performance: 0.222 ns/day, 108.296 hours/ns, 2.565 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 | 389.36 | 389.36 | 389.36 | 0.0 | 99.87 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.077114 | 0.077114 | 0.077114 | 0.0 | 0.02 Output | 0.00017779 | 0.00017779 | 0.00017779 | 0.0 | 0.00 Modify | 0.36745 | 0.36745 | 0.36745 | 0.0 | 0.09 Other | | 0.05824 | | | 0.01 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.35 ghost atom cutoff = 5.35 binsize = 2.675, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, ghost, cut 5.35 pair build: full/bin/ghost stencil: full/ghost/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.082 | 4.082 | 4.082 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 2000 -8186.4405 -8186.4405 -8274.9791 -8274.9791 342.65386 342.65386 24097.872 24097.872 823.90599 823.90599 3000 -8192.7403 -8192.7403 -8275.2042 -8275.2042 319.14383 319.14383 24102.015 24102.015 536.68052 536.68052 Loop time of 389.4 on 1 procs for 1000 steps with 2000 atoms Performance: 0.222 ns/day, 108.167 hours/ns, 2.568 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 | 388.9 | 388.9 | 388.9 | 0.0 | 99.87 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.077226 | 0.077226 | 0.077226 | 0.0 | 0.02 Output | 0.00019013 | 0.00019013 | 0.00019013 | 0.0 | 0.00 Modify | 0.36487 | 0.36487 | 0.36487 | 0.0 | 0.09 Other | | 0.05832 | | | 0.01 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.35 ghost atom cutoff = 5.35 binsize = 2.675, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, ghost, cut 5.35 pair build: full/bin/ghost stencil: full/ghost/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.082 | 4.082 | 4.082 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 3000 -8192.7403 -8192.7403 -8275.2042 -8275.2042 319.14383 319.14383 24102.015 24102.015 536.68052 536.68052 4000 -8187.311 -8187.311 -8271.9887 -8271.9887 327.71147 327.71147 24125.926 24125.926 -1664.6975 -1664.6975 Loop time of 391.459 on 1 procs for 1000 steps with 2000 atoms Performance: 0.221 ns/day, 108.739 hours/ns, 2.555 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 | 390.95 | 390.95 | 390.95 | 0.0 | 99.87 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.077053 | 0.077053 | 0.077053 | 0.0 | 0.02 Output | 0.00027714 | 0.00027714 | 0.00027714 | 0.0 | 0.00 Modify | 0.3686 | 0.3686 | 0.3686 | 0.0 | 0.09 Other | | 0.0592 | | | 0.02 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.35 ghost atom cutoff = 5.35 binsize = 2.675, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, ghost, cut 5.35 pair build: full/bin/ghost stencil: full/ghost/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.082 | 4.082 | 4.082 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 4000 -8187.311 -8187.311 -8271.9887 -8271.9887 327.71147 327.71147 24125.926 24125.926 -1664.6975 -1664.6975 5000 -8190.3377 -8190.3377 -8275.5325 -8275.5325 329.71266 329.71266 24090.56 24090.56 1308.2116 1308.2116 Loop time of 407.437 on 1 procs for 1000 steps with 2000 atoms Performance: 0.212 ns/day, 113.177 hours/ns, 2.454 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 | 406.91 | 406.91 | 406.91 | 0.0 | 99.87 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.078637 | 0.078637 | 0.078637 | 0.0 | 0.02 Output | 0.00015351 | 0.00015351 | 0.00015351 | 0.0 | 0.00 Modify | 0.38809 | 0.38809 | 0.38809 | 0.0 | 0.10 Other | | 0.06044 | | | 0.01 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 = 331.232953098472, Press = -1363.28475874877 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.35 ghost atom cutoff = 5.35 binsize = 2.675, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, ghost, cut 5.35 pair build: full/bin/ghost stencil: full/ghost/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.082 | 4.082 | 4.082 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 5000 -8190.3377 -8190.3377 -8275.5325 -8275.5325 329.71266 329.71266 24090.56 24090.56 1308.2116 1308.2116 6000 -8188.3374 -8188.3374 -8276.5778 -8276.5778 341.49954 341.49954 24071.557 24071.557 2764.8283 2764.8283 Loop time of 405.568 on 1 procs for 1000 steps with 2000 atoms Performance: 0.213 ns/day, 112.658 hours/ns, 2.466 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 | 405.03 | 405.03 | 405.03 | 0.0 | 99.87 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.079022 | 0.079022 | 0.079022 | 0.0 | 0.02 Output | 0.00017779 | 0.00017779 | 0.00017779 | 0.0 | 0.00 Modify | 0.40188 | 0.40188 | 0.40188 | 0.0 | 0.10 Other | | 0.05985 | | | 0.01 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 = 334.030121707144, Press = 56.7334718292653 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.35 ghost atom cutoff = 5.35 binsize = 2.675, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, ghost, cut 5.35 pair build: full/bin/ghost stencil: full/ghost/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.082 | 4.082 | 4.082 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 6000 -8188.3374 -8188.3374 -8276.5778 -8276.5778 341.49954 341.49954 24071.557 24071.557 2764.8283 2764.8283 7000 -8190.5785 -8190.5785 -8273.8258 -8273.8258 322.1754 322.1754 24146.994 24146.994 -3642.4472 -3642.4472 Loop time of 407.201 on 1 procs for 1000 steps with 2000 atoms Performance: 0.212 ns/day, 113.111 hours/ns, 2.456 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 | 406.66 | 406.66 | 406.66 | 0.0 | 99.87 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.07834 | 0.07834 | 0.07834 | 0.0 | 0.02 Output | 0.00024311 | 0.00024311 | 0.00024311 | 0.0 | 0.00 Modify | 0.40302 | 0.40302 | 0.40302 | 0.0 | 0.10 Other | | 0.06025 | | | 0.01 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 = 333.507572245315, Press = -33.7181950686494 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.35 ghost atom cutoff = 5.35 binsize = 2.675, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, ghost, cut 5.35 pair build: full/bin/ghost stencil: full/ghost/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.082 | 4.082 | 4.082 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 7000 -8190.5785 -8190.5785 -8273.8258 -8273.8258 322.1754 322.1754 24146.994 24146.994 -3642.4472 -3642.4472 8000 -8184.1284 -8184.1284 -8273.678 -8273.678 346.56622 346.56622 24113.768 24113.768 -300.76175 -300.76175 Loop time of 390.39 on 1 procs for 1000 steps with 2000 atoms Performance: 0.221 ns/day, 108.442 hours/ns, 2.562 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 | 389.87 | 389.87 | 389.87 | 0.0 | 99.87 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.075851 | 0.075851 | 0.075851 | 0.0 | 0.02 Output | 0.00018658 | 0.00018658 | 0.00018658 | 0.0 | 0.00 Modify | 0.3832 | 0.3832 | 0.3832 | 0.0 | 0.10 Other | | 0.05817 | | | 0.01 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 = 334.146432941083, Press = -22.7477877686432 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.35 ghost atom cutoff = 5.35 binsize = 2.675, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, ghost, cut 5.35 pair build: full/bin/ghost stencil: full/ghost/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.082 | 4.082 | 4.082 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 8000 -8184.1284 -8184.1284 -8273.678 -8273.678 346.56622 346.56622 24113.768 24113.768 -300.76175 -300.76175 9000 -8190.3556 -8190.3556 -8275.1129 -8275.1129 328.01964 328.01964 24091.228 24091.228 1455.6669 1455.6669 Loop time of 411.603 on 1 procs for 1000 steps with 2000 atoms Performance: 0.210 ns/day, 114.334 hours/ns, 2.430 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 | 411.05 | 411.05 | 411.05 | 0.0 | 99.87 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.079846 | 0.079846 | 0.079846 | 0.0 | 0.02 Output | 0.000182 | 0.000182 | 0.000182 | 0.0 | 0.00 Modify | 0.4089 | 0.4089 | 0.4089 | 0.0 | 0.10 Other | | 0.06045 | | | 0.01 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 = 334.325594066276, Press = -25.5354236418228 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.35 ghost atom cutoff = 5.35 binsize = 2.675, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, ghost, cut 5.35 pair build: full/bin/ghost stencil: full/ghost/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.082 | 4.082 | 4.082 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 9000 -8190.3556 -8190.3556 -8275.1129 -8275.1129 328.01964 328.01964 24091.228 24091.228 1455.6669 1455.6669 10000 -8189.0855 -8189.0855 -8275.451 -8275.451 334.24366 334.24366 24108.476 24108.476 477.7328 477.7328 Loop time of 405.908 on 1 procs for 1000 steps with 2000 atoms Performance: 0.213 ns/day, 112.752 hours/ns, 2.464 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 | 405.36 | 405.36 | 405.36 | 0.0 | 99.87 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.078801 | 0.078801 | 0.078801 | 0.0 | 0.02 Output | 0.00015637 | 0.00015637 | 0.00015637 | 0.0 | 0.00 Modify | 0.40449 | 0.40449 | 0.40449 | 0.0 | 0.10 Other | | 0.06018 | | | 0.01 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 = 334.55536227997, Press = -0.283829958511007 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.35 ghost atom cutoff = 5.35 binsize = 2.675, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, ghost, cut 5.35 pair build: full/bin/ghost stencil: full/ghost/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.082 | 4.082 | 4.082 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 10000 -8189.0855 -8189.0855 -8275.451 -8275.451 334.24366 334.24366 24108.476 24108.476 477.7328 477.7328 11000 -8189.8618 -8189.8618 -8277.8398 -8277.8398 340.48391 340.48391 24165.072 24165.072 -4479.7992 -4479.7992 Loop time of 407.24 on 1 procs for 1000 steps with 2000 atoms Performance: 0.212 ns/day, 113.122 hours/ns, 2.456 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 | 406.7 | 406.7 | 406.7 | 0.0 | 99.87 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.078874 | 0.078874 | 0.078874 | 0.0 | 0.02 Output | 0.00015336 | 0.00015336 | 0.00015336 | 0.0 | 0.00 Modify | 0.40274 | 0.40274 | 0.40274 | 0.0 | 0.10 Other | | 0.06036 | | | 0.01 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 = 334.316770413371, Press = -13.3739740565999 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.35 ghost atom cutoff = 5.35 binsize = 2.675, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, ghost, cut 5.35 pair build: full/bin/ghost stencil: full/ghost/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.082 | 4.082 | 4.082 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 11000 -8189.8618 -8189.8618 -8277.8398 -8277.8398 340.48391 340.48391 24165.072 24165.072 -4479.7992 -4479.7992 12000 -8189.3312 -8189.3312 -8277.1573 -8277.1573 339.89634 339.89634 24088.1 24088.1 2226.3193 2226.3193 Loop time of 409.361 on 1 procs for 1000 steps with 2000 atoms Performance: 0.211 ns/day, 113.711 hours/ns, 2.443 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 | 408.81 | 408.81 | 408.81 | 0.0 | 99.87 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.079534 | 0.079534 | 0.079534 | 0.0 | 0.02 Output | 0.00015352 | 0.00015352 | 0.00015352 | 0.0 | 0.00 Modify | 0.40873 | 0.40873 | 0.40873 | 0.0 | 0.10 Other | | 0.06042 | | | 0.01 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 = 334.174079057124, Press = -18.4187523402845 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.35 ghost atom cutoff = 5.35 binsize = 2.675, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, ghost, cut 5.35 pair build: full/bin/ghost stencil: full/ghost/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.082 | 4.082 | 4.082 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 12000 -8189.3312 -8189.3312 -8277.1573 -8277.1573 339.89634 339.89634 24088.1 24088.1 2226.3193 2226.3193 13000 -8191.4036 -8191.4036 -8276.5952 -8276.5952 329.70041 329.70041 24094.782 24094.782 1359.0886 1359.0886 Loop time of 405.561 on 1 procs for 1000 steps with 2000 atoms Performance: 0.213 ns/day, 112.656 hours/ns, 2.466 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 | 405.02 | 405.02 | 405.02 | 0.0 | 99.87 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.07837 | 0.07837 | 0.07837 | 0.0 | 0.02 Output | 0.00015278 | 0.00015278 | 0.00015278 | 0.0 | 0.00 Modify | 0.40097 | 0.40097 | 0.40097 | 0.0 | 0.10 Other | | 0.05992 | | | 0.01 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 = 333.86726339787, Press = 1.3446998745073 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.35 ghost atom cutoff = 5.35 binsize = 2.675, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, ghost, cut 5.35 pair build: full/bin/ghost stencil: full/ghost/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.082 | 4.082 | 4.082 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 13000 -8191.4036 -8191.4036 -8276.5952 -8276.5952 329.70041 329.70041 24094.782 24094.782 1359.0886 1359.0886 14000 -8194.0739 -8194.0739 -8278.1936 -8278.1936 325.5521 325.5521 24109.772 24109.772 388.27663 388.27663 Loop time of 410.279 on 1 procs for 1000 steps with 2000 atoms Performance: 0.211 ns/day, 113.966 hours/ns, 2.437 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 | 409.74 | 409.74 | 409.74 | 0.0 | 99.87 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.07879 | 0.07879 | 0.07879 | 0.0 | 0.02 Output | 0.00015134 | 0.00015134 | 0.00015134 | 0.0 | 0.00 Modify | 0.40163 | 0.40163 | 0.40163 | 0.0 | 0.10 Other | | 0.05992 | | | 0.01 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 = 333.765457669982, Press = -4.38181109251149 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.35 ghost atom cutoff = 5.35 binsize = 2.675, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, ghost, cut 5.35 pair build: full/bin/ghost stencil: full/ghost/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.082 | 4.082 | 4.082 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 14000 -8194.0739 -8194.0739 -8278.1936 -8278.1936 325.5521 325.5521 24109.772 24109.772 388.27663 388.27663 15000 -8189.1514 -8189.1514 -8274.2374 -8274.2374 329.29144 329.29144 24129.523 24129.523 -2103.8758 -2103.8758 Loop time of 407.658 on 1 procs for 1000 steps with 2000 atoms Performance: 0.212 ns/day, 113.238 hours/ns, 2.453 timesteps/s 100.0% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 407.11 | 407.11 | 407.11 | 0.0 | 99.87 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.080019 | 0.080019 | 0.080019 | 0.0 | 0.02 Output | 0.00015569 | 0.00015569 | 0.00015569 | 0.0 | 0.00 Modify | 0.4091 | 0.4091 | 0.4091 | 0.0 | 0.10 Other | | 0.06032 | | | 0.01 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 = 333.405502701112, Press = -9.58429867291999 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.35 ghost atom cutoff = 5.35 binsize = 2.675, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, ghost, cut 5.35 pair build: full/bin/ghost stencil: full/ghost/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.082 | 4.082 | 4.082 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 15000 -8189.1514 -8189.1514 -8274.2374 -8274.2374 329.29144 329.29144 24129.523 24129.523 -2103.8758 -2103.8758 16000 -8195.298 -8195.298 -8276.9024 -8276.9024 315.81762 315.81762 24102.867 24102.867 -168.33519 -168.33519 Loop time of 404.972 on 1 procs for 1000 steps with 2000 atoms Performance: 0.213 ns/day, 112.492 hours/ns, 2.469 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 | 404.43 | 404.43 | 404.43 | 0.0 | 99.87 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.078278 | 0.078278 | 0.078278 | 0.0 | 0.02 Output | 0.00015599 | 0.00015599 | 0.00015599 | 0.0 | 0.00 Modify | 0.39899 | 0.39899 | 0.39899 | 0.0 | 0.10 Other | | 0.05949 | | | 0.01 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 = 333.236563824655, Press = -2.39657940241327 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.35 ghost atom cutoff = 5.35 binsize = 2.675, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, ghost, cut 5.35 pair build: full/bin/ghost stencil: full/ghost/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.082 | 4.082 | 4.082 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 16000 -8195.298 -8195.298 -8276.9024 -8276.9024 315.81762 315.81762 24102.867 24102.867 -168.33519 -168.33519 17000 -8188.6078 -8188.6078 -8275.0216 -8275.0216 334.43037 334.43037 24109.707 24109.707 -347.91797 -347.91797 Loop time of 413.05 on 1 procs for 1000 steps with 2000 atoms Performance: 0.209 ns/day, 114.736 hours/ns, 2.421 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 | 412.5 | 412.5 | 412.5 | 0.0 | 99.87 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.079822 | 0.079822 | 0.079822 | 0.0 | 0.02 Output | 0.00015314 | 0.00015314 | 0.00015314 | 0.0 | 0.00 Modify | 0.41209 | 0.41209 | 0.41209 | 0.0 | 0.10 Other | | 0.06003 | | | 0.01 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 = 333.13441796445, Press = -9.56009749885033 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.35 ghost atom cutoff = 5.35 binsize = 2.675, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, ghost, cut 5.35 pair build: full/bin/ghost stencil: full/ghost/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.082 | 4.082 | 4.082 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 17000 -8188.6078 -8188.6078 -8275.0216 -8275.0216 334.43037 334.43037 24109.707 24109.707 -347.91797 -347.91797 18000 -8185.3778 -8185.3778 -8273.439 -8273.439 340.80607 340.80607 24072.493 24072.493 2636.4075 2636.4075 Loop time of 413.277 on 1 procs for 1000 steps with 2000 atoms Performance: 0.209 ns/day, 114.799 hours/ns, 2.420 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 | 412.72 | 412.72 | 412.72 | 0.0 | 99.86 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.080547 | 0.080547 | 0.080547 | 0.0 | 0.02 Output | 0.00015589 | 0.00015589 | 0.00015589 | 0.0 | 0.00 Modify | 0.41814 | 0.41814 | 0.41814 | 0.0 | 0.10 Other | | 0.06092 | | | 0.01 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 = 333.141962785027, Press = 3.60333907324603 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.35 ghost atom cutoff = 5.35 binsize = 2.675, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, ghost, cut 5.35 pair build: full/bin/ghost stencil: full/ghost/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.082 | 4.082 | 4.082 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 18000 -8185.3778 -8185.3778 -8273.439 -8273.439 340.80607 340.80607 24072.493 24072.493 2636.4075 2636.4075 19000 -8190.5292 -8190.5292 -8276.5123 -8276.5123 332.76332 332.76332 24149.104 24149.104 -3191.6768 -3191.6768 Loop time of 414.5 on 1 procs for 1000 steps with 2000 atoms Performance: 0.208 ns/day, 115.139 hours/ns, 2.413 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 | 413.94 | 413.94 | 413.94 | 0.0 | 99.87 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.080011 | 0.080011 | 0.080011 | 0.0 | 0.02 Output | 0.00015415 | 0.00015415 | 0.00015415 | 0.0 | 0.00 Modify | 0.41662 | 0.41662 | 0.41662 | 0.0 | 0.10 Other | | 0.06175 | | | 0.01 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 = 333.279264149599, Press = -4.81155303932463 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.35 ghost atom cutoff = 5.35 binsize = 2.675, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, ghost, cut 5.35 pair build: full/bin/ghost stencil: full/ghost/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.082 | 4.082 | 4.082 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 19000 -8190.5292 -8190.5292 -8276.5123 -8276.5123 332.76332 332.76332 24149.104 24149.104 -3191.6768 -3191.6768 20000 -8188.6298 -8188.6298 -8275.062 -8275.062 334.50127 334.50127 24082.436 24082.436 2065.7697 2065.7697 Loop time of 406.248 on 1 procs for 1000 steps with 2000 atoms Performance: 0.213 ns/day, 112.847 hours/ns, 2.462 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 | 405.7 | 405.7 | 405.7 | 0.0 | 99.87 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.078795 | 0.078795 | 0.078795 | 0.0 | 0.02 Output | 0.00015448 | 0.00015448 | 0.00015448 | 0.0 | 0.00 Modify | 0.40607 | 0.40607 | 0.40607 | 0.0 | 0.10 Other | | 0.06035 | | | 0.01 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 = 333.333240549974, Press = -4.92743879361674 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.35 ghost atom cutoff = 5.35 binsize = 2.675, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, ghost, cut 5.35 pair build: full/bin/ghost stencil: full/ghost/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.082 | 4.082 | 4.082 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 20000 -8188.6298 -8188.6298 -8275.062 -8275.062 334.50127 334.50127 24082.436 24082.436 2065.7697 2065.7697 21000 -8194.0927 -8194.0927 -8277.971 -8277.971 324.61775 324.61775 24120.831 24120.831 -1666.352 -1666.352 Loop time of 407.385 on 1 procs for 1000 steps with 2000 atoms Performance: 0.212 ns/day, 113.163 hours/ns, 2.455 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 | 406.84 | 406.84 | 406.84 | 0.0 | 99.87 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.078735 | 0.078735 | 0.078735 | 0.0 | 0.02 Output | 0.000155 | 0.000155 | 0.000155 | 0.0 | 0.00 Modify | 0.40577 | 0.40577 | 0.40577 | 0.0 | 0.10 Other | | 0.06023 | | | 0.01 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 = 333.25413733737, Press = -1.97913154491983 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.35 ghost atom cutoff = 5.35 binsize = 2.675, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, ghost, cut 5.35 pair build: full/bin/ghost stencil: full/ghost/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.082 | 4.082 | 4.082 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 21000 -8194.0927 -8194.0927 -8277.971 -8277.971 324.61775 324.61775 24120.831 24120.831 -1666.352 -1666.352 22000 -8188.7282 -8188.7282 -8272.0785 -8272.0785 322.57439 322.57439 24115.547 24115.547 -697.41352 -697.41352 Loop time of 409.895 on 1 procs for 1000 steps with 2000 atoms Performance: 0.211 ns/day, 113.860 hours/ns, 2.440 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 | 409.34 | 409.34 | 409.34 | 0.0 | 99.87 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.07965 | 0.07965 | 0.07965 | 0.0 | 0.02 Output | 0.00015257 | 0.00015257 | 0.00015257 | 0.0 | 0.00 Modify | 0.4097 | 0.4097 | 0.4097 | 0.0 | 0.10 Other | | 0.06068 | | | 0.01 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 = 333.139076237123, Press = -5.46295225233058 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.35 ghost atom cutoff = 5.35 binsize = 2.675, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, ghost, cut 5.35 pair build: full/bin/ghost stencil: full/ghost/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.082 | 4.082 | 4.082 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 22000 -8188.7282 -8188.7282 -8272.0785 -8272.0785 322.57439 322.57439 24115.547 24115.547 -697.41352 -697.41352 23000 -8192.8127 -8192.8127 -8276.0104 -8276.0104 321.98367 321.98367 24106.307 24106.307 236.50755 236.50755 Loop time of 403.028 on 1 procs for 1000 steps with 2000 atoms Performance: 0.214 ns/day, 111.952 hours/ns, 2.481 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 | 402.49 | 402.49 | 402.49 | 0.0 | 99.87 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.078887 | 0.078887 | 0.078887 | 0.0 | 0.02 Output | 0.00015756 | 0.00015756 | 0.00015756 | 0.0 | 0.00 Modify | 0.39646 | 0.39646 | 0.39646 | 0.0 | 0.10 Other | | 0.05934 | | | 0.01 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 = 332.99003084664, Press = -5.23005645520914 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.35 ghost atom cutoff = 5.35 binsize = 2.675, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, ghost, cut 5.35 pair build: full/bin/ghost stencil: full/ghost/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.082 | 4.082 | 4.082 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 23000 -8192.8127 -8192.8127 -8276.0104 -8276.0104 321.98367 321.98367 24106.307 24106.307 236.50755 236.50755 24000 -8191.4799 -8191.4799 -8276.9795 -8276.9795 330.89225 330.89225 24095.979 24095.979 1289.8733 1289.8733 Loop time of 399.3 on 1 procs for 1000 steps with 2000 atoms Performance: 0.216 ns/day, 110.917 hours/ns, 2.504 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 | 398.77 | 398.77 | 398.77 | 0.0 | 99.87 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.078137 | 0.078137 | 0.078137 | 0.0 | 0.02 Output | 0.00015416 | 0.00015416 | 0.00015416 | 0.0 | 0.00 Modify | 0.39577 | 0.39577 | 0.39577 | 0.0 | 0.10 Other | | 0.05943 | | | 0.01 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_T333.15.out" else "print 'not_converged' file output/vol_T333.15.out" print '${V}' file output/vol_T333.15.out 24109.2632622817 print "LAMMPS calculation completed" LAMMPS calculation completed quit 0