Conformational memory in photodissociation of formic acid

The impact of molecular configuration on the bond breaking rates of hydrocarbons: a computational study

The dissociation of hydrocarbon bonds plays a pivotal role in their utilization, whether through fuel combustion or the thermo-cracking of large hydrocarbons in petroleum refinement. Previous studies have primarily focused on the effects of temperature, pressure, and chemical environment on hydrocarbon reactions. However, the influence of molecular configuration on bond breaking rates has not been thoroughly explored. In this study, we propose an approach to compute bond dissociation rates, and apply it to the reactive molecular dynamics simulation (ReaxFF) trajectories of three molecules: n-tridecane, n-pentane, and 1,3-propanediol. Our results reveal that the bond dissociation rate depends not only on the bond position in the chain, but also on the molecular configuration. Stretched configurations exhibit higher dissociation rates, particularly favoring the breaking of central bonds. Conversely, when the molecule is coiled, resulting in a reduced size, terminal bonds exhibit higher dissociation rates. This research contributes to a deeper understanding of molecular dissociation properties in the oxidation of hydrocarbons, and provides a way to explore the bond breaking properties of other molecules.

The Intrinsic Barrier Width and Its Role in Chemical Reactivity

Chemical reactions are in virtually all cases understood and explained on the basis of depicting the molecular potential energy landscape, i.e., the change in atomic positions vs the free-energy change. With such landscapes, the features of the reaction barriers solely determine chemical reactivities. The Marcus dissection of the barrier height (activation energy) on such a potential into the thermodynamically independent (intrinsic) and the thermodynamically dependent (Bell-Evans-Polanyi) contributions successfully models the interplay of reaction rate and driving force. This has led to the well-known and ubiquitously used reactivity paradigm of "kinetic versus thermodynamic control". However, an analogous dissection concept regarding the barrier width is absent. Here we define and outline the concept of intrinsic barrier width and the driving force effect on the barrier width and report experimental as well as theoretical studies to demonstrate their distinct roles. We present the idea of changing the barrier widths of conformational isomerizations of some simple aromatic carboxylic acids as models and use quantum mechanical tunneling (QMT) half-lives as a read-out for these changes because QMT is particularly sensitive to barrier widths. We demonstrate the distinct roles of the intrinsic and the thermodynamic contributions of the barrier width on QMT half-lives. This sheds light on resolving conflicting trends in chemical reactivities where barrier widths are relevant and allows us to draw some important conclusions about the general relevance of barrier widths, their qualitative definition, and the consequences for more complete descriptions of chemical reactions.

Conformational Structures of Neutral and Cationic Pivaldehyde Revealed by IR-Resonant VUV-MATI Mass Spectroscopy

Pivaldehyde, which is an unwanted by-product released with engine exhaust, has received considerable research attention because of its hydrocarbon oxidations at atmospheric temperature. To gain insight into the conformer-specific reaction dynamics, we investigated the conformational structures of the pivaldehyde molecule in neutral (S₀) and cationic (D₀) states using the recently invented IR-resonant VUV-MATI mass spectroscopy. Additionally, we constructed the two-dimensional potential energy surfaces (2D PESs) associated with the conformational transformations in the S₀ and D₀ states to deduce the conformations corresponding to the measured vibrational spectra. The 2D PESs indicated the presence of only the eclipsed conformation in the global minima of both states, unlike those in propanal and isobutanal. However, comparing the IR-dip VUV-MATI spectra from two intense peaks in the VUV-MATI spectrum with the anharmonic IR simulations revealed the correspondence between the gauche conformer on the S₀ state and the measured IR spectra. Furthermore, Franck-Condon analysis confirmed that most peaks in the VUV-MATI spectrum are attributed to the adiabatic ionic transitions between the neutral gauche and cationic eclipsed conformers in pivaldehyde. Consequently, electron removal from the highest occupied molecular orbital, consisting of the nonbonding orbital of the oxygen atom in pivaldehyde, promoted the formyl-relevant modes in the induced cationic eclipsed conformer.

