The long-range non-additive three-body dispersion interactions for the rare gases, alkali, and alkaline-earth atoms

A neural network potential energy surface of the Li3 system and quantum dynamics studies for the 7Li + 6Li2 → 6Li7Li + 6Li reaction

A high-precision global potential energy surface (PES) of the Li3 system is constructed based on high-level ab initio calculations, and the root-mean-square error is 5.54 cm-1. The short-range of the PES is fitted by the fundamental invariant neural network (FI-NN) method, while the long-range uses a function with an accurate asymptotic potential energy form, and the two regions are connected by a switching function. Based on the new PES, the statistical quantum-mechanical (SQM) and the time-dependent wave packet (TDWP) methods are used to study the dynamics of 7Li + 6Li2 (v = 0, j = 0) → 6Li7Li + 6Li reactions in the low collision energy region (10-11 to 10-3 cm-1) and the high collision energy region (8 to 800 cm-1), respectively. In the high collision energy region, the calculation results of the SQM method and the TDWP method are inconsistent, indicating that the reaction dynamics does not follow the statistical behavior in the high collision energy region. In addition, we found that the Coriolis coupling effect plays an important role in this reaction. The symmetric forward-backward scattering in the total DCS indicates that the reaction follows the complex-forming reaction mechanism.

Three-body potential and third virial coefficients for helium including relativistic and nuclear-motion effects

The non-additive three-body interaction potential for helium was computed using the coupled-cluster theory and the full configuration interaction method. The obtained potential comprises an improved nonrelativistic Born-Oppenheimer energy and the leading relativistic and nuclear-motion corrections. The mean absolute uncertainty of our calculations due to the incompleteness of the orbital basis set was determined employing complete-basis-set extrapolation techniques and was found to be 1.2%. For three helium atoms forming an equilateral triangle with the side length of 5.6 bohr - a geometry close to the minimum of the total potential energy surface - our three-body potential amounts to -90.6 mK, with an estimated uncertainty of 0.5 mK. An analytic function, developed to accurately fit the computed three-body interaction energies, was chosen to correctly describe the asymptotic behavior of the three-body potential for trimer configurations corresponding to both the three-atomic and the atom-diatom fragmentation channels. For large triangles with sides r12, r23, and r31, the potential takes correctly into account all angular terms decaying as r-l12 r-m23 r-n21 with l + m + n ≤ 14 for the nonrelativistic Born-Oppenheimer energy and l + m + n ≤ 9 for the post-Born-Oppenheimer corrections. We also developed a short-range analytic function describing the local behavior of the total uncertainty of the computed three-body interaction energies. Using both fits we calculated the third pressure and acoustic virial coefficients for helium and their uncertainties for a wide range of temperatures. The results of these calculations were compared with available experimental data and with previous theoretical determinations. The estimated uncertainties of present calculations are 3-5 times smaller than those reported in the best previous works.

Variational perturbation theory for dynamic polarizabilities and dispersion coefficients

An efficient method based on the variational perturbation theory (VPT) is proposed to conveniently calculate the atomic real- and imaginary-frequency dynamic polarizabilities and the interatomic dispersion coefficients. The developed method holds the great advantage that only the system ground state wave function and corresponding radial mean values are needed. Verification of the VPT method on one- and two-electron atoms indicates that the present approximation shows good agreement with calculations based on the sophisticated sum-over-states method. We apply the VPT method to examine the approximate Z-scaling laws of polarizabilities and dispersion coefficients in the He isoelectronic sequence, and to investigate the plasma screening effect on these quantities for embedded atoms. Our calculation demonstrates very well that the VPT method is capable of producing reasonably accurate static and dynamic polarizabilities as well as two- and three-atom dispersion coefficients for plasma-embedded atoms in a wide range of screening parameters.

Density Functional Response Calculations of Dispersion Coefficients C6 and C9 of Closed- and Open-Shell Systems

Dipole polarizabilities and C₆ and C₉ dispersion coefficients are computed for closed- and open-shell atoms and molecules, using dynamic (time-dependent) density functional (TD-DFT) linear response theory as implemented in the response module of the NWChem quantum chemistry package. The response module is capable of accurate calculations of these properties, based on spin-restricted and spin-unrestricted formalisms. The calculated static polarizabilities and dispersion coefficients are compared to available experimental and other theoretical data. The behavior of the dynamic polarizability at imaginary frequencies is analyzed for differently sized closed- and open-shell systems. An interpolation method enforcing the monotonic decrease of the polarizability with increasing imaginary frequency is beneficial for the integration used to obtain C₆ and C₉. Scaling of the TD-DFT data by ratios of the static polarizability, which can be calculated with a variety of methods, including highly accurate theories, may be used as a leading-order correction.

Bulk viscosity of gaseous argon from molecular dynamics simulations

