Anharmonic Rovibrational Partition Functions at High Temperatures: Tests of Reduced-Dimensional Models for Systems with up to Three Fluxional Modes

Methanetriol─Formation of an Impossible Molecule

Orthocarboxylic acids─organic molecules carrying three hydroxyl groups at the same carbon atom─have been distinguished as vital reactive intermediates by the atmospheric science and physical (organic) chemistry communities as transients in the atmospheric aerosol cycle. Predicted short lifetimes and their tendency to dehydrate to a carboxylic acid, free orthocarboxylic acids, signify one of the most elusive classes of organic reactive intermediates, with even the simplest representative methanetriol (CH(OH)3)─historically known as orthoformic acid─not previously been detected experimentally. Here, we report the first synthesis of the previously elusive methanetriol molecule in low-temperature mixed methanol (CH3OH) and molecular oxygen (O2) ices subjected to energetic irradiation. Supported by electronic structure calculations, methanetriol was identified in the gas phase upon sublimation via isomer-selective photoionization reflectron time-of-flight mass spectrometry combined with isotopic substitution studies and the detection of photoionization fragments. The first synthesis and detection of methanetriol (CH(OH)3) reveals its gas-phase stability as supported by a significant barrier hindering unimolecular decomposition. These findings progress our fundamental understanding of the chemistry and chemical bonding of methanetriol, hydroxyperoxymethane (CH3OOOH), and hydroxyperoxymethanol (CH2(OH)OOH), which are all prototype molecules in the oxidation chemistry of the atmosphere.

PESPIP: Software to fit complex molecular and many-body potential energy surfaces with permutationally invariant polynomials

We wish to describe a potential energy surface by using a basis of permutationally invariant polynomials whose coefficients will be determined by numerical regression so as to smoothly fit a dataset of electronic energies as well as, perhaps, gradients. The polynomials will be powers of transformed internuclear distances, usually either Morse variables, exp(-r_i,j/λ), where λ is a constant range hyperparameter, or reciprocals of the distances, 1/r_i,j. The question we address is how to create the most efficient basis, including (a) which polynomials to keep or discard, (b) how many polynomials will be needed, (c) how to make sure the polynomials correctly reproduce the zero interaction at a large distance, (d) how to ensure special symmetries, and (e) how to calculate gradients efficiently. This article discusses how these questions can be answered by using a set of programs to choose and manipulate the polynomials as well as to write efficient Fortran programs for the calculation of energies and gradients. A user-friendly interface for access to monomial symmetrization approach results is also described. The software for these programs is now publicly available.

Predicting third-body collision efficiencies for water and other polyatomic baths

Low-pressure-limit microcanonical (collisional activation) and thermal rate constants are predicted using a combination of automated ab initio potential energy surface construction, classical trajectories, transition state theory, and a detailed kinetic...

Spiers Memorial Lecture: theory of unimolecular reactions

One hundred years ago, at an earlier Faraday Discussion meeting, Lindemann presented a mechanism that provides the foundation for contemplating the pressure dependence of unimolecular reactions. Since that time, our...

Coupled Anharmonic Thermochemistry from Stratified Monte Carlo Integration

This study presents configuration integral Monte Carlo integration (CIMCI), a new semiclassical method for handling fully coupled anharmonicity in gas-phase thermodynamics that promises to be black boxable, to be applicable to all kinds of anharmonicity, and to scale better at higher dimensionality than other methods for handling gas-phase molecular anharmonicity. The method does so using automatically and recursively stratified, simultaneous Monte Carlo (MC) integration of multiple functions, following a modified version of the standard MISER scheme that converges at a rate of about the square of naïve MC integration. For the small systems analyzed by this study where proper reference data is available (H₂O and H₂O₂), the method's anharmonic entropy corrections match reference data better than those of other black box anharmonic methods, e.g., vibrational perturbation theory (VPT2) and the McClurg hindered rotor model used with automatic detection of rotors; for H₂O₂ and NH₂OH, the method is also in general agreement with one-dimensional hindered rotor treatments at low temperatures. This holds even when sampling with CIMCI is done with primitive force fields, e.g., UFF, while the competing methods are used with proper, comprehensive potentials, e.g., the M06-2X metahybrid density-functional theory (DFT) functional.

Parsimonious Potential Energy Surface Expansions Using Dictionary Learning with Multipass Greedy Selection