Conformer-Specific Dissociation Dynamics in Dimethyl Methylphosphonate Radical Cation

The dynamics of the dimethyl methylphosphonate (DMMP) radical cation after production by strong field adiabatic ionization have been investigated. Pump-probe experiments using strong field 1300 nm pulses to adiabatically ionize DMMP and a 800 nm non-ionizing probe induce coherent oscillations of the parent ion yield with a period of about 45 fs. The yields of two fragments, PO2C2H7+ and PO2CH4+, oscillate approximately out of phase with the parent ion, but with a slight phase shift relative to each other. We use electronic structure theory and nonadiabatic surface hopping dynamics to understand the underlying dynamics. The results show that while the cation oscillates on the ground state along the P=O bond stretch coordinate, the probe excites population to higher electronic states that can lead to fragments PO2C2H7+ and PO2CH4+. The computational results combined with the experimental observations indicate that the two conformers of DMMP that are populated under experimental conditions exhibit different dynamics after being excited to the higher electronic states of the cation leading to different dissociation products. These results highlight the potential usefulness of these pump-probe measurements as a tool to study conformer-specific dynamics in molecules of biological interest.

IR-induced and tunneling reactions in cryogenic matrices: the (incomplete) story of a successful endeavor

In this article, IR-induced and tunneling-driven reactions observed in cryogenic matrices are described in a historical perspective, the entangling of the two types of processes being highlighted. The story of this still ongoing fascinating scientific endeavor is presented here following closely our own involvement in the field for more than 30 years, and thus focuses mostly on our work. It is, because of this reason, also an incomplete story. Nevertheless, it considers a large range of examples, from very selective IR-induced conformational isomerizations to IR-induced bond-breaking/bond-forming reactions and successful observations of rare heavy atom tunneling processes. As a whole, this article provides a rather general overview of the major progress achieved in the field.

Switching on H-Tunneling through Conformational Control

H-tunneling is a ubiquitous phenomenon, relevant to fields from biochemistry to materials science, but harnessing it for mastering the manipulation of chemical structures still remains nearly illusory. Here, we demonstrate how to switch on H-tunneling by conformational control using external radiation. This is outlined with a triplet 2-hydroxyphenylnitrene generated in an N₂ matrix at 10 K by UV-irradiation of an azide precursor. The anti-orientation of the nitrene's OH moiety was converted to syn by selective vibrational excitation at the 2ν(OH) frequency, thereby moving the H atom closer to the vicinal nitrene center. This triggers spontaneous H-tunneling to a singlet 6-imino-2,4-cyclohexadienone. Computations reveal that such fast H-tunneling occurs through crossing the triplet-to-singlet potential energy surfaces. Our experimental realization provides an exciting novel strategy to attain control over tunneling, opening new avenues for directing chemical transformations.

Roaming Dynamics and Conformational Memory in Photolysis of Formic Acid at 193 nm Using Time-resolved Fourier-transform Infrared Emission Spectroscopy

In photodissociation of trans-formic acid (HCOOH) at 193 nm, we have observed two molecular channels of CO + H₂O and CO₂ + H₂ by using 1 μs-resolved Fourier-transform infrared emission spectroscopy. With the aid of spectral simulation, the CO spectra are rotationally resolved for each vibrational state (v = 1-8). Each of the resulting vibrational and rotational population distributions is characteristic of two Boltzmann profiles with different temperatures, originating from either transition state pathway or OH-roaming to form the same CO + H₂O products. The H₂O roaming co-product is also spectrally simulated to understand the interplay with the CO product in the internal energy partitioning. Accordingly, this work has evaluated the internal energy disposal for the CO and H₂O roaming products; especially the vibrational-state dependence of the roaming signature is reported for the first time. Further, given a 1 μs resolution, the temporal dependence of the CO/CO₂ product ratio at v ≥ 1 rises from 3 to 10 of study, thereby characterizing the effect of conformational memory and well reconciling with the disputed results reported previously between absorption and emission methods.

