Gas-Phase Oxidation of Atmospherically Relevant Unsaturated Hydrocarbons by Acyl Peroxy Radicals.

This study assesses the atmospheric impact of reactions between unsaturated hydrocarbons such as isoprene and monoterpenes and peroxy radicals containing various functional groups. We find that reactions between alkenes and acyl peroxy radicals have reaction rates high enough to be feasible in the atmosphere and lead to high molar mass accretion products. Moreover, the reaction between unsaturated hydrocarbons and acyl peroxy radicals leads to an alkyl radical, to which molecular oxygen rapidly adds. This finding is confirmed by both theoretical calculations and experiments. The formed perester peroxy radical may either undergo further H-shift reactions or react bimolecularly. The multifunctional oxygenated compounds formed through acyl peroxy radical + alkene reactions are potentially important contributors to particle formation and growth. Thus, acyl peroxy radical-initiated oxidation chemistry may need to be included in atmospheric models.

Theoretical analysis of the OH-initiated atmospheric oxidation reactions of imidazole.

Imidazoles are present in Earth's atmosphere in both the gas-phase and in aerosol particles, and have been implicated in the formation of brown carbon aerosols. The gas-phase oxidation of imidazole (C3N2H4) by hydroxyl radicals has been shown to be preferentially initiated via OH-addition to position C5, producing the 5-hydroxyimidazolyl radical adduct. However, the fate of this adduct upon reaction with O2 in the atmospheric gas-phase is currently unknown. We employed an automated approach to investigate the reaction mechanism and kinetics of imidazole's OH-initiated gas-phase oxidation, in the presence of O2 and NOx. The explored mechanism included reactions available to first-generation RO2 radicals, as well as alkoxyl radicals produced from RO2 + NO reactions. Product distributions were obtained by assembling and solving a master equation, under conditions relevant to the Earth's atmosphere. Our calculations show a complex, branched reaction mechanism, which nevertheless converges to yield two major closed-shell products: 4H-imidazol-4-ol (4H-4ol) and N,N'-diformylformamidine (FMF). At 298 K and 1 atm, we estimate the yields of 4H-4ol and FMF from imidazole oxidation initiated via OH-addition to position C5 to be 34 : 66, 12 : 85 and 2 : 95 under 10 ppt, 100 ppt and 1 ppb of NO respectively. This work also revealed O2-migration pathways between the α-N-imino peroxyl radical isomers. This reaction channel is fast for the first-generation RO2 radicals, and may be important during the atmospheric oxidation of other unsaturated organic nitrogen compounds as well.

Effect of an Electric Field on the Structure and Stability of Atmospheric Clusters.

We study the influence of an applied electric field on the structure and stability of some common bimolecular clusters that are found in the atmosphere. These clusters play an important role in new particle formation (NPF). For low values of the electric field (i.e., |E| ≤ 0.01 V Å-1), we demonstrate that the field response of the clusters can be predicted from simply calculating the dipole moment of the cluster and the dipole moments of the constituent molecules and that the influence on the association energy of the cluster is minimal (i.e., <0.5 kcal mol-1). For higher field strengths |E| > 0.2 V Å-1, there can be more dramatic effects on both structure and energetics, as the induced dipole, charge transfer, and geometric distortion play a larger role. Although such large fields are not very relevant in the atmosphere, they do exist in some situations of experimental interest, such as near interfaces and in intense laser fields.

Multireference and Coupled-Cluster Study of Dimethyltetroxide (MeO4Me) Formation and Decomposition.

Peroxyl radicals (RO2) are important intermediates in the atmospheric oxidation processes. The RO2 can react with other RO2 to form reactive intermediates known as tetroxides, RO4R. The reaction mechanisms of RO4R formation and its various decomposition channels have been investigated in multiple computational studies, but previous approaches have not been able to provide a unified methodology that is able to connect all relevant reactions together. An apparent difficulty in modeling the RO4R formation and its decomposition is the involvement of open-shell singlet electronic states along the reaction pathway. Modeling such electronic states requires multireference (MR) methods, which we use in the present study to investigate in detail a model reaction of MeO2 + MeO2 → MeO4Me, and its decomposition, MeO4Me → MeO + O2 + MeO, as well as the two-body product complexes MeO···O2 + MeO and MeO···MeO + O2. The used MR methods are benchmarked against relative energies of MeO2 + MeO2, MeO4Me, and MeO + MeO + O2, calculated with CCSD(T)/CBS, W2X, and W3X-L methods. We found that the calculated relative energies of the overall MeO2 + MeO2 → MeO4Me → MeO + O2 + MeO reaction are very sensitive to the chosen MR method and that the CASPT2(22e,14o)-IPEA method is able to reproduce the relative energies obtained by the various coupled-cluster methods. Furthermore, CASPT2(22e,14o)-IPEA and W3X-L results show that the MeO···O2 product complex + MeO is lower in energy than the MeO···MeO complex + O2. The formation of MeO4Me is exothermic, and its decomposition is endothermic, relative to the tetroxide. Both the formation and decomposition reactions appear to follow pathways with no saddle points. According to potential energy surface scans of the decomposition pathway, the concerted cleavage of both MeO···O bonds in MeO4Me is energetically preferred over the corresponding sequential decomposition.

