# 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.1465589553117748*${_u_distance} variable latticeconst_converted equal 3.1465589553117748*1 lattice bcc ${latticeconst_converted} lattice bcc 3.14655895531177 Lattice spacing in x,y,z = 3.1465590 3.1465590 3.1465590 region simbox block 0 10 0 10 0 10 units lattice create_box 1 simbox Created orthogonal box = (0.0000000 0.0000000 0.0000000) to (31.465590 31.465590 31.465590) 1 by 1 by 1 MPI processor grid create_atoms 1 box Created 2000 atoms using lattice units in orthogonal box = (0.0000000 0.0000000 0.0000000) to (31.465590 31.465590 31.465590) create_atoms CPU = 0.002 seconds variable mass_converted equal 95.94*${_u_mass} variable mass_converted equal 95.94*1 kim_interactions Mo WARNING: 'kim_' has been renamed to 'kim '. Please update your input. kim interactions Mo #=== BEGIN kim interactions ================================== pair_style kim MEAM_LAMMPS_LeeBaskesKim_2001_Mo__MO_805823015127_000 pair_coeff * * Mo #=== END kim interactions ==================================== mass 1 ${mass_converted} mass 1 95.94 # initial volume variable v equal vol # assign formula variable V0 equal ${v} # evaluate initial value variable V0 equal 31153.5555569502 variable V0_metal equal ${V0}/(${_u_distance}*${_u_distance}*${_u_distance}) variable V0_metal equal 31153.5555569502/(${_u_distance}*${_u_distance}*${_u_distance}) variable V0_metal equal 31153.5555569502/(1*${_u_distance}*${_u_distance}) variable V0_metal equal 31153.5555569502/(1*1*${_u_distance}) variable V0_metal equal 31153.5555569502/(1*1*1) variable V0_metal_times1000 equal ${V0_metal}*1000 variable V0_metal_times1000 equal 31153.5555569502*1000 print "Initial system volume: ${V0_metal} Angstroms^3" Initial system volume: 31153.5555569502 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 273.15*${_u_temperature} variable temp_converted equal 273.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 273.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 273.15 ${temp_converted} ${Tdamp_converted} iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 273.15 273.15 ${Tdamp_converted} iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 273.15 273.15 0.1 iso ${press_converted} ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 273.15 273.15 0.1 iso 0 ${press_converted} ${Pdamp_converted} fix ensemble all npt temp 273.15 273.15 0.1 iso 0 0 ${Pdamp_converted} fix ensemble all npt temp 273.15 273.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 "273.15 - 0.2" variable T_up equal "273.15 + 0.2" variable P_low equal "0.0 - 0.2" variable P_up equal "0.0 + 0.2" # print to logfile every 1000 timesteps thermo_style custom step etotal v_etotal_metal pe v_pe_metal temp v_T_metal vol v_V_metal press v_P_metal thermo 1000 # Run a simulation for at most 2000*1000 timesteps. At each 1000th time step, check # whether the temperature and pressure have converged. If yes, break. label top variable a loop 2000 run 1000 CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE Your simulation uses code contributions which should be cited: - OpenKIM: https://doi.org/10.1007/s11837-011-0102-6 @Article{tadmor:elliott:2011, author = {E. B. Tadmor and R. S. Elliott and J. P. Sethna and R. E. Miller and C. A. Becker}, title = {The potential of atomistic simulations and the {K}nowledgebase of {I}nteratomic {M}odels}, journal = {{JOM}}, year = 2011, volume = 63, number = 17, pages = {17}, doi = {10.1007/s11837-011-0102-6} } - OpenKIM potential: https://openkim.org/cite/MO_805823015127_000#item-citation CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 0 -13549.42 -13549.42 -13620 -13620 273.15 273.15 31153.556 31153.556 2419.8966 2419.8966 1000 -13475.788 -13475.788 -13544.094 -13544.094 264.35103 264.35103 31239.995 31239.995 3850.6849 3850.6849 Loop time of 97.5499 on 1 procs for 1000 steps with 2000 atoms Performance: 0.886 ns/day, 27.097 hours/ns, 10.251 timesteps/s 99.9% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 97.054 | 97.054 | 97.054 | 0.0 | 99.49 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.077826 | 0.077826 | 0.077826 | 0.0 | 0.08 Output | 0.00020155 | 0.00020155 | 0.00020155 | 0.0 | 0.00 Modify | 0.35945 | 0.35945 | 0.35945 | 0.0 | 0.37 Other | | 0.05804 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.000000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 0, Press = 0 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 1000 -13475.788 -13475.788 -13544.094 -13544.094 264.35103 264.35103 31239.995 31239.995 3850.6849 3850.6849 2000 -13476 -13476 -13547.251 -13547.251 275.74874 275.74874 31289.903 31289.903 -671.74196 -671.74196 Loop time of 99.7345 on 1 procs for 1000 steps with 2000 atoms Performance: 0.866 ns/day, 27.704 hours/ns, 10.027 timesteps/s 99.9% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 99.252 | 99.252 | 99.252 | 0.0 | 99.52 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.074925 | 0.074925 | 0.074925 | 0.0 | 0.08 Output | 0.00018889 | 0.00018889 | 0.00018889 | 0.0 | 0.00 Modify | 0.34949 | 0.34949 | 0.34949 | 0.0 | 0.35 Other | | 0.05753 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.000000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 0, Press = 0 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 2000 -13476 -13476 -13547.251 -13547.251 275.74874 275.74874 31289.903 31289.903 -671.74196 -671.74196 3000 -13477.36 -13477.36 -13547.919 -13547.919 273.07065 273.07065 31299.662 31299.662 -1679.1472 -1679.1472 Loop time of 101.96 on 1 procs for 1000 steps with 2000 atoms Performance: 0.847 ns/day, 28.322 hours/ns, 9.808 timesteps/s 99.9% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 101.46 | 101.46 | 101.46 | 0.0 | 99.51 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.075243 | 0.075243 | 0.075243 | 0.0 | 0.07 Output | 0.00019188 | 0.00019188 | 0.00019188 | 0.0 | 0.00 Modify | 0.36184 | 0.36184 | 0.36184 | 0.0 | 0.35 Other | | 0.05893 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.000000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 0, Press = 0 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 3000 -13477.36 -13477.36 -13547.919 -13547.919 273.07065 273.07065 31299.662 31299.662 -1679.1472 -1679.1472 4000 -13475.435 -13475.435 -13545.869 -13545.869 272.58535 272.58535 31279.974 31279.974 317.35734 317.35734 Loop time of 100.355 on 1 procs for 1000 steps with 2000 atoms Performance: 0.861 ns/day, 27.877 hours/ns, 9.965 timesteps/s 99.9% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 99.865 | 99.865 | 99.865 | 0.0 | 99.51 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.075243 | 0.075243 | 0.075243 | 0.0 | 0.07 Output | 0.00021863 | 0.00021863 | 0.00021863 | 0.0 | 0.00 Modify | 0.35697 | 0.35697 | 0.35697 | 0.0 | 0.36 Other | | 0.05772 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.000000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 0, Press = 0 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 4000 -13475.435 -13475.435 -13545.869 -13545.869 272.58535 272.58535 31279.974 31279.974 317.35734 317.35734 5000 -13477.753 -13477.753 -13548.538 -13548.538 273.94224 273.94224 31305.724 31305.724 -2263.1888 -2263.1888 Loop time of 99.3349 on 1 procs for 1000 steps with 2000 atoms Performance: 0.870 ns/day, 27.593 hours/ns, 10.067 timesteps/s 99.9% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 98.856 | 98.856 | 98.856 | 0.0 | 99.52 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.074204 | 0.074204 | 0.074204 | 0.0 | 0.07 Output | 0.0001544 | 0.0001544 | 0.0001544 | 0.0 | 0.00 Modify | 0.348 | 0.348 | 0.348 | 0.0 | 0.35 Other | | 0.05695 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.000000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 274.478255108884, Press = 117.560378427542 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 5000 -13477.753 -13477.753 -13548.538 -13548.538 273.94224 273.94224 31305.724 31305.724 -2263.1888 -2263.1888 6000 -13477.119 -13477.119 -13541.742 -13541.742 250.10045 250.10045 31306.861 31306.861 -1550.2132 -1550.2132 Loop time of 101.526 on 1 procs for 1000 steps with 2000 atoms Performance: 0.851 ns/day, 28.202 hours/ns, 9.850 timesteps/s 99.9% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 101.01 | 101.01 | 101.01 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.075974 | 0.075974 | 0.075974 | 0.0 | 0.07 Output | 0.0001554 | 0.0001554 | 0.0001554 | 0.0 | 0.00 Modify | 0.37743 | 0.37743 | 0.37743 | 0.0 | 0.37 Other | | 0.05866 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.000000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 272.038670207192, Press = 8.32993432790081 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 6000 -13477.119 -13477.119 -13541.742 -13541.742 250.10045 250.10045 31306.861 31306.861 -1550.2132 -1550.2132 7000 -13474.206 -13474.206 -13543.972 -13543.972 269.99899 269.99899 31285.281 31285.281 -203.44537 -203.44537 Loop time of 101.291 on 1 procs for 1000 steps with 2000 atoms Performance: 0.853 ns/day, 28.136 hours/ns, 9.873 timesteps/s 100.0% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 100.78 | 100.78 | 100.78 | 0.0 | 99.49 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.07563 | 0.07563 | 0.07563 | 0.0 | 0.07 Output | 0.00020703 | 0.00020703 | 0.00020703 | 0.0 | 0.00 Modify | 0.37875 | 0.37875 | 0.37875 | 0.0 | 0.37 Other | | 0.05828 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.000000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 272.868706816387, Press = 22.478189692007 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 7000 -13474.206 -13474.206 -13543.972 -13543.972 269.99899 269.99899 31285.281 31285.281 -203.44537 -203.44537 8000 -13477.871 -13477.871 -13550.178 -13550.178 279.83793 279.83793 31275.738 31275.738 170.58125 170.58125 Loop time of 100.933 on 1 procs for 1000 steps with 2000 atoms Performance: 0.856 ns/day, 28.037 hours/ns, 9.908 timesteps/s 100.0% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 100.42 | 100.42 | 100.42 | 0.0 | 99.49 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.075623 | 0.075623 | 0.075623 | 0.0 | 0.07 Output | 0.00015648 | 0.00015648 | 0.00015648 | 0.0 | 0.00 Modify | 0.3766 | 0.3766 | 0.3766 | 0.0 | 0.37 Other | | 0.05766 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.000000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 273.905710906177, Press = -0.396180893891652 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 8000 -13477.871 -13477.871 -13550.178 -13550.178 279.83793 279.83793 31275.738 31275.738 170.58125 170.58125 9000 -13475.567 -13475.567 -13550.981 -13550.981 291.86127 291.86127 31290.444 31290.444 -963.00823 -963.00823 Loop time of 100.755 on 1 procs for 1000 steps with 2000 atoms Performance: 0.858 ns/day, 27.988 hours/ns, 9.925 timesteps/s 99.9% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 100.25 | 100.25 | 100.25 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.075448 | 0.075448 | 0.075448 | 0.0 | 0.07 Output | 0.00015784 | 0.00015784 | 0.00015784 | 0.0 | 0.00 Modify | 0.37298 | 0.37298 | 0.37298 | 0.0 | 0.37 Other | | 0.05786 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.000000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 273.479679250448, Press = -1.8135001168359 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 9000 -13475.567 -13475.567 -13550.981 -13550.981 291.86127 291.86127 31290.444 31290.444 -963.00823 -963.00823 10000 -13479.986 -13479.986 -13551.617 -13551.617 277.21962 277.21962 31279.786 31279.786 -320.56781 -320.56781 Loop time of 102.114 on 1 procs for 1000 steps with 2000 atoms Performance: 0.846 ns/day, 28.365 hours/ns, 9.793 timesteps/s 99.9% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 101.6 | 101.6 | 101.6 | 0.0 | 99.49 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.076205 | 0.076205 | 0.076205 | 0.0 | 0.07 Output | 0.00016086 | 0.00016086 | 0.00016086 | 0.0 | 0.00 Modify | 0.38108 | 0.38108 | 0.38108 | 0.0 | 0.37 Other | | 0.05874 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.000000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 273.167452851149, Press = 0.239377493700772 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 10000 -13479.986 -13479.986 -13551.617 -13551.617 277.21962 277.21962 31279.786 31279.786 -320.56781 -320.56781 11000 -13477.208 -13477.208 -13548.554 -13548.554 276.11314 276.11314 31275.908 31275.908 269.6638 269.6638 Loop time of 101.361 on 1 procs for 1000 steps with 2000 atoms Performance: 0.852 ns/day, 28.156 hours/ns, 9.866 timesteps/s 99.9% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 100.85 | 100.85 | 100.85 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.075252 | 0.075252 | 0.075252 | 0.0 | 0.07 Output | 0.00015588 | 0.00015588 | 0.00015588 | 0.0 | 0.00 Modify | 0.37653 | 0.37653 | 0.37653 | 0.0 | 0.37 Other | | 0.05786 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.000000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 273.116311427098, Press = -0.281899848688962 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 11000 -13477.208 -13477.208 -13548.554 -13548.554 276.11314 276.11314 31275.908 31275.908 269.6638 269.6638 12000 -13472.292 -13472.292 -13545.997 -13545.997 285.24724 285.24724 31280.644 31280.644 357.63238 357.63238 Loop time of 100.785 on 1 procs for 1000 steps with 2000 atoms Performance: 0.857 ns/day, 27.996 hours/ns, 9.922 timesteps/s 99.9% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 100.28 | 100.28 | 100.28 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.075127 | 0.075127 | 0.075127 | 0.0 | 0.07 Output | 0.00021283 | 0.00021283 | 0.00021283 | 0.0 | 0.00 Modify | 0.37306 | 0.37306 | 0.37306 | 0.0 | 0.37 Other | | 0.05752 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.000000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 273.254275553123, Press = -1.79830047380356 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 12000 -13472.292 -13472.292 -13545.997 -13545.997 285.24724 285.24724 31280.644 31280.644 357.63238 357.63238 13000 -13477.004 -13477.004 -13552.676 -13552.676 292.85613 292.85613 31302.521 31302.521 -2251.0244 -2251.0244 Loop time of 100.757 on 1 procs for 1000 steps with 2000 atoms Performance: 0.858 ns/day, 27.988 hours/ns, 9.925 timesteps/s 99.9% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 100.25 | 100.25 | 100.25 | 0.0 | 99.49 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.075929 | 0.075929 | 0.075929 | 0.0 | 0.08 Output | 0.00019928 | 0.00019928 | 0.00019928 | 0.0 | 0.00 Modify | 0.37479 | 0.37479 | 0.37479 | 0.0 | 0.37 Other | | 0.05792 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.000000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 273.433644005451, Press = -2.98119175630314 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 13000 -13477.004 -13477.004 -13552.676 -13552.676 292.85613 292.85613 31302.521 31302.521 -2251.0244 -2251.0244 14000 -13475.931 -13475.931 -13545.676 -13545.676 269.92147 269.92147 31280.031 31280.031 286.38294 286.38294 Loop time of 102.98 on 1 procs for 1000 steps with 2000 atoms Performance: 0.839 ns/day, 28.605 hours/ns, 9.711 timesteps/s 99.9% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 102.46 | 102.46 | 102.46 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.076647 | 0.076647 | 0.076647 | 0.0 | 0.07 Output | 0.00016045 | 0.00016045 | 0.00016045 | 0.0 | 0.00 Modify | 0.383 | 0.383 | 0.383 | 0.0 | 0.37 Other | | 0.05865 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.000000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 273.572873174435, Press = 0.288112004005871 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 14000 -13475.931 -13475.931 -13545.676 -13545.676 269.92147 269.92147 31280.031 31280.031 286.38294 286.38294 15000 -13477.526 -13477.526 -13542.544 -13542.544 251.62459 251.62459 31233.425 31233.425 4362.8912 4362.8912 Loop time of 101.649 on 1 procs for 1000 steps with 2000 atoms Performance: 0.850 ns/day, 28.236 hours/ns, 9.838 timesteps/s 99.9% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 101.14 | 101.14 | 101.14 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.07551 | 0.07551 | 0.07551 | 0.0 | 0.07 Output | 0.00041207 | 0.00041207 | 0.00041207 | 0.0 | 0.00 Modify | 0.37914 | 0.37914 | 0.37914 | 0.0 | 0.37 Other | | 0.05789 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.000000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 273.599367120018, Press = -1.03171724308701 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 15000 -13477.526 -13477.526 -13542.544 -13542.544 251.62459 251.62459 31233.425 31233.425 4362.8912 4362.8912 16000 -13478.251 -13478.251 -13547.374 -13547.374 267.51249 267.51249 31283.978 31283.978 -234.57546 -234.57546 Loop time of 99.1723 on 1 procs for 1000 steps with 2000 atoms Performance: 0.871 ns/day, 27.548 hours/ns, 10.083 timesteps/s 99.9% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 98.672 | 98.672 | 98.672 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.073875 | 0.073875 | 0.073875 | 0.0 | 0.07 Output | 0.00015792 | 0.00015792 | 0.00015792 | 0.0 | 0.00 Modify | 0.36836 | 0.36836 | 0.36836 | 0.0 | 0.37 Other | | 0.05769 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.000000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 273.607022729405, Press = 1.4691393560608 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 16000 -13478.251 -13478.251 -13547.374 -13547.374 267.51249 267.51249 31283.978 31283.978 -234.57546 -234.57546 17000 -13479.552 -13479.552 -13547.409 -13547.409 262.61383 262.61383 31289.56 31289.56 -717.64245 -717.64245 Loop time of 100.873 on 1 procs for 1000 steps with 2000 atoms Performance: 0.857 ns/day, 28.020 hours/ns, 9.913 timesteps/s 99.9% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 100.37 | 100.37 | 100.37 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.074791 | 0.074791 | 0.074791 | 0.0 | 0.07 Output | 0.00016232 | 0.00016232 | 0.00016232 | 0.0 | 0.00 Modify | 0.37434 | 0.37434 | 0.37434 | 0.0 | 0.37 Other | | 0.05761 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.000000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 273.56786015704, Press = -0.0161634893899602 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 17000 -13479.552 -13479.552 -13547.409 -13547.409 262.61383 262.61383 31289.56 31289.56 -717.64245 -717.64245 18000 -13475.526 -13475.526 -13545.372 -13545.372 270.31074 270.31074 31240.478 31240.478 3786.3435 3786.3435 Loop time of 99.6454 on 1 procs for 1000 steps with 2000 atoms Performance: 0.867 ns/day, 27.679 hours/ns, 10.036 timesteps/s 99.9% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 99.144 | 99.144 | 99.144 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.075001 | 0.075001 | 0.075001 | 0.0 | 0.08 Output | 0.00015457 | 0.00015457 | 0.00015457 | 0.0 | 0.00 Modify | 0.36935 | 0.36935 | 0.36935 | 0.0 | 0.37 Other | | 0.057 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.000000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 273.489341303406, Press = 1.57847999244884 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 18000 -13475.526 -13475.526 -13545.372 -13545.372 270.31074 270.31074 31240.478 31240.478 3786.3435 3786.3435 19000 -13477.287 -13477.287 -13550.078 -13550.078 281.70905 281.70905 31246.914 31246.914 2802.3579 2802.3579 Loop time of 99.1888 on 1 procs for 1000 steps with 2000 atoms Performance: 0.871 ns/day, 27.552 hours/ns, 10.082 timesteps/s 99.9% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 98.688 | 98.688 | 98.688 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.07475 | 0.07475 | 0.07475 | 0.0 | 0.08 Output | 0.00015517 | 0.00015517 | 0.00015517 | 0.0 | 0.00 Modify | 0.36858 | 0.36858 | 0.36858 | 0.0 | 0.37 Other | | 0.05726 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.000000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 273.458699485778, Press = -1.03649692799097 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 19000 -13477.287 -13477.287 -13550.078 -13550.078 281.70905 281.70905 31246.914 31246.914 2802.3579 2802.3579 20000 -13476.949 -13476.949 -13545.813 -13545.813 266.50932 266.50932 31251.221 31251.221 2792.3617 2792.3617 Loop time of 100.284 on 1 procs for 1000 steps with 2000 atoms Performance: 0.862 ns/day, 27.857 hours/ns, 9.972 timesteps/s 99.9% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 99.78 | 99.78 | 99.78 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.075173 | 0.075173 | 0.075173 | 0.0 | 0.07 Output | 0.00016162 | 0.00016162 | 0.00016162 | 0.0 | 0.00 Modify | 0.3722 | 0.3722 | 0.3722 | 0.0 | 0.37 Other | | 0.05712 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.000000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 273.463179368787, Press = 1.83220507979132 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 20000 -13476.949 -13476.949 -13545.813 -13545.813 266.50932 266.50932 31251.221 31251.221 2792.3617 2792.3617 21000 -13475.68 -13475.68 -13543.127 -13543.127 261.0239 261.0239 31268.39 31268.39 1467.8901 1467.8901 Loop time of 101.753 on 1 procs for 1000 steps with 2000 atoms Performance: 0.849 ns/day, 28.265 hours/ns, 9.828 timesteps/s 99.9% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 101.24 | 101.24 | 101.24 | 0.0 | 99.49 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.075414 | 0.075414 | 0.075414 | 0.0 | 0.07 Output | 0.00015562 | 0.00015562 | 0.00015562 | 0.0 | 0.00 Modify | 0.38171 | 0.38171 | 0.38171 | 0.0 | 0.38 Other | | 0.05835 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.000000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 273.494866920745, Press = 2.59518784277694 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 21000 -13475.68 -13475.68 -13543.127 -13543.127 261.0239 261.0239 31268.39 31268.39 1467.8901 1467.8901 22000 -13477.923 -13477.923 -13548.51 -13548.51 273.17971 273.17971 31277.34 31277.34 214.24018 214.24018 Loop time of 101.229 on 1 procs for 1000 steps with 2000 atoms Performance: 0.854 ns/day, 28.119 hours/ns, 9.879 timesteps/s 99.9% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 100.72 | 100.72 | 100.72 | 0.0 | 99.49 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.075609 | 0.075609 | 0.075609 | 0.0 | 0.07 Output | 0.00015637 | 0.00015637 | 0.00015637 | 0.0 | 0.00 Modify | 0.37814 | 0.37814 | 0.37814 | 0.0 | 0.37 Other | | 0.05823 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.000000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 273.42042620148, Press = 3.06010568855333 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 22000 -13477.923 -13477.923 -13548.51 -13548.51 273.17971 273.17971 31277.34 31277.34 214.24018 214.24018 23000 -13476.203 -13476.203 -13547.256 -13547.256 274.98504 274.98504 31263.286 31263.286 1596.5842 1596.5842 Loop time of 100.828 on 1 procs for 1000 steps with 2000 atoms Performance: 0.857 ns/day, 28.008 hours/ns, 9.918 timesteps/s 99.9% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 100.32 | 100.32 | 100.32 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.07498 | 0.07498 | 0.07498 | 0.0 | 0.07 Output | 0.00015698 | 0.00015698 | 0.00015698 | 0.0 | 0.00 Modify | 0.37422 | 0.37422 | 0.37422 | 0.0 | 0.37 Other | | 0.05769 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.000000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 273.387585468858, Press = 3.07032283104547 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 23000 -13476.203 -13476.203 -13547.256 -13547.256 274.98504 274.98504 31263.286 31263.286 1596.5842 1596.5842 24000 -13477.89 -13477.89 -13548.165 -13548.165 271.97287 271.97287 31267.538 31267.538 1098.8632 1098.8632 Loop time of 99.7342 on 1 procs for 1000 steps with 2000 atoms Performance: 0.866 ns/day, 27.704 hours/ns, 10.027 timesteps/s 99.9% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 99.233 | 99.233 | 99.233 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.074708 | 0.074708 | 0.074708 | 0.0 | 0.07 Output | 0.00019208 | 0.00019208 | 0.00019208 | 0.0 | 0.00 Modify | 0.36846 | 0.36846 | 0.36846 | 0.0 | 0.37 Other | | 0.05746 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.000000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 273.298816093865, Press = 2.25136532565099 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 24000 -13477.89 -13477.89 -13548.165 -13548.165 271.97287 271.97287 31267.538 31267.538 1098.8632 1098.8632 25000 -13476.638 -13476.638 -13548.696 -13548.696 278.86947 278.86947 31249.506 31249.506 2647.5576 2647.5576 Loop time of 100.785 on 1 procs for 1000 steps with 2000 atoms Performance: 0.857 ns/day, 27.996 hours/ns, 9.922 timesteps/s 99.9% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 100.28 | 100.28 | 100.28 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.075343 | 0.075343 | 0.075343 | 0.0 | 0.07 Output | 0.00019562 | 0.00019562 | 0.00019562 | 0.0 | 0.00 Modify | 0.37254 | 0.37254 | 0.37254 | 0.0 | 0.37 Other | | 0.05815 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 115998.0 ave 115998 max 115998 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 115998 Ave neighs/atom = 57.999000 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 = 273.265974989678, Press = 5.03865178795797 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 25000 -13476.638 -13476.638 -13548.696 -13548.696 278.86947 278.86947 31249.506 31249.506 2647.5576 2647.5576 26000 -13477.924 -13477.924 -13545.774 -13545.774 262.58928 262.58928 31255.662 31255.662 2288.4124 2288.4124 Loop time of 99.9314 on 1 procs for 1000 steps with 2000 atoms Performance: 0.865 ns/day, 27.759 hours/ns, 10.007 timesteps/s 100.0% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 99.43 | 99.43 | 99.43 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.074337 | 0.074337 | 0.074337 | 0.0 | 0.07 Output | 0.00015413 | 0.00015413 | 0.00015413 | 0.0 | 0.00 Modify | 0.3688 | 0.3688 | 0.3688 | 0.0 | 0.37 Other | | 0.0577 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.000000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 273.118252802619, Press = 4.88788611834392 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 26000 -13477.924 -13477.924 -13545.774 -13545.774 262.58928 262.58928 31255.662 31255.662 2288.4124 2288.4124 27000 -13476.644 -13476.644 -13546.516 -13546.516 270.4109 270.4109 31282.362 31282.362 71.082538 71.082538 Loop time of 99.3034 on 1 procs for 1000 steps with 2000 atoms Performance: 0.870 ns/day, 27.584 hours/ns, 10.070 timesteps/s 99.9% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 98.806 | 98.806 | 98.806 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.074384 | 0.074384 | 0.074384 | 0.0 | 0.07 Output | 0.0001536 | 0.0001536 | 0.0001536 | 0.0 | 0.00 Modify | 0.36653 | 0.36653 | 0.36653 | 0.0 | 0.37 Other | | 0.05668 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.000000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 273.099582995863, Press = 2.86265973911928 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 27000 -13476.644 -13476.644 -13546.516 -13546.516 270.4109 270.4109 31282.362 31282.362 71.082538 71.082538 28000 -13477.739 -13477.739 -13546.155 -13546.155 264.77611 264.77611 31290.851 31290.851 -583.06984 -583.06984 Loop time of 98.2391 on 1 procs for 1000 steps with 2000 atoms Performance: 0.879 ns/day, 27.289 hours/ns, 10.179 timesteps/s 99.9% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 97.747 | 97.747 | 97.747 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.073029 | 0.073029 | 0.073029 | 0.0 | 0.07 Output | 0.00015442 | 0.00015442 | 0.00015442 | 0.0 | 0.00 Modify | 0.36179 | 0.36179 | 0.36179 | 0.0 | 0.37 Other | | 0.05704 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.000000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 273.110775141061, Press = 2.06595140386486 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 28000 -13477.739 -13477.739 -13546.155 -13546.155 264.77611 264.77611 31290.851 31290.851 -583.06984 -583.06984 29000 -13474.451 -13474.451 -13544.417 -13544.417 270.77691 270.77691 31306.799 31306.799 -1755.0594 -1755.0594 Loop time of 102.006 on 1 procs for 1000 steps with 2000 atoms Performance: 0.847 ns/day, 28.335 hours/ns, 9.803 timesteps/s 99.9% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 101.49 | 101.49 | 101.49 | 0.0 | 99.49 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.076449 | 0.076449 | 0.076449 | 0.0 | 0.07 Output | 0.00019453 | 0.00019453 | 0.00019453 | 0.0 | 0.00 Modify | 0.38211 | 0.38211 | 0.38211 | 0.0 | 0.37 Other | | 0.0588 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.000000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 273.127288686473, Press = 2.04697313494176 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 29000 -13474.451 -13474.451 -13544.417 -13544.417 270.77691 270.77691 31306.799 31306.799 -1755.0594 -1755.0594 30000 -13476.776 -13476.776 -13547.73 -13547.73 274.60227 274.60227 31319.976 31319.976 -3163.4859 -3163.4859 Loop time of 101.476 on 1 procs for 1000 steps with 2000 atoms Performance: 0.851 ns/day, 28.188 hours/ns, 9.855 timesteps/s 99.9% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 100.96 | 100.96 | 100.96 | 0.0 | 99.49 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.076631 | 0.076631 | 0.076631 | 0.0 | 0.08 Output | 0.00019454 | 0.00019454 | 0.00019454 | 0.0 | 0.00 Modify | 0.38085 | 0.38085 | 0.38085 | 0.0 | 0.38 Other | | 0.05901 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.000000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 273.060346972729, Press = 1.68649502543947 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 30000 -13476.776 -13476.776 -13547.73 -13547.73 274.60227 274.60227 31319.976 31319.976 -3163.4859 -3163.4859 31000 -13477.064 -13477.064 -13545.292 -13545.292 264.05028 264.05028 31292.752 31292.752 -774.67876 -774.67876 Loop time of 100.4 on 1 procs for 1000 steps with 2000 atoms Performance: 0.861 ns/day, 27.889 hours/ns, 9.960 timesteps/s 99.9% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 99.892 | 99.892 | 99.892 | 0.0 | 99.49 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.075311 | 0.075311 | 0.075311 | 0.0 | 0.08 Output | 0.00018992 | 0.00018992 | 0.00018992 | 0.0 | 0.00 Modify | 0.37446 | 0.37446 | 0.37446 | 0.0 | 0.37 Other | | 0.0582 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.000000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 273.028967413956, Press = 0.661237569602157 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 31000 -13477.064 -13477.064 -13545.292 -13545.292 264.05028 264.05028 31292.752 31292.752 -774.67876 -774.67876 32000 -13476.582 -13476.582 -13547.497 -13547.497 274.45035 274.45035 31285.985 31285.985 -384.91058 -384.91058 Loop time of 100.995 on 1 procs for 1000 steps with 2000 atoms Performance: 0.855 ns/day, 28.054 hours/ns, 9.902 timesteps/s 100.0% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 100.49 | 100.49 | 100.49 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.075802 | 0.075802 | 0.075802 | 0.0 | 0.08 Output | 0.00015787 | 0.00015787 | 0.00015787 | 0.0 | 0.00 Modify | 0.37469 | 0.37469 | 0.37469 | 0.0 | 0.37 Other | | 0.05842 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.000000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 273.021155834029, Press = 0.981148694889127 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 32000 -13476.582 -13476.582 -13547.497 -13547.497 274.45035 274.45035 31285.985 31285.985 -384.91058 -384.91058 33000 -13472.361 -13472.361 -13546.353 -13546.353 286.35685 286.35685 31294.828 31294.828 -853.16119 -853.16119 Loop time of 99.4703 on 1 procs for 1000 steps with 2000 atoms Performance: 0.869 ns/day, 27.631 hours/ns, 10.053 timesteps/s 99.9% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 98.967 | 98.967 | 98.967 | 0.0 | 99.49 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.074692 | 0.074692 | 0.074692 | 0.0 | 0.08 Output | 0.00015251 | 0.00015251 | 0.00015251 | 0.0 | 0.00 Modify | 0.37075 | 0.37075 | 0.37075 | 0.0 | 0.37 Other | | 0.05791 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.000000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" print "flag: Temp = ${T}, Press = ${P}" flag: Temp = 273.125849438458, Press = 0.584468630844962 next a jump SELF top variable a loop 2000 run 1000 Neighbor list info ... update every 1 steps, delay 10 steps, check yes max neighbors/atom: 2000, page size: 100000 master list distance cutoff = 5.8 ghost atom cutoff = 5.8 binsize = 2.9, bins = 11 11 11 1 neighbor lists, perpetual/occasional/extra = 1 0 0 (1) pair kim, perpetual attributes: full, newton off, cut 5.8 pair build: full/bin/atomonly stencil: full/bin/3d bin: standard Per MPI rank memory allocation (min/avg/max) = 3.701 | 3.701 | 3.701 Mbytes Step TotEng v_etotal_metal PotEng v_pe_metal Temp v_T_metal Volume v_V_metal Press v_P_metal 33000 -13472.361 -13472.361 -13546.353 -13546.353 286.35685 286.35685 31294.828 31294.828 -853.16119 -853.16119 34000 -13477.298 -13477.298 -13547.573 -13547.573 271.97172 271.97172 31251.575 31251.575 2510.936 2510.936 Loop time of 99.5637 on 1 procs for 1000 steps with 2000 atoms Performance: 0.868 ns/day, 27.657 hours/ns, 10.044 timesteps/s 100.0% CPU use with 1 MPI tasks x 1 OpenMP threads MPI task timing breakdown: Section | min time | avg time | max time |%varavg| %total --------------------------------------------------------------- Pair | 99.066 | 99.066 | 99.066 | 0.0 | 99.50 Neigh | 0 | 0 | 0 | 0.0 | 0.00 Comm | 0.074353 | 0.074353 | 0.074353 | 0.0 | 0.07 Output | 0.00015354 | 0.00015354 | 0.00015354 | 0.0 | 0.00 Modify | 0.36604 | 0.36604 | 0.36604 | 0.0 | 0.37 Other | | 0.05729 | | | 0.06 Nlocal: 2000.00 ave 2000 max 2000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Nghost: 2941.00 ave 2941 max 2941 min Histogram: 1 0 0 0 0 0 0 0 0 0 Neighs: 0.00000 ave 0 max 0 min Histogram: 1 0 0 0 0 0 0 0 0 0 FullNghs: 116000.0 ave 116000 max 116000 min Histogram: 1 0 0 0 0 0 0 0 0 0 Total # of neighbors = 116000 Ave neighs/atom = 58.000000 Neighbor list builds = 0 Dangerous builds = 0 if "${V_metal}>${V0_metal_times1000}" then "jump SELF unstable" if "${T}>${T_low} && ${T}<${T_up} && ${P}>${P_low} && ${P}<${P_up}" then "jump SELF break" 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_T273.15.out" else "print 'not_converged' file output/vol_T273.15.out" print '${V}' file output/vol_T273.15.out 31282.0287503745 print "LAMMPS calculation completed" LAMMPS calculation completed quit 0