The structure of hydrogen bonded networks is intimately intertwined with their dynamics. Despite the incredibly wide range of hydrogen bond strengths encountered in water clusters, ion−water clusters, and liquid water, we demonstrate that the previously reported correlation between the change in the equilibrium bond length of the hydrogen bonded OH covalent bond and the corresponding shift in its harmonic frequency in water clusters is much more broadly applicable. Surprisingly, this correlation describes the ratios for both the equilibrium OH bond length/ harmonic frequency and the vibrationally averaged bond length/anharmonic frequency in water, hydronium water, and halide water clusters. Consideration of harmonic and anaharmonic data leads to a correlation of −19 ± 1 cm −1 /0.001 Å. The fundamental nature of this correlation is further confirmed through the analysis of ab initio Molecular Dynamics (AIMD) trajectories for liquid water. We demonstrate that this simple correlation for both harmonic and anharmonic systems can be modeled by the response of an OH bond to an external field. Treating the OH bond as a Morse oscillator, we develop analytic expressions, which relate the ratio of the shift in the vibrational frequency of the hydrogen-bonded OH bond to the shift in OH bond length, to parameters in the Morse potential and the ratio of the first and second derivatives of the field-dependent projection of the dipole moment of water onto the hydrogen-bonded OH bond. Based on our analysis, we develop a protocol for reconstructing the AIMD spectra of liquid water from the sampled distribution of the OH bond lengths. Our findings elucidate the origins of the relationship between the molecular structure of the fleeting hydrogen-bonded network and the ensuing dynamics, which can be probed by vibrational spectroscopy.
The ground-state structures of water hexamer and several deuterated variants are studied by using the diffusion Monte Carlo (DMC) method. We demonstrate that a recently developed guided DMC approach allows us to study these systems using substantially smaller ensembles than are required for standard DMC approaches. DMC calculations of the ground states of (H 2 O) 6 and (D 2 O) 6 using the MB-pol potential with 50 000 walkers and a 1 au time step show that for (H 2 O) 6 the cage structure is 51 ± 7 cm −1 lower in energy than the prism structure. In the case of (D 2 O) 6 , the two structures have nearly equal energies (ΔE = 9 ± 11 cm −1 ). The structures of the singly substituted, (H 2 O)(D 2 O) 5 and (H 2 O) 5 (D 2 O), variants of water hexamer are also explored to identify the impact of the location of the unique water molecule on the relative stability of the possible structure. The effect on the stability of the hydrogen bond on whether a hydrogen bond has an OD or OH bond as the donor is also explored. Article pubs.acs.org/JPCA
Infrared (IR) action spectroscopy is utilized to characterize a prototypical carbon-centered hydroperoxyalkyl radical (•QOOH) transiently formed in the oxidation of volatile organic compounds. The •QOOH radical formed in isobutane oxidation, 2-hydroperoxy-2-methylprop-1-yl, •CH2(CH3)2COOH, is generated in the laboratory by H-atom abstraction from tert-butyl hydroperoxide (TBHP). IR spectral features of jet-cooled and stabilized •QOOH radicals are observed from 2950 to 7050 cm−1 at energies that lie below and above the transition state barrier leading to OH radical and cyclic ether products. The observed •QOOH features include overtone OH and CH stretch transitions, combination bands involving OH or CH stretch and a lower frequency mode, and fundamental OH and CH stretch transitions. Most features arise from a single vibrational transition with band contours well simulated at a rotational temperature of 10 K. In each case, the OH products resulting from unimolecular decay of vibrationally activated •QOOH are detected by UV laser-induced fluorescence. Assignments of observed •QOOH IR transitions are guided by anharmonic frequencies computed using second order vibrational perturbation theory, a 2 + 1 model that focuses on the coupling of the OH stretch with two low-frequency torsions, as well as recently predicted statistical •QOOH unimolecular decay rates that include heavy-atom tunneling. Most of the observed vibrational transitions of •QOOH are readily distinguished from those of the TBHP precursor. The distinctive IR transitions of •QOOH, including the strong fundamental OH stretch, provide a general means for detection of •QOOH under controlled laboratory and real-world conditions.
A sparse linear algebra based implementation of Rayleigh–Schrödinger vibrational perturbation theory is presented. This implementation allows for flexibility in the coordinates used to expand the vibrational Hamiltonian as well as the order to which the perturbation theory is performed. It also provides a powerful tool for investigating the origin of spectral intensity and transition frequencies. Specifically, this flexibility allows for the analysis of which terms in the expansions of the Hamiltonian and dipole surface lead to the largest corrections to the energies and transition intensities, and how these conclusions depend on the coordinates used for these expansions. Comparisons of corrections to transition frequencies are reported for the Morse oscillator when the potential is expanded in Δr and Morse coordinates as well as for water, water dimer, and peroxynitrous acid when the molecular Hamiltonians and dipole surfaces are expanded in Cartesian displacement coordinates and in the displacements of the bond-angle-dihedral internal coordinates. Further comparisons of the corrections to the transitions moments are made for H2O and (H2O)2. It is found that while the transition frequencies and intensities are independent of coordinate choice, a good choice of coordinates leads to a cleaner interpretation of the origins of the anharmonicities in these systems.
An approach for idenifying resonances in vibrational perturbation theory calculations is introduced. This approach makes use of the corrections to the wave functions that are obtained from non-degenerate perturbation theory calculations to identify spaces of states that must be treated with degenerate perturbation theory. Pairs of states are considered to be in resonance if the magnitude of expansion coefficients in the corrections to the wave functions in the non-degenerate perturbation theory calculation are greater than a specified threshold, χmax. This approach is applied to calculations of the vibrational spectra of CH4, H2CO, HNO3 and cc-HOONO. The question of how the identified resonances depend on the value of χmax and how the choice of the resonance spaces affects the calculated vibrational spectrum is further explored for \formaldehyde. The approach is also compared to the Martin test [J. Chem. Phys. 103, 2589-2602 (1995)] for calculations of the vibrational spectra of H2CO and cc-HOONO.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.