Large Gas-Phase Source of Esters and Other Accretion Products in the Atmosphere.

Dimeric accretion products have been observed both in atmospheric aerosol particles and in the gas phase. With their low volatilities, they are key contributors to the formation of new aerosol particles, acting as seeds for more volatile organic vapors to partition onto. Many particle-phase accretion products have been identified as esters. Various gas- and particle-phase formation pathways have been suggested for them, yet evidence remains inconclusive. In contrast, peroxide accretion products have been shown to form via gas-phase peroxy radical (RO2) cross reactions. Here, we show that these reactions can also be a major source of esters and other types of accretion products. We studied α-pinene ozonolysis using state-of-the-art chemical ionization mass spectrometry together with different isotopic labeling approaches and quantum chemical calculations, finding strong evidence for fast radical isomerization before accretion. Specifically, this isomerization seems to happen within the intermediate complex of two alkoxy (RO) radicals, which generally determines the branching of all RO2-RO2 reactions. Accretion products are formed when the radicals in the complex recombine. We found that RO with suitable structures can undergo extremely rapid C-C ß scissions before recombination, often resulting in ester products. We also found evidence of this previously overlooked RO2-RO2 reaction pathway forming alkyl accretion products and speculate that some earlier peroxide identifications may in fact be hemiacetals or ethers. Our findings help answer several outstanding questions on the sources of accretion products in organic aerosol and bridge our knowledge of the gas phase formation and particle phase detection of accretion products. As esters are inherently more stable than peroxides, this also impacts their further reactivity in the aerosol.

Influence of plasmons on the luminescence properties of solvatochromic merocyanine dyes with different solvatochromism.

The effect of localized surface plasmon resonance (LSPR) of a system consisting of a highly dipolar merocyanine dye and a silver nanoparticle (NP) was studied experimentally and theoretically. A theoretical model for estimating the fluorescence quantum yield (φfl) using quantum chemical calculations of intramolecular and intermolecular electronic transition rate constants was developed. Calculations show that the main deactivation channels of the lowest excited singlet state of the studied merocyanines are internal conversion (kIC(S1 → S0)) and fluorescence (kr(S1 → S0)). The intersystem-crossing transition has a low probability due to the large energy difference between the singlet and triplet levels. In the presence of plasmonic NPs, the fluorescence quantum yield is increased by a factor of two according to both experiment and computations. The calculated values of φfl, when considering changes in kr(S1 → S0) and the energy-transfer rate constant (ktransfer) from the dye to the NP was also twice as large at distances of 6-8 nm between the NP and the dye molecule. We also found that the LSPR effect can be increased or decreased depending on the value of the dielectric constant (εm) of the environment.

An Experimental and Master Equation Investigation of Kinetics of the CH2OO + RCN Reactions (R = H, CH3, C2H5) and Their Atmospheric Relevance.

