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.

Interaction-induced electric (hyper)polarizability in the dihydrogen-neon pair: basis set and electron correlation effects

We investigated the interaction (hyper)polarizability of neon-dihydrogen pairs by performing high-level ab initio calculations with atom/molecule-specific, purpose-oriented Gaussian basis sets. We obtained interaction-induced electric properties at the SCF, MP2, and CCSD levels of theory. At the CCSD level, for the T-shaped configuration, around the respective potential minimum of 6.437 a₀, the interaction-induced mean first hyperpolarizability varies for 5 <  R/a₀ < 10 as[Formula: see text]Again, at the CCSD level, but for the L-shaped configuration around the respective potential minimum of 6.572 a₀, this property varies for 5 <  R/a₀ < 10 as[Formula: see text] Graphical Abstract Interaction-induced mean dipole polarizability ([Formula: see text]) for the T-shaped configuration of H₂-Ne calculated at the SCF, MP2, and CCSD levels of theory.

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).

Diverging effects of isotopic fractionation upon molecular diffusion of noble gases in water: mechanistic insights through ab initio molecular dynamics simulations

Atmospheric noble gases are routinely used as natural tracers to analyze gas transfer processes in aquatic systems. Their isotopic ratios can be employed to discriminate between different physical transport mechanisms by comparison to the unfractionated atmospheric isotope composition. In many applications of aquatic systems molecular diffusion was thought to cause a mass dependent fractionation of noble gases and their isotopes according to the square root ratio of their masses. However, recent experiments focusing on isotopic fractionation within a single element challenged this broadly accepted assumption. The determined fractionation factors of Ne, Ar, Kr and Xe isotopes revealed that only Ar follows the prediction of the so-called square root relation, whereas within the Ne, Kr and Xe elements no mass-dependence was found. The reason for this unexpected divergence of Ar is not yet understood. The aim of our computational exercise is to establish the molecular-resolved mechanisms behind molecular diffusion of noble gases in water. We make the hypothesis that weak intermolecular interactions are relevant for the dynamical properties of noble gases dissolved in water. Therefore, we used ab initio molecular dynamics to explicitly account for the electronic degrees of freedom. Depending on the size and polarizability of the hydrophobic particles such as noble gases, their motion in dense and polar liquids like water is subject to different diffusive regimes: the inter-cavity hopping mechanism of small particles (He, Ne) breaks down if a critical particle size achieved. For the case of large particles (Kr, Xe), the motion through the water solvent is governed by mass-independent viscous friction leading to hydrodynamical diffusion. Finally, Ar falls in between the two diffusive regimes, where particle dispersion is propagated at the molecular collision time scale of the surrounding water molecules.

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.