Thilanga P. Liyana-Arachchi

Molecular Simulation Study of the Adsorption of Naphthalene and Ozone on Atmospheric Air/Ice Interfaces

The adsorption of gas-phase naphthalene and ozone on atmospheric air/ice interfaces was investigated using classical molecular dynamics (MD) simulations and potential of mean force (PMF) calculations. Naphthalene and ozone exhibit a strong preference to be adsorbed at the air/ice interface, rather than being dissolved into the bulk of the quasi-liquid layer (QLL) or incorporated into the ice crystals. When the air/ice interface is coated with increasing concentrations of naphthalene molecules, the QLL becomes thinner and surface adsorption of ozone is enhanced. Furthermore, ozone tends to adsorb on top of the naphthalene film, although significant penetration of ozone into this film is also observed. Naphthalene molecules tend to adopt a flat orientation on the air/ice interface. Less variation in the orientation was observed for lower concentrations of naphthalene, whereas variations in the ozone concentration do not affect the orientation of naphthalene molecules. However, as the concentration of ozone increases, most of the naphthalene molecules still prefer to stay close to the mobile water molecules in the QLL, but a significant fraction of the naphthalene molecules spends a considerable amount of time inside the thicker layer of ozone. We also monitored the number of contacts between naphthalene and ozone at the air/ice interface upon variations in the concentrations of these two species. These contacts were assumed to be proportional to the reaction rate between these two species. When the number of ozone molecules was held constant, the number of contacts showed a linear relationship to the number of naphthalene molecules. However, when the naphthalene concentration was held constant, for all systems we observed a linear relationship at low ozone concentrations and a plateau at high ozone concentrations.

Adsorption of naphthalene and ozone on atmospheric air/ice interfaces coated with surfactants: a molecular simulation study.

The adsorption of gas-phase naphthalene and ozone molecules onto air/ice interfaces coated with different surfactant species (1-octanol, 1-hexadecanol, or 1-octanal) was investigated using classical molecular dynamics (MD) simulations. Naphthalene and ozone exhibit a strong preference to be adsorbed at the surfactant-coated air/ice interfaces, as opposed to either being dissolved into the bulk of the quasi-liquid layer (QLL) or being incorporated into the ice crystals. The QLL becomes thinner when the air/ice interface is coated with surfactant molecules. The adsorption of both naphthalene and ozone onto surfactant-coated air/ice interfaces is enhanced when compared to bare air/ice interface. Both naphthalene and ozone tend to stay dissolved in the surfactant layer and close to the QLL, rather than adsorbing on top of the surfactant molecules and close to the air region of our systems. Surfactants prefer to orient at a tilted angle with respect to the air/ice interface; the angular distribution and the most preferred angle vary depending on the hydrophilic end group, the length of the hydrophobic tail, and the surfactant concentration at the air/ice interface. Naphthalene prefers to have a flat orientation on the surfactant coated air/ice interface, except at high concentrations of 1-hexadecanol at the air/ice interface; the angular distribution of naphthalene depends on the specific surfactant and its concentration at the air/ice interface. The dynamics of naphthalene molecules at the surfactant-coated air/ice interface slow down as compared to those observed at bare air/ice interfaces. The presence of surfactants does not seem to affect the self-association of naphthalene molecules at the air/ice interface, at least for the specific surfactants and the range of concentrations considered in this study.

Adsorption of Gas-Phase Phenanthrene on Atmospheric Water and Ice Films

The adsorption of gas-phase phenanthrene on atmospheric water and ice films was investigated in a flow-tube reactor with a view to understanding the processing of semi-volatile organic compounds by fog and snow. Air-water (ice) interface partition constants were obtained by measuring the mass uptake of phenanthrene vapor on thin water and ice films with varying thickness. Adsorption enthalpies and entropies were obtained from the temperature dependence of the interfacial partition constants. The surface adsorption is the predominant mechanism for the uptake of phenanthrene in water and ice films with small film thickness or at low temperature. The adsorption of phenanthrene to ice resembles that to sub-cooled water and theres no significant difference between the adsorption of phenanthrene to water and that to quasi liquid layer (QLL) if we take into account the uncertainties on the thermodynamic quantities measured. Molecular dynamics simulations of phenanthrene at air/water and air/ice interfaces support these experimental observations. The interfacial air-water and air-ice partition constants of phenanthrene increased greatly in the presence of surface-active substances, indicating that surface active materials effectively enhanced the uptake of organic compounds by atmospheric water and ice films.

Molecular Modeling of the Green Leaf Volatile Methyl Salicylate on Atmospheric Air/Water Interfaces

