Predicting liquid-liquid phase separation in ternary organic-organic-water mixtures.

Liquid-liquid phase separation (LLPS) affects the water uptake of aerosol particles in the atmosphere through Kelvin and Raoult effects. This study investigates LLPS in ternary mixtures containing water and two organic compounds, using a conductor-like screening model for real solvents (COSMO-RS). COSMO-RS found LLPS in all of the studied mixtures containing water and proxies for primary and secondary organic aerosol (POA and SOA, respectively), due to the limited solubility of the hydrophobic POA proxies in water. The computations predict additional three-phase states in some of the SOA-POA-water mixtures at relative humidity (RH) close to 100%, which was not observed in experiments, likely due to the relatively low RH (90%) used in the experiments. A computational method, such as COSMO-RS, allows for the estimation of new information on mixing states and mixtures that cannot be accessed experimentally. Comparison with experiments can also provide insight into which types of compounds may be present in SOA. Additionally, the possibility of LLPS can be assessed faster with rough estimates rather than by computing the whole phase diagram.

Machine Learning for Predicting Chemical Potentials of Multifunctional Organic Compounds in Atmospherically Relevant Solutions.

We have trained the Extreme Minimum Learning Machine (EMLM) machine learning model to predict chemical potentials of individual conformers of multifunctional organic compounds containing carbon, hydrogen, and oxygen. The model is able to predict chemical potentials of molecules that are in the size range of the training data with a root-mean-square error (RMSE) of 0.5 kcal/mol. There is also a linear correlation between calculated and predicted chemical potentials of molecules that are larger than those included in the training set. Finding the lowest chemical potential conformers is useful in condensed phase thermodynamic property calculations, in order to reduce the number of computationally demanding density functional theory calculations.

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.

New Particle Formation and Growth from Dimethyl Sulfide Oxidation by Hydroxyl Radicals.

Dimethyl sulfide (DMS) is produced by plankton in oceans and constitutes the largest natural emission of sulfur to the atmosphere. In this work, we examine new particle formation from the primary pathway of oxidation of gas-phase DMS by OH radicals. We particularly focus on particle growth and mass yield as studied experimentally under dry conditions using the atmospheric simulation chamber AURA. Experimentally, we show that aerosol mass yields from oxidation of 50-200 ppb of DMS are low (2-7%) and that particle growth rates (8.2-24.4 nm/h) are comparable with ambient observations. An HR-ToF-AMS was calibrated using methanesulfonic acid (MSA) to account for fragments distributed across both the organic and sulfate fragmentation table. AMS-derived chemical compositions revealed that MSA was always more dominant than sulfate in the secondary aerosols formed. Modeling using the Aerosol Dynamics, gas- and particle-phase chemistry kinetic multilayer model for laboratory CHAMber studies (ADCHAM) indicates that the Master Chemical Mechanism gas-phase chemistry alone underestimates experimentally observed particle formation and that DMS multiphase and autoxidation chemistry is needed to explain observations. Based on quantum chemical calculations, we conclude that particle formation from DMS oxidation in the ambient atmosphere will most likely be driven by mixed sulfuric acid/MSA clusters clustering with both amines and ammonia.

Gas-to-Particle Partitioning of Cyclohexene- and α-Pinene-Derived Highly Oxygenated Dimers Evaluated Using COSMOtherm.

Oxidized organic compounds are expected to contribute to secondary organic aerosol (SOA) if they have sufficiently low volatilities. We estimated saturation vapor pressures and activity coefficients (at infinite dilution in water and a model water-insoluble organic phase) of cyclohexene- and α-pinene-derived accretion products, "dimers", using the COSMOtherm19 program. We found that these two property estimates correlate with the number of hydrogen bond-donating functional groups and oxygen atoms in the compound. In contrast, when the number of H-bond donors is fixed, no clear differences are seen either between functional group types (e.g., OH or OOH as H-bond donors) or the formation mechanisms (e.g., gas-phase radical recombination vs liquid-phase closed-shell esterification). For the cyclohexene-derived dimers studied here, COSMOtherm19 predicts lower vapor pressures than the SIMPOL.1 group-contribution method in contrast to previous COSMOtherm estimates using older parameterizations and nonsystematic conformer sampling. The studied dimers can be classified as low, extremely low, or ultra-low-volatility organic compounds based on their estimated saturation mass concentrations. In the presence of aqueous and organic aerosol particles, all of the studied dimers are likely to partition into the particle phase and thereby contribute to SOA formation.

