# periodic boundary conditions along all three dimensions boundary p p p # Set neighbor skin variable neigh_skin equal 2.0*${_u_distance} variable neigh_skin equal 2.0*1 neighbor ${neigh_skin} bin neighbor 2 bin # create a supercell with cubic lattice (fcc, bcc, sc, or diamond) # using 10*10*10 conventional (orthogonal) unit cells variable latticeconst_converted equal 5.428877472877502*${_u_distance} variable latticeconst_converted equal 5.428877472877502*1 lattice diamond ${latticeconst_converted} lattice diamond 5.4288774728775 Lattice spacing in x,y,z = 5.4288775 5.4288775 5.4288775 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 (54.288775 54.288775 54.288775) 1 by 1 by 1 MPI processor grid create_atoms 1 box Created 8000 atoms using lattice units in orthogonal box = (0.0000000 0.0000000 0.0000000) to (54.288775 54.288775 54.288775) create_atoms CPU = 0.003 seconds variable mass_converted equal 28.0855*${_u_mass} variable mass_converted equal 28.0855*1 kim_interactions Si WARNING: 'kim_' has been renamed to 'kim '. Please update your input. kim interactions Si #=== BEGIN kim interactions ================================== pair_style kim Tersoff_LAMMPS_ErhartAlbe_2005_SiC__MO_903987585848_005 pair_coeff * * Si #=== END kim interactions ==================================== mass 1 ${mass_converted} mass 1 28.0855 # initial volume variable v equal vol # assign formula variable V0 equal ${v} # evaluate initial value variable V0 equal 160003.734725217 variable V0_metal equal ${V0}/(${_u_distance}*${_u_distance}*${_u_distance}) variable V0_metal equal 160003.734725217/(${_u_distance}*${_u_distance}*${_u_distance}) variable V0_metal equal 160003.734725217/(1*${_u_distance}*${_u_distance}) variable V0_metal equal 160003.734725217/(1*1*${_u_distance}) variable V0_metal equal 160003.734725217/(1*1*1) variable V0_metal_times1000 equal ${V0_metal}*1000 variable V0_metal_times1000 equal 160003.734725217*1000 print "Initial system volume: ${V0_metal} Angstroms^3" Initial system volume: 160003.734725217 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_903987585848_005#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 = 4.96 ghost atom cutoff = 4.96 binsize = 2.48, bins = 22 22 22 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, ghost, cut 4.96 pair build: full/bin/ghost stencil: full/ghost/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.617 | 4.617 | 4.617 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 0 -36724.214 -36724.214 -37027.317 -37027.317 293.15 293.15 160003.73 160003.73 2023.3843 2023.3843 1000 -36393.166 -36393.166 -36708.063 -36708.063 304.55633 304.55633 160801.53 160801.53 872.88706 872.88706 Loop time of 54.3365 on 1 procs for 1000 steps with 8000 atoms Performance: 1.590 ns/day, 15.093 hours/ns, 18.404 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.625 | 53.625 | 53.625 | 0.0 | 98.69 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.04668 | 0.04668 | 0.04668 | 0.0 | 0.09 Output | 0.00027889 | 0.00027889 | 0.00027889 | 0.0 | 0.00 Modify | 0.5952 | 0.5952 | 0.5952 | 0.0 | 1.10 Other | | 0.069 | | | 0.13 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4995.00 ave 4995 max 4995 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: 224000.0 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 28.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 = 4.96 ghost atom cutoff = 4.96 binsize = 2.48, bins = 22 22 22 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, ghost, cut 4.96 pair build: full/bin/ghost stencil: full/ghost/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.617 | 4.617 | 4.617 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 1000 -36393.166 -36393.166 -36708.063 -36708.063 304.55633 304.55633 160801.53 160801.53 872.88706 872.88706 2000 -36425.065 -36425.065 -36728.911 -36728.911 293.86879 293.86879 160809.79 160809.79 654.2549 654.2549 Loop time of 55.059 on 1 procs for 1000 steps with 8000 atoms Performance: 1.569 ns/day, 15.294 hours/ns, 18.162 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.349 | 54.349 | 54.349 | 0.0 | 98.71 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.046916 | 0.046916 | 0.046916 | 0.0 | 0.09 Output | 0.00017668 | 0.00017668 | 0.00017668 | 0.0 | 0.00 Modify | 0.59754 | 0.59754 | 0.59754 | 0.0 | 1.09 Other | | 0.06576 | | | 0.12 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4995.00 ave 4995 max 4995 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: 224000.0 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 28.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 = 4.96 ghost atom cutoff = 4.96 binsize = 2.48, bins = 22 22 22 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, ghost, cut 4.96 pair build: full/bin/ghost stencil: full/ghost/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.617 | 4.617 | 4.617 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 2000 -36425.065 -36425.065 -36728.911 -36728.911 293.86879 293.86879 160809.79 160809.79 654.2549 654.2549 3000 -36411.191 -36411.191 -36715.56 -36715.56 294.37505 294.37505 160847.83 160847.83 524.65429 524.65429 Loop time of 54.9379 on 1 procs for 1000 steps with 8000 atoms Performance: 1.573 ns/day, 15.261 hours/ns, 18.202 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.221 | 54.221 | 54.221 | 0.0 | 98.69 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.047308 | 0.047308 | 0.047308 | 0.0 | 0.09 Output | 0.00015264 | 0.00015264 | 0.00015264 | 0.0 | 0.00 Modify | 0.60144 | 0.60144 | 0.60144 | 0.0 | 1.09 Other | | 0.06836 | | | 0.12 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4995.00 ave 4995 max 4995 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: 224000.0 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 28.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 = 4.96 ghost atom cutoff = 4.96 binsize = 2.48, bins = 22 22 22 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, ghost, cut 4.96 pair build: full/bin/ghost stencil: full/ghost/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.617 | 4.617 | 4.617 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 3000 -36411.191 -36411.191 -36715.56 -36715.56 294.37505 294.37505 160847.83 160847.83 524.65429 524.65429 4000 -36418.398 -36418.398 -36723.408 -36723.408 294.99452 294.99452 160871.62 160871.62 194.5681 194.5681 Loop time of 56.5019 on 1 procs for 1000 steps with 8000 atoms Performance: 1.529 ns/day, 15.695 hours/ns, 17.699 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 | 55.757 | 55.757 | 55.757 | 0.0 | 98.68 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.048605 | 0.048605 | 0.048605 | 0.0 | 0.09 Output | 0.00015967 | 0.00015967 | 0.00015967 | 0.0 | 0.00 Modify | 0.62434 | 0.62434 | 0.62434 | 0.0 | 1.10 Other | | 0.07223 | | | 0.13 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4995.00 ave 4995 max 4995 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: 224000.0 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 28.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 = 4.96 ghost atom cutoff = 4.96 binsize = 2.48, bins = 22 22 22 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, ghost, cut 4.96 pair build: full/bin/ghost stencil: full/ghost/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.617 | 4.617 | 4.617 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 4000 -36418.398 -36418.398 -36723.408 -36723.408 294.99452 294.99452 160871.62 160871.62 194.5681 194.5681 5000 -36417.434 -36417.434 -36710.479 -36710.479 283.42229 283.42229 160913.75 160913.75 39.600241 39.600241 Loop time of 53.9595 on 1 procs for 1000 steps with 8000 atoms Performance: 1.601 ns/day, 14.989 hours/ns, 18.532 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.259 | 53.259 | 53.259 | 0.0 | 98.70 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.046774 | 0.046774 | 0.046774 | 0.0 | 0.09 Output | 0.00013797 | 0.00013797 | 0.00013797 | 0.0 | 0.00 Modify | 0.58876 | 0.58876 | 0.58876 | 0.0 | 1.09 Other | | 0.06528 | | | 0.12 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4995.00 ave 4995 max 4995 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: 224000.0 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 28.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 = 294.335221908485, Press = 184.979322852748 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 = 4.96 ghost atom cutoff = 4.96 binsize = 2.48, bins = 22 22 22 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, ghost, cut 4.96 pair build: full/bin/ghost stencil: full/ghost/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.617 | 4.617 | 4.617 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 5000 -36417.434 -36417.434 -36710.479 -36710.479 283.42229 283.42229 160913.75 160913.75 39.600241 39.600241 6000 -36412.256 -36412.256 -36710.258 -36710.258 288.2168 288.2168 160914.35 160914.35 -110.08429 -110.08429 Loop time of 52.3767 on 1 procs for 1000 steps with 8000 atoms Performance: 1.650 ns/day, 14.549 hours/ns, 19.092 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 | 51.695 | 51.695 | 51.695 | 0.0 | 98.70 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.045545 | 0.045545 | 0.045545 | 0.0 | 0.09 Output | 0.0001296 | 0.0001296 | 0.0001296 | 0.0 | 0.00 Modify | 0.57533 | 0.57533 | 0.57533 | 0.0 | 1.10 Other | | 0.06039 | | | 0.12 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4995.00 ave 4995 max 4995 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: 224000.0 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 28.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 = 292.79911457405, Press = 11.764647521104 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 = 4.96 ghost atom cutoff = 4.96 binsize = 2.48, bins = 22 22 22 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, ghost, cut 4.96 pair build: full/bin/ghost stencil: full/ghost/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.617 | 4.617 | 4.617 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 6000 -36412.256 -36412.256 -36710.258 -36710.258 288.2168 288.2168 160914.35 160914.35 -110.08429 -110.08429 7000 -36419.131 -36419.131 -36726.383 -36726.383 297.16321 297.16321 160894.66 160894.66 64.316425 64.316425 Loop time of 55.3476 on 1 procs for 1000 steps with 8000 atoms Performance: 1.561 ns/day, 15.374 hours/ns, 18.068 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.612 | 54.612 | 54.612 | 0.0 | 98.67 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.047174 | 0.047174 | 0.047174 | 0.0 | 0.09 Output | 9.8475e-05 | 9.8475e-05 | 9.8475e-05 | 0.0 | 0.00 Modify | 0.61926 | 0.61926 | 0.61926 | 0.0 | 1.12 Other | | 0.0695 | | | 0.13 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4995.00 ave 4995 max 4995 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: 224000.0 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 28.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 = 293.120946823687, Press = 6.84639136669055 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 = 4.96 ghost atom cutoff = 4.96 binsize = 2.48, bins = 22 22 22 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, ghost, cut 4.96 pair build: full/bin/ghost stencil: full/ghost/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.617 | 4.617 | 4.617 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 7000 -36419.131 -36419.131 -36726.383 -36726.383 297.16321 297.16321 160894.66 160894.66 64.316425 64.316425 8000 -36411.315 -36411.315 -36710.937 -36710.937 289.78279 289.78279 160948.55 160948.55 -353.59902 -353.59902 Loop time of 54.8802 on 1 procs for 1000 steps with 8000 atoms Performance: 1.574 ns/day, 15.245 hours/ns, 18.222 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.158 | 54.158 | 54.158 | 0.0 | 98.68 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.046503 | 0.046503 | 0.046503 | 0.0 | 0.08 Output | 0.00013333 | 0.00013333 | 0.00013333 | 0.0 | 0.00 Modify | 0.6096 | 0.6096 | 0.6096 | 0.0 | 1.11 Other | | 0.06636 | | | 0.12 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4995.00 ave 4995 max 4995 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: 224000.0 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 28.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 = 293.062813179363, Press = 2.7672411739741 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 = 4.96 ghost atom cutoff = 4.96 binsize = 2.48, bins = 22 22 22 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, ghost, cut 4.96 pair build: full/bin/ghost stencil: full/ghost/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.617 | 4.617 | 4.617 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 8000 -36411.315 -36411.315 -36710.937 -36710.937 289.78279 289.78279 160948.55 160948.55 -353.59902 -353.59902 9000 -36419.791 -36419.791 -36722.455 -36722.455 292.72549 292.72549 160956.15 160956.15 -300.21656 -300.21656 Loop time of 57.0108 on 1 procs for 1000 steps with 8000 atoms Performance: 1.516 ns/day, 15.836 hours/ns, 17.541 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 | 56.24 | 56.24 | 56.24 | 0.0 | 98.65 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.048538 | 0.048538 | 0.048538 | 0.0 | 0.09 Output | 0.0001436 | 0.0001436 | 0.0001436 | 0.0 | 0.00 Modify | 0.647 | 0.647 | 0.647 | 0.0 | 1.13 Other | | 0.0756 | | | 0.13 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4995.00 ave 4995 max 4995 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: 224000.0 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 28.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 = 293.085748186934, Press = 0.643118653030338 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 = 4.96 ghost atom cutoff = 4.96 binsize = 2.48, bins = 22 22 22 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, ghost, cut 4.96 pair build: full/bin/ghost stencil: full/ghost/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.617 | 4.617 | 4.617 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 9000 -36419.791 -36419.791 -36722.455 -36722.455 292.72549 292.72549 160956.15 160956.15 -300.21656 -300.21656 10000 -36414.285 -36414.285 -36720.694 -36720.694 296.34746 296.34746 160975.85 160975.85 -356.91212 -356.91212 Loop time of 57.0234 on 1 procs for 1000 steps with 8000 atoms Performance: 1.515 ns/day, 15.840 hours/ns, 17.537 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 | 56.248 | 56.248 | 56.248 | 0.0 | 98.64 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.048507 | 0.048507 | 0.048507 | 0.0 | 0.09 Output | 0.00015425 | 0.00015425 | 0.00015425 | 0.0 | 0.00 Modify | 0.6486 | 0.6486 | 0.6486 | 0.0 | 1.14 Other | | 0.07806 | | | 0.14 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4995.00 ave 4995 max 4995 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: 224000.0 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 28.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 = 293.144835330427, Press = -2.18184585628892 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 = 4.96 ghost atom cutoff = 4.96 binsize = 2.48, bins = 22 22 22 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, ghost, cut 4.96 pair build: full/bin/ghost stencil: full/ghost/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.617 | 4.617 | 4.617 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 10000 -36414.285 -36414.285 -36720.694 -36720.694 296.34746 296.34746 160975.85 160975.85 -356.91212 -356.91212 11000 -36423.261 -36423.261 -36717.909 -36717.909 284.97239 284.97239 160986 160986 -557.13919 -557.13919 Loop time of 56.9909 on 1 procs for 1000 steps with 8000 atoms Performance: 1.516 ns/day, 15.831 hours/ns, 17.547 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 | 56.214 | 56.214 | 56.214 | 0.0 | 98.64 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.048906 | 0.048906 | 0.048906 | 0.0 | 0.09 Output | 0.00022457 | 0.00022457 | 0.00022457 | 0.0 | 0.00 Modify | 0.64963 | 0.64963 | 0.64963 | 0.0 | 1.14 Other | | 0.07837 | | | 0.14 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4995.00 ave 4995 max 4995 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: 224000.0 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 28.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 = 293.107884750974, Press = -1.97996841044126 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 = 4.96 ghost atom cutoff = 4.96 binsize = 2.48, bins = 22 22 22 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, ghost, cut 4.96 pair build: full/bin/ghost stencil: full/ghost/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.617 | 4.617 | 4.617 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 11000 -36423.261 -36423.261 -36717.909 -36717.909 284.97239 284.97239 160986 160986 -557.13919 -557.13919 12000 -36418.522 -36418.522 -36722.291 -36722.291 293.79362 293.79362 160952.71 160952.71 -285.98298 -285.98298 Loop time of 57.1067 on 1 procs for 1000 steps with 8000 atoms Performance: 1.513 ns/day, 15.863 hours/ns, 17.511 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 | 56.327 | 56.327 | 56.327 | 0.0 | 98.63 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.048791 | 0.048791 | 0.048791 | 0.0 | 0.09 Output | 0.00014365 | 0.00014365 | 0.00014365 | 0.0 | 0.00 Modify | 0.65217 | 0.65217 | 0.65217 | 0.0 | 1.14 Other | | 0.07873 | | | 0.14 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4995.00 ave 4995 max 4995 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: 224000.0 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 28.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 = 293.088597203492, Press = -2.23598294200227 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 = 4.96 ghost atom cutoff = 4.96 binsize = 2.48, bins = 22 22 22 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, ghost, cut 4.96 pair build: full/bin/ghost stencil: full/ghost/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.617 | 4.617 | 4.617 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 12000 -36418.522 -36418.522 -36722.291 -36722.291 293.79362 293.79362 160952.71 160952.71 -285.98298 -285.98298 13000 -36414.954 -36414.954 -36712.432 -36712.432 287.7095 287.7095 160911.25 160911.25 -17.577058 -17.577058 Loop time of 57.4658 on 1 procs for 1000 steps with 8000 atoms Performance: 1.504 ns/day, 15.963 hours/ns, 17.402 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 | 56.68 | 56.68 | 56.68 | 0.0 | 98.63 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.048825 | 0.048825 | 0.048825 | 0.0 | 0.08 Output | 0.00023402 | 0.00023402 | 0.00023402 | 0.0 | 0.00 Modify | 0.6578 | 0.6578 | 0.6578 | 0.0 | 1.14 Other | | 0.0793 | | | 0.14 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4995.00 ave 4995 max 4995 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: 224000.0 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 28.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 = 293.096070640782, Press = 0.280278940710547 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 = 4.96 ghost atom cutoff = 4.96 binsize = 2.48, bins = 22 22 22 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, ghost, cut 4.96 pair build: full/bin/ghost stencil: full/ghost/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.617 | 4.617 | 4.617 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 13000 -36414.954 -36414.954 -36712.432 -36712.432 287.7095 287.7095 160911.25 160911.25 -17.577058 -17.577058 14000 -36413.873 -36413.873 -36719.23 -36719.23 295.33025 295.33025 160892.86 160892.86 213.69458 213.69458 Loop time of 57.0753 on 1 procs for 1000 steps with 8000 atoms Performance: 1.514 ns/day, 15.854 hours/ns, 17.521 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 | 56.295 | 56.295 | 56.295 | 0.0 | 98.63 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.048324 | 0.048324 | 0.048324 | 0.0 | 0.08 Output | 0.0001327 | 0.0001327 | 0.0001327 | 0.0 | 0.00 Modify | 0.65286 | 0.65286 | 0.65286 | 0.0 | 1.14 Other | | 0.07859 | | | 0.14 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4995.00 ave 4995 max 4995 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: 224000.0 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 28.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 = 293.129158083649, Press = 1.36966888703017 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 = 4.96 ghost atom cutoff = 4.96 binsize = 2.48, bins = 22 22 22 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, ghost, cut 4.96 pair build: full/bin/ghost stencil: full/ghost/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.617 | 4.617 | 4.617 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 14000 -36413.873 -36413.873 -36719.23 -36719.23 295.33025 295.33025 160892.86 160892.86 213.69458 213.69458 15000 -36416.507 -36416.507 -36721.921 -36721.921 295.38449 295.38449 160890.38 160890.38 50.247948 50.247948 Loop time of 54.7428 on 1 procs for 1000 steps with 8000 atoms Performance: 1.578 ns/day, 15.206 hours/ns, 18.267 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.011 | 54.011 | 54.011 | 0.0 | 98.66 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.046633 | 0.046633 | 0.046633 | 0.0 | 0.09 Output | 0.00010712 | 0.00010712 | 0.00010712 | 0.0 | 0.00 Modify | 0.61447 | 0.61447 | 0.61447 | 0.0 | 1.12 Other | | 0.07039 | | | 0.13 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4995.00 ave 4995 max 4995 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: 224000.0 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 28.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 = 293.296005541094, Press = 1.59250796317072 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 = 4.96 ghost atom cutoff = 4.96 binsize = 2.48, bins = 22 22 22 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, ghost, cut 4.96 pair build: full/bin/ghost stencil: full/ghost/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.617 | 4.617 | 4.617 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 15000 -36416.507 -36416.507 -36721.921 -36721.921 295.38449 295.38449 160890.38 160890.38 50.247948 50.247948 16000 -36413.648 -36413.648 -36717.147 -36717.147 293.53255 293.53255 160866.72 160866.72 294.37901 294.37901 Loop time of 56.8294 on 1 procs for 1000 steps with 8000 atoms Performance: 1.520 ns/day, 15.786 hours/ns, 17.597 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 | 56.056 | 56.056 | 56.056 | 0.0 | 98.64 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.048278 | 0.048278 | 0.048278 | 0.0 | 0.08 Output | 0.00014492 | 0.00014492 | 0.00014492 | 0.0 | 0.00 Modify | 0.64721 | 0.64721 | 0.64721 | 0.0 | 1.14 Other | | 0.07801 | | | 0.14 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4995.00 ave 4995 max 4995 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: 224000.0 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 28.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 = 293.268387333416, Press = 2.53480466930454 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 = 4.96 ghost atom cutoff = 4.96 binsize = 2.48, bins = 22 22 22 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, ghost, cut 4.96 pair build: full/bin/ghost stencil: full/ghost/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.617 | 4.617 | 4.617 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 16000 -36413.648 -36413.648 -36717.147 -36717.147 293.53255 293.53255 160866.72 160866.72 294.37901 294.37901 17000 -36421.036 -36421.036 -36718.791 -36718.791 287.97797 287.97797 160861.39 160861.39 232.95436 232.95436 Loop time of 56.1701 on 1 procs for 1000 steps with 8000 atoms Performance: 1.538 ns/day, 15.603 hours/ns, 17.803 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 | 55.41 | 55.41 | 55.41 | 0.0 | 98.65 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.047588 | 0.047588 | 0.047588 | 0.0 | 0.08 Output | 0.00013578 | 0.00013578 | 0.00013578 | 0.0 | 0.00 Modify | 0.63689 | 0.63689 | 0.63689 | 0.0 | 1.13 Other | | 0.07527 | | | 0.13 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4995.00 ave 4995 max 4995 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: 224000.0 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 28.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 = 293.298505723236, Press = 2.99997392015803 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 = 4.96 ghost atom cutoff = 4.96 binsize = 2.48, bins = 22 22 22 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, ghost, cut 4.96 pair build: full/bin/ghost stencil: full/ghost/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.617 | 4.617 | 4.617 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 17000 -36421.036 -36421.036 -36718.791 -36718.791 287.97797 287.97797 160861.39 160861.39 232.95436 232.95436 18000 -36413.815 -36413.815 -36719.268 -36719.268 295.42319 295.42319 160936.05 160936.05 -152.82584 -152.82584 Loop time of 50.5542 on 1 procs for 1000 steps with 8000 atoms Performance: 1.709 ns/day, 14.043 hours/ns, 19.781 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 | 49.898 | 49.898 | 49.898 | 0.0 | 98.70 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.043619 | 0.043619 | 0.043619 | 0.0 | 0.09 Output | 0.00021114 | 0.00021114 | 0.00021114 | 0.0 | 0.00 Modify | 0.55381 | 0.55381 | 0.55381 | 0.0 | 1.10 Other | | 0.05836 | | | 0.12 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4995.00 ave 4995 max 4995 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: 224000.0 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 28.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 = 293.243306515084, Press = 2.37445369284117 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 = 4.96 ghost atom cutoff = 4.96 binsize = 2.48, bins = 22 22 22 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, ghost, cut 4.96 pair build: full/bin/ghost stencil: full/ghost/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.617 | 4.617 | 4.617 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 18000 -36413.815 -36413.815 -36719.268 -36719.268 295.42319 295.42319 160936.05 160936.05 -152.82584 -152.82584 19000 -36422.947 -36422.947 -36723.887 -36723.887 291.05813 291.05813 160939.86 160939.86 -199.62775 -199.62775 Loop time of 50.5505 on 1 procs for 1000 steps with 8000 atoms Performance: 1.709 ns/day, 14.042 hours/ns, 19.782 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 | 49.895 | 49.895 | 49.895 | 0.0 | 98.70 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.043662 | 0.043662 | 0.043662 | 0.0 | 0.09 Output | 9.7008e-05 | 9.7008e-05 | 9.7008e-05 | 0.0 | 0.00 Modify | 0.55363 | 0.55363 | 0.55363 | 0.0 | 1.10 Other | | 0.05828 | | | 0.12 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4995.00 ave 4995 max 4995 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: 224000.0 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 28.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 = 293.130836335768, Press = 1.54972365432019 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 = 4.96 ghost atom cutoff = 4.96 binsize = 2.48, bins = 22 22 22 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, ghost, cut 4.96 pair build: full/bin/ghost stencil: full/ghost/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.617 | 4.617 | 4.617 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 19000 -36422.947 -36422.947 -36723.887 -36723.887 291.05813 291.05813 160939.86 160939.86 -199.62775 -199.62775 20000 -36418.426 -36418.426 -36715.807 -36715.807 287.61586 287.61586 160935.17 160935.17 -176.50314 -176.50314 Loop time of 50.5697 on 1 procs for 1000 steps with 8000 atoms Performance: 1.709 ns/day, 14.047 hours/ns, 19.775 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 | 49.912 | 49.912 | 49.912 | 0.0 | 98.70 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.043682 | 0.043682 | 0.043682 | 0.0 | 0.09 Output | 0.00013577 | 0.00013577 | 0.00013577 | 0.0 | 0.00 Modify | 0.55499 | 0.55499 | 0.55499 | 0.0 | 1.10 Other | | 0.05876 | | | 0.12 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4995.00 ave 4995 max 4995 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: 224000.0 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 28.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 = 293.134399492867, Press = 0.982774584654763 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 = 4.96 ghost atom cutoff = 4.96 binsize = 2.48, bins = 22 22 22 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, ghost, cut 4.96 pair build: full/bin/ghost stencil: full/ghost/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.617 | 4.617 | 4.617 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 20000 -36418.426 -36418.426 -36715.807 -36715.807 287.61586 287.61586 160935.17 160935.17 -176.50314 -176.50314 21000 -36414.504 -36414.504 -36720.231 -36720.231 295.68706 295.68706 160966.52 160966.52 -390.97419 -390.97419 Loop time of 50.5476 on 1 procs for 1000 steps with 8000 atoms Performance: 1.709 ns/day, 14.041 hours/ns, 19.783 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 | 49.891 | 49.891 | 49.891 | 0.0 | 98.70 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.043497 | 0.043497 | 0.043497 | 0.0 | 0.09 Output | 0.00012813 | 0.00012813 | 0.00012813 | 0.0 | 0.00 Modify | 0.55443 | 0.55443 | 0.55443 | 0.0 | 1.10 Other | | 0.05848 | | | 0.12 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4995.00 ave 4995 max 4995 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: 224000.0 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 28.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 = 293.064271288901, Press = 0.770721727122356 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 = 4.96 ghost atom cutoff = 4.96 binsize = 2.48, bins = 22 22 22 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, ghost, cut 4.96 pair build: full/bin/ghost stencil: full/ghost/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.617 | 4.617 | 4.617 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 21000 -36414.504 -36414.504 -36720.231 -36720.231 295.68706 295.68706 160966.52 160966.52 -390.97419 -390.97419 22000 -36415.709 -36415.709 -36718.894 -36718.894 293.22835 293.22835 160997.71 160997.71 -518.38133 -518.38133 Loop time of 50.5497 on 1 procs for 1000 steps with 8000 atoms Performance: 1.709 ns/day, 14.042 hours/ns, 19.783 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 | 49.893 | 49.893 | 49.893 | 0.0 | 98.70 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.043879 | 0.043879 | 0.043879 | 0.0 | 0.09 Output | 0.00012187 | 0.00012187 | 0.00012187 | 0.0 | 0.00 Modify | 0.55461 | 0.55461 | 0.55461 | 0.0 | 1.10 Other | | 0.05851 | | | 0.12 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4995.00 ave 4995 max 4995 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: 224000.0 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 28.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 = 293.033312712007, Press = 0.3553744785194 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 = 4.96 ghost atom cutoff = 4.96 binsize = 2.48, bins = 22 22 22 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, ghost, cut 4.96 pair build: full/bin/ghost stencil: full/ghost/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.617 | 4.617 | 4.617 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 22000 -36415.709 -36415.709 -36718.894 -36718.894 293.22835 293.22835 160997.71 160997.71 -518.38133 -518.38133 23000 -36420.781 -36420.781 -36724.413 -36724.413 293.66133 293.66133 161047.73 161047.73 -838.44817 -838.44817 Loop time of 50.536 on 1 procs for 1000 steps with 8000 atoms Performance: 1.710 ns/day, 14.038 hours/ns, 19.788 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 | 49.88 | 49.88 | 49.88 | 0.0 | 98.70 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.043759 | 0.043759 | 0.043759 | 0.0 | 0.09 Output | 9.4399e-05 | 9.4399e-05 | 9.4399e-05 | 0.0 | 0.00 Modify | 0.55423 | 0.55423 | 0.55423 | 0.0 | 1.10 Other | | 0.05819 | | | 0.12 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4995.00 ave 4995 max 4995 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: 224000.0 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 28.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 = 293.063112990095, Press = -0.324377990342751 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 = 4.96 ghost atom cutoff = 4.96 binsize = 2.48, bins = 22 22 22 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, ghost, cut 4.96 pair build: full/bin/ghost stencil: full/ghost/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.617 | 4.617 | 4.617 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 23000 -36420.781 -36420.781 -36724.413 -36724.413 293.66133 293.66133 161047.73 161047.73 -838.44817 -838.44817 24000 -36414.257 -36414.257 -36719.214 -36719.214 294.94312 294.94312 161041.62 161041.62 -817.83555 -817.83555 Loop time of 50.5537 on 1 procs for 1000 steps with 8000 atoms Performance: 1.709 ns/day, 14.043 hours/ns, 19.781 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 | 49.897 | 49.897 | 49.897 | 0.0 | 98.70 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.043812 | 0.043812 | 0.043812 | 0.0 | 0.09 Output | 0.00012828 | 0.00012828 | 0.00012828 | 0.0 | 0.00 Modify | 0.55434 | 0.55434 | 0.55434 | 0.0 | 1.10 Other | | 0.05843 | | | 0.12 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4995.00 ave 4995 max 4995 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: 224000.0 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 28.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 = 293.047318881742, Press = -0.622739852693009 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 = 4.96 ghost atom cutoff = 4.96 binsize = 2.48, bins = 22 22 22 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, ghost, cut 4.96 pair build: full/bin/ghost stencil: full/ghost/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.617 | 4.617 | 4.617 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 24000 -36414.257 -36414.257 -36719.214 -36719.214 294.94312 294.94312 161041.62 161041.62 -817.83555 -817.83555 25000 -36416.69 -36416.69 -36717.728 -36717.728 291.15306 291.15306 160972.16 160972.16 -409.15657 -409.15657 Loop time of 50.57 on 1 procs for 1000 steps with 8000 atoms Performance: 1.709 ns/day, 14.047 hours/ns, 19.775 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 | 49.913 | 49.913 | 49.913 | 0.0 | 98.70 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.043743 | 0.043743 | 0.043743 | 0.0 | 0.09 Output | 0.0001243 | 0.0001243 | 0.0001243 | 0.0 | 0.00 Modify | 0.55427 | 0.55427 | 0.55427 | 0.0 | 1.10 Other | | 0.05839 | | | 0.12 Nlocal: 8000.00 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4995.00 ave 4995 max 4995 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: 224000.0 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 28.000000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" jump SELF break # Write final averaged volume to file if temperature and volume have converged; otherwise wirte a # flag to indicate non-convergence. variable myStep equal step if "${myStep} < 2000000" then "print '${V}' file output/vol_T293.15.out" else "print 'not_converged' file output/vol_T293.15.out" print '${V}' file output/vol_T293.15.out 160908.064084862 print "LAMMPS calculation completed" LAMMPS calculation completed quit 0