Search NASA⌕ Search

SEARCH · Search NASA

Results for “materials and molecular modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
NASA NTRS1986

Brown dwarfs and Jovian planets: A comparison

The recent detection of a subluminous companion to the M dwarf star VB8 has renewed interest in the characteristics of objects spanning the mass range from Jupiter to hydrogen burning stars. Atmospheric and interior models were constructed for objects in this mass regime, up to 30 Jupiter masses, with emphasis on understanding the relationship of brown dwarfs such as the VB8 companion to the better-studied Jovian planets. The atmospheric model solves the equation of radiative transfer assuming frequency dependent molecular opacity sources H2, He, H2O, CO, and CH4 which are important by virtue of the high cosmic abundance of their constituent atoms. Condensation of cosmochemically important materials, iron and silicates, in the atmosphere is possible, and the effect of such grains as opacity sources is assessed. The luminosity of the object is presumed due to degenerate cooling following a collapse phase and possibly deuterium burning and an interior model is constructed using as an outer boundary condition the temperature and pressure level at which the atmosphere becomes convective. The interior model is analogous to Jupiter, with a large liquid metallic-hydrogen core and a thinner molecular-hydrogen envelope. The oxidation state of carbon in the outer envelope of a brown dwarf of similar age to Jupiter is a function of the object's mass. This makes the wavelength dependence of the atmospheric opacity sensitive to the carbon to oxygen ratio, since the abundance of the primary source of molecular opacity, H2O, decreases as more oxygen is tied up as CO.

Lunine, J. I.↗
DOE OSTI2023

Sensitivity of Dislocation-GB interactions to simulation setups in atomistic models

Dislocation-grain boundary (GB) interactions play an integral role in strengthening of crystalline materials, and are dependent on external loading conditions and atomic arrangements at the grain boundary. While molecular dynamics (MD) simulations can provide critical insights into relationships between these parameters and the exact interaction, the outcomes are sensitive to the computational setup. Here, in this work, we explore the effect of computational setup (system size, GB orientation and boundary conditions) on stress-controlled dislocation-grain boundary (DGB) reactions using MD simulations for a well-studied copper symmetric $\langle$1 1 2$\rangle$ tilt boundary. The study demonstrates that the DGB reactions (pinning, absorption, transmission, etc.) are sensitive to the system size and the orientation of the GB to the applied stress. Additionally, the current work shows that there is a critical system size for the dislocation-grain boundary reaction stresses to converge, one which was found to be much larger than the setups in comparable MD studies. This was attributed to the specific stress states and the influence of surfaces, fixed or free, as these computational setups are modified. Finally, a stress-controlled setup with minimal model artifacts is examined, and the effect of coupled Schmid and non-Schmid stress components on the dislocation-GB reactions are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Machine learning prediction of glass transition temperature of conjugated polymers from chemical structure

Predicting the glass transition temperature (T g ) is of critical importance as it governs the thermomechanical performance of conjugated polymers (CPs). Here, we report a predictive modeling framework to predict T g of CPs through the integration of machine learning (ML), molecular dynamics (MD) simulations, and experiments. With 154 T g data collected, an ML model is developed by taking simplified “geometry” of six chemical building blocks as molecular features, where side-chain fraction, isolated rings, fused rings, and bridged rings features are identified as the dominant ones for T g . MD simulations further unravel the fundamental roles of those chemical building blocks in dynamical heterogeneity and local mobility of CPs at a molecular level. The developed ML model is demonstrated for its capability of predicting T g of several new high-performance solar cell materials to a good approximation. The established predictive framework facilitates the design and prediction of T g of complex CPs, paving the way for addressing device stability issues that have hampered the field from developing stable organic electronics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Closing the Gap Between Modeling and Experiments in the Self-assembly of Biomolecules at Interfaces and in Solution

