Spectroscopic Study of the Ion−Radical H-Bond in H4O2+

Spontaneous Oxidation in Aqueous Microdroplets: Water Radical Cation as Primary Oxidizing Agent

Exploration of the unique chemical properties of interfaces can unlock new understanding. A striking example is the finding of accelerated reactions, particularly spontaneous oxidation reactions, that occur without assistance of catalysts or external oxidants at the air interface of both aqueous and organic solutions (provided they contain some water). This finding opened a new area of interfacial chemistry but also caused heated debate regarding the primary chemical species responsible for the observed oxidation. An overview of the literature covering oxidation in microdroplets with air interfaces is provided, together with a critical examination of previous findings and hypotheses. The water radical cation/radical anion pair, formed spontaneously and responsible for the electric field at or near the droplet/air interface, is suggested to constitute the primary redox species. Mechanisms of accelerated microdroplet reactions are critically discussed and it is shown that hydroxyl radical/hydrogen peroxide formation in microdroplets does not require that these species be the primary oxidant. Instead, we suggest that hydroxyl radical and hydrogen peroxide are the products of water radical cation decay in water. The importance of microdroplet chemistry in the prebiotic environment is sketched briefly and the role of partial solvation in reaction acceleration is noted.

Correlation of the Charge Resonance Interaction with Cluster Conformations Probed by Electronic Spectroscopy of Dimer Radical Cations of CO2 and CS2 in a Cryogenic Ion Trap

Radical cations of dimeric clusters of carbon dioxide/disulfide, [(CX2)2]+• (X = O and S), form strong intracluster bonds through charge resonance (CR) interactions. We herein performed electronic photodissociation spectroscopy of [(CX2)2]+• while regulating the temperature under ambient and cryogenic conditions using a quadrupole ion trap. Both ions exhibited broad band absorption in the near-infrared-visible light region; it is called the "CR band", as a measure of the strength of the CR interaction. Strikingly, this band underwent a noticeable blue shift upon cryogenic cooling for [(CS2)2]+• while not for [(CO2)2]+•. On the basis of quantum chemical calculations with a coupled cluster method, the band shift was attributed to the variations in the relative population of two energetically close conformers found for [(CS2)2]+•. This study highlights a strong correlation between CR interactions and conformation of the radical dimer cations, demonstrating the exceptional significance of cryogenic cooling in the chemistry of ionic molecular clusters.

Isolation and Infrared Spectroscopic Characterization of Hemibonded Water Dimer Cation in Superfluid Helium Nanodroplets

The structure of the minimum unit of the radical cationic water clusters, the (H2O)2+ dimer, has attracted much attention because of its importance for the radiation chemistry of water. Previous spectroscopic studies indicated that the dimers have a proton-transferred structure (H3O+·OH), though the alternate metastable hemibonded structure (H2O·OH2)+ was also predicted based on theoretical calculations. Here, we produce (H2O)2+ dimers in superfluid helium nanodroplets and study their infrared spectra in the range of OH stretching vibrations. The observed spectra indicate the coexistence of the two structures in the droplets, supported by density functional theory calculations. This is the first spectroscopic identification of the hemibonded isomer of water radical cation dimers. The observation of the higher-energy isomer reveals efficient kinetic trapping for metastable ionic clusters due to the rapid cooling in helium droplets.

Exploring the Origins of the Intensity of the OH Stretch-HOH Bend Combination Band in Water

While the intensity of the OH stretching fundamental transition is strongly correlated to hydrogen-bond strength, the intensity of the corresponding transition to the state with one quantum of excitation in both the OH stretching and HOH bending vibrations in the same water molecule shows a much weaker sensitivity to the hydrogen-bonding environment. The origins of this difference are explored through analyses of the contributions of terms in the expansion of the dipole moment to the calculated intensity. It is found that the leading contribution to the stretch-bend intensity involves the second derivative of the dipole moment with respect to the OH bond length and HOH angle. While this is not surprising, the insensitivity of this derivative to the hydrogen-bonding environment is unexpected. Possible contributions of mode mixing are also explored. While mode mixing leads to splittings of the energies of nearly degenerate excited states, it does not result in significant changes in the sum of the intensities of these transitions. Analysis of changes in the partial charges on the hydrogen atoms upon displacement of the HOH angles shows that these charges generally increase with increasing HOH angle. This effect is partially canceled by a decrease in the charge of the hydrogen atom when a hydrogen bond is broken. The extent of this cancellation increases with the hydrogen bond strength, which is reflected in the observed insensitivity of the intensity of the stretch-bend transition to hydrogen-bond strength.

