The water dimer cation, (H2O)2+, has long served as a prototypical reference system for water oxidation chemistry. In spite of this status, a definitive explanation for the anomalous—and dominant—features in the experimental vibrational spectrum [Gardenier, G. H.; Johnson, M. A.; McCoy, A. B. J. Phys. Chem. A, 2009, 113, 4772–4779] has not been determined, and harmonic analyses qualitatively fail to reproduce these features. In this computational study, accurate quantum chemistry methods are combined with a fully coupled, six-dimensional anharmonic model to show that the unassigned bands are the result of resonant mode interactions and strong anharmonic coupling. Such coupling is fundamentally due to the unique electronic structure of this open-shell ion and the manner in which auxiliary modes affect the natural charge-transfer properties of the shared-proton stretch. These unique vibrational signatures provide a key reference point for modern spectroscopic and mechanistic analyses of water-oxidation catalysts.

This work investigates the use of multiple-timestep schemes in imaginary time for computationally efficient ab initio equilibrium path integral simulations of quantum molecular motion. In the simplest formulation, only every nth path integral replica is computed at the target level of electronic structure theory, whereas the remaining low-level replicas still account for nuclear motion quantum effects with a more computationally economical theory. Motivated by recent developments for multiple-timestep techniques in real-time classical molecular dynamics, both 1-electron (atomic-orbital basis set) and 2-electron (electron correlation) truncations are shown to be effective. Structural distributions and thermodynamic averages are tested for representative analytic potentials and ab initio molecular examples. Target quantum chemistry methods include density functional theory and second-order Møller–Plesset perturbation theory, although any level of theory is formally amenable to this framework. For a standard two-level splitting, computational speedups of 1.6–4.0x are observed when using a 4-fold reduction in time slices; an 8-fold reduction is feasible in some cases. Multitiered options further reduce computational requirements and suggest that quantum mechanical motion could potentially be obtained at a cost not significantly different from the cost of classical simulations.

The isomers of a hydrated Cu(I) ion with n = 1–10 water molecules were investigated by using ab initio quantum chemistry and an automated isomer-search algorithm. The electronic structure and vibrational spectra of the hundreds of resulting isomers were used to analyze the source of the observed bonding patterns. A structural evolution from dominantly two-coordinate structures (n = 1–4) toward a mixture of two- and three-coordinate structures was observed at n = 5–6, where the stability provided by expanded hydrogen-bonding was competitive with the dominantly electrostatic interaction between the water ligand and remaining binding sites of the metal ion. Further hydration (n = 7–10) led to a mixture of three- and four-coordinate structures. The metal ion was found, through spectroscopic signatures, to appreciably perturb the O–H bonds of even third-shell water molecules, which highlighted the ability of this nominally simple ion to partially activate the surrounding water network.

The structures, properties, and spectroscopic signatures of oxidized water clusters,(H2O)+n=6–21, are examined in this work, to provide fundamental insight into renewable energy and radiological processes. Computational quantum chemistry approaches are employed to sample cluster morphologies, yielding hundreds of low-lying isomers with low barriers to interconversion. The ion–radical pair-separation trend, however, which was observed in previous computational studies and in small-cluster spectroscopy experiments, is shown to continue in this larger cluster size regime. The source of this trend is preferential solvation of the hydronium ion by water, including effects beyond the first solvation shell. The fundamental conclusion of this work, therefore, is that the initially formed ion–radical dimer, which has served as a prototypical model of oxidized water, is a nascent species in large, oxidized water clusters and, very likely, bulk water.

Because of both experimental and computational challenges, protonated tryptophan has remained the last aromatic amino acid for which the intrinsic structures of low-energy conformers have not been unambiguously solved. The IR–IR–UV hole-burning spectroscopy technique has been applied to overcome the limitations of the commonly used IR–UV double resonance technique and to measure conformer-specific vibrational spectra of TrpH+, cooled to T = 10 K. Anharmonic ab initio vibrational spectroscopy simulations unambiguously assign the dominant conformers to the two lowest-energy geometries from benchmark coupled-cluster structure computations. The match between experimental and ab initio spectra provides an unbiased validation of the calculated structures of the two experimentally observed conformers of this benchmark ion. Furthermore, the vibrational spectra provide conformer-specific signatures of the stabilizing interactions, including hydrogen bonding and an intramolecular cation-π interaction.

