Search

The copper hydroxide ion, CuOH+, serves as the catalytic core in several recently developed water-splitting catalysts, and an understanding of its chemistry is critical to determining viable cataly...

Chemical activation of water by a single chlorine atom was examined computationally for clusters of chlorine radicals and water in a size regime just prior to internal hydration of water/ions, Cl·(H2O)n=4-8,17. This investigation follows a recent analysis of this radical-molecule interaction [Christensen et al. J. Phys. Chem. A 2019, 123, 8657] for n = 1-4, which demonstrated that n = 4 marked a transition in which an oxidized-water structural motif became viable, albeit high in energy. Thousands of unique isomers were computed in the present analysis, which resulted in three structural classes of isomers, including intact hydrated chlorine, hydrogen-transferred (HCl)(OH·)(H2O)n-1, and charge-transferred (Cl-)(H3O+)(OH·)(H2O)n-2 configurations. The electronic structures of these classes were investigated, along with harmonic vibrational signatures that probed the degree of water-network perturbations and generated experimentally verifiable computational predictions. The main outcome of this analysis is that the charge-transferred isomers were stabilized considerably upon increased hydration-leading to an energetic crossover with the hydrogen-transferred forms-but the degree of hydration was surprisingly still not sufficient to achieve crossover between the intact chlorine-water complexes and these charge-separated configurations. Internal hydration of the ions appears to be necessary in order to achieve this separation, which will likely occur at larger cluster sizes.

The gas-phase vibrational spectrum of the isolated iodide–water cluster ion (I–·H2O), first reported in 1996, presents one of the most difficult, long-standing spectroscopic puzzles involving ion m...

This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange–correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods. Methods highlighted in Q-Chem 5 include a suite of tools for modeling core-level spectroscopy, methods for describing metastable resonances, methods for computing vibronic spectra, the nuclear–electronic orbital method, and several different energy decomposition analysis techniques. High-performance capabilities including multithreaded parallelism and support for calculations on graphics processing units are described. Q-Chem boasts a community of well over 100 active academic developers, and the continuing evolution of the software is supported by an “open teamware” model and an increasingly modular design. This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange-correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods. Methods highlighted in Q-Chem 5 include a suite of tools for modeling core-level spectroscopy, methods for describing metastable resonances, methods for computing vibronic spectra, the nuclear-electronic orbital method, and several different energy decomposition analysis techniques. High-performance capabilities including multithreaded parallelism and support for calculations on graphics processing units are described. Q-Chem boasts a community of well over 100 active academic developers, and the continuing evolution of the software is supported by an "open teamware" model and an increasingly modular design.

The partial chemical activation of water by reactive radicals was examined computationally for small clusters of chlorine and water, Cl•(H2O)n=1–4. Using an automated isomer-search procedure, dozen...

The 1,3-diaza-2,4-diborobutane (NBNB) molecule serves as the smallest model complex of an intramolecular "dihydrogen bond," which involves a nominally hydrogen-bonding interaction between amine and borane hydrogen atoms. In the present study, the role of this dihydrogen bond in influencing the inherent molecular dynamics of NBNB is investigated computationally with ab initio molecular dynamics and path integral molecular dynamics techniques, as well as vibrational spectra analysis and static quantum chemistry computations. These simulations indicate that the dihydrogen-bonding interaction impacts both the high- and low-frequency motions of the molecule, with the dominant motions involving low-frequency backbone isomerization and terminal amine rotation. Geometric isotope effects were found to be modest. The analysis also addresses the paradoxical fostering of amine rotation via a relatively strong dihydrogen bond interaction. Electrostatic and geometric factors most directly explain this effect, and although some orbital evidence was found for a small covalent component of this interaction, the dynamics and electronic structure suggest that electrostatic contributions are the controlling factors for molecular motion in NBNB.