Stretching of cis-formic acid: warm-up and cool-down as molecular work-out

A new technique to rotationally simplify and Raman-probe conformationally and vibrationally excited small molecules is applied to the cis-trans isomerism of formic acid. It quintuples the previously available gas phase vibrational data base on this excited form of a strongly anharmonic planar molecule despite its limited spectral resolution. The newly determined cis-formic acid fundamentals allow for a balanced vibrational benchmark on both rotamers of formic acid. Assuming the adequacy of vibrational perturbation theory, it reveals weaknesses of standard methods for these systems like B3LYP-D3(BJ)/aVQZ VPT2 or PBE0-D3(BJ)/aVQZ VPT2. The functionals ωB97-XD and M06-2X additionally suffer from severe integration grid size and symmetry dependencies. The vibrational benchmark suggests B2PLYP-D3(BJ)/aVQZ VPT2 and MP2/aVQZ VPT2 as partially competitive and in any case efficient alternatives to computationally demanding coupled cluster vibrational configuration interaction calculations. Whether this is due to fortuitous compensation between electronic structure and vibrational perturbation error remains to be explored.

Roaming Dynamics in the Photodissociation of Formic Acid at 230 nm

Roaming dynamics is observed in the photodissociation of formic acid (HCOOH) at 230 nm by using the slice imaging method. In combination with rotational state selective (2 + 1) resonance-enhanced multiphoton ionization of the CO fragments, the speed distributions of the CO fragments exhibit a low recoil velocity at low rotational levels of J = 9 and 20, while the velocity distributions of CO at high rotational levels of J = 30 and 48 show a relatively large recoil velocity. The experimental results indicate that the roaming of OH radical should be related with the formation of CO + H₂O channel at the present photolysis energy. Unlike the roaming pathways occurring in H₂CO that can be described by loose flat potential, our CO speed distribution analysis suggests the presence of a "tight" flat potential in the roaming dynamics of HCOOH molecules.

Radiation-induced transformations of methanol molecules in low-temperature solids: a matrix isolation study

Radiation-induced transformations of methanol in inert solids at 6 K reveal remarkable matrix effects, and mechanisms and astrochemical implications are discussed.

Radiation-induced transformations of matrix-isolated formic acid: evidence for the HCOOH → HOCO + H channel

HOCO radicals are produced under X-ray irradiation of isolated formic acid at low temperatures, and possible mechanisms and astrochemical implications are discussed.

Automated Exploration of Photolytic Channels of HCOOH: Conformational Memory via Excited-State Roaming

To elucidate the photodissociation mechanism of HCOOH, we systematically explored reaction pathways starting from the first excited singlet state (S1) by using automated reaction path search methods. All critical points, that is, minima, transition states, minimum energy conical intersections, and minima on seam of crossing, for the S0, T1, and S1 potential energy surfaces (PESs) obtained in the present search were optimized at the CASPT2 level. The structure list obtained by the search explained all experimentally reported photolytic channels. A new mechanism for the previously suggested but unexplained conformational memory in the 193 nm photolysis is proposed, which involves two steps: partial dissociation and succeeding roaming of one of H atoms on the S1 PES, followed by intramolecular recombination on the S0 PES after radiationless transition through a conical intersection at a partially dissociated geometry. This is partially similar to the excited-state roaming recently discovered for the NO3 radical.

Transient H2O Infrared Satellite Peaks Produced in UV Irradiated Formic Acid Doped Solid Parahydrogen

We report newly identified satellite features of the R(0) rovibrational transition of all the fundamental modes of HDO and the ν3 mode of H2O measured via FTIR spectroscopy immediately after the 193 nm in situ photolysis of formic acid (HCOOH and DCOOD) in solid parahydrogen. The intensities of these satellite features decay slowly with a time constant of τ = 121(7) min after photolysis, even when the sample is maintained below 2 K. We propose that the van der Waals complex H···H2O (H···HDO) is the carrier of the satellite peaks and that these metastable complexes are produced after the low-temperature tunneling reaction of the OH (OD) photoproduct with the parahydrogen host.

