The water–hydroxyl radical complex: A matrix isolation study

The stability of oxygen-centered radicals and its response to hydrogen bonding interactions

The stability of various alkoxy/aryloxy/peroxy radicals, as well as TEMPO and triplet dioxygen (3 O2 ) has been explored at a variety of theoretical levels. Good correlations between RSEtheor and RSEexp are found for hybrid DFT methods, for compound schemes such as G3B3-D3, and also for DLPNO-CCSD(T) calculations. The effects of hydrogen bonding interactions on the stability of oxygen-centered radicals have been probed by addition of a single solvating water molecule. While this water molecule always acts as a H-bond donor to the oxygen-centered radical itself, it can act as a H-bond donor or acceptor to the respective closed-shell parent.

3D Characterization of the Molecular Neighborhood of •OH Radical in High Temperature Water by MD Simulation and Voronoi Polyhedra

Understanding the properties of the ^•OH radical in aqueous environments is essential for biochemistry, atmospheric chemistry, and the development of green chemistry technologies. In particular, the technological applications involve knowledge of microsolvation of the ^•OH radical in high temperature water. In this study, the classical molecular dynamics (MD) simulation and the technique based on the construction of Voronoi polyhedra were used to provide 3D characteristics of the molecular vicinity of the aqueous hydroxyl radical (^•OH_aq). The statistical distribution functions of metric and topological features of solvation shells represented by the constructed Voronoi polyhedra are reported for several thermodynamic states of water, including the pressurized high-temperature liquid and supercritical fluid. Calculations showed a decisive influence of the water density on the geometrical properties of the ^•OH solvation shell in the sub- and supercritical region: with the decreasing density, the span and asymmetry of the solvation shell increase. We also showed that the 1D analysis based on the oxygen-oxygen radial distribution functions (RDFs) overestimates the solvation number of ^•OH and insufficiently reflects the influence of transformations in the hydrogen-bonded network of water on the structure of the solvation shell.

Formation of Para-H2O by Vacuum-UV Photolysis of O2 in Solid Hydrogen: Implication for Astrochemistry

The observation that the ortho to para ratio (OPR) of interstellar H₂O is smaller than 3 is an important yet unresolved subject in astronomy. We irradiated O₂ embedded in solid H₂ at 3 K with vacuum-ultraviolet (VUV) light and observed IR lines associated with para-H₂O (denoted as pH₂O) and nonrotating H₂O-(oH₂)_n (where oH₂ denotes ortho-H₂) but no lines associated with ortho-H₂O (denoted as oH₂O). After maintaining the matrix in darkness for ∼30 h, the amount of pH₂O decreased, accompanied by an increase in H₂O-(oH₂)_n via diffusion of oH₂. After that, the continuous nuclear-spin conversion from oH₂ to para-H₂ (denoted as pH₂) in solid H₂ over time resulted in the conversion of nonrotating H₂O-(oH₂)_n to rotating pH₂O in solid pH₂. The observation of the formation and conversion of pH₂O in our experiment suggests a plausible route in which VUV irradiation of O₂ and H₂ adsorbed on grain surfaces might be responsible for the smaller OPR of interstellar H₂O.

VUV photochemistry and nuclear spin conversion of water and water-orthohydrogen complexes in parahydrogen crystals at 4 K

Samples of H2O, HDO, and D2O were isolated in solid parahydrogen (pH2) matrices and irradiated by vacuum ultraviolet (VUV) radiation at 147 nm. Fourier-Transform Infrared (FTIR) spectra showed a clear depletion of D2O and an enrichment of both HDO and H2O by 147 nm irradiation. These irradiation-dependent changes are attributed to the production of OH and/or OD radicals through photodissociations of H2O, HDO, and D2O. The radicals subsequently react with the hydrogen matrix, leading to the observed enrichment of H2O. No trace of isolated OH or OD was detected in the FTIR spectra, indicating that the OH/OD radicals react with the surrounding matrix hydrogen molecules via quantum tunneling within our experimental timescale. The observed temporal changes in concentrations, especially the increase of HDO concentration during VUV irradiation, can be interpreted by a model with a rapid conversion from orthohydrogen (oH2) to pH2 in water-oH2 complexes upon VUV photodissociation, indicating either the acceleration of the nuclear spin conversion (NSC) of H2 due to the magnetic moment of the intermediate OH/OD radical, or the preferential reaction of the OH/OD radical with a nearby oH2 molecule over other pH2 molecules. We have also identified and quantified an anomalously slow NSC of H2O and D2O complexed with oH2 in solid pH2.

