Coupled cluster calculations of the ground state potential and interaction induced electric properties of the mixed dimers of helium, neon and argon

A Simple Model for the Pauli Repulsion with Possible Utility in QM, MM and Chemical Education

The Pauli repulsion is the intermolecular force responsible for the volume and low compressibility of condensed-phase matter at normal conditions. A simple model for this force is presented, wherein per-atom electron densities are represented as spherical charge distributions that are prevented from significantly overlapping. In the example of two noble gas atoms approaching one another beyond their van der Waals radii, the distance between the centers of the electronic charge distributions becomes larger than the distance between the nuclei, giving rise to an unfavorable electrostatic interaction. For the purpose of calculating this interaction, the model is further simplified by representing the per-atom electron density as a negative point charge, loosely inspired by the classical Drude oscillator. The dispersion interaction is simplified to an R-6 term, centered on the aforementioned point charges. Despite the gross simplicity of the resulting formalism, near-quantitative agreement with high-level QM interaction energies of noble gas dimers is achieved. Accordingly, the present model is thought to have utility in force fields, in post-HF and post-DFT dispersion corrections, and in chemical education.

Impact of Combination Rules, Level of Theory, and Potential Function on the Modeling of Gas- and Condensed-Phase Properties of Noble Gases

The systems of noble gases are particularly instructive for molecular modeling due to the elemental nature of their interactions. They do not normally form bonds nor possess a (permanent) dipole moment, and the only forces determining their bonding/clustering stems from van der Waals forces─dispersion and Pauli repulsion, which can be modeled by empirical potential functions. Combination rules, that is, formulas to derive parameters for pair potentials of heterodimers from parameters of corresponding homodimers, have been studied at length for the Lennard-Jones 12-6 potentials but not in great detail for other, more accurate, potentials. In this work, we examine the usefulness of nine empirical potentials in their ability to reproduce quantum mechanical (QM) benchmark dissociation curves of noble gas dimers (He, Ne, Ar, Kr, and Xe homo- and heterodimers), and we systematically study the efficacy of different permutations of combination relations for each parameter of the potentials. Our QM benchmark comprises dissociation curves computed by several different coupled cluster implementations as well as symmetry-adapted perturbation theory. The two-parameter Lennard-Jones potentials were decisively outperformed by more elaborate potentials that sport a 25-30 times lower root-mean-square error (RMSE) when fitted to QM dissociation curves. Very good fits to the QM dissociation curves can be achieved with relatively inexpensive four- or even three-parameter potentials, for instance, the damped 14-7 potential (Halgren, J. Am. Chem. Soc. 1992, 114, 7827-7843), a four-parameter Buckingham potential (Werhahn et al., Chem. Phys. Lett. 2015, 619, 133-138), or the three-parameter Morse potential (Morse, Phys. Rev. 1929, 34, 57-64). Potentials for heterodimers that are generated from combination rules have an RMSE that is up to 20 times higher than potentials that are directly fitted to the QM dissociation curves. This means that the RMSE, in particular, for light atoms, is comparable in magnitude to the well-depth of the potential. Based on a systematic permutation of combination rules, we present one or more combination rules for each potential tested that yield a relatively low RMSE. Two new combination rules are introduced that perform well, one for the van der Waals radius σij as ( 1 2 ( σ i 3 + σ j 3 ) ) 1 / 3 and one for the well-depth ϵij as ( 1 2 ( ϵ i - 2 + ϵ j - 2 ) ) - 1 / 2 . The QM data and the fitted potentials were evaluated in the gas phase against experimental second virial coefficients for homo- and heterodimers, the latter of which allowed evaluation of the combination rules. The fitted models were used to perform condensed phase molecular dynamics simulations to verify the melting points, liquid densities at the melting point, and the enthalpies of vaporization produced by the models for pure substances. Subtle differences in the benchmark potentials, in particular, the well-depth, due to the level of theory used were found here to have a profound effect on the macroscopic properties of noble gases: second virial coefficients or the bulk properties in simulations. By explicitly including three-body dispersion in molecular simulations employing the best pair potential, we were able to obtain accurate melting points as well as satisfactory densities and enthalpies of vaporization.

Universal Pairwise Interatomic van der Waals Potentials Based on Quantum Drude Oscillators

