Taking the plunge: chemical reaction dynamics in liquids

More than just smoke and mirrors: Gas-phase polaritons for optical control of chemistry

Gas-phase molecules are a promising platform to elucidate the mechanisms of action and scope of polaritons for optical control of chemistry. Polaritons arise from the strong coupling of a dipole-allowed molecular transition with the photonic mode of an optical cavity. There is mounting evidence of modified reactivity under polaritonic conditions; however, the complex condensed-phase environment of most experimental demonstrations impedes mechanistic understanding of this phenomenon. While the gas phase was the playground of early efforts in atomic cavity quantum electrodynamics, we have only recently demonstrated the formation of molecular polaritons under these conditions. Studying the reactivity of isolated gas-phase molecules under strong coupling would eliminate solvent interactions and enable quantum state resolution of reaction progress. In this Perspective, we contextualize recent gas-phase efforts in the field of polariton chemistry and offer a practical guide for experimental design moving forward.

Modelling chemical processes in explicit solvents with machine learning potentials

Solvent effects influence all stages of the chemical processes, modulating the stability of intermediates and transition states, as well as altering reaction rates and product ratios. However, accurately modelling these effects remains challenging. Here, we present a general strategy for generating reactive machine learning potentials to model chemical processes in solution. Our approach combines active learning with descriptor-based selectors and automation, enabling the construction of data-efficient training sets that span the relevant chemical and conformational space. We apply this strategy to investigate a Diels-Alder reaction in water and methanol. The generated machine learning potentials enable us to obtain reaction rates that are in agreement with experimental data and analyse the influence of these solvents on the reaction mechanism. Our strategy offers an efficient approach to the routine modelling of chemical reactions in solution, opening up avenues for studying complex chemical processes in an efficient manner.

Factors that Impact Photochemical Cage Escape Yields

The utilization of visible light to mediate chemical reactions in fluid solutions has applications that range from solar fuel production to medicine and organic synthesis. These reactions are typically initiated by electron transfer between a photoexcited dye molecule (a photosensitizer) and a redox-active quencher to yield radical pairs that are intimately associated within a solvent cage. Many of these radicals undergo rapid thermodynamically favored "geminate" recombination and do not diffuse out of the solvent cage that surrounds them. Those that do escape the cage are useful reagents that may undergo subsequent reactions important to the above-mentioned applications. The cage escape process and the factors that determine the yields remain poorly understood despite decades of research motivated by their practical and fundamental importance. Herein, state-of-the-art research on light-induced electron transfer and cage escape that has appeared since the seminal 1972 review by J. P. Lorand entitled "The Cage Effect" is reviewed. This review also provides some background for those new to the field and discusses the cage escape process of both homolytic bond photodissociation and bimolecular light induced electron transfer reactions. The review concludes with some key goals and directions for future research that promise to elevate this very vibrant field to even greater heights.

Laser induced temperature-jump time resolved IR spectroscopy of zeolites

Combining pulsed laser heating and time-resolved infrared (TR-IR) absorption spectroscopy provides a means of initiating and studying thermally activated chemical reactions and diffusion processes in heterogeneous catalysts on timescales from nanoseconds to seconds. To this end, we investigated single pulse and burst laser heating in zeolite catalysts under realistic conditions using TR-IR spectroscopy. 1 ns, 70 μJ, 2.8 μm laser pulses from a Nd:YAG-pumped optical parametric oscillator were observed to induce temperature-jumps (T-jumps) in zeolite pellets in nanoseconds, with the sample cooling over 1-3 ms. By adopting a tightly focused beam geometry, T-jumps as large as 145 °C from the starting temperature were achieved, demonstrated through comparison of the TR-IR spectra with temperature dependent IR absorption spectra and three dimensional heat transfer modelling using realistic experimental parameters. The simulations provide a detailed understanding of the temperature distribution within the sample and its evolution over the cooling period, which we observe to be bi-exponential. These results provide foundations for determining the magnitude of a T-jump in a catalyst/adsorbate system from its absorption spectrum and physical properties, and for applying T-jump TR-IR spectroscopy to the study of reactive chemistry in heterogeneous catalysts.

