Photochemical reactions of the low-lying excited states of formaldehyde: T1/S0 intersystem crossings, characteristics of the S1 and T1 potential energy surfaces, and a global T1 potential energy surface

Theoretical studies on triplet-state driven dissociation of formaldehyde by quasi-classical molecular dynamics simulation on machine-learning potential energy surface

The H-atom dissociation of formaldehyde on the lowest triplet state (T1) is studied by quasi-classical molecular dynamic simulations on the high-dimensional machine-learning potential energy surface (PES) model. An atomic-energy based deep-learning neural network (NN) is used to represent the PES function, and the weighted atom-centered symmetry functions are employed as inputs of the NN model to satisfy the translational, rotational, and permutational symmetries, and to capture the geometry features of each atom and its individual chemical environment. Several standard technical tricks are used in the construction of NN-PES, which includes the application of clustering algorithm in the formation of the training dataset, the examination of the reliability of the NN-PES model by different fitted NN models, and the detection of the out-of-confidence region by the confidence interval of the training dataset. The accuracy of the full-dimensional NN-PES model is examined by two benchmark calculations with respect to ab initio data. Both the NN and electronic-structure calculations give a similar H-atom dissociation reaction pathway on the T1 state in the intrinsic reaction coordinate analysis. The small-scaled trial dynamics simulations based on NN-PES and ab initio PES give highly consistent results. After confirming the accuracy of the NN-PES, a large number of trajectories are calculated in the quasi-classical dynamics, which allows us to get a better understanding of the T1-driven H-atom dissociation dynamics efficiently. Particularly, the dynamics simulations from different initial conditions can be easily simulated with a rather low computational cost. The influence of the mode-specific vibrational excitations on the H-atom dissociation dynamics driven by the T1 state is explored. The results show that the vibrational excitations on symmetric C-H stretching, asymmetric C-H stretching, and C=O stretching motions always enhance the H-atom dissociation probability obviously.

High-fidelity first principles nonadiabaticity: diabatization, analytic representation of global diabatic potential energy matrices, and quantum dynamics

Nonadiabatic dynamics, which goes beyond the Born-Oppenheimer approximation, has increasingly been shown to play an important role in chemical processes, particularly those involving electronically excited states. Understanding multistate dynamics requires rigorous quantum characterization of both electronic and nuclear motion. However, such first principles treatments of multi-dimensional systems have so far been rather limited due to the lack of accurate coupled potential energy surfaces and difficulties associated with quantum dynamics. In this Perspective, we review recent advances in developing high-fidelity analytical diabatic potential energy matrices for quantum dynamical investigations of polyatomic uni- and bi-molecular nonadiabatic processes, by machine learning of high-level ab initio data. Special attention is paid to methods of diabatization, high fidelity construction of multi-state coupled potential energy surfaces and property surfaces, as well as quantum mechanical characterization of nonadiabatic nuclear dynamics. To illustrate the tremendous progress made by these new developments, several examples are discussed, in which direct comparison with quantum state resolved measurements led to either confirmation of the observation or sometimes reinterpretation of the experimental data. The insights gained in these prototypical systems greatly advance our understanding of nonadiabatic dynamics in chemical systems.

Carbon Atom Reactivity with Amorphous Solid Water: H2O-Catalyzed Formation of H2CO

We report new computational and experimental evidence of an efficient and astrochemically relevant formation route to formaldehyde (H₂CO). This simplest carbonylic compound is central to the formation of complex organics in cold interstellar clouds and is generally regarded to be formed by the hydrogenation of solid-state carbon monoxide. We demonstrate H₂CO formation via the reaction of carbon atoms with amorphous solid water. Crucial to our proposed mechanism is a concerted proton transfer catalyzed by the water hydrogen bonding network. Consequently, the reactions ³C + H₂O → ³HCOH and ¹HCOH → ¹H₂CO can take place with low or without barriers, contrary to the high-barrier traditional internal hydrogen migration. These low barriers (or the absence thereof) explain the very small kinetic isotope effect in our experiments when comparing the formation of H₂CO to D₂CO. Our results reconcile the disagreement found in the literature on the reaction route C + H₂O → H₂CO.

Enabling a Unified Description of Both Internal Conversion and Intersystem Crossing in Formaldehyde: A Global Coupled Quasi-Diabatic Hamiltonian for Its S0, S1, and T1 States

In our recent work, a diabatic Hamiltonian that couples the S₀ and S₁ states of formaldehyde was constructed using a robust fitting-and-diabatizing procedure with artificial neural networks, which is capable of representing adiabatic energies, energy gradients, and derivative couplings over a wide range of geometries including seams of conical intersection. In this work, based on the diabatization of S₀ and S₁, the spin-orbit couplings between singlet states (S₀, S₁) and triplet state T₁ are also determined in the same diabatic representation. The diabatized spin-orbit couplings are then fit with a symmetrized neural-network functional form. The ab initio spin-orbit couplings are well reproduced in large configuration space. Together with the neural-network-based potential energy surface for T₁, the full quasi-diabatic Hamiltonian for the S₀, S₁, and T₁ states is completed, enabling a unified description of both internal conversion and intersystem crossing in formaldehyde. The vibrational levels on the three adiabatic states are found to be in good agreement with known experimental band origins.

Roaming Reactions and Dynamics in the van der Waals Region

