The effect of functional groups on the glass transition temperature of atmospheric organic compounds: a molecular dynamics study

Organic compounds constitute a substantial part of atmospheric particulate matter not only in terms of mass concentration but also in terms of distinct functional groups. The glass transition temperature provides an indirect way to investigate the phase state of the organic compounds, playing a crucial role in understanding their behavior and influence on aerosol processes. Molecular dynamics (MD) simulations were implemented here to predict the glass transition temperature (Tg) of atmospherically relevant organic compounds as well as the influence of their functional groups and length of their carbon chain. The cooling step used in the simulations was chosen to be neither too low (to supress crystallization) nor too high (to avoid Tg overprediction). According to the MD simulations, the predicted Tg is sensitive to the functional groups as follows: carboxylic acid (-COOH) > hydroxyl (-OH) and (-COOH) > carbonyls (-CO). Increasing the number of carbon atoms leads to higher Tg for the linearly structured compounds. Linear compounds with lower molecular weight were found to exhibit a lower Tg. No clear correlation between O : C and Tg was observed. The architecture of the carbon chain (linear, or branched, or ring) was also found to impact the glass transition temperature. Compounds containing a non-aromatic carbon ring are characterized by a higher Tg compared to linear and branched ones with the same number of carbon atoms.

Phase Transitions in Organic and Organic/Inorganic Aerosol Particles

The phase state of aerosol particles can impact numerous atmospheric processes, including new particle growth, heterogeneous chemistry, cloud condensation nucleus formation, and ice nucleation. In this article, the phase transitions of inorganic, organic, and organic/inorganic aerosol particles are discussed, with particular focus on liquid-liquid phase separation (LLPS). The physical chemistry that determines whether LLPS occurs, at what relative humidity it occurs, and the resultant particle morphology is explained using both theoretical and experimental methods. The known impacts of LLPS on aerosol processes in the atmosphere are discussed. Finally, potential evidence for LLPS from field and chamber studies is presented. By understanding the physical chemistry of the phase transitions of aerosol particles, we will acquire a better understanding of aerosol processes, which in turn impact human health and climate.

Geometric Analysis of Free and Accessible Volume in Atmospheric Nanoparticles

Computer-generated atomistic microstructures of atmospheric nanoparticles are geometrically analyzed using Delaunay tessellation followed by Monte Carlo integration to compute their free and accessible volume. The nanoparticles studied consist of cis-pinonic acid (a biogenic organic aerosol component), inorganic ions (sulfate and ammonium), and water. Results are presented for the free or unoccupied volume in different domains of the nanoparticles and its dependence on relative humidity and organic content. We also compute the accessible volume to small penetrants such as water molecules. Most of the free volume or volume accessible to a penetrant as large as a water molecule is located in the domains occupied by organics. In contrast, regions dominated by inorganics do not have any cavities with sizes larger than 1 Å. Solid inorganic domains inside the particle are practically impermeable to any small molecule, thereby offering practically infinite resistance to diffusion. A guest molecule can find diffusive channels to wander around within the nanoparticle only through the aqueous and organic-rich domains. The largest pores are observed in nanoparticles with high levels of organic mass and low relative humidity. At high relative humidity, the presence of more water molecules reduces the empty space in the inner domains of the nanoparticle, since areas rich in organic molecules (which are the only ones where appreciable pores are found) are pushed to the outer area of the particle. This, however, should not be expected to affect the diffusive process as transport through the aqueous phase inside the particle will be, by default, fast due to its fluid-like nature.

Inversion model for extracting chemically resolved depth profiles across liquid interfaces of various configurations from XPS data: PROPHESY

PROPHESY, a technique for the reconstruction of surface-depth profiles from X-ray photoelectron spectroscopy data, is introduced. The inversion methodology is based on a Bayesian framework and primal-dual convex optimization. The acquisition model is developed for several geometries representing different sample types: plane (bulk sample), cylinder (liquid microjet) and sphere (droplet). The methodology is tested and characterized with respect to simulated data as a proof of concept. Possible limitations of the method due to uncertainty in the attenuation length of the photo-emitted electron are illustrated.

Revealing the Impacts of Chemical Complexity on Submicrometer Sea Spray Aerosol Morphology

Sea spray aerosol (SSA) ejected through bursting bubbles at the ocean surface is a complex mixture of salts and organic species. Submicrometer SSA particles have long atmospheric lifetimes and play a critical role in the climate system. Composition impacts their ability to form marine clouds, yet their cloud-forming potential is difficult to study due to their small size. Here, we use large-scale molecular dynamics (MD) simulations as a "computational microscope" to provide never-before-seen views of 40 nm model aerosol particles and their molecular morphologies. We investigate how increasing chemical complexity impacts the distribution of organic material throughout individual particles for a range of organic constituents with varying chemical properties. Our simulations show that common organic marine surfactants readily partition between both the surface and interior of the aerosol, indicating that nascent SSA may be more heterogeneous than traditional morphological models suggest. We support our computational observations of SSA surface heterogeneity with Brewster angle microscopy on model interfaces. These observations indicate that increased chemical complexity in submicrometer SSA leads to a reduced surface coverage by marine organics, which may facilitate water uptake in the atmosphere. Our work thus establishes large-scale MD simulations as a novel technique for interrogating aerosols at the single-particle level.