The hydrogen-bonding dynamics of water to a nitrile-functionalized electrode is modulated by voltage according to ultrafast 2D IR spectroscopy

We report the hydrogen-bonding dynamics of water to a nitrile-functionalized and plasmonic electrode surface as a function of applied voltage. The surface-enhanced two-dimensional infrared spectra exhibit hydrogen-bonded and non-hydrogen-bonded nitrile features in similar proportions, plus cross peaks between the two. Isotopic dilution experiments show that the cross peaks arise predominantly from chemical exchange between hydrogen-bonded and non-hydrogen-bonded nitriles. The chemical exchange rate depends upon voltage, with the hydrogen bond of the water to the nitriles breaking 2 to 3 times slower (>63 vs. 25 ps) under a positive as compared to a negative potential. Spectral diffusion created by hydrogen-bond fluctuations occurs on a ~1 ps timescale and is moderately potential-dependent. Timescales from molecular dynamics simulations agree qualitatively with the experiment and show that a negative voltage causes a small net displacement of water away from the surface. These results show that the voltage applied to an electrode can alter the timescales of solvent motion at its interface, which has implications for electrochemically driven reactions.

Ultrafast dynamics of CN radical reactions with chloroform solvent under vibrational strong coupling

Polariton chemistry may provide a new means to control molecular reactivity, permitting remote, reversible modification of reaction energetics, kinetics, and product yields. A considerable body of experimental and theoretical work has already demonstrated that strong coupling between a molecular vibrational mode and the confined electromagnetic field of an optical cavity can alter chemical reactivity without external illumination. However, the mechanisms underlying cavity-altered chemistry remain unclear in large part because the experimental systems examined previously are too complex for detailed analysis of their reaction dynamics. Here, we experimentally investigate photolysis-induced reactions of cyanide radicals with strongly-coupled chloroform (CHCl3) solvent molecules and examine the intracavity rates of photofragment recombination, solvent complexation, and hydrogen abstraction. We use a microfluidic optical cavity fitted with dichroic mirrors to facilitate vibrational strong coupling (VSC) of the C-H stretching mode of CHCl3 while simultaneously permitting optical access at visible wavelengths. Ultrafast transient absorption experiments performed with cavities tuned on- and off-resonance reveal that VSC of the CHCl3 C-H stretching transition does not significantly modify any measured rate constants, including those associated with the hydrogen abstraction reaction. This work represents, to the best of our knowledge, the first experimental study of an elementary bimolecular reaction under VSC. We discuss how the conspicuous absence of cavity-altered effects in this system may provide insights into the mechanisms of modified ground state reactivity under VSC and help bridge the divide between experimental results and theoretical predictions in vibrational polariton chemistry.

Spectrophotometric Concentration Analysis Without Molar Absorption Coefficients by Two-Dimensional-Infrared and Fourier Transform Infrared Spectroscopy

A spectrophotometric method for determining relative concentrations of infrared (IR)-active analytes with unknown concentration and unknown molar absorption coefficient is explored. This type of method may be useful for the characterization of complex/heterogeneous liquids or solids, the study of transient species, and for other scenarios where it might be difficult to gain concentration information by other means. Concentration ratios of two species are obtained from their IR absorption and two-dimensional (2D)-IR diagonal bleach signals using simple ratiometric calculations. A simple calculation framework for deriving concentration ratios from spectral data is developed, extended to IR-pump-probe signals, and applied to the calculation of transition dipole ratios. Corrections to account for the attenuation of the 2D-IR signal caused by population relaxation, spectral overlap, wavelength-dependent pump absorption, inhomogeneous broadening, and laser intensity variations are described. A simple formula for calculating the attenuation of the 2D-IR signal due to sample absorption is deduced and by comparison with 2D-IR signals at varying total sample absorbance found to be quantitatively accurate. 2D-IR and Fourier transform infrared spectroscopy of two carbonyl containing species acetone and N-methyl-acetamide dissolved in D2O are used to experimentally confirm the validity of the ratiometric calculations. Finally, to address ambiguities over units and scaling of 2D-IR signals, a physical unit of 2D-IR spectral amplitude in mOD/c m - 1 is proposed.