SO2 formation and peroxy radical isomerization in the atmospheric reaction of OH radicals with dimethyl disulfide.

The atmospheric reaction of OH radicals with dimethyl disulfide, CH3SSCH3, proceeds primarily via OH addition forming CH3S and CH3SOH as reactive intermediates, and to a lesser extent via H-abstraction resulting in the peroxy radical CH3SSCH2OO in the presence of O2. The latter undergoes a fast two-step isomerization process leading to HOOCH2SSCHO. CH3S and CH3SOH are both converted to SO2 and CH3O2 with near unity yields under atmospheric conditions.

Improving Solubility and Activity Estimates of Multifunctional Atmospheric Organics by Selecting Conformers in COSMOtherm.

We estimated aqueous solubilities and activity coefficients of atmospherically relevant highly oxidized multifunctional organic compounds in binary mixtures with water at temperatures between 278.15 and 338.15 K, using the COSMOtherm program. Physicochemical properties of organic aerosol constituents are needed in the modeling of atmospheric aerosol processes. As experimental data are often impossible to obtain, reliable estimates from theoretical approaches are a promising path to fill this gap. We investigated the effect of intramolecular hydrogen bonds on the estimation of these condensed-phase properties, attempting to improve the agreement between experimental and estimated values. Citric, tartaric, malic, and maleic acids, which are often used in atmospheric models as representatives of oxidized compounds, were selected to benchmark our calculations. In addition, we estimated aqueous solubilities and activity coefficients of α-pinene-derived organosulfates and highly oxidized isoprene-derived organic compounds, for which no experimental data are available. Our results indicate that the absolute aqueous solubility and activity coefficient estimates of citric, tartaric, malic, and maleic acids, and likely other multifunctional organics, can be improved significantly by selecting conformers on the basis of their intramolecular hydrogen bonding in COSMOtherm calculations.

Solubility and Activity Coefficients of Atmospheric Surfactants in Aqueous Solution Evaluated Using COSMOtherm.

Fatty acids (CH3(CH2)n-2COOH) and their salts are an important class of atmospheric surfactants. Here, we use COSMOtherm to predict solubility and activity coefficients for C2-C12 fatty acids with even number of carbon atoms and their sodium salts in binary water solutions and also in ternary water-inorganic salt solutions. COSMOtherm is a continuum solvent model implementation which can calculate properties of complex systems using quantum chemistry and thermodynamics. Calculated solubility values of the organic acids in pure water are in good agreement with reported experimental values. The comparison of the COSMOtherm-derived Setschenow constants for ternary solutions comprising NaCl with the corresponding experimental values from the literature shows that COSMOtherm overpredicts the salting out effect in all cases except for the solutions of acetic acid. The calculated activity and mean activity coefficients of fatty acids and fatty acid sodium salts, respectively, show deviation of the systems from ideal solution. The computed mean activity coefficients of the fatty acid salts in binary systems are in better agreement with experimentally derived values for the organic salts with longer aliphatic chain (C8-C10). The deviation of the solutions from ideality could lead to biased estimations of cloud condensation nuclei number concentrations if not considered in Köhler calculations and cloud microphysics.

Strong Even/Odd Pattern in the Computed Gas-Phase Stability of Dicarboxylic Acid Dimers: Implications for Condensation Thermodynamics.