We have performed direct kinetic measurements of the CH2OO + RCN reactions (R = H, CH3, C2H5) in the temperature range 233-360 K and pressure range 10-250 Torr using time-resolved UV-absorption spectroscopy. We have utilized a new photolytic precursor, chloroiodomethane (CH2ICl), whose photolysis at 193 nm in the presence of O2 produces CH2OO. Observed bimolecular rate coefficients for CH2OO + HCN, CH2OO + CH3CN, and CH2OO + C2H5CN reactions at 296 K are (2.22 ± 0.65) × 10-14 cm3 molecule-1 s-1, (1.02 ± 0.10) × 10-14 cm3 molecule-1 s-1, and (2.55 ± 0.13) × 10-14 cm3 molecule-1 s-1, respectively, suggesting that reaction with CH2OO is a potential atmospheric degradation pathway for nitriles. All the reactions have negligible temperature and pressure dependence in the studied regions. Quantum chemical calculations (ωB97X-D/aug-cc-pVTZ optimization with CCSD(T)-F12a/VDZ-F12 electronic energy correction) of the CH2OO + RCN reactions indicate that the barrierless lowest-energy reaction path leads to a ring closure, resulting in the formation of a 1,2,4-dioxazole compound. Master equation modeling results suggest that following the ring closure, chemical activation in the case of CH2OO + HCN and CH2OO + CH3CN reactions leads to a rapid decomposition of 1,2,4-dioxazole into a CH2O + RNCO pair, or by a rearrangement, into a formyl amide (RC(O)NHC(O)H), followed by decomposition into CO and an imidic acid (RC(NH)OH). The 1,2,4-dioxazole, the CH2O + RNCO pair, and the CO + RC(NH)OH pair are atmospherically significant end products to varying degrees.

Computational Investigation of Substituent Effects on the Alcohol + Carbonyl Channel of Peroxy Radical Self- and Cross-Reactions.

Organic peroxy radicals (RO2) are key intermediates in atmospheric chemistry and can undergo a large variety of both uni- and bimolecular reactions. One of the least understood reaction classes of RO2 are their self- and cross-reactions: RO2 + R'O2. In our previous work, we have investigated how RO2 + R'O2 reactions can lead to the formation of ROOR' accretion products through intersystem crossings and subsequent recombination of a triplet intermediate complex 3(RO···OR'). Accretion products can potentially have very low saturation vapor pressures, and may therefore participate in the formation of aerosol particles. In this work, we investigate the competing H-shift channel, which leads to the formation of more volatile carbonyl and alcohol products. This is one of the main, and sometimes the dominant, RO2 + R'O2 reaction channels for small RO2. We investigate how substituents (R and R' groups) affect the H-shift barriers and rates for a set of 3(RO···OR') complexes. The variation in barrier heights and rates is found to be surprisingly small, and most computed H-shift rates are fast: around 108-109 s-1. We find that the barrier height is affected by three competing factors: (1) the weakening of the breaking C-H bond due to interactions with adjacent functional groups; (2) the overall binding energy of the 3(RO···OR'), which tends to increase the barrier height; and (3) the thermodynamic stability of the reaction products. We also calculated intersystem crossing rate coefficients (ISC) for the same systems and found that most of them were of the same order of magnitude as the H-shift rates. This suggests that both studied channels are competitive for small and medium-sized RO2. However, for complex enough R or R' groups, the binding energy effect may render the H-shift channel uncompetitive with intersystem crossings (and thus ROOR' formation), as the rate of the latter, while variable, seems to be largely independent of system size. This may help explain the experimental observation that accretion product formation becomes highly effective for large and multifunctional RO2.

Pathways to Highly Oxidized Products in the Δ3-Carene + OH System.

Oxidation of the monoterpene Δ3-carene (C10H16) is a potentially important and understudied source of atmospheric secondary organic aerosol (SOA). We present chamber-based measurements of speciated gas and particle phases during photochemical oxidation of Δ3-carene. We find evidence of highly oxidized organic molecules (HOMs) in the gas phase and relatively low-volatility SOA dominated by C7-C10 species. We then use computational methods to develop the first stages of a Δ3-carene photochemical oxidation mechanism and explain some of our measured compositions. We find that alkoxy bond scission of the cyclohexyl ring likely leads to efficient HOM formation, in line with previous studies. We also find a surprising role for the abstraction of primary hydrogens from methyl groups, which has been calculated to be rapid in the α-pinene system, and suggest more research is required to determine if this is more general to other systems and a feature of autoxidation. This work develops a more comprehensive view of Δ3-carene photochemical oxidation products via measurements and lays out a suggested mechanism of oxidation via computationally derived rate coefficients.

Critical Role of Iodous Acid in Neutral Iodine Oxoacid Nucleation.