A multiple-timestep ab initio molecular dynamics scheme based on varying the two-electron integral screening method used in Hartree–Fock or density functional theory calculations is presented. Although screening is motivated by numerical considerations, it is also related to separations in the length- and timescales characterizing forces in a molecular system: Loose thresholds are sufficient to describe fast motions over short distances, while tight thresholds may be employed for larger length scales and longer times, leading to a practical acceleration of ab initio molecular dynamics simulations. Standard screening approaches can lead, however, to significant discontinuities in (and inconsistencies between) the energy and gradient when the screening threshold is loose, making them inappropriate for use in dynamics. To remedy this problem, a consistent window-screening method that smooths these discontinuities is devised. Further algorithmic improvements reuse electronic-structure information within the dynamics step and enhance efficiency relative to a naı̈ve multiple-timestepping protocol. The resulting scheme is shown to realize meaningful reductions in the cost of Hartree–Fock and B3LYP simulations of a moderately large system, the protonated sarcosine/glycine dipeptide embedded in a 19-water cluster.

This work describes an approach to accelerate ab initio Born–Oppenheimer molecular dynamics (MD) simulations by exploiting the inherent timescale separation between contributions from different atom-centered Gaussian basis sets. Several MD steps are propagated with a cost-efficient, low-level basis set, after which a dynamical correction accounts for large basis set relaxation effects in a time-reversible fashion. This multiple-timestep scheme is shown to generate valid MD trajectories, on the basis of rigorous testing for water clusters, the methanol dimer, an alanine polypeptide, protonated hydrazine, and the oxidized water dimer. This new approach generates observables that are consistent with those of target basis set trajectories, including MD-based vibrational spectra. This protocol is shown to be valid for Hartree–Fock, density functional theory, and second-order Møller–Plesset perturbation theory approaches. Recommended pairings include 6-31G as a low-level basis set for 6-31G** or 6-311G**, as well as cc-pVDZ as the subset for accurate dynamics with aug-cc-pVTZ. Demonstrated cost savings include factors of 2.6–7.3 on the systems tested and are expected to remain valid across system sizes.

Ionized water clusters serve as a model of water-splitting chemistry for energetic purposes, as well as postradiolytic events in condensed-phase systems. Structures, properties, and relative energies are presented for oxidized water clusters, (H2O)n=1–5+, using equation-of-motion coupled-cluster theory approaches. In small clusters, an ion–radical contact pair OH···H3O+ is known to form upon ionization. The transition from n = 4 to n = 5 molecules in the cluster, however, is found to demarcate a size regime in which a propensity for the ion and radical to separate exists. This trend is consistent with recent experimental vibrational analyses. Decomposition of the cluster energetics reveals that preferential solvation of the hydronium cation by water serves as the dominant driving force for this pair separation, which should persist in larger clusters and bulk water ionization.

The dynamic, quantum structure of [Mg(H2)n=1–6]2+complexes is investigated via ab initio path integral molecular dynamics simulations. These complexes represent the strong, σ-complex regime of metal–H2 interactions and are representative of bonding motifs found in metal–organic frameworks. Significant nuclear motion within the coordination sphere is observed, even though the ligands remain largely intact. Quantum effects are found to be important in the H–H and metal–H2 stretch coordinates, but the remaining motion in the molecule is well represented by classical simulations. Nearly free rotation of the dihydrogen moiety is observed in all complexes. Statistical averages and distributions of structural parameters are found to deviate nontrivially from the same parameters in static, equilibrium structures.

The cation dimer of water and hydrogen sulfide, [(H2O)(H2S)]+, serves as a fundamental model for the oxidation chemistry of H2S. The known oxidative metabolism of H2S by biological species in sulfur-rich environments has motivated the study of the inherent properties of this benchmark complex, with possible mechanistic implications for modern water oxidation chemistry. The low-energy isomer of this open-shell ion is a proton-transferred (PT) structure, H3O+⋯SH˙. An alternative PT structure, H3S+⋯OH˙, and a hemibonded (HB) isomer, [H2O·SH2]+, are also stable isomers, placing this complex intermediate to known (H2O)2+ (PT) and (H2S)2+ (HB) limiting regimes. This intermediate character suggested the possibility of unique molecular motion, even in the vibrational ground state. Path integral molecular dynamics and anharmonic vibrational spectroscopy simulations have been performed in this study, in order to understand the inherent quantum molecular motion of this complex. The resulting structural distributions were found to deviate significantly from both classical and harmonic analyses, including the observation of large-amplitude anharmonic motion of the central proton and nearly free rotation of the terminal hydrogens. The predicted vibrational spectra are particularly unique and suggest characteristic signatures of the strong electronic interactions and anharmonic vibrational mode couplings in this radical cation.