Evidence for the co-existence of isomers of water dimer radical cations and their inter-conversion in a linear ion trap

Water dimer radical cations are regarded as key intermediates in many aqueous reactions and biochemical processes. However, the structure of the water dimer radical cations, and particularly the inter-conversion between their isomers, remain difficult to investigate experimentally due to their short lifetime and low abundance under ambient conditions. Furthermore, the isomers cannot be distinguished in a full mass spectra. In this study, we report the experimental evidence for the hemi-bonded and proton-transferred isomers of gas-phase water dimer radical cations, and the inter-conversion process between them in a linear ion trap at low pressure and near room temperature. Multiple collisions of isolated water dimer radical cations with He inside the ion trap were systematically investigated; first, under different trapping times (i.e., reaction times) ranging from 0.03 to 800 ms, and then at a very low collision energies ranging from 0.1% to 10% normalized collision energy. The proton-transferred isomers were dominant at shorter trapping times (≤250 ms), while the hemi-bonded isomers were dominant at longer trapping times (250-800 ms). Moreover, the difference in symmetry of the shapes of the H₂O^•+ signal profiles and the H₃O⁺ signal profiles implied the existence of two kinds of isomers and there were small potential differences between them. Our results also suggested that by tuning the experimental parameters the hemi-bonded isomers would become dominant, which could allow the study of novel chemical reactions involving the hemi-bonded two-center-three-electron (2c-3e) structure in a linear ion trap.

An 241Am Plasma Desorption Ionization (AmDI) Source Scavenged from Smoke Detectors for Ambient Mass Spectrometry Sampling

The development of efficient, low-cost, easy-to-use ambient ionization methods has been a major goal of modern mass spectrometry. In this Letter, we present a gas-free, voltage-free, economic, and safe desorption ionization method using the plasma generated by a radioactive element, americium-241, scavenged from smoke detectors that equip almost every household. No other energy sources, such as laser, discharge, fast-moving carrier gas, solvent droplet, ultrasound, or heat are needed. We name this new method as americium-241 desorption ionization (AmDI). AmDI is tested for the detection of more than 20 volatile and nonvolatile chemicals under different sampling conditions, and the detection limit can be in the range of tens of picograms for some analytes. Mechanistically, we provide evidence that the α particles emitted from radioactive decay ionize ambient air, and the resulting plasma further energizes and ionizes the surface analytes for mass spectrometry detection. We anticipate wide applications of AmDI in mass spectrometric sampling in the near future because of the plethora of merits.

Imaging the short-lived hydroxyl-hydronium pair in ionized liquid water

[Figure: see text].

Infrared Spectroscopy and Anharmonic Vibrational Analysis of (H2O-Krn)+ (n = 1-3): Hemibond Formation of the Water Radical Cation

The hemibond is a nonclassical covalent bond formed between a radical (cation) and a closed shell molecule. The hemibond formation ability of water has attracted great interest, concerning its role in ionization of water. While many computational studies on the water hemibond have been performed, clear experimental evidence has been hardly reported because the hydrogen bond formation overwhelms the hemibond formation. In the present study, infrared photodissociation spectroscopy is applied to (H₂O-Kr_n)⁺ (n = 1-3) radical cation clusters. The observed spectra of (H₂O-Kr_n)⁺ are well reproduced by the anharmonic vibrational simulations based on the hemibonded isomer structures. The firm evidence of the hemibond formation ability of water is revealed.

Photoelectron Spectroscopy of the Water Dimer Reveals Unpredicted Vibrational Structure