Contribution of Different Molecules and Moieties to the Surface Tension in Aqueous Surfactant Solutions. II: Role of the Size and Charge Sign of the Counterions

Understanding the role of the counterion species in surfactant solutions is a complicated task, made harder by the fact that, experimentally, it is not possible to vary independently bulk and surface quantities. Here, we perform molecular dynamics simulations at constant surface coverage of the liquid/vapor interface of lithium, sodium, potassium, rubidium, and cesium dodecyl sulfate aqueous solutions. We investigate the effect of counterion type and charge sign on the surface tension of the solution, analyzing the contribution of different species and moieties to the lateral pressure profile. The observed trends are qualitatively compatible with the Hofmeister series, with the notable exception of sodium. We point out a possible shortcoming of what is at the moment, in our experience, the most realistic nonpolarizable force field (CHARMM36) that includes the parametrization for the whole series of alkali counterions. In the artificial system where the counterion and surfactant charges are inverted in sign, the counterions become considerably harder. This charge inversion changes considerably the surface tension contributions of the counterions, surfactant headgroups, and water molecules, stressing the key role of the hardness of the counterions in this respect. However, the hydration free energy gain of the counterions, occurring upon charge inversion, is compensated by the concomitant free energy loss of the headgroups and water molecules, leading to a negligible change in the surface tension of the entire system.

Molecular dynamics simulations demonstrate that non-ideal mixing dominates subsaturation organic aerosol hygroscopicity

The microscopic properties that determine hygroscopic behavior are complex. The importance of hygroscopicity to many areas, and particularly atmospheric chemistry, in terms of aerosol growth and cloud nucleation, mandate the need for robust models to understand this behavior. Toward this end, we have employed molecular dynamics simulations to calculate hygroscopicity from atomistic models using free energy perturbation. We find that currently available force fields may not be well-suited to modeling the extreme environments of aerosol particles. Nonetheless, the results illuminate some shortcomings in our current understanding of hygroscopic growth and cloud nucleation. The most widely used model of hygroscopicity, κ-Köhler Theory (κKT), breaks down in the case of deviations from ideal solution behavior and empirical adjustments within the simplified framework cannot account for non-ideal behavior. A revised model that incorporates non-ideal mixing rescues the general framework of κKT and allows us to understand our simulation results as well as the behavior of atmospheric aerosols over the full range of humidity. The revised model shows that non-ideal mixing dominates hygroscopic growth at subsaturation humidity. Thus, a model based on ideal mixing will fail to predict subsaturation growth from cloud condensation nucleus (CCN) activation or vice versa; a single parameter model for hygroscopicity will generally be insufficient to extrapolate across wide ranges of humidity. We argue that in many cases, when data are limited to subsaturation humidity, an empirical model for non-ideal mixing may be more successful than one for ideal mixing.

Seed-Adsorbate Interactions as the Key of Heterogeneous Butanol and Diethylene Glycol Nucleation on Ammonium Bisulfate and Tetramethylammonium Bromide

Condensation particle counter (CPC) instruments are commonly used to detect atmospheric nanoparticles. They operate on the basis of condensing an organic working fluid on the nanoparticle seeds to grow the particles to a detectable size, and at the size of few nanometers, their efficiency depends on how well the working fluid interacts with the seeds under the measurement conditions. This study models the first steps of heterogeneous nucleation of two working fluids commonly used in CPCs (diethylene glycol (DEG) and n-butanol) onto two positively charged seeds, ammonium bisulfate and tetramethylammonium bromide. The nucleation process is modeled on a molecular level using a combination of systematic configurational sampling and density functional theory (DFT). We take into account the conformational flexibility of DEG and n-butanol and determine the key factors that can improve the efficiency of nanoparticle measurements by CPCs. The results show that hydrogen bonding between the seed and the working fluid molecules is central to the adsorption of the first DEG/n-butanol molecules onto the seeds. However, intermolecular hydrogen bonding between the adsorbed molecules can also enhance the nucleation process for the weakly adsorbing vapor molecules. Accordingly, the heterogeneous nucleation probability is higher for working fluid-nanoparticle combinations with a higher potential for hydrogen bonding; in this case, DEG and ammonium bisulfate. Moreover, conformational analysis and methodology evaluations indicate that the consideration of adsorbate conformers and step-wise addition of the vapor molecules to the seeds is not essential for qualitative modeling of heterogeneous nucleation systems, at least for systems where the adsorbate and seed chemical properties are clearly different. This is the first molecular-level modeling study reporting detailed chemical reasons for experimentally observed seed and working fluid preferences in CPCs and reproducing the experimental observations. Our presented approach can be likely used for predicting preferences in similar nucleating systems.