Characterizing structural and phase transformations of water at the molecular level is key to understanding a variety of multiphase processes ranging from ice nucleation in the atmosphere to hydration of biomolecules and wetting of solid surfaces. In this study, state-of-the-art quantum simulations with a many-body water potential energy surface, which exhibits chemical and spectroscopic accuracy, are carried out to monitor the microscopic melting of the water hexamer through the analysis of vibrational spectra and appropriate structural order parameters as a function of temperature. The water hexamer is specifically chosen as a case study due to the central role of this cluster in the molecular-level understanding of hydrogen bonding in water. Besides being in agreement with the experimental data available for selected isomers at very low temperature, the present results provide quantitative insights into the interplay between energetic, entropic, and nuclear quantum effects on the evolution of water clusters from "solid-like" to "liquid-like" structures. This study thus demonstrates that computer simulations can now bridge the gap between measurements currently possible for individual isomers at very low temperature and observations of isomer mixtures at ambient conditions.

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.

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.

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

In combined quantum mechanical/molecular mechanical (QM/MM) free energy calculations, it is often advantageous to have a frozen geometry for the quantum mechanical (QM) region. For such multiple-environment single-system (MESS) cases, two schemes are proposed here for estimating the polarization energy: the first scheme, termed MESS-E, involves a Roothaan step extrapolation of the self-consistent field (SCF) energy; whereas the other scheme, termed MESS-H, employs a Newton-Raphson correction using an approximate inverse electronic Hessian of the QM region (which is constructed only once). Both schemes are extremely efficient, because the expensive Fock updates and SCF iterations in standard QM/MM calculations are completely avoided at each configuration. They produce reasonably accurate QM/MM polarization energies: MESS-E can predict the polarization energy within 0.25 kcal/mol in terms of the mean signed error for two of our test cases, solvated methanol and solvated β-alanine, using the M06-2X or ωB97X-D functionals; MESS-H can reproduce the polarization energy within 0.2 kcal/mol for these two cases and for the oxyluciferin-luciferase complex, if the approximate inverse electronic Hessians are constructed with sufficient accuracy.

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

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.

Polynomial extrapolation of response (‘z-vector’) elements is shown to reduce the cost of correlated-wavefunction, classical ab initio molecular dynamics. Demonstrations with resolution-of-the-identity Moller–Plesset perturbation theory (RI-MP2) show that the number of response equation iterations is reduced by a factor of 2–3, thereby enabling accelerated MP2 dynamics. Coupled with previously demonstrated Fock matrix extrapolation, the combined iterative SCF and z-vector calculations are reduced to a minority share of the total calculation. Though demonstrated for MP2, these methods can also be generalized to higher levels of electron correlation or excited state methods.

It's the Network Numerous reactions of small molecules and ions in the atmosphere take place in the confines of watery aerosols. Relph et al. (p. 308 ; see the Perspective by Siefermann and Abel ) explored the specific influence of a water cluster's geometry on the transformation of solvated nitrosonium (NO + ) to nitrous acid (HONO). The reaction involves (O)N–O(H) bond formation with one water molecule, concomitant with proton transfer to additional, surrounding water molecules. Vibrational spectroscopy and theoretical simulations suggest that certain arrangements of the surrounding water network are much more effective than others in accommodating this charge transfer, and thus facilitating the reaction.

Ground state structures with side-on nitrosyl ({eta}{sup 2}-NO) and isonitrosyl (ON) ligands have been observed in a variety of transition-metal complexes. In contrast, excited state structures with bent-NO ligands have been proposed for years but never directly observed. Here we use picosecond time-resolved infrared spectroscopy and density functional theory (DFT) modeling to study the photochemistry of Co(CO){sub 3}(NO), a model transition-metal-NO compound. Surprisingly, we have observed no evidence for ON and {eta}{sup 2}-NO structural isomers, but have observed two bent-NO complexes. DFT modeling of the ground and excited state potentials indicates that the bent-NO complexes correspond to triplet excited states. Photolysis of Co(CO){sub 3}(NO) with a 400-nm pump pulse leads to population of a manifold of excited states which decay to form an excited state triplet bent-NO complex within 1 ps. This structure relaxes to the ground triplet state in ca. 350 ps to form a second bent-NO structure.

