# 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 3.521391771733761*${_u_distance} variable latticeconst_converted equal 3.521391771733761*1 lattice fcc ${latticeconst_converted} lattice fcc 3.52139177173376 Lattice spacing in x,y,z = 3.5213918 3.5213918 3.5213918 region simbox block 0 10 0 10 0 10 units lattice create_box 1 simbox Created orthogonal box = (0 0 0) to (35.213918 35.213918 35.213918) 1 by 1 by 1 MPI processor grid create_atoms 1 box Created 4000 atoms using lattice units in orthogonal box = (0 0 0) to (35.213918 35.213918 35.213918) create_atoms CPU = 0.001 seconds variable mass_converted equal 58.6934*${_u_mass} variable mass_converted equal 58.6934*1 kim_interactions Ni WARNING: 'kim_' has been renamed to 'kim '. Please update your input. kim interactions Ni #=== BEGIN kim interactions ================================== pair_style kim MEAM_LAMMPS_ChoiKimSeol_2017_NiMn__MO_348689608050_002 pair_coeff * * Ni #=== END kim interactions ==================================== mass 1 ${mass_converted} mass 1 58.6934 # initial volume variable v equal vol # assign formula variable V0 equal ${v} # evaluate initial value variable V0 equal 43665.9622831874 variable V0_metal equal ${V0}/(${_u_distance}*${_u_distance}*${_u_distance}) variable V0_metal equal 43665.9622831874/(${_u_distance}*${_u_distance}*${_u_distance}) variable V0_metal equal 43665.9622831874/(1*${_u_distance}*${_u_distance}) variable V0_metal equal 43665.9622831874/(1*1*${_u_distance}) variable V0_metal equal 43665.9622831874/(1*1*1) variable V0_metal_times1000 equal ${V0_metal}*1000 variable V0_metal_times1000 equal 43665.9622831874*1000 print "Initial system volume: ${V0_metal} Angstroms^3" Initial system volume: 43665.9622831874 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 253.15*${_u_temperature} variable temp_converted equal 253.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 253.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 253.15 ${temp_converted} ${Tdamp_converted} iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 253.15 253.15 ${Tdamp_converted} iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 253.15 253.15 0.1 iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 253.15 253.15 0.1 iso 0 ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 253.15 253.15 0.1 iso 0 0 ${Pdamp_converted} fix ensemble all npt temp 253.15 253.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 "253.15 - 0.2" variable T_up equal "253.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_348689608050_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 ghost atom cutoff = 6 binsize = 3, bins = 12 12 12 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.61 | 4.61 | 4.61 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 0 -17669.144 -17669.144 -17800 -17800 253.15 253.15 43665.962 43665.962 3200.8637 3200.8637 1000 -17530.093 -17530.093 -17668.356 -17668.356 267.47983 267.47983 44119.697 44119.697 -2026.3015 -2026.3015 Loop time of 49.5601 on 1 procs for 1000 steps with 4000 atoms Performance: 1.743 ns/day, 13.767 hours/ns, 20.178 timesteps/s 99.6% 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.409 | 49.409 | 49.409 | 0.0 | 99.70 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.021759 | 0.021759 | 0.021759 | 0.0 | 0.04 Output | 0.00013947 | 0.00013947 | 0.00013947 | 0.0 | 0.00 Modify | 0.11844 | 0.11844 | 0.11844 | 0.0 | 0.24 Other | | 0.01058 | | | 0.02 Nlocal: 4000 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841 ave 5841 max 5841 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: 312000 ave 312000 max 312000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 312000 Ave neighs/atom = 78 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 ghost atom cutoff = 6 binsize = 3, bins = 12 12 12 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.61 | 4.61 | 4.61 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 1000 -17530.093 -17530.093 -17668.356 -17668.356 267.47983 267.47983 44119.697 44119.697 -2026.3015 -2026.3015 2000 -17537.851 -17537.851 -17665.863 -17665.863 247.64814 247.64814 44046.438 44046.438 1072.7777 1072.7777 Loop time of 56.869 on 1 procs for 1000 steps with 4000 atoms Performance: 1.519 ns/day, 15.797 hours/ns, 17.584 timesteps/s 99.7% 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.72 | 56.72 | 56.72 | 0.0 | 99.74 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.021683 | 0.021683 | 0.021683 | 0.0 | 0.04 Output | 4.9784e-05 | 4.9784e-05 | 4.9784e-05 | 0.0 | 0.00 Modify | 0.11733 | 0.11733 | 0.11733 | 0.0 | 0.21 Other | | 0.01025 | | | 0.02 Nlocal: 4000 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841 ave 5841 max 5841 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: 313488 ave 313488 max 313488 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 313488 Ave neighs/atom = 78.372 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 ghost atom cutoff = 6 binsize = 3, bins = 12 12 12 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.61 | 4.61 | 4.61 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 2000 -17537.851 -17537.851 -17665.863 -17665.863 247.64814 247.64814 44046.438 44046.438 1072.7777 1072.7777 3000 -17537.617 -17537.617 -17668.667 -17668.667 253.52538 253.52538 44068.269 44068.269 -111.50552 -111.50552 Loop time of 57.2229 on 1 procs for 1000 steps with 4000 atoms Performance: 1.510 ns/day, 15.895 hours/ns, 17.476 timesteps/s 99.7% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 57.073 | 57.073 | 57.073 | 0.0 | 99.74 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.021621 | 0.021621 | 0.021621 | 0.0 | 0.04 Output | 5.0415e-05 | 5.0415e-05 | 5.0415e-05 | 0.0 | 0.00 Modify | 0.11749 | 0.11749 | 0.11749 | 0.0 | 0.21 Other | | 0.0103 | | | 0.02 Nlocal: 4000 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841 ave 5841 max 5841 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: 313800 ave 313800 max 313800 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 313800 Ave neighs/atom = 78.45 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 ghost atom cutoff = 6 binsize = 3, bins = 12 12 12 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.61 | 4.61 | 4.61 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 3000 -17537.617 -17537.617 -17668.667 -17668.667 253.52538 253.52538 44068.269 44068.269 -111.50552 -111.50552 4000 -17535.019 -17535.019 -17664.216 -17664.216 249.94055 249.94055 44083.053 44083.053 -268.1606 -268.1606 Loop time of 57.3076 on 1 procs for 1000 steps with 4000 atoms Performance: 1.508 ns/day, 15.919 hours/ns, 17.450 timesteps/s 99.7% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 57.158 | 57.158 | 57.158 | 0.0 | 99.74 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.021658 | 0.021658 | 0.021658 | 0.0 | 0.04 Output | 5.3752e-05 | 5.3752e-05 | 5.3752e-05 | 0.0 | 0.00 Modify | 0.11739 | 0.11739 | 0.11739 | 0.0 | 0.20 Other | | 0.01033 | | | 0.02 Nlocal: 4000 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841 ave 5841 max 5841 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: 313746 ave 313746 max 313746 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 313746 Ave neighs/atom = 78.4365 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 ghost atom cutoff = 6 binsize = 3, bins = 12 12 12 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.61 | 4.61 | 4.61 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 4000 -17535.019 -17535.019 -17664.216 -17664.216 249.94055 249.94055 44083.053 44083.053 -268.1606 -268.1606 5000 -17538.817 -17538.817 -17671.232 -17671.232 256.16704 256.16704 44032.977 44032.977 1228.552 1228.552 Loop time of 57.2194 on 1 procs for 1000 steps with 4000 atoms Performance: 1.510 ns/day, 15.894 hours/ns, 17.477 timesteps/s 99.7% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 57.07 | 57.07 | 57.07 | 0.0 | 99.74 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.021638 | 0.021638 | 0.021638 | 0.0 | 0.04 Output | 3.9395e-05 | 3.9395e-05 | 3.9395e-05 | 0.0 | 0.00 Modify | 0.11771 | 0.11771 | 0.11771 | 0.0 | 0.21 Other | | 0.01031 | | | 0.02 Nlocal: 4000 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841 ave 5841 max 5841 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: 313686 ave 313686 max 313686 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 313686 Ave neighs/atom = 78.4215 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 = 250.9134755731, Press = -532.831960493775 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 ghost atom cutoff = 6 binsize = 3, bins = 12 12 12 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.61 | 4.61 | 4.61 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 5000 -17538.817 -17538.817 -17671.232 -17671.232 256.16704 256.16704 44032.977 44032.977 1228.552 1228.552 6000 -17534.041 -17534.041 -17664.645 -17664.645 252.66261 252.66261 44110.352 44110.352 -1397.0525 -1397.0525 Loop time of 57.2888 on 1 procs for 1000 steps with 4000 atoms Performance: 1.508 ns/day, 15.914 hours/ns, 17.455 timesteps/s 99.7% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 57.135 | 57.135 | 57.135 | 0.0 | 99.73 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.021719 | 0.021719 | 0.021719 | 0.0 | 0.04 Output | 4.0286e-05 | 4.0286e-05 | 4.0286e-05 | 0.0 | 0.00 Modify | 0.12126 | 0.12126 | 0.12126 | 0.0 | 0.21 Other | | 0.01036 | | | 0.02 Nlocal: 4000 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841 ave 5841 max 5841 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: 313686 ave 313686 max 313686 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 313686 Ave neighs/atom = 78.4215 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 = 252.808174911893, Press = -42.2789794010795 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 ghost atom cutoff = 6 binsize = 3, bins = 12 12 12 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.61 | 4.61 | 4.61 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 6000 -17534.041 -17534.041 -17664.645 -17664.645 252.66261 252.66261 44110.352 44110.352 -1397.0525 -1397.0525 7000 -17538.784 -17538.784 -17668.817 -17668.817 251.55639 251.55639 44033.259 44033.259 1429.1677 1429.1677 Loop time of 57.2955 on 1 procs for 1000 steps with 4000 atoms Performance: 1.508 ns/day, 15.915 hours/ns, 17.453 timesteps/s 99.7% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 57.142 | 57.142 | 57.142 | 0.0 | 99.73 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.021828 | 0.021828 | 0.021828 | 0.0 | 0.04 Output | 3.8863e-05 | 3.8863e-05 | 3.8863e-05 | 0.0 | 0.00 Modify | 0.1213 | 0.1213 | 0.1213 | 0.0 | 0.21 Other | | 0.01034 | | | 0.02 Nlocal: 4000 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841 ave 5841 max 5841 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: 313710 ave 313710 max 313710 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 313710 Ave neighs/atom = 78.4275 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 = 253.02878829138, Press = -5.41002935989565 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 ghost atom cutoff = 6 binsize = 3, bins = 12 12 12 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.61 | 4.61 | 4.61 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 7000 -17538.784 -17538.784 -17668.817 -17668.817 251.55639 251.55639 44033.259 44033.259 1429.1677 1429.1677 8000 -17536.264 -17536.264 -17670.456 -17670.456 259.60338 259.60338 44076.851 44076.851 -566.33379 -566.33379 Loop time of 57.3251 on 1 procs for 1000 steps with 4000 atoms Performance: 1.507 ns/day, 15.924 hours/ns, 17.444 timesteps/s 99.7% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 57.172 | 57.172 | 57.172 | 0.0 | 99.73 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.021739 | 0.021739 | 0.021739 | 0.0 | 0.04 Output | 3.8603e-05 | 3.8603e-05 | 3.8603e-05 | 0.0 | 0.00 Modify | 0.1212 | 0.1212 | 0.1212 | 0.0 | 0.21 Other | | 0.01037 | | | 0.02 Nlocal: 4000 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841 ave 5841 max 5841 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: 313740 ave 313740 max 313740 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 313740 Ave neighs/atom = 78.435 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 = 253.00707729478, Press = -29.3109557003675 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 ghost atom cutoff = 6 binsize = 3, bins = 12 12 12 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.61 | 4.61 | 4.61 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 8000 -17536.264 -17536.264 -17670.456 -17670.456 259.60338 259.60338 44076.851 44076.851 -566.33379 -566.33379 9000 -17537.79 -17537.79 -17668.187 -17668.187 252.26117 252.26117 44093.367 44093.367 -1040.9415 -1040.9415 Loop time of 57.2765 on 1 procs for 1000 steps with 4000 atoms Performance: 1.508 ns/day, 15.910 hours/ns, 17.459 timesteps/s 99.7% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 57.123 | 57.123 | 57.123 | 0.0 | 99.73 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.021638 | 0.021638 | 0.021638 | 0.0 | 0.04 Output | 3.9354e-05 | 3.9354e-05 | 3.9354e-05 | 0.0 | 0.00 Modify | 0.12125 | 0.12125 | 0.12125 | 0.0 | 0.21 Other | | 0.01034 | | | 0.02 Nlocal: 4000 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841 ave 5841 max 5841 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: 313558 ave 313558 max 313558 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 313558 Ave neighs/atom = 78.3895 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 = 252.578876359855, Press = -2.44059132063483 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 ghost atom cutoff = 6 binsize = 3, bins = 12 12 12 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.61 | 4.61 | 4.61 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 9000 -17537.79 -17537.79 -17668.187 -17668.187 252.26117 252.26117 44093.367 44093.367 -1040.9415 -1040.9415 10000 -17536.061 -17536.061 -17667.104 -17667.104 253.51173 253.51173 44030.477 44030.477 1715.3786 1715.3786 Loop time of 57.2999 on 1 procs for 1000 steps with 4000 atoms Performance: 1.508 ns/day, 15.917 hours/ns, 17.452 timesteps/s 99.7% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 57.147 | 57.147 | 57.147 | 0.0 | 99.73 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.021652 | 0.021652 | 0.021652 | 0.0 | 0.04 Output | 3.8733e-05 | 3.8733e-05 | 3.8733e-05 | 0.0 | 0.00 Modify | 0.12112 | 0.12112 | 0.12112 | 0.0 | 0.21 Other | | 0.01035 | | | 0.02 Nlocal: 4000 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841 ave 5841 max 5841 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: 313562 ave 313562 max 313562 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 313562 Ave neighs/atom = 78.3905 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 = 252.499440288546, Press = -9.90854631825337 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 ghost atom cutoff = 6 binsize = 3, bins = 12 12 12 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.61 | 4.61 | 4.61 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 10000 -17536.061 -17536.061 -17667.104 -17667.104 253.51173 253.51173 44030.477 44030.477 1715.3786 1715.3786 11000 -17536.499 -17536.499 -17665.641 -17665.641 249.83424 249.83424 44110.879 44110.879 -1599.8296 -1599.8296 Loop time of 57.3213 on 1 procs for 1000 steps with 4000 atoms Performance: 1.507 ns/day, 15.923 hours/ns, 17.446 timesteps/s 99.7% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 57.168 | 57.168 | 57.168 | 0.0 | 99.73 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.021612 | 0.021612 | 0.021612 | 0.0 | 0.04 Output | 3.8202e-05 | 3.8202e-05 | 3.8202e-05 | 0.0 | 0.00 Modify | 0.12121 | 0.12121 | 0.12121 | 0.0 | 0.21 Other | | 0.01032 | | | 0.02 Nlocal: 4000 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841 ave 5841 max 5841 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: 313754 ave 313754 max 313754 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 313754 Ave neighs/atom = 78.4385 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 = 252.498985733799, Press = -9.4414617755967 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 ghost atom cutoff = 6 binsize = 3, bins = 12 12 12 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.61 | 4.61 | 4.61 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 11000 -17536.499 -17536.499 -17665.641 -17665.641 249.83424 249.83424 44110.879 44110.879 -1599.8296 -1599.8296 12000 -17542.273 -17542.273 -17669.749 -17669.749 246.61147 246.61147 44038.672 44038.672 1017.7026 1017.7026 Loop time of 57.3267 on 1 procs for 1000 steps with 4000 atoms Performance: 1.507 ns/day, 15.924 hours/ns, 17.444 timesteps/s 99.7% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 57.173 | 57.173 | 57.173 | 0.0 | 99.73 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.021791 | 0.021791 | 0.021791 | 0.0 | 0.04 Output | 3.7962e-05 | 3.7962e-05 | 3.7962e-05 | 0.0 | 0.00 Modify | 0.12121 | 0.12121 | 0.12121 | 0.0 | 0.21 Other | | 0.01037 | | | 0.02 Nlocal: 4000 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841 ave 5841 max 5841 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: 313640 ave 313640 max 313640 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 313640 Ave neighs/atom = 78.41 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 = 252.509911343097, Press = -1.18160366865634 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 ghost atom cutoff = 6 binsize = 3, bins = 12 12 12 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.61 | 4.61 | 4.61 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 12000 -17542.273 -17542.273 -17669.749 -17669.749 246.61147 246.61147 44038.672 44038.672 1017.7026 1017.7026 13000 -17535.199 -17535.199 -17668.989 -17668.989 258.82456 258.82456 44075.227 44075.227 -283.59325 -283.59325 Loop time of 57.3506 on 1 procs for 1000 steps with 4000 atoms Performance: 1.507 ns/day, 15.931 hours/ns, 17.437 timesteps/s 99.7% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 57.197 | 57.197 | 57.197 | 0.0 | 99.73 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.021632 | 0.021632 | 0.021632 | 0.0 | 0.04 Output | 2.9586e-05 | 2.9586e-05 | 2.9586e-05 | 0.0 | 0.00 Modify | 0.12118 | 0.12118 | 0.12118 | 0.0 | 0.21 Other | | 0.01034 | | | 0.02 Nlocal: 4000 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841 ave 5841 max 5841 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: 313660 ave 313660 max 313660 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 313660 Ave neighs/atom = 78.415 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 = 252.47056091564, Press = -8.51146623195887 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 ghost atom cutoff = 6 binsize = 3, bins = 12 12 12 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.61 | 4.61 | 4.61 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 13000 -17535.199 -17535.199 -17668.989 -17668.989 258.82456 258.82456 44075.227 44075.227 -283.59325 -283.59325 14000 -17538.312 -17538.312 -17670.219 -17670.219 255.18247 255.18247 44077.559 44077.559 -569.34288 -569.34288 Loop time of 57.2806 on 1 procs for 1000 steps with 4000 atoms Performance: 1.508 ns/day, 15.911 hours/ns, 17.458 timesteps/s 99.7% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 57.127 | 57.127 | 57.127 | 0.0 | 99.73 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.021639 | 0.021639 | 0.021639 | 0.0 | 0.04 Output | 7.3209e-05 | 7.3209e-05 | 7.3209e-05 | 0.0 | 0.00 Modify | 0.12118 | 0.12118 | 0.12118 | 0.0 | 0.21 Other | | 0.01035 | | | 0.02 Nlocal: 4000 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841 ave 5841 max 5841 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: 313670 ave 313670 max 313670 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 313670 Ave neighs/atom = 78.4175 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 = 252.42449079862, Press = -2.9867514831822 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 ghost atom cutoff = 6 binsize = 3, bins = 12 12 12 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.61 | 4.61 | 4.61 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 14000 -17538.312 -17538.312 -17670.219 -17670.219 255.18247 255.18247 44077.559 44077.559 -569.34288 -569.34288 15000 -17536.276 -17536.276 -17667.159 -17667.159 253.20155 253.20155 44044.501 44044.501 1118.6824 1118.6824 Loop time of 57.3065 on 1 procs for 1000 steps with 4000 atoms Performance: 1.508 ns/day, 15.918 hours/ns, 17.450 timesteps/s 99.7% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 57.153 | 57.153 | 57.153 | 0.0 | 99.73 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.02171 | 0.02171 | 0.02171 | 0.0 | 0.04 Output | 4.0567e-05 | 4.0567e-05 | 4.0567e-05 | 0.0 | 0.00 Modify | 0.12124 | 0.12124 | 0.12124 | 0.0 | 0.21 Other | | 0.01035 | | | 0.02 Nlocal: 4000 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841 ave 5841 max 5841 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: 313604 ave 313604 max 313604 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 313604 Ave neighs/atom = 78.401 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 = 252.350525362619, Press = -4.91404627893523 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 ghost atom cutoff = 6 binsize = 3, bins = 12 12 12 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.61 | 4.61 | 4.61 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 15000 -17536.276 -17536.276 -17667.159 -17667.159 253.20155 253.20155 44044.501 44044.501 1118.6824 1118.6824 16000 -17539.643 -17539.643 -17670.712 -17670.712 253.56072 253.56072 44129.337 44129.337 -2829.475 -2829.475 Loop time of 57.3982 on 1 procs for 1000 steps with 4000 atoms Performance: 1.505 ns/day, 15.944 hours/ns, 17.422 timesteps/s 99.7% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 57.245 | 57.245 | 57.245 | 0.0 | 99.73 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.021676 | 0.021676 | 0.021676 | 0.0 | 0.04 Output | 4.0476e-05 | 4.0476e-05 | 4.0476e-05 | 0.0 | 0.00 Modify | 0.12129 | 0.12129 | 0.12129 | 0.0 | 0.21 Other | | 0.01035 | | | 0.02 Nlocal: 4000 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841 ave 5841 max 5841 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: 313702 ave 313702 max 313702 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 313702 Ave neighs/atom = 78.4255 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 = 252.423825860421, Press = -5.40613906038276 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 ghost atom cutoff = 6 binsize = 3, bins = 12 12 12 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.61 | 4.61 | 4.61 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 16000 -17539.643 -17539.643 -17670.712 -17670.712 253.56072 253.56072 44129.337 44129.337 -2829.475 -2829.475 17000 -17535.412 -17535.412 -17668.399 -17668.399 257.27226 257.27226 44020.054 44020.054 2009.7248 2009.7248 Loop time of 57.3149 on 1 procs for 1000 steps with 4000 atoms Performance: 1.507 ns/day, 15.921 hours/ns, 17.447 timesteps/s 99.7% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 57.162 | 57.162 | 57.162 | 0.0 | 99.73 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.021649 | 0.021649 | 0.021649 | 0.0 | 0.04 Output | 4.0076e-05 | 4.0076e-05 | 4.0076e-05 | 0.0 | 0.00 Modify | 0.12129 | 0.12129 | 0.12129 | 0.0 | 0.21 Other | | 0.01033 | | | 0.02 Nlocal: 4000 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841 ave 5841 max 5841 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: 313414 ave 313414 max 313414 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 313414 Ave neighs/atom = 78.3535 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 = 252.572458790754, Press = -0.695403992221019 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 ghost atom cutoff = 6 binsize = 3, bins = 12 12 12 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.61 | 4.61 | 4.61 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 17000 -17535.412 -17535.412 -17668.399 -17668.399 257.27226 257.27226 44020.054 44020.054 2009.7248 2009.7248 18000 -17534.305 -17534.305 -17664.059 -17664.059 251.01856 251.01856 44084.38 44084.38 -348.27103 -348.27103 Loop time of 57.5302 on 1 procs for 1000 steps with 4000 atoms Performance: 1.502 ns/day, 15.981 hours/ns, 17.382 timesteps/s 99.7% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 57.377 | 57.377 | 57.377 | 0.0 | 99.73 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.021691 | 0.021691 | 0.021691 | 0.0 | 0.04 Output | 4.1208e-05 | 4.1208e-05 | 4.1208e-05 | 0.0 | 0.00 Modify | 0.12148 | 0.12148 | 0.12148 | 0.0 | 0.21 Other | | 0.01036 | | | 0.02 Nlocal: 4000 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841 ave 5841 max 5841 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: 313828 ave 313828 max 313828 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 313828 Ave neighs/atom = 78.457 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 = 252.662467864716, Press = -5.25893813970224 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 ghost atom cutoff = 6 binsize = 3, bins = 12 12 12 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.61 | 4.61 | 4.61 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 18000 -17534.305 -17534.305 -17664.059 -17664.059 251.01856 251.01856 44084.38 44084.38 -348.27103 -348.27103 19000 -17539.278 -17539.278 -17669.886 -17669.886 252.66961 252.66961 44064.003 44064.003 -49.778569 -49.778569 Loop time of 57.3945 on 1 procs for 1000 steps with 4000 atoms Performance: 1.505 ns/day, 15.943 hours/ns, 17.423 timesteps/s 99.7% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 57.241 | 57.241 | 57.241 | 0.0 | 99.73 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.021655 | 0.021655 | 0.021655 | 0.0 | 0.04 Output | 3.9084e-05 | 3.9084e-05 | 3.9084e-05 | 0.0 | 0.00 Modify | 0.12123 | 0.12123 | 0.12123 | 0.0 | 0.21 Other | | 0.01034 | | | 0.02 Nlocal: 4000 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841 ave 5841 max 5841 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: 313708 ave 313708 max 313708 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 313708 Ave neighs/atom = 78.427 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 = 252.787463824164, Press = -1.77317437042269 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 ghost atom cutoff = 6 binsize = 3, bins = 12 12 12 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.61 | 4.61 | 4.61 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 19000 -17539.278 -17539.278 -17669.886 -17669.886 252.66961 252.66961 44064.003 44064.003 -49.778569 -49.778569 20000 -17540.112 -17540.112 -17668.388 -17668.388 248.15844 248.15844 44061.611 44061.611 182.2083 182.2083 Loop time of 57.3306 on 1 procs for 1000 steps with 4000 atoms Performance: 1.507 ns/day, 15.925 hours/ns, 17.443 timesteps/s 99.7% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 57.178 | 57.178 | 57.178 | 0.0 | 99.73 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.021466 | 0.021466 | 0.021466 | 0.0 | 0.04 Output | 3.9695e-05 | 3.9695e-05 | 3.9695e-05 | 0.0 | 0.00 Modify | 0.12102 | 0.12102 | 0.12102 | 0.0 | 0.21 Other | | 0.01028 | | | 0.02 Nlocal: 4000 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841 ave 5841 max 5841 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: 313656 ave 313656 max 313656 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 313656 Ave neighs/atom = 78.414 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 = 252.807657896009, Press = -2.72314884502105 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 ghost atom cutoff = 6 binsize = 3, bins = 12 12 12 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.61 | 4.61 | 4.61 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 20000 -17540.112 -17540.112 -17668.388 -17668.388 248.15844 248.15844 44061.611 44061.611 182.2083 182.2083 21000 -17538.442 -17538.442 -17667.905 -17667.905 250.45535 250.45535 44083.015 44083.015 -611.87885 -611.87885 Loop time of 57.4349 on 1 procs for 1000 steps with 4000 atoms Performance: 1.504 ns/day, 15.954 hours/ns, 17.411 timesteps/s 99.7% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 57.282 | 57.282 | 57.282 | 0.0 | 99.73 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.021633 | 0.021633 | 0.021633 | 0.0 | 0.04 Output | 3.698e-05 | 3.698e-05 | 3.698e-05 | 0.0 | 0.00 Modify | 0.12105 | 0.12105 | 0.12105 | 0.0 | 0.21 Other | | 0.01015 | | | 0.02 Nlocal: 4000 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841 ave 5841 max 5841 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: 313646 ave 313646 max 313646 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 313646 Ave neighs/atom = 78.4115 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 = 252.847672587776, Press = -2.76030629381302 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 ghost atom cutoff = 6 binsize = 3, bins = 12 12 12 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.61 | 4.61 | 4.61 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 21000 -17538.442 -17538.442 -17667.905 -17667.905 250.45535 250.45535 44083.015 44083.015 -611.87885 -611.87885 22000 -17533.212 -17533.212 -17667.193 -17667.193 259.19617 259.19617 44061.241 44061.241 432.60885 432.60885 Loop time of 57.4405 on 1 procs for 1000 steps with 4000 atoms Performance: 1.504 ns/day, 15.956 hours/ns, 17.409 timesteps/s 99.7% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 57.288 | 57.288 | 57.288 | 0.0 | 99.73 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.021549 | 0.021549 | 0.021549 | 0.0 | 0.04 Output | 3.7611e-05 | 3.7611e-05 | 3.7611e-05 | 0.0 | 0.00 Modify | 0.12124 | 0.12124 | 0.12124 | 0.0 | 0.21 Other | | 0.01007 | | | 0.02 Nlocal: 4000 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841 ave 5841 max 5841 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: 313602 ave 313602 max 313602 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 313602 Ave neighs/atom = 78.4005 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 = 252.915793062524, Press = -2.3844973865662 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 ghost atom cutoff = 6 binsize = 3, bins = 12 12 12 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.61 | 4.61 | 4.61 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 22000 -17533.212 -17533.212 -17667.193 -17667.193 259.19617 259.19617 44061.241 44061.241 432.60885 432.60885 23000 -17538.709 -17538.709 -17666.609 -17666.609 247.43136 247.43136 44093.604 44093.604 -1000.4281 -1000.4281 Loop time of 57.4142 on 1 procs for 1000 steps with 4000 atoms Performance: 1.505 ns/day, 15.948 hours/ns, 17.417 timesteps/s 99.7% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 57.261 | 57.261 | 57.261 | 0.0 | 99.73 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.021638 | 0.021638 | 0.021638 | 0.0 | 0.04 Output | 3.723e-05 | 3.723e-05 | 3.723e-05 | 0.0 | 0.00 Modify | 0.12111 | 0.12111 | 0.12111 | 0.0 | 0.21 Other | | 0.0101 | | | 0.02 Nlocal: 4000 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841 ave 5841 max 5841 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: 313716 ave 313716 max 313716 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 313716 Ave neighs/atom = 78.429 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 = 252.986654772505, Press = -3.35269724801844 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 ghost atom cutoff = 6 binsize = 3, bins = 12 12 12 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 6 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.61 | 4.61 | 4.61 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 23000 -17538.709 -17538.709 -17666.609 -17666.609 247.43136 247.43136 44093.604 44093.604 -1000.4281 -1000.4281 24000 -17536.972 -17536.972 -17667.519 -17667.519 252.55243 252.55243 44047.518 44047.518 888.00658 888.00658 Loop time of 57.3069 on 1 procs for 1000 steps with 4000 atoms Performance: 1.508 ns/day, 15.919 hours/ns, 17.450 timesteps/s 99.7% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 57.154 | 57.154 | 57.154 | 0.0 | 99.73 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.021564 | 0.021564 | 0.021564 | 0.0 | 0.04 Output | 3.8282e-05 | 3.8282e-05 | 3.8282e-05 | 0.0 | 0.00 Modify | 0.12123 | 0.12123 | 0.12123 | 0.0 | 0.21 Other | | 0.01011 | | | 0.02 Nlocal: 4000 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841 ave 5841 max 5841 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: 313636 ave 313636 max 313636 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 313636 Ave neighs/atom = 78.409 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_T253.15.out" else "print 'not_converged' file output/vol_T253.15.out" print '${V}' file output/vol_T253.15.out 44068.7708480297 print "LAMMPS calculation completed" LAMMPS calculation completed quit 0