# 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.863657310605049*${_u_distance} variable latticeconst_converted equal 2.863657310605049*1 lattice bcc ${latticeconst_converted} lattice bcc 2.86365731060505 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.001 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_LeeShimPark_2001_FeCr__MO_150993986463_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 23483.5168277961 variable V0_metal equal ${V0}/(${_u_distance}*${_u_distance}*${_u_distance}) variable V0_metal equal 23483.5168277961/(${_u_distance}*${_u_distance}*${_u_distance}) variable V0_metal equal 23483.5168277961/(1*${_u_distance}*${_u_distance}) variable V0_metal equal 23483.5168277961/(1*1*${_u_distance}) variable V0_metal equal 23483.5168277961/(1*1*1) variable V0_metal_times1000 equal ${V0_metal}*1000 variable V0_metal_times1000 equal 23483.5168277961*1000 print "Initial system volume: ${V0_metal} Angstroms^3" Initial system volume: 23483.5168277961 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_150993986463_000#item-citation CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 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.2866 3445.2866 1000 -8425.0239 -8425.0239 -8497.9604 -8497.9604 282.27172 282.27172 23692.276 23692.276 2408.0418 2408.0418 Loop time of 53.0315 on 1 procs for 1000 steps with 2000 atoms Performance: 1.629 ns/day, 14.731 hours/ns, 18.857 timesteps/s 99.2% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 52.796 | 52.796 | 52.796 | 0.0 | 99.56 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.036937 | 0.036937 | 0.036937 | 0.0 | 0.07 Output | 0.00018643 | 0.00018643 | 0.00018643 | 0.0 | 0.00 Modify | 0.16948 | 0.16948 | 0.16948 | 0.0 | 0.32 Other | | 0.02864 | | | 0.05 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.276 23692.276 2408.0418 2408.0418 2000 -8425.4004 -8425.4004 -8500.2416 -8500.2416 289.64312 289.64312 23725.037 23725.037 -178.5232 -178.5232 Loop time of 53.8883 on 1 procs for 1000 steps with 2000 atoms Performance: 1.603 ns/day, 14.969 hours/ns, 18.557 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 | 53.671 | 53.671 | 53.671 | 0.0 | 99.60 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.035796 | 0.035796 | 0.035796 | 0.0 | 0.07 Output | 0.00011402 | 0.00011402 | 0.00011402 | 0.0 | 0.00 Modify | 0.15483 | 0.15483 | 0.15483 | 0.0 | 0.29 Other | | 0.02704 | | | 0.05 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.5232 -178.5232 3000 -8426.809 -8426.809 -8497.7632 -8497.7632 274.60014 274.60014 23739.586 23739.586 -1024.9432 -1024.9432 Loop time of 54.0347 on 1 procs for 1000 steps with 2000 atoms Performance: 1.599 ns/day, 15.010 hours/ns, 18.507 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 | 53.821 | 53.821 | 53.821 | 0.0 | 99.60 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.035204 | 0.035204 | 0.035204 | 0.0 | 0.07 Output | 0.00011201 | 0.00011201 | 0.00011201 | 0.0 | 0.00 Modify | 0.15258 | 0.15258 | 0.15258 | 0.0 | 0.28 Other | | 0.02577 | | | 0.05 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.9432 -1024.9432 4000 -8424.0945 -8424.0945 -8502.4446 -8502.4446 303.22271 303.22271 23759.323 23759.323 -2785.8961 -2785.8961 Loop time of 54.3091 on 1 procs for 1000 steps with 2000 atoms Performance: 1.591 ns/day, 15.086 hours/ns, 18.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 | 54.09 | 54.09 | 54.09 | 0.0 | 99.60 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.036109 | 0.036109 | 0.036109 | 0.0 | 0.07 Output | 0.00012083 | 0.00012083 | 0.00012083 | 0.0 | 0.00 Modify | 0.156 | 0.156 | 0.156 | 0.0 | 0.29 Other | | 0.02702 | | | 0.05 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.22271 303.22271 23759.323 23759.323 -2785.8961 -2785.8961 5000 -8427.0408 -8427.0408 -8501.1089 -8501.1089 286.65155 286.65155 23754.254 23754.254 -2415.6515 -2415.6515 Loop time of 54.0082 on 1 procs for 1000 steps with 2000 atoms Performance: 1.600 ns/day, 15.002 hours/ns, 18.516 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 | 53.792 | 53.792 | 53.792 | 0.0 | 99.60 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.035193 | 0.035193 | 0.035193 | 0.0 | 0.07 Output | 5.9984e-05 | 5.9984e-05 | 5.9984e-05 | 0.0 | 0.00 Modify | 0.15378 | 0.15378 | 0.15378 | 0.0 | 0.28 Other | | 0.02677 | | | 0.05 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.418488573909, Press = -145.778488308888 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.6515 -2415.6515 6000 -8425.4554 -8425.4554 -8501.9076 -8501.9076 295.87768 295.87768 23740.732 23740.732 -1536.6526 -1536.6526 Loop time of 53.862 on 1 procs for 1000 steps with 2000 atoms Performance: 1.604 ns/day, 14.962 hours/ns, 18.566 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 | 53.638 | 53.638 | 53.638 | 0.0 | 99.58 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.03545 | 0.03545 | 0.03545 | 0.0 | 0.07 Output | 8.4142e-05 | 8.4142e-05 | 8.4142e-05 | 0.0 | 0.00 Modify | 0.16128 | 0.16128 | 0.16128 | 0.0 | 0.30 Other | | 0.02734 | | | 0.05 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.01260216549, Press = -19.8351688883026 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.87768 295.87768 23740.732 23740.732 -1536.6526 -1536.6526 7000 -8422.887 -8422.887 -8498.9469 -8498.9469 294.35987 294.35987 23721.445 23721.445 286.42315 286.42315 Loop time of 53.883 on 1 procs for 1000 steps with 2000 atoms Performance: 1.603 ns/day, 14.968 hours/ns, 18.559 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 | 53.658 | 53.658 | 53.658 | 0.0 | 99.58 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.035575 | 0.035575 | 0.035575 | 0.0 | 0.07 Output | 5.5685e-05 | 5.5685e-05 | 5.5685e-05 | 0.0 | 0.00 Modify | 0.16137 | 0.16137 | 0.16137 | 0.0 | 0.30 Other | | 0.02767 | | | 0.05 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.820591321661, Press = -14.660988407761 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.887 -8422.887 -8498.9469 -8498.9469 294.35987 294.35987 23721.445 23721.445 286.42315 286.42315 8000 -8426.7698 -8426.7698 -8502.4698 -8502.4698 292.96648 292.96648 23701.017 23701.017 1239.758 1239.758 Loop time of 54.0068 on 1 procs for 1000 steps with 2000 atoms Performance: 1.600 ns/day, 15.002 hours/ns, 18.516 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 | 53.784 | 53.784 | 53.784 | 0.0 | 99.59 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.035176 | 0.035176 | 0.035176 | 0.0 | 0.07 Output | 0.000101 | 0.000101 | 0.000101 | 0.0 | 0.00 Modify | 0.15998 | 0.15998 | 0.15998 | 0.0 | 0.30 Other | | 0.02717 | | | 0.05 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.160436262301, Press = -13.1960206594905 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.7698 -8426.7698 -8502.4698 -8502.4698 292.96648 292.96648 23701.017 23701.017 1239.758 1239.758 9000 -8423.0518 -8423.0518 -8501.8866 -8501.8866 305.09876 305.09876 23672.344 23672.344 3499.2638 3499.2638 Loop time of 54.1436 on 1 procs for 1000 steps with 2000 atoms Performance: 1.596 ns/day, 15.040 hours/ns, 18.469 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 | 53.916 | 53.916 | 53.916 | 0.0 | 99.58 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.036504 | 0.036504 | 0.036504 | 0.0 | 0.07 Output | 8.4737e-05 | 8.4737e-05 | 8.4737e-05 | 0.0 | 0.00 Modify | 0.162 | 0.162 | 0.162 | 0.0 | 0.30 Other | | 0.02856 | | | 0.05 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.542530045409, Press = -3.92103096826358 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.0518 -8423.0518 -8501.8866 -8501.8866 305.09876 305.09876 23672.344 23672.344 3499.2638 3499.2638 10000 -8426.7252 -8426.7252 -8501.6701 -8501.6701 290.04449 290.04449 23682.882 23682.882 2586.9922 2586.9922 Loop time of 54.2666 on 1 procs for 1000 steps with 2000 atoms Performance: 1.592 ns/day, 15.074 hours/ns, 18.428 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 | 54.037 | 54.037 | 54.037 | 0.0 | 99.58 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.036489 | 0.036489 | 0.036489 | 0.0 | 0.07 Output | 0.00018418 | 0.00018418 | 0.00018418 | 0.0 | 0.00 Modify | 0.16613 | 0.16613 | 0.16613 | 0.0 | 0.31 Other | | 0.027 | | | 0.05 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.816279220524, Press = 7.53014383030853 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.7252 -8426.7252 -8501.6701 -8501.6701 290.04449 290.04449 23682.882 23682.882 2586.9922 2586.9922 11000 -8422.0423 -8422.0423 -8499.8227 -8499.8227 301.01822 301.01822 23710.964 23710.964 1134.5101 1134.5101 Loop time of 54.1422 on 1 procs for 1000 steps with 2000 atoms Performance: 1.596 ns/day, 15.039 hours/ns, 18.470 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 | 53.915 | 53.915 | 53.915 | 0.0 | 99.58 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.036381 | 0.036381 | 0.036381 | 0.0 | 0.07 Output | 9.889e-05 | 9.889e-05 | 9.889e-05 | 0.0 | 0.00 Modify | 0.16424 | 0.16424 | 0.16424 | 0.0 | 0.30 Other | | 0.02613 | | | 0.05 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.029500215212, Press = 6.74871515349487 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.0423 -8422.0423 -8499.8227 -8499.8227 301.01822 301.01822 23710.964 23710.964 1134.5101 1134.5101 12000 -8426.1649 -8426.1649 -8500.4486 -8500.4486 287.48547 287.48547 23722.22 23722.22 -91.301447 -91.301447 Loop time of 53.9414 on 1 procs for 1000 steps with 2000 atoms Performance: 1.602 ns/day, 14.984 hours/ns, 18.539 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 | 53.719 | 53.719 | 53.719 | 0.0 | 99.59 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.035516 | 0.035516 | 0.035516 | 0.0 | 0.07 Output | 7.4502e-05 | 7.4502e-05 | 7.4502e-05 | 0.0 | 0.00 Modify | 0.16144 | 0.16144 | 0.16144 | 0.0 | 0.30 Other | | 0.02546 | | | 0.05 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.131637882459, Press = 9.46566573801482 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.1649 -8426.1649 -8500.4486 -8500.4486 287.48547 287.48547 23722.22 23722.22 -91.301447 -91.301447 13000 -8420.6507 -8420.6507 -8499.5574 -8499.5574 305.37704 305.37704 23732.837 23732.837 -289.23817 -289.23817 Loop time of 53.8866 on 1 procs for 1000 steps with 2000 atoms Performance: 1.603 ns/day, 14.969 hours/ns, 18.557 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 | 53.664 | 53.664 | 53.664 | 0.0 | 99.59 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.035473 | 0.035473 | 0.035473 | 0.0 | 0.07 Output | 0.00010274 | 0.00010274 | 0.00010274 | 0.0 | 0.00 Modify | 0.16162 | 0.16162 | 0.16162 | 0.0 | 0.30 Other | | 0.02536 | | | 0.05 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.269640444561, Press = 6.90228858891914 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.6507 -8420.6507 -8499.5574 -8499.5574 305.37704 305.37704 23732.837 23732.837 -289.23817 -289.23817 14000 -8425.929 -8425.929 -8500.2078 -8500.2078 287.46683 287.46683 23738.818 23738.818 -1013.7762 -1013.7762 Loop time of 53.7771 on 1 procs for 1000 steps with 2000 atoms Performance: 1.607 ns/day, 14.938 hours/ns, 18.595 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 | 53.554 | 53.554 | 53.554 | 0.0 | 99.59 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.03549 | 0.03549 | 0.03549 | 0.0 | 0.07 Output | 9.1607e-05 | 9.1607e-05 | 9.1607e-05 | 0.0 | 0.00 Modify | 0.16189 | 0.16189 | 0.16189 | 0.0 | 0.30 Other | | 0.02553 | | | 0.05 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.254650679389, Press = 4.67993015460134 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.929 -8425.929 -8500.2078 -8500.2078 287.46683 287.46683 23738.818 23738.818 -1013.7762 -1013.7762 15000 -8425.7448 -8425.7448 -8501.3077 -8501.3077 292.43621 292.43621 23754.326 23754.326 -2297.7484 -2297.7484 Loop time of 53.5848 on 1 procs for 1000 steps with 2000 atoms Performance: 1.612 ns/day, 14.885 hours/ns, 18.662 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 | 53.365 | 53.365 | 53.365 | 0.0 | 99.59 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.034714 | 0.034714 | 0.034714 | 0.0 | 0.06 Output | 8.7233e-05 | 8.7233e-05 | 8.7233e-05 | 0.0 | 0.00 Modify | 0.15991 | 0.15991 | 0.15991 | 0.0 | 0.30 Other | | 0.02554 | | | 0.05 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.191642102872, Press = 1.41229783868615 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.7448 -8425.7448 -8501.3077 -8501.3077 292.43621 292.43621 23754.326 23754.326 -2297.7484 -2297.7484 16000 -8428.5226 -8428.5226 -8502.4001 -8502.4001 285.91346 285.91346 23742.401 23742.401 -1812.0266 -1812.0266 Loop time of 53.5009 on 1 procs for 1000 steps with 2000 atoms Performance: 1.615 ns/day, 14.861 hours/ns, 18.691 timesteps/s 100.0% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 53.282 | 53.282 | 53.282 | 0.0 | 99.59 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.034687 | 0.034687 | 0.034687 | 0.0 | 0.06 Output | 6.3334e-05 | 6.3334e-05 | 6.3334e-05 | 0.0 | 0.00 Modify | 0.1593 | 0.1593 | 0.1593 | 0.0 | 0.30 Other | | 0.02534 | | | 0.05 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: 128020.0 ave 128020 max 128020 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 128020 Ave neighs/atom = 64.010000 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.889433627271, Press = -1.84168680911592 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.5226 -8428.5226 -8502.4001 -8502.4001 285.91346 285.91346 23742.401 23742.401 -1812.0266 -1812.0266 17000 -8425.9191 -8425.9191 -8502.8484 -8502.8484 297.7243 297.7243 23732.167 23732.167 -999.57226 -999.57226 Loop time of 53.6281 on 1 procs for 1000 steps with 2000 atoms Performance: 1.611 ns/day, 14.897 hours/ns, 18.647 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 | 53.408 | 53.408 | 53.408 | 0.0 | 99.59 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.034504 | 0.034504 | 0.034504 | 0.0 | 0.06 Output | 8.8381e-05 | 8.8381e-05 | 8.8381e-05 | 0.0 | 0.00 Modify | 0.16028 | 0.16028 | 0.16028 | 0.0 | 0.30 Other | | 0.02569 | | | 0.05 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.839042800818, Press = -2.29796908201378 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.9191 -8425.9191 -8502.8484 -8502.8484 297.7243 297.7243 23732.167 23732.167 -999.57226 -999.57226 18000 -8423.7666 -8423.7666 -8501.4594 -8501.4594 300.67939 300.67939 23711.28 23711.28 766.23559 766.23559 Loop time of 53.6428 on 1 procs for 1000 steps with 2000 atoms Performance: 1.611 ns/day, 14.901 hours/ns, 18.642 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 | 53.423 | 53.423 | 53.423 | 0.0 | 99.59 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.034687 | 0.034687 | 0.034687 | 0.0 | 0.06 Output | 6.6356e-05 | 6.6356e-05 | 6.6356e-05 | 0.0 | 0.00 Modify | 0.15888 | 0.15888 | 0.15888 | 0.0 | 0.30 Other | | 0.02586 | | | 0.05 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: 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.644093532862, Press = -2.453725921639 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 -8423.7666 -8423.7666 -8501.4594 -8501.4594 300.67939 300.67939 23711.28 23711.28 766.23559 766.23559 19000 -8428.3826 -8428.3826 -8503.553 -8503.553 290.91711 290.91711 23653.776 23653.776 4411.425 4411.425 Loop time of 53.5862 on 1 procs for 1000 steps with 2000 atoms Performance: 1.612 ns/day, 14.885 hours/ns, 18.662 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 | 53.367 | 53.367 | 53.367 | 0.0 | 99.59 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.034653 | 0.034653 | 0.034653 | 0.0 | 0.06 Output | 0.00020332 | 0.00020332 | 0.00020332 | 0.0 | 0.00 Modify | 0.15865 | 0.15865 | 0.15865 | 0.0 | 0.30 Other | | 0.02554 | | | 0.05 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.48714346137, Press = -1.43363384602608 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 -8428.3826 -8428.3826 -8503.553 -8503.553 290.91711 290.91711 23653.776 23653.776 4411.425 4411.425 20000 -8425.6073 -8425.6073 -8501.8177 -8501.8177 294.94207 294.94207 23685.177 23685.177 2569.1585 2569.1585 Loop time of 53.647 on 1 procs for 1000 steps with 2000 atoms Performance: 1.611 ns/day, 14.902 hours/ns, 18.640 timesteps/s 100.0% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 53.421 | 53.421 | 53.421 | 0.0 | 99.58 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.035171 | 0.035171 | 0.035171 | 0.0 | 0.07 Output | 9.4964e-05 | 9.4964e-05 | 9.4964e-05 | 0.0 | 0.00 Modify | 0.16198 | 0.16198 | 0.16198 | 0.0 | 0.30 Other | | 0.0287 | | | 0.05 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.286657910698, Press = 2.74225445249418 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.6073 -8425.6073 -8501.8177 -8501.8177 294.94207 294.94207 23685.177 23685.177 2569.1585 2569.1585 21000 -8427.2433 -8427.2433 -8502.9985 -8502.9985 293.1802 293.1802 23704.618 23704.618 1012.8165 1012.8165 Loop time of 53.5676 on 1 procs for 1000 steps with 2000 atoms Performance: 1.613 ns/day, 14.880 hours/ns, 18.668 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 | 53.348 | 53.348 | 53.348 | 0.0 | 99.59 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.034781 | 0.034781 | 0.034781 | 0.0 | 0.06 Output | 9.5301e-05 | 9.5301e-05 | 9.5301e-05 | 0.0 | 0.00 Modify | 0.15903 | 0.15903 | 0.15903 | 0.0 | 0.30 Other | | 0.02562 | | | 0.05 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: 128032.0 ave 128032 max 128032 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 128032 Ave neighs/atom = 64.016000 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.252094027662, Press = 3.2931683159916 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 -8427.2433 -8427.2433 -8502.9985 -8502.9985 293.1802 293.1802 23704.618 23704.618 1012.8165 1012.8165 22000 -8423.2448 -8423.2448 -8499.5668 -8499.5668 295.37387 295.37387 23727.561 23727.561 -69.467639 -69.467639 Loop time of 53.3008 on 1 procs for 1000 steps with 2000 atoms Performance: 1.621 ns/day, 14.806 hours/ns, 18.761 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 | 53.082 | 53.082 | 53.082 | 0.0 | 99.59 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.03487 | 0.03487 | 0.03487 | 0.0 | 0.07 Output | 8.9387e-05 | 8.9387e-05 | 8.9387e-05 | 0.0 | 0.00 Modify | 0.15872 | 0.15872 | 0.15872 | 0.0 | 0.30 Other | | 0.02555 | | | 0.05 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.121521604355, Press = 3.49937117422787 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 -8423.2448 -8423.2448 -8499.5668 -8499.5668 295.37387 295.37387 23727.561 23727.561 -69.467639 -69.467639 23000 -8427.9157 -8427.9157 -8502.6887 -8502.6887 289.37928 289.37928 23744.876 23744.876 -1990.7042 -1990.7042 Loop time of 53.4358 on 1 procs for 1000 steps with 2000 atoms Performance: 1.617 ns/day, 14.843 hours/ns, 18.714 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 | 53.221 | 53.221 | 53.221 | 0.0 | 99.60 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.034323 | 0.034323 | 0.034323 | 0.0 | 0.06 Output | 9.0971e-05 | 9.0971e-05 | 9.0971e-05 | 0.0 | 0.00 Modify | 0.15632 | 0.15632 | 0.15632 | 0.0 | 0.29 Other | | 0.02413 | | | 0.05 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: 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 = 293.059355352867, Press = 2.81092017476944 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 -8427.9157 -8427.9157 -8502.6887 -8502.6887 289.37928 289.37928 23744.876 23744.876 -1990.7042 -1990.7042 24000 -8423.5413 -8423.5413 -8500.1442 -8500.1442 296.4612 296.4612 23783.837 23783.837 -4056.8484 -4056.8484 Loop time of 53.4368 on 1 procs for 1000 steps with 2000 atoms Performance: 1.617 ns/day, 14.844 hours/ns, 18.714 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 | 53.224 | 53.224 | 53.224 | 0.0 | 99.60 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.034131 | 0.034131 | 0.034131 | 0.0 | 0.06 Output | 8.7383e-05 | 8.7383e-05 | 8.7383e-05 | 0.0 | 0.00 Modify | 0.15486 | 0.15486 | 0.15486 | 0.0 | 0.29 Other | | 0.02394 | | | 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 = 292.966659483934, Press = 0.605730911960306 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 -8423.5413 -8423.5413 -8500.1442 -8500.1442 296.4612 296.4612 23783.837 23783.837 -4056.8484 -4056.8484 25000 -8423.8467 -8423.8467 -8499.1488 -8499.1488 291.42661 291.42661 23759.436 23759.436 -2368.7252 -2368.7252 Loop time of 53.5415 on 1 procs for 1000 steps with 2000 atoms Performance: 1.614 ns/day, 14.873 hours/ns, 18.677 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 | 53.328 | 53.328 | 53.328 | 0.0 | 99.60 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.034277 | 0.034277 | 0.034277 | 0.0 | 0.06 Output | 6.2282e-05 | 6.2282e-05 | 6.2282e-05 | 0.0 | 0.00 Modify | 0.15559 | 0.15559 | 0.15559 | 0.0 | 0.29 Other | | 0.02398 | | | 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: 128050.0 ave 128050 max 128050 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 128050 Ave neighs/atom = 64.025000 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.974839743572, Press = -0.950442221894222 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 -8423.8467 -8423.8467 -8499.1488 -8499.1488 291.42661 291.42661 23759.436 23759.436 -2368.7252 -2368.7252 26000 -8425.0585 -8425.0585 -8501.5102 -8501.5102 295.87597 295.87597 23725.336 23725.336 -262.42813 -262.42813 Loop time of 53.4009 on 1 procs for 1000 steps with 2000 atoms Performance: 1.618 ns/day, 14.834 hours/ns, 18.726 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 | 53.183 | 53.183 | 53.183 | 0.0 | 99.59 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.034452 | 0.034452 | 0.034452 | 0.0 | 0.06 Output | 8.6219e-05 | 8.6219e-05 | 8.6219e-05 | 0.0 | 0.00 Modify | 0.15817 | 0.15817 | 0.15817 | 0.0 | 0.30 Other | | 0.02472 | | | 0.05 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: 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" 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 23721.9238688783 print "LAMMPS calculation completed" LAMMPS calculation completed quit 0