The physical properties of small straight-chain dicarboxylic acids are well known to exhibit even/odd alternations with respect to the carbon chain length. For example, odd numbered diacids have lower melting points and higher saturation vapor pressures than adjacent even numbered diacids. This alternation has previously been explained in terms of solid-state properties, such as higher torsional strain of odd number diacids. Using quantum chemical methods, we demonstrate an additional contribution to this alternation in properties resulting from gas-phase dimer formation. Due to a combination of hydrogen bond strength and torsional strain, dimer formation in the gas phase occurs efficiently for glutaric acid (C5) and pimelic acid (C7) but is unfavorable for succinic acid (C4) and adipic acid (C6). Our results indicate that a significant fraction of the total atmospheric gas-phase concentration of glutaric and pimelic acid may consist of dimers.

The role of highly oxygenated organic molecules in the Boreal aerosol-cloud-climate system.

Over Boreal regions, monoterpenes emitted from the forest are the main precursors for secondary organic aerosol (SOA) formation and the primary driver of the growth of new aerosol particles to climatically important cloud condensation nuclei (CCN). Autoxidation of monoterpenes leads to rapid formation of Highly Oxygenated organic Molecules (HOM). We have developed the first model with near-explicit representation of atmospheric new particle formation (NPF) and HOM formation. The model can reproduce the observed NPF, HOM gas-phase composition and SOA formation over the Boreal forest. During the spring, HOM SOA formation increases the CCN concentration by ~10 % and causes a direct aerosol radiative forcing of -0.10 W/m2. In contrast, NPF reduces the number of CCN at updraft velocities < 0.2 m/s, and causes a direct aerosol radiative forcing of +0.15 W/m2. Hence, while HOM SOA contributes to climate cooling, NPF can result in climate warming over the Boreal forest.

NH 4 + Association and Proton Transfer Reactions With a Series of Organic Molecules.

In this study, we present reactions of NH 4 + with a series of analytes (A): acetone (C3H6O), methyl vinyl ketone (C4H6O), methyl ethyl ketone (C4H8O), and eight monoterpene isomers (C10H16) using a Selective Reagent Ionization Time-of-Flight Mass Spectrometer (SRI-ToF-MS). We studied the ion-molecule reactions at collision energies of 55 and 80 meV. The ketones, having a substantially lower proton affinity than NH3, produce only cluster ions NH 4 + (A) in detectable amounts at 55 meV. At 80 meV, no cluster ions were detected meaning that these adduct ions are formed by strongly temperature dependent association reactions. Bond energies of cluster ions and proton affinities for most monoterpenes are not known and were estimated by high level quantum chemical calculations. The calculations reveal monoterpene proton affinities, which range from slightly smaller to substantially higher than the proton affinity of NH3. Proton affinities and cluster bond energies allow to group the monoterpenes as a function of the enthalpy for the dissociation reaction N H 4 + A → A H + + N H 3 . We find that this enthalpy can be used to predict the NH 4 + (A) cluster ion yield. The present study explains product ion formation involving NH 4 + ion chemistry. This is of importance for chemical ionization mass spectrometry (CIMS) utilizing NH 4 + as well as NH 4 + (H2O) as reagent ions to quantitatively detect atmospherically important organic compounds in real-time.

Computational Comparison of Different Reagent Ions in the Chemical Ionization of Oxidized Multifunctional Compounds.

