The hydration of the OH radical: Microsolvation modeling and statistical mechanics simulation

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.

Role of Hydroxyl Radical in Degradation of NTO: DFT Study

Hydroxyl radicals are important reactive oxygen species produced in the aquatic environment under sunlight irradiation. Many organic pollutants may be decomposed as they encounter hydroxyl radicals, due to their high oxidative ability. NTO (5-nitro-1,2,4-triazol-3-one), an energetic material used in military applications, may be released to the environment and dissolved in surface water and groundwater due to its good water solubility. A detailed investigation of the possible mechanism for NTO decomposition in water induced by hydroxyl radical as one of the pathways for NTO environmental degradation was performed by computational study at the PCM/M06-2X/6-311++G(d,p) level. Decomposition of NTO was found to be a multistep process that may begin with an addition of hydroxyl radical to the carbon atom of C═N double bond and consequent release of a nitrite radical. The formed intermediate undergoes a series of chemical transformations that include the attachments of hydroxyl radical to carbon atoms, the transfer of hydrogen to hydroxyl radical, isomerization, and bond cleavage, leading to low-weight inorganic compounds, such as ammonia, nitrogen gas, nitrous acid, nitric acid, and carbon(IV) oxide. The anionic form of NTO is more reactive toward interaction with the hydroxyl radical as compared with its neutral form. Calculated activation energies and high exergonicity of the studied process support the significant contribution of the hydroxyl radical to NTO mineralization in environment.

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.

Electrostatics and Chemical Reactivity at the Air-Water Interface

It has been recently discovered that chemical reactions at aqueous interfaces can be orders of magnitude faster compared to conventional bulk phase reactions, but despite its wide-ranging implications, which extend from atmospheric to synthetic chemistry or technological applications, the phenomenon is still incompletely understood. The role of strong electric fields due to space asymmetry and the accumulation of ions at the interface has been claimed as a possible cause from some experiments, but the reorganization of the solvent around the reactive system should provide even greater additional electrostatic contributions that have not yet been analyzed. In this study, with the help of first-principles molecular dynamics simulations, we go deeper into this issue by a careful assessment of solvation electrostatics at the air-water interface. Our simulations confirm that electrostatic forces can indeed be a key factor in rate acceleration compared to bulk solution. Remarkably, the study reveals that the effect cannot simply be attributed to the magnitude of the local electric field and that the fluctuations of the full electrostatic potential resulting from unique dynamical behavior of the solvation shells at the interface must be accounted for. This finding paves the way for future applications of the phenomenon in organic synthesis, especially for charge transfer or redox reactions in thin films and microdroplets.

The water trimer reaction OH + (H2O)3→ (H2O)2OH + H2O

All important stationary points on the potential energy surface (PES) for the reaction OH + (H₂O)₃→ (H₂O)₂OH + H₂O have been fully optimized using the "gold standard" CCSD(T) method with the large Dunning correlation-consistent cc-pVQZ basis sets. Three types of pathways were found. For the pathway without hydrogen abstraction, the barrier height of the transition state (TS1) is predicted to lie 5.9 kcal mol^-1 below the reactants. The two major complexes (H₂O)₃OH (CP1 and CP2a) are found to lie 6.3 and 11.0 kcal mol^-1, respectively, below the reactants [OH + (H₂O)₃]. For one of the H-abstraction pathways the lowest classical barrier height is predicted to be much higher, 6.1 kcal mol^-1 (TS2a) above the reactants. For the other H-abstraction pathway the barrier height is even higher, 15.0 (TS3) kcal mol^-1. Vibrational frequencies and the zero-point vibrational energies connected to the PES are also reported. The energy barriers for the H-abstraction pathways are compared with those for the OH + (H₂O)₂ and OH + H₂O reactions, and the effects of the third water on the energetics are usually minor (0.2 kcal mol^-1).

Structure, stability, infrared spectra, and bonding of OHm(H2O)7 (m = 0, ±1) clusters: ab initio study combining the particle swarm optimization algorithm

The various structural candidates of anionic, neutral, and cationic water clusters OHm(H2O)7 (m = 0, ±1) have been globally predicted by combining the particle swarm optimization method and quantum chemical calculations. Geometry optimization and vibrational analysis for the optimal structures were performed with the MP2/aug-cc-pVDZ method, and the energy profile was further refined at the CCSD(T)/CBS level. Special attention was paid to the relationships between configurations and energies, particularly the first solvation shell coordination number of OH- and OH. For OH-(H2O)7, OH(H2O)7, and OH+(H2O)7 clusters, the most stable species at room temperature are predicted to be the tetra-solvated multi-ring structure A6, the tri-solvated hemibond cage structure N1, and the single five-membered ring structure C2, respectively. The temperature effects on the stability of these three systems were also explored via Gibbs free energies. Furthermore, for the OH-(H2O)7 clusters, the assignments of vibrational transitions in the OH stretching region are in good agreement with the studies of small hydroxide ion-water clusters, and the IR spectra of two isomers (tetra-solvated multi-ring A6 and penta-solvated cage A3) may match future experimental observation well. By topological analysis and reduced density gradient analysis, the structural characteristics and bonding strengths of the studied clusters were investigated. This work indicates the excellent performance of the PSO search algorithm and CALYPSO on water clusters, and may further provide extensive insights into the chemical behavior such as the transport mechanism of OH- ions and OH radicals in the aqueous phase.

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.

One- or Two-Electron Water Oxidation, Hydroxyl Radical, or H2O2 Evolution

