Glass-Forming Properties of 3-Methylbutane-1,2,3-tricarboxylic Acid and Its Mixtures with Water and Pinonic Acid

Volatility Basis Set Distributions and Viscosity of Organic Aerosol Mixtures: Insights from Chemical Characterization Using Temperature-Programmed Desorption-Direct Analysis in Real-Time High-Resolution Mass Spectrometry

Quantitative assessment of gas-particle partitioning of individual components within complex atmospheric organic aerosol (OA) mixtures is critical for predicting and comprehending the formation and evolution of OA particles in the atmosphere. This investigation leverages previously documented data obtained through a temperature-programmed desorption-direct analysis in real-time, high-resolution mass spectrometry (TPD-DART-HRMS) platform. This methodology facilitates the bottom-up construction of volatility basis set (VBS) distributions for constituents found in three biogenic secondary organic aerosol (SOA) mixtures produced through the ozonolysis of α-pinene, limonene, and ocimene. The apparent enthalpies (ΔH*, kJ mol-1) and saturation mass concentrations (CT*, μg·m-3) of individual SOA components, determined as a function of temperature (T, K), facilitated an assessment of changes in VBS distributions and gas-particle partitioning with respect to T and atmospheric total organic mass loadings (tOM, μg·m-3). The VBS distributions reveal distinct differences in volatilities among monomers, dimers, and trimers, categorized into separate volatility bins. At the ambient temperature of T = 298 K, only monomers efficiently partition between gas and particle phases across a broad range of atmospherically relevant tOM values of 1-100 μg·m-3. Partitioning of dimers and trimers becomes notable only at T > 360 K and T > 420 K, respectively. The viscosity of SOA mixtures is assessed using a bottom-up calculation approach, incorporating the input of elemental formulas, ΔH*, CT*, and particle-phase mass fractions of the SOA components. Through this approach, we are able to accurately estimate the variations in SOA viscosity that result from the evaporation of its components. These variations are, in turn, influenced by atmospherically relevant changes in tOM and T. Comparison of the calculated SOA viscosity and diffusivity values with literature reported experimental results shows close agreement, thereby validating the employed calculation approach. These findings underscore the significant potential for TPD-DART-HRMS measurements in enabling the untargeted analysis of organic molecules within OA mixtures. This approach facilitates quantitative assessment of their gas-particle partitioning and allows for the estimation of their viscosity and condensed-phase diffusion, thereby contributing valuable insights to atmospheric models.

Modeling the molecular composition of secondary organic aerosol under highly polluted conditions: A case study in the Yangtze River Delta Region in China

A near-explicit mechanism, the master chemical mechanism (MCMv3.3.1), coupled with the Community Multiscale Air Quality (CMAQ) model (CMAQ-MCM-SOA), was applied to investigate the characteristics of secondary organic aerosol (SOA) during a pollution event in the Yangtze River Delta (YRD) region in summer 2018. Model performances in predicting explicit volatile organic compounds (VOCs), organic aerosol (OA), secondary organic carbon (SOC), and other related pollutants in Taizhou, as well as ozone (O3) and fine particulate matter (PM2.5) in multiple cities in this region, were evaluated against observations and model predictions by the CMAQ model coupled with a lumped photochemical mechanism (SAPRC07tic, S07). MCM and S07 exhibited similar performances in predicting gaseous species, while MCM better captured the observed PM2.5 and inorganic aerosols. Both models underpredicted OA concentrations. When excluding data during biomass burning events, SOC concentrations were underpredicted by the CMAQ-MCM-SOA model (-28.4 %) and overpredicted by the CMAQ-S07 model (134.4 %), with better agreement with observations in the trend captured by the CMAQ-MCM-SOA model. Dicarbonyl SOA accounted for a significant fraction of total SOA in the YRD, while organic nitrates originating from aromatics were the most abundant species contributing to the SOA formation from gas-particle partitioning. The oxygen-to‑carbon ratio (O/C) for SOA and OA were 0.68-0.75 and 0.20-0.65, respectively, indicating a higher oxidation state in the areas influenced by biogenic emissions. Finally, the phase state of SOA was examined by calculating the glass transition temperature (Tg) based on its molecular composition. It was found that semi-solid state characterized SOA in the YRD, which could potentially impact their chemical transformation and lifetimes along with those of their precursors.