Nucleation of neutral iodine particles has recently been found to involve both iodic acid (HIO3) and iodous acid (HIO2). However, the precise role of HIO2 in iodine oxoacid nucleation remains unclear. Herein, we probe such a role by investigating the cluster formation mechanisms and kinetics of (HIO3)m(HIO2)n (m = 0-4, n = 0-4) clusters with quantum chemical calculations and atmospheric cluster dynamics modeling. When compared with HIO3, we find that HIO2 binds more strongly with HIO3 and also more strongly with HIO2. After accounting for ambient vapor concentrations, the fastest nucleation rate is predicted for mixed HIO3-HIO2 clusters rather than for pure HIO3 or HIO2 ones. Our calculations reveal that the strong binding results from HIO2 exhibiting a base behavior (accepting a proton from HIO3) and forming stronger halogen bonds. Moreover, the binding energies of (HIO3)m(HIO2)n clusters show a far more tolerant choice of growth paths when compared with the strict stoichiometry required for sulfuric acid-base nucleation. Our predicted cluster formation rates and dimer concentrations are acceptably consistent with those measured by the Cosmic Leaving Outdoor Droplets (CLOUD) experiment. This study suggests that HIO2 could facilitate the nucleation of other acids beyond HIO3 in regions where base vapors such as ammonia or amines are scarce.

Energy transfer, pre-reactive complex formation and recombination reactions during the collision of peroxy radicals.

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 (S0) and a triplet (T1) state. The method predicts differences between collisions on S0 and T1 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.

Magnetically induced ring currents in metallocenothiaporphyrins.

The magnetically induced current-density susceptibility tensor (CDT) of the lowest singlet and triplet states of the metallocenothiaporphyrins, where the metal is V, Cr, Mn, Fe, Co, Ni, Mo, Tc, Ru, or Rh, have been studied with the gauge-including magnetically induced currents (GIMIC) method. The compounds containing V, Mn, Co, Tc or Rh were studied as cations because the neutral molecules have an odd number of electrons. The calculations show that the aromatic nature of most of the studied molecules follows the Hückel and Baird rules of aromaticity. CDT calculations on the high-spin states of the neutral metallocenothiaporphyrins with V, Mn, Co, Tc or Rh also shows that these molecules follow a unified extended Hückel and Baird aromaticity orbital-count rule stating that molecules with an odd number of occupied conjugated valence orbitals are aromatic, whereas molecules with an even number of occupied conjugated orbitals are antiaromatic.

Gas-Phase Peroxyl Radical Recombination Reactions: A Computational Study of Formation and Decomposition of Tetroxides.

