The calculated infrared spectrum of Cl-H2O using a new full dimensional ab initio potential surface and dipole moment surface

Efficient vibrationally correlated calculations using n-mode expansion-based kinetic energy operators

Due to its efficiency and flexibility, the n-mode expansion is a frequently used tool for representing molecular potential energy surfaces in quantum chemical simulations. In this work, we investigate the performance of n-mode expansion-based models of kinetic energy operators in general polyspherical coordinate systems. In particular, we assess the operators with respect to accuracy in vibrationally correlated calculations and their effect on potential energy surface construction with the adaptive density guided approach. Our results show that the n-mode expansion-based operator variants are reliable and systematically improvable approximations of the full kinetic energy operator. Moreover, we introduce a workflow to generate the n-mode expanded kinetic energy operators on-the-fly within the adaptive density guided approach. This scheme can be applied in studies of species and coordinate systems, for which an analytical form of the kinetic energy operator is not available.

Quantum molecular simulations of micro-hydrated halogen anions

We report the results of a detailed and accurate investigation focused on structures and energetics of poly-hydrated halides employing first-principles polarizable halide-water potentials to describe the underlying forces. Following a bottom-up data-driven potential approach, we initially looked into the classical behavior of higher-order X^-(H₂O)_N clusters. We have located several low-lying energies, such as global and local minima, structures for each cluster, with various water molecules (up to N = 8) surrounding the halide anion (X^- = F^-, Cl^-, Br^-, I^-), employing an evolutionary programming method. It is found that the F^--water clusters exhibit different structural configurations than the heavier halides, however independently of the halide anion, all clusters show in general a selective growth with the anion preferring to be connected to the outer shell of the water molecule arrangements. In turn, path-integral molecular dynamics simulations are performed to incorporate explicitly nuclear quantum and thermal effects in describing the nature of halide ion microsolvation in such prototypical model systems. Our data reveal that at low finite temperatures, nuclear quantum effects affect certain structural properties, such as weakening hydrogen bonding between the halide anion and water molecules, with minor distortions in the water network beyond the first hydration shell, indicating local structure rearrangements. Such structural characteristics and the promising cluster size trends observed in the single-ion solvation energies motivated us to draw connections of small size cluster data to the limits of continuum bulk values, toward the investigation of the challenging computational modeling of bulk single ion hydration.

A Benchmark Protocol for DFT Approaches and Data-Driven Models for Halide-Water Clusters

Dissolved ions in aqueous media are ubiquitous in many physicochemical processes, with a direct impact on research fields, such as chemistry, climate, biology, and industry. Ions play a crucial role in the structure of the surrounding network of water molecules as they can either weaken or strengthen it. Gaining a thorough understanding of the underlying forces from small clusters to bulk solutions is still challenging, which motivates further investigations. Through a systematic analysis of the interaction energies obtained from high-level electronic structure methodologies, we assessed various dispersion-corrected density functional approaches, as well as ab initio-based data-driven potential models for halide ion-water clusters. We introduced an active learning scheme to automate the generation of optimally weighted datasets, required for the development of efficient bottom-up anion-water models. Using an evolutionary programming procedure, we determined optimized and reference configurations for such polarizable and first-principles-based representation of the potentials, and we analyzed their structural characteristics and energetics in comparison with estimates from DF-MP2 and DFT+D quantum chemistry computations. Moreover, we presented new benchmark datasets, considering both equilibrium and non-equilibrium configurations of higher-order species with an increasing number of water molecules up to 54 for each F, Cl, Br, and I anions, and we proposed a validation protocol to cross-check methods and approaches. In this way, we aim to improve the predictive ability of future molecular computer simulations for determining the ongoing conflicting distribution of different ions in aqueous environments, as well as the transition from nanoscale clusters to macroscopic condensed phases.

Permutationally Invariant Polynomial Expansions with Unrestricted Complexity