High pressure anion chemical ionization is commonly used for the detection of neutral molecules in the gas phase. The detection efficiency in these measurements depends on how strongly the reagent ion binds to the neutral target molecule. We have calculated the binding strength of nitrate (NO3-), acetate (CH3C(O)O-), lactate (CH3CH(OH)C(O)O-), trifluoroacetate (CF3C(O)O-), trifluoromethanolate (CF3O-), bromide (Br-), and iodide (I-) reagent ions to ten different products derived from the OH radical-initiated oxidation of butadiene. We found that the binding of these oxidation products to the reagent ions depends almost linearly on the number of oxygen atoms in the target molecule, with the precise chemical identity of the compound (e.g., the number and relative position of hydroxyl or hydroperoxy groups) playing a more minor role. For acetate, the formation free energy decreases on average by around 4 kcal/mol when the number of oxygen atoms in the sample molecule increases by one. For the other reagent ions the corresponding decrease is around 3 kcal/mol. For all of the molecules studied, acetate forms the most stable clusters and I- the least stable. We also investigated the effect of humidity on the chemical ionization by calculating how strongly water molecules bind to both the reagent ions and the ion-molecule clusters. Water binds much more strongly to the reagent ion monomers compared to the reagent ion "dimers" (defined here as a cluster of the reagent anion with the corresponding neutral conjugate acid, e.g., HNO3(NO3-)) or the ion-molecule clusters. This likely leads to a stronger humidity dependence when using reagent ions that are not able to form reagent ion dimers (such as CF3C(O)O-, CF3O-, Br-, and I-).

Computational and Experimental Investigation of the Detection of HO2 Radical and the Products of Its Reaction with Cyclohexene Ozonolysis Derived RO2 Radicals by an Iodide-Based Chemical Ionization Mass Spectrometer.

The HO2 radical is an important atmospheric molecule that can potentially influence the termination of autoxidation processes of volatile organic compounds (VOCs) that lead to the formation of highly oxygenated multifunctional compounds (HOMs). In this work, we demonstrate the direct detection of the HO2 radical using an iodide-based chemical ionization mass spectrometer (iodide-CIMS). Expanding on the previously established correlation between molecule-iodide binding enthalpy and iodide-CIMS instrument sensitivity, the experimental detection of the HO2 radical was preceded by the quantum chemical calculation of the HO2*I- cluster (PBE/aug-cc-pVTZ-PP level), which showed a reasonably strong binding enthalpy of 21.60 kcal/mol. Cyclohexene ozonolysis intermediates and closed-shell products were next detected by the iodide-CIMS. The ozone-initiated cyclohexene oxidation mechanism was perturbed by the introduction of the HO2 radical, leading to the formation of closed-shell hydroperoxides. The experimental investigation once again followed the initial computational molecule-iodide binding enthalpy calculations. The quantum chemical calculations were performed at the PBE/aug-cc-pVTZ-PP level for radicals and DLPNO-CCSD(T)/def2-QZVPP//PBE/aug-cc-pVTZ-PP level for the closed-shell products. A comparison between the iodide-CIMS and nitrate-CIMS spectra with identical measurement steps revealed that the iodide-CIMS was able to detect the low-oxidized (O/C ratio 0.5 and 0.66) cyclohexene ozonolysis monomer products more efficiently than nitrate-CIMS. Higher-oxidized monomers (O/C ratio 1 to 1.5) were detected equally well by both methods. An investigation of dimers showed that both iodide- and nitrate-CIMS were able to detect the dimer compositions possibly formed from reactions between the peroxy radical monomers considered in this study. Additionally, iodide-CIMS detected organic ions that were formed by a previously suggested mechanism of dehydroxylation of peroxy acids (and deoxygenation of acyl peroxy radicals) by H2O*I- clusters. These mechanisms were computationally verified.

Computational Comparison of Acetate and Nitrate Chemical Ionization of Highly Oxidized Cyclohexene Ozonolysis Intermediates and Products.

During the past few years nitrate chemical ionization has been used to detect highly oxidized products from OH- and O3-initiated alkene autoxidation. These have been speculated to play a significant role in atmospheric aerosol formation. As less oxidized autoxidation products have not been detected using nitrate chemical ionization, and the absolute concentrations of the highly oxidized species are as yet unknown, other reagent ions, such as acetate, are needed both to verify the detection efficiency of nitrate chemical ionization and to measure the less oxidized compounds. Here we compare the formation free energies of the acetate and nitrate clusters of several atmospherically relevant RO2 intermediates and products derived from cyclohexene ozonolysis, calculated at the ωB97xD/aug-cc-pVTZ level of theory. We found that, for the molecules with one hydrogen bonding peroxy acid group, the binding with nitrate is on average 7.5 kcal/mol weaker than with acetate and the binding is on average 10.5 kcal/mol weaker for molecules with two hydrogen bonding peroxy acid groups. We also calculated the deprotonation energies of the RO2 intermediates and the closed-shell products and found that acetate is able to deprotonate almost all of these molecules, while deprotonation with nitrate is (as expected for the conjugate base of a strong acid) not favorable.