Electrochemical or photoelectrochemcial oxidation of water to form hydrogen peroxide (H₂O₂) or hydroxyl radicals (^•OH) offers a very attractive route to water disinfection, and the first process could be the basis for a clean way to produce hydrogen peroxide. A major obstacle in the development of effective catalysts for these reactions is that the electrocatalyst must suppress the thermodynamically favored four-electron pathway leading to O₂ evolution. We develop a thermochemical picture of the catalyst properties that determine selectivity toward the one, two, and four electron processes leading to ^•OH, H₂O₂, and O₂.

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.

Study of electronic structure and dynamics of interacting free radicals influenced by water

We present a study of electronic structure, stability, and dynamics of interaction and recombination of free radicals such as HO(2) and OH influenced by water. As simple model calculations, we performed ab initio and density functional calculations for the interaction of HO(2) and OH in the presence of water cluster. Results indicate that a significant interaction, overcoming the repulsive Columbic barrier, occurs at a separation distance between the radicals of 5.7 A. This confirms early predictions of the minimum size of molecular dianions stable in the gas phase. It is well known that atomic dianions are unstable in the gas phase but molecular dianions are stable when the size of the molecule is larger than 5.7 A. Ab initio molecular dynamics calculations with Car-Parrinello scheme show that the reaction is very fast and occurs on a time scale of about 1.5 ps. The difference in stability and dynamics of the interacting free radicals on singlet and triplet potential energy surfaces is also discussed.

Computational observation of enhanced solvation of the hydroxyl radical with increased NaCl concentration

Classical molecular dynamics simulations with many-body potentials were carried out to quantitatively determine the effect of NaCl salt concentration on the aqueous solvation and surface concentration of hydroxyl radicals. The potential of mean force technique was used to track the incremental free energy of the hydroxyl radical from the vapor, crossing the air-water interface into the aqueous bulk. Results showed increased NaCl salt concentration significantly enhanced hydroxyl radical solvation, which should significantly increase its accommodation on water droplets. This has been experimentally observed for ozone aqueous accommodation with increased NaI concentration, but, to our knowledge, no experimental study has probed this for hydroxyl radicals. The origin for this effect was found to be very favorable hydroxyl radical-chloride ion interactions, being stronger than those for water-chloride.

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.

The Kohn-Sham density of states and band gap of water: from small clusters to liquid water

Electronic properties of water clusters (H2O)(n), with n=2, 4, 8, 10, 15, 20, and 30 molecules were investigated by sequential Monte Carlo/density-functional theory (DFT) calculations. DFT calculations were carried out over uncorrelated configurations generated by Monte Carlo simulations of liquid water with a reparametrized exchange-correlation functional that reproduces the experimental information on the electronic properties (first ionization energy and highest occupied molecular orbital-lowest unoccupied molecular orbital gap) of the water dimer. The dependence of electronic properties on the cluster size (n) shows that the density of states (DOS) of small water clusters (n>10) exhibits the same basic features that are typical of larger aggregates, such as the mixing of the 3a1 and 1b1 valence bands. When long-ranged polarization effects are taken into account by the introduction of embedding charges, the DOS associated with 3a1 orbitals is significantly enhanced. In agreement with valence-band photoelectron spectra of liquid water, the 1b1, 3a1, and 1b2 electron binding energies in water aggregates are redshifted by approximately 1 eV relative to the isolated molecule. By extrapolating the results for larger clusters the threshold energy for photoelectron emission is 9.6+/-0.15 eV (free clusters) and 10.58+/-0.10 eV (embedded clusters). Our results for the electron affinity (V0=-0.17+/-0.05 eV) and adiabatic band gap (E(G,Ad)=6.83+/-0.05 eV) of liquid water are in excellent agreement with recent information from theoretical and experimental works.

Ab initio and analytical intermolecular potential for ClO–H2O

In recent years, the ClO free radical has been found to play an important role in the ozone removal processes in the atmosphere. In this work, the authors present a potential energy surface scan of the ClO.H2O system with high-level ab initio methods. Because of the existence of low-lying excited states of the ClO.H2O complex and their potential impact on the chemical behavior of the ClO radical in the atmosphere, the authors perform the potential energy surface scan at the CCSD(T)/aug-cc-pVTZ level of theory of both the first excited and ground states. Analytical potentials for both ground and excited states, with the ClO and H2O units held fixed at their optimized geometries and with anisotropic atomic polarizabilities modeling the physics of the unpaired electron in the ClO radical, were built based on a Thole-type model. The two minima of the ClO.H2O complex are recovered by the analytical potential.

Bond-dissociation enthalpies in the gas phase and in organic solvents: Making ends meet

Solvent effects are responsible for the difference between gas- and solution-phase bond-dissociation enthalpies (BDEs), and are thus crucial for understanding reactivity in solution. While solvation effects can be negligible (e.g., in reactions involving carbon-centered radicals), they may be rather significant (e.g., when oxygen-centered radicals are formed). This paper reviews a number of models which have been proposed to deal with the difference between the solvation energetics of a radical and its parent molecule. It is concluded that the radical-solvent interaction may be larger than previously anticipated.

Do Hydroxyl Radical−Water Clusters, OH(H2O)n, n = 1−5, Exist in the Atmosphere?

It has been speculated that the presence of OH(H2O)n clusters in the troposphere could have significant effects on the solar absorption balance and the reactivity of the hydroxyl radical. We have used the G3 and G3B3 model chemistries to model the structures and predict the frequencies of hydroxyl radical/water clusters containing one to five water molecules. The reaction between hydroxyl radical clusters and methane was examined as a function of water cluster size to gain an understanding of how cluster size affects the hydroxyl radical reactivity.