A general strategy is presented for constructing and validating permutationally invariant polynomial (PIP) expansions for chemical systems of any stoichiometry. Demonstrations are made for three categories of gas-phase dynamics and kinetics: collisional energy-transfer trajectories for predicting pressure-dependent kinetics, three-body collisions for describing transient van der Waals adducts relevant to atmospheric chemistry, and nonthermal reactivity via quasiclassical trajectories. In total, 30 systems are considered with up to 15 atoms and 39 degrees of freedom. Permutational invariance is enforced in PIP expansions with as many as 13 million terms and 13 permutationally distinct atom types by taking advantage of petascale computational resources. The quality of the PIP expansions is demonstrated through the systematic convergence of in-sample and out-of-sample errors with respect to both the number of training data and the order of the expansion, and these errors are shown to predict errors in the dynamics for both reactive and nonreactive applications. The parallelized code distributed as part of this work enables the automation of PIP generation for complex systems with multiple channels and flexible user-defined symmetry constraints and for automatically removing unphysical unconnected terms from the basis set expansions, all of which are required for simulating complex reactive systems.

Discrete Oligomers and Polymers of Chloride Monohydrate Can Form in Encapsulated Environments: Structures and Infrared Spectra of [Cl4 (H2 O)4 ]4- and {[Cl(H2 O)]- }∞

A discrete tetrachloride tetrahydrate cluster, [Cl₄ (H₂ O)₄ ]^4- , was obtained with a partially-fluorinated triaminocyclopropenium cation, [C₃ (N(CH₂ CF₃ )₂ )(NEt₂ )(NPr₂ )]⁺ . The cluster consists of a [Cl₂ (H₂ O)₂ ]^2- square with each Cl^- coordinated by another H₂ O bridged to another Cl^- . A 1D polymer of chloride monohydrate, {[Cl(H₂ O)]^- }_∞ , was obtained with [C₃ (N(CH₂ CF₃ )₂ )₂ (NBuMe)]⁺ . The tetrameric and polymeric structures were found to be computationally-unstable in the gas phase which indicates that an encapsulated environment is essential for their isolation. DFT and DFTB calculations were carried out on gas-phase [Cl₄ (H₂ O)₄ ]^4- to assist the infrared assignments. Anharmonically-corrected B3LYP transition frequencies were in close agreement with experiment, but DFTB models were only appropriate for qualitative interpretation. Solid-state DFTB calculations allowed the vibrational modes to be assigned. The results found are consistent with "discrete" chloride hydrates.

Electrostatics, Charge Transfer, and the Nature of the Halide-Water Hydrogen Bond

Binary halide-water complexes X^-(H₂O) are examined by means of symmetry-adapted perturbation theory, using charge-constrained promolecular reference densities to extract a meaningful charge-transfer component from the induction energy. As is known, the X^-(H₂O) potential energy surface (for X = F, Cl, Br, or I) is characterized by symmetric left and right hydrogen bonds separated by a C_2v-symmetric saddle point, with a tunneling barrier height that is <2 kcal/mol except in the case of F^-(H₂O). Our analysis demonstrates that the charge-transfer energy is correspondingly small (<2 kcal/mol except for X = F), considerably smaller than the electrostatic interaction energy. Nevertheless, charge transfer plays a crucial role determining the conformational preferences of X^-(H₂O) and provides a driving force for the formation of quasi-linear X··· H-O hydrogen bonds. Charge-transfer energies correlate well with measured O-H vibrational redshifts for the halide-water complexes and also for OH^-(H₂O) and NO₂^-(H₂O), providing some indication of a general mechanism.

A Series of Discrete Dichloride Dihydrates: Characterisation and Symmetry Effects

A series of three discrete dichloride dihydrates [Cl₂ (H₂ O)₂ ]^2- have been isolated with different triaminocyclopropenium (TAC) cations and with different crystallographic symmetries. The cluster exhibits D_2h symmetry with the tris(dimethylamino)cyclopropenium cation [C₃ (NMe₂ )₃ ]⁺ , C_2h symmetry with the fluorinated cation [C₃ (N(CH₂ CF₃ )₂ )(NBu₂ )₂ ]⁺ (containing two 2,2,2-trifluoroethyl substituents) and C_2v symmetry with the more fluorinated [C₃ (N(CH₂ CF₃ )₂ )₂ (NBu₂ )]⁺ cation. The effect of symmetry on the infrared spectra of the dichloride ion-pair clusters, as well as deuterated analogues, has been investigated. The D_2h - and C_2h -symmetric clusters each exhibit two stretching bands in the infrared at 3427 and 3368 cm^-1 for D_2h symmetry and 3444 and 3392 cm^-1 for C_2h symmetry, whereas the C_2v -symmetric cluster exhibits three bands at 3475, 3426 and 3373 cm^-1 . Computational studies were carried out on a [Cl₂ (H₂ O)₂ ]^2- cluster with C_2v symmetry to aid the infrared band assignments.