Reaction dynamics of Diels-Alder reactions from machine learned potentials

Recent advances in the development of reactive machine-learned potentials (MLPs) promise to transform reaction modelling. However, such methods have remained computationally expensive and limited to experts. Here, we employ different MLP methods (ACE, NequIP, GAP), combined with automated fitting and active learning, to study the reaction dynamics of representative Diels-Alder reactions. We demonstrate that the ACE and NequIP MLPs can consistently achieve chemical accuracy (±1 kcal mol^-1) to the ground-truth surface with only a few hundred reference calculations. These strategies are shown to enable routine ab initio-quality classical and quantum dynamics, and obtain dynamical quantities such as product ratios and free energies from non-static methods. For ambimodal reactions, product distributions were found to be strongly dependent on the QM method and less so on the type of dynamics propagated.

Heavy-Atom Kinetic Isotope Effects: Primary Interest or Zero Point?

Chemists have many options for elucidating reaction mechanisms. Global kinetic analysis and classic transition-state probes (e.g., LFERs, Eyring) inevitably form the cornerstone of any strategy, yet their application to increasingly sophisticated synthetic methodologies often leads to a wide range of indistinguishable mechanistic proposals. Computational chemistry provides powerful tools for narrowing the field in such cases, yet wholly simulated mechanisms must be interpreted with great caution. Heavy-atom kinetic isotope effects (KIEs) offer an exquisite but underutilized method for reconciling the two approaches, anchoring the theoretician in the world of calculable observables and providing the experimentalist with atomistic insights. This Perspective provides a personal outlook on this synergy. It surveys the computation of heavy-atom KIEs and their measurement by NMR spectroscopy, discusses recent case studies, highlights the intellectual reward that lies in alignment of experiment and theory, and reflects on the changes required in chemical education in the area.

Effect of solvent motions on the dynamics of the Diels-Alder reaction

How solvent motions affect the dynamics of chemical reactions in which the solute undergoes a substantial shape change is a fundamental but elusive issue. This work utilizes reactive simulation and Grote-Hynes theory to explore the effect of solvent motions on the dynamics of the Diels-Alder reaction (in the reverse direction, this reaction involves very substantial solute expansion) in aprotic solvents. The results reveal that the solvent environment is not sufficiently constraining to influence transition state passage dynamics, with the calculated transmission coefficients being close to unity. Even when solvent motions are suppressed or artificially slowed down, the solvent only affects the reaction dynamics in the transition state region to a very small extent. The only notable effect of solvent occurs far from the transition state region and corresponds to caging of the reactants within the reactant well.

Solvent Effects on Ultrafast Photochemical Pathways