Molecular self-assembly is a powerful tool in materials design, wherein non-covalent interactions like electrostatic, hydrophobic, hydrogen bonding, and van der Waals can be exploited to produce supramolecular nanostructures that are functional and highly tunable. Biomolecules are attractive building blocks, as they are biocompatible, biodegradable and adopt a wide array of higher order structures. Moreover, naturally occurring protein systems display a manifold of structures and interactions that can be replicated in synthetic biomolecules. In this perspective, we highlight advances in multiscale simulation techniques across broad spatiotemporal scales that can aid in characterizing self-assembly of hybrid and hierarchical bionanomaterial systems, with an emphasis on physics-based simulation approaches currently employed to study biomolecules at mineral interfaces. The power of these approaches is highlighted across a few recent areas where molecular simulations have advanced our understanding of self-assembly spanning peptides to protein self-assembly. Looking forward, we discuss how in the near future emerging methods in statistical and machine learning will advance this research field in all areas from expanding the capabilities of physics-based simulation methods to enabling new analyses of high throughput experiments. These advances will pave the way for understanding the molecular recognition patterns in systems that are dictated by self-assembly - biomineralizing peptides, hierarchical peptoids, and large protein assemblies, and will aid in the development of a new synthesis science for achieving precise molecular control in materials design

Sampath, Janani↗
DOE OSTI2023

Response of Sulfonated Polystyrene Melts to Nonlinear Elongation Flows

Ionizable polymers form dynamic networks with domains controlled by two distinct energy scales, ionic interactions and van der Waals forces; both evolve under elongational flows during their processing into viable materials. A molecular level insight of their nonlinear response, paramount to controlling their structure, is attained by fully atomistic molecular dynamics simulations of a model ionizable polymer, polystyrene sulfonate. As a function of increasing elongational flow rate, the systems display an initial elastic response, followed by an ionic fraction-dependent strain hardening, stress overshoot, and eventually strain-thinning. As the sulfonation fraction increases, the chain elongation becomes more heterogeneous. Finally, flow-driven ionic assembly dynamics that continuously break and reform control the response of the system.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Development of neural network force fields for corrosion studies

To fully understand the chemistry and physics of corrosion, novel methods of simulation must be developed. One approach is designing machine learning (ML) algorithms integrated with density functional theory to develop adaptive force fields to gain insight into corrosion behavior namely at the surface of metal oxides. Current methods of modeling corrosion are slow due to the computational cost of resolving both reaction mechanics and mass transport processes. Machine learning methods can be implemented to obtain structure-activity relationships at both the molecular and bulk scale while still retaining the accuracy of density functional theory (DFT) and significantly decreasing the time needed for simulations of complex chemical processes in the various environments of corrosion. Multiscale models are needed for corrosion studies to fully understand its processes not only at the atomic length scale (chemical bonding, energies, and forces), but also at the nano and meso length scales (solid-state physics and material science processes). Current methods of study include DFT, molecular dynamics, and Monte Carlo. The limitation of DFT is that only a small number of atoms or molecules can be simulated at that level of theory. Density functional theory is used to study the electronic structure of atoms and molecules, and calculate the force component of each atom. However, these calculations are limited to about 1000 atoms. Custom periodic boundary conditions (PBC) can be used to describe the various environments and defects that affect the atomic forces to produce a large data set from which a training set can be derived. Machine learning can be utilized to overcome the barrier of modeling macroscopic and multi-scale processes from ab initio calculations through the development of adaptive force fields. Local environments determine the atomic forces of a given system, therefore adaptive force fields must be created to produce reliable quantum mechanical calculations. This can be achieved by developing a learning algorithm that uses the mapped atomic forces or fingerprint as an input to produce energies and magnetic moments as output. A systematic approach was used to begin to build a data set in order to accurately describe the atomic forces in various environments. In Figure 4 below, a simple PBC cell of Fe{sub 2}O{sub 3} was first optimized. A surface optimization was performed next, followed by a hydroxylated surface optimization. Once this calculation has converged, the adsorption of halide species to the hydroxylated surface will be investigated. TensorFlow is an open source platform for machine learning developed by Google. Using a high level application program interface (API) such as Keras allows for building and training ML models easily in a number of different environments and languages. For this project, a neural network was developed within Anaconda in Python. Future Work: Further development of reference data set; Refining neural network and learning algorithm; Fingerprinting atomic environment to enable mapping of atomic force components; Choosing appropriate training set from reference data; Learning from training set and enabling non-linear mapping of training set fingerprints and the atomic forces; Estimation of uncertainty to identify ranges of outside applicability; Testing and analysis of molecular dynamic simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Investigation of fracture in porous materials: a phase-field fracture study informed by ReaxFF

