# 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.863657325506211*${_u_distance} variable latticeconst_converted equal 2.863657325506211*1 lattice bcc ${latticeconst_converted} lattice bcc 2.86365732550621 Lattice spacing in x,y,z = 2.8636573 2.8636573 2.8636573 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.636573 28.636573 28.636573) 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.636573 28.636573 28.636573) create_atoms CPU = 0.004 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 MEAM_LAMMPS_WuLeeSu_2017_NiCrFe__MO_912636107108_001 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 23483.5171943885 variable V0_metal equal ${V0}/(${_u_distance}*${_u_distance}*${_u_distance}) variable V0_metal equal 23483.5171943885/(${_u_distance}*${_u_distance}*${_u_distance}) variable V0_metal equal 23483.5171943885/(1*${_u_distance}*${_u_distance}) variable V0_metal equal 23483.5171943885/(1*1*${_u_distance}) variable V0_metal equal 23483.5171943885/(1*1*1) variable V0_metal_times1000 equal ${V0_metal}*1000 variable V0_metal_times1000 equal 23483.5171943885*1000 print "Initial system volume: ${V0_metal} Angstroms^3" Initial system volume: 23483.5171943885 Angstroms^3 # set the time step to 0.001 picoseconds variable timestep_converted equal 0.001*${_u_time} variable timestep_converted equal 0.001*1 timestep ${timestep_converted} timestep 0.001 variable temp_converted equal 293.15*${_u_temperature} variable temp_converted equal 293.15*1 variable Tdamp_converted equal 0.1*${_u_time} variable Tdamp_converted equal 0.1*1 variable press_converted equal 0.0*${_u_pressure} variable press_converted equal 0.0*1 variable Pdamp_converted equal 1*${_u_time} variable Pdamp_converted equal 1*1 # create initial velocities consistent with the chosen temperature velocity all create ${temp_converted} 17 mom yes rot yes velocity all create 293.15 17 mom yes rot yes # set NPT ensemble for all atoms fix ensemble all npt temp ${temp_converted} ${temp_converted} ${Tdamp_converted} iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 293.15 ${temp_converted} ${Tdamp_converted} iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 293.15 293.15 ${Tdamp_converted} iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 293.15 293.15 0.1 iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 293.15 293.15 0.1 iso 0 ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 293.15 293.15 0.1 iso 0 0 ${Pdamp_converted} fix ensemble all npt temp 293.15 293.15 0.1 iso 0 0 1 # compute the time averages of pressure, temperature, and volume, respectively # ignore the first 5000 timesteps variable etotal_metal equal etotal/${_u_energy} variable etotal_metal equal etotal/1 variable pe_metal equal pe/${_u_energy} variable pe_metal equal pe/1 variable T_metal equal temp/${_u_temperature} variable T_metal equal temp/1 variable V_metal equal vol/(${_u_distance}*${_u_distance}*${_u_distance}) variable V_metal equal vol/(1*${_u_distance}*${_u_distance}) variable V_metal equal vol/(1*1*${_u_distance}) variable V_metal equal vol/(1*1*1) variable P_metal equal press/${_u_pressure} variable P_metal equal press/1 fix avgmyTemp all ave/time 5 20 100 v_T_metal ave running start 5000 fix avgmyPress all ave/time 5 20 100 v_P_metal ave running start 5000 fix avgmyVol all ave/time 5 20 100 v_V_metal ave running start 5000 # extract fix quantities into variables so they can be used in if-else logic later. variable T equal f_avgmyTemp variable P equal f_avgmyPress variable V equal f_avgmyVol # set error bounds for temperature and pressure in original metal units (K and bar) variable T_low equal "293.15 - 0.2" variable T_up equal "293.15 + 0.2" variable P_low equal "0.0 - 0.2" variable P_up equal "0.0 + 0.2" # print to logfile every 1000 timesteps thermo_style custom step etotal v_etotal_metal pe v_pe_metal temp v_T_metal vol v_V_metal press v_P_metal thermo 1000 # Run a simulation for at most 2000*1000 timesteps. At each 1000th time step, check # whether the temperature and pressure have converged. If yes, break. label top variable a loop 2000 run 1000 CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE Your simulation uses code contributions which should be cited: - OpenKIM: https://doi.org/10.1007/s11837-011-0102-6 @Article{tadmor:elliott:2011, author = {E. B. Tadmor and R. S. Elliott and J. P. Sethna and R. E. Miller and C. A. Becker}, title = {The potential of atomistic simulations and the {K}nowledgebase of {I}nteratomic {M}odels}, journal = {{JOM}}, year = 2011, volume = 63, number = 17, pages = {17}, doi = {10.1007/s11837-011-0102-6} } - OpenKIM potential: https://openkim.org/cite/MO_912636107108_001#item-citation CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 10 10 10 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 0 -8504.2527 -8504.2527 -8580 -8580 293.15 293.15 23483.517 23483.517 3445.2596 3445.2596 1000 -8425.0239 -8425.0239 -8497.9604 -8497.9604 282.27172 282.27172 23692.277 23692.277 2408.0405 2408.0405 Loop time of 144.547 on 1 procs for 1000 steps with 2000 atoms Performance: 0.598 ns/day, 40.152 hours/ns, 6.918 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 | 144.03 | 144.03 | 144.03 | 0.0 | 99.64 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.098164 | 0.098164 | 0.098164 | 0.0 | 0.07 Output | 0.00019451 | 0.00019451 | 0.00019451 | 0.0 | 0.00 Modify | 0.35556 | 0.35556 | 0.35556 | 0.0 | 0.25 Other | | 0.06043 | | | 0.04 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119.00 ave 4119 max 4119 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: 128000.0 ave 128000 max 128000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 128000 Ave neighs/atom = 64.000000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 0, Press = 0 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 10 10 10 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 1000 -8425.0239 -8425.0239 -8497.9604 -8497.9604 282.27172 282.27172 23692.277 23692.277 2408.0405 2408.0405 2000 -8425.4004 -8425.4004 -8500.2416 -8500.2416 289.64312 289.64312 23725.037 23725.037 -178.52369 -178.52369 Loop time of 151.319 on 1 procs for 1000 steps with 2000 atoms Performance: 0.571 ns/day, 42.033 hours/ns, 6.609 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 | 150.81 | 150.81 | 150.81 | 0.0 | 99.66 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.097232 | 0.097232 | 0.097232 | 0.0 | 0.06 Output | 0.00018464 | 0.00018464 | 0.00018464 | 0.0 | 0.00 Modify | 0.35142 | 0.35142 | 0.35142 | 0.0 | 0.23 Other | | 0.0604 | | | 0.04 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119.00 ave 4119 max 4119 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: 128062.0 ave 128062 max 128062 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 128062 Ave neighs/atom = 64.031000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 0, Press = 0 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 10 10 10 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 2000 -8425.4004 -8425.4004 -8500.2416 -8500.2416 289.64312 289.64312 23725.037 23725.037 -178.52369 -178.52369 3000 -8426.809 -8426.809 -8497.7632 -8497.7632 274.60014 274.60014 23739.586 23739.586 -1024.946 -1024.946 Loop time of 154.341 on 1 procs for 1000 steps with 2000 atoms Performance: 0.560 ns/day, 42.873 hours/ns, 6.479 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 | 153.82 | 153.82 | 153.82 | 0.0 | 99.66 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.099329 | 0.099329 | 0.099329 | 0.0 | 0.06 Output | 0.00023005 | 0.00023005 | 0.00023005 | 0.0 | 0.00 Modify | 0.36053 | 0.36053 | 0.36053 | 0.0 | 0.23 Other | | 0.06128 | | | 0.04 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119.00 ave 4119 max 4119 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: 128078.0 ave 128078 max 128078 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 128078 Ave neighs/atom = 64.039000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 0, Press = 0 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 10 10 10 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 3000 -8426.809 -8426.809 -8497.7632 -8497.7632 274.60014 274.60014 23739.586 23739.586 -1024.946 -1024.946 4000 -8424.0945 -8424.0945 -8502.4446 -8502.4446 303.2227 303.2227 23759.323 23759.323 -2785.9023 -2785.9023 Loop time of 153.52 on 1 procs for 1000 steps with 2000 atoms Performance: 0.563 ns/day, 42.644 hours/ns, 6.514 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 | 153 | 153 | 153 | 0.0 | 99.66 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.098768 | 0.098768 | 0.098768 | 0.0 | 0.06 Output | 0.00019287 | 0.00019287 | 0.00019287 | 0.0 | 0.00 Modify | 0.36219 | 0.36219 | 0.36219 | 0.0 | 0.24 Other | | 0.06145 | | | 0.04 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119.00 ave 4119 max 4119 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: 128040.0 ave 128040 max 128040 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 128040 Ave neighs/atom = 64.020000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 0, Press = 0 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 10 10 10 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 4000 -8424.0945 -8424.0945 -8502.4446 -8502.4446 303.2227 303.2227 23759.323 23759.323 -2785.9023 -2785.9023 5000 -8427.0408 -8427.0408 -8501.1089 -8501.1089 286.65155 286.65155 23754.254 23754.254 -2415.6535 -2415.6535 Loop time of 154.932 on 1 procs for 1000 steps with 2000 atoms Performance: 0.558 ns/day, 43.037 hours/ns, 6.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 | 154.4 | 154.4 | 154.4 | 0.0 | 99.66 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.099511 | 0.099511 | 0.099511 | 0.0 | 0.06 Output | 0.0001545 | 0.0001545 | 0.0001545 | 0.0 | 0.00 Modify | 0.3706 | 0.3706 | 0.3706 | 0.0 | 0.24 Other | | 0.06187 | | | 0.04 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119.00 ave 4119 max 4119 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: 128008.0 ave 128008 max 128008 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 128008 Ave neighs/atom = 64.004000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 293.418487576325, Press = -145.778826926322 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 10 10 10 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 5000 -8427.0408 -8427.0408 -8501.1089 -8501.1089 286.65155 286.65155 23754.254 23754.254 -2415.6535 -2415.6535 6000 -8425.4554 -8425.4554 -8501.9076 -8501.9076 295.87769 295.87769 23740.732 23740.732 -1536.6454 -1536.6454 Loop time of 152.942 on 1 procs for 1000 steps with 2000 atoms Performance: 0.565 ns/day, 42.484 hours/ns, 6.538 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 | 152.4 | 152.4 | 152.4 | 0.0 | 99.65 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.098437 | 0.098437 | 0.098437 | 0.0 | 0.06 Output | 0.00015298 | 0.00015298 | 0.00015298 | 0.0 | 0.00 Modify | 0.38151 | 0.38151 | 0.38151 | 0.0 | 0.25 Other | | 0.06081 | | | 0.04 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119.00 ave 4119 max 4119 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: 128024.0 ave 128024 max 128024 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 128024 Ave neighs/atom = 64.012000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 292.012596898415, Press = -19.8350405041784 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 10 10 10 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 6000 -8425.4554 -8425.4554 -8501.9076 -8501.9076 295.87769 295.87769 23740.732 23740.732 -1536.6454 -1536.6454 7000 -8422.8871 -8422.8871 -8498.947 -8498.947 294.35961 294.35961 23721.445 23721.445 286.40608 286.40608 Loop time of 153.936 on 1 procs for 1000 steps with 2000 atoms Performance: 0.561 ns/day, 42.760 hours/ns, 6.496 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 | 153.39 | 153.39 | 153.39 | 0.0 | 99.65 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.099163 | 0.099163 | 0.099163 | 0.0 | 0.06 Output | 0.0001547 | 0.0001547 | 0.0001547 | 0.0 | 0.00 Modify | 0.38263 | 0.38263 | 0.38263 | 0.0 | 0.25 Other | | 0.06145 | | | 0.04 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4115.00 ave 4115 max 4115 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: 128030.0 ave 128030 max 128030 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 128030 Ave neighs/atom = 64.015000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 292.82059838539, Press = -14.6603800117828 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 10 10 10 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 7000 -8422.8871 -8422.8871 -8498.947 -8498.947 294.35961 294.35961 23721.445 23721.445 286.40608 286.40608 8000 -8426.7697 -8426.7697 -8502.4697 -8502.4697 292.96668 292.96668 23701.019 23701.019 1239.6316 1239.6316 Loop time of 154.045 on 1 procs for 1000 steps with 2000 atoms Performance: 0.561 ns/day, 42.790 hours/ns, 6.492 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 | 153.5 | 153.5 | 153.5 | 0.0 | 99.65 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.098643 | 0.098643 | 0.098643 | 0.0 | 0.06 Output | 0.00018631 | 0.00018631 | 0.00018631 | 0.0 | 0.00 Modify | 0.38001 | 0.38001 | 0.38001 | 0.0 | 0.25 Other | | 0.0613 | | | 0.04 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4113.00 ave 4113 max 4113 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: 128026.0 ave 128026 max 128026 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 128026 Ave neighs/atom = 64.013000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 293.160490454784, Press = -13.1959297594132 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 10 10 10 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 8000 -8426.7697 -8426.7697 -8502.4697 -8502.4697 292.96668 292.96668 23701.019 23701.019 1239.6316 1239.6316 9000 -8423.0508 -8423.0508 -8501.886 -8501.886 305.10066 305.10066 23672.344 23672.344 3499.3852 3499.3852 Loop time of 152.291 on 1 procs for 1000 steps with 2000 atoms Performance: 0.567 ns/day, 42.303 hours/ns, 6.566 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 | 151.76 | 151.76 | 151.76 | 0.0 | 99.65 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.097974 | 0.097974 | 0.097974 | 0.0 | 0.06 Output | 0.0001573 | 0.0001573 | 0.0001573 | 0.0 | 0.00 Modify | 0.37461 | 0.37461 | 0.37461 | 0.0 | 0.25 Other | | 0.06069 | | | 0.04 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119.00 ave 4119 max 4119 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: 128024.0 ave 128024 max 128024 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 128024 Ave neighs/atom = 64.012000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 293.542756769594, Press = -3.92188932226494 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 10 10 10 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 9000 -8423.0508 -8423.0508 -8501.886 -8501.886 305.10066 305.10066 23672.344 23672.344 3499.3852 3499.3852 10000 -8426.722 -8426.722 -8501.6684 -8501.6684 290.05008 290.05008 23682.886 23682.886 2587.1381 2587.1381 Loop time of 152.097 on 1 procs for 1000 steps with 2000 atoms Performance: 0.568 ns/day, 42.249 hours/ns, 6.575 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 | 151.56 | 151.56 | 151.56 | 0.0 | 99.65 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.097201 | 0.097201 | 0.097201 | 0.0 | 0.06 Output | 0.0001522 | 0.0001522 | 0.0001522 | 0.0 | 0.00 Modify | 0.37316 | 0.37316 | 0.37316 | 0.0 | 0.25 Other | | 0.06173 | | | 0.04 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119.00 ave 4119 max 4119 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: 128044.0 ave 128044 max 128044 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 128044 Ave neighs/atom = 64.022000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 293.817224692678, Press = 7.52927974243659 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 10 10 10 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 10000 -8426.722 -8426.722 -8501.6684 -8501.6684 290.05008 290.05008 23682.886 23682.886 2587.1381 2587.1381 11000 -8422.0396 -8422.0396 -8499.8211 -8499.8211 301.02269 301.02269 23710.976 23710.976 1133.9934 1133.9934 Loop time of 153.015 on 1 procs for 1000 steps with 2000 atoms Performance: 0.565 ns/day, 42.504 hours/ns, 6.535 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 | 152.48 | 152.48 | 152.48 | 0.0 | 99.65 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.098929 | 0.098929 | 0.098929 | 0.0 | 0.06 Output | 0.00015124 | 0.00015124 | 0.00015124 | 0.0 | 0.00 Modify | 0.37777 | 0.37777 | 0.37777 | 0.0 | 0.25 Other | | 0.06073 | | | 0.04 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119.00 ave 4119 max 4119 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: 128056.0 ave 128056 max 128056 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 128056 Ave neighs/atom = 64.028000 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 = 294.030239867756, Press = 6.74982819092394 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 10 10 10 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 11000 -8422.0396 -8422.0396 -8499.8211 -8499.8211 301.02269 301.02269 23710.976 23710.976 1133.9934 1133.9934 12000 -8426.1699 -8426.1699 -8500.451 -8500.451 287.47565 287.47565 23722.217 23722.217 -91.624173 -91.624173 Loop time of 151.477 on 1 procs for 1000 steps with 2000 atoms Performance: 0.570 ns/day, 42.077 hours/ns, 6.602 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 | 150.95 | 150.95 | 150.95 | 0.0 | 99.65 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.096712 | 0.096712 | 0.096712 | 0.0 | 0.06 Output | 0.00015489 | 0.00015489 | 0.00015489 | 0.0 | 0.00 Modify | 0.37268 | 0.37268 | 0.37268 | 0.0 | 0.25 Other | | 0.06073 | | | 0.04 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119.00 ave 4119 max 4119 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: 128080.0 ave 128080 max 128080 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 128080 Ave neighs/atom = 64.040000 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 = 294.131445752861, Press = 9.46603413193528 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 10 10 10 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 12000 -8426.1699 -8426.1699 -8500.451 -8500.451 287.47565 287.47565 23722.217 23722.217 -91.624173 -91.624173 13000 -8420.6681 -8420.6681 -8499.5671 -8499.5671 305.34744 305.34744 23732.81 23732.81 -289.45182 -289.45182 Loop time of 150.021 on 1 procs for 1000 steps with 2000 atoms Performance: 0.576 ns/day, 41.672 hours/ns, 6.666 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 | 149.5 | 149.5 | 149.5 | 0.0 | 99.65 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.096319 | 0.096319 | 0.096319 | 0.0 | 0.06 Output | 0.0001546 | 0.0001546 | 0.0001546 | 0.0 | 0.00 Modify | 0.3658 | 0.3658 | 0.3658 | 0.0 | 0.24 Other | | 0.06019 | | | 0.04 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4118.00 ave 4118 max 4118 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: 128034.0 ave 128034 max 128034 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 128034 Ave neighs/atom = 64.017000 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 = 294.270037078439, Press = 6.89801629270594 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 10 10 10 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 13000 -8420.6681 -8420.6681 -8499.5671 -8499.5671 305.34744 305.34744 23732.81 23732.81 -289.45182 -289.45182 14000 -8425.9226 -8425.9226 -8500.2041 -8500.2041 287.47684 287.47684 23738.865 23738.865 -1016.3577 -1016.3577 Loop time of 149.731 on 1 procs for 1000 steps with 2000 atoms Performance: 0.577 ns/day, 41.592 hours/ns, 6.679 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 | 149.21 | 149.21 | 149.21 | 0.0 | 99.65 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.096111 | 0.096111 | 0.096111 | 0.0 | 0.06 Output | 0.00015346 | 0.00015346 | 0.00015346 | 0.0 | 0.00 Modify | 0.36479 | 0.36479 | 0.36479 | 0.0 | 0.24 Other | | 0.06019 | | | 0.04 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119.00 ave 4119 max 4119 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: 128046.0 ave 128046 max 128046 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 128046 Ave neighs/atom = 64.023000 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 = 294.254518329013, Press = 4.68665651007902 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 10 10 10 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 14000 -8425.9226 -8425.9226 -8500.2041 -8500.2041 287.47684 287.47684 23738.865 23738.865 -1016.3577 -1016.3577 15000 -8425.7644 -8425.7644 -8501.3179 -8501.3179 292.3997 292.3997 23754.341 23754.341 -2301.0505 -2301.0505 Loop time of 151.948 on 1 procs for 1000 steps with 2000 atoms Performance: 0.569 ns/day, 42.208 hours/ns, 6.581 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 | 151.42 | 151.42 | 151.42 | 0.0 | 99.65 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.09759 | 0.09759 | 0.09759 | 0.0 | 0.06 Output | 0.00015703 | 0.00015703 | 0.00015703 | 0.0 | 0.00 Modify | 0.37537 | 0.37537 | 0.37537 | 0.0 | 0.25 Other | | 0.06026 | | | 0.04 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119.00 ave 4119 max 4119 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: 128028.0 ave 128028 max 128028 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 128028 Ave neighs/atom = 64.014000 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 = 294.191448717104, Press = 1.42799428020963 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 10 10 10 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 15000 -8425.7644 -8425.7644 -8501.3179 -8501.3179 292.3997 292.3997 23754.341 23754.341 -2301.0505 -2301.0505 16000 -8428.5517 -8428.5517 -8502.4116 -8502.4116 285.84528 285.84528 23742.459 23742.459 -1819.1268 -1819.1268 Loop time of 152.211 on 1 procs for 1000 steps with 2000 atoms Performance: 0.568 ns/day, 42.281 hours/ns, 6.570 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 | 151.68 | 151.68 | 151.68 | 0.0 | 99.65 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.099905 | 0.099905 | 0.099905 | 0.0 | 0.07 Output | 0.00015605 | 0.00015605 | 0.00015605 | 0.0 | 0.00 Modify | 0.37305 | 0.37305 | 0.37305 | 0.0 | 0.25 Other | | 0.06133 | | | 0.04 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119.00 ave 4119 max 4119 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: 128018.0 ave 128018 max 128018 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 128018 Ave neighs/atom = 64.009000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 293.886312495119, Press = -1.82622584187789 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 10 10 10 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 16000 -8428.5517 -8428.5517 -8502.4116 -8502.4116 285.84528 285.84528 23742.459 23742.459 -1819.1268 -1819.1268 17000 -8425.9913 -8425.9913 -8502.8931 -8502.8931 297.6178 297.6178 23732.039 23732.039 -999.21309 -999.21309 Loop time of 153.496 on 1 procs for 1000 steps with 2000 atoms Performance: 0.563 ns/day, 42.638 hours/ns, 6.515 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 | 152.96 | 152.96 | 152.96 | 0.0 | 99.65 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.098706 | 0.098706 | 0.098706 | 0.0 | 0.06 Output | 0.00015241 | 0.00015241 | 0.00015241 | 0.0 | 0.00 Modify | 0.37971 | 0.37971 | 0.37971 | 0.0 | 0.25 Other | | 0.06184 | | | 0.04 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119.00 ave 4119 max 4119 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: 128052.0 ave 128052 max 128052 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 128052 Ave neighs/atom = 64.026000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 293.840081847425, Press = -2.23216959942243 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 10 10 10 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 17000 -8425.9913 -8425.9913 -8502.8931 -8502.8931 297.6178 297.6178 23732.039 23732.039 -999.21309 -999.21309 18000 -8428.7138 -8428.7138 -8504.127 -8504.127 291.8569 291.8569 23705.057 23705.057 624.44245 624.44245 Loop time of 150.887 on 1 procs for 1000 steps with 2000 atoms Performance: 0.573 ns/day, 41.913 hours/ns, 6.627 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 | 150.36 | 150.36 | 150.36 | 0.0 | 99.65 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.097041 | 0.097041 | 0.097041 | 0.0 | 0.06 Output | 0.00025932 | 0.00025932 | 0.00025932 | 0.0 | 0.00 Modify | 0.37024 | 0.37024 | 0.37024 | 0.0 | 0.25 Other | | 0.06115 | | | 0.04 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119.00 ave 4119 max 4119 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: 128028.0 ave 128028 max 128028 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 128028 Ave neighs/atom = 64.014000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 293.670792922751, Press = -2.59894088882807 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 10 10 10 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 18000 -8428.7138 -8428.7138 -8504.127 -8504.127 291.8569 291.8569 23705.057 23705.057 624.44245 624.44245 19000 -8424.9051 -8424.9051 -8501.6778 -8501.6778 297.11832 297.11832 23663.843 23663.843 4101.3862 4101.3862 Loop time of 151.009 on 1 procs for 1000 steps with 2000 atoms Performance: 0.572 ns/day, 41.947 hours/ns, 6.622 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 | 150.48 | 150.48 | 150.48 | 0.0 | 99.65 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.097041 | 0.097041 | 0.097041 | 0.0 | 0.06 Output | 0.00019213 | 0.00019213 | 0.00019213 | 0.0 | 0.00 Modify | 0.37272 | 0.37272 | 0.37272 | 0.0 | 0.25 Other | | 0.06074 | | | 0.04 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119.00 ave 4119 max 4119 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: 128014.0 ave 128014 max 128014 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 128014 Ave neighs/atom = 64.007000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 293.595989286398, Press = -2.76891189501532 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 10 10 10 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 19000 -8424.9051 -8424.9051 -8501.6778 -8501.6778 297.11832 297.11832 23663.843 23663.843 4101.3862 4101.3862 20000 -8425.0755 -8425.0755 -8501.5142 -8501.5142 295.82563 295.82563 23690.577 23690.577 2244.6148 2244.6148 Loop time of 152.519 on 1 procs for 1000 steps with 2000 atoms Performance: 0.566 ns/day, 42.366 hours/ns, 6.557 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 | 151.98 | 151.98 | 151.98 | 0.0 | 99.65 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.099181 | 0.099181 | 0.099181 | 0.0 | 0.07 Output | 0.00015434 | 0.00015434 | 0.00015434 | 0.0 | 0.00 Modify | 0.37813 | 0.37813 | 0.37813 | 0.0 | 0.25 Other | | 0.06105 | | | 0.04 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119.00 ave 4119 max 4119 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: 128070.0 ave 128070 max 128070 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 128070 Ave neighs/atom = 64.035000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 293.430533899205, Press = 2.35276788366058 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 10 10 10 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 20000 -8425.0755 -8425.0755 -8501.5142 -8501.5142 295.82563 295.82563 23690.577 23690.577 2244.6148 2244.6148 21000 -8423.7617 -8423.7617 -8501.1245 -8501.1245 299.40188 299.40188 23708.752 23708.752 1130.5263 1130.5263 Loop time of 152.388 on 1 procs for 1000 steps with 2000 atoms Performance: 0.567 ns/day, 42.330 hours/ns, 6.562 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 | 151.85 | 151.85 | 151.85 | 0.0 | 99.65 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.098274 | 0.098274 | 0.098274 | 0.0 | 0.06 Output | 0.00015303 | 0.00015303 | 0.00015303 | 0.0 | 0.00 Modify | 0.37731 | 0.37731 | 0.37731 | 0.0 | 0.25 Other | | 0.06128 | | | 0.04 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119.00 ave 4119 max 4119 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: 128028.0 ave 128028 max 128028 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 128028 Ave neighs/atom = 64.014000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 293.338095758912, Press = 2.98306715717043 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 10 10 10 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 21000 -8423.7617 -8423.7617 -8501.1245 -8501.1245 299.40188 299.40188 23708.752 23708.752 1130.5263 1130.5263 22000 -8429.4523 -8429.4523 -8502.9404 -8502.9404 284.40657 284.40657 23712.05 23712.05 295.17247 295.17247 Loop time of 150.133 on 1 procs for 1000 steps with 2000 atoms Performance: 0.575 ns/day, 41.704 hours/ns, 6.661 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 | 149.6 | 149.6 | 149.6 | 0.0 | 99.65 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.098055 | 0.098055 | 0.098055 | 0.0 | 0.07 Output | 0.00015107 | 0.00015107 | 0.00015107 | 0.0 | 0.00 Modify | 0.3734 | 0.3734 | 0.3734 | 0.0 | 0.25 Other | | 0.06093 | | | 0.04 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119.00 ave 4119 max 4119 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: 128058.0 ave 128058 max 128058 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 128058 Ave neighs/atom = 64.029000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 293.280086096165, Press = 3.83142257544669 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 10 10 10 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 22000 -8429.4523 -8429.4523 -8502.9404 -8502.9404 284.40657 284.40657 23712.05 23712.05 295.17247 295.17247 23000 -8424.6513 -8424.6513 -8500.9637 -8500.9637 295.33668 295.33668 23746.089 23746.089 -1695.0024 -1695.0024 Loop time of 151.202 on 1 procs for 1000 steps with 2000 atoms Performance: 0.571 ns/day, 42.000 hours/ns, 6.614 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 | 150.66 | 150.66 | 150.66 | 0.0 | 99.64 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.097664 | 0.097664 | 0.097664 | 0.0 | 0.06 Output | 0.00015651 | 0.00015651 | 0.00015651 | 0.0 | 0.00 Modify | 0.37825 | 0.37825 | 0.37825 | 0.0 | 0.25 Other | | 0.06119 | | | 0.04 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4118.00 ave 4118 max 4118 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: 128060.0 ave 128060 max 128060 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 128060 Ave neighs/atom = 64.030000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 293.166125221264, Press = 3.58955256371124 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 10 10 10 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 23000 -8424.6513 -8424.6513 -8500.9637 -8500.9637 295.33668 295.33668 23746.089 23746.089 -1695.0024 -1695.0024 24000 -8428.0003 -8428.0003 -8502.5668 -8502.5668 288.57992 288.57992 23776.265 23776.265 -4057.2027 -4057.2027 Loop time of 152.765 on 1 procs for 1000 steps with 2000 atoms Performance: 0.566 ns/day, 42.435 hours/ns, 6.546 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 | 152.22 | 152.22 | 152.22 | 0.0 | 99.65 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.099025 | 0.099025 | 0.099025 | 0.0 | 0.06 Output | 0.00019916 | 0.00019916 | 0.00019916 | 0.0 | 0.00 Modify | 0.38048 | 0.38048 | 0.38048 | 0.0 | 0.25 Other | | 0.06124 | | | 0.04 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119.00 ave 4119 max 4119 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: 128046.0 ave 128046 max 128046 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 128046 Ave neighs/atom = 64.023000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 293.091347383311, Press = 2.54940274336024 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 10 10 10 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 24000 -8428.0003 -8428.0003 -8502.5668 -8502.5668 288.57992 288.57992 23776.265 23776.265 -4057.2027 -4057.2027 25000 -8424.3252 -8424.3252 -8499.5893 -8499.5893 291.27985 291.27985 23761.373 23761.373 -2576.6913 -2576.6913 Loop time of 150.895 on 1 procs for 1000 steps with 2000 atoms Performance: 0.573 ns/day, 41.915 hours/ns, 6.627 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 | 150.37 | 150.37 | 150.37 | 0.0 | 99.65 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.096801 | 0.096801 | 0.096801 | 0.0 | 0.06 Output | 0.00015605 | 0.00015605 | 0.00015605 | 0.0 | 0.00 Modify | 0.37093 | 0.37093 | 0.37093 | 0.0 | 0.25 Other | | 0.06054 | | | 0.04 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119.00 ave 4119 max 4119 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: 128054.0 ave 128054 max 128054 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 128054 Ave neighs/atom = 64.027000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 293.089710098702, Press = 0.451494002698893 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 10 10 10 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 25000 -8424.3252 -8424.3252 -8499.5893 -8499.5893 291.27985 291.27985 23761.373 23761.373 -2576.6913 -2576.6913 26000 -8426.4251 -8426.4251 -8502.1861 -8502.1861 293.20316 293.20316 23728.64 23728.64 -646.27095 -646.27095 Loop time of 151.925 on 1 procs for 1000 steps with 2000 atoms Performance: 0.569 ns/day, 42.201 hours/ns, 6.582 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 | 151.39 | 151.39 | 151.39 | 0.0 | 99.65 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.098053 | 0.098053 | 0.098053 | 0.0 | 0.06 Output | 0.00015571 | 0.00015571 | 0.00015571 | 0.0 | 0.00 Modify | 0.37399 | 0.37399 | 0.37399 | 0.0 | 0.25 Other | | 0.06119 | | | 0.04 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119.00 ave 4119 max 4119 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: 128038.0 ave 128038 max 128038 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 128038 Ave neighs/atom = 64.019000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 293.075629594633, Press = 0.28460594667375 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 10 10 10 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 26000 -8426.4251 -8426.4251 -8502.1861 -8502.1861 293.20316 293.20316 23728.64 23728.64 -646.27095 -646.27095 27000 -8421.5917 -8421.5917 -8498.5399 -8498.5399 297.7977 297.7977 23726.303 23726.303 304.35455 304.35455 Loop time of 153.542 on 1 procs for 1000 steps with 2000 atoms Performance: 0.563 ns/day, 42.651 hours/ns, 6.513 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 | 152.99 | 152.99 | 152.99 | 0.0 | 99.64 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.099239 | 0.099239 | 0.099239 | 0.0 | 0.06 Output | 0.00015623 | 0.00015623 | 0.00015623 | 0.0 | 0.00 Modify | 0.38627 | 0.38627 | 0.38627 | 0.0 | 0.25 Other | | 0.06218 | | | 0.04 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4115.00 ave 4115 max 4115 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: 128036.0 ave 128036 max 128036 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 128036 Ave neighs/atom = 64.018000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 293.172260495641, Press = 0.412687830213166 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 10 10 10 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.755 | 3.755 | 3.755 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 27000 -8421.5917 -8421.5917 -8498.5399 -8498.5399 297.7977 297.7977 23726.303 23726.303 304.35455 304.35455 28000 -8424.7254 -8424.7254 -8502.1317 -8502.1317 299.57021 299.57021 23705.013 23705.013 1096.9819 1096.9819 Loop time of 151.647 on 1 procs for 1000 steps with 2000 atoms Performance: 0.570 ns/day, 42.124 hours/ns, 6.594 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 | 151.11 | 151.11 | 151.11 | 0.0 | 99.65 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.097826 | 0.097826 | 0.097826 | 0.0 | 0.06 Output | 0.00015298 | 0.00015298 | 0.00015298 | 0.0 | 0.00 Modify | 0.37563 | 0.37563 | 0.37563 | 0.0 | 0.25 Other | | 0.0607 | | | 0.04 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119.00 ave 4119 max 4119 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: 128048.0 ave 128048 max 128048 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 128048 Ave neighs/atom = 64.024000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" jump SELF break # Write final averaged volume to file if temperature and volume have converged; otherwise wirte a # flag to indicate non-convergence. variable myStep equal step if "${myStep} < 2000000" then "print '${V}' file output/vol_T293.15.out" else "print 'not_converged' file output/vol_T293.15.out" print '${V}' file output/vol_T293.15.out 23722.1298916476 print "LAMMPS calculation completed" LAMMPS calculation completed quit 0