Photochemical reactions are increasingly being used for chemical and materials synthesis, for example, in photoredox catalysis, and generally involve photoexcitation of molecular chromophores dissolved in a liquid solvent. The choice of solvent influences the outcomes of the photochemistry because solute-solvent interactions modify the energies of and crossings between electronic states of the chromophores, and they affect the evolving structures of the photoexcited molecules. Ultrafast laser spectroscopy methods with femtosecond to picosecond time resolution can resolve the dynamics of these photoexcited molecules as they undergo structural and electronic changes, relax back to the ground state, dissipate their excess internal energy to the surrounding solvent, or undergo photochemical reactions. In this Account, we illustrate how experimental studies using ultrafast lasers can reveal the influences that different solvents or cosolutes exert on the photoinduced nonadiabatic dynamics of internal conversion and intersystem crossing in nonradiative relaxation pathways. Although the environment surrounding a solute molecule is rapidly changing, with fluctuations in the coordination to neighboring solvent molecules occurring on femtosecond or picosecond time scales, we show that it is possible to photoexcite selectively only those molecular chromophores transiently experiencing specific solute-solvent interactions such as intermolecular hydrogen bonding.The effects of different solvation environments on the photodynamics are illustrated using four selected examples of photochemical processes in which the solvent has a marked effect on the outcomes. We first consider two aromatic carbonyl compounds, benzophenone and acetophenone, which are known to undergo fast intersystem crossing to populate the first excited triplet state on time scales of a few picoseconds. We show that the nonadiabatic excited-state dynamics are modified by transient hydrogen bonding of the carbonyl group to a protic solvent or by coordination to a metal cation cosolute. We then examine how different solvents modify the competition between two alternative relaxation pathways in a photoexcited UVA-sunscreen molecule, diethylamino hydroxybenzoyl hexyl benzoate (DHHB). This relaxation back to the ground electronic state is an essential part of the effective operation of the sunscreen compound, but the dynamics are sensitive to the surrounding environment. Finally, we consider how solvents of different polarity affect the energies and lifetimes of excited states with locally excited or charge-transfer character in heterocyclic organic compounds used as excited-state electron donors for photoredox catalysis. With these and other examples, we seek to develop a molecular level understanding of how the choice of solution environment might be used to control the outcomes of photochemical reactions.

Mapping the structural dynamics of water dissociation

The formation of transient water photolysis products is followed with ultrafast electron diffraction.

Molecules, the Ultimate Nanomotor: Linking Chemical Reaction Intermediates to their Molecular Diffusivity

The intellectual community focused on nanomotors has recently become interested in extending these concepts to individual molecules. Here, we study a chemical reaction according to whose mechanism some intermediate species should speed up while others slow down in predictable ways, if the nanomotor hypothesis of boosted diffusion holds. Accordingly, we scrutinize the absolute diffusion coefficient (D) during intermediate steps of the catalytic cycle for the CuAAC reaction (copper-catalyzed azide-alkyne cycloaddition click reaction), using proton pulsed field-gradient nuclear magnetic resonance to discriminate between the diffusion of various reaction intermediates. We observe time-dependent diffusion that is enhanced for some intermediate molecular species and depressed for those whose size increases owing to complex formation. These findings point to the failure of the conventional Stokes-Einstein equation to fully explain diffusivity during chemical reaction. Without attempting a firm explanation, this paper highlights aspects of the physics of chemical reactions that are imperfectly understood and presents systematic data that can be used to assess hypotheses.

Impact of dimensionality and confinement on reaction dynamics and thermodynamics within 1D and 2D nanostructures

Confinement has been shown to contribute to the dynamics of small molecules within nanoscale hydrophobic or hydrophilic cavities. Enclosure within a confined space can also influence energy transfer pathways, such as the enhancement of fluorescence over thermal relaxation. In this paper, the effect of confinement on the thermodynamic properties and reaction kinetics of small hydrophobic molecules confined in a soft polymeric template is detailed. A quasi-elastic neutron scattering experiment identified a substantial decrease in translational diffusion of pyrrole after solubilization within a hydrophobic cavity. This decrease in mobility is due to pyrrole's closer packing and increased density under confinement vs the bulk liquid. The decreased mobility and increased density explain the spontaneous polymerization reaction of pyrrole observed within the cavity. The precise characterization of the polymerization kinetics under confinement found that the reaction is independent of pyrrole concentration, consistent with the close packing density. Kinetic data also show that confinement dimensionality finds a thermodynamic expression in the transition state entropy. The dynamics and kinetics experiments reported here offer rare empirical insight into the important influence that cavity geometry places on the reactions they host.

Vibrational Energy Relaxation of Deuterium Fluoride in d-Dichloromethane: Insights from Different Potentials