Microscopic features (e.g, pore shapes, sizes, and distribution) in porous material substantially affect the overall mechanical properties such as stiffness and strength. In turn, these material properties determine the macroscopic behaviors of fracture in the porous material. In certain cases, macroscopic properties can be derived from the porous skeleton and void ratio (i.e., porosity), but in many other cases, such derivation is a challenging task. This paper presents a numerical investigation of microporosity and micropore shapes effect on the macrofracture behavior in porous amorphous silica. For this study, we extend the recently-proposed combined molecular dynamic (MD) and phase-field (PF) fracture modeling approach by including different pore shapes in the atomistic domain. In the MD simulations, we adopt ReaxFF to evaluate the material properties, where four different micropore cases are considered. Based on the material properties derived from MD simulations, the macrofracture propagation of porous media is studied using hybrid PF simulation. In the characterization of the pore structure, the concept of pore ligament is proposed to relate the pore shape and the critical energy release rate. Two classical fracture problems were used to evaluate the effect of pore shape on the macrofracture behavior. The results of the case study show that although the micropore shapes change the macrofracture behaviors, these effects vary with the geometry and loading conditions of macroscopic boundary value problems. Furthermore, the case study also shows that the influence of micropore structure can be captured at the macroscopic level through the material properties derived from the MD simulations.

42 ENGINEERING↗
DOE OSTI2021

Understanding geopolymer binder-aggregate interfacial characteristics at molecular level

The interfacial characteristics of geopolymer binder to aggregate composites are poorly understood, especially at molecular level. Herein, molecular models are developed to study, for the first time, the geopolymer-aggregate interface. Chemically, various forms of interfacial bonding are characterized, including Al-O-Si bonding through condensation reactions, NaO and H-bonding. An atomic-level interfacial transition zone (ITZ) is identified, attributed to the concentration of –OH groups. Increasing the Si/Al ratio of geopolymer is found to decrease the ITZ density, but have limited effect on the ITZ width. A heterogeneous diffusion characteristic occurs in geopolymer, due to the weak interfacial interaction. Mechanically, lowering the Si/Al ratio promotes the interfacial strength due to the stronger interfacial interaction and higher cross-linking degree in geopolymer. Under loading the interfacial fracture undergoes three stages: crack propagation, chain bridging (including aluminosilicate and ionic bridging) and breakage. The above atomic-level findings may facilitate a better design of geopolymer concrete in engineering.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Systems-To-Atoms (S2A): enabling hydrogen for climate security

The project addresses a critical gap in hydrogen infrastructure by integrating system-level energy models with atomic-scale material simulations in a unified Systems-to-Atoms (S2A) framework. The motivation stems from the need to develop efficient, cost-effective, and durable hydrogen transport and utilization technologies to support decarbonization of hard-to-electrify sectors such as heavy-duty transportation. Current system models lack awareness of material performance mechanisms, while material-scale models do not account for system-level usage and variability. To bridge this divide, the team developed a co-simulation capability linking techno-economic analyses, reactor/process-flow modeling, and molecular-scale catalysis simulations. Applied to hydrogen delivery in California, the framework enabled comparative evaluations of compressed, cryogenic, and liquid organic hydrogen carrier (LOHC) pathways, highlighting how catalyst operation and unit process efficiency influence overall performance. The results demonstrate that no single material or transport mode is universally optimal; instead, heterogeneous solutions tuned to specific operational contexts deliver better performance. The project delivers a new capability for cross-scale material co-design, advancing hydrogen infrastructure readiness and informing DOE and LLNL missions in climate and energy resilience.

organic↗
DOE OSTI2023

Crystallographic and temperature effects in low-energy collisions for plasma–material interactions