Glass Transition Temperatures of Organic Mixtures from Isoprene Epoxydiol-Derived Secondary Organic Aerosol

The phase states and glass transition temperatures (Tg) of secondary organic aerosol (SOA) particles are important to resolve for understanding the formation, growth, and fate of SOA as well as their cloud formation properties. Currently, there is a limited understanding of how Tg changes with the composition of organic and inorganic components of atmospheric aerosol. Using broadband dielectric spectroscopy, we measured the Tg of organic mixtures containing isoprene epoxydiol (IEPOX)-derived SOA components, including 2-methyltetrols (2-MT), 2-methyltetrol-sulfate (2-MTS), and 3-methyltetrol-sulfate (3-MTS). The results demonstrate that the Tg of mixtures depends on their composition. The Kwei equation, a modified Gordon-Taylor equation with an added quadratic term and a fitting parameter representing strong intermolecular interactions, provides a good fit for the Tg-composition relationship of complex mixtures. By combining Raman spectroscopy with geometry optimization simulations obtained using density functional theory, we demonstrate that the non-linear deviation of Tg as a function of composition may be caused by changes in the extent of hydrogen bonding in the mixture.

Machine-Learning-Based Prediction of the Glass Transition Temperature of Organic Compounds Using Experimental Data

Knowledge of the glass transition temperature of molecular compounds that occur in atmospheric aerosol particles is important for estimating their viscosity, as it directly influences the kinetics of chemical reactions and particle phase state. While there is a great diversity of organic compounds present in aerosol particles, for only a minor fraction of them experimental glass transition temperatures are known. Therefore, we have developed a machine learning model designed to predict the glass transition temperature of organic molecular compounds based on molecule-derived input variables. The extremely randomized trees (extra trees) procedure was chosen for this purpose. Two approaches using different sets of input variables were followed. The first one uses the number of selected functional groups present in the compound, while the second one generates descriptors from a SMILES (Simplified Molecular Input Line Entry System) string. Organic compounds containing carbon, hydrogen, oxygen, nitrogen, and halogen atoms are included. For improved results, both approaches can be combined with the melting temperature of the compound as an additional input variable. The results show that the predictions of both approaches show a similar mean absolute error of about 12-13 K, with the SMILES-based predictions performing slightly better. In general, the model shows good predictive power considering the diversity of the experimental input data. Furthermore, we also show that its performance exceeds that of previous parameterizations developed for this purpose and also performs better than existing machine learning models. In order to provide user-friendly versions of the model for applications, we have developed a web site where the model can be run by interested scientists via a web-based interface without prior technical knowledge. We also provide Python code of the model. Additionally, all experimental input data are provided in form of the Bielefeld Molecular Organic Glasses (BIMOG) database. We believe that this model is a powerful tool for many applications in atmospheric aerosol science and material science.

Molecular Understanding of the Enhancement in Organic Aerosol Mass at High Relative Humidity

The mechanistic pathway by which high relative humidity (RH) affects gas-particle partitioning remains poorly understood, although many studies report increased secondary organic aerosol (SOA) yields at high RH. Here, we use real-time, molecular measurements of both the gas and particle phase to provide a mechanistic understanding of the effect of RH on the partitioning of biogenic oxidized organic molecules (from α-pinene and isoprene) at low temperatures (243 and 263 K) at the CLOUD chamber at CERN. We observe increases in SOA mass of 45 and 85% with increasing RH from 10-20 to 60-80% at 243 and 263 K, respectively, and attribute it to the increased partitioning of semi-volatile compounds. At 263 K, we measure an increase of a factor 2-4 in the concentration of C10H16O2-3, while the particle-phase concentrations of low-volatility species, such as C10H16O6-8, remain almost constant. This results in a substantial shift in the chemical composition and volatility distribution toward less oxygenated and more volatile species at higher RH (e.g., at 263 K, O/C ratio = 0.55 and 0.40, at RH = 10 and 80%, respectively). By modeling particle growth using an aerosol growth model, which accounts for kinetic limitations, we can explain the enhancement in the semi-volatile fraction through the complementary effect of decreased compound activity and increased bulk-phase diffusivity. Our results highlight the importance of particle water content as a diluting agent and a plasticizer for organic aerosol growth.