Vibrationally excited deuterium fluoride (DF) formed by fluorine atom reaction with a solvent was found (Science, 2015, 347, 530) to relax rapidly (less than 10 ps) in acetonitrile-d₃ (CD₃CN) and dichloromethane-d₂ (CD₂Cl₂). However, insights into how CD₂Cl₂ facilitates this energy relaxation have so far been lacking, given the weak interaction between DF and a single CD₂Cl₂. In this work, we report the results of reactive simulations with a two-state reactive empirical valence bond (EVB) potential to study the energy deposited into nascent DF after transition-state passage and of nonequilibrium molecular dynamics simulations using multiple different potential energy functions to model the relaxation dynamics. For these second simulations, we used the standard Merck molecular force field (MMFF) potential, an MMFF-based covalent-ionic empirical valence bond (EVB) potential (EVB^CI), a newly developed potential [referred to as MMFF(rDF)] which extends upon the MMFF potential by making the DF/CD₂Cl₂ interaction depend on the value of the D-F bond stretching coordinate and by taking the anisotropic charge distribution of the solvent molecules into account, the polarizable atomic multipole optimized energetics for biomolecular applications (AMOEBA) potential, and the quantum mechanics/molecular mechanics (QM/MM) potential. The relaxation is revealed to be highly sensitive to the potential used. Neither standard MMFF nor EVB^CI reproduces the experimentally observed rapid relaxation dynamics, and they also fail to provide a good description of the interaction potential between DF and CD₂Cl₂ as calculated using CCSD(T)-F12. This is attributed to the use of a point-charge model for the solute and to failing to model the anisotropic electrostatic properties of CD₂Cl₂. The MMFF(rDF), AMOEBA, and QM/MM potentials all reproduce the CCSD(T)-F12 two-body DF---CD₂Cl₂ interaction potential rather well but only with the QM/MM approach is fast vibrational relaxation obtained (lifetimes of ∼288, ∼186, and ∼8 ps, respectively), which we attribute to differences in the solute-solvent local structure. With QM/MM, a unique "many-body" interaction pattern in which DF is in close contact with two solvent Cl atoms and more than three solvent D atoms is found, but this structure is not seen with other potentials. The QM/MM dynamics also display enhanced solute-solvent interactions with vibrationally excited DF that induce a DF band redshift and hence a resonant overlap with solvent C-D modes, which facilitate the intermolecular energy transfer. Our work also suggests that potentials used to model energy relaxation need to capture the fine structure of solute-solvent interactions and not just the two-body part.

The electronic structure of the aqueous permanganate ion: aqueous-phase energetics and molecular bonding studied using liquid jet photoelectron spectroscopy

Permanganate aqueous solutions, MnO4-(aq.), were studied using liquid-micro-jet-based soft X-ray non-resonant and resonant photoelectron spectroscopy to determine valence and core-level binding energies. To identify possible differences in the energetics between the aqueous bulk and the solution-gas interface, non-resonant spectra were recorded at two different probing depths. Similar experiments were performed with different counter ions, Na+ and K+, with the two solutions yielding indistinguishable anion electron binding energies. Our resonant photoelectron spectroscopy measurements, performed near the Mn LII,III- and O K-edges, selectively probed valence charge distributions between the Mn metal center, O ligands, and first solvation shell in the aqueous bulk. Associated resonantly-enhanced solute ionisation signals revealed hybridisation of the solute constituents' atomic orbitals, including the inner valence Mn 3p and O 2s. We identified intermolecular coulombic decay relaxation processes following resonant X-ray excitation of the solute that highlight valence MnO4-(aq.)-H2O(l) electronic couplings. Furthermore, our results allowed us to infer oxidative reorganisation energies of MnO4˙(aq.) and adiabatic valence ionisation energies of MnO4-(aq.), revealing the Gibbs free energy of oxidation and permitting estimation of the vertical electron affinity of MnO4˙(aq.). Finally, the Gibbs free energy of hydration of isolated MnO4- was determined. Our results and analysis allowed a near-complete binding-energy-scaled MnO4-(aq.) molecular orbital and a valence energy level diagram to be produced for the MnO4-(aq.)/MnO4˙(aq.) system. Cumulatively, our mapping of the aqueous-phase electronic structure of MnO4- is expected to contribute to a deeper understanding of the exceptional redox properties of this widely applied aqueous transition-metal complex ion.

