# Variables that can be adjusted by kim-lammps-preprocessor to switch unit sets for # Simulator Models variable _u_distance equal 1.0 variable _u_energy equal 1.0 variable _u_mass equal 1.0 variable _u_time equal 1.0 variable _u_pressure equal 1.0 variable _u_temperature equal 1.0 # This line may be swapped out by kim-lammps-preprocessor if running against a Simulator # Model whose atom_style is not 'atomic' atom_style atomic # 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 2.855324655771256*${_u_distance} variable latticeconst_converted equal 2.855324655771256*1 lattice bcc ${latticeconst_converted} lattice bcc 2.85532465577126 Lattice spacing in x,y,z = 2.85532 2.85532 2.85532 region simbox block 0 10 0 10 0 10 units lattice create_box 1 simbox Created orthogonal box = (0 0 0) to (28.5532 28.5532 28.5532) 1 by 1 by 1 MPI processor grid create_atoms 1 box Created 2000 atoms create_atoms CPU = 0.000348091 secs variable mass_converted equal 55.845*${_u_mass} variable mass_converted equal 55.845*1 # specify which KIM Model to use pair_style kim EAM_Dynamo_Mendelev_2003_Fe__MO_546673549085_000 pair_coeff * * Fe mass 1 ${mass_converted} mass 1 55.845 # initial volume variable v equal vol # assign formula variable V0 equal ${v} # evaluate initial value variable V0 equal 23279.1161097207 variable V0_metal equal ${V0}/(${_u_distance}*${_u_distance}*${_u_distance}) variable V0_metal equal 23279.1161097207/(${_u_distance}*${_u_distance}*${_u_distance}) variable V0_metal equal 23279.1161097207/(1*${_u_distance}*${_u_distance}) variable V0_metal equal 23279.1161097207/(1*1*${_u_distance}) variable V0_metal equal 23279.1161097207/(1*1*1) variable V0_metal_times1000 equal ${V0_metal}*1000 variable V0_metal_times1000 equal 23279.1161097207*1000 print "Initial system volume: ${V0_metal} Angstroms^3" Initial system volume: 23279.1161097207 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 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 7.3 ghost atom cutoff = 7.3 binsize = 3.65, bins = 8 8 8 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7.3 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.195 | 4.195 | 4.195 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 0 -8158.7872 -8158.7872 -8244.8702 -8244.8702 333.15 333.15 23279.116 23279.116 3949.6836 3949.6836 1000 -8066.9646 -8066.9646 -8147.5701 -8147.5701 311.95137 311.95137 23388.021 23388.021 -138.72209 -138.72209 Loop time of 23.437 on 1 procs for 1000 steps with 2000 atoms Performance: 3.686 ns/day, 6.510 hours/ns, 42.668 timesteps/s 26.5% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 22.816 | 22.816 | 22.816 | 0.0 | 97.35 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.22796 | 0.22796 | 0.22796 | 0.0 | 0.97 Output | 4.5061e-05 | 4.5061e-05 | 4.5061e-05 | 0.0 | 0.00 Modify | 0.37896 | 0.37896 | 0.37896 | 0.0 | 1.62 Other | | 0.01383 | | | 0.06 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5471 ave 5471 max 5471 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 272000 ave 272000 max 272000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 272000 Ave neighs/atom = 136 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 Per MPI rank memory allocation (min/avg/max) = 4.195 | 4.195 | 4.195 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 1000 -8066.9646 -8066.9646 -8147.5701 -8147.5701 311.95137 311.95137 23388.021 23388.021 -138.72209 -138.72209 2000 -8067.406 -8067.406 -8151.8256 -8151.8256 326.71267 326.71267 23373.826 23373.826 1005.7786 1005.7786 Loop time of 23.8234 on 1 procs for 1000 steps with 2000 atoms Performance: 3.627 ns/day, 6.618 hours/ns, 41.976 timesteps/s 27.5% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 23.373 | 23.373 | 23.373 | 0.0 | 98.11 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.067348 | 0.067348 | 0.067348 | 0.0 | 0.28 Output | 4.4107e-05 | 4.4107e-05 | 4.4107e-05 | 0.0 | 0.00 Modify | 0.339 | 0.339 | 0.339 | 0.0 | 1.42 Other | | 0.04413 | | | 0.19 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5419 ave 5419 max 5419 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 278974 ave 278974 max 278974 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 278974 Ave neighs/atom = 139.487 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 Per MPI rank memory allocation (min/avg/max) = 4.195 | 4.195 | 4.195 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 2000 -8067.406 -8067.406 -8151.8256 -8151.8256 326.71267 326.71267 23373.826 23373.826 1005.7786 1005.7786 3000 -8067.6418 -8067.6418 -8151.4449 -8151.4449 324.3265 324.3265 23397.429 23397.429 -898.51408 -898.51408 Loop time of 23.6733 on 1 procs for 1000 steps with 2000 atoms Performance: 3.650 ns/day, 6.576 hours/ns, 42.242 timesteps/s 27.3% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 23.133 | 23.133 | 23.133 | 0.0 | 97.72 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.17714 | 0.17714 | 0.17714 | 0.0 | 0.75 Output | 3.5048e-05 | 3.5048e-05 | 3.5048e-05 | 0.0 | 0.00 Modify | 0.34789 | 0.34789 | 0.34789 | 0.0 | 1.47 Other | | 0.01478 | | | 0.06 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5451 ave 5451 max 5451 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 278138 ave 278138 max 278138 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 278138 Ave neighs/atom = 139.069 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 Per MPI rank memory allocation (min/avg/max) = 4.195 | 4.195 | 4.195 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 3000 -8067.6418 -8067.6418 -8151.4449 -8151.4449 324.3265 324.3265 23397.429 23397.429 -898.51408 -898.51408 4000 -8066.5247 -8066.5247 -8154.0573 -8154.0573 338.76001 338.76001 23352.225 23352.225 2803.5237 2803.5237 Loop time of 23.3329 on 1 procs for 1000 steps with 2000 atoms Performance: 3.703 ns/day, 6.481 hours/ns, 42.858 timesteps/s 27.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 | 22.873 | 22.873 | 22.873 | 0.0 | 98.03 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.14796 | 0.14796 | 0.14796 | 0.0 | 0.63 Output | 3.6001e-05 | 3.6001e-05 | 3.6001e-05 | 0.0 | 0.00 Modify | 0.2683 | 0.2683 | 0.2683 | 0.0 | 1.15 Other | | 0.04384 | | | 0.19 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5440 ave 5440 max 5440 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 278210 ave 278210 max 278210 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 278210 Ave neighs/atom = 139.105 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 Per MPI rank memory allocation (min/avg/max) = 4.195 | 4.195 | 4.195 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 4000 -8066.5247 -8066.5247 -8154.0573 -8154.0573 338.76001 338.76001 23352.225 23352.225 2803.5237 2803.5237 5000 -8068.1393 -8068.1393 -8155.4632 -8155.4632 337.95223 337.95223 23435.806 23435.806 -3212.1205 -3212.1205 Loop time of 22.2354 on 1 procs for 1000 steps with 2000 atoms Performance: 3.886 ns/day, 6.176 hours/ns, 44.973 timesteps/s 29.2% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 21.816 | 21.816 | 21.816 | 0.0 | 98.12 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.1277 | 0.1277 | 0.1277 | 0.0 | 0.57 Output | 2.408e-05 | 2.408e-05 | 2.408e-05 | 0.0 | 0.00 Modify | 0.23698 | 0.23698 | 0.23698 | 0.0 | 1.07 Other | | 0.0544 | | | 0.24 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5451 ave 5451 max 5451 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 278390 ave 278390 max 278390 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 278390 Ave neighs/atom = 139.195 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.563131970498, Press = -122.379388963097 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.195 | 4.195 | 4.195 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 5000 -8068.1393 -8068.1393 -8155.4632 -8155.4632 337.95223 337.95223 23435.806 23435.806 -3212.1205 -3212.1205 6000 -8065.912 -8065.912 -8152.29 -8152.29 334.29179 334.29179 23333.255 23333.255 3413.3052 3413.3052 Loop time of 23.5438 on 1 procs for 1000 steps with 2000 atoms Performance: 3.670 ns/day, 6.540 hours/ns, 42.474 timesteps/s 27.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 | 23.061 | 23.061 | 23.061 | 0.0 | 97.95 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.11657 | 0.11657 | 0.11657 | 0.0 | 0.50 Output | 2.3127e-05 | 2.3127e-05 | 2.3127e-05 | 0.0 | 0.00 Modify | 0.32241 | 0.32241 | 0.32241 | 0.0 | 1.37 Other | | 0.04419 | | | 0.19 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5453 ave 5453 max 5453 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 277062 ave 277062 max 277062 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 277062 Ave neighs/atom = 138.531 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.358239076938, Press = -16.4334054160244 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.195 | 4.195 | 4.195 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 6000 -8065.912 -8065.912 -8152.29 -8152.29 334.29179 334.29179 23333.255 23333.255 3413.3052 3413.3052 7000 -8064.0742 -8064.0742 -8147.9904 -8147.9904 324.76449 324.76449 23412.496 23412.496 -1423.3399 -1423.3399 Loop time of 24.2221 on 1 procs for 1000 steps with 2000 atoms Performance: 3.567 ns/day, 6.728 hours/ns, 41.285 timesteps/s 28.4% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 23.422 | 23.422 | 23.422 | 0.0 | 96.70 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.21005 | 0.21005 | 0.21005 | 0.0 | 0.87 Output | 2.8849e-05 | 2.8849e-05 | 2.8849e-05 | 0.0 | 0.00 Modify | 0.5441 | 0.5441 | 0.5441 | 0.0 | 2.25 Other | | 0.04564 | | | 0.19 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5439 ave 5439 max 5439 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 279148 ave 279148 max 279148 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 279148 Ave neighs/atom = 139.574 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.39885414634, Press = 6.35554218370662 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.195 | 4.195 | 4.195 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 7000 -8064.0742 -8064.0742 -8147.9904 -8147.9904 324.76449 324.76449 23412.496 23412.496 -1423.3399 -1423.3399 8000 -8068.6485 -8068.6485 -8152.047 -8152.047 322.7609 322.7609 23382.01 23382.01 172.51458 172.51458 Loop time of 23.6409 on 1 procs for 1000 steps with 2000 atoms Performance: 3.655 ns/day, 6.567 hours/ns, 42.300 timesteps/s 28.5% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 23.073 | 23.073 | 23.073 | 0.0 | 97.60 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.17867 | 0.17867 | 0.17867 | 0.0 | 0.76 Output | 3.0994e-05 | 3.0994e-05 | 3.0994e-05 | 0.0 | 0.00 Modify | 0.35505 | 0.35505 | 0.35505 | 0.0 | 1.50 Other | | 0.03454 | | | 0.15 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5432 ave 5432 max 5432 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 277930 ave 277930 max 277930 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 277930 Ave neighs/atom = 138.965 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.381578089578, Press = 2.91509633588534 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.195 | 4.195 | 4.195 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 8000 -8068.6485 -8068.6485 -8152.047 -8152.047 322.7609 322.7609 23382.01 23382.01 172.51458 172.51458 9000 -8061.8305 -8061.8305 -8149.3242 -8149.3242 338.6093 338.6093 23409.794 23409.794 -945.35336 -945.35336 Loop time of 24.5101 on 1 procs for 1000 steps with 2000 atoms Performance: 3.525 ns/day, 6.808 hours/ns, 40.799 timesteps/s 27.1% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 23.973 | 23.973 | 23.973 | 0.0 | 97.81 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.14831 | 0.14831 | 0.14831 | 0.0 | 0.61 Output | 2.408e-05 | 2.408e-05 | 2.408e-05 | 0.0 | 0.00 Modify | 0.35417 | 0.35417 | 0.35417 | 0.0 | 1.45 Other | | 0.03471 | | | 0.14 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5444 ave 5444 max 5444 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 278398 ave 278398 max 278398 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 278398 Ave neighs/atom = 139.199 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.839557580917, Press = 0.0960758961706994 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.195 | 4.195 | 4.195 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 9000 -8061.8305 -8061.8305 -8149.3242 -8149.3242 338.6093 338.6093 23409.794 23409.794 -945.35336 -945.35336 10000 -8069.0222 -8069.0222 -8153.7929 -8153.7929 328.07136 328.07136 23370.91 23370.91 833.29881 833.29881 Loop time of 22.4645 on 1 procs for 1000 steps with 2000 atoms Performance: 3.846 ns/day, 6.240 hours/ns, 44.515 timesteps/s 29.3% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 21.93 | 21.93 | 21.93 | 0.0 | 97.62 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.18895 | 0.18895 | 0.18895 | 0.0 | 0.84 Output | 5.1975e-05 | 5.1975e-05 | 5.1975e-05 | 0.0 | 0.00 Modify | 0.33133 | 0.33133 | 0.33133 | 0.0 | 1.47 Other | | 0.01467 | | | 0.07 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5443 ave 5443 max 5443 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 277796 ave 277796 max 277796 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 277796 Ave neighs/atom = 138.898 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 = 334.023327708374, Press = -0.925202090701491 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.195 | 4.195 | 4.195 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 10000 -8069.0222 -8069.0222 -8153.7929 -8153.7929 328.07136 328.07136 23370.91 23370.91 833.29881 833.29881 11000 -8067.289 -8067.289 -8150.8598 -8150.8598 323.42732 323.42732 23432.61 23432.61 -2774.9529 -2774.9529 Loop time of 21.7187 on 1 procs for 1000 steps with 2000 atoms Performance: 3.978 ns/day, 6.033 hours/ns, 46.043 timesteps/s 30.3% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 21.221 | 21.221 | 21.221 | 0.0 | 97.71 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.10812 | 0.10812 | 0.10812 | 0.0 | 0.50 Output | 3.0041e-05 | 3.0041e-05 | 3.0041e-05 | 0.0 | 0.00 Modify | 0.31474 | 0.31474 | 0.31474 | 0.0 | 1.45 Other | | 0.0745 | | | 0.34 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5462 ave 5462 max 5462 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 277996 ave 277996 max 277996 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 277996 Ave neighs/atom = 138.998 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.83224950256, Press = 4.00455935526068 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.195 | 4.195 | 4.195 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 11000 -8067.289 -8067.289 -8150.8598 -8150.8598 323.42732 323.42732 23432.61 23432.61 -2774.9529 -2774.9529 12000 -8070.681 -8070.681 -8157.3819 -8157.3819 335.54149 335.54149 23353.505 23353.505 1893.9116 1893.9116 Loop time of 21.3266 on 1 procs for 1000 steps with 2000 atoms Performance: 4.051 ns/day, 5.924 hours/ns, 46.890 timesteps/s 30.7% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 20.749 | 20.749 | 20.749 | 0.0 | 97.29 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.27794 | 0.27794 | 0.27794 | 0.0 | 1.30 Output | 3.0041e-05 | 3.0041e-05 | 3.0041e-05 | 0.0 | 0.00 Modify | 0.22405 | 0.22405 | 0.22405 | 0.0 | 1.05 Other | | 0.07546 | | | 0.35 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5444 ave 5444 max 5444 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 277900 ave 277900 max 277900 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 277900 Ave neighs/atom = 138.95 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.290137648472, Press = 0.346235014439888 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.195 | 4.195 | 4.195 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 12000 -8070.681 -8070.681 -8157.3819 -8157.3819 335.54149 335.54149 23353.505 23353.505 1893.9116 1893.9116 13000 -8066.7267 -8066.7267 -8153.1562 -8153.1562 334.49107 334.49107 23421.652 23421.652 -2156.2969 -2156.2969 Loop time of 21.5638 on 1 procs for 1000 steps with 2000 atoms Performance: 4.007 ns/day, 5.990 hours/ns, 46.374 timesteps/s 30.4% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 21.009 | 21.009 | 21.009 | 0.0 | 97.43 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.11813 | 0.11813 | 0.11813 | 0.0 | 0.55 Output | 3.0041e-05 | 3.0041e-05 | 3.0041e-05 | 0.0 | 0.00 Modify | 0.37215 | 0.37215 | 0.37215 | 0.0 | 1.73 Other | | 0.06452 | | | 0.30 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5458 ave 5458 max 5458 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 277730 ave 277730 max 277730 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 277730 Ave neighs/atom = 138.865 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.182390070904, Press = 3.11332399026088 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.195 | 4.195 | 4.195 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 13000 -8066.7267 -8066.7267 -8153.1562 -8153.1562 334.49107 334.49107 23421.652 23421.652 -2156.2969 -2156.2969 14000 -8070.5278 -8070.5278 -8155.8324 -8155.8324 330.13773 330.13773 23340.278 23340.278 2694.653 2694.653 Loop time of 21.3936 on 1 procs for 1000 steps with 2000 atoms Performance: 4.039 ns/day, 5.943 hours/ns, 46.743 timesteps/s 30.3% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 21.059 | 21.059 | 21.059 | 0.0 | 98.44 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.12768 | 0.12768 | 0.12768 | 0.0 | 0.60 Output | 2.9802e-05 | 2.9802e-05 | 2.9802e-05 | 0.0 | 0.00 Modify | 0.14275 | 0.14275 | 0.14275 | 0.0 | 0.67 Other | | 0.06413 | | | 0.30 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5445 ave 5445 max 5445 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 277708 ave 277708 max 277708 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 277708 Ave neighs/atom = 138.854 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.899248692738, Press = 0.989654099597374 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.195 | 4.195 | 4.195 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 14000 -8070.5278 -8070.5278 -8155.8324 -8155.8324 330.13773 330.13773 23340.278 23340.278 2694.653 2694.653 15000 -8066.7478 -8066.7478 -8151.2611 -8151.2611 327.07523 327.07523 23427.646 23427.646 -2827.7543 -2827.7543 Loop time of 20.6864 on 1 procs for 1000 steps with 2000 atoms Performance: 4.177 ns/day, 5.746 hours/ns, 48.341 timesteps/s 31.7% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 20.248 | 20.248 | 20.248 | 0.0 | 97.88 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.15784 | 0.15784 | 0.15784 | 0.0 | 0.76 Output | 2.8849e-05 | 2.8849e-05 | 2.8849e-05 | 0.0 | 0.00 Modify | 0.26363 | 0.26363 | 0.26363 | 0.0 | 1.27 Other | | 0.01686 | | | 0.08 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5446 ave 5446 max 5446 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 278170 ave 278170 max 278170 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 278170 Ave neighs/atom = 139.085 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.773662970965, Press = 2.74841909576337 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.195 | 4.195 | 4.195 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 15000 -8066.7478 -8066.7478 -8151.2611 -8151.2611 327.07523 327.07523 23427.646 23427.646 -2827.7543 -2827.7543 16000 -8070.4693 -8070.4693 -8155.1433 -8155.1433 327.69702 327.69702 23339.111 23339.111 3311.5944 3311.5944 Loop time of 20.5061 on 1 procs for 1000 steps with 2000 atoms Performance: 4.213 ns/day, 5.696 hours/ns, 48.766 timesteps/s 32.1% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 20.111 | 20.111 | 20.111 | 0.0 | 98.07 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.077715 | 0.077715 | 0.077715 | 0.0 | 0.38 Output | 2.9087e-05 | 2.9087e-05 | 2.9087e-05 | 0.0 | 0.00 Modify | 0.24338 | 0.24338 | 0.24338 | 0.0 | 1.19 Other | | 0.07434 | | | 0.36 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5437 ave 5437 max 5437 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 277808 ave 277808 max 277808 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 277808 Ave neighs/atom = 138.904 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.596470889198, Press = -3.1429889131338 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.195 | 4.195 | 4.195 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 16000 -8070.4693 -8070.4693 -8155.1433 -8155.1433 327.69702 327.69702 23339.111 23339.111 3311.5944 3311.5944 17000 -8067.6533 -8067.6533 -8153.6632 -8153.6632 332.86712 332.86712 23410.093 23410.093 -1405.8301 -1405.8301 Loop time of 19.0362 on 1 procs for 1000 steps with 2000 atoms Performance: 4.539 ns/day, 5.288 hours/ns, 52.532 timesteps/s 34.4% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 18.671 | 18.671 | 18.671 | 0.0 | 98.08 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.05751 | 0.05751 | 0.05751 | 0.0 | 0.30 Output | 3.0994e-05 | 3.0994e-05 | 3.0994e-05 | 0.0 | 0.00 Modify | 0.29375 | 0.29375 | 0.29375 | 0.0 | 1.54 Other | | 0.0142 | | | 0.07 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5442 ave 5442 max 5442 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 278578 ave 278578 max 278578 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 278578 Ave neighs/atom = 139.289 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.775856958394, Press = 5.02610241728431 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.195 | 4.195 | 4.195 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 17000 -8067.6533 -8067.6533 -8153.6632 -8153.6632 332.86712 332.86712 23410.093 23410.093 -1405.8301 -1405.8301 18000 -8067.1879 -8067.1879 -8152.4863 -8152.4863 330.11368 330.11368 23364.183 23364.183 1612.6207 1612.6207 Loop time of 18.988 on 1 procs for 1000 steps with 2000 atoms Performance: 4.550 ns/day, 5.274 hours/ns, 52.665 timesteps/s 35.4% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 18.554 | 18.554 | 18.554 | 0.0 | 97.72 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.079428 | 0.079428 | 0.079428 | 0.0 | 0.42 Output | 2.8849e-05 | 2.8849e-05 | 2.8849e-05 | 0.0 | 0.00 Modify | 0.23154 | 0.23154 | 0.23154 | 0.0 | 1.22 Other | | 0.1227 | | | 0.65 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5448 ave 5448 max 5448 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 277544 ave 277544 max 277544 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 277544 Ave neighs/atom = 138.772 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.885641320938, Press = -0.847749903120894 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.195 | 4.195 | 4.195 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 18000 -8067.1879 -8067.1879 -8152.4863 -8152.4863 330.11368 330.11368 23364.183 23364.183 1612.6207 1612.6207 19000 -8066.78 -8066.78 -8151.5876 -8151.5876 328.21412 328.21412 23403.395 23403.395 -1071.5975 -1071.5975 Loop time of 20.8029 on 1 procs for 1000 steps with 2000 atoms Performance: 4.153 ns/day, 5.779 hours/ns, 48.070 timesteps/s 34.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 | 20.178 | 20.178 | 20.178 | 0.0 | 97.00 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.22141 | 0.22141 | 0.22141 | 0.0 | 1.06 Output | 3.0041e-05 | 3.0041e-05 | 3.0041e-05 | 0.0 | 0.00 Modify | 0.35808 | 0.35808 | 0.35808 | 0.0 | 1.72 Other | | 0.04539 | | | 0.22 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5461 ave 5461 max 5461 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 278054 ave 278054 max 278054 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 278054 Ave neighs/atom = 139.027 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.846209213741, Press = 1.99899555575496 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.195 | 4.195 | 4.195 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 19000 -8066.78 -8066.78 -8151.5876 -8151.5876 328.21412 328.21412 23403.395 23403.395 -1071.5975 -1071.5975 20000 -8070.7772 -8070.7772 -8155.4622 -8155.4622 327.73995 327.73995 23366.275 23366.275 1193.691 1193.691 Loop time of 20.1083 on 1 procs for 1000 steps with 2000 atoms Performance: 4.297 ns/day, 5.586 hours/ns, 49.731 timesteps/s 34.2% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 19.515 | 19.515 | 19.515 | 0.0 | 97.05 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.18065 | 0.18065 | 0.18065 | 0.0 | 0.90 Output | 3.0041e-05 | 3.0041e-05 | 3.0041e-05 | 0.0 | 0.00 Modify | 0.37744 | 0.37744 | 0.37744 | 0.0 | 1.88 Other | | 0.03475 | | | 0.17 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5450 ave 5450 max 5450 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 277908 ave 277908 max 277908 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 277908 Ave neighs/atom = 138.954 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.7750479749, Press = -0.481349636417936 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.195 | 4.195 | 4.195 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 20000 -8070.7772 -8070.7772 -8155.4622 -8155.4622 327.73995 327.73995 23366.275 23366.275 1193.691 1193.691 21000 -8066.8257 -8066.8257 -8153.3809 -8153.3809 334.97752 334.97752 23421.106 23421.106 -2082.4653 -2082.4653 Loop time of 19.6703 on 1 procs for 1000 steps with 2000 atoms Performance: 4.392 ns/day, 5.464 hours/ns, 50.838 timesteps/s 34.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 | 19.357 | 19.357 | 19.357 | 0.0 | 98.41 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.062224 | 0.062224 | 0.062224 | 0.0 | 0.32 Output | 2.1935e-05 | 2.1935e-05 | 2.1935e-05 | 0.0 | 0.00 Modify | 0.23606 | 0.23606 | 0.23606 | 0.0 | 1.20 Other | | 0.01488 | | | 0.08 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5450 ave 5450 max 5450 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 277954 ave 277954 max 277954 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 277954 Ave neighs/atom = 138.977 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.620414212998, Press = 2.23126141922721 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.195 | 4.195 | 4.195 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 21000 -8066.8257 -8066.8257 -8153.3809 -8153.3809 334.97752 334.97752 23421.106 23421.106 -2082.4653 -2082.4653 22000 -8070.423 -8070.423 -8152.7636 -8152.7636 318.66658 318.66658 23331.198 23331.198 3804.6126 3804.6126 Loop time of 17.3377 on 1 procs for 1000 steps with 2000 atoms Performance: 4.983 ns/day, 4.816 hours/ns, 57.678 timesteps/s 38.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 | 16.95 | 16.95 | 16.95 | 0.0 | 97.76 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.079267 | 0.079267 | 0.079267 | 0.0 | 0.46 Output | 3.4809e-05 | 3.4809e-05 | 3.4809e-05 | 0.0 | 0.00 Modify | 0.29431 | 0.29431 | 0.29431 | 0.0 | 1.70 Other | | 0.01408 | | | 0.08 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5433 ave 5433 max 5433 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 277554 ave 277554 max 277554 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 277554 Ave neighs/atom = 138.777 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.573594379102, Press = -2.3823390020138 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.195 | 4.195 | 4.195 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 22000 -8070.423 -8070.423 -8152.7636 -8152.7636 318.66658 318.66658 23331.198 23331.198 3804.6126 3804.6126 23000 -8061.4578 -8061.4578 -8150.4552 -8150.4552 344.42928 344.42928 23437.588 23437.588 -2836.4598 -2836.4598 Loop time of 17.3454 on 1 procs for 1000 steps with 2000 atoms Performance: 4.981 ns/day, 4.818 hours/ns, 57.652 timesteps/s 38.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 | 16.976 | 16.976 | 16.976 | 0.0 | 97.87 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.099416 | 0.099416 | 0.099416 | 0.0 | 0.57 Output | 2.3842e-05 | 2.3842e-05 | 2.3842e-05 | 0.0 | 0.00 Modify | 0.21571 | 0.21571 | 0.21571 | 0.0 | 1.24 Other | | 0.05445 | | | 0.31 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5436 ave 5436 max 5436 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 278612 ave 278612 max 278612 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 278612 Ave neighs/atom = 139.306 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.703129745073, Press = 1.7653280614047 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.195 | 4.195 | 4.195 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 23000 -8061.4578 -8061.4578 -8150.4552 -8150.4552 344.42928 344.42928 23437.588 23437.588 -2836.4598 -2836.4598 24000 -8067.8275 -8067.8275 -8153.9318 -8153.9318 333.23274 333.23274 23358.137 23358.137 2137.0125 2137.0125 Loop time of 18.3638 on 1 procs for 1000 steps with 2000 atoms Performance: 4.705 ns/day, 5.101 hours/ns, 54.455 timesteps/s 39.2% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 18.025 | 18.025 | 18.025 | 0.0 | 98.16 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.082345 | 0.082345 | 0.082345 | 0.0 | 0.45 Output | 3.0994e-05 | 3.0994e-05 | 3.0994e-05 | 0.0 | 0.00 Modify | 0.22146 | 0.22146 | 0.22146 | 0.0 | 1.21 Other | | 0.03491 | | | 0.19 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5436 ave 5436 max 5436 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 277592 ave 277592 max 277592 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 277592 Ave neighs/atom = 138.796 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.842317251828, Press = 0.158552255543147 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.195 | 4.195 | 4.195 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 24000 -8067.8275 -8067.8275 -8153.9318 -8153.9318 333.23274 333.23274 23358.137 23358.137 2137.0125 2137.0125 25000 -8066.4762 -8066.4762 -8155.7154 -8155.7154 345.36519 345.36519 23405.046 23405.046 -1076.2518 -1076.2518 Loop time of 17.3935 on 1 procs for 1000 steps with 2000 atoms Performance: 4.967 ns/day, 4.832 hours/ns, 57.493 timesteps/s 37.5% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 16.997 | 16.997 | 16.997 | 0.0 | 97.72 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.088246 | 0.088246 | 0.088246 | 0.0 | 0.51 Output | 2.9087e-05 | 2.9087e-05 | 2.9087e-05 | 0.0 | 0.00 Modify | 0.25391 | 0.25391 | 0.25391 | 0.0 | 1.46 Other | | 0.05475 | | | 0.31 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5450 ave 5450 max 5450 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 278374 ave 278374 max 278374 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 278374 Ave neighs/atom = 139.187 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.844734517742, Press = 2.17552047799384 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.195 | 4.195 | 4.195 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 25000 -8066.4762 -8066.4762 -8155.7154 -8155.7154 345.36519 345.36519 23405.046 23405.046 -1076.2518 -1076.2518 26000 -8069.5135 -8069.5135 -8155.5043 -8155.5043 332.79321 332.79321 23386.648 23386.648 -138.51553 -138.51553 Loop time of 18.5912 on 1 procs for 1000 steps with 2000 atoms Performance: 4.647 ns/day, 5.164 hours/ns, 53.789 timesteps/s 34.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 | 18.025 | 18.025 | 18.025 | 0.0 | 96.95 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.14071 | 0.14071 | 0.14071 | 0.0 | 0.76 Output | 3.1233e-05 | 3.1233e-05 | 3.1233e-05 | 0.0 | 0.00 Modify | 0.41171 | 0.41171 | 0.41171 | 0.0 | 2.21 Other | | 0.01386 | | | 0.07 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5428 ave 5428 max 5428 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 277602 ave 277602 max 277602 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 277602 Ave neighs/atom = 138.801 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.818176000126, Press = -0.634799198003051 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.195 | 4.195 | 4.195 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 26000 -8069.5135 -8069.5135 -8155.5043 -8155.5043 332.79321 332.79321 23386.648 23386.648 -138.51553 -138.51553 27000 -8065.0508 -8065.0508 -8150.9453 -8150.9453 332.42051 332.42051 23406.835 23406.835 -975.37548 -975.37548 Loop time of 20.8707 on 1 procs for 1000 steps with 2000 atoms Performance: 4.140 ns/day, 5.797 hours/ns, 47.914 timesteps/s 30.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 | 20.306 | 20.306 | 20.306 | 0.0 | 97.29 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.107 | 0.107 | 0.107 | 0.0 | 0.51 Output | 3.0041e-05 | 3.0041e-05 | 3.0041e-05 | 0.0 | 0.00 Modify | 0.41367 | 0.41367 | 0.41367 | 0.0 | 1.98 Other | | 0.04403 | | | 0.21 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5445 ave 5445 max 5445 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 277754 ave 277754 max 277754 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 277754 Ave neighs/atom = 138.877 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.839330401841, Press = 0.564810485911629 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.195 | 4.195 | 4.195 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 27000 -8065.0508 -8065.0508 -8150.9453 -8150.9453 332.42051 332.42051 23406.835 23406.835 -975.37548 -975.37548 28000 -8067.3357 -8067.3357 -8154.4793 -8154.4793 337.2545 337.2545 23366.841 23366.841 1448.7901 1448.7901 Loop time of 20.8847 on 1 procs for 1000 steps with 2000 atoms Performance: 4.137 ns/day, 5.801 hours/ns, 47.882 timesteps/s 31.2% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 20.48 | 20.48 | 20.48 | 0.0 | 98.06 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.097747 | 0.097747 | 0.097747 | 0.0 | 0.47 Output | 2.9802e-05 | 2.9802e-05 | 2.9802e-05 | 0.0 | 0.00 Modify | 0.27326 | 0.27326 | 0.27326 | 0.0 | 1.31 Other | | 0.03389 | | | 0.16 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5447 ave 5447 max 5447 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 277898 ave 277898 max 277898 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 277898 Ave neighs/atom = 138.949 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.892978290645, Press = -0.287302356614405 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.195 | 4.195 | 4.195 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 28000 -8067.3357 -8067.3357 -8154.4793 -8154.4793 337.2545 337.2545 23366.841 23366.841 1448.7901 1448.7901 29000 -8066.9201 -8066.9201 -8151.8734 -8151.8734 328.77803 328.77803 23422.427 23422.427 -2062.0823 -2062.0823 Loop time of 19.8534 on 1 procs for 1000 steps with 2000 atoms Performance: 4.352 ns/day, 5.515 hours/ns, 50.369 timesteps/s 32.5% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 19.402 | 19.402 | 19.402 | 0.0 | 97.73 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.1169 | 0.1169 | 0.1169 | 0.0 | 0.59 Output | 2.8849e-05 | 2.8849e-05 | 2.8849e-05 | 0.0 | 0.00 Modify | 0.25991 | 0.25991 | 0.25991 | 0.0 | 1.31 Other | | 0.07413 | | | 0.37 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 5448 ave 5448 max 5448 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 278462 ave 278462 max 278462 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 278462 Ave neighs/atom = 139.231 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 23387.7274656289 print "LAMMPS calculation completed" LAMMPS calculation completed quit 0