One water to bind a chloride-chloride ion pair: isolation of discrete [Cl2(H2O)]2- in the solid state

A discrete dichloride ion pair in the form of a monohydrate, [Cl₂(H₂O)]^2-, was isolated using the triaminocyclopropenium cation [C₃(NHex₂)(N(CH₂CF₃)₂)₂]⁺. Although this ion pair is calculated to be unstable in the gas phase, the ionic lattice and weak CH-Cl hydrogen bonds assist the stabilization of the cluster. The D₂O and HDO isotopomers were also prepared and characterized.

A Discrete Dichloride Tetrahydrate Trapped by a Cyclopropenium Cation: Structure and Spectroscopic Properties

A discrete dichloride tetrahydrate cluster, [Cl₂ (H₂ O)₄ ]^2- , was obtained as a salt of the bis(diphenylamino)diethylamino cyclopropenium cation [C₃ (NPh₂ )₂ (NEt₂ )]⁺ and characterized by single-crystal X-ray diffraction and infrared spectroscopy. This chloride-chloride ion-pair cluster consists of a [Cl₂ (H₂ O)₂ ]^2- square with opposite edges bridged by water molecules to give a chair-like structure of the non-hydrogen atoms. The solid-state structure is essentially the same as the calculated gas-phase structure. Infrared spectra were also collected on the deuterium analogue [Cl₂ (D₂ O)₄ ]^2- . Computational studies were carried out on gas-phase [Cl₂ (H₂ O)₄ ]^2- to confirm the infrared band assignments in the solid state. The structure and infrared spectrum are consistent with the discrete nature of the cluster.

Spectroscopic Signatures of Mode-Dependent Tunnel Splitting in the Iodide-Water Binary Complex

The gas-phase vibrational spectrum of the isolated iodide-water cluster ion (I^-·H₂O), first reported in 1996, presents one of the most difficult, long-standing spectroscopic puzzles involving ion microhydration. Although the spectra of the smaller halides are well described in the context of an asymmetrical ground-state structure in which only one OH group is hydrogen-bonded to the ion, the I^-·H₂O spectrum displays multiplet structures with partially resolved rotational patterns that are additionally influenced by quantum nuclear spin statistics. In this study, this complex behavior is unraveled with a combination of experimental methods, including ion preparation in a temperature-controlled ion trap and spectral simplification through applications of tag-free, two-color IR-IR double-resonance spectroscopy. Analysis of the double-resonance spectra reveals a vibrational ground-state tunneling splitting of about 20 cm^-1, which is on the same order as the spacing between the peaks that comprise the multiplet structure. These findings are further supported by the results obtained from a fully coupled, six-dimensional calculation of the vibrational spectrum. The underlying level structure can then be understood as a consequence of experimentally measurable, vibrational mode-dependent tunneling splittings (which, in the case of the ground vibrational state, is comparable to the rotational energy spacing between levels with K_a = 0 and 1), as well as Fermi resonance interactions. The latter include the hydrogen-bonded OH stretches and combination bands that involve the HOH bend overtones and soft-mode excitations of frustrated translation and rotation displacements of the water molecule relative to the ion. These anharmonic couplings yield closely spaced bands that are activated in the IR by borrowing intensity from the OH stretch fundamentals.

Improving the description of solvent pairwise interactions using local solute/solvent three-body functions. The case of halides and carboxylates in aqueous environment

We propose a general strategy to remediate force-field artifacts in describing pairwise interactions among similar molecules M in the vicinity of another chemical species, C, like water molecules interacting at short distance from a monoatomic ion. This strategy is based on introducing a three-body potential energy term that alters the pairwise interactions among M-type molecules when they lie at short range from the species C. In other words the species C is the center of a space domain where the pairwise interactions among the molecules M is altered. Here, we apply it to improve the description of the water interactions provided by the polarizable water model TCPE/2013 in the vicinity of halides, from F^- to At^- , and of the prototypical carboxylate anion CH₃ COO^- . We show the accuracy and the transferability of such an approach to investigate not only the hydration process of single anions but also of a salt solution NH 4 + / Cl - in aqueous phase. This strategy can be used to remediate the drawbacks of any kind of force fields. © 2019 Wiley Periodicals, Inc.

Calculation of vibrational spectroscopic and geometrical characteristics of the [F(HF)2]- and [F(DF)2]- complexes using the second-order vibrational perturbation theory and a 6D variational method

