# 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.430545687675476*${_u_distance} variable latticeconst_converted equal 5.430545687675476*1 lattice diamond ${latticeconst_converted} lattice diamond 5.43054568767548 Lattice spacing in x,y,z = 5.4305457 5.4305457 5.4305457 region simbox block 0 10 0 10 0 10 units lattice create_box 1 simbox Created orthogonal box = (0 0 0) to (54.305457 54.305457 54.305457) 1 by 1 by 1 MPI processor grid create_atoms 1 box Created 8000 atoms using lattice units in orthogonal box = (0 0 0) to (54.305457 54.305457 54.305457) create_atoms CPU = 0.012 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 MEAM_LAMMPS_LenoskySadighAlonso_2000_Si__MO_533426548156_002 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 160151.280490546 variable V0_metal equal ${V0}/(${_u_distance}*${_u_distance}*${_u_distance}) variable V0_metal equal 160151.280490546/(${_u_distance}*${_u_distance}*${_u_distance}) variable V0_metal equal 160151.280490546/(1*${_u_distance}*${_u_distance}) variable V0_metal equal 160151.280490546/(1*1*${_u_distance}) variable V0_metal equal 160151.280490546/(1*1*1) variable V0_metal_times1000 equal ${V0_metal}*1000 variable V0_metal_times1000 equal 160151.280490546*1000 print "Initial system volume: ${V0_metal} Angstroms^3" Initial system volume: 160151.280490546 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 273.15*${_u_temperature} variable temp_converted equal 273.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 273.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 273.15 ${temp_converted} ${Tdamp_converted} iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 273.15 273.15 ${Tdamp_converted} iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 273.15 273.15 0.1 iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 273.15 273.15 0.1 iso 0 ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 273.15 273.15 0.1 iso 0 0 ${Pdamp_converted} fix ensemble all npt temp 273.15 273.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 "273.15 - 0.2" variable T_up equal "273.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_533426548156_002#item-citation CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 17 17 17 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.677 | 4.677 | 4.677 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 0 -36571.936 -36571.936 -36854.36 -36854.36 273.15 273.15 160151.28 160151.28 1883.5973 1883.5973 1000 -36251.693 -36251.693 -36535.743 -36535.743 274.72296 274.72296 158697.48 158697.48 532.3892 532.3892 Loop time of 364.885 on 1 procs for 1000 steps with 8000 atoms Performance: 0.237 ns/day, 101.357 hours/ns, 2.741 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 | 363.17 | 363.17 | 363.17 | 0.0 | 99.53 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.21382 | 0.21382 | 0.21382 | 0.0 | 0.06 Output | 0.000337 | 0.000337 | 0.000337 | 0.0 | 0.00 Modify | 1.3367 | 1.3367 | 1.3367 | 0.0 | 0.37 Other | | 0.163 | | | 0.04 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 368000 ave 368000 max 368000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 368000 Ave neighs/atom = 46 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 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 17 17 17 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.679 | 4.679 | 4.679 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 1000 -36251.693 -36251.693 -36535.743 -36535.743 274.72296 274.72296 158697.48 158697.48 532.3892 532.3892 2000 -36284.989 -36284.989 -36565.013 -36565.013 270.82911 270.82911 158932.05 158932.05 -224.47693 -224.47693 Loop time of 370.668 on 1 procs for 1000 steps with 8000 atoms Performance: 0.233 ns/day, 102.963 hours/ns, 2.698 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 | 368.95 | 368.95 | 368.95 | 0.0 | 99.54 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.21405 | 0.21405 | 0.21405 | 0.0 | 0.06 Output | 0.00027513 | 0.00027513 | 0.00027513 | 0.0 | 0.00 Modify | 1.3432 | 1.3432 | 1.3432 | 0.0 | 0.36 Other | | 0.1644 | | | 0.04 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 384454 ave 384454 max 384454 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 384454 Ave neighs/atom = 48.05675 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 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 17 17 17 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.679 | 4.679 | 4.679 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 2000 -36284.989 -36284.989 -36565.013 -36565.013 270.82911 270.82911 158932.05 158932.05 -224.47693 -224.47693 3000 -36261.526 -36261.526 -36547.322 -36547.322 276.41091 276.41091 158884.14 158884.14 -615.39575 -615.39575 Loop time of 370.234 on 1 procs for 1000 steps with 8000 atoms Performance: 0.233 ns/day, 102.843 hours/ns, 2.701 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 | 368.52 | 368.52 | 368.52 | 0.0 | 99.54 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.2148 | 0.2148 | 0.2148 | 0.0 | 0.06 Output | 0.0002742 | 0.0002742 | 0.0002742 | 0.0 | 0.00 Modify | 1.3391 | 1.3391 | 1.3391 | 0.0 | 0.36 Other | | 0.1626 | | | 0.04 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 382204 ave 382204 max 382204 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 382204 Ave neighs/atom = 47.7755 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 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 17 17 17 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.679 | 4.679 | 4.679 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 3000 -36261.526 -36261.526 -36547.322 -36547.322 276.41091 276.41091 158884.14 158884.14 -615.39575 -615.39575 4000 -36280.03 -36280.03 -36557.684 -36557.684 268.53667 268.53667 158890.67 158890.67 -185.63102 -185.63102 Loop time of 368.9 on 1 procs for 1000 steps with 8000 atoms Performance: 0.234 ns/day, 102.472 hours/ns, 2.711 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 | 367.19 | 367.19 | 367.19 | 0.0 | 99.54 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.21166 | 0.21166 | 0.21166 | 0.0 | 0.06 Output | 0.00035705 | 0.00035705 | 0.00035705 | 0.0 | 0.00 Modify | 1.333 | 1.333 | 1.333 | 0.0 | 0.36 Other | | 0.1624 | | | 0.04 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 383670 ave 383670 max 383670 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 383670 Ave neighs/atom = 47.95875 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 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 17 17 17 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.679 | 4.679 | 4.679 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 4000 -36280.03 -36280.03 -36557.684 -36557.684 268.53667 268.53667 158890.67 158890.67 -185.63102 -185.63102 5000 -36266.192 -36266.192 -36548.254 -36548.254 272.80003 272.80003 158805.37 158805.37 -58.318968 -58.318968 Loop time of 369.745 on 1 procs for 1000 steps with 8000 atoms Performance: 0.234 ns/day, 102.707 hours/ns, 2.705 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 | 368.03 | 368.03 | 368.03 | 0.0 | 99.54 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.21573 | 0.21573 | 0.21573 | 0.0 | 0.06 Output | 0.00023038 | 0.00023038 | 0.00023038 | 0.0 | 0.00 Modify | 1.3406 | 1.3406 | 1.3406 | 0.0 | 0.36 Other | | 0.1626 | | | 0.04 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 382670 ave 382670 max 382670 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 382670 Ave neighs/atom = 47.83375 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 = 276.462339385153, Press = -266.647676560377 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 17 17 17 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.679 | 4.679 | 4.679 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 5000 -36266.192 -36266.192 -36548.254 -36548.254 272.80003 272.80003 158805.37 158805.37 -58.318968 -58.318968 6000 -36276.212 -36276.212 -36554.843 -36554.843 269.48169 269.48169 158889.68 158889.68 -343.71216 -343.71216 Loop time of 369.448 on 1 procs for 1000 steps with 8000 atoms Performance: 0.234 ns/day, 102.624 hours/ns, 2.707 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 | 367.7 | 367.7 | 367.7 | 0.0 | 99.53 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.21266 | 0.21266 | 0.21266 | 0.0 | 0.06 Output | 0.00028366 | 0.00028366 | 0.00028366 | 0.0 | 0.00 Modify | 1.3748 | 1.3748 | 1.3748 | 0.0 | 0.37 Other | | 0.1623 | | | 0.04 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 383476 ave 383476 max 383476 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 383476 Ave neighs/atom = 47.9345 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 = 272.994271106161, Press = -29.8168063805911 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 17 17 17 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.679 | 4.679 | 4.679 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 6000 -36276.212 -36276.212 -36554.843 -36554.843 269.48169 269.48169 158889.68 158889.68 -343.71216 -343.71216 7000 -36268.945 -36268.945 -36552.467 -36552.467 274.21194 274.21194 158845.48 158845.48 -49.261186 -49.261186 Loop time of 369.409 on 1 procs for 1000 steps with 8000 atoms Performance: 0.234 ns/day, 102.614 hours/ns, 2.707 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 | 367.67 | 367.67 | 367.67 | 0.0 | 99.53 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.21162 | 0.21162 | 0.21162 | 0.0 | 0.06 Output | 0.0002875 | 0.0002875 | 0.0002875 | 0.0 | 0.00 Modify | 1.3681 | 1.3681 | 1.3681 | 0.0 | 0.37 Other | | 0.1632 | | | 0.04 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 382994 ave 382994 max 382994 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 382994 Ave neighs/atom = 47.87425 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 = 273.308932886114, Press = 0.45499059295584 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 17 17 17 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.679 | 4.679 | 4.679 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 7000 -36268.945 -36268.945 -36552.467 -36552.467 274.21194 274.21194 158845.48 158845.48 -49.261186 -49.261186 8000 -36273.579 -36273.579 -36550.929 -36550.929 268.24251 268.24251 158647.22 158647.22 613.69442 613.69442 Loop time of 369.837 on 1 procs for 1000 steps with 8000 atoms Performance: 0.234 ns/day, 102.733 hours/ns, 2.704 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 | 368.09 | 368.09 | 368.09 | 0.0 | 99.53 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.21196 | 0.21196 | 0.21196 | 0.0 | 0.06 Output | 0.0002327 | 0.0002327 | 0.0002327 | 0.0 | 0.00 Modify | 1.3714 | 1.3714 | 1.3714 | 0.0 | 0.37 Other | | 0.162 | | | 0.04 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 383006 ave 383006 max 383006 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 383006 Ave neighs/atom = 47.87575 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 = 273.110132530826, Press = -8.36140815851527 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 17 17 17 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.679 | 4.679 | 4.679 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 8000 -36273.579 -36273.579 -36550.929 -36550.929 268.24251 268.24251 158647.22 158647.22 613.69442 613.69442 9000 -36272.069 -36272.069 -36555.209 -36555.209 273.84315 273.84315 158843.87 158843.87 -125.66895 -125.66895 Loop time of 368.693 on 1 procs for 1000 steps with 8000 atoms Performance: 0.234 ns/day, 102.415 hours/ns, 2.712 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 | 366.95 | 366.95 | 366.95 | 0.0 | 99.53 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.21324 | 0.21324 | 0.21324 | 0.0 | 0.06 Output | 0.00023781 | 0.00023781 | 0.00023781 | 0.0 | 0.00 Modify | 1.3713 | 1.3713 | 1.3713 | 0.0 | 0.37 Other | | 0.1626 | | | 0.04 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 384218 ave 384218 max 384218 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 384218 Ave neighs/atom = 48.02725 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 = 273.208808180703, Press = -3.0061283070099 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 17 17 17 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.679 | 4.679 | 4.679 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 9000 -36272.069 -36272.069 -36555.209 -36555.209 273.84315 273.84315 158843.87 158843.87 -125.66895 -125.66895 10000 -36270.703 -36270.703 -36553.15 -36553.15 273.17203 273.17203 158816.54 158816.54 -136.51756 -136.51756 Loop time of 369.144 on 1 procs for 1000 steps with 8000 atoms Performance: 0.234 ns/day, 102.540 hours/ns, 2.709 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 | 367.4 | 367.4 | 367.4 | 0.0 | 99.53 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.2121 | 0.2121 | 0.2121 | 0.0 | 0.06 Output | 0.00034667 | 0.00034667 | 0.00034667 | 0.0 | 0.00 Modify | 1.3712 | 1.3712 | 1.3712 | 0.0 | 0.37 Other | | 0.1634 | | | 0.04 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 383050 ave 383050 max 383050 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 383050 Ave neighs/atom = 47.88125 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 = 273.126675300096, Press = -8.22580882586744 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 17 17 17 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.679 | 4.679 | 4.679 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 10000 -36270.703 -36270.703 -36553.15 -36553.15 273.17203 273.17203 158816.54 158816.54 -136.51756 -136.51756 11000 -36271.429 -36271.429 -36554.216 -36554.216 273.50047 273.50047 158803.21 158803.21 4.3959858 4.3959858 Loop time of 369.599 on 1 procs for 1000 steps with 8000 atoms Performance: 0.234 ns/day, 102.666 hours/ns, 2.706 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 | 367.86 | 367.86 | 367.86 | 0.0 | 99.53 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.21315 | 0.21315 | 0.21315 | 0.0 | 0.06 Output | 0.00022618 | 0.00022618 | 0.00022618 | 0.0 | 0.00 Modify | 1.359 | 1.359 | 1.359 | 0.0 | 0.37 Other | | 0.1625 | | | 0.04 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 383534 ave 383534 max 383534 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 383534 Ave neighs/atom = 47.94175 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 = 273.235604617161, Press = 1.95995401694425 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 17 17 17 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.679 | 4.679 | 4.679 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 11000 -36271.429 -36271.429 -36554.216 -36554.216 273.50047 273.50047 158803.21 158803.21 4.3959858 4.3959858 12000 -36269.821 -36269.821 -36554.029 -36554.029 274.87574 274.87574 158743.47 158743.47 252.37459 252.37459 Loop time of 369.655 on 1 procs for 1000 steps with 8000 atoms Performance: 0.234 ns/day, 102.682 hours/ns, 2.705 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 | 367.91 | 367.91 | 367.91 | 0.0 | 99.53 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.21241 | 0.21241 | 0.21241 | 0.0 | 0.06 Output | 0.0002356 | 0.0002356 | 0.0002356 | 0.0 | 0.00 Modify | 1.368 | 1.368 | 1.368 | 0.0 | 0.37 Other | | 0.1618 | | | 0.04 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 383452 ave 383452 max 383452 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 383452 Ave neighs/atom = 47.9315 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 = 273.214348632995, Press = -5.86415966491388 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 17 17 17 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.679 | 4.679 | 4.679 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 12000 -36269.821 -36269.821 -36554.029 -36554.029 274.87574 274.87574 158743.47 158743.47 252.37459 252.37459 13000 -36270.653 -36270.653 -36552.617 -36552.617 272.70468 272.70468 158726.01 158726.01 139.5204 139.5204 Loop time of 370.745 on 1 procs for 1000 steps with 8000 atoms Performance: 0.233 ns/day, 102.985 hours/ns, 2.697 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 | 369 | 369 | 369 | 0.0 | 99.53 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.21007 | 0.21007 | 0.21007 | 0.0 | 0.06 Output | 0.0002312 | 0.0002312 | 0.0002312 | 0.0 | 0.00 Modify | 1.3699 | 1.3699 | 1.3699 | 0.0 | 0.37 Other | | 0.1614 | | | 0.04 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 383734 ave 383734 max 383734 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 383734 Ave neighs/atom = 47.96675 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 = 273.371354145941, Press = -1.29428368225591 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 17 17 17 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.679 | 4.679 | 4.679 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 13000 -36270.653 -36270.653 -36552.617 -36552.617 272.70468 272.70468 158726.01 158726.01 139.5204 139.5204 14000 -36275.196 -36275.196 -36557.889 -36557.889 273.40974 273.40974 158837.23 158837.23 -186.92182 -186.92182 Loop time of 370.812 on 1 procs for 1000 steps with 8000 atoms Performance: 0.233 ns/day, 103.003 hours/ns, 2.697 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 | 369.06 | 369.06 | 369.06 | 0.0 | 99.53 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.21283 | 0.21283 | 0.21283 | 0.0 | 0.06 Output | 0.00023458 | 0.00023458 | 0.00023458 | 0.0 | 0.00 Modify | 1.3742 | 1.3742 | 1.3742 | 0.0 | 0.37 Other | | 0.1628 | | | 0.04 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 384016 ave 384016 max 384016 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 384016 Ave neighs/atom = 48.002 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 = 273.332704665437, Press = -2.35384844938902 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 17 17 17 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.679 | 4.679 | 4.679 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 14000 -36275.196 -36275.196 -36557.889 -36557.889 273.40974 273.40974 158837.23 158837.23 -186.92182 -186.92182 15000 -36269.393 -36269.393 -36560.06 -36560.06 281.12295 281.12295 158846.16 158846.16 -113.33758 -113.33758 Loop time of 370.299 on 1 procs for 1000 steps with 8000 atoms Performance: 0.233 ns/day, 102.861 hours/ns, 2.701 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 | 368.55 | 368.55 | 368.55 | 0.0 | 99.53 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.21566 | 0.21566 | 0.21566 | 0.0 | 0.06 Output | 0.00023879 | 0.00023879 | 0.00023879 | 0.0 | 0.00 Modify | 1.375 | 1.375 | 1.375 | 0.0 | 0.37 Other | | 0.1624 | | | 0.04 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 383054 ave 383054 max 383054 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 383054 Ave neighs/atom = 47.88175 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 = 273.272325758867, Press = -1.10775043010717 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 17 17 17 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.679 | 4.679 | 4.679 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 15000 -36269.393 -36269.393 -36560.06 -36560.06 281.12295 281.12295 158846.16 158846.16 -113.33758 -113.33758 16000 -36277.218 -36277.218 -36553.911 -36553.911 267.6077 267.6077 158799.7 158799.7 -235.35338 -235.35338 Loop time of 371.693 on 1 procs for 1000 steps with 8000 atoms Performance: 0.232 ns/day, 103.248 hours/ns, 2.690 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 | 369.93 | 369.93 | 369.93 | 0.0 | 99.53 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.21171 | 0.21171 | 0.21171 | 0.0 | 0.06 Output | 0.00023482 | 0.00023482 | 0.00023482 | 0.0 | 0.00 Modify | 1.383 | 1.383 | 1.383 | 0.0 | 0.37 Other | | 0.1644 | | | 0.04 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 382810 ave 382810 max 382810 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 382810 Ave neighs/atom = 47.85125 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 = 273.189322611923, Press = -1.31777328049207 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 17 17 17 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.679 | 4.679 | 4.679 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 16000 -36277.218 -36277.218 -36553.911 -36553.911 267.6077 267.6077 158799.7 158799.7 -235.35338 -235.35338 17000 -36270.634 -36270.634 -36551.821 -36551.821 271.95359 271.95359 158870.49 158870.49 -501.10952 -501.10952 Loop time of 370.309 on 1 procs for 1000 steps with 8000 atoms Performance: 0.233 ns/day, 102.864 hours/ns, 2.700 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 | 368.56 | 368.56 | 368.56 | 0.0 | 99.53 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.21235 | 0.21235 | 0.21235 | 0.0 | 0.06 Output | 0.00023324 | 0.00023324 | 0.00023324 | 0.0 | 0.00 Modify | 1.3759 | 1.3759 | 1.3759 | 0.0 | 0.37 Other | | 0.1642 | | | 0.04 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 383994 ave 383994 max 383994 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 383994 Ave neighs/atom = 47.99925 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 = 273.18634913539, Press = -3.20440313177325 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 17 17 17 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.679 | 4.679 | 4.679 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 17000 -36270.634 -36270.634 -36551.821 -36551.821 271.95359 271.95359 158870.49 158870.49 -501.10952 -501.10952 18000 -36271.024 -36271.024 -36553.792 -36553.792 273.48272 273.48272 158839.16 158839.16 -104.34565 -104.34565 Loop time of 368.523 on 1 procs for 1000 steps with 8000 atoms Performance: 0.234 ns/day, 102.368 hours/ns, 2.714 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 | 366.79 | 366.79 | 366.79 | 0.0 | 99.53 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.21075 | 0.21075 | 0.21075 | 0.0 | 0.06 Output | 0.00024143 | 0.00024143 | 0.00024143 | 0.0 | 0.00 Modify | 1.3611 | 1.3611 | 1.3611 | 0.0 | 0.37 Other | | 0.1623 | | | 0.04 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 383422 ave 383422 max 383422 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 383422 Ave neighs/atom = 47.92775 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 = 273.169697886443, Press = 0.576887673679553 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 17 17 17 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.679 | 4.679 | 4.679 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 18000 -36271.024 -36271.024 -36553.792 -36553.792 273.48272 273.48272 158839.16 158839.16 -104.34565 -104.34565 19000 -36270.879 -36270.879 -36553.434 -36553.434 273.27631 273.27631 158811.32 158811.32 -16.71501 -16.71501 Loop time of 369.391 on 1 procs for 1000 steps with 8000 atoms Performance: 0.234 ns/day, 102.609 hours/ns, 2.707 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 | 367.64 | 367.64 | 367.64 | 0.0 | 99.53 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.21303 | 0.21303 | 0.21303 | 0.0 | 0.06 Output | 0.00023494 | 0.00023494 | 0.00023494 | 0.0 | 0.00 Modify | 1.3772 | 1.3772 | 1.3772 | 0.0 | 0.37 Other | | 0.1623 | | | 0.04 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6726 ave 6726 max 6726 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 383332 ave 383332 max 383332 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 383332 Ave neighs/atom = 47.9165 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 = 273.210709382268, Press = -2.49195273684127 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 17 17 17 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.679 | 4.679 | 4.679 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 19000 -36270.879 -36270.879 -36553.434 -36553.434 273.27631 273.27631 158811.32 158811.32 -16.71501 -16.71501 20000 -36270.786 -36270.786 -36554.421 -36554.421 274.32175 274.32175 158791.89 158791.89 -80.83165 -80.83165 Loop time of 370.531 on 1 procs for 1000 steps with 8000 atoms Performance: 0.233 ns/day, 102.925 hours/ns, 2.699 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 | 368.77 | 368.77 | 368.77 | 0.0 | 99.53 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.21328 | 0.21328 | 0.21328 | 0.0 | 0.06 Output | 0.00031819 | 0.00031819 | 0.00031819 | 0.0 | 0.00 Modify | 1.3829 | 1.3829 | 1.3829 | 0.0 | 0.37 Other | | 0.1632 | | | 0.04 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6727 ave 6727 max 6727 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 383618 ave 383618 max 383618 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 383618 Ave neighs/atom = 47.95225 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 = 273.212120395523, Press = -0.263172747280541 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 17 17 17 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.679 | 4.679 | 4.679 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 20000 -36270.786 -36270.786 -36554.421 -36554.421 274.32175 274.32175 158791.89 158791.89 -80.83165 -80.83165 21000 -36271.407 -36271.407 -36559.439 -36559.439 278.5741 278.5741 158816.03 158816.03 162.95966 162.95966 Loop time of 369.158 on 1 procs for 1000 steps with 8000 atoms Performance: 0.234 ns/day, 102.544 hours/ns, 2.709 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 | 367.42 | 367.42 | 367.42 | 0.0 | 99.53 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.21183 | 0.21183 | 0.21183 | 0.0 | 0.06 Output | 0.00023559 | 0.00023559 | 0.00023559 | 0.0 | 0.00 Modify | 1.3647 | 1.3647 | 1.3647 | 0.0 | 0.37 Other | | 0.1614 | | | 0.04 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 383786 ave 383786 max 383786 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 383786 Ave neighs/atom = 47.97325 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 = 273.25102633216, Press = -1.67667411511727 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 17 17 17 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.679 | 4.679 | 4.679 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 21000 -36271.407 -36271.407 -36559.439 -36559.439 278.5741 278.5741 158816.03 158816.03 162.95966 162.95966 22000 -36270.807 -36270.807 -36555.092 -36555.092 274.95028 274.95028 158795.96 158795.96 131.77202 131.77202 Loop time of 371.009 on 1 procs for 1000 steps with 8000 atoms Performance: 0.233 ns/day, 103.058 hours/ns, 2.695 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 | 369.25 | 369.25 | 369.25 | 0.0 | 99.53 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.21397 | 0.21397 | 0.21397 | 0.0 | 0.06 Output | 0.00024225 | 0.00024225 | 0.00024225 | 0.0 | 0.00 Modify | 1.3814 | 1.3814 | 1.3814 | 0.0 | 0.37 Other | | 0.1628 | | | 0.04 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 382772 ave 382772 max 382772 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 382772 Ave neighs/atom = 47.8465 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 = 273.231070031074, Press = -0.320071840918573 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 17 17 17 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.679 | 4.679 | 4.679 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 22000 -36270.807 -36270.807 -36555.092 -36555.092 274.95028 274.95028 158795.96 158795.96 131.77202 131.77202 23000 -36268.863 -36268.863 -36551.775 -36551.775 273.62215 273.62215 158661.63 158661.63 649.59647 649.59647 Loop time of 370.055 on 1 procs for 1000 steps with 8000 atoms Performance: 0.233 ns/day, 102.793 hours/ns, 2.702 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 | 368.3 | 368.3 | 368.3 | 0.0 | 99.53 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.21265 | 0.21265 | 0.21265 | 0.0 | 0.06 Output | 0.00023295 | 0.00023295 | 0.00023295 | 0.0 | 0.00 Modify | 1.3771 | 1.3771 | 1.3771 | 0.0 | 0.37 Other | | 0.1631 | | | 0.04 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 383250 ave 383250 max 383250 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 383250 Ave neighs/atom = 47.90625 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 = 273.306855151543, Press = -1.5161011809559 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 17 17 17 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.679 | 4.679 | 4.679 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 23000 -36268.863 -36268.863 -36551.775 -36551.775 273.62215 273.62215 158661.63 158661.63 649.59647 649.59647 24000 -36267.496 -36267.496 -36555.718 -36555.718 278.75743 278.75743 158950.03 158950.03 -584.59988 -584.59988 Loop time of 370.739 on 1 procs for 1000 steps with 8000 atoms Performance: 0.233 ns/day, 102.983 hours/ns, 2.697 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 | 368.98 | 368.98 | 368.98 | 0.0 | 99.52 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.21302 | 0.21302 | 0.21302 | 0.0 | 0.06 Output | 0.00023558 | 0.00023558 | 0.00023558 | 0.0 | 0.00 Modify | 1.3837 | 1.3837 | 1.3837 | 0.0 | 0.37 Other | | 0.1647 | | | 0.04 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 384110 ave 384110 max 384110 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 384110 Ave neighs/atom = 48.01375 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 = 273.299952497581, Press = -0.484025835165537 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 17 17 17 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.679 | 4.679 | 4.679 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 24000 -36267.496 -36267.496 -36555.718 -36555.718 278.75743 278.75743 158950.03 158950.03 -584.59988 -584.59988 25000 -36268.372 -36268.372 -36551.616 -36551.616 273.94284 273.94284 158696.25 158696.25 581.26084 581.26084 Loop time of 369.215 on 1 procs for 1000 steps with 8000 atoms Performance: 0.234 ns/day, 102.560 hours/ns, 2.708 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 | 367.48 | 367.48 | 367.48 | 0.0 | 99.53 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.21223 | 0.21223 | 0.21223 | 0.0 | 0.06 Output | 0.00024257 | 0.00024257 | 0.00024257 | 0.0 | 0.00 Modify | 1.3611 | 1.3611 | 1.3611 | 0.0 | 0.37 Other | | 0.1613 | | | 0.04 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 382366 ave 382366 max 382366 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 382366 Ave neighs/atom = 47.79575 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 = 273.348938150793, Press = -0.9228851863801 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 17 17 17 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.679 | 4.679 | 4.679 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 25000 -36268.372 -36268.372 -36551.616 -36551.616 273.94284 273.94284 158696.25 158696.25 581.26084 581.26084 26000 -36270.529 -36270.529 -36554.628 -36554.628 274.76966 274.76966 159028.22 159028.22 -1006.1726 -1006.1726 Loop time of 368.262 on 1 procs for 1000 steps with 8000 atoms Performance: 0.235 ns/day, 102.295 hours/ns, 2.715 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 | 366.53 | 366.53 | 366.53 | 0.0 | 99.53 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.2119 | 0.2119 | 0.2119 | 0.0 | 0.06 Output | 0.00029415 | 0.00029415 | 0.00029415 | 0.0 | 0.00 Modify | 1.3576 | 1.3576 | 1.3576 | 0.0 | 0.37 Other | | 0.1618 | | | 0.04 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6727 ave 6727 max 6727 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 383774 ave 383774 max 383774 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 383774 Ave neighs/atom = 47.97175 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 = 273.355758268838, Press = -0.235101985530751 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 17 17 17 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.679 | 4.679 | 4.679 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 26000 -36270.529 -36270.529 -36554.628 -36554.628 274.76966 274.76966 159028.22 159028.22 -1006.1726 -1006.1726 27000 -36272.653 -36272.653 -36553.605 -36553.605 271.72563 271.72563 158820.88 158820.88 -113.95119 -113.95119 Loop time of 370.121 on 1 procs for 1000 steps with 8000 atoms Performance: 0.233 ns/day, 102.811 hours/ns, 2.702 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 | 368.37 | 368.37 | 368.37 | 0.0 | 99.53 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.21276 | 0.21276 | 0.21276 | 0.0 | 0.06 Output | 0.00023179 | 0.00023179 | 0.00023179 | 0.0 | 0.00 Modify | 1.3723 | 1.3723 | 1.3723 | 0.0 | 0.37 Other | | 0.1625 | | | 0.04 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 382414 ave 382414 max 382414 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 382414 Ave neighs/atom = 47.80175 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 = 273.399412686055, Press = -0.877288322568025 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 17 17 17 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.679 | 4.679 | 4.679 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 27000 -36272.653 -36272.653 -36553.605 -36553.605 271.72563 271.72563 158820.88 158820.88 -113.95119 -113.95119 28000 -36268.751 -36268.751 -36554.47 -36554.47 276.33718 276.33718 158768.06 158768.06 242.40373 242.40373 Loop time of 370.28 on 1 procs for 1000 steps with 8000 atoms Performance: 0.233 ns/day, 102.856 hours/ns, 2.701 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 | 368.53 | 368.53 | 368.53 | 0.0 | 99.53 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.21033 | 0.21033 | 0.21033 | 0.0 | 0.06 Output | 0.00023522 | 0.00023522 | 0.00023522 | 0.0 | 0.00 Modify | 1.3757 | 1.3757 | 1.3757 | 0.0 | 0.37 Other | | 0.1635 | | | 0.04 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 383422 ave 383422 max 383422 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 383422 Ave neighs/atom = 47.92775 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 = 273.377824594472, Press = -0.241039749836448 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 17 17 17 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.679 | 4.679 | 4.679 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 28000 -36268.751 -36268.751 -36554.47 -36554.47 276.33718 276.33718 158768.06 158768.06 242.40373 242.40373 29000 -36275.617 -36275.617 -36557.41 -36557.41 272.54045 272.54045 158855.2 158855.2 -223.98653 -223.98653 Loop time of 372.469 on 1 procs for 1000 steps with 8000 atoms Performance: 0.232 ns/day, 103.464 hours/ns, 2.685 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 | 370.7 | 370.7 | 370.7 | 0.0 | 99.53 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.21266 | 0.21266 | 0.21266 | 0.0 | 0.06 Output | 0.0002864 | 0.0002864 | 0.0002864 | 0.0 | 0.00 Modify | 1.3908 | 1.3908 | 1.3908 | 0.0 | 0.37 Other | | 0.1634 | | | 0.04 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 383286 ave 383286 max 383286 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 383286 Ave neighs/atom = 47.91075 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 = 273.33041251111, Press = -1.29656117299147 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 17 17 17 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.679 | 4.679 | 4.679 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 29000 -36275.617 -36275.617 -36557.41 -36557.41 272.54045 272.54045 158855.2 158855.2 -223.98653 -223.98653 30000 -36269.702 -36269.702 -36554.589 -36554.589 275.53232 275.53232 158707.88 158707.88 508.8522 508.8522 Loop time of 374.709 on 1 procs for 1000 steps with 8000 atoms Performance: 0.231 ns/day, 104.086 hours/ns, 2.669 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 | 372.92 | 372.92 | 372.92 | 0.0 | 99.52 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.21335 | 0.21335 | 0.21335 | 0.0 | 0.06 Output | 0.00023277 | 0.00023277 | 0.00023277 | 0.0 | 0.00 Modify | 1.4102 | 1.4102 | 1.4102 | 0.0 | 0.38 Other | | 0.1656 | | | 0.04 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 382966 ave 382966 max 382966 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 382966 Ave neighs/atom = 47.87075 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 = 273.290985888865, Press = -0.263109359256337 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 17 17 17 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.679 | 4.679 | 4.679 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 30000 -36269.702 -36269.702 -36554.589 -36554.589 275.53232 275.53232 158707.88 158707.88 508.8522 508.8522 31000 -36272.563 -36272.563 -36554.696 -36554.696 272.86844 272.86844 159014.67 159014.67 -1116.7366 -1116.7366 Loop time of 372.105 on 1 procs for 1000 steps with 8000 atoms Performance: 0.232 ns/day, 103.363 hours/ns, 2.687 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 | 370.34 | 370.34 | 370.34 | 0.0 | 99.53 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.21013 | 0.21013 | 0.21013 | 0.0 | 0.06 Output | 0.00023375 | 0.00023375 | 0.00023375 | 0.0 | 0.00 Modify | 1.3891 | 1.3891 | 1.3891 | 0.0 | 0.37 Other | | 0.1641 | | | 0.04 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 383812 ave 383812 max 383812 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 383812 Ave neighs/atom = 47.9765 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 = 273.266065818342, Press = -0.852767693044262 next a jump SELF top variable a loop 2000 run 1000 Generated 0 of 0 mixed pair_coeff terms from geometric mixing rule Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 17 17 17 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6.5 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.679 | 4.679 | 4.679 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 31000 -36272.563 -36272.563 -36554.696 -36554.696 272.86844 272.86844 159014.67 159014.67 -1116.7366 -1116.7366 32000 -36271.411 -36271.411 -36556.25 -36556.25 275.48585 275.48585 158657.11 158657.11 817.70711 817.70711 Loop time of 372.376 on 1 procs for 1000 steps with 8000 atoms Performance: 0.232 ns/day, 103.438 hours/ns, 2.685 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 | 370.61 | 370.61 | 370.61 | 0.0 | 99.52 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.21149 | 0.21149 | 0.21149 | 0.0 | 0.06 Output | 0.00023372 | 0.00023372 | 0.00023372 | 0.0 | 0.00 Modify | 1.3938 | 1.3938 | 1.3938 | 0.0 | 0.37 Other | | 0.1648 | | | 0.04 Nlocal: 8000 ave 8000 max 8000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 6725 ave 6725 max 6725 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 382752 ave 382752 max 382752 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 382752 Ave neighs/atom = 47.844 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_T273.15.out" else "print 'not_converged' file output/vol_T273.15.out" print '${V}' file output/vol_T273.15.out 158802.001992138 print "LAMMPS calculation completed" LAMMPS calculation completed quit 0