Hydrogen bonds and proton transfer reactions can be considered as being at the very heart of aqueous chemistry and of utmost importance for many processes of biological relevance. Nevertheless, these processes are not yet well understood, even in seemingly simple model systems like small water clusters. We present a study of the photoelectron spectrum of the water dimer, revealing previously unresolved vibrational structure with 10-30 meV (80-242 cm^-1) typical splitting, in disagreement with a previous theoretical photoionization study predicting an apparent main vibrational progression with an ∼130 meV spacing [Kamarchik et al.; J. Chem. Phys. 2010, 132, 194311]. The observed vibrational structure and its deviation from the theoretical prediction is discussed in terms of known difficulties with calculations of strongly coupled anharmonic systems involving large amplitude motions. Potential contributions of the nonzero vibrational energy of the neutral water dimer at a finite experimental internal temperature are addressed. The internal temperature is estimated from the breakdown diagram associated with the dissociative ionization of the water dimer to be around to 130 K. This analysis also provides two additional, independently measured values for the 0 K appearance energy of the hydronium ion (H₃O⁺) from dissociative ionization of the water dimer.

Proton Transfer in Nitromethane-Ammonia Clusters under VUV Single-Photon Ionization Explored by Infrared Spectroscopy and Theoretical Calculations

It is known that the acidity and reactivity of the CH bond can be enhanced after ionization. Also, this property plays a pivotal role in proton transfer reaction and in the formation of new molecules. Herein, infrared spectroscopy and high-precision quantum chemical calculations are used to study the neutral and cationic clusters of nitromethane-ammonia (CH₃NO₂-NH₃). It is found that in the neutral cluster, CH₃NO₂ and NH₃ are mainly bonded by three intermolecular hydrogen bonds, in which electrostatic contribution plays a major role. After vacuum ultraviolet (VUV) single-photon ionization of CH₃NO₂-NH₃, the positive charge redistributes from the ionized nitrogen atom of NH₃ to the CH₃NO₂ molecule immediately. Then, the proton of CH₃NO₂ transfers to NH₃ to form a proton-transferred type structure CH₂NO₂-NH₄⁺, without any effective energy barrier, due to the positive hyperconjugation of cationic nitromethane. A closed loop of positive charge transfer takes place in the CH₃NO₂-NH₃ cluster after VUV ionization. The present work demonstrates that both the proton transfer reaction and charge transfer process have occurred in the ionized CH₃NO₂-NH₃ cluster. Moreover, it is found that the proton transfer reaction is a result of the highly acidic CH bond caused by hyperconjugation between the σ (CH) bond and π orbital.

Evaluation of Ar tagging toward the vibrational spectra and zero point energy of X-HOH, X-DOH, and X-HOD, for X = F, Cl, Br

In this study, we theoretically evaluated the effect of argon tagging toward the binding energy and vibrational spectra of water halide anion complexes Ar.X-HOH, Ar.X-HOD, and Ar.X-DOH (X = F, Cl, Br). The ionic hydrogen bond (IHB) OH stretching mode was calculated to have a strong peak in the vibrational spectra, and coupling to intermolecular modes as well as bending modes was observed as combination bands and Fermi resonances. We found that the argon tagging affected the IHB OH stretching peak position in Ar.F-H2O, but not in Ar.Cl-H2O and Ar.Br-H2O. Furthermore, D-binding is favored for Cl and Br based on zero point energies, but for F our calculated zero point energies did not show a preference between H- and D-binding. We show that the competition of the energy lowering in the zero point energy of the anharmonic IHB OH (OD) stretching mode versus the intermolecular out-of-plane IHB OH (OD) wagging mode is important for determining the preference between H- and D-binding for these monohydrated halide clusters. We also found that for X-HOD the HOD bending fundamental peak is blue shifted compared to bare HOD, but is redshifted for F-DOH.

Water Radical Cations in the Gas Phase: Methods and Mechanisms of Formation, Structure and Chemical Properties