We report the geometries of two distorted T-shaped benzene dimer structures optimized at the RI-MP2/aug-cc-pVTZ level of theory. At the extrapolated RI-MP2/aug-cc-pV(TQ)Z level, the Cs over atom and Cs over bond configurations were found to be lower in energy than the conventionally accepted C2v T-shaped structure by 0.146 and 0.163 kcal/mol, respectively. When DCCSD(T)/6-311+G(2df,p) corrections were included, these structures remained lower in energy than the C2v reference by 0.127 and 0.132 kcal/mol, respectively, with Cs over bond as the minimum energy T-shaped structure. While not the focus of this Letter, we also report that the C2v T-shaped configuration is stabilized by 0.31 kcal/mol over the C2h parallel-displaced configuration at the DCCSD(T)/aug-cc-pVTZ approximation to the CCSD(T)/CBS limit. � 2007 Elsevier B.V. All rights reserved.

We present the local charge state modification at room temperature of small insulator-supported molecular ensembles formed by 1,1'-ferrocenedicarboxylic acid on calcite. Single electron tunnelling between the conducting tip of a noncontact atomic force microscope (NC-AFM) and the molecular islands is observed. By joining NC-AFM with Kelvin probe force microscopy, successive charge build-up in the sample is observed from consecutive experiments. Charge transfer within the islands and structural relaxation of the adsorbate/surface system is suggested by the experimental data.

Analytic gradients of dual-basis Hartree-Fock and density functional theory energies have been derived and implemented, which provide the opportunity for capturing large basis-set gradient effects at reduced cost. Suggested pairings for gradient calculations are 6-31G/6-31G**, dual[-f,-d]/cc-pVTZ, and 6-311G*/6-311 + +G(3df,3pd). Equilibrium geometries are produced within 0.0005 A of large-basis results for the latter two pairings. Though a single, iterative SCF response equation must be solved (unlike standard SCF gradients), it may be obtained in the smaller basis set, and integral screening further reduces the cost for well-chosen subsets. Total nuclear force calculations exhibit up to 75% savings, relative to large-basis calculations.

A “tiered” approach to Monte Carlo sampling of nuclear configurations is presented for ab initio, self-consistent field (SCF)-based potentials, including Hartree-Fock and density functional theory. Rather than Metropolis testing only the final SCF energy, individual cycle energies are tested in a tiered fashion, without approximation. Accordingly, rejected configurations are terminated early in the SCF procedure. The method is shown to properly obey detailed balance, and effective modifications are presented for cases in which the initial SCF guess is particularly poor. Demonstrations on simple systems are provided, including an assessment of the thermal properties of the neutral water dimer with B3LYP/6-31++G**. Cost analysis indicates a factor-of-two reduction in SCF cycles, which makes the method competitive with accelerated molecular dynamics sampling techniques, without the need for forces.

On-the-fly, ab initio classical molecular dynamics are demonstrated with an underlying dual basis set potential energy surface. Dual-basis self-consistent field (Hartree-Fock and density functional theory) and resolution-of-the-identity second-order Moller-Plesset perturbation theory (RI-MP2) dynamics are tested for small systems, including the water dimer. The resulting dynamics are shown to be faithful representations of their single-basis analogues for individual trajectories, as well as vibrational spectra. Computational cost savings of 58% are demonstrated for SCF methods, even relative to Fock-extrapolated dynamics, and savings are further increased to 71% with RI-MP2. Notably, these timings outperform an idealized estimate of extended-Lagrangian molecular dynamics. The method is subsequently demonstrated on the vibrational absorption spectrum of two NO(+)(H₂O)₃ isomers and is shown to recover the significant width of the shared-proton bands observed experimentally.