Single photon ionization of hydrogen bonded clusters with a soft x-ray laser: (HCOOH)x and (HCOOH)y(H2O)z

Pure, neutral formic acid (HCOOH)n+1 clusters and mixed (HCOOH)(H2O) clusters are investigated employing time of flight mass spectroscopy and single photon ionization at 26.5 eV using a very compact, capillary discharge, soft x-ray laser. During the ionization process, neutral clusters suffer little fragmentation because almost all excess energy above the vertical ionization energy is taken away by the photoelectron, leaving only a small part of the photon energy deposited into the (HCOOH)n+1+ cluster. The vertical ionization energy minus the adiabatic ionization energy is enough excess energy in the clusters to surmount the proton transfer energy barrier and induce the reaction (HCOOH)n+1+-->(HCOOH)nH+ +HCOO making the protonated (HCOOH)nH+ series dominant in all data obtained. The distribution of pure (HCOOH)nH+ clusters is dependent on experimental conditions. Under certain conditions, a magic number is found at n=5. Metastable dissociation rate constants of (HCOOH)nH+ are measured in the range (0.1-0.8)x10(4) s(-1) for cluster sizes 4<n<9. The rate constants display an odd/even alternating behavior between monomer and dimer loss that can be attributed to the structure of the cluster. When small amounts of water are added to the formic acid, the predominant signals in the mass spectrum are still (HCOOH)nH+ cluster ions. Also observed are the protonated mixed cluster series (HCOOH)n(H2O)mH+ for n=1-8 and m=0-4. A magic number in the cluster series n=5, m=1 is observed. The mechanisms and dynamics of formation of these neutral and ionic clusters are discussed.

Conformationally controlled chemistry: excited-state dynamics dictate ground-state reaction

Ion imaging reveals distinct photodissociation dynamics for propanal cations initially prepared in either the cis or gauche conformation, even though these isomers differ only slightly in energy and face a small interconversion barrier. The product kinetic energy distributions for the hydrogen atom elimination channels are bimodal, and the two peaks are readily assigned to propanoyl cation or hydroxyallyl cation coproducts. Ab initio multiple spawning dynamical calculations suggest that distinct ultrafast dynamics in the excited state deposit each conformer in isolated regions of the ground-state potential energy surface, and, from these distinct regions, conformer interconversion does not effectively compete with dissociation.

Photo-induced cold vapor generation with low molecular weight alcohol, aldehyde, or carboxylic acid for atomic fluorescence spectrometric determination of mercury

With UV irradiation, Hg(2+) in aqueous solution can be converted into Hg(0) cold vapor by low molecular weight alcohols, aldehydes, or carboxylic acids, e.g., methanol, formaldehyde, acetaldehyde, glycol, 1,2-propanediol, glycerol, acetic acid, oxalic acid, or malonic acid. It was found that the presence of nano-TiO(2) more or less improved the efficiency of the photo-induced chemical/cold vapor generation (photo-CVG) with most of the organic reductants. The nano-TiO(2)-enhanced photo-CVG systems can be coupled to various analytical atomic spectrometric techniques for the determination of ultratrace mercury. In this work, we evaluated the application of this method to the atomic fluorescence spectrometric (AFS) determination of mercury in cold vapor mode. Under the optimized experimental conditions, the instrumental limits of detection (based on three times the standard deviation of 11 measurements of a blank solution) were around 0.02-0.04 microg L(-1), with linear dynamic ranges up to 15 microg L(-1). The interference of transition metals and the mechanism of the photo-CVG are briefly discussed. Real sample analysis using the photo-CVG-AFS method revealed that it was promising for water and geological analysis of ultralow levels of mercury.