Nanomechanics and Morphology of Simulated Respiratory Particles

The impact of respiratory particle composition on the equilibrium morphology and phase is not well understood. Furthermore, the effects of these different phases and morphologies on the viability of viruses embedded within these particles are equally unknown. Physiologically relevant respiratory fluid analogues were constructed, and their hygroscopic behavior was measured using an ensemble technique. A relationship between hygroscopicity and protein concentration was determined, providing additional validation to the high protein content of respiratory aerosol measured in prior works (>90%). It was found that the salt component of the respiratory particles could crystallize as a single crystal, multiple crystals, or would not crystallize at all. It was found that dried protein particles at indoor-relevant climatic conditions could exist separately in a glassy (∼77% of particles) or viscoelastic state (∼23% of particles). The phase state and morphology of respiratory particles may influence the viability of embedded pathogens. We recommend that pathogen research aiming to mimic the native composition of respiratory fluid should use a protein concentration of at least 90% by solute volume to improve the representativity of the pathogen's microenvironment.

Emerging investigator series: chemical and physical properties of organic mixtures on indoor surfaces during HOMEChem

Organic films on indoor surfaces serve as a medium for reactions and for partitioning of semi-volatile organic compounds and thus play an important role in indoor chemistry. However, the chemical and physical properties of these films are poorly characterized. Here, we investigate the chemical composition of an organic film collected during the HOMEChem campaign, over three cumulative weeks in the kitchen, using both Fourier Transform Ion Cyclotron Resonance Mass Spectrometry (FT-ICR MS) and offline Aerosol Mass Spectrometry (AMS). We also characterize the viscosity of this film using a model based on molecular formulas as well as poke-flow measurements. We find that the film contains organic material similar to cooking organic aerosol (COA) measured during the campaign using on-line AMS. However, the average molecular formula observed using FT-ICR MS is ∼C50H90O11, which is larger and more oxidized than fresh COA. Solvent extracted film material is a low viscous semisolid, with a measured viscosity <104 Pa s. This is much lower than the viscosity model predicts, which is parametrized with atmospherically relevant organic molecules, but sensitivity tests demonstrate that including unsaturation can explain the differences. The presence of unsaturation is supported by reactions of film material with ozone. In contrast to the solvent extract, manually removed material appears to be highly viscous, highlighting the need for continued work understanding both viscosity measurements as well as parameterizations for modeled viscosity of indoor organic films.

Indoor aerosol water content and phase state in U.S. residences: impacts of relative humidity, aerosol mass and composition, and mechanical system operation

