# 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.302443556487561*${_u_distance} variable latticeconst_converted equal 3.302443556487561*1 lattice bcc ${latticeconst_converted} lattice bcc 3.30244355648756 Lattice spacing in x,y,z = 3.3024436 3.3024436 3.3024436 region simbox block 0 10 0 10 0 10 units lattice create_box 1 simbox Created orthogonal box = (0.0000000 0.0000000 0.0000000) to (33.024436 33.024436 33.024436) 1 by 1 by 1 MPI processor grid create_atoms 1 box Created 2000 atoms using lattice units in orthogonal box = (0.0000000 0.0000000 0.0000000) to (33.024436 33.024436 33.024436) create_atoms CPU = 0.004 seconds variable mass_converted equal 180.9479*${_u_mass} variable mass_converted equal 180.9479*1 kim_interactions Ta WARNING: 'kim_' has been renamed to 'kim '. Please update your input. kim interactions Ta #=== BEGIN kim interactions ================================== pair_style kim MEAM_LAMMPS_LeeBaskesKim_2001_Ta__MO_644143102837_000 pair_coeff * * Ta #=== END kim interactions ==================================== mass 1 ${mass_converted} mass 1 180.9479 # initial volume variable v equal vol # assign formula variable V0 equal ${v} # evaluate initial value variable V0 equal 36016.8901176252 variable V0_metal equal ${V0}/(${_u_distance}*${_u_distance}*${_u_distance}) variable V0_metal equal 36016.8901176252/(${_u_distance}*${_u_distance}*${_u_distance}) variable V0_metal equal 36016.8901176252/(1*${_u_distance}*${_u_distance}) variable V0_metal equal 36016.8901176252/(1*1*${_u_distance}) variable V0_metal equal 36016.8901176252/(1*1*1) variable V0_metal_times1000 equal ${V0_metal}*1000 variable V0_metal_times1000 equal 36016.8901176252*1000 print "Initial system volume: ${V0_metal} Angstroms^3" Initial system volume: 36016.8901176252 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 313.15*${_u_temperature} variable temp_converted equal 313.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 313.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 313.15 ${temp_converted} ${Tdamp_converted} iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 313.15 313.15 ${Tdamp_converted} iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 313.15 313.15 0.1 iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 313.15 313.15 0.1 iso 0 ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 313.15 313.15 0.1 iso 0 0 ${Pdamp_converted} fix ensemble all npt temp 313.15 313.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 "313.15 - 0.2" variable T_up equal "313.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_644143102837_000#item-citation CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 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) = 3.704 | 3.704 | 3.704 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 0 -16099.085 -16099.085 -16180 -16180 313.15 313.15 36016.89 36016.89 2399.6044 2399.6044 1000 -16014.194 -16014.194 -16085.744 -16085.744 276.90552 276.90552 36106.729 36106.729 6233.0658 6233.0658 Loop time of 99.0344 on 1 procs for 1000 steps with 2000 atoms Performance: 0.872 ns/day, 27.510 hours/ns, 10.098 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 | 98.535 | 98.535 | 98.535 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.077653 | 0.077653 | 0.077653 | 0.0 | 0.08 Output | 0.00022961 | 0.00022961 | 0.00022961 | 0.0 | 0.00 Modify | 0.36326 | 0.36326 | 0.36326 | 0.0 | 0.37 Other | | 0.058 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.000000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 0, Press = 0 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 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) = 3.704 | 3.704 | 3.704 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 1000 -16014.194 -16014.194 -16085.744 -16085.744 276.90552 276.90552 36106.729 36106.729 6233.0658 6233.0658 2000 -16013.771 -16013.771 -16097.994 -16097.994 325.95101 325.95101 36217.675 36217.675 -850.94592 -850.94592 Loop time of 99.7047 on 1 procs for 1000 steps with 2000 atoms Performance: 0.867 ns/day, 27.696 hours/ns, 10.030 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 | 99.217 | 99.217 | 99.217 | 0.0 | 99.51 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.07622 | 0.07622 | 0.07622 | 0.0 | 0.08 Output | 0.00019425 | 0.00019425 | 0.00019425 | 0.0 | 0.00 Modify | 0.35363 | 0.35363 | 0.35363 | 0.0 | 0.35 Other | | 0.05798 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 115978.0 ave 115978 max 115978 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 115978 Ave neighs/atom = 57.989000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 0, Press = 0 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 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) = 3.704 | 3.704 | 3.704 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 2000 -16013.771 -16013.771 -16097.994 -16097.994 325.95101 325.95101 36217.675 36217.675 -850.94592 -850.94592 3000 -16016.389 -16016.389 -16101.222 -16101.222 328.31206 328.31206 36178.725 36178.725 1066.2696 1066.2696 Loop time of 103.952 on 1 procs for 1000 steps with 2000 atoms Performance: 0.831 ns/day, 28.876 hours/ns, 9.620 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 | 103.44 | 103.44 | 103.44 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.078781 | 0.078781 | 0.078781 | 0.0 | 0.08 Output | 0.00019213 | 0.00019213 | 0.00019213 | 0.0 | 0.00 Modify | 0.37589 | 0.37589 | 0.37589 | 0.0 | 0.36 Other | | 0.0598 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 115966.0 ave 115966 max 115966 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 115966 Ave neighs/atom = 57.983000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 0, Press = 0 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 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) = 3.704 | 3.704 | 3.704 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 3000 -16016.389 -16016.389 -16101.222 -16101.222 328.31206 328.31206 36178.725 36178.725 1066.2696 1066.2696 4000 -16015.804 -16015.804 -16094.518 -16094.518 304.63215 304.63215 36200.882 36200.882 326.08253 326.08253 Loop time of 107.49 on 1 procs for 1000 steps with 2000 atoms Performance: 0.804 ns/day, 29.858 hours/ns, 9.303 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 | 106.95 | 106.95 | 106.95 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.080428 | 0.080428 | 0.080428 | 0.0 | 0.07 Output | 0.00019267 | 0.00019267 | 0.00019267 | 0.0 | 0.00 Modify | 0.39483 | 0.39483 | 0.39483 | 0.0 | 0.37 Other | | 0.06159 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 115974.0 ave 115974 max 115974 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 115974 Ave neighs/atom = 57.987000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 0, Press = 0 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 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) = 3.704 | 3.704 | 3.704 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 4000 -16015.804 -16015.804 -16094.518 -16094.518 304.63215 304.63215 36200.882 36200.882 326.08253 326.08253 5000 -16016.045 -16016.045 -16098.858 -16098.858 320.49431 320.49431 36210.733 36210.733 -641.25926 -641.25926 Loop time of 107.342 on 1 procs for 1000 steps with 2000 atoms Performance: 0.805 ns/day, 29.817 hours/ns, 9.316 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 | 106.81 | 106.81 | 106.81 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.080555 | 0.080555 | 0.080555 | 0.0 | 0.08 Output | 0.0001854 | 0.0001854 | 0.0001854 | 0.0 | 0.00 Modify | 0.39315 | 0.39315 | 0.39315 | 0.0 | 0.37 Other | | 0.06118 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 115966.0 ave 115966 max 115966 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 115966 Ave neighs/atom = 57.983000 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 = 303.089409542722, Press = -209.019905769434 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 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) = 3.704 | 3.704 | 3.704 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 5000 -16016.045 -16016.045 -16098.858 -16098.858 320.49431 320.49431 36210.733 36210.733 -641.25926 -641.25926 6000 -16016.899 -16016.899 -16095.941 -16095.941 305.90064 305.90064 36208.438 36208.438 -184.94237 -184.94237 Loop time of 105.735 on 1 procs for 1000 steps with 2000 atoms Performance: 0.817 ns/day, 29.371 hours/ns, 9.458 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 | 105.18 | 105.18 | 105.18 | 0.0 | 99.48 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.081481 | 0.081481 | 0.081481 | 0.0 | 0.08 Output | 0.00019536 | 0.00019536 | 0.00019536 | 0.0 | 0.00 Modify | 0.40745 | 0.40745 | 0.40745 | 0.0 | 0.39 Other | | 0.06117 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 115976.0 ave 115976 max 115976 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 115976 Ave neighs/atom = 57.988000 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 = 312.197107320621, Press = -40.1466712128268 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 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) = 3.704 | 3.704 | 3.704 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 6000 -16016.899 -16016.899 -16095.941 -16095.941 305.90064 305.90064 36208.438 36208.438 -184.94237 -184.94237 7000 -16017.102 -16017.102 -16093.126 -16093.126 294.2232 294.2232 36247.25 36247.25 -2091.4323 -2091.4323 Loop time of 105.479 on 1 procs for 1000 steps with 2000 atoms Performance: 0.819 ns/day, 29.300 hours/ns, 9.481 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 | 104.94 | 104.94 | 104.94 | 0.0 | 99.49 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.080188 | 0.080188 | 0.080188 | 0.0 | 0.08 Output | 0.00019844 | 0.00019844 | 0.00019844 | 0.0 | 0.00 Modify | 0.40232 | 0.40232 | 0.40232 | 0.0 | 0.38 Other | | 0.06051 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 115970.0 ave 115970 max 115970 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 115970 Ave neighs/atom = 57.985000 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 = 312.590521292136, Press = 21.328097631395 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 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) = 3.704 | 3.704 | 3.704 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 7000 -16017.102 -16017.102 -16093.126 -16093.126 294.2232 294.2232 36247.25 36247.25 -2091.4323 -2091.4323 8000 -16015.696 -16015.696 -16101.889 -16101.889 333.57758 333.57758 36225.402 36225.402 -1559.9312 -1559.9312 Loop time of 108.329 on 1 procs for 1000 steps with 2000 atoms Performance: 0.798 ns/day, 30.091 hours/ns, 9.231 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 | 107.77 | 107.77 | 107.77 | 0.0 | 99.49 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.08126 | 0.08126 | 0.08126 | 0.0 | 0.08 Output | 0.0001529 | 0.0001529 | 0.0001529 | 0.0 | 0.00 Modify | 0.41507 | 0.41507 | 0.41507 | 0.0 | 0.38 Other | | 0.06118 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 115954.0 ave 115954 max 115954 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 115954 Ave neighs/atom = 57.977000 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 = 312.545562827576, Press = 4.73084617594917 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 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) = 3.704 | 3.704 | 3.704 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 8000 -16015.696 -16015.696 -16101.889 -16101.889 333.57758 333.57758 36225.402 36225.402 -1559.9312 -1559.9312 9000 -16018.883 -16018.883 -16101.035 -16101.035 317.93739 317.93739 36223.136 36223.136 -1323.1834 -1323.1834 Loop time of 104.92 on 1 procs for 1000 steps with 2000 atoms Performance: 0.823 ns/day, 29.145 hours/ns, 9.531 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 | 104.38 | 104.38 | 104.38 | 0.0 | 99.49 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.079585 | 0.079585 | 0.079585 | 0.0 | 0.08 Output | 0.00015571 | 0.00015571 | 0.00015571 | 0.0 | 0.00 Modify | 0.39988 | 0.39988 | 0.39988 | 0.0 | 0.38 Other | | 0.05989 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 115972.0 ave 115972 max 115972 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 115972 Ave neighs/atom = 57.986000 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 = 312.422787975239, Press = 5.58204631892851 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 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) = 3.704 | 3.704 | 3.704 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 9000 -16018.883 -16018.883 -16101.035 -16101.035 317.93739 317.93739 36223.136 36223.136 -1323.1834 -1323.1834 10000 -16016.889 -16016.889 -16093.655 -16093.655 297.09304 297.09304 36241.966 36241.966 -2073.3044 -2073.3044 Loop time of 107.311 on 1 procs for 1000 steps with 2000 atoms Performance: 0.805 ns/day, 29.808 hours/ns, 9.319 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 | 106.75 | 106.75 | 106.75 | 0.0 | 99.48 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.0814 | 0.0814 | 0.0814 | 0.0 | 0.08 Output | 0.00015432 | 0.00015432 | 0.00015432 | 0.0 | 0.00 Modify | 0.41404 | 0.41404 | 0.41404 | 0.0 | 0.39 Other | | 0.06091 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 115980.0 ave 115980 max 115980 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 115980 Ave neighs/atom = 57.990000 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 = 312.33295977899, Press = 9.32967778856425 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 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) = 3.704 | 3.704 | 3.704 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 10000 -16016.889 -16016.889 -16093.655 -16093.655 297.09304 297.09304 36241.966 36241.966 -2073.3044 -2073.3044 11000 -16014.479 -16014.479 -16097.812 -16097.812 322.50889 322.50889 36173.314 36173.314 1792.4593 1792.4593 Loop time of 106.528 on 1 procs for 1000 steps with 2000 atoms Performance: 0.811 ns/day, 29.591 hours/ns, 9.387 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 | 105.98 | 105.98 | 105.98 | 0.0 | 99.49 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.080699 | 0.080699 | 0.080699 | 0.0 | 0.08 Output | 0.00020046 | 0.00020046 | 0.00020046 | 0.0 | 0.00 Modify | 0.40674 | 0.40674 | 0.40674 | 0.0 | 0.38 Other | | 0.0603 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 115952.0 ave 115952 max 115952 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 115952 Ave neighs/atom = 57.976000 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 = 312.499643094891, Press = -0.765697140494696 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 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) = 3.704 | 3.704 | 3.704 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 11000 -16014.479 -16014.479 -16097.812 -16097.812 322.50889 322.50889 36173.314 36173.314 1792.4593 1792.4593 12000 -16017.31 -16017.31 -16102.344 -16102.344 329.0904 329.0904 36210.86 36210.86 -672.32063 -672.32063 Loop time of 105.507 on 1 procs for 1000 steps with 2000 atoms Performance: 0.819 ns/day, 29.308 hours/ns, 9.478 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 | 104.96 | 104.96 | 104.96 | 0.0 | 99.49 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.080113 | 0.080113 | 0.080113 | 0.0 | 0.08 Output | 0.00015466 | 0.00015466 | 0.00015466 | 0.0 | 0.00 Modify | 0.40241 | 0.40241 | 0.40241 | 0.0 | 0.38 Other | | 0.06056 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 115966.0 ave 115966 max 115966 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 115966 Ave neighs/atom = 57.983000 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 = 312.518659284209, Press = 6.33112017405448 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 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) = 3.704 | 3.704 | 3.704 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 12000 -16017.31 -16017.31 -16102.344 -16102.344 329.0904 329.0904 36210.86 36210.86 -672.32063 -672.32063 13000 -16016.959 -16016.959 -16095.846 -16095.846 305.30183 305.30183 36220.672 36220.672 -833.25923 -833.25923 Loop time of 107.701 on 1 procs for 1000 steps with 2000 atoms Performance: 0.802 ns/day, 29.917 hours/ns, 9.285 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 | 107.14 | 107.14 | 107.14 | 0.0 | 99.48 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.080962 | 0.080962 | 0.080962 | 0.0 | 0.08 Output | 0.00018742 | 0.00018742 | 0.00018742 | 0.0 | 0.00 Modify | 0.41335 | 0.41335 | 0.41335 | 0.0 | 0.38 Other | | 0.06154 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 115972.0 ave 115972 max 115972 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 115972 Ave neighs/atom = 57.986000 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 = 312.803327704475, Press = -2.46451775332707 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 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) = 3.704 | 3.704 | 3.704 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 13000 -16016.959 -16016.959 -16095.846 -16095.846 305.30183 305.30183 36220.672 36220.672 -833.25923 -833.25923 14000 -16015.312 -16015.312 -16099.268 -16099.268 324.91886 324.91886 36224.626 36224.626 -1095.1244 -1095.1244 Loop time of 104.052 on 1 procs for 1000 steps with 2000 atoms Performance: 0.830 ns/day, 28.903 hours/ns, 9.611 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 | 103.52 | 103.52 | 103.52 | 0.0 | 99.49 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.079275 | 0.079275 | 0.079275 | 0.0 | 0.08 Output | 0.00015712 | 0.00015712 | 0.00015712 | 0.0 | 0.00 Modify | 0.39699 | 0.39699 | 0.39699 | 0.0 | 0.38 Other | | 0.05936 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 115966.0 ave 115966 max 115966 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 115966 Ave neighs/atom = 57.983000 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 = 313.022273349253, Press = 2.52129801180619 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 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) = 3.704 | 3.704 | 3.704 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 14000 -16015.312 -16015.312 -16099.268 -16099.268 324.91886 324.91886 36224.626 36224.626 -1095.1244 -1095.1244 15000 -16018.721 -16018.721 -16095.485 -16095.485 297.08409 297.08409 36207.74 36207.74 -164.68221 -164.68221 Loop time of 103.214 on 1 procs for 1000 steps with 2000 atoms Performance: 0.837 ns/day, 28.670 hours/ns, 9.689 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 | 102.68 | 102.68 | 102.68 | 0.0 | 99.48 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.079546 | 0.079546 | 0.079546 | 0.0 | 0.08 Output | 0.00019939 | 0.00019939 | 0.00019939 | 0.0 | 0.00 Modify | 0.39524 | 0.39524 | 0.39524 | 0.0 | 0.38 Other | | 0.06041 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 115958.0 ave 115958 max 115958 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 115958 Ave neighs/atom = 57.979000 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 = 313.023627684844, Press = 7.1942057245276 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 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) = 3.704 | 3.704 | 3.704 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 15000 -16018.721 -16018.721 -16095.485 -16095.485 297.08409 297.08409 36207.74 36207.74 -164.68221 -164.68221 16000 -16014.26 -16014.26 -16095.932 -16095.932 316.08254 316.08254 36146.059 36146.059 3148.8899 3148.8899 Loop time of 106.53 on 1 procs for 1000 steps with 2000 atoms Performance: 0.811 ns/day, 29.592 hours/ns, 9.387 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 | 105.98 | 105.98 | 105.98 | 0.0 | 99.49 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.081069 | 0.081069 | 0.081069 | 0.0 | 0.08 Output | 0.00015399 | 0.00015399 | 0.00015399 | 0.0 | 0.00 Modify | 0.40617 | 0.40617 | 0.40617 | 0.0 | 0.38 Other | | 0.06075 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 115972.0 ave 115972 max 115972 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 115972 Ave neighs/atom = 57.986000 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 = 313.04998629346, Press = 2.19600797231884 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 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) = 3.704 | 3.704 | 3.704 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 16000 -16014.26 -16014.26 -16095.932 -16095.932 316.08254 316.08254 36146.059 36146.059 3148.8899 3148.8899 17000 -16020.557 -16020.557 -16098.625 -16098.625 302.1317 302.1317 36207.566 36207.566 -508.23141 -508.23141 Loop time of 105.539 on 1 procs for 1000 steps with 2000 atoms Performance: 0.819 ns/day, 29.316 hours/ns, 9.475 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 | 104.99 | 104.99 | 104.99 | 0.0 | 99.48 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.080604 | 0.080604 | 0.080604 | 0.0 | 0.08 Output | 0.00015461 | 0.00015461 | 0.00015461 | 0.0 | 0.00 Modify | 0.40623 | 0.40623 | 0.40623 | 0.0 | 0.38 Other | | 0.06037 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 115970.0 ave 115970 max 115970 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 115970 Ave neighs/atom = 57.985000 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 = 312.898136351347, Press = -3.47168098317435 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 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) = 3.704 | 3.704 | 3.704 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 17000 -16020.557 -16020.557 -16098.625 -16098.625 302.1317 302.1317 36207.566 36207.566 -508.23141 -508.23141 18000 -16015.575 -16015.575 -16097.895 -16097.895 318.58927 318.58927 36257.766 36257.766 -2950.1579 -2950.1579 Loop time of 106.657 on 1 procs for 1000 steps with 2000 atoms Performance: 0.810 ns/day, 29.627 hours/ns, 9.376 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 | 106.11 | 106.11 | 106.11 | 0.0 | 99.48 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.081145 | 0.081145 | 0.081145 | 0.0 | 0.08 Output | 0.00020646 | 0.00020646 | 0.00020646 | 0.0 | 0.00 Modify | 0.4089 | 0.4089 | 0.4089 | 0.0 | 0.38 Other | | 0.06129 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 115982.0 ave 115982 max 115982 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 115982 Ave neighs/atom = 57.991000 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 = 312.755114426338, Press = 2.58659379205239 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 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) = 3.704 | 3.704 | 3.704 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 18000 -16015.575 -16015.575 -16097.895 -16097.895 318.58927 318.58927 36257.766 36257.766 -2950.1579 -2950.1579 19000 -16014.181 -16014.181 -16097.008 -16097.008 320.54863 320.54863 36181.976 36181.976 1346.1713 1346.1713 Loop time of 107.381 on 1 procs for 1000 steps with 2000 atoms Performance: 0.805 ns/day, 29.828 hours/ns, 9.313 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 | 106.82 | 106.82 | 106.82 | 0.0 | 99.48 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.082656 | 0.082656 | 0.082656 | 0.0 | 0.08 Output | 0.00015375 | 0.00015375 | 0.00015375 | 0.0 | 0.00 Modify | 0.41486 | 0.41486 | 0.41486 | 0.0 | 0.39 Other | | 0.06151 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 115980.0 ave 115980 max 115980 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 115980 Ave neighs/atom = 57.990000 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 = 312.877233058052, Press = 3.38754810828689 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 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) = 3.704 | 3.704 | 3.704 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 19000 -16014.181 -16014.181 -16097.008 -16097.008 320.54863 320.54863 36181.976 36181.976 1346.1713 1346.1713 20000 -16016.849 -16016.849 -16098.72 -16098.72 316.85092 316.85092 36187.897 36187.897 767.34083 767.34083 Loop time of 102.756 on 1 procs for 1000 steps with 2000 atoms Performance: 0.841 ns/day, 28.543 hours/ns, 9.732 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 | 102.23 | 102.23 | 102.23 | 0.0 | 99.49 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.078533 | 0.078533 | 0.078533 | 0.0 | 0.08 Output | 0.00020342 | 0.00020342 | 0.00020342 | 0.0 | 0.00 Modify | 0.38742 | 0.38742 | 0.38742 | 0.0 | 0.38 Other | | 0.06069 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 115986.0 ave 115986 max 115986 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 115986 Ave neighs/atom = 57.993000 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 = 312.967875798418, Press = 0.356721248005224 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 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) = 3.704 | 3.704 | 3.704 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 20000 -16016.849 -16016.849 -16098.72 -16098.72 316.85092 316.85092 36187.897 36187.897 767.34083 767.34083 21000 -16011.937 -16011.937 -16095.839 -16095.839 324.70913 324.70913 36193.5 36193.5 840.08073 840.08073 Loop time of 105.44 on 1 procs for 1000 steps with 2000 atoms Performance: 0.819 ns/day, 29.289 hours/ns, 9.484 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 | 104.9 | 104.9 | 104.9 | 0.0 | 99.48 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.080197 | 0.080197 | 0.080197 | 0.0 | 0.08 Output | 0.00025809 | 0.00025809 | 0.00025809 | 0.0 | 0.00 Modify | 0.40345 | 0.40345 | 0.40345 | 0.0 | 0.38 Other | | 0.05999 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 115980.0 ave 115980 max 115980 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 115980 Ave neighs/atom = 57.990000 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 = 313.067889718896, Press = -0.299843420630641 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 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) = 3.704 | 3.704 | 3.704 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 21000 -16011.937 -16011.937 -16095.839 -16095.839 324.70913 324.70913 36193.5 36193.5 840.08073 840.08073 22000 -16017.979 -16017.979 -16094.816 -16094.816 297.36709 297.36709 36188.661 36188.661 874.86435 874.86435 Loop time of 105.148 on 1 procs for 1000 steps with 2000 atoms Performance: 0.822 ns/day, 29.208 hours/ns, 9.510 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 | 104.61 | 104.61 | 104.61 | 0.0 | 99.48 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.08043 | 0.08043 | 0.08043 | 0.0 | 0.08 Output | 0.00015332 | 0.00015332 | 0.00015332 | 0.0 | 0.00 Modify | 0.40121 | 0.40121 | 0.40121 | 0.0 | 0.38 Other | | 0.06005 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 115970.0 ave 115970 max 115970 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 115970 Ave neighs/atom = 57.985000 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 = 313.061905708067, Press = 4.00796404754847 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 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) = 3.704 | 3.704 | 3.704 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 22000 -16017.979 -16017.979 -16094.816 -16094.816 297.36709 297.36709 36188.661 36188.661 874.86435 874.86435 23000 -16016.069 -16016.069 -16094.311 -16094.311 302.80315 302.80315 36143.855 36143.855 3354.7786 3354.7786 Loop time of 108.26 on 1 procs for 1000 steps with 2000 atoms Performance: 0.798 ns/day, 30.072 hours/ns, 9.237 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 | 107.7 | 107.7 | 107.7 | 0.0 | 99.48 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.081992 | 0.081992 | 0.081992 | 0.0 | 0.08 Output | 0.00019275 | 0.00019275 | 0.00019275 | 0.0 | 0.00 Modify | 0.41544 | 0.41544 | 0.41544 | 0.0 | 0.38 Other | | 0.0617 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 115970.0 ave 115970 max 115970 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 115970 Ave neighs/atom = 57.985000 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 = 313.108134228569, Press = 0.640379177387611 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 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) = 3.704 | 3.704 | 3.704 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 23000 -16016.069 -16016.069 -16094.311 -16094.311 302.80315 302.80315 36143.855 36143.855 3354.7786 3354.7786 24000 -16018.993 -16018.993 -16095.817 -16095.817 297.31589 297.31589 36237.172 36237.172 -1770.1851 -1770.1851 Loop time of 103.332 on 1 procs for 1000 steps with 2000 atoms Performance: 0.836 ns/day, 28.703 hours/ns, 9.678 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 | 102.8 | 102.8 | 102.8 | 0.0 | 99.49 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.078612 | 0.078612 | 0.078612 | 0.0 | 0.08 Output | 0.00015571 | 0.00015571 | 0.00015571 | 0.0 | 0.00 Modify | 0.39084 | 0.39084 | 0.39084 | 0.0 | 0.38 Other | | 0.05976 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 115968.0 ave 115968 max 115968 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 115968 Ave neighs/atom = 57.984000 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 = 312.890251017187, Press = -1.90272956769879 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 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) = 3.704 | 3.704 | 3.704 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 24000 -16018.993 -16018.993 -16095.817 -16095.817 297.31589 297.31589 36237.172 36237.172 -1770.1851 -1770.1851 25000 -16015.56 -16015.56 -16095.736 -16095.736 310.2891 310.2891 36227.578 36227.578 -1089.7073 -1089.7073 Loop time of 104.884 on 1 procs for 1000 steps with 2000 atoms Performance: 0.824 ns/day, 29.135 hours/ns, 9.534 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 | 104.35 | 104.35 | 104.35 | 0.0 | 99.49 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.079528 | 0.079528 | 0.079528 | 0.0 | 0.08 Output | 0.00020344 | 0.00020344 | 0.00020344 | 0.0 | 0.00 Modify | 0.39917 | 0.39917 | 0.39917 | 0.0 | 0.38 Other | | 0.06051 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 115962.0 ave 115962 max 115962 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 115962 Ave neighs/atom = 57.981000 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 = 312.88115140369, Press = 0.784144821159072 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 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) = 3.704 | 3.704 | 3.704 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 25000 -16015.56 -16015.56 -16095.736 -16095.736 310.2891 310.2891 36227.578 36227.578 -1089.7073 -1089.7073 26000 -16016.726 -16016.726 -16096.551 -16096.551 308.93349 308.93349 36218.024 36218.024 -652.70517 -652.70517 Loop time of 103.165 on 1 procs for 1000 steps with 2000 atoms Performance: 0.837 ns/day, 28.657 hours/ns, 9.693 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 | 102.64 | 102.64 | 102.64 | 0.0 | 99.49 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.078663 | 0.078663 | 0.078663 | 0.0 | 0.08 Output | 0.00019664 | 0.00019664 | 0.00019664 | 0.0 | 0.00 Modify | 0.38956 | 0.38956 | 0.38956 | 0.0 | 0.38 Other | | 0.05922 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 115974.0 ave 115974 max 115974 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 115974 Ave neighs/atom = 57.987000 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 = 312.900587012967, Press = 2.04637312753449 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 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) = 3.704 | 3.704 | 3.704 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 26000 -16016.726 -16016.726 -16096.551 -16096.551 308.93349 308.93349 36218.024 36218.024 -652.70517 -652.70517 27000 -16015.082 -16015.082 -16096.385 -16096.385 314.64887 314.64887 36173.151 36173.151 1762.6777 1762.6777 Loop time of 104.602 on 1 procs for 1000 steps with 2000 atoms Performance: 0.826 ns/day, 29.056 hours/ns, 9.560 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 | 104.07 | 104.07 | 104.07 | 0.0 | 99.49 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.078768 | 0.078768 | 0.078768 | 0.0 | 0.08 Output | 0.00019111 | 0.00019111 | 0.00019111 | 0.0 | 0.00 Modify | 0.39682 | 0.39682 | 0.39682 | 0.0 | 0.38 Other | | 0.06 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 115970.0 ave 115970 max 115970 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 115970 Ave neighs/atom = 57.985000 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 = 312.876875614583, Press = 2.38026094197945 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 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) = 3.704 | 3.704 | 3.704 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 27000 -16015.082 -16015.082 -16096.385 -16096.385 314.64887 314.64887 36173.151 36173.151 1762.6777 1762.6777 28000 -16018.278 -16018.278 -16099.932 -16099.932 316.01054 316.01054 36185.399 36185.399 817.20117 817.20117 Loop time of 103.171 on 1 procs for 1000 steps with 2000 atoms Performance: 0.837 ns/day, 28.659 hours/ns, 9.693 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 | 102.64 | 102.64 | 102.64 | 0.0 | 99.49 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.078522 | 0.078522 | 0.078522 | 0.0 | 0.08 Output | 0.00019926 | 0.00019926 | 0.00019926 | 0.0 | 0.00 Modify | 0.39155 | 0.39155 | 0.39155 | 0.0 | 0.38 Other | | 0.05986 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 115980.0 ave 115980 max 115980 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 115980 Ave neighs/atom = 57.990000 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 = 312.812106489488, Press = 0.866024383896713 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 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) = 3.704 | 3.704 | 3.704 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 28000 -16018.278 -16018.278 -16099.932 -16099.932 316.01054 316.01054 36185.399 36185.399 817.20117 817.20117 29000 -16011.354 -16011.354 -16094.456 -16094.456 321.61089 321.61089 36195.803 36195.803 708.12123 708.12123 Loop time of 104.203 on 1 procs for 1000 steps with 2000 atoms Performance: 0.829 ns/day, 28.945 hours/ns, 9.597 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 | 103.67 | 103.67 | 103.67 | 0.0 | 99.49 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.078494 | 0.078494 | 0.078494 | 0.0 | 0.08 Output | 0.00015283 | 0.00015283 | 0.00015283 | 0.0 | 0.00 Modify | 0.39248 | 0.39248 | 0.39248 | 0.0 | 0.38 Other | | 0.05949 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 115978.0 ave 115978 max 115978 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 115978 Ave neighs/atom = 57.989000 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 = 312.902228853418, Press = -0.704363275471234 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 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) = 3.704 | 3.704 | 3.704 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 29000 -16011.354 -16011.354 -16094.456 -16094.456 321.61089 321.61089 36195.803 36195.803 708.12123 708.12123 30000 -16013.944 -16013.944 -16094.46 -16094.46 311.60755 311.60755 36229.868 36229.868 -1130.0342 -1130.0342 Loop time of 105.271 on 1 procs for 1000 steps with 2000 atoms Performance: 0.821 ns/day, 29.242 hours/ns, 9.499 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 | 104.73 | 104.73 | 104.73 | 0.0 | 99.49 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.079219 | 0.079219 | 0.079219 | 0.0 | 0.08 Output | 0.00038319 | 0.00038319 | 0.00038319 | 0.0 | 0.00 Modify | 0.39895 | 0.39895 | 0.39895 | 0.0 | 0.38 Other | | 0.06026 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 115976.0 ave 115976 max 115976 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 115976 Ave neighs/atom = 57.988000 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 = 312.916274740248, Press = -1.13995194593712 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 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) = 3.704 | 3.704 | 3.704 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 30000 -16013.944 -16013.944 -16094.46 -16094.46 311.60755 311.60755 36229.868 36229.868 -1130.0342 -1130.0342 31000 -16018.259 -16018.259 -16098.388 -16098.388 310.10882 310.10882 36215.857 36215.857 -713.84381 -713.84381 Loop time of 103.786 on 1 procs for 1000 steps with 2000 atoms Performance: 0.832 ns/day, 28.829 hours/ns, 9.635 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 | 103.25 | 103.25 | 103.25 | 0.0 | 99.49 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.078897 | 0.078897 | 0.078897 | 0.0 | 0.08 Output | 0.00019746 | 0.00019746 | 0.00019746 | 0.0 | 0.00 Modify | 0.39223 | 0.39223 | 0.39223 | 0.0 | 0.38 Other | | 0.05959 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 115978.0 ave 115978 max 115978 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 115978 Ave neighs/atom = 57.989000 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 = 312.886920073363, Press = 1.78909358753308 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 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) = 3.704 | 3.704 | 3.704 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 31000 -16018.259 -16018.259 -16098.388 -16098.388 310.10882 310.10882 36215.857 36215.857 -713.84381 -713.84381 32000 -16017.411 -16017.411 -16095.774 -16095.774 303.26948 303.26948 36204.088 36204.088 -8.0624638 -8.0624638 Loop time of 103.056 on 1 procs for 1000 steps with 2000 atoms Performance: 0.838 ns/day, 28.627 hours/ns, 9.703 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 | 102.52 | 102.52 | 102.52 | 0.0 | 99.48 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.078583 | 0.078583 | 0.078583 | 0.0 | 0.08 Output | 0.00015431 | 0.00015431 | 0.00015431 | 0.0 | 0.00 Modify | 0.39288 | 0.39288 | 0.39288 | 0.0 | 0.38 Other | | 0.05938 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 115982.0 ave 115982 max 115982 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 115982 Ave neighs/atom = 57.991000 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 = 312.905034682452, Press = -0.16986698786467 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 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) = 3.704 | 3.704 | 3.704 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 32000 -16017.411 -16017.411 -16095.774 -16095.774 303.26948 303.26948 36204.088 36204.088 -8.0624638 -8.0624638 33000 -16015.822 -16015.822 -16099.701 -16099.701 324.61899 324.61899 36187.446 36187.446 848.65869 848.65869 Loop time of 106.151 on 1 procs for 1000 steps with 2000 atoms Performance: 0.814 ns/day, 29.486 hours/ns, 9.421 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 | 105.6 | 105.6 | 105.6 | 0.0 | 99.49 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.079904 | 0.079904 | 0.079904 | 0.0 | 0.08 Output | 0.00015903 | 0.00015903 | 0.00015903 | 0.0 | 0.00 Modify | 0.40557 | 0.40557 | 0.40557 | 0.0 | 0.38 Other | | 0.06101 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 115970.0 ave 115970 max 115970 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 115970 Ave neighs/atom = 57.985000 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 = 312.787903992732, Press = 0.244898367678001 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 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) = 3.704 | 3.704 | 3.704 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 33000 -16015.822 -16015.822 -16099.701 -16099.701 324.61899 324.61899 36187.446 36187.446 848.65869 848.65869 34000 -16014.817 -16014.817 -16096.149 -16096.149 314.7637 314.7637 36209.414 36209.414 -120.93468 -120.93468 Loop time of 104.946 on 1 procs for 1000 steps with 2000 atoms Performance: 0.823 ns/day, 29.152 hours/ns, 9.529 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 | 104.41 | 104.41 | 104.41 | 0.0 | 99.49 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.079234 | 0.079234 | 0.079234 | 0.0 | 0.08 Output | 0.00018595 | 0.00018595 | 0.00018595 | 0.0 | 0.00 Modify | 0.39547 | 0.39547 | 0.39547 | 0.0 | 0.38 Other | | 0.05986 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 115984.0 ave 115984 max 115984 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 115984 Ave neighs/atom = 57.992000 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 = 312.824194483963, Press = 0.992194923035761 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 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) = 3.704 | 3.704 | 3.704 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 34000 -16014.817 -16014.817 -16096.149 -16096.149 314.7637 314.7637 36209.414 36209.414 -120.93468 -120.93468 35000 -16017.637 -16017.637 -16096.129 -16096.129 303.76989 303.76989 36174.815 36174.815 1536.6966 1536.6966 Loop time of 106.527 on 1 procs for 1000 steps with 2000 atoms Performance: 0.811 ns/day, 29.591 hours/ns, 9.387 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 | 105.98 | 105.98 | 105.98 | 0.0 | 99.49 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.080653 | 0.080653 | 0.080653 | 0.0 | 0.08 Output | 0.00019972 | 0.00019972 | 0.00019972 | 0.0 | 0.00 Modify | 0.40482 | 0.40482 | 0.40482 | 0.0 | 0.38 Other | | 0.06102 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 115972.0 ave 115972 max 115972 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 115972 Ave neighs/atom = 57.986000 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 = 312.822106463268, Press = 1.6567444912213 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 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) = 3.704 | 3.704 | 3.704 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 35000 -16017.637 -16017.637 -16096.129 -16096.129 303.76989 303.76989 36174.815 36174.815 1536.6966 1536.6966 36000 -16015.208 -16015.208 -16094.098 -16094.098 305.31472 305.31472 36178.587 36178.587 1480.7359 1480.7359 Loop time of 102.771 on 1 procs for 1000 steps with 2000 atoms Performance: 0.841 ns/day, 28.548 hours/ns, 9.730 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 | 102.24 | 102.24 | 102.24 | 0.0 | 99.49 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.077898 | 0.077898 | 0.077898 | 0.0 | 0.08 Output | 0.00015242 | 0.00015242 | 0.00015242 | 0.0 | 0.00 Modify | 0.39053 | 0.39053 | 0.39053 | 0.0 | 0.38 Other | | 0.0591 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 115976.0 ave 115976 max 115976 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 115976 Ave neighs/atom = 57.988000 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 = 312.828912007421, Press = -1.23196050672661 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 6 ghost atom cutoff = 6 binsize = 3, bins = 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) = 3.704 | 3.704 | 3.704 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 36000 -16015.208 -16015.208 -16094.098 -16094.098 305.31472 305.31472 36178.587 36178.587 1480.7359 1480.7359 37000 -16014.294 -16014.294 -16095.047 -16095.047 312.5213 312.5213 36231.92 36231.92 -1324.1586 -1324.1586 Loop time of 102.275 on 1 procs for 1000 steps with 2000 atoms Performance: 0.845 ns/day, 28.410 hours/ns, 9.778 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 | 101.75 | 101.75 | 101.75 | 0.0 | 99.49 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.077463 | 0.077463 | 0.077463 | 0.0 | 0.08 Output | 0.00015653 | 0.00015653 | 0.00015653 | 0.0 | 0.00 Modify | 0.38398 | 0.38398 | 0.38398 | 0.0 | 0.38 Other | | 0.059 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 115978.0 ave 115978 max 115978 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 115978 Ave neighs/atom = 57.989000 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_T313.15.out" else "print 'not_converged' file output/vol_T313.15.out" print '${V}' file output/vol_T313.15.out 36204.0932898592 print "LAMMPS calculation completed" LAMMPS calculation completed quit 0