Vibrational spectroscopic and average geometrical parameters of the strong H-bonded complexes [F(HF)₂]^- and [F(DF)₂]^- are determined for the first time from nine-dimensional (9D) perturbative and 6D variational calculations. The frequencies and intensities for all fundamental and some combination and overtone transitions obtained by the method of second-order vibrational perturbation theory (VPT2) are reported. A two-fold decrease in the H-F (D-F) stretching band frequency and a more than ten-fold increase in the intensity of this band upon complexation are predicted. The theoretical frequencies for both isolated isotopologues are in satisfactory agreement (to better than 70 cm^-1) with the scarce experimental data obtained in condensed phases. The main purpose of variational calculations is to analyze the intermode anharmonic coupling and the changes in the geometrical parameters upon vibrational excitation and H/D isotopic substitution. The equilibrium nuclear configuration and the 2D potential energy surface (PES) of [F(HF)₂]^- for H-F stretches are calculated in the MP2/6-311++G(3df,3pd), CCSD(T)/6-311++G(3df,3pd), CCSD(T)/aug-cc-pVTZ, and CCSD(T)/d-aug-cc-pVTZ approximations with the basis set superposition error taken into account. Anharmonic vibrational problems are solved by the variational method for 2D, 4D, and 6D systems of H-bond and H-F (D-F) stretches and in-plane bends. The VPT2 calculations and calculations of the PESs for 4D and 6D systems are performed in the MP2/6-311++G(3df,3pd) approximation. Comparison of variational anharmonic solutions for different vibrational subsystems demonstrates the influence of intermode anharmonic coupling on the mixing of wave functions and spectroscopic and geometrical characteristics. The inverse Ubbelohde effect is predicted and substantiated.

Quantum dynamics of ClH2O- photodetachment: Isotope effect and impact of anion vibrational excitation

Photodetachment of the ClH₂O^- anion is investigated using full-dimensional quantum mechanics on accurate potential energy surfaces of both the anion and neutral species. Detailed analysis of the photoelectron spectrum and the corresponding wavefunctions reveals that the photodetachment leads to, in the product channel of the exothermic HCl + OH → Cl + H₂O reaction, the formation of numerous Feshbach resonances due apparently to slow energy transfer from H₂O vibrational modes to the dissociation coordinate. These long-lived resonances can be grouped into two broad peaks in the low-resolution photoelectron spectrum, which is in good agreement with available experiments, and they are assigned to the ground and first excited OH stretching vibrational manifolds of H₂O complexed with Cl. In addition, effects of isotope substitution on the photoelectron spectrum were small. Finally, photodetachment of the vibrationally excited ClH₂O^- in the ionic hydrogen bond mode is found to lead to Feshbach resonances with higher stretching vibrational excitations in H₂O.

Effects of vibrational excitation on the F + H2O → HF + OH reaction: dissociative photodetachment of overtone-excited [F-H-OH]. Chem Sci 2017;8:7821-7833. [PMID: 29163919 PMCID: PMC5674243 DOI: 10.1039/c7sc03364h]

Photodetaching vibrationally excited FH₂O^– channels energy into the reaction coordinate of the F + H₂O reaction, as shown in this joint experimental-theoretical study.

The reaction F + H₂O → HF + OH is a four-atom system that provides an important benchmark for reaction dynamics. Hydrogen atom transfer at the transition state for this reaction is expected to exhibit a strong dependence on reactant vibrational excitation. In the present study, the vibrational effects are examined by photodetachment of vibrationally excited F^–(H₂O) precursor anions using photoelectron-photofragment coincidence (PPC) spectroscopy and compared with full six-dimensional quantum dynamical calculations on ab initio potential energy surfaces. Prior to photodetachment at hν_UV = 4.80 eV, the overtone of the ionic hydrogen bond mode in the precursor F^–(H₂O), 2ν_IHB at 2885 cm^–1, was excited using a tunable IR laser. Experiment and theory show that vibrational energy in the anion can be effectively carried away by the photoelectron upon a Franck–Condon photodetachment, and also show evidence for an increase of branching into the F + H₂O reactant channel. The experimental results suggest a greater role for product rotational excitation than theory. Improved potential energy surfaces and longer wavepacket propagation times would be helpful to further examine the nature of the discrepancy.

Modeling Molecular Interactions in Water: From Pairwise to Many-Body Potential Energy Functions