Water radical cations, (H₂O)_n^+•, are of great research interest in both fundamental and applied sciences. Fundamental studies of water radical reactions are important to better understand the mechanisms of natural processes, such as proton transfer in aqueous solutions, the formation of hydrogen bonds and DNA damage, as well as for the discovery of new gas-phase reactions and products. In applied science, the interest in water radicals is prompted by their potential in radiobiology and as a source of primary ions for selective and sensitive chemical ionization. However, in contrast to protonated water clusters, (H₂O)_nH⁺, which are relatively easy to generate and isolate in experiments, the generation and isolation of radical water clusters, (H₂O)_n^+•, is tremendously difficult due to their ultra-high reactivity. This review focuses on the current knowledge and unknowns regarding (H₂O)_n^+• species, including the methods and mechanisms of their formation, structure and chemical properties.

Probing the Partial Activation of Water by Open-Shell Interactions, Cl(H2O)1-4

The partial chemical activation of water by reactive radicals was examined computationally for small clusters of chlorine and water, Cl^•(H₂O)_n=1-4. Using an automated isomer-search procedure, dozens of unique, stable structures were computed. Among the resulting structural classes were intact, hydrated-chlorine isomers, as well as hydrogen-abstracted (HCl)(OH)(H₂O)_n-1 configurations. The latter showed increased stability as the degree of hydration increased, until n = 4, where a new class of structures was discovered with a chloride ion bound to an oxidized water network. The electronic structure of these three structural classes was investigated, and spectral signatures of this hydration-based evolution were connected to these electronic properties. An ancillary outcome of this detailed computational analysis, including coupled-cluster benchmarks, was the calibration of cost-effective quantum chemistry methods for future studies of these radical-water complexes.

Activation of CO2 in Photoirradiated CO2-H2O Clusters: Direct Ab Initio Molecular Dynamics (MD) Study

Carbon dioxide (CO₂) is one of the stable and inactive molecules that contribute to greenhouse gases. The development of new reactions of CO₂ activation, chemical fixation, and conversion is a very important issue. In this report, the reactions of CO₂-H₂O binary clusters were investigated using a direct ab initio molecular dynamics (AIMD) method to find a new reaction of CO₂ activation. Clusters composed of carbon dioxide and water molecules, CO₂(H₂O) _n ( n = 2-5), were utilized as a model of the binary cluster. The reaction dynamics of [CO₂(H₂O) _n]⁺ following the ionization of parent neutral clusters were also investigated. Two electronic states of [CO₂(H₂O) _n]⁺ were examined for direct AIMD surfaces: CO₂[(H₂O) _n]⁺ (ground state) and (CO₂)⁺(H₂O) _n (excited charge transfer (CT) state). After the ionization of the clusters, a proton-transfer (PT) reaction occurred within the (H₂O) _n⁺ moiety at the ground state, whereas the reactive HCO₃ radical was formed at the CT state for OH addition to CO₂⁺: CO₂⁺(H₂O) _n → HCO₃ + H⁺(H₂O) _n-1. The mechanisms of the PT process and the HCO₃ radical formation were discussed based on the theoretical results.

Proton transfer from pinene stabilizes water clusters

Molecular beams experiments and ab initio theory reveal indirect formation of protonated water clusters by ionization of pinene.

Infrared photodissociation spectroscopy of ion-radical networks in cationic dimethylamine complexes

Infrared photodissociation spectroscopy was employed to establish the general trends in the stepwise growth motif of cationic dimethylamine (DMA)n+ (n = 4-13) complexes. Electronic structure calculations were performed to identify the structure of the low-lying isomers and to assign the observed spectral features. The results showed the preference of the formation of the proton-transferred (CH3)2NH2+ ion core. The (CH3)2NH2+-[(CH3)2N] ion-radical pair contact and the ion-radical separated pair could coexist at n = 4. The [(CH3)2N] radical is separated from the (CH3)2NH2+ ion core by one DMA molecule at n = 4-6 and by two or more DMA molecules in the larger clusters. This suggests that the (CH3)2NH2+-[(CH3)2N] ion-radical contact pair is not stable in the subsequent radiation-induced processes of DMA, and the [(CH3)2N] radical is released from the charged site in the cationic DMA networks.