Mode-specific quantum dynamics and kinetics of the hydrogen abstraction reaction OH + H2O → H2O + OH

The synergistic effect between the reactant stretching and bending modes on promoting the reaction.

Tracking the energy flow in the hydrogen exchange reaction OH + H2O → H2O + OH

The prototypical hydrogen exchange reaction OH + H2O → H2O + OH has attracted considerable interest due to its importance in a wide range of chemically active environments. In this work, an accurate global potential energy surface (PES) for the ground electronic state was developed based on ∼44 000 ab initio points at the level of UCCSD(T)-F12a/aug-cc-pVTZ. The PES was fitted using the fundamental invariant-neural network method with a root mean squared error of 4.37 meV. The mode specific dynamics was then studied by the quasi-classical trajectory method on the PES. Furthermore, the normal mode analysis approach was employed to calculate the final vibrational state distribution of the product H2O, in which a new scheme to acquire the Cartesian coordinates and momenta of each atom in the product molecule from the trajectories was proposed. It was found that, on one hand, excitation of either the symmetric stretching mode or the asymmetric stretching mode of the reactant H2O promotes the reaction more than the translational energy, which can be rationalized by the sudden vector projection model. On the other hand, the relatively higher efficacy of exciting the symmetric stretching mode than that of the asymmetric stretching mode is caused by the prevalence of the indirect mechanism at low collision energies and the stripping mechanism at high collision energies. In addition, the initial collision energy turns ineffectively into the vibrational energy of the products H2O and OH while a fraction of the energy transforms into the rotational energy of the product H2O. Fundamental excitation of the stretching modes of H2O results in the product H2O having the highest population in the fundamental state of the asymmetric stretching mode, followed by the ground state and the fundamental state of the symmetric stretching mode.

VUV photochemistry of the H2O⋯CO complex in noble-gas matrices: formation of the OH⋯CO complex and the HOCO radical

VUV photolysis of the H₂O⋯CO complexes leads to the formation of the OH⋯CO radical–molecule complexes and trans-HOCO radicals.

Quasi-classical trajectory studies on the full-dimensional accurate potential energy surface for the OH + H2O = H2O + OH reaction

The first accurate PES for the OH + H₂O reaction is developed by using the permutation invariant polynomial-neural network method to fit ∼48 000 CCSD(T)-F12a/AVTZ calculated points.

I + (H2O)2 → HI + (H2O)OH Forward and Reverse Reactions. CCSD(T) Studies Including Spin-Orbit Coupling

The potential energy profile for the atomic iodine plus water dimer reaction I + (H2O)2 → HI + (H2O)OH has been explored using the "Gold Standard" CCSD(T) method with quadruple-ζ correlation-consistent basis sets. The corresponding information for the reverse reaction HI + (H2O)OH → I + (H2O)2 is also derived. Both zero-point vibrational energies (ZPVEs) and spin-orbit (SO) coupling are considered, and these notably alter the classical energetics. On the basis of the CCSD(T)/cc-pVQZ-PP results, including ZPVE and SO coupling, the forward reaction is found to be endothermic by 47.4 kcal/mol, implying a significant exothermicity for the reverse reaction. The entrance complex I···(H2O)2 is bound by 1.8 kcal/mol, and this dissociation energy is significantly affected by SO coupling. The reaction barrier lies 45.1 kcal/mol higher than the reactants. The exit complex HI···(H2O)OH is bound by 3.0 kcal/mol relative to the asymptotic limit. At every level of theory, the reverse reaction HI + (H2O)OH → I + (H2O)2 proceeds without a barrier. Compared with the analogous water monomer reaction I + H2O → HI + OH, the additional water molecule reduces the relative energies of the entrance stationary point, transition state, and exit complex by 3-5 kcal/mol. The I + (H2O)2 reaction is related to the valence isoelectronic bromine and chlorine reactions but is distinctly different from the F + (H2O)2 system.