Cost-Effective Implementation of Multiconformer Transition State Theory for Peroxy Radical Hydrogen Shift Reactions.

Based on a small test system, (R)-CH(OH)(OO·)CH2CHO, we have developed a cost-effective approach to the practical implementation of multiconformer transition state theory for peroxy radical hydrogen shift reactions at atmospherically relevant temperatures. While conformer searching is crucial for accurate reaction rates, an energy cutoff can be used to significantly reduce the computational cost with little loss of accuracy. For the reaction barrier, high-level calculations are needed, but the highest level of electronic structure theory is not necessary for the relative energy between conformers. Improving the approach to both transition state theory and electronic structure theory decreases the calculated reaction rate significantly, so low-level calculations can be used to rule out slow reactions. Further computational time can be saved by approximating the tunneling coefficients for each transition state by only that of the lowest-energy transition state. Finally, we test and validate our approach using higher-level theoretical values for our test system and existing experimental results for additional peroxy radical hydrogen shift reactions in three slightly larger systems.

Hydroxyl radical-induced formation of highly oxidized organic compounds.

Explaining the formation of secondary organic aerosol is an intriguing question in atmospheric sciences because of its importance for Earth's radiation budget and the associated effects on health and ecosystems. A breakthrough was recently achieved in the understanding of secondary organic aerosol formation from ozone reactions of biogenic emissions by the rapid formation of highly oxidized multifunctional organic compounds via autoxidation. However, the important daytime hydroxyl radical reactions have been considered to be less important in this process. Here we report measurements on the reaction of hydroxyl radicals with α- and ß-pinene applying improved mass spectrometric methods. Our laboratory results prove that the formation of highly oxidized products from hydroxyl radical reactions proceeds with considerably higher yields than previously reported. Field measurements support these findings. Our results allow for a better description of the diurnal behaviour of the highly oxidized product formation and subsequent secondary organic aerosol formation in the atmosphere.

Unimolecular HO2 Loss from Peroxy Radicals Formed in Autoxidation Is Unlikely under Atmospheric Conditions.

A concerted HO2 loss reaction from a peroxy radical (RO2), formed from the addition of O2 to an alkyl radical, has been proposed as a mechanism to form closed-shell products in the atmospheric oxidation of organic molecules. We investigate this reaction computationally with four progressively oxidized radicals. Potential energy surfaces of the O2 addition and HO2 loss reactions were calculated at ROHF-RCCSD(T)-F12a/VDZ-F12//ωB97xD/aug-cc-pVTZ level of theory and the master equation solver for multienergy well reactions (MESMER) was used to calculate Bartis-Widom phenomenological rate coefficients. The rate coefficients were also compared with the unimolecular rate coefficients of the HO2 loss reaction calculated with transition state theory at atmospheric temperature and pressure. On the basis of our calculations, the unimolecular concerted HO2 loss is unlikely to be a major pathway in the formation of highly oxidized closed-shell molecules in the atmosphere.

Computational Study of Hydrogen Shifts and Ring-Opening Mechanisms in α-Pinene Ozonolysis Products.

