Anisotropy of long range interactions between linear molecules: H2-H2 and H2-He

R-8 Dispersion Interaction: Derivation and Application to the Effective Fragment Potential Method

The anisotropic and isotropic R-8 dispersion contributions (disp8) are derived and implemented within the framework of the effective fragment potential (EFP) method formulated with imaginary frequency-dependent Cartesian polarizability tensors distributed at the centroids of the localized molecular orbitals (LMOs). Two forms of damping functions, intermolecular overlap-based and Tang-Toennies, are extended for disp8. To obtain LMO polarizability tensors centered at LMO centroids, an origin-shifting transformation is derived and implemented for the dipole-octopole polarizability tensor and the quadrupole-quadrupole polarizability tensor. The analytic gradient is derived and implemented for the isotropic disp8 contribution. Relative to the previously implemented empirical EFP disp8 energy, the isotropic disp8 component of the interaction energy improves the overall agreement of the EFP dispersion energies with the symmetry-adapted perturbation theory (SAPT) benchmarks, reducing the mean absolute errors (MAEs) and mean absolute percentage errors for most of the databases examined in this work. While the anisotropic disp8 can further enhance the accuracy of the EFP dispersion energy and yield smaller MAEs, significantly overbound dispersion energies are predicted by the anisotropic disp8 when the maximum element in the intermolecular overlap matrix is greater than 0.1, possibly due to the breakdown of the approximations made in the EFP dispersion derivation at a short range. For potential energy scan databases, the newly developed EFP dispersion model with isotropic disp8 yields the overall correct curvature and good agreement with SAPT benchmarks around equilibrium and longer but overestimates the dispersion interactions at a short range. While the overlap-based dispersion-damping functions produce better MAEs than Tang-Toennies damping functions, further improvement is needed to better screen the large attractive dispersion energies at a short range (overlap >0.1).

Dependence of the multipole moments, static polarizabilities, and static hyperpolarizabilities of the hydrogen molecule on the H-H separation in the ground singlet state

In this work, we provide values for the quadrupole moment Θ, the hexadecapole moment Φ, the dipole polarizability α, the quadrupole polarizability C, the dipole-octopole polarizability E, the second dipole hyperpolarizability γ, and the dipole-dipole-quadrupole hyperpolarizability B for the hydrogen molecule in the ground singlet state, evaluated by finite-field configuration interaction singles and doubles (CISD) and coupled-cluster singles and doubles (CCSD) methods for 26 different H-H separations r, ranging from 0.567 a.u. to 10.0 a.u. Results obtained with various large correlation-consistent basis sets are compared at the vibrationally averaged bond length r₀ in the ground state. Results over the full range of r values are presented at the CISD/d-aug-cc-pV6Z level for all of the independent components of the property tensors. In general, our values agree well with previous ab initio results of high accuracy for the ranges of H-H distances that have been treated in common. To our knowledge, for H₂ in the ground state, our results are the first to be reported in the literature for Φ for r > 7.0 a.u., γ and B for r > 6.0 a.u., and C and E for any H-H separation outside a narrow range around the potential minimum. Quantum Monte Carlo values of Θ have been given previously for H-H distances out to 10.0 a.u., but the statistical error is relatively large for r > 7.0 a.u. At the larger r values in this work, α_xx and α_zz show the expected functional forms, to leading order in r^-1. As r increases further, Θ and Φ vanish, while α, γ, and the components of B converge to twice the isolated-atom values. Components of C and E diverge as r increases. Vibrationally averaged values of the properties are reported for all of the bound states (vibrational quantum numbers υ = 0-14) with rotational quantum numbers J = 0-3.

Controllable binding of polar molecules and metastability of one-dimensional gases with attractive dipole forces

We explore one-dimensional samples of ultracold polar molecules with attractive dipole-dipole interactions and show the existence of a repulsive barrier caused by a strong quadrupole interaction between molecules. This barrier can stabilize a gas of ultracold KRb molecules and even lead to long-range wells supporting bound states between the molecules. The properties of these wells can be controlled by external electric fields, allowing the formation of long polymerlike chains of KRb and studies of quantum phase transitions by varying the effective interaction between molecules. We discuss the generalization of those results to other systems.

Determination of shear viscosity of molecular nitrogen (N2): molecular dynamic hard rotor methodology and the results

Determination of shear viscosity of molecular nitrogen (N(2)) by molecular dynamics (MD) in the high density range needs explicit incorporation of the rotational motion and therefore precise knowledge of angular dependence of N(2)-N(2) intermolecular potential. Newly designed Couette flow nonequilibrium molecular dynamic (NEMD) simulation procedure employs corrugated moving boundary, coupling the moving walls to translational and rotational motion exactly. Low density data on nitrogen viscosity show good agreement between MD data and experiment, confirming the radial dependence of the potential derived from quantum mechanical (QM) high precision calculations (coupled-cluster singles-and-doubles with a perturbative triples corrections [CCSD(T)]). Additionally, the angular dependence of the potential is verified using shear viscosity data for high density region, obtained from newly developed molecular dynamics (MD) simulations. It was demonstrated that the corrugated wall flow simulations provide results that are independent of the details of wall potential, fulfilling a basic requirement for application of MD simulations. These results are in good agreement with the equilibrium molecular dynamics (EMD) viscosity, derived from the Green-Kubo formula. Derived analytical dependence of the shear viscosity on the density and temperature shows that the MD data are in good agreement with experiment. Thus, MD simulations indicate that the CCSD(T) potential angular form is sufficiently precise for determination of the viscosity in a wide range of temperature and pressure.

Molecular nitrogen-N2 properties: The intermolecular potential and the equation of state

Quantum mechanical (QM) high precision calculations were used to determine N(2)-N(2) intermolecular interaction potential. Using QM numerical data the anisotropic potential energy surface was obtained for all orientations of the pair of the nitrogen molecules in the rotation invariant form. The new N(2)-N(2) potential is in reasonably good agreement with the scaled potential obtained by van der Avoird et al. using the results of Hartree-Fock calculations [J. Chem. Phys. 84, 1629 (1986)]. The molecular dynamics (MD) of the N(2) molecules has been used to determine nitrogen equation of state. The classical motion of N(2) molecules was integrated in rigid rotor approximation, i.e., it accounted only translational and rotational degrees of freedom. Fincham [Mol. Simul. 11, 79 (1993)] algorithm was shown to be superior in terms of precision and energy stability to other algorithms, including Singer [Mol. Phys. 33, 1757 (1977)], fifth order predictor-corrector, or Runge-Kutta, and was therefore used in the MD modeling of the nitrogen pressure [S. Krukowski and P. Strak, J. Chem. Phys. 124, 134501 (2006)]. Nitrogen equation of state at pressures up to 30 GPa (300 kbars) and temperatures from the room temperature to 2000 K was obtained using MD simulation results. Results of MD simulations are in very good agreement (the error below 1%) with the experimental data on nitrogen equation of state at pressures below 1 GPa (10 kbars) for temperatures below 1800 K [R. T. Jacobsen et al., J. Phys. Chem. Ref. Data 15, 735 (1986)]. For higher temperatures, the deviation is slightly larger, about 2.5% which still is a very good agreement. The slightly larger difference may be attributed to the vibrational motion not accounted explicitly by rigid rotor approximation, which may be especially important at high temperatures. These results allow to obtain reliable equation of state of nitrogen for pressures up to 30 GPa (300 kbars), i.e., close to molecular nitrogen stability limit, determined by Nellis et al. [Phys. Rev. Lett. 53, 1661 (1984)].