The bulk viscosity of dilute argon gas is calculated using molecular dynamics simulations in the temperature range 150-500 K and is found to be proportional to density squared in the investigated range of densities 0.001-1 kg m^{-3}. A comparison of the results obtained using Lennard-Jones and Tang-Toennies models of pair interaction potential reveals that the value of the bulk viscosity coefficient is sensitive to the choice of the pair interaction model. The inclusion of the Axilrod-Teller-Muto three-body interaction in the model does not noticeably affect the values of the bulk viscosity in dilute states, contrary to the previously investigated case of dense fluids.

Three-Body Dispersion Potentials Involving Electric Octupole Coupling

Non-pairwise additive three-body dispersion potentials dependent upon one or more electric octupole moments are evaluated using the theory of molecular quantum electrodynamics. To simplify the perturbation theory calculations, an effective two-photon interaction Hamiltonian operator is employed. This leads to only third-order theory being required to evaluate energy shifts instead of the usual sixth-order formula, and the summation over six time-ordered sequences of virtual photon creation and annihilation events. Specific energy shifts computed include DD-DD-DO, DD-DO-DO, DO-DO-DO, and DD-DO-OO terms, where D and O are electric dipole and octupole moments, respectively. The formulae obtained are applicable to an arbitrary arrangement of the three particles, and we present explicit results for the equilateral triangle and collinear configurations, which complements the recently published DD-DD-OO potential. In this last case it was found that the contribution from the octupole weight-1 term could be viewed as a higher-order correction to the triple-dipole dispersion potential DD-DD-DD. In a similar fashion the octupole moment is decomposed into its irreducible components of weights-1 and -3, enabling insight to be gained into the potentials obtained in this study. Dispersion interaction energies proportional to mixed dipole-octupole polarisabilities, for example, are found to depend only on the weight-1 octupole moment for isotropic species and are retarded. Additional approximations are necessary in the evaluation of wave vector integrals for these cases in order to yield energy shifts that are valid in the near-zone.

Hyperspherical coupled channel calculations of energy and structure of (4)He-(4)He-Li(+) and its isotopic combinations

The ground state vibrational energy and spatial features of (4)He-(4)He-Li(+) and its triatomic isotopic complexes are studied using the slow variable discretization (SVD) method in the hyperspherical coordinates for the zero total angular momentum. Our results show that the dominant structure of the system is an isosceles triangle with the shorter side associated with the two Li(+)-He distances using the sum-of-potential approximation. Corrections caused by the induced dipole-induced dipole interactions on the He atoms are also investigated. The effects are seen to be small and have a minor influence on the binding energy and the structure of present system. The results are also compared with the full ab initio calculations including all the three-body interactions and information of three-body corrections is obtained.

Evaluation of the grand-canonical partition function using expanded Wang-Landau simulations. IV. Performance of many-body force fields and tight-binding schemes for the fluid phases of silicon

We extend Expanded Wang-Landau (EWL) simulations beyond classical systems and develop the EWL method for systems modeled with a tight-binding Hamiltonian. We then apply the method to determine the partition function and thus all thermodynamic properties, including the Gibbs free energy and entropy, of the fluid phases of Si. We compare the results from quantum many-body (QMB) tight binding models, which explicitly calculate the overlap between the atomic orbitals of neighboring atoms, to those obtained with classical many-body (CMB) force fields, which allow to recover the tetrahedral organization in condensed phases of Si through, e.g., a repulsive 3-body term that favors the ideal tetrahedral angle. Along the vapor-liquid coexistence, between 3000 K and 6000 K, the densities for the two coexisting phases are found to vary significantly (by 5 orders of magnitude for the vapor and by up to 25% for the liquid) and to provide a stringent test of the models. Transitions from vapor to liquid are predicted to occur for chemical potentials that are 10%-15% higher for CMB models than for QMB models, and a ranking of the force fields is provided by comparing the predictions for the vapor pressure to the experimental data. QMB models also reveal the formation of a gap in the electronic density of states of the coexisting liquid at high temperatures. Subjecting Si to a nanoscopic confinement has a dramatic effect on the phase diagram with, e.g. at 6000 K, a decrease in liquid densities by about 50% for both CMB and QMB models and an increase in vapor densities between 90% (CMB) and 170% (QMB). The results presented here provide a full picture of the impact of the strategy (CMB or QMB) chosen to model many-body effects on the thermodynamic properties of the fluid phases of Si.

Density-functional approach to the three-body dispersion interaction based on the exchange dipole moment

We implement and compute the density functional nonadditive three-body dispersion interaction using a combination of Tang-Karplus formalism and the exchange-dipole moment model of Becke and Johnson. The computation of the C9 dispersion coefficients is done in a non-empirical fashion. The obtained C9 values of a series of noble atom triplets agree well with highly accurate values in the literature. We also calculate the C9 values for a series of benzene trimers and find a good agreement with high-level ab initio values reported recently in the literature. For the question of damping of the three-body dispersion at short distances, we propose two damping schemes and optimize them based on the benzene trimers data, and the fitted analytic potentials of He3 and Ar3 trimers fitted to the results of high-level wavefunction theories available from the literature. Both damping schemes respond well to the optimization of two parameters.