Hygroscopic particulate matter (PM) constituents promote uptake of aerosol water (AW), depending on relative humidity (RH), which can constrain qualities such as organic aerosol (OA) phase state and inorganic aerosol (IA) deliquescence and efflorescence. This work provides a first incorporation of AW predictions into residential indoor PM simulations. The indoor model, IMAGES, which simulates factored OA concentrations and thermodynamics using the 2D-volatility basis set, was expanded to predict speciated IA concentrations, AW with κ-Köhler theory of hygroscopic growth, and OA phase state with glass transition temperatures. Since RH is the largest driver of AW and varies with meteorology, simulations were conducted using a database of historical ambient weather and pollution records spanning the sixteen U.S. climate zones, facilitating assessment of seasonal and regional trends. Over this diverse simulation set, the residential indoor AW mass was ∼10 to 100 times smaller than dry PM mass. This relative AW amount indoors was about ∼10 times smaller than outdoors, since indoor-emitted aerosol is likely less hygroscopic. The indoor OA phase state was typically semisolid, suggesting kinetic limitations might inhibit thermodynamic OA partitioning equilibrium from being established indoors. Residences in hot and humid climates during the summertime may have liquid indoor OA, while amorphous solid indoor OA can exist in cold climates. Deliquescence and efflorescence of recirculated IA within HVAC systems during cooling or heating, respectively, was also modeled. Oftentimes, two IA populations with different histories existing as wet or dry aerosol were generated by HVAC operation depending on indoor and outdoor environmental conditions and the HVAC operating mode.

Physicochemical Properties of Terebic Acid, MBTCA, Diaterpenylic Acid Acetate, and Pinanediol as Relevant α-Pinene Oxidation Products

The physicochemical properties and the synthesis of four α-pinene oxidation products, terebic acid, 3-methyl-1,2,3-butanetricarboxylic acid (MBTCA), diaterpenylic acid acetate (DTAA), and pinanediol, are presented in this study. The physicochemical properties encompass thermal properties, solubility in water, and dissociation constant (pK _a) for the investigated compounds. It was found that terebic acid exhibits a relatively high melting temperature of 449.29 K, whereas pinanediol revealed a low melting temperature of 329.26 K. The solubility in water was determined with the dynamic method and the experimental results were correlated using three different mathematical models: Wilson, NRTL, and UNIQUAC equations. The results of the correlation indicate that the Wilson equation appears to work the best for terebic acid and pinanediol. The calculated standard deviation was for 3.79 for terebic acid and 1.25 for pinanediol. In contrast, UNIQUAC was the best mathematical model for DTAA and MBTCA. The calculated standard deviation was 0.57 for DTAA and 2.21 for MBTCA. The measured water solubility increased in the following order: pinanediol > DTAA ≥ MBTCA > terebic acid, which affects their multiphase aging chemistry in the atmosphere. Moreover, acidity constants (pK _a) at 298, 303, and 308 K were determined for DTAA with the Bates-Schwarzenbach spectrophotometric method. The pK _a values obtained at 298, 303, and 308 K were found to be 3.76, 3.85, and 3.88, respectively.

The Cooling Rate- and Volatility-Dependent Glass-Forming Properties of Organic Aerosols Measured by Broadband Dielectric Spectroscopy

Glass transitions of secondary organic aerosols (SOA) from liquid/semisolid to solid phase states have important implications for aerosol reactivity, growth, and cloud formation properties. In the present study, glass transition temperatures (T_g) of isoprene SOA components, including isoprene hydroxy hydroperoxide (ISOPOOH), isoprene-derived epoxydiols (IEPOX), 2-methyltetrols, and 2-methyltetrol sulfates, were measured at atmospherically relevant cooling rates (2-10 K/min) by thin film broadband dielectric spectroscopy. The results indicate that 2-methyltetrol sulfates have the highest glass transition temperature, while ISOPOOH has the lowest glass transition temperature. By varying the cooling rate of the same compound from 2 to 10 K/min, the T_g of these compounds increased by 4-5 K. This temperature difference leads to a height difference of 400-800 m in the atmosphere for the corresponding updraft induced cooling rates, assuming a hygroscopicity value (κ) of 0.1 and relative humidity less than 95%. The T_g of the organic compounds was found to be strongly correlated with volatility, and a semiempirical formula between glass transition temperatures and volatility was derived. The Gordon-Taylor equation was applied to calculate the effect of relative humidity (RH) and water content at five mixing ratios on the T_g of organic aerosols. The model shows that T_g could drop by 15-40 K as the RH changes from <5 to 90%, whereas the mixing ratio of water in the particle increases from 0 to 0.5. These results underscore the importance of chemical composition, updraft rates, and water content (RH) in determining the phase states and hygroscopic properties of organic particles.

