Proton transfer reaction (PTR) is a commonly applied ionization technique for mass spectrometers, in which hydronium ions (H₃O⁺) transfer a proton to analytes with higher proton affinities than the water molecule. This method has most commonly been used to quantify volatile hydrocarbons, but later-generation PTR instruments have been designed for better throughput of less volatile species, allowing detection of more functionalized molecules as well. For example, the recently developed Vocus PTR time-of-flight mass spectrometer (PTR-TOF) has been shown to agree well with an iodide-adduct-based chemical ionization mass spectrometer (CIMS) for products with 3–5 O atoms from oxidation of monoterpenes (C₁₀H₁₆). However, while several different types of CIMS instruments (including those using iodide) detect abundant signals also at “dimeric” species, believed to be primarily ROOR peroxides, no such signals have been observed in the Vocus PTR even though these compounds fulfil the condition of having higher proton affinity than water. More traditional PTR instruments have been limited to volatile molecules as the inlets have not been designed for transmission of easily condensable species. Some newer instruments, like the Vocus PTR, have overcome this limitation but are still not able to detect the full range of functionalized products, suggesting that other limitations need to be considered. One such limitation, well-documented in PTR literature, is the tendency of protonation to lead to fragmentation of some analytes. In this work, we evaluate the potential for PTR to detect dimers and the most oxygenated compounds as these have been shown to be crucial for forming atmospheric aerosol particles. We studied the detection of dimers using a Vocus PTR-TOF in laboratory experiments, as well as through quantum chemical calculations. Only noisy signals of potential dimers were observed during experiments on the ozonolysis of the monoterpene α-pinene, while a few small signals of dimeric compounds were detected during the ozonolysis of cyclohexene. During the latter experiments, we also tested varying the pressures and electric fields in the ionization region of the Vocus PTR-TOF, finding that only small improvements were possible in the relative dimer contributions. Calculations for model ROOR and ROOH systems showed that most of these peroxides should fragment partially following protonation. With the inclusion of additional energy from the ion–molecule collisions driven by the electric fields in the ionization source, computational results suggest substantial or nearly complete fragmentation of dimers. Our study thus suggests that while the improved versions of PTR-based mass spectrometers are very powerful tools for measuring hydrocarbons and their moderately oxidized products, other types of CIMS are likely more suitable for the detection of ROOR and ROOH species.

The oxidation of aromatics contributes significantly to the formation of atmospheric aerosol. Using toluene as an example, we demonstrate the existence of a molecular rearrangement channel in the oxidation mechanism. Based on both flow reactor experiments and quantum chemical calculations, we show that the bicyclic peroxy radicals (BPRs) formed in OH-initiated aromatic oxidation are much less stable than previously thought, and in the case of the toluene derived ipso-BPRs, lead to aerosol-forming low-volatility products with up to 9 oxygen atoms on sub-second timescales. Similar results are predicted for ipso-BPRs formed from many other aromatic compounds. This reaction class is likely a key route for atmospheric aerosol formation, and including the molecular rearrangement of BPRs may be vital for accurate chemical modeling of the atmosphere.

We study the pressure-driven flow of aqueous NaCl in amorphous silica nanotubes using nonequilibrium molecular dynamics simulations featuring both polarizable and non-polarizable molecular models. Different pressures, electrolyte concentrations and pore sizes are examined. Our results indicate a flow that deviates considerably from the predictions of Poiseuille fluid mechanics. Due to preferential adsorption of the different ionic species by surface SiO⁻ or SiOH groups, we find that a significant electric current is generated, but with opposite polarities using polarizable vs. fixed charge models for water and ions, emphasizing the need for careful parameterization in such complex systems. We also examine the influence of partial deprotonation of the silica surface, and we find that much more current is generated in a dehydrogenated nanopore, even though the overall efficiency remains low. These findings indicate that different methods of nanopore preparation, which can produce a range of surface properties, should be examined more closely in the related experimental methods to generate electrokinetic current.

In this paper we study collisions between polyatomic radicals – an important process in fields ranging from biology to combustion. Energy transfer, formation of intermediate complexes and recombination reactions are treated, with applications to peroxy radicals in atmospheric chemistry. Multi-reference perturbation theory, supplemented by coupled-cluster calculations, describes the potential energy surfaces with high accuracy, including the interaction of singlet and triplet spin states during radical recombination. Our multi-reference molecular dynamics (MD) trajectories on methyl peroxy radicals confirm the reaction mechanism postulated in earlier studies. Specifically, they show that if suitable pre-reactive complexes are formed, they will rapidly lead to the formation and subsequent decomposition of tetroxide intermediates. However, generating multi-reference MD trajectories is exceedingly computationally demanding, and we cannot adequately sample the whole conformational space. To answer this challenge, we promote the use of a novel simplified semi-empirical MD methodology. It assumes the collision is governed by two states, a singlet (S₀) and a triplet (T₁) state. The method predicts differences between collisions on S₀ and T₁ surfaces, and qualitatively includes not only pre-reactive complex formation, but also recombination processes such as tetroxide formation. Finally, classical MD simulations using force-fields for non-reactive collisions are employed to generate thousands of collision trajectories, to verify that the semi-empirical method is sampling collisions adequately, and to carry out preliminary investigations of larger systems. For systems with low activation energies, the experimental rate coefficient is surprisingly well reproduced by simply multiplying the gas-kinetic collision rate by the simulated probability for long-lived complex formation.

We propose to analyze molecular dynamics (MD) output via a supervised machine learning (ML) algorithm, the decision tree. The approach aims to identify the predominant geometric features which correlate with trajectories that transition between two arbitrarily defined states. The data-driven algorithm aims to identify these features without the bias of human “chemical intuition”. We demonstrate the method by analyzing the proton exchange reactions in formic acid solvated in small water clusters. The simulations were performed with ab initio MD combined with a method to efficiently sample the rare event, path sampling. Our ML analysis identified relevant geometric variables involved in the proton transfer reaction and how they may change as the number of solvating water molecules changes.