Polarizable Molecular Dynamics and Experiments of 1,2-Dimethoxyethane Electrolytes with Lithium and Sodium Salts: Structure and Transport Properties.

The structure and transport properties of electrolyte solutions of 1,2-dimethoxyethane (DME) having salts of Li+ with bis(trifluoromethanesulfonyl)imide ([TFSI]-) or Na+ with [TFSI]- are investigated with polarizable molecular dynamics and experiments. Polarizable force fields for Li+ and Na+ with DME and [TFSI]- were developed based on quantum chemistry calculations, ab initio molecular dynamics simulations, and thermodynamic liquid-state properties. Simulation results for density, viscosity, self-diffusion coefficient, and conductivity of the electrolytes all agree well with the trends and magnitudes of available experimental data for a wide range of salt concentrations. As the concentration of salt increases, the electrolytes become more viscous and molecular species become less mobile. Ionic conductivity does not change monotonically with salt concentration and exhibits a maximum between 0.5 and 1.0 M for both Li[TFSI] and Na[TFSI] electrolytes. Comparatively, both cations are solvated by 5-6 DME or [TFSI]- oxygen atoms and exhibit similar diffusivities and conductivities. The solvation shell of Na+ is found to be more weakly bound and to have a lower binding residence time than that of Li+. The transport of Li+ therefore is more vehicular, through the motion of the solvation shell, while the transport of Na+ is based more on exchange, through the replacement of solvating species. The atomistic insight provided by this work can be used as the basis for future rational design of improved electrolyte solvents for lithium-oxygen, sodium-oxygen, and lithium-sulfur batteries.

Ionic-Functionalized Polymers of Intrinsic Microporosity for Gas Separation Applications.

Ionic-functionalized microporous materials are attractive for energy-efficient gas adsorption and separation processes and have shown promising results in gas mixtures at pressure ranges and compositions that are relevant for industrial applications. In this work, we studied the influence of different counterions (Li+, Na+, K+, Rb+, and Mg2+) on the porosity, carbon dioxide (CO2) gas adsorption, and selectivity in ionic-functionalized PIM-1 (IonomIMs), a polymer belonging to the class of linear and amorphous microporous polymers known as polymers of intrinsic microporosity (PIMs). It was found that an increase in the concentration of ionic groups led to a decrease in the free volume, resulting in a less porous polymer framework, and Mg2+-functionalized IonomIMs exhibited a relatively larger porosity compared to other IonomIMs. The CO2 adsorption capacity was affected by the different counterions for IonomIM-1, and a higher loading capacity for pure CO2 was observed for Mg2+. Furthermore, the IonomIMs showed an enhanced CO2 selectivity in CO2/CH4 and CO2/N2 gas mixtures at conditions used in pressure swing adsorption and vacuum swing adsorption applications. It was also observed that the concentration of ionic groups plays a vital role in changing the CO2 gas adsorption and selectivity.

NLDFT Pore Size Distribution in Amorphous Microporous Materials.

The pore size distribution (PSD) is one of the most important properties when characterizing and designing materials for gas storage and separation applications. Experimentally, one of the current standards for determining microscopic PSD is using indirect molecular adsorption methods such as nonlocal density functional theory (NLDFT) and N2 isotherms at 77 K. Because determining the PSD from NLDFT is an indirect method, the validation can be a nontrivial task for amorphous microporous materials. This is especially crucial since this method is known to produce artifacts. In this work, the accuracy of NLDFT PSD was compared against the exact geometric PSD for 11 different simulated amorphous microporous materials. The geometric surface area and micropore volumes of these materials were between 5 and 1698 m2/g and 0.039 and 0.55 cm3/g, respectively. N2 isotherms at 77 K were constructed using Gibbs ensemble Monte Carlo (GEMC) simulations. Our results show that the discrepancies between NLDFT and geometric PSD are significant. NLDFT PSD produced several artificial gaps and peaks that were further confirmed by the coordinates of inserted particles of a specific size. We found that dominant peaks from NLDFT typically reported in the literature do not necessarily represent the truly dominant pore size within the system. The confirmation provides concrete evidence for artifacts that arise from the NLDFT method. Furthermore, a sensitivity analysis was performed to show the high dependency of PSD as a function of the regularization parameter, λ. A higher value of λ produced a broader and smoother PSD that closely resembles geometric PSD. As an alternative, a new criterion for choosing λ, called here the smooth-shift method (SSNLDFT), is proposed that tuned the NLDFT PSD to better match the true geometric PSD. Using the geometric pore size distribution as our reference, the smooth-shift method reduced the root-mean-square deviation by ∼70% when the geometric surface area of the material is greater than 100 m2/g.