Repulsive short-range and attractive long-range van der Waals (vdW) forces play an appreciable role in the behavior of extended molecular systems. When using empirical force fields, the most popular computational methods applied to such systems, vdW forces are typically described by Lennard-Jones-like potentials, which unfortunately have a limited predictive power. Here, we present a universal parameterization of a quantum-mechanical vdW potential, which requires only two free-atom properties─the static dipole polarizability α1 and the dipole-dipole C6 dispersion coefficient. This is achieved by deriving the functional form of the potential from the quantum Drude oscillator (QDO) model, employing scaling laws for the equilibrium distance and the binding energy, and applying the microscopic law of corresponding states. The vdW-QDO potential is shown to be accurate for vdW binding energy curves, as demonstrated by comparing to the ab initio binding curves of 21 noble-gas dimers. The functional form of the vdW-QDO potential has the correct asymptotic behavior at both zero and infinite distances. In addition, it is shown that the damped vdW-QDO potential can accurately describe vdW interactions in dimers consisting of group II elements. Finally, we demonstrate the applicability of the atom-in-molecule vdW-QDO model for predicting accurate dispersion energies for molecular systems. The present work makes an important step toward constructing universal vdW potentials, which could benefit (bio)molecular computational studies.

Basis sets dependency in constructing spectroscopy-accuracy Ab Initio global electric dipole moment functions

Recently, more attention have been paid on the construction of dipole moment functions (DMF) using theoretical methods. However, the computational methods to construct DMFs are not validated as much as those for potential energy surfaces do. In this letter, using Ar ⋯ He as an example, we tested how spectroscopy-accuracy DMFs can be constructed using ab initio methods. We especially focused on the basis set dependency in this scenario, i.e., the convergence of DMF with the sizes of basis sets, basis set superposition error, and mid-bond functions. We also tested the explicitly correlated method, which converges with smaller basis sets than the conventional methods do. This work can serve as a pictorial sample of all these computational technologies behaving in the context of constructing DMFs.

Conformal Analytical Potential for All the Rare Gas Dimers over the Full Range of Internuclear Distances

An analytical model for the potential between two rare gas atoms at distances between R=0 to R→∞ is assumed to be conformal with the previously published potential for He_{2} [J. Chem. Phys. 142, 131102 (2015)JCPSA60021-960610.1063/1.4916740]. The potential curves of the rare gas dimers all have the same shape and only depend on the well parameters D_{e} and R_{e}. The potentials and the vibrational levels for the 11 homonuclear and heteronuclear dimers for which recent ab initio calculations are available agree, within several percent, with the ab initio results. For the other rare gas dimers, the new potential provides the first realistic estimates for the potentials.

Can Density Functional Theory Be Trusted for High-Order Electric Properties? The Case of Hydrogen-Bonded Complexes

This work reports on an extensive assessment of the performance of a wide palette of density functional approximations in predicting the (high-order) electric properties of hydrogen-bonded complexes. To this end, we compute the electronic and vibrational contributions to the electric polarizability and the first and second hyperpolarizabilities, using the CCSD(T)/aug-cc-pVTZ level of theory as reference. For all the studied properties, the average absolute errors below 20% can only be obtained using the CAM-B3LYP functional, while LC-BLYP and MN15 are shown to be only slightly less accurate (average absolute errors not exceeding 30%). Among Minnesota density functionals, i.e., M06, M06-2X, and MN15, we only recommend the latter one, which quite accurately predicts the electronic and vibrational (hyper)polarizabilities. We also analyze the optimal tuning of the range-separation parameter μ for the LC-BLYP functional, finding that this approach does not bring any systematic improvement in the predictions of electronic and vibrational (hyper)polarizabilities and the accuracy of computed properties is largely system-dependent. Finally, we report huge errors in predicting the vibrational second hyperpolarizability by ωB97X, M06, and M06-2X functionals. Based on the explicit evaluation of anharmonic terms contributing to the second hyperpolarizability, this failure is traced down to a poor determination of third- and fourth-order energy derivatives with respect to normal modes. These results reveal serious flaws of some density functional approximations and suggest caution in selecting the appropriate functional to calculate not only electronic and vibrational (hyper)polarizabilities but also other molecular properties that contain vibrational anharmonic contributions.

Thermodynamic Properties of Vapor-Liquid Equilibria from Monte-Carlo Simulation using ab initio Intermolecular Potentials of Systems H2-H2 and F2-F2

In this work, we have been carried out GEMC-NVT simulations in the temperature range 18 K–32 K for fluid hydrogen and in range 60 K–140 K for fluid fluorine using four our developed ab initio 5-site intermolecular potentials for dimers H2-H2 and F2-F2, respectively. The thermodynamic properties of vapor-liquid equilibria and the critical points of fluids hydrogen and fluorine were calculated with the obtained densities of coexisting phases and vapor pressures. The simulation results drived from ab initio pair potentials were compared with those from ab initio potential plus three-body Axilrod-Teller potential and experimental data as well as those from Monte Carlo simulation using Lennard-Jones potentials, Deiters equation of state (D1-EOS) and Benedict-Webb-Rubin equation of state (EOS) reported in the literature.