Accessing the Activation Mechanisms of Ethylene Photo-Polymerization under Pressure by Transient Infrared Absorption Spectroscopy

The ambient temperature photoinduced polymerization of compressed (P < 1 GPa) fluid ethylene was characterized by transient infrared absorption spectroscopy with a resolution of few nanoseconds, 3 orders of magnitude higher than previously reported. The reaction has been studied under both one- and two-photon excitation evidencing in the latter case its occurrence only in the presence of different transition metal oxides. Their photocatalytic activity is ascribed to the stabilization of the excited biradicals through electron density exchange between the d orbitals of the metal and the π antibonding orbitals of ethylene which lengthens the lifetime of the biradicals. In both one- and two-photon activation cases the polymerization is characterized by an initial step distinguished by a molecularity of 0.15 ± 0.02 identified as the activation step of the reaction lasting, in the one-photon excitation case, a few hundreds of nanoseconds. Using pulsed excitation the reaction evolves toward a free radical polymerization only under one-photon excitation whereas the critical concentration of radicals required to propagate the reaction is never achieved in the two-photon excitation case. Comparison with continuous wave excitation unambiguously identifies in the average power released to the sample the key factor to drive quantitatively and qualitatively the polymerization.

Direct observation of o-benzyne formation in photochemical hexadehydro-Diels-Alder (hν-HDDA) reactions

Reactive ortho-benzyne derivatives are believed to be the initial products of liquid-phase [4 + 2]-cycloadditions between a 1,3-diyne and an alkyne via what is known as a hexadehydro-Diels–Alder (HDDA) reaction. The UV/VIS spectroscopic observation of o-benzyne derivatives and their photochemical dynamics in solution, however, have not been reported previously. Herein, we report direct UV/VIS spectroscopic evidence for the existence of an o-benzyne in solution, and establish the dynamics of its formation in a photoinduced reaction. For this purpose, we investigated a bis-diyne compound using femtosecond transient absorption spectroscopy in the ultraviolet/visible region. In the first step, we observe excited-state isomerization on a sub-10 ps time scale. For identification of the o-benzyne species formed within 50–70 ps, and the corresponding photochemical hexadehydro-Diels–Alder (hν-HDDA) reactions, we employed two intermolecular trapping strategies. In the first case, the o-benzyne was trapped by a second bis-diyne, i.e., self-trapping. The self-trapping products were then identified in the transient absorption experiments by comparing their spectral features to those of the isolated products. In the second case, we used perylene for trapping and reconstructed the spectrum of the trapping product by removing the contribution of irrelevant species from the experimentally observed spectra. Taken together, the UV/VIS spectroscopic data provide a consistent picture for o-benzyne derivatives in solution as the products of photo-initiated HDDA reactions, and we deduce the time scales for their formation.

We report the transient ultraviolet/visible absorption spectrum of an o-benzyne species in solution for the first time.

Intramolecular vibrational energy redistribution and the quantum ergodicity transition: a phase space perspective

The onset of facile intramolecular vibrational energy flow can be related to features in the connected network of anharmonic resonances in the classical phase space.

When classical trajectories get to quantum accuracy: II. The scattering of rotationally excited H2 on Pd(111)

The classical trajectory method in a quantum spirit assigns statistical weights to classical paths on the basis of two semiclassical corrections: Gaussian binning and the adiabaticity correction.

