Ethanol is a molecule of fundamental interest in combustion, astrochemistry, and condensed phase as a solvent. It is characterized by two methyl rotors and trans (anti) and gauche conformers, which are known to be very close in energy. Here we show that based on rigorous quantum calculations of the vibrational zero-point state, using a new ab initio potential energy surface (PES), the ground state resembles the trans conformer, but substantial delocalization to the gauche conformer is present. This explains experimental issues about identification and isolation of the two conformers. This "leak" effect is partially quenched when deuterating the OH group, which further demonstrates the need for a quantum mechanical approach. Diffusion Monte Carlo and full-dimensional semiclassical dynamics calculations are employed. The new PES is obtained by means of a Δ-machine learning approach starting from a pre-existing low level density functional theory surface. This surface is brought to the CCSD(T) level of theory using a relatively small number of ab initio CCSD(T) energies. Agreement between the corrected PES and direct ab initio results for standard tests is excellent. One- and two-dimensional discrete variable representation calculations focusing on the trans-gauche torsional motion are also reported, in reasonable agreement with experiment.
The vibrational strong coupling (VSC) between molecular vibrations and cavity photon modes has recently emerged as a promising tool for influencing chemical reactivities. Despite numerous experimental and theoretical efforts, the underlying mechanism of VSC effects remains elusive. In this study, we combine state-of-art quantum cavity vibrational self-consistent field/configuration interaction theory (cav-VSCF/VCI), quasi-classical trajectory method, along with the quantum-chemical CCSD(T)-level machine learning potential, to simulate the hydrogen bond dissociation dynamics of water dimer under VSC. We observe that manipulating the light-matter coupling strength and cavity frequencies can either inhibit or accelerate the dissociation rate. Furthermore, we discover that the cavity surprisingly modifies the vibrational dissociation channels, with a pathway involving both water fragments in their ground vibrational states becoming the major channel, which is a minor one when the water dimer is outside the cavity. We elucidate the mechanisms behind these effects by investigating the critical role of the optical cavity in modifying the intramolecular and intermolecular coupling patterns. While our work focuses on single water dimer system, it provides direct and statistically significant evidence of VSC effects on molecular reaction dynamics.
Tropolone, a 15-atom cyclic molecule, has received much interest both experimentally and theoretically due to its H-transfer tunneling dynamics. An accurate theoretical description is challenging owing to the need to develop a high-level potential energy surface (PES) and then to simulate quantum-mechanical tunneling on this PES in full dimensionality. Here, we tackle both aspects of this challenge and make detailed comparisons with experiments for numerous isotopomers. The PES, of near CCSD(T)-quality, is obtained using a Δ-machine learning approach starting from a pre-existing low-level DFT PES and corrected by a small number of approximate CCSD(T) energies obtained using the fragmentation-based molecular tailoring approach. The resulting PES is benchmarked against DF-FNO-CCSD(T) and CCSD(T)-F12 calculations. Ring-polymer instanton calculations of the splittings, obtained with the Δ-corrected PES are in good agreement with previously reported experiments and a significant improvement over those obtained using the low-level DFT PES. The instanton path includes heavy-atom tunneling effects and cuts the corner, thereby avoiding passing through the conventional saddle-point transition state. This is in contradistinction with typical approaches based on the minimum-energy reaction path. Finally, the subtle changes in the splittings for some of the heavy-atom isotopomers seen experimentally are reproduced and explained.
A recent full-dimensional Δ-Machine learning potential energy surface (PES) for ethanol is employed in semiclassical and vibrational self-consistent field (VSCF) and virtual-state configuration interaction (VCI) calculations, using MULTIMODE, to determine the anharmonic vibrational frequencies of vibration for both the trans and gauche conformers of ethanol. Both semiclassical and VSCF/VCI energies agree well with the experimental data. We find significant mixing between the VSCF basis states due to Fermi resonances between bending and stretching modes. The same effects are also accurately described by the full-dimensional semiclassical calculations. These are the first high-level anharmonic calculations using a PES, in particular a "gold-standard" CCSD(T) one.
Machine-learned potential energy surfaces (PESs) for molecules with more than 10 atoms are typically forced to use lower-level electronic structure methods such as density functional theory (DFT) and second-order Møller-Plesset perturbation theory (MP2). While these are efficient and realistic, they fall short of the accuracy of the "gold standard"coupled-cluster method, especially with respect to reaction and isomerization barriers. We report a major step forward in applying a Δ-machine learning method to the challenging case of acetylacetone, whose MP2 barrier height for H-atom transfer is low by roughly 1.1 kcal/mol relative to the benchmark CCSD(T) barrier of 3.2 kcal/mol. From a database of 2151 local CCSD(T) energies and training with as few as 430 energies, we obtain a new PES with a barrier of 3.5 kcal/mol in agreement with the LCCSD(T) barrier of 3.5 kcal/mol and close to the benchmark value. Tunneling splittings due to H-atom transfer are calculated using this new PES, providing improved estimates over previous ones obtained using an MP2-based PES.
Roaming is a novel mechanism in reaction dynamics. It describes an unusual pathway, which can be quite different from the conventional minimum-energy path, leading to products. While roaming has been reported or suggested in a number of unimolecular reactions, it has been rarely reported for bimolecular reactions. Here, we report a high-level computational study of roaming dynamics in the important bimolecular combustion reaction H + C2H4 → H2 + C2H3, using a new, high-level machine learning-based potential energy surface. In addition to the complex-mediated roaming mechanism, a non-complex forming roaming mechanism is found. It can be described as a direct inelastic collision where the departing H atom roams and then abstracts an H atom. We denoted this as "collision-induced" roaming. These two roaming mechanisms have different angular distributions; however, both produce highly internally excited C2H3. The roaming pathway leads to remarkably different dynamics as compared with the direct abstraction pathway. A clear signature of the roaming mechanism is highly internally excited C2H3, which could be observed experimentally.
In this issue of PNAS, Houston and Kable describe experiments that implicate a new mechanism for chemical reactions in the photodissociation of acetaldehyde. The “roaming” mechanism was discovered in a recent combined experimental/theoretical study of formaldehyde photodissociation. (Formaldehyde is acetaldehyde with the methyl group replaced by hydrogen.) This mechanism is a reaction pathway that bypasses the conventional transition state of the reaction. It has major consequences for experimental observations in the case of formaldehyde photodissociation, and its explication solved a puzzle brought to the attention of researchers in the field by an important paper on formaldehyde photodissociation by van Zee et al. in 1993.
An efficient method is described for generating a fragmented, permutationally invariant polynomial basis to fit electronic energies and, if available, gradients for large molecules. The method presented rests on the fragmentation of a large molecule into any number of fragments while maintaining the permutational invariance and uniqueness of the polynomials. The new approach improves on a previous one reported by Qu and Bowman by avoiding repetition of polynomials in the fitting basis set and speeding up gradient evaluations while keeping the accuracy of the PES. The method is demonstrated for CH3-NH-CO-CH3 (N-methylacetamide) and NH2-CH2-COOH (glycine).
by
Brianna R. Heazlewood;
Meredith J. T. Jordan;
Scott H. Kable;
Talitha M. Selby;
David L. Osborn;
Benjamin C. Shepler;
Bastiaan J. Braams;
Joel M Bowman
Reaction pathways that bypass the conventional saddle-point transition state (TS) are of considerable interest and importance. An example of such a pathway, termed “roaming,” has been described in the photodissociation of H2CO. In a combined experimental and theoretical study, we show that roaming pathways are important in the 308-nm photodissociation of CH3CHO to CH4 + CO. The CH4 product is found to have extreme vibrational excitation, with the vibrational distribution peaked at ≈95% of the total available energy. Quasiclassical trajectory calculations on full-dimensional potential energy surfaces reproduce these results and are used to infer that the major route to CH4 + CO products is via a roaming pathway where a CH3 fragment abstracts an H from HCO. The conventional saddle-point TS pathway to CH4 + CO formation plays only a minor role. H-atom roaming is also observed, but this is also a minor pathway. The dominance of the CH3 roaming mechanism is attributed to the fact that the CH3 + HCO radical asymptote and the TS saddle-point barrier to CH4 + CO are nearly isoenergetic. Roaming dynamics are therefore not restricted to small molecules such as H2CO, nor are they limited to H atoms being the roaming fragment. The observed dominance of the roaming mechanism over the conventional TS mechanism presents a significant challenge to current reaction rate theory.
by
Dr. Raffael Schwan;
Dr. Chen Qu;
Dr. Devendra Mani;
Nitish Pal;
Dr. Lex van der Meer;
Dr. Britta Redlich;
Prof. Claude Leforestier;
Joel Bowman;
Gerhard Schwaab;
Prof. Martina Havenith
Using the helium nanodroplet isolation setup at the ultrabright free-electron laser source FELIX in Nijmegen (BoHeNDI@FELIX), the intermolecular modes of water dimer in the frequency region from 70 to 550 cm−1 were recorded. Observed bands were assigned to donor torsion, acceptor wag, acceptor twist, intermolecular stretch, donor torsion overtone, and in-plane and out-of-plane librational modes. This experimental data set provides a sensitive test for state-of-the-art water potentials and dipole moment surfaces. Theoretical calculations of the IR spectrum are presented using high-level quantum and approximate quasiclassical molecular dynamics approaches. These calculations use the full-dimensional ab initio WHHB potential and dipole moment surfaces. Based on the experimental data, a considerable increase of the acceptor switch and a bifurcation tunneling splitting in the librational mode is deduced, which is a consequence of the effective decrease in the tunneling barrier.