# 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.917371571063995*${_u_distance} variable latticeconst_converted equal 3.917371571063995*1 lattice fcc ${latticeconst_converted} lattice fcc 3.91737157106399 Lattice spacing in x,y,z = 3.9173716 3.9173716 3.9173716 region simbox block 0 10 0 10 0 10 units lattice create_box 1 simbox Created orthogonal box = (0 0 0) to (39.173716 39.173716 39.173716) 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 (39.173716 39.173716 39.173716) create_atoms CPU = 0.004 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_KimSeolJi_2017_PtNi__MO_020840179467_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 60115.200756225 variable V0_metal equal ${V0}/(${_u_distance}*${_u_distance}*${_u_distance}) variable V0_metal equal 60115.200756225/(${_u_distance}*${_u_distance}*${_u_distance}) variable V0_metal equal 60115.200756225/(1*${_u_distance}*${_u_distance}) variable V0_metal equal 60115.200756225/(1*1*${_u_distance}) variable V0_metal equal 60115.200756225/(1*1*1) variable V0_metal_times1000 equal ${V0_metal}*1000 variable V0_metal_times1000 equal 60115.200756225*1000 print "Initial system volume: ${V0_metal} Angstroms^3" Initial system volume: 60115.200756225 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 333.15*${_u_temperature} variable temp_converted equal 333.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 333.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 333.15 ${temp_converted} ${Tdamp_converted} iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 333.15 333.15 ${Tdamp_converted} iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 333.15 333.15 0.1 iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 333.15 333.15 0.1 iso 0 ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 333.15 333.15 0.1 iso 0 0 ${Pdamp_converted} fix ensemble all npt temp 333.15 333.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 "333.15 - 0.2" variable T_up equal "333.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_020840179467_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 = 13 13 13 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.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 0 -22907.791 -22907.791 -23080 -23080 333.15 333.15 60115.201 60115.201 3059.7675 3059.7675 1000 -22727.237 -22727.237 -22905.516 -22905.516 344.89383 344.89383 60741.06 60741.06 -1796.6215 -1796.6215 Loop time of 189.324 on 1 procs for 1000 steps with 4000 atoms Performance: 0.456 ns/day, 52.590 hours/ns, 5.282 timesteps/s 100.0% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 188.34 | 188.34 | 188.34 | 0.0 | 99.48 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.15746 | 0.15746 | 0.15746 | 0.0 | 0.08 Output | 0.0002391 | 0.0002391 | 0.0002391 | 0.0 | 0.00 Modify | 0.72939 | 0.72939 | 0.72939 | 0.0 | 0.39 Other | | 0.09997 | | | 0.05 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.5 ghost atom cutoff = 6.5 binsize = 3.25, bins = 13 13 13 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.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 1000 -22727.237 -22727.237 -22905.516 -22905.516 344.89383 344.89383 60741.06 60741.06 -1796.6215 -1796.6215 2000 -22735.815 -22735.815 -22907.923 -22907.923 332.95456 332.95456 60648.314 60648.314 515.99209 515.99209 Loop time of 188.258 on 1 procs for 1000 steps with 4000 atoms Performance: 0.459 ns/day, 52.294 hours/ns, 5.312 timesteps/s 100.0% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 187.35 | 187.35 | 187.35 | 0.0 | 99.52 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.14784 | 0.14784 | 0.14784 | 0.0 | 0.08 Output | 0.00028442 | 0.00028442 | 0.00028442 | 0.0 | 0.00 Modify | 0.65847 | 0.65847 | 0.65847 | 0.0 | 0.35 Other | | 0.09664 | | | 0.05 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: 311968 ave 311968 max 311968 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311968 Ave neighs/atom = 77.992 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 = 13 13 13 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.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 2000 -22735.815 -22735.815 -22907.923 -22907.923 332.95456 332.95456 60648.314 60648.314 515.99209 515.99209 3000 -22737.168 -22737.168 -22910.598 -22910.598 335.51325 335.51325 60663.091 60663.091 -116.91973 -116.91973 Loop time of 219.118 on 1 procs for 1000 steps with 4000 atoms Performance: 0.394 ns/day, 60.866 hours/ns, 4.564 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 | 218.04 | 218.04 | 218.04 | 0.0 | 99.51 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16636 | 0.16636 | 0.16636 | 0.0 | 0.08 Output | 0.00028786 | 0.00028786 | 0.00028786 | 0.0 | 0.00 Modify | 0.80566 | 0.80566 | 0.80566 | 0.0 | 0.37 Other | | 0.1087 | | | 0.05 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: 311950 ave 311950 max 311950 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311950 Ave neighs/atom = 77.9875 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 = 13 13 13 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.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 3000 -22737.168 -22737.168 -22910.598 -22910.598 335.51325 335.51325 60663.091 60663.091 -116.91973 -116.91973 4000 -22731.129 -22731.129 -22904.971 -22904.971 336.30948 336.30948 60668.005 60668.005 185.68305 185.68305 Loop time of 218.898 on 1 procs for 1000 steps with 4000 atoms Performance: 0.395 ns/day, 60.805 hours/ns, 4.568 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 | 217.82 | 217.82 | 217.82 | 0.0 | 99.51 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.1665 | 0.1665 | 0.1665 | 0.0 | 0.08 Output | 0.00023112 | 0.00023112 | 0.00023112 | 0.0 | 0.00 Modify | 0.80617 | 0.80617 | 0.80617 | 0.0 | 0.37 Other | | 0.1073 | | | 0.05 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: 311976 ave 311976 max 311976 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311976 Ave neighs/atom = 77.994 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 = 13 13 13 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.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 4000 -22731.129 -22731.129 -22904.971 -22904.971 336.30948 336.30948 60668.005 60668.005 185.68305 185.68305 5000 -22739.04 -22739.04 -22908.072 -22908.072 327.00331 327.00331 60662.931 60662.931 82.448927 82.448927 Loop time of 191.768 on 1 procs for 1000 steps with 4000 atoms Performance: 0.451 ns/day, 53.269 hours/ns, 5.215 timesteps/s 100.0% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 190.85 | 190.85 | 190.85 | 0.0 | 99.52 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.15037 | 0.15037 | 0.15037 | 0.0 | 0.08 Output | 0.00019193 | 0.00019193 | 0.00019193 | 0.0 | 0.00 Modify | 0.67579 | 0.67579 | 0.67579 | 0.0 | 0.35 Other | | 0.09618 | | | 0.05 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: 311986 ave 311986 max 311986 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311986 Ave neighs/atom = 77.9965 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 = 329.669121966326, Press = -144.201419186194 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 = 13 13 13 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.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 5000 -22739.04 -22739.04 -22908.072 -22908.072 327.00331 327.00331 60662.931 60662.931 82.448927 82.448927 6000 -22731.871 -22731.871 -22905.809 -22905.809 336.4946 336.4946 60687.848 60687.848 -529.70543 -529.70543 Loop time of 186.276 on 1 procs for 1000 steps with 4000 atoms Performance: 0.464 ns/day, 51.743 hours/ns, 5.368 timesteps/s 100.0% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 185.36 | 185.36 | 185.36 | 0.0 | 99.51 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.14751 | 0.14751 | 0.14751 | 0.0 | 0.08 Output | 0.00023485 | 0.00023485 | 0.00023485 | 0.0 | 0.00 Modify | 0.67625 | 0.67625 | 0.67625 | 0.0 | 0.36 Other | | 0.09459 | | | 0.05 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: 311970 ave 311970 max 311970 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311970 Ave neighs/atom = 77.9925 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 = 332.905288297028, Press = 10.3933265803032 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 = 13 13 13 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.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 6000 -22731.871 -22731.871 -22905.809 -22905.809 336.4946 336.4946 60687.848 60687.848 -529.70543 -529.70543 7000 -22738.584 -22738.584 -22906.355 -22906.355 324.56611 324.56611 60629.392 60629.392 1072.7001 1072.7001 Loop time of 209.414 on 1 procs for 1000 steps with 4000 atoms Performance: 0.413 ns/day, 58.171 hours/ns, 4.775 timesteps/s 100.0% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 208.36 | 208.36 | 208.36 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16103 | 0.16103 | 0.16103 | 0.0 | 0.08 Output | 0.00018942 | 0.00018942 | 0.00018942 | 0.0 | 0.00 Modify | 0.79173 | 0.79173 | 0.79173 | 0.0 | 0.38 Other | | 0.1029 | | | 0.05 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: 311974 ave 311974 max 311974 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311974 Ave neighs/atom = 77.9935 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 = 332.79033608566, Press = -18.8306339301346 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 = 13 13 13 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.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 7000 -22738.584 -22738.584 -22906.355 -22906.355 324.56611 324.56611 60629.392 60629.392 1072.7001 1072.7001 8000 -22734.696 -22734.696 -22909.252 -22909.252 337.69037 337.69037 60730.062 60730.062 -1871.9672 -1871.9672 Loop time of 218.887 on 1 procs for 1000 steps with 4000 atoms Performance: 0.395 ns/day, 60.802 hours/ns, 4.569 timesteps/s 99.8% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 217.78 | 217.78 | 217.78 | 0.0 | 99.49 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16555 | 0.16555 | 0.16555 | 0.0 | 0.08 Output | 0.00018922 | 0.00018922 | 0.00018922 | 0.0 | 0.00 Modify | 0.83891 | 0.83891 | 0.83891 | 0.0 | 0.38 Other | | 0.107 | | | 0.05 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: 311982 ave 311982 max 311982 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311982 Ave neighs/atom = 77.9955 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 = 332.940373049177, Press = 10.2124863256836 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 = 13 13 13 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.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 8000 -22734.696 -22734.696 -22909.252 -22909.252 337.69037 337.69037 60730.062 60730.062 -1871.9672 -1871.9672 9000 -22738.9 -22738.9 -22906.361 -22906.361 323.9654 323.9654 60607.118 60607.118 1646.7728 1646.7728 Loop time of 208.056 on 1 procs for 1000 steps with 4000 atoms Performance: 0.415 ns/day, 57.793 hours/ns, 4.806 timesteps/s 99.8% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 207 | 207 | 207 | 0.0 | 99.49 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16071 | 0.16071 | 0.16071 | 0.0 | 0.08 Output | 0.00019049 | 0.00019049 | 0.00019049 | 0.0 | 0.00 Modify | 0.78685 | 0.78685 | 0.78685 | 0.0 | 0.38 Other | | 0.1034 | | | 0.05 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: 311978 ave 311978 max 311978 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311978 Ave neighs/atom = 77.9945 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 = 332.925935636756, Press = -4.98167945739705 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 = 13 13 13 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.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 9000 -22738.9 -22738.9 -22906.361 -22906.361 323.9654 323.9654 60607.118 60607.118 1646.7728 1646.7728 10000 -22732.48 -22732.48 -22905.026 -22905.026 333.80104 333.80104 60692.364 60692.364 -590.38289 -590.38289 Loop time of 193.159 on 1 procs for 1000 steps with 4000 atoms Performance: 0.447 ns/day, 53.655 hours/ns, 5.177 timesteps/s 100.0% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 192.2 | 192.2 | 192.2 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.14893 | 0.14893 | 0.14893 | 0.0 | 0.08 Output | 0.00018859 | 0.00018859 | 0.00018859 | 0.0 | 0.00 Modify | 0.71349 | 0.71349 | 0.71349 | 0.0 | 0.37 Other | | 0.1002 | | | 0.05 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: 311966 ave 311966 max 311966 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311966 Ave neighs/atom = 77.9915 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 = 332.711073296549, Press = -2.27002786416188 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 = 13 13 13 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.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 10000 -22732.48 -22732.48 -22905.026 -22905.026 333.80104 333.80104 60692.364 60692.364 -590.38289 -590.38289 11000 -22738.937 -22738.937 -22912.745 -22912.745 336.2425 336.2425 60642.178 60642.178 339.81184 339.81184 Loop time of 184.204 on 1 procs for 1000 steps with 4000 atoms Performance: 0.469 ns/day, 51.168 hours/ns, 5.429 timesteps/s 100.0% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 183.29 | 183.29 | 183.29 | 0.0 | 99.51 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.14214 | 0.14214 | 0.14214 | 0.0 | 0.08 Output | 0.00019662 | 0.00019662 | 0.00019662 | 0.0 | 0.00 Modify | 0.67108 | 0.67108 | 0.67108 | 0.0 | 0.36 Other | | 0.09769 | | | 0.05 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: 311964 ave 311964 max 311964 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311964 Ave neighs/atom = 77.991 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 = 332.548813748526, Press = -0.21990115200365 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 = 13 13 13 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.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 11000 -22738.937 -22738.937 -22912.745 -22912.745 336.2425 336.2425 60642.178 60642.178 339.81184 339.81184 12000 -22733.001 -22733.001 -22906.675 -22906.675 335.98399 335.98399 60680.023 60680.023 -329.57434 -329.57434 Loop time of 200.002 on 1 procs for 1000 steps with 4000 atoms Performance: 0.432 ns/day, 55.556 hours/ns, 5.000 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 | 199 | 199 | 199 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.15409 | 0.15409 | 0.15409 | 0.0 | 0.08 Output | 0.00023288 | 0.00023288 | 0.00023288 | 0.0 | 0.00 Modify | 0.74597 | 0.74597 | 0.74597 | 0.0 | 0.37 Other | | 0.1007 | | | 0.05 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: 311988 ave 311988 max 311988 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311988 Ave neighs/atom = 77.997 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 = 332.688605278555, Press = -1.79657880645885 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 = 13 13 13 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.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 12000 -22733.001 -22733.001 -22906.675 -22906.675 335.98399 335.98399 60680.023 60680.023 -329.57434 -329.57434 13000 -22737.221 -22737.221 -22912.73 -22912.73 339.53506 339.53506 60634.102 60634.102 579.65928 579.65928 Loop time of 218.812 on 1 procs for 1000 steps with 4000 atoms Performance: 0.395 ns/day, 60.781 hours/ns, 4.570 timesteps/s 99.9% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 217.71 | 217.71 | 217.71 | 0.0 | 99.49 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16721 | 0.16721 | 0.16721 | 0.0 | 0.08 Output | 0.00024992 | 0.00024992 | 0.00024992 | 0.0 | 0.00 Modify | 0.83164 | 0.83164 | 0.83164 | 0.0 | 0.38 Other | | 0.106 | | | 0.05 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: 311978 ave 311978 max 311978 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311978 Ave neighs/atom = 77.9945 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 = 332.734629135333, Press = 0.43145036006743 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 = 13 13 13 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.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 13000 -22737.221 -22737.221 -22912.73 -22912.73 339.53506 339.53506 60634.102 60634.102 579.65928 579.65928 14000 -22726.801 -22726.801 -22902.253 -22902.253 339.42377 339.42377 60721.777 60721.777 -1181.9988 -1181.9988 Loop time of 218.822 on 1 procs for 1000 steps with 4000 atoms Performance: 0.395 ns/day, 60.784 hours/ns, 4.570 timesteps/s 100.0% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 217.71 | 217.71 | 217.71 | 0.0 | 99.49 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16733 | 0.16733 | 0.16733 | 0.0 | 0.08 Output | 0.00019153 | 0.00019153 | 0.00019153 | 0.0 | 0.00 Modify | 0.83922 | 0.83922 | 0.83922 | 0.0 | 0.38 Other | | 0.1069 | | | 0.05 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: 311978 ave 311978 max 311978 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311978 Ave neighs/atom = 77.9945 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 = 332.946562557522, Press = -3.53627056774667 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 = 13 13 13 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.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 14000 -22726.801 -22726.801 -22902.253 -22902.253 339.42377 339.42377 60721.777 60721.777 -1181.9988 -1181.9988 15000 -22733.457 -22733.457 -22905.486 -22905.486 332.80123 332.80123 60627.251 60627.251 1245.1385 1245.1385 Loop time of 219.09 on 1 procs for 1000 steps with 4000 atoms Performance: 0.394 ns/day, 60.858 hours/ns, 4.564 timesteps/s 100.0% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 217.97 | 217.97 | 217.97 | 0.0 | 99.49 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16748 | 0.16748 | 0.16748 | 0.0 | 0.08 Output | 0.00019694 | 0.00019694 | 0.00019694 | 0.0 | 0.00 Modify | 0.84527 | 0.84527 | 0.84527 | 0.0 | 0.39 Other | | 0.107 | | | 0.05 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: 311978 ave 311978 max 311978 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311978 Ave neighs/atom = 77.9945 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 = 333.015847513453, Press = 3.47108193736627 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 = 13 13 13 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.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 15000 -22733.457 -22733.457 -22905.486 -22905.486 332.80123 332.80123 60627.251 60627.251 1245.1385 1245.1385 16000 -22743.528 -22743.528 -22912.476 -22912.476 326.84171 326.84171 60687.529 60687.529 -1021.2862 -1021.2862 Loop time of 219.007 on 1 procs for 1000 steps with 4000 atoms Performance: 0.395 ns/day, 60.835 hours/ns, 4.566 timesteps/s 100.0% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 217.89 | 217.89 | 217.89 | 0.0 | 99.49 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16781 | 0.16781 | 0.16781 | 0.0 | 0.08 Output | 0.00025746 | 0.00025746 | 0.00025746 | 0.0 | 0.00 Modify | 0.8428 | 0.8428 | 0.8428 | 0.0 | 0.38 Other | | 0.106 | | | 0.05 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: 311974 ave 311974 max 311974 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311974 Ave neighs/atom = 77.9935 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 = 332.957236989178, Press = -5.35734132997796 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 = 13 13 13 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.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 16000 -22743.528 -22743.528 -22912.476 -22912.476 326.84171 326.84171 60687.529 60687.529 -1021.2862 -1021.2862 17000 -22734.173 -22734.173 -22906.044 -22906.044 332.49651 332.49651 60659.312 60659.312 352.15437 352.15437 Loop time of 194.54 on 1 procs for 1000 steps with 4000 atoms Performance: 0.444 ns/day, 54.039 hours/ns, 5.140 timesteps/s 100.0% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 193.57 | 193.57 | 193.57 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.1503 | 0.1503 | 0.1503 | 0.0 | 0.08 Output | 0.00019226 | 0.00019226 | 0.00019226 | 0.0 | 0.00 Modify | 0.72162 | 0.72162 | 0.72162 | 0.0 | 0.37 Other | | 0.1006 | | | 0.05 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: 311984 ave 311984 max 311984 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311984 Ave neighs/atom = 77.996 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 = 332.886977439011, Press = 3.01608770119449 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 = 13 13 13 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.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 17000 -22734.173 -22734.173 -22906.044 -22906.044 332.49651 332.49651 60659.312 60659.312 352.15437 352.15437 18000 -22740.753 -22740.753 -22909.658 -22909.658 326.75757 326.75757 60648.37 60648.37 304.45561 304.45561 Loop time of 184.484 on 1 procs for 1000 steps with 4000 atoms Performance: 0.468 ns/day, 51.246 hours/ns, 5.421 timesteps/s 100.0% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 183.57 | 183.57 | 183.57 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.14212 | 0.14212 | 0.14212 | 0.0 | 0.08 Output | 0.00019103 | 0.00019103 | 0.00019103 | 0.0 | 0.00 Modify | 0.67219 | 0.67219 | 0.67219 | 0.0 | 0.36 Other | | 0.09913 | | | 0.05 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: 311984 ave 311984 max 311984 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311984 Ave neighs/atom = 77.996 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 = 332.739567825419, Press = -2.48424831738521 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 = 13 13 13 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.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 18000 -22740.753 -22740.753 -22909.658 -22909.658 326.75757 326.75757 60648.37 60648.37 304.45561 304.45561 19000 -22734.541 -22734.541 -22906.137 -22906.137 331.9643 331.9643 60674.527 60674.527 -74.250698 -74.250698 Loop time of 184.024 on 1 procs for 1000 steps with 4000 atoms Performance: 0.470 ns/day, 51.118 hours/ns, 5.434 timesteps/s 100.0% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 183.11 | 183.11 | 183.11 | 0.0 | 99.51 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.14211 | 0.14211 | 0.14211 | 0.0 | 0.08 Output | 0.00019236 | 0.00019236 | 0.00019236 | 0.0 | 0.00 Modify | 0.66914 | 0.66914 | 0.66914 | 0.0 | 0.36 Other | | 0.09828 | | | 0.05 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: 311984 ave 311984 max 311984 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311984 Ave neighs/atom = 77.996 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 = 332.700983882972, Press = 1.00693197967183 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 = 13 13 13 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.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 19000 -22734.541 -22734.541 -22906.137 -22906.137 331.9643 331.9643 60674.527 60674.527 -74.250698 -74.250698 20000 -22730.842 -22730.842 -22903.519 -22903.519 334.05481 334.05481 60656.787 60656.787 606.36545 606.36545 Loop time of 196.293 on 1 procs for 1000 steps with 4000 atoms Performance: 0.440 ns/day, 54.526 hours/ns, 5.094 timesteps/s 100.0% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 195.32 | 195.32 | 195.32 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.1506 | 0.1506 | 0.1506 | 0.0 | 0.08 Output | 0.00019082 | 0.00019082 | 0.00019082 | 0.0 | 0.00 Modify | 0.72633 | 0.72633 | 0.72633 | 0.0 | 0.37 Other | | 0.1005 | | | 0.05 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: 311976 ave 311976 max 311976 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311976 Ave neighs/atom = 77.994 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 = 332.697019375008, Press = -1.32107295813457 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 = 13 13 13 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.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 20000 -22730.842 -22730.842 -22903.519 -22903.519 334.05481 334.05481 60656.787 60656.787 606.36545 606.36545 21000 -22736.965 -22736.965 -22907.638 -22907.638 330.17902 330.17902 60730.863 60730.863 -1858.1249 -1858.1249 Loop time of 215.883 on 1 procs for 1000 steps with 4000 atoms Performance: 0.400 ns/day, 59.967 hours/ns, 4.632 timesteps/s 99.8% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 214.79 | 214.79 | 214.79 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16453 | 0.16453 | 0.16453 | 0.0 | 0.08 Output | 0.00019258 | 0.00019258 | 0.00019258 | 0.0 | 0.00 Modify | 0.8176 | 0.8176 | 0.8176 | 0.0 | 0.38 Other | | 0.1064 | | | 0.05 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: 311968 ave 311968 max 311968 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311968 Ave neighs/atom = 77.992 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 = 332.823462667515, Press = 1.03262064128614 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 = 13 13 13 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.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 21000 -22736.965 -22736.965 -22907.638 -22907.638 330.17902 330.17902 60730.863 60730.863 -1858.1249 -1858.1249 22000 -22729.277 -22729.277 -22904.244 -22904.244 338.48558 338.48558 60563.993 60563.993 3171.0635 3171.0635 Loop time of 219.329 on 1 procs for 1000 steps with 4000 atoms Performance: 0.394 ns/day, 60.925 hours/ns, 4.559 timesteps/s 99.8% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 218.22 | 218.22 | 218.22 | 0.0 | 99.49 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16741 | 0.16741 | 0.16741 | 0.0 | 0.08 Output | 0.00019734 | 0.00019734 | 0.00019734 | 0.0 | 0.00 Modify | 0.83633 | 0.83633 | 0.83633 | 0.0 | 0.38 Other | | 0.1083 | | | 0.05 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: 311984 ave 311984 max 311984 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311984 Ave neighs/atom = 77.996 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 = 332.907889307709, Press = -1.04217233436677 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 = 13 13 13 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.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 22000 -22729.277 -22729.277 -22904.244 -22904.244 338.48558 338.48558 60563.993 60563.993 3171.0635 3171.0635 23000 -22737.622 -22737.622 -22905.652 -22905.652 325.06568 325.06568 60722.95 60722.95 -1610.6139 -1610.6139 Loop time of 218.49 on 1 procs for 1000 steps with 4000 atoms Performance: 0.395 ns/day, 60.692 hours/ns, 4.577 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 | 217.39 | 217.39 | 217.39 | 0.0 | 99.49 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16579 | 0.16579 | 0.16579 | 0.0 | 0.08 Output | 0.00018947 | 0.00018947 | 0.00018947 | 0.0 | 0.00 Modify | 0.83048 | 0.83048 | 0.83048 | 0.0 | 0.38 Other | | 0.1076 | | | 0.05 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: 311966 ave 311966 max 311966 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311966 Ave neighs/atom = 77.9915 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 = 332.943380400481, Press = -0.509657337617675 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 = 13 13 13 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.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 23000 -22737.622 -22737.622 -22905.652 -22905.652 325.06568 325.06568 60722.95 60722.95 -1610.6139 -1610.6139 24000 -22734.249 -22734.249 -22906.087 -22906.087 332.43254 332.43254 60630.366 60630.366 1152.0838 1152.0838 Loop time of 217.291 on 1 procs for 1000 steps with 4000 atoms Performance: 0.398 ns/day, 60.359 hours/ns, 4.602 timesteps/s 100.0% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 216.19 | 216.19 | 216.19 | 0.0 | 99.49 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16498 | 0.16498 | 0.16498 | 0.0 | 0.08 Output | 0.00019078 | 0.00019078 | 0.00019078 | 0.0 | 0.00 Modify | 0.82948 | 0.82948 | 0.82948 | 0.0 | 0.38 Other | | 0.1065 | | | 0.05 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: 311962 ave 311962 max 311962 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311962 Ave neighs/atom = 77.9905 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 = 332.939230386741, Press = 0.464417204735799 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 = 13 13 13 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.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 24000 -22734.249 -22734.249 -22906.087 -22906.087 332.43254 332.43254 60630.366 60630.366 1152.0838 1152.0838 25000 -22738.895 -22738.895 -22911.839 -22911.839 334.57209 334.57209 60683.066 60683.066 -754.12935 -754.12935 Loop time of 217.241 on 1 procs for 1000 steps with 4000 atoms Performance: 0.398 ns/day, 60.345 hours/ns, 4.603 timesteps/s 100.0% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 216.15 | 216.15 | 216.15 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.1631 | 0.1631 | 0.1631 | 0.0 | 0.08 Output | 0.0002503 | 0.0002503 | 0.0002503 | 0.0 | 0.00 Modify | 0.81745 | 0.81745 | 0.81745 | 0.0 | 0.38 Other | | 0.1056 | | | 0.05 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: 311954 ave 311954 max 311954 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311954 Ave neighs/atom = 77.9885 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 = 332.872096559234, Press = -1.17839117482021 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 = 13 13 13 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.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 25000 -22738.895 -22738.895 -22911.839 -22911.839 334.57209 334.57209 60683.066 60683.066 -754.12935 -754.12935 26000 -22734.436 -22734.436 -22904.799 -22904.799 329.57963 329.57963 60655.563 60655.563 493.11419 493.11419 Loop time of 217.477 on 1 procs for 1000 steps with 4000 atoms Performance: 0.397 ns/day, 60.410 hours/ns, 4.598 timesteps/s 100.0% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 216.36 | 216.36 | 216.36 | 0.0 | 99.49 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16745 | 0.16745 | 0.16745 | 0.0 | 0.08 Output | 0.00019246 | 0.00019246 | 0.00019246 | 0.0 | 0.00 Modify | 0.83758 | 0.83758 | 0.83758 | 0.0 | 0.39 Other | | 0.1066 | | | 0.05 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: 311978 ave 311978 max 311978 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311978 Ave neighs/atom = 77.9945 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 = 332.911783050652, Press = 0.911542503795297 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 = 13 13 13 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.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 26000 -22734.436 -22734.436 -22904.799 -22904.799 329.57963 329.57963 60655.563 60655.563 493.11419 493.11419 27000 -22737.485 -22737.485 -22906.844 -22906.844 327.63539 327.63539 60664.69 60664.69 59.111912 59.111912 Loop time of 217.353 on 1 procs for 1000 steps with 4000 atoms Performance: 0.398 ns/day, 60.376 hours/ns, 4.601 timesteps/s 100.0% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 216.27 | 216.27 | 216.27 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16388 | 0.16388 | 0.16388 | 0.0 | 0.08 Output | 0.00027257 | 0.00027257 | 0.00027257 | 0.0 | 0.00 Modify | 0.81567 | 0.81567 | 0.81567 | 0.0 | 0.38 Other | | 0.1043 | | | 0.05 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: 311956 ave 311956 max 311956 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311956 Ave neighs/atom = 77.989 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 = 332.975356087206, Press = -3.0686777627929 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 = 13 13 13 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.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 27000 -22737.485 -22737.485 -22906.844 -22906.844 327.63539 327.63539 60664.69 60664.69 59.111912 59.111912 28000 -22732.553 -22732.553 -22907.833 -22907.833 339.09082 339.09082 60708.426 60708.426 -1150.8622 -1150.8622 Loop time of 217.031 on 1 procs for 1000 steps with 4000 atoms Performance: 0.398 ns/day, 60.286 hours/ns, 4.608 timesteps/s 100.1% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 215.92 | 215.92 | 215.92 | 0.0 | 99.49 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16611 | 0.16611 | 0.16611 | 0.0 | 0.08 Output | 0.00027214 | 0.00027214 | 0.00027214 | 0.0 | 0.00 Modify | 0.83698 | 0.83698 | 0.83698 | 0.0 | 0.39 Other | | 0.1069 | | | 0.05 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: 311966 ave 311966 max 311966 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311966 Ave neighs/atom = 77.9915 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 = 333.068523490281, Press = 2.17558309468703 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 = 13 13 13 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.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 28000 -22732.553 -22732.553 -22907.833 -22907.833 339.09082 339.09082 60708.426 60708.426 -1150.8622 -1150.8622 29000 -22734.757 -22734.757 -22906.44 -22906.44 332.13329 332.13329 60625.627 60625.627 1269.8154 1269.8154 Loop time of 217.515 on 1 procs for 1000 steps with 4000 atoms Performance: 0.397 ns/day, 60.421 hours/ns, 4.597 timesteps/s 100.0% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 216.41 | 216.41 | 216.41 | 0.0 | 99.49 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16517 | 0.16517 | 0.16517 | 0.0 | 0.08 Output | 0.00027276 | 0.00027276 | 0.00027276 | 0.0 | 0.00 Modify | 0.83158 | 0.83158 | 0.83158 | 0.0 | 0.38 Other | | 0.1064 | | | 0.05 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: 311968 ave 311968 max 311968 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311968 Ave neighs/atom = 77.992 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 = 333.10356466134, Press = -1.20413729222365 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 = 13 13 13 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.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 29000 -22734.757 -22734.757 -22906.44 -22906.44 332.13329 332.13329 60625.627 60625.627 1269.8154 1269.8154 30000 -22734.611 -22734.611 -22907.945 -22907.945 335.3266 335.3266 60699.683 60699.683 -953.42887 -953.42887 Loop time of 217.693 on 1 procs for 1000 steps with 4000 atoms Performance: 0.397 ns/day, 60.470 hours/ns, 4.594 timesteps/s 100.0% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 216.58 | 216.58 | 216.58 | 0.0 | 99.49 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16668 | 0.16668 | 0.16668 | 0.0 | 0.08 Output | 0.00019137 | 0.00019137 | 0.00019137 | 0.0 | 0.00 Modify | 0.83981 | 0.83981 | 0.83981 | 0.0 | 0.39 Other | | 0.1072 | | | 0.05 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: 311982 ave 311982 max 311982 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311982 Ave neighs/atom = 77.9955 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 = 333.137174326002, Press = 0.730811032671893 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 = 13 13 13 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.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 30000 -22734.611 -22734.611 -22907.945 -22907.945 335.3266 335.3266 60699.683 60699.683 -953.42887 -953.42887 31000 -22734.133 -22734.133 -22908.823 -22908.823 337.95132 337.95132 60627.799 60627.799 1049.1465 1049.1465 Loop time of 217.766 on 1 procs for 1000 steps with 4000 atoms Performance: 0.397 ns/day, 60.490 hours/ns, 4.592 timesteps/s 100.0% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 216.66 | 216.66 | 216.66 | 0.0 | 99.49 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16413 | 0.16413 | 0.16413 | 0.0 | 0.08 Output | 0.00031549 | 0.00031549 | 0.00031549 | 0.0 | 0.00 Modify | 0.8359 | 0.8359 | 0.8359 | 0.0 | 0.38 Other | | 0.1076 | | | 0.05 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: 311978 ave 311978 max 311978 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311978 Ave neighs/atom = 77.9945 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 = 333.111022553702, Press = -0.897391688621584 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 = 13 13 13 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.613 | 4.613 | 4.613 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 31000 -22734.133 -22734.133 -22908.823 -22908.823 337.95132 337.95132 60627.799 60627.799 1049.1465 1049.1465 32000 -22737.691 -22737.691 -22905.878 -22905.878 325.37024 325.37024 60719.813 60719.813 -1458.505 -1458.505 Loop time of 209.188 on 1 procs for 1000 steps with 4000 atoms Performance: 0.413 ns/day, 58.108 hours/ns, 4.780 timesteps/s 100.0% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 208.12 | 208.12 | 208.12 | 0.0 | 99.49 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.16047 | 0.16047 | 0.16047 | 0.0 | 0.08 Output | 0.00019154 | 0.00019154 | 0.00019154 | 0.0 | 0.00 Modify | 0.80529 | 0.80529 | 0.80529 | 0.0 | 0.38 Other | | 0.1062 | | | 0.05 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: 311972 ave 311972 max 311972 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311972 Ave neighs/atom = 77.993 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_T333.15.out" else "print 'not_converged' file output/vol_T333.15.out" print '${V}' file output/vol_T333.15.out 60667.3818046979 print "LAMMPS calculation completed" LAMMPS calculation completed quit 0