Matrix-isolation studies of noncovalent interactions: more sophisticated approaches

Noncovalent interactions are crucial for many physical, chemical, and biological phenomena. Matrix isolation is a powerful method to study noncovalent interactions, including hydrogen-bonded species, and it has been extensively used in this field. However, there are difficult situations, such as in the case of species that are impossible to prepare in the gas phase. In this article, we describe some advanced approaches allowing studies of complexes that are problematic for the traditional methods. Photolysis of a suitable precursor in a matrix can lead to a large concentration of 1:1 complexes, which are otherwise very difficult to prepare (e.g., the H2O···O complex). Photolysis of species combined with annealing can lead to complexes of molecules with mobile atoms (e.g., the same H2O···O complex). Simultaneous photolysis of two species combined with annealing can produce complexes of radicals via reactions of the photogenerated complexes with mobile atoms (e.g., the H2O···HCO complex). Interaction of noble-gas (Ng) hydrides with other species is another topic (e.g., the N2···HArF complex) and very large blue shifts of the H-Ng stretching modes are normally observed for these systems. Complexes and dimers of the higher-energy conformer of formic acid have been prepared by using selective vibrational excitation of the ground-state conformer. The higher-energy conformer of formic acid can be efficiently stabilized in the complexes with strong hydrogen bonding. We also consider some problematic cases when the changes in the vibrational frequencies of the 1:1 complexes are very small (e.g., the phenol···Xe complex) and when the complex formation is prevented by strong solvation in the matrix (e.g., species in solid xenon).

The exothermic HCl + OH·(H2O) reaction: removal of the HCl + OH barrier by a single water molecule

The entrance complex, transition state, and exit complex for the title reaction have been investigated using the CCSD(T) method with correlation consistent basis sets up to cc-pVQZ. The stationary point geometries for the reaction are related to but different from those for the water monomer reaction HCl + OH → Cl + H2O. Our most important conclusion is that the hydrogen-bonded water molecule removes the classical barrier entirely. For the endothermic reverse reaction Cl + (H2O)2, the second water molecule lowers the relative energies of the entrance complex, transition state, and exit complex by about 4 kcal/mol. The title reaction is exothermic by 17.7 kcal/mol. The entrance complex HCl⋯OH·(H2O) is bound by 6.9 kcal/mol relative to the separated reactants. The classical barrier height for the reverse reaction is predicted to be 16.5 kcal/mol. The exit complex Cl⋯(H2O)2 is found to lie 6.8 kcal/mol below the separated products. The potential energy surface for the Cl + (H2O)2 reaction is radically different from that for the valence isoelectronic F + (H2O)2 system.

The hemibond as an alternative condensed phase structure for the hydroxyl radical

Despite the critical importance of the hydroxyl radical in major scientific fields, there are still open questions on the behavior of this species in the aqueous phase. In particular, there has been much debate on the existence of a hemibonded interaction between the hydroxyl radical and water molecules. While some reports indicate that the hemibonded radical might explain some experimental data, others have claimed that this interaction is simply a density functional theory (DFT) artifact. Here, we provide results from high level (basis set limit of coupled-cluster levels up to single, double, triple excitations (CCSD(T)) and beyond) ab initio calculations of different OH•(H2O)n clusters in the gas phase to accurately explore the existence of the hemibonded interaction and its energy difference with respect to other well-defined hydrogen bond interactions. Additional comparisons with second order perturbation theory (MP2) and DFT are also presented. Constrained molecular dynamics was applied to determine the free energy for the formation/disruption and ice systems. Overall, our findings confirm that the hemibond can be an alternative structure for the hydroxyl radical in the condensed phase when the formation of hydrogen bonds is impeded. These results will aid the understanding of theoretical and experimental data and help future experimental designs for the detection of this important species.

Effects of a single water molecule on the OH + H2O2 reaction