Partitioning of interaction-induced nonlinear optical properties of molecular complexes. I. Hydrogen-bonded systems

Interaction-induced electronic and vibrational (hyper)polarizabilities were decomposed into different interaction types (electrostatic, exchange, induction and dispersion).

On the physical origins of interaction-induced vibrational (hyper)polarizabilities

This paper presents the results of a pioneering exploration of the physical origins of vibrational contributions to the interaction-induced electric properties of molecular complexes. In order to analyze the excess nuclear relaxation (hyper)polarizabilities, a new scheme was proposed which relies on the computationally efficient Bishop-Hasan-Kirtman method for determining the nuclear relaxation contributions to electric properties. The extension presented herein is general and can be used with any interaction-energy partitioning method. As an example, in this study we employed the variational-perturbational interaction-energy decomposition scheme (at the MP2/aug-cc-pVQZ level) and the extended transition state method by employing three exchange-correlation functionals (BLYP, LC-BLYP, and LC-BLYP-dDsC) to study the excess properties of the HCN dimer. It was observed that the first-order electrostatic contribution to the excess nuclear relaxation polarizability cancels with the negative exchange repulsion term out to a large extent, resulting in a positive value of Δα(nr) due to the contributions from the delocalization and the dispersion terms. In the case of the excess nuclear relaxation first hyperpolarizability, the pattern of interaction contributions is very similar to that for Δα(nr), both in terms of their sign as well as relative magnitude. Finally, our results show that the LC-BLYP and LC-BLYP-dDsC functionals, which yield smaller values of the orbital relaxation term than BLYP, are more successful in predicting excess properties.

Heavy rare-gas atomic pairs and the "double penalty" issue: Isotropic Raman lineshapes by Kr2, Xe2, and KrXe at room temperature

We report absolutely calibrated isotropic Raman lineshapes for Kr2 and Xe2 and for KrXe at 294.5 K and compare them to quantum-mechanically generated lineshapes by using state-of-the-art second-order Møller-Plesset and DFT/B3LYP data sets for the induced mean dipole polarizability ᾱ. A very good agreement between the numerical and the experimental data was observed but the large uncertainty margins and the short Raman frequency interval probed in our experiment prevented us from rating on a more refined scale the performance of the tested ᾱ models. These drawbacks are inherent in isotropic Raman spectrum measurements and amplified for dissimilar pairs because, for such systems and spectra, the unreliable operation of subtracting optical signals of comparable magnitude occurs twice per Raman frequency shift value, thus penalizing twice the quality of the measured data. In light of our findings and of previously reported evidence about related electric properties in Kr2 and Xe2 and in KrXe, we are left with no doubt as to the consistency of the induced-polarizability and interatomic-potential data used for these three systems at the reported level of accuracy.

Transport coefficients of helium-argon mixture based on ab initio potential

The viscosity, thermal conductivity, diffusion coefficient, and thermal diffusion factor of helium-argon mixtures are calculated for a wide range of temperature and for various mole fractions up to the 12th order of the Sonine polynomial expansion with an ab initio intermolecular potential. The calculated values for these transport coefficients are compared with other data available in the open literature. The comparison shows that the obtained transport coefficients of helium-argon mixture have the best accuracy for the moment.

Ab initiostudy of interaction-induced NMR shielding constants in mixed rare gas dimers

The interaction-induced contribution to the NMR shielding constants in homonuclear A2 and heteronuclear AB (A,B=He,Ne,Ar) dimers is obtained ab initio by employing a coupled cluster singles and doubles with perturbative treatment of triples wave function model and extended correlation-consistent basis sets. The second virial coefficients entering the expansion of the property with the density are then computed in a fully quantum mechanical approach, for temperatures ranging from the limit of dissociation of the dimer to well above standard conditions. The results can be used to describe the density and temperature dependence of the shielding constants in binary mixtures of helium, neon, and argon. The predicted effects should be observable for the interaction of 21Ne with other rare gases.

Hyper-Rayleigh spectral intensities of gaseous Kr–Xe mixture

Binary collision-induced hyper-Rayleigh (CIHR) spectra of Kr-Xe gaseous system are computed quantum mechanically and classically within the frequency range up to 380 cm(-1). The intensities are expressed in absolute units. The details of the theory developed for the CIHR spectra are given and the properties of the profiles as well as the depolarization ratio frequency dependence are discussed. The contributions to the spectra related to the vector b10(r) and the septor b30(r) components of the hyperpolarizability tensor are evaluated.