# 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.0000000 0.0000000 0.0000000) 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.0000000 0.0000000 0.0000000) to (39.173716 39.173716 39.173716) create_atoms CPU = 0.005 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_001 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_001#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.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 180.975 on 1 procs for 1000 steps with 4000 atoms Performance: 0.477 ns/day, 50.271 hours/ns, 5.526 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 | 180.05 | 180.05 | 180.05 | 0.0 | 99.49 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.15014 | 0.15014 | 0.15014 | 0.0 | 0.08 Output | 0.00027755 | 0.00027755 | 0.00027755 | 0.0 | 0.00 Modify | 0.67716 | 0.67716 | 0.67716 | 0.0 | 0.37 Other | | 0.09824 | | | 0.05 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 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: 312000.0 ave 312000 max 312000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 312000 Ave neighs/atom = 78.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.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 191.201 on 1 procs for 1000 steps with 4000 atoms Performance: 0.452 ns/day, 53.111 hours/ns, 5.230 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 | 190.28 | 190.28 | 190.28 | 0.0 | 99.52 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.14741 | 0.14741 | 0.14741 | 0.0 | 0.08 Output | 0.00025331 | 0.00025331 | 0.00025331 | 0.0 | 0.00 Modify | 0.67701 | 0.67701 | 0.67701 | 0.0 | 0.35 Other | | 0.09941 | | | 0.05 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 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: 311968.0 ave 311968 max 311968 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311968 Ave neighs/atom = 77.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 = 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.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 190.153 on 1 procs for 1000 steps with 4000 atoms Performance: 0.454 ns/day, 52.820 hours/ns, 5.259 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 | 189.23 | 189.23 | 189.23 | 0.0 | 99.52 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.14703 | 0.14703 | 0.14703 | 0.0 | 0.08 Output | 0.00026701 | 0.00026701 | 0.00026701 | 0.0 | 0.00 Modify | 0.67284 | 0.67284 | 0.67284 | 0.0 | 0.35 Other | | 0.09902 | | | 0.05 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 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: 311950.0 ave 311950 max 311950 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311950 Ave neighs/atom = 77.987500 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.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 190.281 on 1 procs for 1000 steps with 4000 atoms Performance: 0.454 ns/day, 52.856 hours/ns, 5.255 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 | 189.37 | 189.37 | 189.37 | 0.0 | 99.52 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.14644 | 0.14644 | 0.14644 | 0.0 | 0.08 Output | 0.00033732 | 0.00033732 | 0.00033732 | 0.0 | 0.00 Modify | 0.66439 | 0.66439 | 0.66439 | 0.0 | 0.35 Other | | 0.09859 | | | 0.05 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 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: 311976.0 ave 311976 max 311976 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311976 Ave neighs/atom = 77.994000 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.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 193.984 on 1 procs for 1000 steps with 4000 atoms Performance: 0.445 ns/day, 53.884 hours/ns, 5.155 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 | 193.05 | 193.05 | 193.05 | 0.0 | 99.52 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.14803 | 0.14803 | 0.14803 | 0.0 | 0.08 Output | 0.00019969 | 0.00019969 | 0.00019969 | 0.0 | 0.00 Modify | 0.68519 | 0.68519 | 0.68519 | 0.0 | 0.35 Other | | 0.0998 | | | 0.05 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 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: 311986.0 ave 311986 max 311986 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311986 Ave neighs/atom = 77.996500 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 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 192.709 on 1 procs for 1000 steps with 4000 atoms Performance: 0.448 ns/day, 53.530 hours/ns, 5.189 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 | 191.75 | 191.75 | 191.75 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.1473 | 0.1473 | 0.1473 | 0.0 | 0.08 Output | 0.0001749 | 0.0001749 | 0.0001749 | 0.0 | 0.00 Modify | 0.70942 | 0.70942 | 0.70942 | 0.0 | 0.37 Other | | 0.09942 | | | 0.05 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 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: 311970.0 ave 311970 max 311970 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311970 Ave neighs/atom = 77.992500 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 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 192.01 on 1 procs for 1000 steps with 4000 atoms Performance: 0.450 ns/day, 53.336 hours/ns, 5.208 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 | 191.07 | 191.07 | 191.07 | 0.0 | 99.51 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.14749 | 0.14749 | 0.14749 | 0.0 | 0.08 Output | 0.00021872 | 0.00021872 | 0.00021872 | 0.0 | 0.00 Modify | 0.69577 | 0.69577 | 0.69577 | 0.0 | 0.36 Other | | 0.09916 | | | 0.05 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 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: 311974.0 ave 311974 max 311974 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311974 Ave neighs/atom = 77.993500 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 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 190.166 on 1 procs for 1000 steps with 4000 atoms Performance: 0.454 ns/day, 52.824 hours/ns, 5.259 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 | 189.22 | 189.22 | 189.22 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.14712 | 0.14712 | 0.14712 | 0.0 | 0.08 Output | 0.00017765 | 0.00017765 | 0.00017765 | 0.0 | 0.00 Modify | 0.69548 | 0.69548 | 0.69548 | 0.0 | 0.37 Other | | 0.09954 | | | 0.05 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 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: 311982.0 ave 311982 max 311982 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311982 Ave neighs/atom = 77.995500 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 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 190.842 on 1 procs for 1000 steps with 4000 atoms Performance: 0.453 ns/day, 53.012 hours/ns, 5.240 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 | 189.91 | 189.91 | 189.91 | 0.0 | 99.51 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.14652 | 0.14652 | 0.14652 | 0.0 | 0.08 Output | 0.00022089 | 0.00022089 | 0.00022089 | 0.0 | 0.00 Modify | 0.68851 | 0.68851 | 0.68851 | 0.0 | 0.36 Other | | 0.0981 | | | 0.05 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 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: 311978.0 ave 311978 max 311978 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311978 Ave neighs/atom = 77.994500 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 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 191.223 on 1 procs for 1000 steps with 4000 atoms Performance: 0.452 ns/day, 53.117 hours/ns, 5.230 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 | 190.28 | 190.28 | 190.28 | 0.0 | 99.51 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.14716 | 0.14716 | 0.14716 | 0.0 | 0.08 Output | 0.00017959 | 0.00017959 | 0.00017959 | 0.0 | 0.00 Modify | 0.69665 | 0.69665 | 0.69665 | 0.0 | 0.36 Other | | 0.09964 | | | 0.05 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 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: 311966.0 ave 311966 max 311966 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311966 Ave neighs/atom = 77.991500 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 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 191.547 on 1 procs for 1000 steps with 4000 atoms Performance: 0.451 ns/day, 53.208 hours/ns, 5.221 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 | 190.6 | 190.6 | 190.6 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.14833 | 0.14833 | 0.14833 | 0.0 | 0.08 Output | 0.0001775 | 0.0001775 | 0.0001775 | 0.0 | 0.00 Modify | 0.70103 | 0.70103 | 0.70103 | 0.0 | 0.37 Other | | 0.09925 | | | 0.05 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 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: 311964.0 ave 311964 max 311964 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311964 Ave neighs/atom = 77.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 = 332.548813748526, Press = -0.21990115200365 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.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 191.016 on 1 procs for 1000 steps with 4000 atoms Performance: 0.452 ns/day, 53.060 hours/ns, 5.235 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 | 190.07 | 190.07 | 190.07 | 0.0 | 99.51 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.14798 | 0.14798 | 0.14798 | 0.0 | 0.08 Output | 0.00017701 | 0.00017701 | 0.00017701 | 0.0 | 0.00 Modify | 0.69657 | 0.69657 | 0.69657 | 0.0 | 0.36 Other | | 0.09933 | | | 0.05 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 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: 311988.0 ave 311988 max 311988 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311988 Ave neighs/atom = 77.997000 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 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 190.79 on 1 procs for 1000 steps with 4000 atoms Performance: 0.453 ns/day, 52.997 hours/ns, 5.241 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 | 189.85 | 189.85 | 189.85 | 0.0 | 99.51 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.14744 | 0.14744 | 0.14744 | 0.0 | 0.08 Output | 0.00017822 | 0.00017822 | 0.00017822 | 0.0 | 0.00 Modify | 0.69228 | 0.69228 | 0.69228 | 0.0 | 0.36 Other | | 0.09857 | | | 0.05 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 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: 311978.0 ave 311978 max 311978 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311978 Ave neighs/atom = 77.994500 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 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 192.407 on 1 procs for 1000 steps with 4000 atoms Performance: 0.449 ns/day, 53.446 hours/ns, 5.197 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 | 191.45 | 191.45 | 191.45 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.14791 | 0.14791 | 0.14791 | 0.0 | 0.08 Output | 0.00017575 | 0.00017575 | 0.00017575 | 0.0 | 0.00 Modify | 0.70806 | 0.70806 | 0.70806 | 0.0 | 0.37 Other | | 0.09942 | | | 0.05 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 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: 311978.0 ave 311978 max 311978 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311978 Ave neighs/atom = 77.994500 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 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 192.489 on 1 procs for 1000 steps with 4000 atoms Performance: 0.449 ns/day, 53.469 hours/ns, 5.195 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 | 191.54 | 191.54 | 191.54 | 0.0 | 99.51 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.14726 | 0.14726 | 0.14726 | 0.0 | 0.08 Output | 0.00049448 | 0.00049448 | 0.00049448 | 0.0 | 0.00 Modify | 0.70371 | 0.70371 | 0.70371 | 0.0 | 0.37 Other | | 0.09942 | | | 0.05 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 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: 311978.0 ave 311978 max 311978 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311978 Ave neighs/atom = 77.994500 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 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 191.421 on 1 procs for 1000 steps with 4000 atoms Performance: 0.451 ns/day, 53.172 hours/ns, 5.224 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.48 | 190.48 | 190.48 | 0.0 | 99.51 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.14621 | 0.14621 | 0.14621 | 0.0 | 0.08 Output | 0.00017943 | 0.00017943 | 0.00017943 | 0.0 | 0.00 Modify | 0.69395 | 0.69395 | 0.69395 | 0.0 | 0.36 Other | | 0.09974 | | | 0.05 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 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: 311974.0 ave 311974 max 311974 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311974 Ave neighs/atom = 77.993500 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 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 192.165 on 1 procs for 1000 steps with 4000 atoms Performance: 0.450 ns/day, 53.379 hours/ns, 5.204 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 | 191.21 | 191.21 | 191.21 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.14695 | 0.14695 | 0.14695 | 0.0 | 0.08 Output | 0.00018116 | 0.00018116 | 0.00018116 | 0.0 | 0.00 Modify | 0.70498 | 0.70498 | 0.70498 | 0.0 | 0.37 Other | | 0.09958 | | | 0.05 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 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: 311984.0 ave 311984 max 311984 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311984 Ave neighs/atom = 77.996000 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 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 190.755 on 1 procs for 1000 steps with 4000 atoms Performance: 0.453 ns/day, 52.988 hours/ns, 5.242 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 | 189.81 | 189.81 | 189.81 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.14721 | 0.14721 | 0.14721 | 0.0 | 0.08 Output | 0.00023826 | 0.00023826 | 0.00023826 | 0.0 | 0.00 Modify | 0.69808 | 0.69808 | 0.69808 | 0.0 | 0.37 Other | | 0.09947 | | | 0.05 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 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: 311984.0 ave 311984 max 311984 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311984 Ave neighs/atom = 77.996000 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 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 192.43 on 1 procs for 1000 steps with 4000 atoms Performance: 0.449 ns/day, 53.453 hours/ns, 5.197 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 | 191.47 | 191.47 | 191.47 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.14841 | 0.14841 | 0.14841 | 0.0 | 0.08 Output | 0.00025258 | 0.00025258 | 0.00025258 | 0.0 | 0.00 Modify | 0.711 | 0.711 | 0.711 | 0.0 | 0.37 Other | | 0.09989 | | | 0.05 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 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: 311984.0 ave 311984 max 311984 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311984 Ave neighs/atom = 77.996000 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 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 193.544 on 1 procs for 1000 steps with 4000 atoms Performance: 0.446 ns/day, 53.762 hours/ns, 5.167 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 | 192.59 | 192.59 | 192.59 | 0.0 | 99.51 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.14748 | 0.14748 | 0.14748 | 0.0 | 0.08 Output | 0.00021781 | 0.00021781 | 0.00021781 | 0.0 | 0.00 Modify | 0.70598 | 0.70598 | 0.70598 | 0.0 | 0.36 Other | | 0.1009 | | | 0.05 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 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: 311976.0 ave 311976 max 311976 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311976 Ave neighs/atom = 77.994000 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 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 194.824 on 1 procs for 1000 steps with 4000 atoms Performance: 0.443 ns/day, 54.118 hours/ns, 5.133 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 | 193.86 | 193.86 | 193.86 | 0.0 | 99.51 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.14941 | 0.14941 | 0.14941 | 0.0 | 0.08 Output | 0.00017743 | 0.00017743 | 0.00017743 | 0.0 | 0.00 Modify | 0.71199 | 0.71199 | 0.71199 | 0.0 | 0.37 Other | | 0.1006 | | | 0.05 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 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: 311968.0 ave 311968 max 311968 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311968 Ave neighs/atom = 77.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 = 332.823462667515, Press = 1.03262064128614 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.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 194.128 on 1 procs for 1000 steps with 4000 atoms Performance: 0.445 ns/day, 53.924 hours/ns, 5.151 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 | 193.16 | 193.16 | 193.16 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.14913 | 0.14913 | 0.14913 | 0.0 | 0.08 Output | 0.0001788 | 0.0001788 | 0.0001788 | 0.0 | 0.00 Modify | 0.7171 | 0.7171 | 0.7171 | 0.0 | 0.37 Other | | 0.09992 | | | 0.05 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 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: 311984.0 ave 311984 max 311984 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311984 Ave neighs/atom = 77.996000 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 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 191.778 on 1 procs for 1000 steps with 4000 atoms Performance: 0.451 ns/day, 53.272 hours/ns, 5.214 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 | 190.82 | 190.82 | 190.82 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.14827 | 0.14827 | 0.14827 | 0.0 | 0.08 Output | 0.00017752 | 0.00017752 | 0.00017752 | 0.0 | 0.00 Modify | 0.70851 | 0.70851 | 0.70851 | 0.0 | 0.37 Other | | 0.1 | | | 0.05 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 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: 311966.0 ave 311966 max 311966 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311966 Ave neighs/atom = 77.991500 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 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 194.899 on 1 procs for 1000 steps with 4000 atoms Performance: 0.443 ns/day, 54.139 hours/ns, 5.131 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.94 | 193.94 | 193.94 | 0.0 | 99.51 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.14927 | 0.14927 | 0.14927 | 0.0 | 0.08 Output | 0.00023608 | 0.00023608 | 0.00023608 | 0.0 | 0.00 Modify | 0.71103 | 0.71103 | 0.71103 | 0.0 | 0.36 Other | | 0.09911 | | | 0.05 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 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: 311962.0 ave 311962 max 311962 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311962 Ave neighs/atom = 77.990500 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 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 195.459 on 1 procs for 1000 steps with 4000 atoms Performance: 0.442 ns/day, 54.294 hours/ns, 5.116 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 | 194.48 | 194.48 | 194.48 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.15108 | 0.15108 | 0.15108 | 0.0 | 0.08 Output | 0.00017109 | 0.00017109 | 0.00017109 | 0.0 | 0.00 Modify | 0.7281 | 0.7281 | 0.7281 | 0.0 | 0.37 Other | | 0.1014 | | | 0.05 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 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: 311954.0 ave 311954 max 311954 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311954 Ave neighs/atom = 77.988500 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 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 196.29 on 1 procs for 1000 steps with 4000 atoms Performance: 0.440 ns/day, 54.525 hours/ns, 5.095 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 | 195.32 | 195.32 | 195.32 | 0.0 | 99.51 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.14918 | 0.14918 | 0.14918 | 0.0 | 0.08 Output | 0.00017955 | 0.00017955 | 0.00017955 | 0.0 | 0.00 Modify | 0.71829 | 0.71829 | 0.71829 | 0.0 | 0.37 Other | | 0.1004 | | | 0.05 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 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: 311978.0 ave 311978 max 311978 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311978 Ave neighs/atom = 77.994500 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 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 195.532 on 1 procs for 1000 steps with 4000 atoms Performance: 0.442 ns/day, 54.314 hours/ns, 5.114 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 | 194.55 | 194.55 | 194.55 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.15006 | 0.15006 | 0.15006 | 0.0 | 0.08 Output | 0.00017686 | 0.00017686 | 0.00017686 | 0.0 | 0.00 Modify | 0.72818 | 0.72818 | 0.72818 | 0.0 | 0.37 Other | | 0.1011 | | | 0.05 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 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: 311956.0 ave 311956 max 311956 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311956 Ave neighs/atom = 77.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 = 332.975356087206, Press = -3.0686777627929 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.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 190.123 on 1 procs for 1000 steps with 4000 atoms Performance: 0.454 ns/day, 52.812 hours/ns, 5.260 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 | 189.19 | 189.19 | 189.19 | 0.0 | 99.51 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.14589 | 0.14589 | 0.14589 | 0.0 | 0.08 Output | 0.00017338 | 0.00017338 | 0.00017338 | 0.0 | 0.00 Modify | 0.69172 | 0.69172 | 0.69172 | 0.0 | 0.36 Other | | 0.09965 | | | 0.05 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 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: 311966.0 ave 311966 max 311966 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311966 Ave neighs/atom = 77.991500 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 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 192.134 on 1 procs for 1000 steps with 4000 atoms Performance: 0.450 ns/day, 53.371 hours/ns, 5.205 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 | 191.18 | 191.18 | 191.18 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.14877 | 0.14877 | 0.14877 | 0.0 | 0.08 Output | 0.00028838 | 0.00028838 | 0.00028838 | 0.0 | 0.00 Modify | 0.70594 | 0.70594 | 0.70594 | 0.0 | 0.37 Other | | 0.09983 | | | 0.05 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 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: 311968.0 ave 311968 max 311968 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311968 Ave neighs/atom = 77.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 = 333.10356466134, Press = -1.20413729222365 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.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 191.886 on 1 procs for 1000 steps with 4000 atoms Performance: 0.450 ns/day, 53.302 hours/ns, 5.211 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.94 | 190.94 | 190.94 | 0.0 | 99.51 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.1487 | 0.1487 | 0.1487 | 0.0 | 0.08 Output | 0.0001721 | 0.0001721 | 0.0001721 | 0.0 | 0.00 Modify | 0.69953 | 0.69953 | 0.69953 | 0.0 | 0.36 Other | | 0.09983 | | | 0.05 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 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: 311982.0 ave 311982 max 311982 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311982 Ave neighs/atom = 77.995500 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 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 193.264 on 1 procs for 1000 steps with 4000 atoms Performance: 0.447 ns/day, 53.684 hours/ns, 5.174 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.3 | 192.3 | 192.3 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.14959 | 0.14959 | 0.14959 | 0.0 | 0.08 Output | 0.0001788 | 0.0001788 | 0.0001788 | 0.0 | 0.00 Modify | 0.71524 | 0.71524 | 0.71524 | 0.0 | 0.37 Other | | 0.1001 | | | 0.05 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 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: 311978.0 ave 311978 max 311978 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311978 Ave neighs/atom = 77.994500 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 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 193.312 on 1 procs for 1000 steps with 4000 atoms Performance: 0.447 ns/day, 53.698 hours/ns, 5.173 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 | 192.36 | 192.36 | 192.36 | 0.0 | 99.51 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.14843 | 0.14843 | 0.14843 | 0.0 | 0.08 Output | 0.00024484 | 0.00024484 | 0.00024484 | 0.0 | 0.00 Modify | 0.70508 | 0.70508 | 0.70508 | 0.0 | 0.36 Other | | 0.09928 | | | 0.05 Nlocal: 4000.00 ave 4000 max 4000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5841.00 ave 5841 max 5841 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: 311972.0 ave 311972 max 311972 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 311972 Ave neighs/atom = 77.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" 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