Roaming reactions were first clearly identified in photodissociation of formaldehyde 15 years ago, and roaming dynamics are now recognized as a universal aspect of chemical reactivity. These reactions typically involve frustrated near-dissociation of a quasibound system to radical fragments, followed by reorientation at long range and intramolecular abstraction. The consequences can be unexpected formation of molecular products, depletion of the radical pool in chemical systems, and formation of products with unusual internal state distributions. In this review, I examine some current aspects of roaming reactions with an emphasis on experimental results, focusing on possible quantum effects in roaming and roaming dynamics in bimolecular systems. These considerations lead to a more inclusive definition of roaming reactions as those for which key dynamics take place at long range.

Accurate Neural Network Representation of the Ab Initio Determined Spin-Orbit Interaction in the Diabatic Representation Including the Effects of Conical Intersections

A method for fitting ab initio determined spin-orbit coupling interactions, in the Breit-Pauli approximation, based on quasidiabatic representations using neural network fits is reported. The algorithm generalizes our recently reported neural network approach for representing the dipole interaction. The S₀, S₁, and T₁ states of formaldehyde are used as an example. First, the two singlet states S₀ and S₁ are diabatized with a modified Boys Localization diabatization method. Second, the spin-orbit coupling between singlet and triplet states is transformed to the diabatic representation. This removes the discontinuities in the adiabatic representation. The diabatized spin-orbit couplings are then fit with smooth neural network functions. The analytic representation of spin-orbit coupling interactions in a diabatic basis by neural networks will make accurate full-dimensional quantum dynamical treatment of both internal conversion and intersystem crossing possible, which will help us to gain better understanding of both processes.

Chemical Kinetics Approves the Occurrence of C (3Pj) Reaction with H2O

Although both atomic carbon and water are omnipresent in human life, there is a debate about the possibility of carbon reaction with water. Some low-temperature spectroscopic investigations have rejected the reaction, whereas some room-temperature experiments and theoretical studies have accepted the possibility of the reaction by reporting rate coefficients ranging from 10⁵ to 10⁹ L mol^-1 s^-1. This study provides new lines of evidence about the reaction through exploration of the reaction mechanism using the CCSD(T) method and solving the corresponding master equation by following two main approaches. According to the results, the rate coefficient of the reaction is significantly influenced by the tunneling and hindered rotation effects, in addition to the selected total angular momentum (J). Furthermore, the total rate coefficient of the reaction increases dramatically (from 10⁷ to 10¹¹ L mol^-1 s^-1) with the rise of temperature from 100 to 4000 K, while the total rate coefficient is insensitive to pressure (0.1-10 atm). Despite some differences between the results of the two approaches, the rate coefficients of both methods are consistent with the previously reported rate coefficients. Also, in agreement with the previous studies, the major products are ²HOC + ²H and ²HCO + ²H. In general, the findings approve the occurrence of the title reaction and indicate that the mentioned conflict is due to the sensitivity of the reaction to the investigated temperature and J level. The sensitivity does not permit low-temperature spectroscopic studies to detect any products and varies the measured and calculated rate coefficients.

Photoinduced C–H bond fission in prototypical organic molecules and radicals

We survey and assess current knowledge regarding the primary photochemistry of hydrocarbon molecules and radicals.

Kinetics and dynamics of the C(3P) + H2O reaction on a full-dimensional accurate triplet state potential energy surface

A full-dimensional accurate PES for the C(³P) + H₂O reaction is developed using the PIP-NN method.

Exploring Multiple Potential Energy Surfaces: Photochemistry of Small Carbonyl Compounds

In theoretical studies of chemical reactions involving multiple potential energy surfaces (PESs) such as photochemical reactions, seams of intersection among the PESs often complicate the analysis. In this paper, we review our recipe for exploring multiple PESs by using an automated reaction path search method which has previously been applied to single PESs. Although any such methods for single PESs can be employed in the recipe, the global reaction route mapping (GRRM) method was employed in this study. By combining GRRM with the proposed recipe, all critical regions, that is, transition states, conical intersections, intersection seams, and local minima, associated with multiple PESs, can be explored automatically. As illustrative examples, applications to photochemistry of formaldehyde and acetone are described. In these examples as well as in recent applications to other systems, the present approach led to discovery of many unexpected nonadiabatic pathways, by which some complicated experimental data have been explained very clearly.

UV photolysis of α-cyclohexanedione in the gas phase

Ultraviolet absorption spectrum of α-cyclohexanedione (α-CHD) vapor in the wavelength range of 220-320 nm has been recorded in a 1 m long path gas cell at room temperature. With the aid of theoretical calculation, the band has been assigned to the S(2) ← S(0) transition of largely ππ* type. The absorption cross section at the band maximum (∼258 nm) is nearly 3 orders of magnitude larger compared to that for the S(2) ← S(0) transition of a linear α-diketo prototype, 2,3-pentanedione. The photolysis was performed by exciting the sample vapor near this band maximum, using the 253.7 nm line of a mercury vapor lamp, and the products were analyzed by mass spectrometry as well as by infrared spectroscopy. The identified products are cyclopentanone, carbon monoxide, ketene, ethylene, and 4-pentenal. Geometry optimization at the CIS/6-311++G** level predicts that the carbonyl group is pyramidally distorted in the excited S(1) and S(2) states, but the α-CHD ring does not show dissociative character. Potential energy curves with respect to a ring rupture coordinate (C-C bond between two carbonyl groups) for S(0), S(1), S(2), T(1), T(2), and T(3) states have been generated by partially optimizing the ground state geometry at DFT/B3LYP/6-311++G** level and calculating the vertical transition energies to the excited states by TDDFT method. Our analysis reveals that the reactions can take place at higher vibrational levels of S(0) as well as T(1) states.