Molecular Orientation at the Squalene/Air Interface from Sum Frequency Generation Spectroscopy and Atomistic Modeling

Molecular Dynamics Simulations of the Interactions of Organic Compounds at Indoor Relevant Surfaces

With markedly different reaction conditions compared to the chemistry of the outside atmosphere, indoor air chemistry poses new challenges to the scientific community that require combined experimental and computational efforts. Here, we review molecular dynamics simulations that have contributed to the mechanistic understanding of the complex dynamics of organic compounds at indoor surfaces and their interplay with experiments and indoor air models. We highlight the rich interactions between volatile organic compounds and silica and titanium dioxide surfaces, serving as proxies for glasses and paints, as well as the dynamics of skin oil lipids and their oxidation products, which sensitively affect the quality of indoor air in crowded environments. As the studies we review here are pioneering in the rapidly emerging field of indoor chemistry, we provide suggestions for increasing the potentially important role that molecular simulations can continue to play.

Indoor surface chemistry variability: microspectroscopic analysis of deposited particles in dwellings across the United States

Dwellings across the United States range dramatically with respect to numerous variables (e.g., size, ventilation, and proximity to outdoor sources), and there are considerable uncertainties regarding the heterogeneity in chemical composition and physical properties of indoor particles and surfaces. Stay-at-home orders early in the COVID-19 pandemic led to significant portions of the population spending high fractions of their time at their primary dwelling. Stay-at-HomeChem leveraged a network of indoor chemistry researchers to study indoor air quality and surface chemistry in their homes (March-April 2020). Within this effort, glass microscope slides were deployed in kitchens and other rooms in dwellings across the country for time periods ranging from as short as three hours up to three weeks. Overall, results from 10 occupied homes (15 distinct rooms) showed that collected material on this time scale was primarily deposited particles, rather than thick films, based on optical microscopy and profilometry. Raman microspectroscopy and optical photothermal infrared (O-PTIR) spectroscopy showed that organic modes were dominant, including ν(C-H), δ(C-H), and ν(CO), with minimal contributions from inorganic ions commonly observed in outdoor particulate matter (sulfate, nitrate, or ammonium). Spectral variability within the C-H stretching and fingerprint regions demonstrate differing compositions of deposited particles, often related to cooking activities (e.g., organic particles from cooking oils). Differences within a single dwelling, highlighted that particles from cooking were key contributors in some other rooms, but not all, reinforcing that sources and ventilation likely led to quite distinct surfaces in different rooms. Overall, these results demonstrate the need for real-world measurements to assess the representativeness of assumptions regarding exposure to organic material indoors.

Molecular Orientation of Carboxylate Anions at the Water-Air Interface Studied with Heterodyne-Detected Vibrational Sum-Frequency Generation

The carboxylate anion group plays an important role in many (bio)chemical systems and polymeric materials. In this work, we study the orientation of carboxylate anions with various aliphatic and aromatic substituents at the water-air interface by probing the carboxylate stretch vibrations with heterodyne-detected vibrational sum-frequency generation spectroscopy in different polarization configurations. We find that carboxylate groups with small aliphatic substituents show a large tilt angle with respect to the surface normal and that this angle decreases with increasing size of the substituent. We further use the information about the orientation of the carboxylate group to determine the hyperpolarizability components of this group.

Unanticipated Hydrophobicity Increases of Squalene and Human Skin Oil Films Upon Ozone Exposure

The C-H and O-H oscillators on the surfaces of thin films of human-derived skin oil and squalene are probed under ambient conditions (300 K, 1 atm total pressure, 40% RH) using second-order vibrational spectroscopy and contact angle goniometry before and after exposure to ppb amounts of ozone. Skin oil and squalene are found to produce different vibrational sum frequency generation spectra in the C-H stretching region, while exposure to ozone results in surface spectra for both materials that is consistent with a loss of C-H oscillators. The measured contact angles show that the hydrophobicity of the films increases following exposure to ozone, consistent with the reduction in C═C···H₂O ("πH") bonding interactions that is expected from C═C double bond loss due to ozonolysis and indicating that the polar functional groups formed point toward the films' interiors. Implications for heterogeneous indoor chemistry are discussed.

Diffusion of Squalene in Nonaqueous Solvents

Capillary flow techniques have been used to determine the translational diffusion constant, D, of squalene in seven alkanes and five cyclohexanes. The alkanes are n-hexane, n-octane, n-decane, n-dodecane, n-tetradecane, 2,2,4,4,6,8,8-heptamethylnonane (isocetane), and 2,6,10,14-tetramethylpentadecane (pristane). The cyclohexanes are cyclohexane, n-butylcyclohexane, n-hexylcyclohexane, n-octylcyclohexane, and n-dodecylcyclohexane. When combined with published data in CD₂Cl₂, ethyl acetate, n-hexadecane, squalane, n-octane-squalane mixtures, and supercritical CO₂, the 35 diffusion constants and viscosities, η, vary by factors of ∼230 and ∼500, respectively. A fit to the modified Stokes-Einstein equation (MSE, D/T = A _SE/η ^p ) gives an average absolute percentage difference (AAPD) of 7.72% between the experimental and calculated D values where p and A _SE are constants, T is the absolute temperature, and the AAPD is the average value of (10²) (|D _calcd - D _exptl|/D _exptl). Two other MSE fits using subsets of the 35 diffusion constants may be useful for (a) estimating the viscosity of the hydrophobic core of lipid droplets, where squalene is a naturally occurring component, and (b) providing estimates of the D values needed to design extraction processes by which squalene is obtained from plant oils. The Wilke-Chang equation also was considered and found to give larger AAPDs than the corresponding MSE fits.