Ultrafast trans → cis Photoisomerization Dynamics of Alkyl-Substituted Stilbenes in a Supramolecular Capsule

Ultrafast spectroscopy reveals the effects of confinement on the excited-state photoisomerization dynamics for a series of alkyl-substituted trans-stilbenes encapsulated in the hydrophobic cavity of an aqueous supramolecular organic host-guest complex. Compared with the solvated compounds, encapsulated trans-stilbenes have broader excited-state absorption spectra, excited-state lifetimes that are 3-4 times longer, and photoisomerization quantum yields that are 1.7-6.5 times lower in the restricted environment. The organic capsule disrupts the equilibrium structure and restricts torsional rotation around the central C═C double bond in the excited state, which is an important motion for the relaxation of trans-stilbene from S₁ to S₀. The location and identity of alkyl substituents play a significant role in determining the excited-state dynamics and photoisomerization quantum yields by tuning the relative crowding inside the capsule. The results are discussed in terms of distortions of the ground- and excited-state potential energy surfaces, including the topology of the S₁-S₀ conical intersection.

Assessment of reactivities with explicit and implicit solvent models: QM/MM and gas-phase evaluation of three different Ag-catalysed furan ring formation routes

A representative reaction illustrates cases where strong solvent–solute interactions can be sufficiently well captured by continuum solvation model rendering QM/MM unnecessary.

EVB and polarizable MM study of energy relaxation in fluorine–acetonitrile reactions

Many-body effects can impact on rates of energy transfer from a ‘hot’ DF solute to acetonitrile solvent.

Time-delayed mass spectrometry of the low-energy electron impact with a liquid beam surface

We set up an experimental apparatus to investigate the low-energy electron impact with a liquid beam surface, in which a cylindrical liquid beam with a diameter of 25 μm emits as the laminar flow from a microjet and the positively charged ions produced by the electron-impact ionizations are detected with a linear time-of-flight mass spectrometer. We propose a time-delayed mass spectrometry for this apparatus to identify the cationic fragments produced on the liquid surface, in which the application of the ion extracting pulse is delayed with different time intervals after the electron beam pulse. Sensitivity and specificity of the present methodology are demonstrated by the combinational experiments of the gas-phase and liquid ethanol. In comparison with the gas-phase experiments, the ion peaks become much broader in the mass spectra of the liquid beam, primarily due to the molecular evaporation and diffusion. After delaying with about 2 μs, we find that the hydrocarbon ions are ultimately the predominant products in the mass spectra of the liquid ethanol and they are proposed to be produced on the liquid surface. Above observations are in line with the widely accepted picture of the molecular orientation on the liquid surface; namely, the ethanol's CH₃-CH₂- group on the liquid surface prefers to be oriented outside. Therefore, we demonstrate a new mass spectrometry to explore the molecular structures of the liquid surface.

Alkyl hydrogen atom abstraction reactions of the CN radical with ethanol

We present a study of the abstraction of alkyl hydrogen atoms from the β and α positions of ethanol by the CN radical in solution using the Empirical Valence Bond (EVB) method. We have built separate 2 × 2 EVB models for the H_β and H_α reactions, where the atom transfer is parameterized using ab initio calculations. The intra- and intermolecular potentials of the reactant and product molecules were modelled with the General AMBER Force Field, with some modifications. We have carried out the dynamics in water and chloroform, which are solvents of contrasting polarity. We have computed the potential of mean force for both abstractions in each of the solvents. They are found to have a small and early barrier along the reaction coordinate with a large energy release. Analyzing the solvent structure around the reaction system, we have found two solvents to have little effect on either reaction. Simulating the dynamics from the transition state, we also study the fate of the energies in the HCN vibrational modes. The HCN molecule is born vibrationally hot in the CH stretch in both reactions and additionally in the HCN bends for the H_α abstraction reaction. In the early stage of the dynamics, we find that the CN stretch mode gains energy at the expense of the energy in CH stretch mode.