Molecular Dynamics Simulations of Liquid and Polymer Electrolytes for Energy Storage Devices

Advancing beyond current lithium-ion technology is necessary in order to enable energy storage devices for electric airplanes. Electrolyte stability is a key limiting factor, yet the design of improved electrolytes remains a formidable challenge. Molecular dynamics (MD) simulations are a powerful tool for studying electrolytes, since they can be used to evaluate structural, thermodynamic, and transport properties, and can provide molecular-level detail often inaccessible to experimental techniques. Our computational materials groups at the NASA Ames Research Center has developed models and methods to accurately simulate both liquid and polymer electrolytes.We report the results from atomistic MD simulations of several electrolyte materials, with lithium salts dissolved in ionic liquids, dimethoxyethane (DME), and polyethylene oxide (PEO). For improved accuracy, we employ polarizable models, where each atom is given an environment-dependent atomic dipole. The simulations accurately predict bulk transport properties, including viscosity, diffusion, and ionic conductivity, in quantitative agreement with available experimental data. Moreover, the simulations provide important insights into the solvation structure of the lithium ions.We also report the results from coarse-grained MD simulations of polyanion electrolytes. In order to more efficiently capture the longer length- and time-scales of these systems, we employ a generic bead-spring model. These simulations provide important insight into how the polymer chain architecture and ionic interaction strengths affect the ionic aggregation behavior and cation dynamics. Despite the simplicity of the model, the simulations yield qualitative agreement with experimental data for similar systems

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

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

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 N₂ 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 m²/g and 0.039 and 0.55 cm³/g, respectively. N₂ 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 m²/g

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

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

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