The 1,4-phenylenediisocyanide (PDI) dimer serves as an intriguing case of the substituted benzene dimer, as well as a prototype system for self-assembled monolayers of organic isocyanide complexes. Structures and binding energies are explored using recently developed dual-basis second-order Moller–Plesset perturbation theory energies and gradients. The structures are dictated by a combination of dispersion and electrostatics, a combination not properly treated with local or gradient-corrected density functionals. The PDI dimer binds more than twice as strongly as unsubstituted benzene dimers in several configurations, and greater directional specificity between parallel-displaced and T-shaped structures is observed. A rotated-parallel structure is the predicted lowest-energy, gas-phase configuration, in which the isocyanide ligands are staggered on the monomers. Relevant potential energy curves of the dimer are also presented, and insights into PDI monolayer formation on metal surfaces are explored via simple two-body models. Based on the adsorbate interaction alone, a high-coverage configuration and non-vertical tilt are predicted to be favorable, although the total binding for PDI in these configurations is still insufficient to form ordered monolayers, a result consistent with previous experimental findings. Additional phenyl rings (biphenyldiisocyanide, triphenyldiisocyanide) significantly stabilize the interaction and provide the additional dispersion necessary for an ordered monolayer.

Accurate potential energy surfaces for benzene.M complexes (M = Li+, Na+, K+, and NH4+) are obtained using coupled-cluster theory through perturbative triple excitations, CCSD(T). Our computations show that off-axis cation-pi interactions, where the cation is not directly above the aromatic ring, can be favorable and may influence molecular recognition. Even perpendicular, side-on interactions retain 18-32% of their pi-face interaction energy in the gas phase, making their bond strengths comparable to hydrogen bonds in the gas phase. Solvent effects have been explored for each complex using the polarizable continuum model.

This article presents a general computational approach for efficient simulations of anharmonic vibrational spectra in chemical systems. An automated local-mode vibrational approach is presented, which borrows techniques from localized molecular orbitals in electronic structure theory. This approach generates spatially localized vibrational modes, in contrast to the delocalization exhibited by canonical normal modes. The method is rigorously tested across a series of chemical systems, ranging from small molecules to large water clusters and a protonated dipeptide. It is interfaced with exact, grid-based approaches, as well as vibrational self-consistent field methods. Most significantly, this new set of reference coordinates exhibits a well-behaved spatial decay of mode couplings, which allows for a systematic, a priori truncation of mode couplings and increased computational efficiency. Convergence can typically be reached by including modes within only about 4 Å. The local nature of this truncation suggests particular promise for the ab initio simulation of anharmonic vibrational motion in large systems, where connection to experimental spectra is currently most challenging.

A time-reversible, multiple-timestep protocol is presented for ab initio molecular dynamics simulations using correlated, wavefunction-based underlying potentials. The method is motivated by the observation that electron correlation contributions to forces vary on a slower timescale than their Hartree-Fock counterparts. An efficient dynamics algorithm, involving short-timestep Hartree-Fock and long-timestep Moøller-Plesset perturbation theory, is presented and tested. Results indicate stable trajectories and relative speedups comparable to those seen in force field-based multiple-timestep schemes, with the highest efficiency improvement occurring for large systems.

The resolution-of-the-identity (RI) approximation has placed the onus of the cost of a second-order Moller-Plesset (MP2) calculation on the underlying self-consistent field (SCF) calculation for many moderately sized molecules. A dual-basis approach to the SCF calculation, based on previous methods demonstrated for density functional theory, is combined with RI-MP2 calculations, and small basis subsets for cc-pVTZ, cc-pVQZ, and 6-311++G(3df,3pd) are presented. These subsets provide time savings of greater than 90%, with negligible errors in absolute and relative energies, compared to the associated full-basis counterpart. The method is tested with a series of rotational barriers, relative conformational energies of alanine tetrapeptides, as well as the full G3/99 molecular set. RI-MP2 calculations on alanine octapeptides (40 heavy atoms, 3460 basis functions), using cc-pVQZ, are presented. Results improve upon previous methods that diagonalize the virtual space separately.