The effect of a single water molecule on the reaction between H(2)O(2) and HO has been investigated by employing MP2 and CCSD(T) theoretical approaches in connection with the aug-cc-PVDZ, aug-cc-PVTZ, and aug-cc-PVQZ basis sets and extrapolation to an ∞ basis set. The reaction without water has two elementary reaction paths that differ from each other in the orientation of the hydrogen atom of the hydroxyl radical moiety. Our computed rate constant, at 298 K, is 1.56 × 10(-12) cm(3) molecule(-1) s(-1), in excellent agreement with the suggested value by the NASA/JPL evaluation. The influence of water vapor has been investigated by considering either that H(2)O(2) first forms a complex with water that reacts with hydroxyl radical or that H(2)O(2) reacts with a previously formed H(2)O·OH complex. With the addition of water, the reaction mechanism becomes much more complex, yielding four different reaction paths. Two pathways do not undergo the oxidation reaction but an exchange reaction where there is an interchange between H(2)O(2)·H(2)O and H(2)O·OH complexes. The other two pathways oxidize H(2)O(2), with a computed total rate constant of 4.09 × 10(-12) cm(3) molecule(-1) s(-1) at 298 K, 2.6 times the value of the rate constant of the unassisted reaction. However, the true effect of water vapor requires taking into account the concentration of the prereactive bimolecular complex, namely, H(2)O(2)·H(2)O. With this consideration, water can actually slow down the oxidation of H(2)O(2) by OH between 1840 and 20.5 times in the 240-425 K temperature range. This is an example that demonstrates how water could be a catalyst in an atmospheric reaction in the laboratory but is slow under atmospheric conditions.

Excited state properties of liquid water

In this paper, we give an overview of the state of the art in calculations of the electronic band structure and absorption spectra of water. After an introduction to the main theoretical and computational schemes used, we present results for the electronic and optical excitations of water. We focus mainly on liquid water, but spectroscopic properties of ice and vapor phase are also described. The applicability and the accuracy of first-principles methods are discussed, and results are critically presented.

Thioperoxy derivative generated by UV-induced transformation of N-hydroxypyridine-2(1H)-thione isolated in low-temperature matrixes

Photochemical transformations of N-hydroxypyridine-2(1H)-thione and its deuterated isotopologue were studied using the matrix-isolation technique. Low-temperature Ar and N2 matrixes containing monomers of this compound were irradiated with continuous-wave near-UV light. Photogeneration of two products was observed in these experiments. The relative population of these photogenerated species was found to be dependent on the wavelength of the UV light used for irradiation. By comparison of the IR spectra of the photoproducts with the spectra simulated theoretically at the DFT(B3LYP)/6-311++G(d, p) level, the final and the intermediate products were identified as rotameric forms of 2-hydroxysulfanyl-pyridine. This is the first report on generation of this thioperoxy derivative of pyridine. The mechanism of photogeneration of 2-hydroxysulfanyl-pyridine involves a photoinduced cleavage of the N-O bond in N-hydroxypyridine-2(1H)-thione, generation of the .OH radical weakly bound with the remaining pyridylthiyl radical, and recombination of these two radicals by formation of the new -S-O- bond. A theoretical model supporting this interpretation was constructed on the basis of approximate coupled cluster (CC2) calculations of the potential energy surfaces of the ground and first excited singlet electronic states of the system. After electronic excitation of the monomeric N-hydroxypyridine-2(1H)-thione, the molecule evolves to the conical intersection with the potential energy surface of the ground state and then to the global minimum corresponding to 2-hydroxysulfanyl-pyridine.

Electric field effects on aromatic and aliphatic hydrocarbons: a density-functional study

The influence of a uniform static external electric field on some aliphatic and aromatic molecular species is studied within the density functional theory (DFT) employing the 6-311++G(2d,2p) basis set with B3LYP exchange-correlation prescription. The electric field perturbs the molecular geometry but drastically alters the dipole moments and engenders, to a varying degree, the molecular vibrational Stark effect, i.e., shifts in the infrared (IR) vibrational frequencies accompanied by spectral intensity redistribution. For polar molecules, significant negative ("red") and positive ("blue") frequency shifts are observed for field orientations both parallel and antiparallel to their permanent dipole moments. Further, a selective reordering of frontier orbitals is observed to be brought about by moderately intense fields. In particular, molecules having a lowest unoccupied molecular orbital (LUMO) with predominant pi character possess a threshold field beyond which energy gap between the highest occupied molecular orbital (HOMO) and LUMO diminishes rapidly. A time-dependent (TD) DFT analysis reveals that an increase in the applied field strength by and large increases the excitation energies corresponding to significant electronic transitions among frontier MOs with a concomitant decrease in their oscillator strengths.

