Attention:The NSF Public Access Repository (PAR) system and access will be unavailable from 5:00 PM ET until 8:00 PM ET on Friday, September 11 due to maintenance. We apologize for the inconvenience.


This content will become publicly available on November 27, 2026

Title: Diffusion Monte Carlo Study of the Structure and Spectroscopy of H 3 O –
A potential energy surface for H3O– has been developed based on the NN+(MOB-ML) approach we developed for studies of complexes of OH– with two or three water molecules. Unlike those systems, H3O– has two low-energy isomers, H–·H2O and OH–·H2, which differ in energy by less than 2.5 kcal mol–1, and which are separated by a barrier of roughly 4.5 kcal mol–1. We find that by training the NN+(MOB-ML) model using structures based on diffusion Monte Carlo (DMC) simulations initiated in the two potential minima, we are able to obtain a potential surface that describes both isomers. Using these potentials, the structure and spectra of H3O– and its deuterated analogues are investigated using DMC. These calculations show that the ground state wave function for H3O– is mainly localized in the H–·H2O minimum in the potential, with a small amount of the probability amplitude (<5%) in the region of the OH–·H2 minimum. The delocalization of the wave function into the secondary minimum is lowered by deuteration of the water molecule, while replacing H– with D– increases the isomerization due to the shortening of the average distance between the hydride ion and the hydrogen atom in water to which it is bound. Introducing one quantum of excitation in the H–···H2O stretching vibration increases the amount of isomerization, while the isomerization decreases with additional excitation of this mode. Excitation of the free OH stretch in the water molecule also increases the amount of isomerization, while excitation of the out-of-plane bending vibration suppresses the isomerization. The effects of partial deuteration on the frequencies for these vibrations are also explored.  more » « less
Award ID(s):
2517458
PAR ID:
10680880
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
American Chemical Society
Date Published:
Journal Name:
The Journal of Physical Chemistry A
Volume:
129
Issue:
47
ISSN:
1089-5639
Page Range / eLocation ID:
10928 to 10939
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Using a combination of two-dimensional infrared (2D IR) and variable temperature Fourier transform infrared (FTIR) spectroscopies the rapid structural isomerization of a five-coordinate ruthenium complex is investigated. In methylene chloride, three exchanging isomers were observed: (1) square pyramidal equatorial, ( 1 ); (2) trigonal bipyramidal, ( 0 ); and (3) square pyramidal apical, ( 2 ). Exchange between 1 and 0 was found to be an endergonic process (Δ H = 0.84 (0.08) kcal mol −1 , Δ S = 0.6 (0.4) eu) with an isomerization time constant of 4.3 (1.5) picoseconds (ps, 10 −12 s). Exchange between 0 and 2 however was found to be exergonic (Δ H = −2.18 (0.06) kcal mol −1 , Δ S = −5.3 (0.3) eu) and rate limiting with an isomerization time constant of 6.3 (1.6) ps. The trigonal bipyramidal complex was found to be an intermediate, with an activation barrier of 2.2 (0.2) kcal mol −1 and 2.4 (0.2) kcal mol −1 relative to the equatorial and apical square pyramidal isomers respectively. This study provides direct validation of the mechanism of Berry pseudorotation – the pairwise exchange of ligands in a five-coordinate complex – a process that was first described over fifty years ago. This study also clearly demonstrates that the rate of pseudorotation approaches the frequency of molecular vibrations. 
    more » « less
  2. The reactions of thioformaldehyde (H 2 CS) with OH radicals and assisted by a single water molecule have been investigated using high level ab initio quantum chemistry calculations. The H 2 CS + ˙OH reaction can in principle proceed through: (1) abstraction, and (2) addition pathways. The barrier height for the addition reaction in the absence of a catalyst was found to be −0.8 kcal mol −1 , relative to the separated reactants, which has a ∼1.0 kcal mol −1 lower barrier than the abstraction channel. The H 2 CS + ˙OH reaction assisted by a single water molecule reduces the barrier heights significantly for both the addition and abstraction channels, to −5.5 and −6.7 kcal mol −1 respectively, compared to the un-catalyzed H 2 CS + ˙OH reaction. These values suggest that water lowers the barriers by ∼6.0 kcal mol −1 for both reaction paths. The rate constants for the H 2 CS⋯H 2 O + ˙OH and OH⋯H 2 O + H 2 CS bimolecular reaction channels were calculated using Canonical Variational Transition state theory (CVT) in conjunction with the Small Curvature Tunneling (SCT) method over the atmospherically relevant temperatures between 200 and 400 K. Rate constants for the H 2 CS + ˙OH reaction paths for comparison with the H 2 CS + ˙OH + H 2 O reaction in the same temperature range were also computed. The results suggest that the rate of the H 2 CS + ˙OH + H 2 O reaction is slower than that of the H 2 CS + ˙OH reaction by ∼1–4 orders of magnitude in the temperatures between 200 and 400 K. For example, at 300 K, the rates of the H 2 CS + ˙OH + H 2 O and H 2 CS + ˙OH reactions were found to be 2.2 × 10 −8 s −1 and 6.4 × 10 −6 s −1 , respectively, calculated using [OH] = 1.0 × 10 6 molecules cm −3 , and [H 2 O] = 8.2 × 10 17 molecules cm −3 (300 K, RH 100%) atmospheric conditions. Electronic structure calculations on the H 2 C(OH)S˙ product in the presence of 3 O 2 were also performed. The results show that H 2 CS is removed from the atmosphere primarily by reacting with ˙OH and O 2 to form thioformic acid, HO 2 , formaldehyde, and SO 2 as the main end products. 
    more » « less
  3. Methods for evaluating vibrational spectra are discussed. These include the harmonic treatment as well as extensions based on vibrational perturbation theory (VPT2). We also describe approaches for obtaining simulated spectra from properties of the ground state wave function. These approaches include spectral mapping approaches in which we exploit a linear relationship between the vibrationally averaged OH bond length and the anharmonic frequency of a hydrogen-bonded OH stretching vibration. An approximation to excited state wave functions that takes the form of products of the ground state wave function and a set of orthogonal polynomials is also described. For these last two approaches, the ground state wave function is obtained from a diffusion Monte Carlo (DMC) calculation. Requirements of the potential surfaces for all of the above calculations are discussed. These approaches are used to explore H+(H2O)n and OH-(H2O)n complexes with n=2 and 3. 
    more » « less
  4. Nitroxyl (HNO) and hydrogen peroxide have both been implicated in a variety of reactions relevant to environmental and physiological processes and may contribute to a unique, unexplored, pathway for the production of nitrous acid (HONO) in soil. To investigate the potential for this reaction, we report an in-depth investigation of the reaction pathway of H 2 O 2 and HNO forming HONO and water. We find the breaking of the peroxide bond and a coupled proton transfer in the first step leads to hydrogen nitryl (HNO 2 ) and an endogenous water, with an extrapolated NEVPT2 (multireference perturbation theory) barrier of 29.3 kcal mol −1 . The first transition state is shown to possess diradical character linking the far peroxide oxygen to the bridging, reacting, peroxide oxygen. The energy of this first step, when calculated using hybrid density functional theory, is shown to depend heavily on the amount of Hartree–Fock exchange in the functional, with higher amounts leading to a higher barrier and more diradical character. Additionally, high amounts of spin contamination cause CCSD(T) to significantly overestimate the TS1 barrier with a value of 36.2 kcal mol −1 when using the stable UHF wavefunction as the reference wavefunction. However, when using the restricted Hartree–Fock reference wavefunction, the TS1 CCSD(T) energy is lowered to yield a barrier of 31.2 kcal mol −1 , in much better agreement with the NEVPT2 result. The second step in the reaction is the isomerization of HNO 2 to trans -HONO through a Grotthuss-like mechanism accepting a proton from and donating a proton to the endogenous water. This new mechanism for the isomerization of HNO 2 is shown to have an NEVPT2 barrier of 23.3 kcal mol −1 , much lower than previous unimolecular estimates not including an explicit water. Finally, inclusion of an additional explicit water is shown to lower the HNO 2 isomerization barrier even further. 
    more » « less
  5. Herein we report the results of preparing metal compounds (where the metal ions are Co2+, Ni2+, Cu2+, Zn2+) with the cyclic ligand 1,4,8,11-tetraazacyclotetradecane [cyclam] under a variety of conditions of metal-ligand ratios and solvent media. In all cases, we used metal Cl2 . nH2O salts (except for anhydrous CoCl2), as specified. Outcome: we isolated species with a four-coordinate metal in the N4 cavity of the ligand alone, and also with either one or two additional axial ligands. Those axial ligands can be (a) a single chloride, leading to penta-coordinated+ products; (b) two chlorides, leading to octahedral-neutral compounds; (c) two waters, giving rise to hexa-coordinated [(cyclam)metal(H2O)2]2+ species. Finally, in the case of HCl added to the reaction medium, the cyclam can be di-protonated and appears as [(cyclam)H2]2+ in the crystals. With such a variety of products, it is not surprising that since the metal coordination numbers vary, the cyclam ligand stereochemistries are thereby affected. Interestingly, the [(cyclam)metal] species are invariably hydrogen-bonded to one another in infinite strings of two kinds: (1) those for which the crystal’s Z’ = 1 have single strings; (2) when Z’ = 2, there is a pair of homogeneous strings attached to one another by a variety of hydrogen-bonding linkages. Finally, we observed an interesting pair of hydroxonium cations: the first is hydoxonium cations in a pleated 2-D sheet consisting of fused pentagons located between sheets of [(cyclam)metal] moieties; the second one is an infinite string of composition (H3O+)-(H2O)-(H3O+)-(H2O)-(H3O+)-(H2O)-(H3O+). 
    more » « less