# 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 3.1399998813867573*${_u_distance} variable latticeconst_converted equal 3.1399998813867573*1 lattice bcc ${latticeconst_converted} lattice bcc 3.13999988138676 Lattice spacing in x,y,z = 3.14 3.14 3.14 region simbox block 0 10 0 10 0 10 units lattice create_box 1 simbox Created orthogonal box = (0 0 0) to (31.4 31.4 31.4) 1 by 1 by 1 MPI processor grid create_atoms 1 box Created 2000 atoms create_atoms CPU = 0.0102842 secs variable mass_converted equal 183.84*${_u_mass} variable mass_converted equal 183.84*1 # specify which KIM Model to use pair_style kim EAM_Dynamo_BonnyGrigorev1Terentyev_2017_W__MO_234187151804_000 pair_coeff * * W mass 1 ${mass_converted} mass 1 183.84 # initial volume variable v equal vol # assign formula variable V0 equal ${v} # evaluate initial value variable V0 equal 30959.1404915628 variable V0_metal equal ${V0}/(${_u_distance}*${_u_distance}*${_u_distance}) variable V0_metal equal 30959.1404915628/(${_u_distance}*${_u_distance}*${_u_distance}) variable V0_metal equal 30959.1404915628/(1*${_u_distance}*${_u_distance}) variable V0_metal equal 30959.1404915628/(1*1*${_u_distance}) variable V0_metal equal 30959.1404915628/(1*1*1) variable V0_metal_times1000 equal ${V0_metal}*1000 variable V0_metal_times1000 equal 30959.1404915628*1000 print "Initial system volume: ${V0_metal} Angstroms^3" Initial system volume: 30959.1404915628 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 293.15*${_u_temperature} variable temp_converted equal 293.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 293.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 293.15 ${temp_converted} ${Tdamp_converted} iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 293.15 293.15 ${Tdamp_converted} iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 293.15 293.15 0.1 iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 293.15 293.15 0.1 iso 0 ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 293.15 293.15 0.1 iso 0 0 ${Pdamp_converted} fix ensemble all npt temp 293.15 293.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 "293.15 - 0.2" variable T_up equal "293.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.46044 ghost atom cutoff = 7.46044 binsize = 3.73022, bins = 9 9 9 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 7.46044 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 4.133 | 4.133 | 4.133 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 0 -17724.213 -17724.213 -17799.96 -17799.96 293.15 293.15 30959.14 30959.14 2613.352 2613.352 1000 -17665.141 -17665.141 -17738.194 -17738.194 282.72303 282.72303 31204.817 31204.817 -31.732168 -31.732168 Loop time of 13.2093 on 1 procs for 1000 steps with 2000 atoms Performance: 6.541 ns/day, 3.669 hours/ns, 75.704 timesteps/s 37.6% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 12.866 | 12.866 | 12.866 | 0.0 | 97.40 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.11537 | 0.11537 | 0.11537 | 0.0 | 0.87 Output | 3.7909e-05 | 3.7909e-05 | 3.7909e-05 | 0.0 | 0.00 Modify | 0.19638 | 0.19638 | 0.19638 | 0.0 | 1.49 Other | | 0.03144 | | | 0.24 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 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: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 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.133 | 4.133 | 4.133 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 1000 -17665.141 -17665.141 -17738.194 -17738.194 282.72303 282.72303 31204.817 31204.817 -31.732168 -31.732168 2000 -17656.937 -17656.937 -17735.392 -17735.392 303.6296 303.6296 31214.173 31214.173 361.69248 361.69248 Loop time of 13.5644 on 1 procs for 1000 steps with 2000 atoms Performance: 6.370 ns/day, 3.768 hours/ns, 73.723 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 | 13.322 | 13.322 | 13.322 | 0.0 | 98.21 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.064979 | 0.064979 | 0.064979 | 0.0 | 0.48 Output | 3.4094e-05 | 3.4094e-05 | 3.4094e-05 | 0.0 | 0.00 Modify | 0.16593 | 0.16593 | 0.16593 | 0.0 | 1.22 Other | | 0.01119 | | | 0.08 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 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: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 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.133 | 4.133 | 4.133 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 2000 -17656.937 -17656.937 -17735.392 -17735.392 303.6296 303.6296 31214.173 31214.173 361.69248 361.69248 3000 -17661.702 -17661.702 -17736.09 -17736.09 287.8903 287.8903 31196.447 31196.447 2187.1219 2187.1219 Loop time of 12.3617 on 1 procs for 1000 steps with 2000 atoms Performance: 6.989 ns/day, 3.434 hours/ns, 80.895 timesteps/s 41.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 | 12.079 | 12.079 | 12.079 | 0.0 | 97.71 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.044953 | 0.044953 | 0.044953 | 0.0 | 0.36 Output | 3.0994e-05 | 3.0994e-05 | 3.0994e-05 | 0.0 | 0.00 Modify | 0.18647 | 0.18647 | 0.18647 | 0.0 | 1.51 Other | | 0.0514 | | | 0.42 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 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: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 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.133 | 4.133 | 4.133 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 3000 -17661.702 -17661.702 -17736.09 -17736.09 287.8903 287.8903 31196.447 31196.447 2187.1219 2187.1219 4000 -17661.972 -17661.972 -17736.124 -17736.124 286.973 286.973 31213.546 31213.546 -938.01083 -938.01083 Loop time of 13.051 on 1 procs for 1000 steps with 2000 atoms Performance: 6.620 ns/day, 3.625 hours/ns, 76.622 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 | 12.709 | 12.709 | 12.709 | 0.0 | 97.38 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.074848 | 0.074848 | 0.074848 | 0.0 | 0.57 Output | 3.6955e-05 | 3.6955e-05 | 3.6955e-05 | 0.0 | 0.00 Modify | 0.23567 | 0.23567 | 0.23567 | 0.0 | 1.81 Other | | 0.0313 | | | 0.24 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 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: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 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.133 | 4.133 | 4.133 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 4000 -17661.972 -17661.972 -17736.124 -17736.124 286.973 286.973 31213.546 31213.546 -938.01083 -938.01083 5000 -17659.123 -17659.123 -17733.158 -17733.158 286.52294 286.52294 31203.721 31203.721 1296.7136 1296.7136 Loop time of 13.2841 on 1 procs for 1000 steps with 2000 atoms Performance: 6.504 ns/day, 3.690 hours/ns, 75.278 timesteps/s 38.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 | 12.941 | 12.941 | 12.941 | 0.0 | 97.42 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.10513 | 0.10513 | 0.10513 | 0.0 | 0.79 Output | 2.7895e-05 | 2.7895e-05 | 2.7895e-05 | 0.0 | 0.00 Modify | 0.18648 | 0.18648 | 0.18648 | 0.0 | 1.40 Other | | 0.05155 | | | 0.39 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 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: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 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 = 297.958513895458, Press = 165.390273812092 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.133 | 4.133 | 4.133 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 5000 -17659.123 -17659.123 -17733.158 -17733.158 286.52294 286.52294 31203.721 31203.721 1296.7136 1296.7136 6000 -17661.709 -17661.709 -17737.284 -17737.284 292.48235 292.48235 31210.571 31210.571 -966.16067 -966.16067 Loop time of 13.5448 on 1 procs for 1000 steps with 2000 atoms Performance: 6.379 ns/day, 3.762 hours/ns, 73.829 timesteps/s 37.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 | 13.218 | 13.218 | 13.218 | 0.0 | 97.59 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.044931 | 0.044931 | 0.044931 | 0.0 | 0.33 Output | 2.7895e-05 | 2.7895e-05 | 2.7895e-05 | 0.0 | 0.00 Modify | 0.25016 | 0.25016 | 0.25016 | 0.0 | 1.85 Other | | 0.03125 | | | 0.23 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 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: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 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 = 291.120470710331, Press = -48.783960678444 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.133 | 4.133 | 4.133 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 6000 -17661.709 -17661.709 -17737.284 -17737.284 292.48235 292.48235 31210.571 31210.571 -966.16067 -966.16067 7000 -17662.551 -17662.551 -17738.51 -17738.51 293.9693 293.9693 31169.412 31169.412 4464.0407 4464.0407 Loop time of 12.9461 on 1 procs for 1000 steps with 2000 atoms Performance: 6.674 ns/day, 3.596 hours/ns, 77.243 timesteps/s 39.6% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 12.721 | 12.721 | 12.721 | 0.0 | 98.26 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.02461 | 0.02461 | 0.02461 | 0.0 | 0.19 Output | 2.9087e-05 | 2.9087e-05 | 2.9087e-05 | 0.0 | 0.00 Modify | 0.14903 | 0.14903 | 0.14903 | 0.0 | 1.15 Other | | 0.05121 | | | 0.40 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 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: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 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 = 291.984047075918, Press = 0.979517422585126 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.133 | 4.133 | 4.133 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 7000 -17662.551 -17662.551 -17738.51 -17738.51 293.9693 293.9693 31169.412 31169.412 4464.0407 4464.0407 8000 -17660.499 -17660.499 -17735.446 -17735.446 290.05484 290.05484 31214.8 31214.8 -613.31354 -613.31354 Loop time of 12.7605 on 1 procs for 1000 steps with 2000 atoms Performance: 6.771 ns/day, 3.545 hours/ns, 78.367 timesteps/s 40.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 | 12.514 | 12.514 | 12.514 | 0.0 | 98.07 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.06528 | 0.06528 | 0.06528 | 0.0 | 0.51 Output | 2.9087e-05 | 2.9087e-05 | 2.9087e-05 | 0.0 | 0.00 Modify | 0.16952 | 0.16952 | 0.16952 | 0.0 | 1.33 Other | | 0.01128 | | | 0.09 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 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: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 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 = 292.212066472952, Press = 5.08813228687585 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.133 | 4.133 | 4.133 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 8000 -17660.499 -17660.499 -17735.446 -17735.446 290.05484 290.05484 31214.8 31214.8 -613.31354 -613.31354 9000 -17662.117 -17662.117 -17736.94 -17736.94 289.57532 289.57532 31194.022 31194.022 1689.8754 1689.8754 Loop time of 12.0057 on 1 procs for 1000 steps with 2000 atoms Performance: 7.197 ns/day, 3.335 hours/ns, 83.294 timesteps/s 42.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 | 11.74 | 11.74 | 11.74 | 0.0 | 97.79 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.10495 | 0.10495 | 0.10495 | 0.0 | 0.87 Output | 2.8849e-05 | 2.8849e-05 | 2.8849e-05 | 0.0 | 0.00 Modify | 0.1491 | 0.1491 | 0.1491 | 0.0 | 1.24 Other | | 0.01114 | | | 0.09 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 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: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 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 = 293.049355284524, Press = 3.91445736739094 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.133 | 4.133 | 4.133 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 9000 -17662.117 -17662.117 -17736.94 -17736.94 289.57532 289.57532 31194.022 31194.022 1689.8754 1689.8754 10000 -17657.866 -17657.866 -17733.835 -17733.835 294.00938 294.00938 31201.262 31201.262 2291.5455 2291.5455 Loop time of 11.7141 on 1 procs for 1000 steps with 2000 atoms Performance: 7.376 ns/day, 3.254 hours/ns, 85.368 timesteps/s 43.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 | 11.448 | 11.448 | 11.448 | 0.0 | 97.73 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.085743 | 0.085743 | 0.085743 | 0.0 | 0.73 Output | 4.6015e-05 | 4.6015e-05 | 4.6015e-05 | 0.0 | 0.00 Modify | 0.1489 | 0.1489 | 0.1489 | 0.0 | 1.27 Other | | 0.03117 | | | 0.27 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 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: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 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 = 293.380075648854, Press = -14.7761012998577 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.133 | 4.133 | 4.133 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 10000 -17657.866 -17657.866 -17733.835 -17733.835 294.00938 294.00938 31201.262 31201.262 2291.5455 2291.5455 11000 -17660.968 -17660.968 -17734.907 -17734.907 286.15169 286.15169 31194.168 31194.168 1840.2251 1840.2251 Loop time of 12.8979 on 1 procs for 1000 steps with 2000 atoms Performance: 6.699 ns/day, 3.583 hours/ns, 77.532 timesteps/s 40.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 | 12.569 | 12.569 | 12.569 | 0.0 | 97.45 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.10233 | 0.10233 | 0.10233 | 0.0 | 0.79 Output | 2.9087e-05 | 2.9087e-05 | 2.9087e-05 | 0.0 | 0.00 Modify | 0.1764 | 0.1764 | 0.1764 | 0.0 | 1.37 Other | | 0.04981 | | | 0.39 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 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: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 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 = 293.349393070714, Press = -10.7120330579527 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.133 | 4.133 | 4.133 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 11000 -17660.968 -17660.968 -17734.907 -17734.907 286.15169 286.15169 31194.168 31194.168 1840.2251 1840.2251 12000 -17659.116 -17659.116 -17735.569 -17735.569 295.88142 295.88142 31212.958 31212.958 -509.36546 -509.36546 Loop time of 12.2182 on 1 procs for 1000 steps with 2000 atoms Performance: 7.071 ns/day, 3.394 hours/ns, 81.845 timesteps/s 42.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 | 11.953 | 11.953 | 11.953 | 0.0 | 97.83 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.064896 | 0.064896 | 0.064896 | 0.0 | 0.53 Output | 2.8133e-05 | 2.8133e-05 | 2.8133e-05 | 0.0 | 0.00 Modify | 0.16913 | 0.16913 | 0.16913 | 0.0 | 1.38 Other | | 0.03126 | | | 0.26 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 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: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 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 = 293.496532289343, Press = -6.54902071667276 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.133 | 4.133 | 4.133 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 12000 -17659.116 -17659.116 -17735.569 -17735.569 295.88142 295.88142 31212.958 31212.958 -509.36546 -509.36546 13000 -17662.239 -17662.239 -17735.996 -17735.996 285.44858 285.44858 31198.982 31198.982 1163.4994 1163.4994 Loop time of 12.4784 on 1 procs for 1000 steps with 2000 atoms Performance: 6.924 ns/day, 3.466 hours/ns, 80.139 timesteps/s 41.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 | 12.24 | 12.24 | 12.24 | 0.0 | 98.09 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.035066 | 0.035066 | 0.035066 | 0.0 | 0.28 Output | 2.7895e-05 | 2.7895e-05 | 2.7895e-05 | 0.0 | 0.00 Modify | 0.19165 | 0.19165 | 0.19165 | 0.0 | 1.54 Other | | 0.01118 | | | 0.09 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 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: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 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 = 293.705292129094, Press = 2.81390889301672 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.133 | 4.133 | 4.133 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 13000 -17662.239 -17662.239 -17735.996 -17735.996 285.44858 285.44858 31198.982 31198.982 1163.4994 1163.4994 14000 -17660.135 -17660.135 -17735.823 -17735.823 292.9214 292.9214 31206.66 31206.66 524.98191 524.98191 Loop time of 11.0299 on 1 procs for 1000 steps with 2000 atoms Performance: 7.833 ns/day, 3.064 hours/ns, 90.663 timesteps/s 46.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 | 10.831 | 10.831 | 10.831 | 0.0 | 98.20 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.02461 | 0.02461 | 0.02461 | 0.0 | 0.22 Output | 2.9087e-05 | 2.9087e-05 | 2.9087e-05 | 0.0 | 0.00 Modify | 0.16317 | 0.16317 | 0.16317 | 0.0 | 1.48 Other | | 0.01121 | | | 0.10 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 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: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 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 = 293.950642970095, Press = -1.27734695768755 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.133 | 4.133 | 4.133 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 14000 -17660.135 -17660.135 -17735.823 -17735.823 292.9214 292.9214 31206.66 31206.66 524.98191 524.98191 15000 -17658.65 -17658.65 -17736.936 -17736.936 302.97572 302.97572 31191.704 31191.704 2390.9842 2390.9842 Loop time of 11.1909 on 1 procs for 1000 steps with 2000 atoms Performance: 7.721 ns/day, 3.109 hours/ns, 89.358 timesteps/s 45.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 | 10.884 | 10.884 | 10.884 | 0.0 | 97.26 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.024754 | 0.024754 | 0.024754 | 0.0 | 0.22 Output | 2.7895e-05 | 2.7895e-05 | 2.7895e-05 | 0.0 | 0.00 Modify | 0.23048 | 0.23048 | 0.23048 | 0.0 | 2.06 Other | | 0.05142 | | | 0.46 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 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: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 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 = 294.044589037801, Press = 3.13550315039238 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.133 | 4.133 | 4.133 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 15000 -17658.65 -17658.65 -17736.936 -17736.936 302.97572 302.97572 31191.704 31191.704 2390.9842 2390.9842 16000 -17665.659 -17665.659 -17739.595 -17739.595 286.1429 286.1429 31217.997 31217.997 -3308.7596 -3308.7596 Loop time of 11.0413 on 1 procs for 1000 steps with 2000 atoms Performance: 7.825 ns/day, 3.067 hours/ns, 90.569 timesteps/s 46.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 | 10.8 | 10.8 | 10.8 | 0.0 | 97.81 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.033948 | 0.033948 | 0.033948 | 0.0 | 0.31 Output | 2.1935e-05 | 2.1935e-05 | 2.1935e-05 | 0.0 | 0.00 Modify | 0.19663 | 0.19663 | 0.19663 | 0.0 | 1.78 Other | | 0.01111 | | | 0.10 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 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: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 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 = 293.856292200251, Press = -2.02477797940202 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.133 | 4.133 | 4.133 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 16000 -17665.659 -17665.659 -17739.595 -17739.595 286.1429 286.1429 31217.997 31217.997 -3308.7596 -3308.7596 17000 -17660.033 -17660.033 -17735.953 -17735.953 293.81703 293.81703 31197.586 31197.586 1793.4793 1793.4793 Loop time of 11.6862 on 1 procs for 1000 steps with 2000 atoms Performance: 7.393 ns/day, 3.246 hours/ns, 85.571 timesteps/s 43.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 | 11.443 | 11.443 | 11.443 | 0.0 | 97.92 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.044864 | 0.044864 | 0.044864 | 0.0 | 0.38 Output | 2.9087e-05 | 2.9087e-05 | 2.9087e-05 | 0.0 | 0.00 Modify | 0.17051 | 0.17051 | 0.17051 | 0.0 | 1.46 Other | | 0.02756 | | | 0.24 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 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: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 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 = 293.628004321773, Press = -1.12752615813362 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.133 | 4.133 | 4.133 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 17000 -17660.033 -17660.033 -17735.953 -17735.953 293.81703 293.81703 31197.586 31197.586 1793.4793 1793.4793 18000 -17661.343 -17661.343 -17737.988 -17737.988 296.62268 296.62268 31187.271 31187.271 1815.7695 1815.7695 Loop time of 11.0372 on 1 procs for 1000 steps with 2000 atoms Performance: 7.828 ns/day, 3.066 hours/ns, 90.602 timesteps/s 46.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 | 10.782 | 10.782 | 10.782 | 0.0 | 97.69 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.044416 | 0.044416 | 0.044416 | 0.0 | 0.40 Output | 2.6941e-05 | 2.6941e-05 | 2.6941e-05 | 0.0 | 0.00 Modify | 0.19958 | 0.19958 | 0.19958 | 0.0 | 1.81 Other | | 0.01121 | | | 0.10 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 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: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 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 = 293.347976157772, Press = -1.15611374578637 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.133 | 4.133 | 4.133 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 18000 -17661.343 -17661.343 -17737.988 -17737.988 296.62268 296.62268 31187.271 31187.271 1815.7695 1815.7695 19000 -17662.078 -17662.078 -17735.383 -17735.383 283.69543 283.69543 31197.626 31197.626 1043.5486 1043.5486 Loop time of 11.5936 on 1 procs for 1000 steps with 2000 atoms Performance: 7.452 ns/day, 3.220 hours/ns, 86.255 timesteps/s 44.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 | 11.288 | 11.288 | 11.288 | 0.0 | 97.37 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.064711 | 0.064711 | 0.064711 | 0.0 | 0.56 Output | 2.8133e-05 | 2.8133e-05 | 2.8133e-05 | 0.0 | 0.00 Modify | 0.2095 | 0.2095 | 0.2095 | 0.0 | 1.81 Other | | 0.03107 | | | 0.27 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 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: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 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 = 293.393461814289, Press = 0.582490781417774 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.133 | 4.133 | 4.133 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 19000 -17662.078 -17662.078 -17735.383 -17735.383 283.69543 283.69543 31197.626 31197.626 1043.5486 1043.5486 20000 -17655.928 -17655.928 -17733.322 -17733.322 299.52406 299.52406 31221.106 31221.106 -900.79441 -900.79441 Loop time of 11.1275 on 1 procs for 1000 steps with 2000 atoms Performance: 7.765 ns/day, 3.091 hours/ns, 89.868 timesteps/s 46.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 | 10.83 | 10.83 | 10.83 | 0.0 | 97.32 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.046171 | 0.046171 | 0.046171 | 0.0 | 0.41 Output | 2.7895e-05 | 2.7895e-05 | 2.7895e-05 | 0.0 | 0.00 Modify | 0.22041 | 0.22041 | 0.22041 | 0.0 | 1.98 Other | | 0.03111 | | | 0.28 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 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: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 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 = 293.474957835319, Press = -2.16940155685515 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.133 | 4.133 | 4.133 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 20000 -17655.928 -17655.928 -17733.322 -17733.322 299.52406 299.52406 31221.106 31221.106 -900.79441 -900.79441 21000 -17660.524 -17660.524 -17737.134 -17737.134 296.48829 296.48829 31202.067 31202.067 -102.55833 -102.55833 Loop time of 10.7266 on 1 procs for 1000 steps with 2000 atoms Performance: 8.055 ns/day, 2.980 hours/ns, 93.226 timesteps/s 47.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 | 10.541 | 10.541 | 10.541 | 0.0 | 98.27 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.024606 | 0.024606 | 0.024606 | 0.0 | 0.23 Output | 2.7895e-05 | 2.7895e-05 | 2.7895e-05 | 0.0 | 0.00 Modify | 0.12953 | 0.12953 | 0.12953 | 0.0 | 1.21 Other | | 0.03119 | | | 0.29 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 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: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 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 = 293.547032744316, Press = -5.43013687055479 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.133 | 4.133 | 4.133 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 21000 -17660.524 -17660.524 -17737.134 -17737.134 296.48829 296.48829 31202.067 31202.067 -102.55833 -102.55833 22000 -17661.925 -17661.925 -17736.377 -17736.377 288.13794 288.13794 31188.04 31188.04 2997.582 2997.582 Loop time of 11.1197 on 1 procs for 1000 steps with 2000 atoms Performance: 7.770 ns/day, 3.089 hours/ns, 89.931 timesteps/s 45.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 | 10.781 | 10.781 | 10.781 | 0.0 | 96.95 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.044651 | 0.044651 | 0.044651 | 0.0 | 0.40 Output | 2.0981e-05 | 2.0981e-05 | 2.0981e-05 | 0.0 | 0.00 Modify | 0.26307 | 0.26307 | 0.26307 | 0.0 | 2.37 Other | | 0.03134 | | | 0.28 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 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: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 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 = 293.58758940141, Press = -2.06588041630454 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.133 | 4.133 | 4.133 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 22000 -17661.925 -17661.925 -17736.377 -17736.377 288.13794 288.13794 31188.04 31188.04 2997.582 2997.582 23000 -17663.062 -17663.062 -17737.731 -17737.731 288.97671 288.97671 31190.179 31190.179 1714.1764 1714.1764 Loop time of 11.2192 on 1 procs for 1000 steps with 2000 atoms Performance: 7.701 ns/day, 3.116 hours/ns, 89.133 timesteps/s 45.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 | 10.937 | 10.937 | 10.937 | 0.0 | 97.48 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.065573 | 0.065573 | 0.065573 | 0.0 | 0.58 Output | 2.6941e-05 | 2.6941e-05 | 2.6941e-05 | 0.0 | 0.00 Modify | 0.16557 | 0.16557 | 0.16557 | 0.0 | 1.48 Other | | 0.05141 | | | 0.46 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 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: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 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 = 293.483552297213, Press = 2.74710784938171 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.133 | 4.133 | 4.133 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 23000 -17663.062 -17663.062 -17737.731 -17737.731 288.97671 288.97671 31190.179 31190.179 1714.1764 1714.1764 24000 -17659.643 -17659.643 -17734.978 -17734.978 291.55428 291.55428 31197.284 31197.284 1625.2095 1625.2095 Loop time of 11.4995 on 1 procs for 1000 steps with 2000 atoms Performance: 7.513 ns/day, 3.194 hours/ns, 86.960 timesteps/s 44.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 | 11.144 | 11.144 | 11.144 | 0.0 | 96.91 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.085066 | 0.085066 | 0.085066 | 0.0 | 0.74 Output | 2.718e-05 | 2.718e-05 | 2.718e-05 | 0.0 | 0.00 Modify | 0.24958 | 0.24958 | 0.24958 | 0.0 | 2.17 Other | | 0.02123 | | | 0.18 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 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: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 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 = 293.377285103144, Press = -1.62230521672032 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.133 | 4.133 | 4.133 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 24000 -17659.643 -17659.643 -17734.978 -17734.978 291.55428 291.55428 31197.284 31197.284 1625.2095 1625.2095 25000 -17662.026 -17662.026 -17737.359 -17737.359 291.54707 291.54707 31187.547 31187.547 1828.996 1828.996 Loop time of 11.5922 on 1 procs for 1000 steps with 2000 atoms Performance: 7.453 ns/day, 3.220 hours/ns, 86.265 timesteps/s 44.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 | 11.236 | 11.236 | 11.236 | 0.0 | 96.93 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.044803 | 0.044803 | 0.044803 | 0.0 | 0.39 Output | 2.7895e-05 | 2.7895e-05 | 2.7895e-05 | 0.0 | 0.00 Modify | 0.28017 | 0.28017 | 0.28017 | 0.0 | 2.42 Other | | 0.03134 | | | 0.27 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 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: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 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 = 293.252240459822, Press = -1.62709561934174 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.133 | 4.133 | 4.133 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 25000 -17662.026 -17662.026 -17737.359 -17737.359 291.54707 291.54707 31187.547 31187.547 1828.996 1828.996 26000 -17663.205 -17663.205 -17737.519 -17737.519 287.60529 287.60529 31192.197 31192.197 1808.894 1808.894 Loop time of 11.1098 on 1 procs for 1000 steps with 2000 atoms Performance: 7.777 ns/day, 3.086 hours/ns, 90.011 timesteps/s 45.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 | 10.832 | 10.832 | 10.832 | 0.0 | 97.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.044505 | 0.044505 | 0.044505 | 0.0 | 0.40 Output | 2.8133e-05 | 2.8133e-05 | 2.8133e-05 | 0.0 | 0.00 Modify | 0.22181 | 0.22181 | 0.22181 | 0.0 | 2.00 Other | | 0.01107 | | | 0.10 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 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: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 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 = 293.258854569093, Press = -0.34195333742879 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.133 | 4.133 | 4.133 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 26000 -17663.205 -17663.205 -17737.519 -17737.519 287.60529 287.60529 31192.197 31192.197 1808.894 1808.894 27000 -17658.678 -17658.678 -17734.529 -17734.529 293.54765 293.54765 31205.99 31205.99 1603.3942 1603.3942 Loop time of 11.4454 on 1 procs for 1000 steps with 2000 atoms Performance: 7.549 ns/day, 3.179 hours/ns, 87.372 timesteps/s 44.6% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 11.17 | 11.17 | 11.17 | 0.0 | 97.60 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.084702 | 0.084702 | 0.084702 | 0.0 | 0.74 Output | 4.5061e-05 | 4.5061e-05 | 4.5061e-05 | 0.0 | 0.00 Modify | 0.1789 | 0.1789 | 0.1789 | 0.0 | 1.56 Other | | 0.01123 | | | 0.10 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 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: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 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 = 293.22890034339, Press = -0.617786528779984 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.133 | 4.133 | 4.133 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 27000 -17658.678 -17658.678 -17734.529 -17734.529 293.54765 293.54765 31205.99 31205.99 1603.3942 1603.3942 28000 -17664.763 -17664.763 -17738.474 -17738.474 285.26998 285.26998 31216.233 31216.233 -2468.5116 -2468.5116 Loop time of 10.8142 on 1 procs for 1000 steps with 2000 atoms Performance: 7.989 ns/day, 3.004 hours/ns, 92.471 timesteps/s 46.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 | 10.509 | 10.509 | 10.509 | 0.0 | 97.18 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.044509 | 0.044509 | 0.044509 | 0.0 | 0.41 Output | 3.9101e-05 | 3.9101e-05 | 3.9101e-05 | 0.0 | 0.00 Modify | 0.22941 | 0.22941 | 0.22941 | 0.0 | 2.12 Other | | 0.0312 | | | 0.29 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 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: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 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 = 293.160572353412, Press = -1.17764603777529 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.133 | 4.133 | 4.133 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 28000 -17664.763 -17664.763 -17738.474 -17738.474 285.26998 285.26998 31216.233 31216.233 -2468.5116 -2468.5116 29000 -17661.377 -17661.377 -17737.481 -17737.481 294.53151 294.53151 31194.384 31194.384 796.81645 796.81645 Loop time of 10.5411 on 1 procs for 1000 steps with 2000 atoms Performance: 8.197 ns/day, 2.928 hours/ns, 94.867 timesteps/s 48.6% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 10.276 | 10.276 | 10.276 | 0.0 | 97.48 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.064777 | 0.064777 | 0.064777 | 0.0 | 0.61 Output | 2.9087e-05 | 2.9087e-05 | 2.9087e-05 | 0.0 | 0.00 Modify | 0.1694 | 0.1694 | 0.1694 | 0.0 | 1.61 Other | | 0.03122 | | | 0.30 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 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: 224000 ave 224000 max 224000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224000 Ave neighs/atom = 112 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 = 293.119391173008, Press = -1.68330311693074 next a jump SELF top variable a loop 2000 run 1000 Per MPI rank memory allocation (min/avg/max) = 4.133 | 4.133 | 4.133 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 29000 -17661.377 -17661.377 -17737.481 -17737.481 294.53151 294.53151 31194.384 31194.384 796.81645 796.81645 30000 -17658.125 -17658.125 -17738.13 -17738.13 309.62873 309.62873 31212.706 31212.706 -1564.1654 -1564.1654 Loop time of 11.0794 on 1 procs for 1000 steps with 2000 atoms Performance: 7.798 ns/day, 3.078 hours/ns, 90.257 timesteps/s 46.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 | 10.764 | 10.764 | 10.764 | 0.0 | 97.15 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.084792 | 0.084792 | 0.084792 | 0.0 | 0.77 Output | 2.2173e-05 | 2.2173e-05 | 2.2173e-05 | 0.0 | 0.00 Modify | 0.21959 | 0.21959 | 0.21959 | 0.0 | 1.98 Other | | 0.01095 | | | 0.10 Nlocal: 2000 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 4119 ave 4119 max 4119 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: 224002 ave 224002 max 224002 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 224002 Ave neighs/atom = 112.001 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_T293.15.out" else "print 'not_converged' file output/vol_T293.15.out" print '${V}' file output/vol_T293.15.out 31207.2859632104 print "LAMMPS calculation completed" LAMMPS calculation completed quit 0