Highly Viscous States Affect the Browning of Atmospheric Organic Particulate Matter

Initially transparent organic particulate matter (PM) can become shades of light-absorbing brown via atmospheric particle-phase chemical reactions. The production of nitrogen-containing compounds is one important pathway for browning. Semisolid or solid physical states of organic PM might, however, have sufficiently slow diffusion of reactant molecules to inhibit browning reactions. Herein, organic PM of secondary organic material (SOM) derived from toluene, a common SOM precursor in anthropogenically affected environments, was exposed to ammonia at different values of relative humidity (RH). The production of light-absorbing organonitrogen imines from ammonia exposure, detected by mass spectrometry and ultraviolet-visible spectrophotometry, was kinetically inhibited for RH < 20% for exposure times of 6 min to 24 h. By comparison, from 20% to 60% RH organonitrogen production took place, implying ammonia uptake and reaction. Correspondingly, the absorption index k across 280 to 320 nm increased from 0.012 to 0.02, indicative of PM browning. The k value across 380 to 420 nm increased from 0.001 to 0.004. The observed RH-dependent behavior of ammonia uptake and browning was well captured by a model that considered the diffusivities of both the large organic molecules that made up the PM and the small reactant molecules taken up from the gas phase into the PM. Within the model, large-molecule diffusivity was calculated based on observed SOM viscosity and evaporation. Small-molecule diffusivity was represented by the water diffusivity measured by a quartz-crystal microbalance. The model showed that the browning reaction rates at RH < 60% could be controlled by the low diffusivity of the large organic molecules from the interior region of the particle to the reactive surface region. The results of this study have implications for accurate modeling of atmospheric brown carbon production and associated influences on energy balance.

Global distribution of particle phase state in atmospheric secondary organic aerosols

Secondary organic aerosols (SOA) are a large source of uncertainty in our current understanding of climate change and air pollution. The phase state of SOA is important for quantifying their effects on climate and air quality, but its global distribution is poorly characterized. We developed a method to estimate glass transition temperatures based on the molar mass and molecular O:C ratio of SOA components, and we used the global chemistry climate model EMAC with the organic aerosol module ORACLE to predict the phase state of atmospheric SOA. For the planetary boundary layer, global simulations indicate that SOA are mostly liquid in tropical and polar air with high relative humidity, semi-solid in the mid-latitudes and solid over dry lands. We find that in the middle and upper troposphere SOA should be mostly in a glassy solid phase state. Thus, slow diffusion of water, oxidants and organic molecules could kinetically limit gas–particle interactions of SOA in the free and upper troposphere, promote ice nucleation and facilitate long-range transport of reactive and toxic organic pollutants embedded in SOA.

Secondary organic aerosols (SOA) are important for climate and aerosol quality, but the phase state is unclear. Here, the authors show that SOA is liquid in tropical and polar air, semi-solid in the mid-latitudes, solid over dry lands and in a glassy solid phase state in the middle and upper troposphere.

Influence of Functional Groups on the Viscosity of Organic Aerosol

Organic aerosols can exist in highly viscous or glassy phase states. A viscosity database for organic compounds with atmospherically relevant functional groups is compiled and analyzed to quantify the influence of number and location of functional groups on viscosity. For weakly functionalized compounds the trend in viscosity sensitivity to functional group addition is carboxylic acid (COOH) ≈ hydroxyl (OH) > nitrate (ONO₂) > carbonyl (CO) ≈ ester (COO) > methylene (CH₂). Sensitivities to group addition increase with greater levels of prior functionalization and decreasing temperature. For carboxylic acids a sharp increase in sensitivity is likely present already at the second addition at room temperature. Ring structures increase viscosity relative to linear structures. Sensitivities are correlated with analogously derived sensitivities of vapor pressure reduction. This may be exploited in the future to predict viscosity in numerical models by piggybacking on schemes that track the evolution of organic aerosol volatility with age.