Potential energy surfaces fit with basis set expansions have been shown to provide accurate representations of electronic energies and have enabled a variety of high-accuracy dynamics, kinetics, and spectroscopy applications. The number of terms in these expansions scales poorly with system size, a drawback that challenges their use for systems with more than ∼10 atoms. A solution is presented here using dictionary learning. Subsets of the full set of conventional basis functions are optimized using a newly developed multipass greedy regression method inspired by forward and backward selection methods from the statistics, signal processing, and machine learning literatures. The optimized representations have accuracies comparable to the full set but are 1 or more orders of magnitude smaller, and notably, the number of terms in the optimized multipass greedy expansions scales approximately linearly with the number of atoms.

Permutationally Invariant Polynomial Expansions with Unrestricted Complexity

A general strategy is presented for constructing and validating permutationally invariant polynomial (PIP) expansions for chemical systems of any stoichiometry. Demonstrations are made for three categories of gas-phase dynamics and kinetics: collisional energy-transfer trajectories for predicting pressure-dependent kinetics, three-body collisions for describing transient van der Waals adducts relevant to atmospheric chemistry, and nonthermal reactivity via quasiclassical trajectories. In total, 30 systems are considered with up to 15 atoms and 39 degrees of freedom. Permutational invariance is enforced in PIP expansions with as many as 13 million terms and 13 permutationally distinct atom types by taking advantage of petascale computational resources. The quality of the PIP expansions is demonstrated through the systematic convergence of in-sample and out-of-sample errors with respect to both the number of training data and the order of the expansion, and these errors are shown to predict errors in the dynamics for both reactive and nonreactive applications. The parallelized code distributed as part of this work enables the automation of PIP generation for complex systems with multiple channels and flexible user-defined symmetry constraints and for automatically removing unphysical unconnected terms from the basis set expansions, all of which are required for simulating complex reactive systems.

Extreme Low-Temperature Combustion Chemistry: Ozone-Initiated Oxidation of Methyl Hexanoate

The accelerating chemical effect of ozone addition on the oxidation chemistry of methyl hexanoate [CH₃(CH₂)₄C(═O)OCH₃] was investigated over a temperature range from 460 to 940 K. Using an externally heated jet-stirred reactor at p = 700 Torr (residence time τ = 1.3 s, stoichiometry φ = 0.5, 80% argon dilution), we explored the relevant chemical pathways by employing molecular-beam mass spectrometry with electron and single-photon ionization to trace the temperature dependencies of key intermediates, including many hydroperoxides. In the absence of ozone, reactivity is observed in the so-called low-temperature chemistry (LTC) regime between 550 and 700 K, which is governed by hydroperoxides formed from sequential O₂ addition and isomerization reactions. At temperatures above 700 K, we observed the negative temperature coefficient (NTC) regime, in which the reactivity decreases with increasing temperatures, until near 800 K, where the reactivity increases again. Upon addition of ozone (1000 ppm), the overall reactivity of the system is dramatically changed due to the time scale of ozone decomposition in comparison to fuel oxidation time scales of the mixtures at different temperatures. While the LTC regime seems to be only slightly affected by the addition of ozone with respect to the identity and quantity of the observed intermediates, we observed an increased reactivity in the intermediate NTC temperature range. Furthermore, we observed experimental evidence for an additional oxidation regime in the range near 500 K, herein referred to as the extreme low-temperature chemistry (ELTC) regime. Experimental evidence and theoretical rate constant calculations indicate that this ELTC regime is likely to be initiated by H abstraction from methyl hexanoate via O atoms, which originate from thermal O₃ decomposition. The theoretical calculations show that the rate constants for methyl ester initiation via abstraction by O atoms increase dramatically with the size of the methyl ester, suggesting that ELTC is likely not important for the smaller methyl esters. Experimental evidence is provided indicating that, similar to the LTC regime, the chemistry in the ELTC regime is dominated by hydroperoxide chemistry. However, mass spectra recorded at various reactor temperatures and at different photon energies provide experimental evidence of some differences in chemical species between the ELTC and the LTC temperature ranges.

Efficient Generation of Permutationally Invariant Potential Energy Surfaces for Large Molecules

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 CH₃–NH–CO–CH₃ (N-methylacetamide) and NH₂–CH₂–COOH (glycine).