Ionization of Ammonia Nanoices with Adsorbed Methanol Molecules

Large ammonia clusters represent a model system of ices that are omnipresent throughout the space. The interaction of ammonia ices with other hydrogen-boding molecules such as methanol or water and their behavior upon an ionization are thus relevant in the astrochemical context. In this study, ammonia clusters (NH₃) _N with the mean size N̅ ≈ 230 were prepared in molecular beams and passed through a pickup cell in which methanol molecules were adsorbed. At the highest exploited pickup pressures, the average composition of (NH₃) _N(CH₃OH) _M clusters was estimated to be N: M ≈ 210:10. On the other hand, the electron ionization of these clusters yielded about 75% of methanol-containing fragments (NH₃) _n(CH₃OH) _mH⁺ compared to 25% contribution of pure ammonia (NH₃) _nH⁺ ions. On the basis of this substantial disproportion, we propose the following ionization mechanism: The prevailing ammonia is ionized in most cases, resulting in NH₄⁺ core solvated most likely with four ammonia molecules, yielding the well-known "magic number" structure (NH₃)₄NH₄⁺. The methanol molecules exhibit a strong propensity for sticking to the fragment ion. We have also considered mechanisms of intracluster reactions. In most cases, proton transfer between ammonia units take place. The theoretical calculations suggested the proton transfer either from the methyl group or from the hydroxyl group of the ionized methanol molecule to ammonia to be the energetically open channels. However, the experiments with selectively deuterated methanols did not show any evidence for the D⁺ transfer from the CD₃ group. The proton transfer from the hydroxyl group could not be excluded entirely or confirmed unambiguously by the experiment.

Vibrational Characterization of Radical Ion Adducts between Imidazole and CO2

We address the molecular level origins of the dramatic difference in the catalytic mechanisms of CO₂ activation by the seemingly similar molecules pyridine (Py) and imidazole (Im). This is accomplished by comparing the fundamental interactions of CO₂ radical anions with Py and Im in the isolated, gas phase PyCO₂^- and ImCO₂^- complexes. These species are prepared by condensation of the neutral compounds onto a (CO₂) _n^- cluster ion beam by entrainment in a supersonic jet ion source. The structures of the anionic complexes are determined by theoretical analysis of their vibrational spectra, obtained by IR photodissociation of weakly bound CO₂ molecules in a photofragmentation mass spectrometer. Although the radical PyCO₂^- system adopts a carbamate-like configuration corresponding to formation of an N-C covalent bond, the ImCO₂^- species is revealed to be best described as an ion-molecule complex in which an oxygen atom in the CO₂^- radical anion is H-bonded to the NH group. Species that feature a covalent N-C interaction in ImCO₂^- are calculated to be locally stable structures, but are much higher in energy than the largely electrostatically bound ion-molecule complex. These results support the suggestion from solution phase electrochemical studies (Bocarsly et al. ACS Catal. 2012, 2, 1684-1692) that the N atoms are not directly involved in the catalytic activation of CO₂ by Im.

Assessment of Real-Time Time-Dependent Density Functional Theory (RT-TDDFT) in Radiation Chemistry: Ionized Water Dimer

Ionization in the condensed phase and molecular clusters leads to a complicated chain of processes with coupled electron-nuclear dynamics. It is difficult to describe such dynamics with conventional nonadiabatic molecular dynamics schemes since the number of states swiftly increases as the molecular system grows. It is therefore attractive to use a direct electron and nuclear propagation such as the real-time time-dependent density functional theory (RT-TDDFT). Here we report a RT-TDDFT benchmark study on simulations of singly and doubly ionized states of a water monomer and dimer as a prototype for more complex processes in a condensed phase. We employed the RT-TDDFT based Ehrenfest molecular dynamics with a generalized gradient approximate (GGA) functional and compared it with wave-function-based surface hopping (SH) simulations. We found that the initial dynamics of a singly HOMO ionized water dimer is similar for both the RT-TDDFT/GGA and the SH simulations but leads to completely different reaction channels on a longer time scale. This failure is attributed to the self-interaction error in the GGA functionals and it can be avoided by using hybrid functionals with large fraction of exact exchange (represented here by the BHandHLYP functional). The simulations of doubly ionized states are reasonably described already at the GGA level. This suggests that the RT-TDDFT/GGA method could describe processes following the autoionization processes such as Auger emission, while its applicability to more complex processes such as intermolecular Coulombic decay remains limited.