The recombination ("dimerization") of peroxyl radicals (RO2•) is one of the pathways suggested in the literature for the formation of peroxides (ROOR', often referred to as dimers or accretion products in the literature) in the atmosphere. It is generally accepted that these dimers play a major role in the first steps of the formation of submicron aerosol particles. However, the precise reaction pathways and energetics of RO2• + R'O2• reactions are still unknown. In this work, we have studied the formation of tetroxide intermediates (RO4R'): their formation from two peroxyl radicals and their decomposition to triplet molecular oxygen (3O2) and a triplet pair of alkoxyl radicals (RO•). We demonstrate this mechanism for several atmospherically relevant primary and secondary peroxyl radicals. The potential energy surface corresponds to an overall singlet state. The subsequent reaction channels of the alkoxyl radicals include, but are not limited to, their dimerization into ROOR'. Our work considers the multiconfigurational character of the tetroxides and the intermediate phases of the reaction, leading to reliable mechanistic insights for the formation and decomposition of the tetroxides. Despite substantial uncertainties in the computed energetics, our results demonstrate that the barrier heights along the reaction path are invariably small for these systems. This suggests that the reaction mechanism, previously validated at a multireference level only for methyl peroxyl radicals, is a plausible pathway for the formation of aerosol-relevant larger peroxides in the atmosphere.

Isomer-Resolved Mobility-Mass Analysis of α-Pinene Ozonolysis Products.

Highly oxygenated organic molecules (HOMs) are important sources of atmospheric aerosols. Resolving the molecular-level formation mechanisms of these HOMs from freshly emitted hydrocarbons improves the understanding of aerosol properties and their influence on the climate. In this study, we measure the electrical mobility and mass-to-charge ratio of α-pinene oxidation products using a secondary electrospray-differential mobility analyzer-mass spectrometer (SESI-DMA-MS). The mass-mobility spectrum of the oxidation products is measured with seven different reagent ions generated by the electrospray. We analyzed the mobility-mass spectra of the oxidation products C9-10H14-18O2-6. Our results show that acetate and chloride yield the highest charging efficiencies. Analysis of the mobility spectra suggests that the clusters have 1-5 isomeric structures (i.e., ion-molecule cluster structures with distinct mobilities), and the number is affected by the reagent ion. Most of the isomers are likely cluster isomers originating from binding of the reagent ion to different sites of the molecule. By comparing the number of observed isomers and measured mobilities and collision cross sections between standard pinanediol and pinonic acid to the values observed for C10H18O2 and C10H16O3 produced from oxidation of α-pinene, we confirm that pinanediol and pinonic acid are the only isomers for these elemental compositions in our experimental conditions. Our study shows that the SESI-DMA-MS produces new information from the first steps of oxidation of α-pinene.

Odd-Number Cyclo[n]Carbons Sustaining Alternating Aromaticity.

Cyclo[n]carbons (n = 5, 7, 9, ..., 29) composed from an odd number of carbon atoms are studied computationally at density functional theory (DFT) and ab initio complete active space self-consistent field (CASSCF) levels of theory to get insight into their electronic structure and aromaticity. DFT calculations predict a strongly delocalized carbene structure of the cyclo[n]carbons and an aromatic character for all of them. In contrast, calculations at the CASSCF level yield geometrically bent and electronically localized carbene structures leading to an alternating double aromaticity of the odd-number cyclo[n]carbons. CASSCF calculations yield a singlet electronic ground state for the studied cyclo[n]carbons except for C25, whereas at the DFT level the energy difference between the lowest singlet and triplet states depends on the employed functional. The BHandHLYP functional predicts a triplet ground state of the larger odd-number cyclo[n]carbons starting from n = 13. Current-density calculations at the BHandHLYP level using the CASSCF-optimized molecular structures show that there is a through-space delocalization in the cyclo[n]carbons. The current density avoids the carbene carbon atom, leading to an alternating double aromaticity of the odd-number cyclo[n]carbons satisfying the antiaromatic [4k+1] and aromatic [4k+3] rules. C11, C15, and C19 are aromatic and can be prioritized in future synthesis. We predict a bond-shift phenomenon for the triplet state of the cyclo[n]carbons leading to resonance structures that have different reactivity toward dimerization.

Molecular-scale evidence of aerosol particle formation via sequential addition of HIO3.

Homogeneous nucleation and subsequent cluster growth leads to the formation of new aerosol particles in the atmosphere. The nucleation of sulfuric acid and organic vapours is thought to be responsible for the formation of new particles over continents, whereas iodine oxide vapours have been implicated in particle formation over coastal regions. The molecular clustering pathways that are involved in atmospheric particle formation have been elucidated in controlled laboratory studies of chemically simple systems, but direct molecular-level observations of nucleation in atmospheric field conditions that involve sulfuric acid, organic or iodine oxide vapours have yet to be reported. Here we present field data from Mace Head, Ireland, and supporting data from northern Greenland and Queen Maud Land, Antarctica, that enable us to identify the molecular steps involved in new particle formation in an iodine-rich, coastal atmospheric environment. We find that the formation and initial growth process is almost exclusively driven by iodine oxoacids and iodine oxide vapours, with average oxygen-to-iodine ratios of 2.4 found in the clusters. On the basis of this high ratio, together with the high concentrations of iodic acid (HIO3) observed, we suggest that cluster formation primarily proceeds by sequential addition of HIO3, followed by intracluster restructuring to I2O5 and recycling of water either in the atmosphere or on dehydration. Our study provides ambient atmospheric molecular-level observations of nucleation, supporting the previously suggested role of iodine-containing species in the formation of new aerosol particles, and identifies the key nucleating compound.

Unexpected quenching effect on new particle formation from the atmospheric reaction of methanol with SO3.

Despite the high abundance in the atmosphere, alcohols in general and methanol in particular are believed to play a small role in atmospheric new particle formation (NPF) largely due to the weak binding abilities of alcohols with the major nucleation precursors, e.g., sulfuric acid (SA) and dimethylamine (DMA). Herein, we identify a catalytic reaction that was previously overlooked, namely, the reaction between methanol and SO3, catalyzed by SA, DMA, or water. We found that alcohols can have unexpected quenching effects on the NPF process, particularly in dry and highly polluted regions with high concentrations of alcohols. Specifically, the catalytic reaction between methanol and SO3 can convert methanol into a less-volatile species--methyl hydrogen sulfate (MHS). The latter was initially thought to be a good nucleation agent for NPF. However, our simulation results suggest that the formation of MHS consumes an appreciable amount of atmospheric SO3, disfavoring further reactions of SO3 with H2O. Indeed, we found that MHS formation can cause a reduction of SA concentration up to 87%, whereas the nucleation ability of MHS toward new particles is not as good as that of SA. Hence, a high abundance of methanol in the atmosphere can lower the particle nucleation rate by as much as two orders of magnitude. Such a quenching effect suggests that the recently identified catalytic reactions between alcohols and SO3 need to be considered in atmospheric modeling in order to predict SA concentration from SO2, while also account for their potentially negative effect on NPF.

Highly Oxygenated Organic Molecules (HOM) from Gas-Phase Autoxidation Involving Peroxy Radicals: A Key Contributor to Atmospheric Aerosol.

Highly oxygenated organic molecules (HOM) are formed in the atmosphere via autoxidation involving peroxy radicals arising from volatile organic compounds (VOC). HOM condense on pre-existing particles and can be involved in new particle formation. HOM thus contribute to the formation of secondary organic aerosol (SOA), a significant and ubiquitous component of atmospheric aerosol known to affect the Earth's radiation balance. HOM were discovered only very recently, but the interest in these compounds has grown rapidly. In this Review, we define HOM and describe the currently available techniques for their identification/quantification, followed by a summary of the current knowledge on their formation mechanisms and physicochemical properties. A main aim is to provide a common frame for the currently quite fragmented literature on HOM studies. Finally, we highlight the existing gaps in our understanding and suggest directions for future HOM research.

Heterogeneous Nucleation of Butanol on NaCl: A Computational Study of Temperature, Humidity, Seed Charge, and Seed Size Effects.

Using a combination of quantum chemistry and cluster size distribution dynamics, we study the heterogeneous nucleation of n-butanol and water onto sodium chloride (NaCl)10 seeds at different butanol saturation ratios and relative humidities. We also investigate how the heterogeneous nucleation of butanol is affected by the seed size through comparing (NaCl)5, (NaCl)10, and (NaCl)25 seeds and by seed electrical charge through comparing (Na10Cl9)+, (NaCl)10, and (Na9Cl10)- seeds. Butanol is a common working fluid for condensation particle counters used in atmospheric aerosol studies, and NaCl seeds are frequently used for calibration purposes and as model systems, for example, sea spray aerosol. In general, our simulations reproduce the experimentally observed trends for the NaCl-BuOH-H2O system, such as the increase of nucleation rate with relative humidity and with temperature (at constant supersaturation of butanol). Our results also provide molecular-level insights into the vapor-seed interactions driving the first steps of the heterogeneous nucleation process. The main purpose of this work is to show that theoretical studies can provide molecular understanding of initial steps of heterogeneous nucleation and that it is possible to find cost-effective yet accurate-enough combinations of methods for configurational sampling and energy evaluation to successfully model heterogeneous nucleation of multicomponent systems. In the future, we anticipate that such simulations can also be extended to chemically more complex seeds.

Computational Investigation of the Formation of Peroxide (ROOR) Accretion Products in the OH- and NO3-Initiated Oxidation of α-Pinene.

The formation of accretion products ("dimers") from recombination reactions of peroxyl radicals (RO2) is a key step in the gas-phase generation of low-volatility vapors, leading to atmospheric aerosol particles. We have recently demonstrated that this recombination reaction very likely proceeds via an intermediate complex of two alkoxy radicals (RO···OR') and that the accretion product pathway involves an intersystem crossing (ISC) of this complex from the triplet to the singlet surface. However, ISC rates have hitherto not been computed for large and chemically complex RO···OR' systems actually relevant to atmospheric aerosol formation. Here, we carry out systematic conformational sampling and ISC rate calculations on 3(RO···OR') clusters formed in the recombination reactions of different diastereomers of the first-generation peroxyl radicals originating in both OH- and NO3-initiated reactions of α-pinene, a key biogenic hydrocarbon for atmospheric aerosol formation. While we find large differences between the ISC rates of different diastereomer pairs, all systems have ISC rates of at least 106 s-1, and many have rates exceeding 1010 s-1. Especially the latter value demonstrates that accretion product formation via the suggested pathway is a competitive process also for α-pinene-derived RO2 and likely explains the experimentally observed gas-phase formation of C20 compounds in α-pinene oxidation.