Solvent Effects on the Photodissociation of Formic Acid: A Theoretical Study

Photodissociation of aqueous formic acid has been investigated with the CASSCF, DFT, and MR-CI methods. Solvent effects are considered as a combination of the hydrogen-bonding interaction from explicit H2O molecules and the effects from the bulk surrounding H2O molecules using the polarizable continuum model. It is found that the hydrogen-bonding effect from the explicit water in the complex is the major factor to influence properties of aqueous formic acid, while the bulk surrounding H2O molecules has a noticeable influence on the structures of the complex. The direct C-O bond fission along the S1 pathway is predicted to be an important channel upon photolysis of aqueous formic acid at 200 nm, which is consistent with experimental observation that aqueous formic acid dissociates predominantly into fragments of HCO and OH. The existence of a dark channel upon photolysis of aqueous formic acid at 200 nm is assigned as fast relaxation from the S1 Franck-Condon geometry to the T1/S1 intersection and subsequent S1-->T1 intersystem crossing process. S1-->S0 internal conversion followed by molecular elimination to CO+H2O is the most probable primary process for formation of carbon monoxide, which was observed with considerable yield upon photolysis of aqueous formic acid at 253.7 nm.

UV Photolysis Products of Propiolic Acid in Noble-Gas Solids

Photolysis (193 nm) of propiolic acid (HCCCOOH) was studied with Fourier transform infrared spectroscopy in noble-gas (Ar, Kr, and Xe) solid matrixes. The photolysis products were assigned using ab initio quantum chemistry calculations. The novel higher-energy conformer of propiolic acid was efficiently formed upon UV irradiation, and it decayed back to the ground-state conformer on a time scale of approximately 10 min by tunneling of the hydrogen atom through the torsional energy barrier. In addition, the photolysis produced a number of matrix-isolated 1:1 molecular complexes such as HCCH...CO2, HCCOH...CO, and H2O...C3O. The HCCH...CO2 complex dominated among the photolysis products, and the computations suggested a parallel geometry of this complex characterized by an interaction energy of -9.6 kJ/mol. The HCCOH...CO complex also formed efficiently, but its concentration was strongly limited by its light-induced decomposition. In this complex, the most probable geometry was found to feature the interaction of carbon monoxide with the OH group via the carbon atom, and the computational interaction energy was determined to be -18.3 kJ/mol. The formation of the strong H2O...C3O complex (interaction energy -21 kJ/mol) was less efficient, which might be due to the inefficiency of the involved radical reaction.

Ultrasound-Promoted Cold Vapor Generation in the Presence of Formic Acid for Determination of Mercury by Atomic Absorption Spectrometry

A new cold vapor technique within the context of green chemistry is described for determination of mercury in liquid samples following high-intensity ultrasonication. Volatile Hg evolved in a sonoreactor without the use of a chemical reducing agent is carried to a quartz cell kept at room temperature for measurement of the atomic absorption. The mechanism involved lies in the reduction of Hg(II) to Hg(0) by reducing gases formed upon sonication and subsequent volatilization of Hg(0) due to the degassing effect caused by the cavitation phenomenon. Addition of a low molecular weight organic acid such as formic acid favors the process, but vapor generation also occurs from Hg solutions in ultrapure water. The detection limit of Hg was 0.1 microg/L, and the repeatability, expressed as relative standard deviation, was 4.4% (peak height). Addition of small amounts of oxidizing substances such as the permanganate or dichromate anions completely suppressed the formation of Hg(0), which confirms the above mechanism. Effect of other factors such as ultrasound irradiation time, ultrasound amplitude, and the presence of concomitants are also investigated. Some complexing anions such as chloride favored the stabilization of Hg(II) in solution, hence causing an interference effect on the ultrasound-assisted reduction/volatilization process.

Different chemical dynamics for different conformers of biological molecules: Photoionization of glycine