The interaction of plasma with materials is critical for the fundamental understanding of non-equilibrium processes and their wide application. Recent achievements in fusion energy research emphasize the importance of this problem. Because modelling and predicting plasma–material interactions (PMI) require considering a tremendous number of single plasma particle–surface interaction events, this type of modelling was, until recently, possible only within a binary collision approximation (BCA). The BCA approach considers materials as uniform temperature-insensitive media. The research presented here addresses the PMI problem within an atomistic-based approach for the first time at a statistically significant level. Approximately 10 5 molecular dynamics trajectories were generated for 100 eV deuterium ions interacting with a tungsten surface. The research demonstrated the critical importance of incorporating the discrete lattice structure of matter into the model. Deuterium penetration depth and fraction of backscattered deuterium ions strongly depend on the surface orientation, the impact incident ion directions, and the ions’ initial positions. On average, the BCA-based calculations underestimate the penetration depth by a factor of two or more, and the fraction of backscattered atoms is overestimated by a factor of two or more. As a result, increasing the material temperature from 500 to 3000 K reduces penetration depth by about 30% and increases the fraction of backscattered atoms by 250%.

36 MATERIALS SCIENCE↗
NASA NTRS1984

Mass loss from red giants - A simple evolutionary model for NGC 7027

NGC 7027 is a young planetary nebula with the remnants of a red giant circumstellar envelope surrounding the central, ionized region. By comparing the outer molecular envelope with the inner ionized material, it is argued that the mass loss rate has decreased by at least a factor of 3, and more probably by about a factor of 10, during the past 1000 years. From this result, it is argued that the luminosity of the central star has also decreased substantially during the same time, consistent with models for the rapid evolution of stars just after they evolve off the asymptotic giant branch. In this picture, the distance to NGC 7027 is less than 1300 pc. NGC 7027 was the last (and best) example of a star where apparently the momentum in the outflowing mass /M(dot)v/ is considerably greater than the momentum in the radiation field (L/c). With the above description of this object, the evidence is now strong that quite often the mass lost from late-type giants is ultimately driven to infinity by radiation pressure on grains. If M(dot)v is as large as L/c for asymptotic branch stars, then it is expected that the total amount of mass lost during this stage of evolution is of the same magnitude as the initial mass of the star, and therefore this mass loss can profoundly affect the star's ultimate fate.

Jura, M.↗
DOE OSTI2024

Aberration corrected RF flipper for high resolution neutron spectroscopy

Project Summary Company: Adelphi Technology, Inc. Title: Aberration-corrected High Frequency RF Flipper for High-Resolution Neutron Spectroscopy PI: Dr. Jay Theodore Cremer Topic: C55-11 Enhancement of Scattering Instrumentation Technology Used at Pulsed and Continuous Sources Subtopic: d. Other Statement of the problem or situation that is being addressed. The quest to understand heterogeneous and hierarchical materials is gathering momentum, as described in a 2015 report by the Basic Energy Sciences Advisory Committee on Challenges at the Frontiers of Matter and Energy. For the past 40 years a technique called neutron spin echo (NSE) has been used to probe molecular motions in such non-crystalline materials over time scales from 10’s of picoseconds to 100’s of nanoseconds. The method has provided unique information about the dynamics of soft heterogeneous materials, including confirmation of the de Gennes model of polymer reptation and quantitative measurement of bending constants of biologically relevant lipid membranes. However, scientists continue to clamor for even higher resolution than NSE can provide. Biomaterials, polymers, glasses, and artificially nanostructured materials all manifest slow molecular motions because of weak or competing interactions between subunits and are amenable to study with neutrons, provided sufficiently long dynamical correlation times can be achieved. All these materials have important applications to advanced technologies so understanding them is key to technological progress. General statement of how this problem is being addressed. We will address the need for high-resolution neutron spectroscopy by using a technique called Neutron Resonance Spin Echo (NRSE). While similar to NSE in many respects, this method has the potential to exceed the NSE capabilities, if 2 technical hurdles can be overcome. The major impediments to successful high-resolution NRSE are the availability of two technologies: a very high frequency, efficient, radiofrequency (rf) flipper for neutrons and a method to correct certain magnetic aberrations. Based on previous STTR support and follow-on research we have developed a suitable rf flipper and we have invented a method to correct the magnetic aberrations. Both technologies need refinement to make them suitable for implementation at a neutron source such as the Oak Ridge National Laboratory nuclear reactor. In this proposal we seek to perfect the two technologies and to combine them into a single, operationally convenient device. Commercial Applications and Other Benefits In view of the increasing demand for the unique scientific information that high resolution neutron spectroscopy can provide, we expect several major instrumentation upgrades at both U.S. and foreign neutron centers will require make use of the NRSE method over the coming decade, creating a market for the devices we will design. These components will enhance scientists’ abilities to probe the time dependence of density fluctuations in a wide range of hierarchical and heterogeneous materials many of which are vital to existing and future technologies. Key Words – Polarized Neutrons, Neutron Spin Echo, Neutron Scattering, advanced materials. Summary for Members of Congress Neutron beams are a powerful materials-science probe that provide unique information about the structure of matter. The proposed devices will accelerate scientific discoveries required to achieve national goals for new technological materials.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Ab Initio Molecular Dynamics Simulations of Amorphous Calcium Carbonate: Interpretation of Pair Distribution Function and X-Ray Absorption Spectroscopy Data