Autoxidation by sequential peroxy radical hydrogen shifts (H-shifts) and O2 additions has recently emerged as a promising mechanism for the rapid formation of highly oxidized, low-volatility organic compounds in the atmosphere. A key prerequisite for autoxidation is that the H-shifts of the initial peroxy radicals formed by, e.g., OH or O3 oxidation are fast enough to compete with bimolecular sink reactions. In most atmospheric conditions, these restrict the lifetime of peroxy radicals to be on the order of seconds. We have systematically investigated all potentially important (nonmethyl, sterically unhindered) H-shifts of all four peroxy radicals formed in the ozonolysis of α-pinene using density functional (ωB97XD) and coupled cluster [CCSD(T)-F12] theory. In contrast to the related but chemically simpler cyclohexene ozonolysis system, none of the calculated H-shifts have rate constants above 1 s(-1) at 298 K, and most are below 0.01 s(-1). The low rate constants are connected to the presence of the strained cyclobutyl ring in the α-pinene-derived peroxy radicals, which hinders H-shifts both from and across the ring. For autoxidation to yield the experimentally observed highly oxidized products in the α-pinene ozonolysis system, additional ring-opening reaction mechanisms breaking the cyclobutyl ring are therefore needed. We further investigate possible uni- and bimolecular pathways for opening the cyclobutyl ring in the α-pinene ozonolysis system.

Computational Study of the Clustering of a Cyclohexene Autoxidation Product C6H8O7 with Itself and Sulfuric Acid.

We investigate the molecular interactions between sulfuric acid and a recently reported C6H8O7 ketodiperoxy acid formed through autoxidation from cyclohexene and ozone. Structurally similar but larger ELVOC (extremely low volatility organic compound) products formed from autoxidation of monoterpenes are believed to play a major role in the formation and early growth of atmospheric aerosol particles. Utilizing density functional theory geometries, with a DLPNO-CCSD(T)/def2-QZVPP single point energy correction, the stepwise Gibbs free energies of formation have been calculated, and the stabilities of the molecular clusters have been evaluated. C6H8O7 interacts weakly with both itself and sulfuric acid, with standard free energies of formation (ΔG at 298 K and 1 atm) around or above 0 kcal/mol. This is due to the presence of strong intramolecular hydrogen bonds in the peroxyacid groups of C6H8O7. These stabilize the isolated molecule with respect to its clusters, and lead to unfavorable interaction energies. The addition of sulfuric acid to clusters containing C6H8O7 is somewhat more favorable, but the formed clusters are still far more likely to evaporate than to grow further in atmospheric conditions. These findings indicate that the O/C ratio cannot exclusively be used as a proxy for volatility in atmospheric new particle formation involving organic compounds. The specific molecular structure, and especially the number of strong hydrogen binding moieties, are equally important. The interaction between the C6H8O7 compound and aqueous phase sulfate ions indicates that ELVOC-type compounds can contribute to aerosol mass by effectively partitioning into the condensed phase.

Modeling the Charging of Highly Oxidized Cyclohexene Ozonolysis Products Using Nitrate-Based Chemical Ionization.

Several extremely low volatility organic compounds (ELVOCs) formed in the ozonolysis of endocyclic alkenes have recently been detected in laboratory and field studies. These experiments have been carried out with chemical ionization atmospheric pressure interface time-of-flight mass spectrometers (CI-APi-TOF) with nitrate ions as reagent ions. The nitrate ion binds to the detected species through hydrogen bonds, but it also binds very strongly to one or two neutral nitric acid molecules. This makes the measurement highly selective when there is an excess amount of neutral nitric acid in the instrument. In this work, we used quantum-chemical methods to calculate the binding energies between a nitrate ion and several highly oxidized ozonolysis products of cyclohexene. These were then compared with the binding energies of nitrate ion-nitric acid clusters. Systematic configurational sampling of the molecules and clusters was carried out at the B3LYP/6-31+G* and ωB97xD/aug-cc-pVTZ levels, and the final single-point energies were calculated with DLPNO-CCSD(T)/def2-QZVPP. The binding energies were used in a kinetic simulation of the measurement system to determine the relative ratios of the detected signals. Our results indicate that at least two hydrogen bond donor functional groups (in this case, hydroperoxide, OOH) are needed for an ELVOC molecule to be detected in a nitrate ion CI-APi-TOF. Also, a double bond in the carbon backbone makes the nitrate cluster formation less favorable.