General Approach to Coupled Reactive Smoluchowski Equations: Integration and Application of Discrete Variable Representation and Generalized Coordinate Methods to Diffusive Problems

A new and more general approach to diffusion problems with the inclusion of reactivity among different coupled diffusional states is rationalized and presented. The integration of our previous developments in such a field [Phys. Chem. Chem. Phys., 2015, 17, 17362-17374; J. Chem. Theory Comput., 2016, 12, 3482-3490] are implemented in a software package tool allowing the generic user to set up and run diffusional calculations with very low efforts. We show the applicability of the whole framework to a generic diffusional case of chemical interest that is the study case of (N,N-dimethylamino)benzonitrile (DMABN) fluorescence, whose excited state undergoes twisted intramolecular charge transfer (TICT) relaxation. The population dynamics of the excited state coupled to the ground state is followed, and a fluorescence decay spectrum is calculated. The theoretical and numerical background here presented is robust and general enough to complement a wide number of diffusional problems of current interest.

Spectroscopic observation of two-center three-electron bonded (hemi-bonded) structures of (H2S) n+ clusters in the gas phase

A two-center three-electron 2c-3e bond (hemi-bond) is a non-classical chemical bond, and its existence has been supposed in radical cation clusters with lone pairs. Though the nature of the hemi-bond and its role in the reactivity of radical cations have attracted great interest, spectroscopic observations of hemi-bonded structures have been very scarce. In the present study, the presence of a stable hemi-bonded core (H2S∴SH2)+ in (H2S) n+ (n = 3-6) in the gas phase is demonstrated by infrared spectroscopy combined with quantum chemical calculations. The spectral features of the free SH stretch of the ion core show that the hemi-bond motif of the ion core is maintained up to the completion of the first H-bonded solvation shell. All of the observed spectra are well reproduced by the minimum energy hemi-bonded isomers, and no sign of the proton-transferred ion core type H3S+-SH, which is estimated to have a much higher energy, is found. Spin density calculations show that the excess charge is almost equally delocalized over the two H2S molecules in the cluster for n = 3 to 6. This also indicates the hemi-bond nature of the (H2S∴SH2)+ ion core and the small impact of the formation of a solvation shell on the ion core.

Correcting density-driven errors in projection-based embedding

Projection-based embedding provides a simple and numerically robust framework for multiscale wavefunction-in-density-functional-theory (WF-in-DFT) calculations. The approach works well when the approximate DFT is sufficiently accurate to describe the energetics of the low-level subsystem and the coupling between subsystems. It is also necessary that the low-level DFT produces a qualitatively reasonable description of the total density, and in this work, we study model systems where delocalization error prevents this from being the case. We find substantial errors in embedding calculations on open-shell doublet systems in which self-interaction errors cause spurious delocalization of the singly occupied orbital. We propose a solution to this error by evaluating the DFT energy using a more accurate self-consistent density, such as that of Hartree-Fock (HF) theory. These so-called WF-in-(HF-DFT) calculations show excellent convergence towards full-system wavefunction calculations.

Ab initio investigation of structure, stability, thermal behavior, bonding, and infrared spectra of ionized water cluster (H2O)6. J Chem Phys 2016;145:154307. [PMID: 27782468 DOI: 10.1063/1.4964860]