Single-photon ionization dynamics of two conformers of glycine is studied by classical trajectory simulations using the semiempirical PM3 potential surface in "on the fly" calculations. Initial conditions for the trajectories are weighted according to the Wigner distribution function computed for the initial vibrational ground state. Vertical ionization in the spirit of the classical Franck-Condon principle is assumed. The dynamics of the two conformers are compared during the first 10 ps. The comparison shows very different dynamical behavior for the two conformers. In particular, the chemical fragmentation pathways differ in part. Also, one of the conformers gives much higher rates of conformational transitions, while the other conformer gives larger chemical fragmentation yields. The example shows significantly different chemical dynamics for two conformers close in energy and separated by a low barrier.

On the nature of bonding in HCOOH...Ar and HCOOH...Kr complexes

The chemical interaction in HCOOH...Ng (Ng=Ar, Kr) complex was analyzed by topological analysis of the electron density based on Atoms-In-Molecules theory. For all computationally stable equilibrium structures of 1:1 HCOOH...Ng complexes, an intermolecular bond path with a "bond" critical point was found and perturbation of formic acid (FA) atomic basins and electron density was observed. The intermolecular interaction between the two complex subunits can be classified, according to its nature, as a closed-shell van der Waals type of interaction. However, one of the computed structures (complex II), pictures a noble gas atom attached linearly to the acidic O-H tail of FA. In this particular case, the electron density at the intermolecular bond critical point was found to resemble a hydrogen-bonded system and thus, may be termed a hydrogen-bond-like interaction. This change in the nature of the interaction is also shown by large perturbations of the FA properties found for this complex structure. The structural and vibrational perturbations are larger than for the other three structures and they increase for the Kr complexes compared to the Ar complex. [Figure: see text]. Electron density analysis of HCOOH...Ng (Ng=Ar,Kr) complex.

Internal Rotation in Propionic Acid: Near-Infrared-Induced Isomerization in Solid Argon

The conformational system of propionic acid (CH3CH2COOH) is studied in solid argon. It is predicted by the ab initio calculations that this molecule has four stable conformers. These four structures are denoted Tt, Tg+/-, Ct, and Cg+/-, and they differ by the arrangement around the C-O and Calpha-C bonds. The ground-state Tt conformer is the only form present at 8 K after deposition of an argon matrix containing propionic acid. For the CH3CH2COOH and CH3CH2COOD isotopologues, narrow-band excitation of the first hydroxyl stretching overtone of the conformational ground state promotes the Calpha-C and C-O internal rotations producing the Tg+/- and Ct conformers, respectively. A subsequent vibrational excitation of the produced Tg+/- form induces its conversion to the Cg+/- conformer by rotation around the C-O bond. In the dark, all of the produced conformers decay to the conformational ground state at different rates. The decay kinetics and its temperature dependence allow the identification of the conformers by IR absorption spectroscopy, which is supported by ab initio calculations of their vibrational spectra. For the CH3CH2COOD isotopologue, the excitation of molecules isolated in different matrix sites results in site-dependent photoisomerization rates for the Calpha-C and C-O internal rotations, which also confirm the identification of the photoproducts.

On the Conformational Memory in the Photodissociation of Formic Acid

The photodissociation of formic acid at 248 and 193 nm was investigated by classical trajectory and RRKM calculations using an interpolated potential energy surface, iteratively constructed using the B3LYP/aug-cc-pVDZ level of calculation. Several sampling schemes in the ground electronic state were employed to explore the possibility of conformational memory in formic acid. The CO/CO2 branching ratios obtained from trajectories initiated at the cis and at the trans conformers are almost identical to each other and in very good accordance with the RRKM results. In addition, when a specific initial excitation that simulates more rigorously the internal conversion process is used, the calculated branching ratio does not vary with respect to those obtained from cis and trans initializations. This result is at odds with the idea of conformational memory in the ground state proposed recently for the interpretation of the experimental results. It was also found that the calculated CO vibrational distributions after dissociation of the parent molecule at 248 nm are in agreement with the experimental available data.