The structure and transformation of hydrous amorphous calcium carbonate (ACC) are key to understanding biomineralization pathways and their relationship with the properties of the resulting material. Quantitative interpretation of scattering experiments aimed at elucidating the structure of ACC is challenging, due to the amorphous nature of this material, and, therefore, requires models for the structure and the scattering physics. Here, we generate physically realistic ensembles of hydrated ACC structures and their vibrational disorder from ab initio molecular dynamics (AIMD) simulations with an emphasis on enabling the consistent interpretation of the finer details of three complementary structural probes: neutron and x-ray pair distribution experiments and x-ray absorption spectroscopy (XAS). In each case, we simulate the signal directly in reciprocal space and then manipulate it into the real-space pair distribution function (PDF) or spectrum using the same procedures for the experimental and theoretical data. Good agreement with experiment was obtained across the three techniques with the simulations accounting well for all features in the spectra. Remaining small discrepancies pointed to differences between real samples and the idealized simulated systems such as deviations from the nominal CaCO3·nH2O stoichiometry. Additionally, the simulations offered a more accurate description of the local coordination environment of calcium than previous shell-by-shell fits to spectra of synthetic ACC and classical molecular dynamics simulations. This work demonstrates that AIMD is a powerful approach for extracting detailed structural information from neutron PDF, x-ray PDF, and XAS of amorphous carbonate phases.

Prange, Micah P.↗
DOE OSTI2022

Computational determination of a primary diffusion mode in γ U-10Mo under irradiation

Low enriched uranium (< 20 % 235 U)-molybdenum (U-Mo) monolithic fuel is the primary candidate for high performance research and test reactors, and is in the process of being qualified to replace highly enriched uranium (≥ 20 % 23 5U) fuel. As part of the qualification process, it is critical to understand and predict the behavior of fission gas bubbles under irradiation, which affects fuel swelling and fuel failure. Mechanistic fuel models are being developed that can both reproduce the existing experimental data for fuel swelling, and be further applied to irradiation conditions beyond the experimental scope. Diffusion of species under irradiation conditions is an important parameter in the mechanistic fuel models; however, no temperature-relevant experimental diffusion data exists. In the present work, radiation-enhanced diffusion coefficients of U, Mo, and Xe in λU-10wt.%Mo were calculated in the temperature range between 300 K and 1400 K via rate-theory models and molecular dynamics simulations with an embedded-atom method interatomic potential for the U-Mo-Xe system. Accordingly, total diffusion coefficients under relevant irradiation conditions are determined using previously obtained intrinsic thermal diffusion and radiation-driven diffusion coefficients, as well as the newly calculated radiation-enhanced diffusion coefficients presented herein. Radiation-enhanced diffusion of U and Mo was dominant in the intermediate temperature range, whereas radiation-enhanced diffusion of Xe did not significantly contribute to total diffusion of Xe at the relevant fission rate densities. Radiation-enhanced diffusion of Xe became faster than both intrinsic thermal diffusion and radiation-driven diffusion at a fission rate density of 5 x 10 22 fissions/m 3 /s, which is higher than the typical fission rate density range in research reactors. The temperature regime that radiation-enhanced diffusion of each element dominated was dependent on the fission rate density. Finally, the total diffusion coefficients of U, Mo, and Xe, updated in this work, will be utilized as parameters in the mechanistic fuel models to help predict the behavior of fission gas bubbles under irradiation more accurately.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Self-Healing Ceramics and Mechanical Metamaterials