The low-lying isomers of cationic water cluster (H₂O)₆⁺ have been globally explored by using particle swarm optimization algorithm in conjunction with quantum chemical calculations. Compared with previous results, our searching method covers a wide range of structural isomers of (H₂O)₆⁺ and therefore turns out to be more effective. With these local minima, geometry optimization and vibrational analysis are performed for the most interesting clusters at second-order Møller-Plesset (MP2)/aug-cc-pVDZ level, and their energies are further refined at MP2/aug-cc-pVTZ and coupled-cluster theory with single, double, and perturbative triple excitations/aug-cc-pVDZ level. The interaction energies using the complete basis set limits at MP2 level are also reported. The relationships between their structure arrangement and their energies are discussed. Based on the results of thermal simulation, structural change from a four-numbered ring to a tree-like structure occurs at T ≈ 45 K, and the relative population of six lowest-free-energy isomers is found to exceed 4% at some point within the studied temperature range. Studies reveal that, among these six isomers, two new-found isomers constitute 10% of isomer population at 180 K, and the experimental spectra can be better explained with inclusions of the two isomers. The molecular orbitals for six representative cationic water clusters are also studied. Through topological and reduced density gradient analysis, we investigated the structural characteristics and the bonding strengths of these water cluster radical cations.

Quantum molecular motion in the mixed ion-radical complex, [(H2O)(H2S)]. Phys Chem Chem Phys 2016;18:27450-27459. [PMID: 27711703 DOI: 10.1039/c6cp05299a]

The cation dimer of water and hydrogen sulfide, [(H₂O)(H₂S)]⁺, serves as a fundamental model for the oxidation chemistry of H₂S. The known oxidative metabolism of H₂S by biological species in sulfur-rich environments has motivated the study of the inherent properties of this benchmark complex, with possible mechanistic implications for modern water oxidation chemistry. The low-energy isomer of this open-shell ion is a proton-transferred (PT) structure, H₃O⁺SH˙. An alternative PT structure, H₃S⁺OH˙, and a hemibonded (HB) isomer, [H₂O·SH₂]⁺, are also stable isomers, placing this complex intermediate to known (H₂O)₂⁺ (PT) and (H₂S)₂⁺ (HB) limiting regimes. This intermediate character suggested the possibility of unique molecular motion, even in the vibrational ground state. Path integral molecular dynamics and anharmonic vibrational spectroscopy simulations have been performed in this study, in order to understand the inherent quantum molecular motion of this complex. The resulting structural distributions were found to deviate significantly from both classical and harmonic analyses, including the observation of large-amplitude anharmonic motion of the central proton and nearly free rotation of the terminal hydrogens. The predicted vibrational spectra are particularly unique and suggest characteristic signatures of the strong electronic interactions and anharmonic vibrational mode couplings in this radical cation.

Photoelectron Spectroscopy of cis-Nitrous Acid Anion (cis-HONO(-))

We report photoelectron spectra of cis-HONO(-) formed from an association reaction of OH(-) and NO in a pulsed, plasma-entrainment ion source. The experimental data are assigned to the cis-HONO(-) isomer, which is predicted to be the global minimum on the anion potential energy surface. We do not find evidence for a significant contribution from trans-HONO(-). Electron photodetachment of cis-HONO(-) with 1613, 1064, 532, 355, and 301 nm photons accesses the ground X̃ (1)A' (S0) and excited ã (3)A″ (T1) states of neutral HONO. The photoelectron spectrum resulting from detachment forming cis-HONO (S0) exhibits a long vibrational progression, dominated by overtones and combination bands involving the central O-N stretching and ONO bending vibrations. This indicates that there is a significant change in the central O-N bond length between cis-HONO(-) and cis-HONO (S0). The electron affinity (EA) of cis-HONO is determined to be 0.356(8) eV. We also report the dissociation energy (D0) of cis-HONO(-), forming OH(-) + NO, as 0.594(9) eV, which is a factor of 4 decrease in the central O-N bond strength compared to neutral cis-HONO. The T1 state of cis-HONO is shown to be ∼2.3 eV higher in energy than cis-HONO (S0). Electron photodetachment to form cis-HONO (T1) accesses a transition state along the HO-NO bond dissociation coordinate. The resulting photoelectron spectrum exhibits broad peaks spaced by the terminal N═O stretching frequency. Electronic structure calculations and photoelectron spectrum simulations reported here show very good agreement with the experimental data.

The mechanism of excited state proton dissociation in microhydrated hydroxylamine clusters