Reactions and Products of Squalene and Ozone: A Review

This critical review describes the squalene-ozone (SqOz) reaction, or squalene ozonolysis. Ambient ozone penetrates indoors and drives indoor air chemistry. Squalene, a component of human skin oil, contains six carbon-carbon double bonds and is very reactive with ozone. Bioeffluents from people contribute to indoor air chemistry and affect the indoor air quality, resulting in exposures because people spend the majority of their time indoors. The SqOz reaction proceeds through various formation pathways and produces compounds that include aldehydes, ketones, carboxylic acids, and dicarbonyl species, which have a range of volatilities. In this critical review of SqOz chemistry, information on the mechanism of reaction, reaction probability, rate constants, and reaction kinetics are compiled. Characterizations of SqOz reaction products have been done in laboratory experiments and real-world settings. The effect of multiple environmental parameters (ozone concentration, air exchange rate (AER), temperature, and relative humidity (RH)) in indoor settings are summarized. This critical review concludes by identifying the paucity of available exposure, health, and toxicological data for known reaction products. Key knowledge gaps about SqOz reactions leading to indoor exposures and adverse health outcomes are provided as well as an outlook on where the field is headed.

Multiphase Ozonolysis of Oleic Acid-Based Lipids: Quantitation of Major Products and Kinetic Multilayer Modeling

Commonly found in atmospheric aerosols, cooking oils, and human sebum, unsaturated lipids rapidly decay upon exposure to ozone, following the Criegee mechanism. Here, the gas-surface ozonolysis of three oleic acid-based compounds was studied in a reactor and indoors. Under dry conditions, quantitative product analyses by ¹H NMR indicate up to 79% molar yield of stable secondary ozonides (SOZs) in oxidized triolein and methyl oleate coatings. Elevated relative humidity (RH) significantly suppresses the SOZ yields, enhancing the formation of condensed-phase aldehydes and volatile C9 products. Along with kinetic parameters informed by molecular dynamics simulations, these results were used as constraints in a kinetic multilayer model (KM-GAP) simulating triolein ozonolysis. Covering a wide range of coating thicknesses and ozone levels, the model predicts a much faster decay near the gas-lipid interface compared to the bulk. Although the dependence of RH on SOZ yields is well predicted, the model overestimates the production of H₂O₂ and aldehydes. With negligible dependence on RH, the product composition for oxidized oleic acid is substantially affected by a competitive reaction between Criegee intermediates (CIs) and carboxylic acids. The resulting α-acyloxyalkyl hydroperoxides (α-AAHPs) have much higher molar yields (29-38%) than SOZs (12-16%). Overall, the ozone-lipid chemistry could affect the indoor environment through "crust" accumulation on surfaces and volatile organic compound (VOC) emission. In the atmosphere, the peroxide formation and changes in particle hygroscopicity may have effects on climate. The related health impacts are also discussed.

Experimental evidence that halogen bonding catalyzes the heterogeneous chlorination of alkenes in submicron liquid droplets

A key challenge in predicting the multiphase chemistry of aerosols and droplets is connecting reaction probabilities, observed in an experiment, with the kinetics of individual elementary steps that control the chemistry that occurs across a gas/liquid interface. Here we report evidence that oxygenated molecules accelerate the heterogeneous reaction rate of chlorine gas with an alkene (squalene, Sqe) in submicron droplets. The effective reaction probability for Sqe is sensitive to both the aerosol composition and gas phase environment. In binary aerosol mixtures with 2-decyl-1-tetradecanol, linoleic acid and oleic acid, Sqe reacts 12-23× more rapidly than in a pure aerosol. In contrast, the reactivity of Sqe is diminished by 3× when mixed with an alkane. Additionally, small oxygenated molecules in the gas phase (water, ethanol, acetone, and acetic acid) accelerate (up to 10×) the heterogeneous chlorination rate of Sqe. The overall reaction mechanism is not altered by the presence of these aerosol and gas phase additives, suggesting instead that they act as catalysts. Since the largest rate acceleration occurs in the presence of oxygenated molecules, we conclude that halogen bonding enhances reactivity by slowing the desorption kinetics of Cl2 at the interface, in a way that is analogous to decreasing temperature. These results highlight the importance of relatively weak interactions in controlling the speed of multiphase reactions important for atmospheric and indoor environments.

Molecular Structure of the Surface-Immobilized Super Uranyl Binding Protein

Recently, a super uranyl binding protein (SUP) was developed, which exhibits excellent sensitivity/selectivity to bind uranyl ions. It can be immobilized onto a surface in sensing devices to detect uranyl ions. Here, sum frequency generation (SFG) vibrational spectroscopy was applied to probe the interfacial structures of surface-immobilized SUP. The collected SFG spectra were compared to the calculated orientation-dependent SUP SFG spectra using a one-excitonic Hamiltonian approach based on the SUP crystal structures to deduce the most likely surface-immobilized SUP orientation(s). Furthermore, discrete molecular dynamics (DMD) simulation was applied to refine the surface-immobilized SUP conformations and orientations. The immobilized SUP structures calculated from DMD simulations confirmed the SUP orientations obtained from SFG data analyzed based on the crystal structures and were then used for a new round of SFG orientation analysis to more accurately determine the interfacial orientations and conformations of immobilized SUP before and after uranyl ion binding, providing an in-depth understanding of molecular interactions between SUP and the surface and the effect of uranyl ion binding on the SUP interfacial structures. We believe that the developed method of combining SFG measurements, DMD simulation, and Hamiltonian data analysis approach is widely applicable to study biomolecules at solid/liquid interfaces.