Contemporary Formulation Development for Inhaled Pharmaceuticals

Pulmonary delivery has gained increased interests over the past few decades. For respiratory conditions, targeted drug delivery directly to the site of action can achieve a high local concentration for efficacy with reduced systemic exposure and adverse effects. For systemic conditions, the unique physiology of the lung evolutionarily designed for rapid gaseous exchange presents an entry route for systemic drug delivery. Although the development of inhaled formulations has come a long way over the last few decades, many aspects of it remain to be elucidated. In particular, a reliable and well-understood method for in vitro-in vivo correlations remains to be established. With the rapid and ongoing advancement of technology, there is much potential to better utilise computational methods including different types of modelling and simulation approaches to support inhaled formulation development. This review intends to provide an introduction on some fundamental concepts in pulmonary drug delivery and inhaled formulation development followed by discussions on some challenges and opportunities in the translation of inhaled pharmaceuticals from preclinical studies to clinical development. The review concludes with some recent advancements in modelling and simulation approaches that could play an increasingly important role in modern formulation development of inhaled pharmaceuticals.

Investigating the effect of starch/Fe3O4 nanoparticles on biodesulfurization using molecular dynamic simulation

The application of dibenzothiophene (DBT) as a source of energy leads to air pollution. The key solution to overcome this drawback is desulfurization. Magnetic nanoparticles have shown an excellent performance in the desulfurization of dibenzothiophene. In this study, molecular dynamic (MD) simulation was considered for the first time to gain insight about the molecule interactions in the biodesulfurization (BDS) process of DBT using Rhodococcus erythropolis IGTS8, in the presence and absence of starch/magnetic nanoparticles. According to the MD simulation results, the density of the system in the presence of starch/Fe₃O₄ was ascending while in the absence of these nanoparticles, the density was descending. Starch/magnetic nanoparticles caused more rapid equilibrium state in the biodesulfurization process. The energy diagram showed that magnetic nanoparticles decrease the energy fluctuation and increase the difference of non-bounding energy and potential energy (8 times) compared to (BDS) without nanoparticle, which reflects higher bounded energy in the system using starch/magnetic nanoparticles. The height of RDF peak in the presence of starch/Fe₃O₄ was 4 times more than the RDF peak in the absence of nanoparticle. In addition, the nanoparticles decreased the fluctuations around optimal temperature in BDS up to 5% compared to other state.

Size-Dependent Liquid-Liquid Phase Separation in Atmospherically Relevant Complex Systems

Physical properties of aerosol particles, such as liquid-liquid phase separation (LLPS), have the potential to impact the climate system. Model systems have been shown to have size-dependent LLPS in the submicron regime; however, these systems are an extreme simplification of ambient aerosol, which can include myriad organic compounds. We expand the studies of LLPS in particles consisting of ammonium sulfate and more complex organic mixtures from multiple organic compounds to α-pinene secondary organic matter (SOM). All systems display a size-dependent morphology, with small particles remaining homogeneous while large particles phase-separate. Surprisingly, three-phase particles were also observed in some of the systems in addition to a new phase state that we have termed channel morphology, which can arise upon efflorescence. The existence of size-dependent LLPS in complex organic mixtures and SOM provides evidence that this is a relevant phenomenon for ambient aerosol and should be considered when modeling atmospheric aerosol.

Numerical Study of Soft Colloidal Nanoparticles Interaction in Shear Flow

The mechanical behavior of nanoparticle assemblies depends on complex particle interactions that are difficult to study experimentally. Depending on the nanoparticle morphology, these interactions could lead to adhesive and elastic-plastic behavior during contact deformation. The aim of this research is to study the effect of contact interactions between polymer nanoparticles and their impact on the macroscopic properties of formed aggregates. For this purpose, the discrete element method (DEM) was used to develop an interaction model combining elastic-plastic deformation and adhesion to study the behavior of spherical polymeric nanoparticles. Initially, a pair of particles interacting in the normal direction was simulated to evaluate the effect of adhesion and plastic deformation in the pull-off force of the contact. Based on these results, the simulations were extended to a dispersed system of nanoparticles, in which multibody interactions become dominant. Considering the aggregation between the nanoparticles induced by a shear flow, we performed an analysis of the number of aggregates and aggregates size in time to characterize the strength of clusters formed during the process. The simulation results showed that the interaction strength upon breakage of the clusters, correlating with the aggregates size, depends on the nanoparticle's softness. In this way, we verified that the type of contact interaction directly influences the macroscopic mechanical response of nanoparticle assemblies. Therefore, our model represents a new way of predicting the mechanical behavior of polymer nanoparticle systems and of optimizing it by adjusting primary particle properties.