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

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.

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.

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

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.

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.