Methyl salicylate (MeSA) is a green leaf volatile (GLV) compound that is emitted in significant amounts by plants, especially when they are under stress conditions. GLVs can then undergo chemical reactions with atmospheric oxidants, yielding compounds that contribute to the formation of secondary organic aerosols (SOAs). We investigated the adsorption of MeSA on atmospheric air/water interfaces at 298 K using thermodynamic integration (TI), potential of mean force (PMF) calculations, and classical molecular dynamics (MD) simulations. Our molecular models can reproduce experimental results of the 1-octanol/water partition coefficient of MeSA. A deep free energy minimum was found for MeSA at the air/water interface, which is mainly driven by energetic interactions between MeSA and water. At the interface, the oxygenated groups in MeSA tend to point toward the water side of the interface, with the aromatic group of MeSA lying farther away from water. Increases in the concentrations of MeSA lead to reductions in the height of the peaks in the MeSA-MeSA g(r) functions, a slowing down of the dynamics of both MeSA and water at the interface, and a reduction in the interfacial surface tension. Our results indicate that MeSA has a strong thermodynamic preference to remain at the air/water interface, and thus chemical reactions with atmospheric oxidants are more likely to take place at this interface, rather than in the water phase of atmospheric water droplets or in the gas phase.

Ice growth from supercooled aqueous solutions of benzene, naphthalene, and phenanthrene.

Classical molecular dynamics (MD) were performed to investigate the growth of ice from supercooled aqueous solutions of benzene, naphthalene, or phenanthrene. The main objective of this study is to explore the fate of those aromatic molecules after freezing of the supercooled aqueous solutions, i.e., if these molecules become trapped inside the ice lattice or if they are displaced to the QLL or to the interface with air. Ice growth from supercooled aqueous solutions of benzene, naphthalene, or phenanthrene result in the formation of quasi-liquid layers (QLLs) at the air/ice interface that are thicker than those observed when pure supercooled water freezes. Naphthalene and phenanthrene molecules in the supercooled aqueous solutions are displaced to the air/ice interface during the freezing process at both 270 and 260 K; no incorporation of these aromatics into the ice lattice is observed throughout the freezing process. Similar trends were observed during freezing of supercooled aqueous solutions of benzene at 270 K. In contrast, a fraction of the benzene molecules become trapped inside the ice lattice during the freezing process at 260 K, with the rest of the benzene molecules being displaced to the air/ice interface. These results suggest that the size of the aromatic molecule in the supercooled aqueous solution is an important parameter in determining whether these molecules become trapped inside the ice crystals. Finally, we also report potential of mean force (PMF) calculations aimed at studying the adsorption of gas-phase benzene and phenanthrene on atmospheric air/ice interfaces. Our PMF calculations indicate the presence of deep free energy minima for both benzene and phenanthrene at the air/ice interface, with these molecules adopting a flat orientation at the air/ice interface.

Molecular dynamics simulation of single-walled carbon nanotubes inside liquid crystals

Abstract Molecular dynamic simulations of systems of single-walled carbon nanotubes (CNTs) in liquid crystalline solvents were performed, in order to investigate the effect of the molecular structure and phase of the liquid crystal (LC) on the interactions between the CNTs. Three different LC molecules (5CB, 8CB and 5CF) were considered in our study. Our results with 5CB and 8CB suggest that increasing the chain length of the hydrophobic part of the LC molecule by three carbon atoms is insufficient to decrease the tendency for the CNTs to aggregate in the LCs. Additionally, varying the phase of the LC is also insufficient to decrease the aggregation tendency of the CNTs. However, simulations with 5CF (which has fluorine atoms in the head group of the LC molecule) suggest that this LC solvent can decrease the tendency of the CNTs to aggregate. This study is relevant to assist experimentalists with the development of high-quality dispersions of large concentrations of CNTs in the LCs.

Combined Experimental and Molecular Simulation Investigation of the Individual Effects of Corexit Surfactants on the Aerosolization of Oil Spill Matter

We report laboratory aerosolization experiments and classical molecular dynamics (MD) simulations, with the objective of investigating the individual effects of the two Corexit surfactants Span 80 (nonionic) and dioctyl sodium sulfosuccinate (DOSS, ionic), on the aerosolization of oil spill matter to the atmosphere. Our simulation results show that Span 80, DOSS, and the oil alkanes n-pentadecane (C15) and n-triacontane (C30) exhibit deep free energy minima at the air/seawater interface. C15 and C30 exhibit deeper free energy minima at the interface when Span 80 is present, as compared to the situation when DOSS or no surfactants are at the interface. These results suggest that Span 80 makes these oil hydrocarbons more likely to be adsorbed at the surface of seawater droplets and carried out to the atmosphere, relative to DOSS or to the situation where no surfactants are present. These simulation trends are in qualitative agreement with our experimental observations in a bubble-column setup, where larger amounts of oil hydrocarbons are ejected when Span 80 is mixed with oil and injected into the column, as compared to when DOSS is used. Our simulations also indicate that Span 80 has a larger thermodynamic incentive than DOSS to move from the seawater phase and into the air/seawater interface. This observation is also in qualitative agreement with our experimental measurements, which indicate that Span 80 is ejected in larger quantities than DOSS. Our simulations also suggest that DOSS predominantly adopts a perpendicular orientation with respect to the air/seawater interface at a dispersant to oil ratio (DOR) of 1:20, but has a slight preference to lie parallel to the interfaces at a DOR = 1:5; in both cases, DOSS molecules have their tails wide open and stretched. In contrast, Span 80 has a slight preference to align parallel to the interfaces with a coiled conformation at both DOR values.