The structures of ozone and HOx radicals in aqueous solution from combined quantum/classical molecular dynamics simulations

Ozone in aqueous solution decomposes through a complex mechanism that involves initial reaction with a hydroxide ion followed by formation of a variety of oxidizing species such as HO, HO(2), and HO(3) radicals. Though a number of hydrogen-bonded complexes have been described in the gas phase, both theoretically and experimentally, the structures of ozone and HO(x) in liquid water remain uncertain. In this work, combined quantum/classical computer simulations of aqueous solutions of these species have been reported. The results show that ozone undergoes noticeable electron polarization but it does not participate in hydrogen bonds with liquid water. The main contribution of the solvation energy comes from dispersion forces. In contrast, HO(x) radicals form strong hydrogen bonds. They are better proton donors but weaker proton acceptors than water. Their electronic and geometrical structures are significantly modified by the solvent, especially in the case of HO(3). In all cases, fluctuations in amplitudes of electronic properties are considerable, suggesting that solvent effects might play a crucial role on oxidation mechanisms initiated by ozone in liquid water. These mechanisms are important in a broad range of domains, such as atmospheric processes, plant response to ambient ozone, and medical and industrial applications.

Spectroscopic implications of partially quenched orbital angular momentum in the OH-water complex

The OH monomer orbital angular momentum is predicted to be partially quenched in the OH-water complex because of the significant splitting of the OH monomer orbital degeneracy into (2)A' and (2)A' ' electronic states. This orbital angular momentum quenching and the associated decoupling of the electron spin from the a inertial axis are shown to have dramatic effects on the rotational band structure of the microwave and infrared transitions of the OH-water complex. At the ab initio values for the splitting between the (2)A' and (2)A' ' surfaces, simulated spectra of a- and b-type bands, such as those expected for the OH radical stretch and water asymmetric stretch, are predicted to have a noticeably different appearance than the well-established limiting cases associated with fully quenched or completely unquenched orbital angular momentum. Spectral identification of the OH-water complex in the gas phase will require explicit consideration of this quenching phenomenon.

Cell signaling (mechanism and reproductive toxicity): redox chains, radicals, electrons, relays, conduit, electrochemistry, and other medical implications

This article deals with a novel, simple, integrated approach to cell signaling involving basic biochemical principles, and their relationship to reproductive toxicity. Initially, an overview of the biological aspects is presented. According to the hypothetical approach, cell signaling entails interaction of redox chains, involving initiation, propagation, and termination. The messengers are mainly radicals and electrons that are generated during electron transfer (ET) and hydrogen atom abstraction reactions. Termination and initiation processes in the chain occur at relay sites occupied by redox functionalities, including quinones, metal complexes, and imines, as well as redox amino acids. Conduits for the messengers, comprising species with nonbonding electrons, are omnipresent. Details are provided for the various electron transfer processes. In relation to the varying rates of cell communication, rationale is based on electrons and size of radicals. Another fit is similarly seen in inspection of endogenous precursors of reactive oxygen species (ROS); namely, proteins bearing redox moieties, lipid oxidation products, and carbohydrate radicals. A hypothesis is advanced in which electromagnetic fields associated with mobile radicals and electrons play a role. Although radicals have previously been investigated as messengers, the area occupies a minor part of the research, and it has not attracted broad consensus as an important component. For the first time, an integrated framework is presented composed of radicals, electrons, relays, conduits, and electrical fields. The approach is in keeping with the vast majority of experimental observations. Cell signaling also plays an important role in reproductive toxicity. The main classes that cause birth defects, including ROS, radiation, metal compounds, medicinals, abused drugs, and miscellaneous substances, are known to participate in the signaling process. A unifying basis exists, in that both signaling and reproductive toxicity are characterized by the electron transfer-reactive oxygen species-oxidative stress (ET-ROS-OS) scheme. This article also incorporates representative examples of the extensive investigations dealing with various medical implications. There is considerable literature pointing to a role for cell communication in a wide variety of illnesses.