The dynamics and mechanism of excited-state proton dissociation and transfer in microhydrated hydroxylamine clusters are studied using NH2OH(H2O)n (n = 1-4) as model systems and the DFT/B3LYP/aug-cc-pVDZ and TD-DFT/B3LYP/aug-cc-pVDZ methods as model calculations. This investigation is based on the Förster acidity scheme and emphasizes the photoacid dissociation in the ground (S0) and lowest singlet-excited states (S1) and the interplay between the photo and thermal excitations. The quantum chemical results suggest that the intermediate complexes are formed only in the S1 state in a low local-dielectric environment (e.g., ε = 1) and that upon the S0→ S1 transition, the photon energy excites mostly NH2OH, which leads to a homolytic cleavage of the O-H bond and to dynamically stable charge-separated Rydberg-like H-bond complexes (e.g., NH2O˙-H3O(+)˙). The potential energy surfaces for proton displacement in the smallest Rydberg-like H-bond complex support the intersection of the S0 and S1 states in low local-dielectric environments, whereas in a high local-dielectric environment (e.g., ε = 78), these two states are completely separated. Based on the static results, a photoacid-dissociation mechanism that involves Rydberg-like H-bond complex formation, an H-bond chain extension and fluctuations in the local-dielectric environment is proposed. NVT-BOMD simulations confirm the static results and show that the dynamic behavior of the dissociating proton in the S1 state is not different from that of the protonated H-bond systems in the ground state, which consists of the oscillatory shuttling and structural diffusion motions. These findings allow our theoretical methods, which have been used successfully in protonated H-bond systems in the ground state, to be applied in the study of the photoacid-dissociation processes. The current theoretical study suggests effective steps as well as guidelines for the investigation of the dynamics of the photoacid-dissociation and transfer processes in the Förster acidity scheme, provided that the exciting photon does not lead to a significant change in the structure of the intermediate complex in the excited state.

Water network-mediated, electron-induced proton transfer in [C5H5N ⋅ (H2O)n](-) clusters

The role of proton-assisted charge accommodation in electron capture by a heterocyclic electron scavenger is investigated through theoretical analysis of the vibrational spectra of cold, gas phase [Py ⋅ (H2O)n=3-5](-) clusters. These radical anions are formed when an excess electron is attached to water clusters containing a single pyridine (Py) molecule in a supersonic jet ion source. Under these conditions, the cluster ion distribution starts promptly at n = 3, and the photoelectron spectra, combined with vibrational predissociation spectra of the Ar-tagged anions, establish that for n > 3, these species are best described as hydrated hydroxide ions with the neutral pyridinium radical, PyH((0)), occupying one of the primary solvation sites of the OH(-). The n = 3 cluster appears to be a special case where charge localization on Py and hydroxide is nearly isoenergetic, and the nature of this species is explored with ab initio molecular dynamics calculations of the trajectories that start from metastable arrangements of the anion based on a diffuse, essentially dipole-bound electron. These calculations indicate that the reaction proceeds via a relatively slow rearrangement of the water network to create a favorable hydration configuration around the water molecule that eventually donates a proton to the Py nitrogen atom to yield the product hydroxide ion. The correlation between the degree of excess charge localization and the evolving shape of the water network revealed by this approach thus provides a microscopic picture of the "solvent coordinate" at the heart of a prototypical proton-coupled electron transfer reaction.

Nonadiabatic dynamics of floppy hydrogen bonded complexes: the case of the ionized ammonia dimer

Non-adiabatic dynamics of a floppy hydrogen bonded ammonia dimer was studied by ab initio molecular dynamics simulations.

Discrete variable representation of the Smoluchowski equation using a sinc basis set

We present a new general framework for solving the monodimensional Smoluchowski equation using a discrete variable representation (DVR) based on the so called sinc basis set. The reliability of our implementation is assessed by comparing the convergence of diffusive operator eigenvalues calculated using our method and using a simple finite difference scheme for some model diffusive problems. The results here presented open encouraging possibilities for dealing with more complicated systems, where additional coordinate dependent terms in the equation or multidimensional treatments are needed and traditional methods often become unfeasible.