We performed large scale molecular dynamics (MD) simulations and reactive & quantum molecular dynamics simulations to investigate: (1) faceting, grain growth, and crack healing in ceramics; (2) stress-release mechanism and healing in fibers; and (3) mechanical properties of two-dimensional metamaterials. We also designed highly accurate interatomic potentials and forces using artificial neural networks (ANN) and first-principles molecular dynamics simulation data for use in our MD and reactive MD simulations. We demonstrated the efficacy of modeling quantum materials atomistic simulation data using unsupervised deep learning approaches on a quantum computer. These achievements are elaborated in the following sections.

36 MATERIALS SCIENCE↗
DOE OSTI2026

LDRD FY25 Program Overview

As Lawrence Livermore National Laboratory’s (LLNL’s) Laboratory Directed Research and Development (LDRD) program enters its fifth decade of leading-edge research and development, its impact and importance have never been stronger. The program continues to advance strategic investments in pioneering science, technology, and engineering, ensuring LLNL will be ready to deliver on our mission as it evolves over the coming decades. Investing in LDRD research, and the people who perform this critical work, gives LLNL the ability to sustain our role as a leader in the Department of Energy and National Nuclear Security Administration enterprise. The LDRD program enables high-risk, high-payoff research that anticipates emerging threats and future mission needs. By nurturing the ingenuity of the Lab’s greatest asset, its people, LDRD funding advances not only our research but also grows and nurtures our workforce: engaging future innovators with student mentoring, challenging postdoctoral researchers to apply their skills to support national security, and strengthening the leadership skills of early career staff. This annual report documents how LDRD investments advance LLNL’s science, technology, and engineering across our mission space. To assess LDRD’s impact we track both short and long-term metrics such as peer-reviewed publications, number of students, or professional fellows. In addition to reviewing these metrics, I encourage you to delve deeper into the breadth of science and technology that illustrate the strategic value of this research portfolio. For instance, a recent exploratory research project used advanced manufacturing to construct miniaturized three-dimensional ion traps for a quantum computer with reduced quantum error rates to enable applications that address national security missions and support basic science. Another project has delved into studying detonation by examining deflagration to enhance the safety and security of the nuclear weapons stockpile. LDRD researchers are also deploying AI agents on two of the world’s most powerful supercomputers to automate and accelerate inertial confinement fusion experiments. Other teams are delivering more accurate optical constants to enable improved validation for aluminum to advance atomic and molecular physics models. LDRD-driven discoveries of how metals deform under extreme conditions strengthen our ability to model and design materials for demanding national security environments. National security challenges are increasingly complex and continuously evolving. LDRD focuses our most innovative science and technology on these challenges, ensuring the Laboratory is developing creative, forward-leaning solutions for our nation and the world. The following pages feature highlights of published scientific advances, patents, and honors that stem from LDRD investments. As you read this report, I hope you will understand how these investments position the Laboratory, and our partners, to meet the demands of the decades ahead.

36 MATERIALS SCIENCE↗
NASA NTRS2004

Cross-Linked Nanotube Materials with Variable Stiffness Tethers

The constitutive properties of a cross-linked single-walled carbon nanotube material are predicted with a multi-scale model. The material is modeled as a transversely isotropic solid using concepts from equivalent-continuum modeling. The elastic constants are determined using molecular dynamics simulation. Some parameters of the molecular force field are determined specifically for the cross-linker from ab initio calculations. A demonstration of how the cross-linked nanotubes may affect the properties of a nanotube/polyimide composite is included using a micromechanical analysis.

Frankland, Sarah-Jane V.↗
NASA NTRS2018

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.

Abbott, Lauren J.↗