Polymer ultrapermeability from the inefficient packing of 2D chains.

The promise of ultrapermeable polymers, such as poly(trimethylsilylpropyne) (PTMSP), for reducing the size and increasing the efficiency of membranes for gas separations remains unfulfilled due to their poor selectivity. We report an ultrapermeable polymer of intrinsic microporosity (PIM-TMN-Trip) that is substantially more selective than PTMSP. From molecular simulations and experimental measurement we find that the inefficient packing of the two-dimensional (2D) chains of PIM-TMN-Trip generates a high concentration of both small (<0.7 nm) and large (0.7-1.0 nm) micropores, the former enhancing selectivity and the latter permeability. Gas permeability data for PIM-TMN-Trip surpass the 2008 Robeson upper bounds for O2/N2, H2/N2, CO2/N2, H2/CH4 and CO2/CH4, with the potential for biogas purification and carbon capture demonstrated for relevant gas mixtures. Comparisons between PIM-TMN-Trip and structurally similar polymers with three-dimensional (3D) contorted chains confirm that its additional intrinsic microporosity is generated from the awkward packing of its 2D polymer chains in a 3D amorphous solid. This strategy of shape-directed packing of chains of microporous polymers may be applied to other rigid polymers for gas separations.

Ultrathin Molecular-Layer-by-Layer Polyamide Membranes: Insights from Atomistic Molecular Simulations.

In this study, we present an atomistic simulation study of several physicochemical properties of polyamide (PA) membranes formed from interfacial polymerization or from a molecular-layer-by-layer (mLbL) on a silicon wafer. These membranes are composed of meta-phenylenediamine (MPD) and benzene-1,3,5-tricarboxylic acid chloride (TMC) for potential reverse osmosis (RO) applications. The mLbL membrane generation procedure and the force field models were validated, by comparison with available experimental data, for hydrated density, membrane swelling, and pore size distributions of PA membranes formed by interfacial polymerization. Physicochemical properties such as density, free volume, thickness, the degree of cross-linking, atomic compositions, and average molecular orientation (which is relevant for the mLbL membranes) are compared for these different processes. The mLbL membranes are investigated systematically with respect to TMC monomer growth rate per substrate surface area, MPD/TMC ratio, and the number of mLbL deposition cycles. Atomistic simulations show that the mLbL deposition generates membranes with a constant film growth if both the TMC monomer growth rate and MPD/TMC monomer ratio are kept constant. The film growth rate increases with TMC monomer growth rate or MPD/TMC ratio. Furthermore, it was found on one hand that the mLbL membrane density and free volume varies significantly with respect to the TMC monomer growth rate, while on the other hand the degree of cross-linking and the atomic composition varies considerably with the MPD/TMC ratio. Additionally, it was found that both TMC and MPD orient at a tilted angle with respect to the substrate surface, where their angular distribution and average angle orientation depend on both the TMC growth rate and the number of deposition cycles. This study illustrates that molecular simulations can play a crucial role in the understanding of structural properties that can empower the design of the next generation RO membranes created from molecular-layer-by-layer (mLbL) on a silicon wafer.

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.

Liquid-liquid equilibria for soft-repulsive particles: improved equation of state and methodology for representing molecules of different sizes and chemistry in dissipative particle dynamics.