Almost 50 years have passed from the first computer simulations of water, and a large number of molecular models have been proposed since then to elucidate the unique behavior of water across different phases. In this article, we review the recent progress in the development of analytical potential energy functions that aim at correctly representing many-body effects. Starting from the many-body expansion of the interaction energy, specific focus is on different classes of potential energy functions built upon a hierarchy of approximations and on their ability to accurately reproduce reference data obtained from state-of-the-art electronic structure calculations and experimental measurements. We show that most recent potential energy functions, which include explicit short-range representations of two-body and three-body effects along with a physically correct description of many-body effects at all distances, predict the properties of water from the gas to the condensed phase with unprecedented accuracy, thus opening the door to the long-sought "universal model" capable of describing the behavior of water under different conditions and in different environments.

Largely reduced grid densities in a vibrational self-consistent field treatment do not significantly impact the resultingwavenumbers

Especially for larger molecules relevant to life sciences, vibrational self-consistent field (VSCF) calculations can become unmanageably demanding even when only first and second order potential coupling terms are considered. This paper investigates to what extent the grid density of the VSCF’s underlying potential energy surface can be reduced without sacrificing accuracy of the resulting wavenumbers. Including single-mode and pair contributions, a reduction to eight points per mode did not introduce a significant deviation but improved the computational efficiency by a factor of four. A mean unsigned deviation of 1.3% from the experiment could be maintained for the fifteen molecules under investigation and the approach was found to be applicable to rigid, semi-rigid and soft vibrational problems likewise. Deprotonated phosphoserine, stabilized by two intramolecular hydrogen bonds, was investigated as an exemplary application.

The impact of highly correlated potential energy surfaces on the anharmonically corrected IR spectrum of acetonitrile

This paper discusses the quality and feasibility of highly correlated ab initio techniques in a vibrational self-consistent field (VSCF) approach using acetonitrile as a model system. The topical renormalized coupled-cluster technique exploiting the similarity-transformed Hamiltonian's left eigenstates (i.e. CR-CC(2,3)) is investigated alongside the well-known Hartree-Fock (HF), Møller-Plesset second-order perturbation theory (MP2) and coupled cluster (CCSD(T)) methods. The inclusion of mode triple interactions is discussed and it is found that the use of an effective core potential (ECP) serves as a viable compromise during the highly demanding task of computing such contributions, thus enabling a grid-based evaluation of three mode interaction terms with coupled cluster techniques also for larger molecules. In this context, a previously proposed reduced coupling scheme [1] is investigated, confirming the applicability of this technique to a system exhibiting a rather complex electronic structure. A combination of Ahlrichs' triple-ζ valence polarized (TZVP) basis set with Dunning's set of core-valence correlation functions is found to deliver results in good agreement with experiment while being computationally very feasible. Since CH3CN exhibits four degenerate vibrational degrees of freedom, it serves as an ideal model system for critically assessing the qualities of the degenerate second-order perturbation theory corrected (DPT2) VSCF technique. Besides fundamental vibrations, a thorough investigation of overtone transitions and combination bands is conducted by means of comparing the results to both available and newly recorded experimental data.

Calculations of rovibrational energies and dipole transition intensities for polyatomic molecules using MULTIMODE

We report rigorous calculations of rovibrational energies and dipole transition intensities for three molecules using a new version of the code MULTIMODE. The key features of this code which permit, for the first time, such calculations for moderately sized but otherwise general polyatomic molecules are briefly described. Calculations for the triatomic molecule BF(2) are done to validate the code. New calculations for H(2)CO and H(2)CS are reported; these make use of semiempirical potentials but ab initio dipole moment surfaces. The new dipole surface for H(2)CO is a full-dimensional fit to the dipole moment obtained with the coupled-cluster with single and double excitations and a perturbative treatment of triple excitations method with the augmented correlation consistent triple zeta basis set. Detailed comparisons are made with experimental results from a fit to relative data for H(2)CS and absolute intensities from the HITRAN database for H(2)CO.

A new glance at HCl-monohydrate spectroscopy, using on-the-fly dynamics

On-the-fly dynamics is used to analyze the remarkably anharmonic infrared spectroscopy of crystalline HCl monohydrate, an ionic solid composed of H3O+ and Cl-. The dominant intense infrared feature is shown to originate from specific sections of the hydronium trajectory, in which one of the H-atoms interacts strongly with a neighboring Cl-.