Three developments are presented that significantly expand the applicability of dissipative particle dynamics (DPD) simulations for symmetric and non-symmetric mixtures, where the former contain particles with equal repulsive parameter for self-interactions but a different repulsive parameter for cross-interactions, and the latter contain particles with different repulsive parameters also for the self-interactions. Monte Carlo and molecular dynamics simulations for unary phases covering a wide range of repulsive parameters and of densities for single-bead DPD particles point to deficiencies of the Groot and Warren equation of state (GW-EOS) [J. Chem. Phys. 107, 4423 (1997)]. A revised version, called rGW-EOS, is proposed here that is significantly more accurate over a wider range of parameters/densities. The second development is the generalization of the relationship between the Flory-Huggins χ parameter and the repulsive cross-interaction parameter when the two particles involved have different molecular volumes. The third aspect is an investigation of Gibbs ensemble Monte Carlo simulation protocols, which demonstrates the importance of volume fluctuations and excess volumes of mixing even for equimolar symmetric mixtures of DPD particles. As an illustrative example, the novel DPD methodology is applied to the prediction of the liquid-liquid equilibria for acetic anhydride/(n-hexane or n-octane) binary mixtures.

Bubble bursting as an aerosol generation mechanism during an oil spill in the deep-sea environment: molecular dynamics simulations of oil alkanes and dispersants in atmospheric air/salt water interfaces.

Potential of mean force (PMF) calculations and molecular dynamics (MD) simulations were performed to investigate the properties of oil n-alkanes [i.e., n-pentadecane (C15), n-icosane (C20) and n-triacontane (C30)], as well as several surfactant species [i.e., the standard anionic surfactant sodium dodecyl sulfate (SDS), and three model dispersants similar to the Tween and Span species present in Corexit 9500A] at air/salt water interfaces. This study was motivated by the 2010 Deepwater Horizon (DWH) oil spill, and our simulation results show that, from the thermodynamic point of view, the n-alkanes and the model dispersants have a strong preference to remain at the air/salt water interface, as indicated by the presence of deep free energy minima at these interfaces. The free energy minimum of these n-alkanes becomes deeper as their chain length increases, and as the concentration of surfactant species at the interface increases. The n-alkanes tend to adopt a flat orientation and form aggregates at the bare air/salt water interface. When this interface is coated with surfactants, the n-alkanes tend to adopt more tilted orientations with respect to the vector normal to the interface. These simulation results are consistent with the experimental findings reported in the accompanying paper [Ehrenhauser et al., Environ. Sci.: Processes Impacts 2013, in press, (DOI: 10.1039/c3em00390f)]. The fact that these long-chain n-alkanes show a strong thermodynamic preference to remain at the air/salt water interfaces, especially if these interfaces are coated with surfactants, makes these species very likely to adsorb at the surface of bubbles or droplets and be ejected to the atmosphere by sea surface processes such as whitecaps (breaking waves) and bubble bursting. Finally, the experimental finding that more oil hydrocarbons are ejected when Corexit 9500A is present in the system is consistent with the deeper free energy minima observed for the n-alkanes at the air/salt water interface at increasing concentrations of surfactant species.

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.

Molecular simulations of green leaf volatiles and atmospheric oxidants on air/water interfaces.

Green leaf volatiles (GLVs) are oxygenated hydrocarbons that are emitted by plants, especially under stress conditions such as mechanical damage and local weather changes. GLVs can react with photochemically-generated oxidants (e.g., OH radicals) in atmospheric water drops, and contribute to the formation of secondary organic aerosols (SOAs). Here we investigated the adsorption of a gas phase GLV, 2-methyl-3-buten-2-ol (MBO) and OH radicals on atmospheric air/water interfaces using classical molecular dynamics (MD) simulations and potential of mean force (PMF) calculations. Our models can reproduce experimental values of the free energy of hydration of MBO and ˙OH, as well as 1-octanol/water partition coefficients of MBO determined experimentally in this study. Both MBO and ˙OH have a strong thermodynamic incentive to remain at the air/water interface, with their density profiles overlapping significantly at the interface. These results suggest that chemical reactions between MBO and ˙OH are more likely to take place at the interface, rather than inside the bulk of water droplets or in the vapor phase. We found a significant number of contacts between MBO and ˙OH in our simulations, which could lead to reactions between these two species.

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.

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.

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.