Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “molecular dynamic (MD) simulation”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Modes of strain accommodation in Cu-Nb multilayered thin film on indentation and cyclic shear

Two-phase layered thin films with a high density of semi-coherent interfaces exhibit excellent mechanical properties and thermal stability. Here, in this study, a magnetron-sputtered Cu-Nb dual-layered thin film (~500 nm for Cu and ~150 nm for Nb) having an amorphous interface between Cu and Nb with a high density of aligned growth twins in Cu is subjected to severe surface deformation. The material is loaded using indentation and cyclic shear under tribological testing. The strain accommodation in the subsurface microstructure after deformation varies based on the local structure and deformation mode. Grain refinement and crack formations in the stressed region of the Nb layer and localized crystallization of the amorphous interface are observed after indentation and scratch testing. Pronounced detwinning of growth twins in the Cu layer under the cyclic shear strain leaves large dislocations sites and loops which are observed both by high-resolution transmission electron microscopy and experiment-guided molecular dynamic (MD) simulations. Our simulations provided insights into understanding the pathway for the detwinning process under cyclic shear loading.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Analysis of the molecular determinants for furin cleavage of the spike protein S1/S2 site in defined strains of the prototype coronavirus murine hepatitis virus (MHV)

We analyzed the spike protein S1/S2 cleavage of selected strains of a prototype coronavirus, mouse hepatitis virus (MHV) by the cellular protease furin, in order to understand the structural requirements underlying the sequence selectivity of the scissile segment. The probability of cleavage of selected MHV strains was first evaluated from furin cleavage scores predicted by the ProP computer software, and then cleavage was measured experimentally with a fluorogenic peptide cleavage assay consisting of S1/S2 peptide mimics and purified furin. We found that in vitro cleavability varied across MHV strains in line with predicted results—but with the notable exception of MHV-A59, which was not cleaved despite a high score predicted for its sequence. Using the known X-Ray structure of furin in complex with a substrate-like inhibitor as an initial structural reference, we carried out molecular dynamics (MD) simulations to learn the modes of binding of the peptides in the furin active site, and the suitability of the complex for initiation of the enzymatic cleavage. We identified the 3D structural requirements of the furin active site configuration that enable bound peptides to undergo cleavage, and the way in which the various strains tested experimentally are fulfilling these requirements. We find that despite some flexibility in the organization of the peptide bound to the active site of the enzyme, the presence of a histidine at P2 of MHV-A59 fails to properly orient the sidechain of His194 of the furin catalytic triad and therefore produces a distortion that renders the peptide/complex structural configuration in the active site incompatible with requirements for cleavage initiation. The Ser/Thr in P1 of MHV-2 and MHV-S has a similar effect of distorting the conformation of the furin active site residues produced by the elimination of the canonical salt-bridge formed by arginine in P1 position. This work informs a study of coronavirus infection and pathogenesis with respect to the function of the viral spike protein, and suggests an important process of viral adaptation and evolution within the spike S1/S2 structural loop.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Chloride Insertion Enhances the Electrochemical Oxidation of Iron Hydroxide Double-Layer Hydroxide into Oxyhydroxide in Alkaline Iron Batteries

Rechargeable alkaline iron batteries that constitute environmentally benign electrolytes and earth-abundant industrial materials are desirable green solutions for large-scale energy storage. As one of the most abundant metal elements in the earth’s crust, iron (Fe) can satisfy nearly all criteria for low-cost and safe battery electrodes. However, challenges in achieving reversible Fe redox impede their extensive implementation in modern energy supply systems. Here, this study revealed that Cl-anion insertion into Fe(OH) 2 layered double hydroxide (LDH) formed a green rust intermediate phase with the formula [Fe 2 2+ Fe 1 3+ (HO – ) 6 ] + [Cl] – , which assisted a high Fe(OH) 2 /FeOOH conversion reaction (64.7%) and improved cycling stability. This new iron redox chemistry was validated by operando X-ray diffraction, electrochemical testing, X-ray absorption spectroscopy (XAS), X-ray photoelectron spectroscopy (XPS) analysis, scanning transmission electron microscopy–energy-dispersive X-ray spectroscopy (STEM-EDS) mapping, and molecular dynamics (MD) simulations. Our study provides new insight into designing LDH materials for high-capacity alkaline iron batteries.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Structure–Activity Relationships in Ether-Functionalized Solid-State Metal–Organic Framework Electrolytes

The structure–property relationships of metal–organic framework (MOF)-based solid-state electrolytes are not well understood. Herein, a systematic investigation of 12 Zr(IV)-based UiO-66 MOFs with varying ether-chain functional groups was carried out to elucidate the critical microscopic interactions that facilitate improved solid-state electrolyte performance. Enhanced sampling molecular dynamics (MD) simulations were employed and revealed a three-tier ion hopping mechanism: linker–linker hopping, linker-counterion hopping, and counterion-counterion hopping. Detailed structural analysis of the MD trajectories revealed that the chemistry and morphology of the linker groups affect the relative stability and population distribution of the electrolyte components, such that crown-ether-based linker groups enhance the probability of extended, low-barrier ion percolation pathways. As a result, we were able to tune the ionic conductivities by rationally manipulating the counterion distributions, linker binding strengths, and the configurational entropy (multivariability of the linkers). The resulting performance of these MOF-based solid-state electrolytes was significantly enhanced, with a methoxy-functionalized framework (UiO-66-L1 100 ) achieving high ionic conductivities of 2.32 × 10 –4 S/cm and 2.07 × 10 –3 S/cm at 30 °C and 90 °C, respectively, an order of magnitude greater than other all-solid-state MOF electrolyte systems. The electrolyte stability was evaluated with LiIn|LPSCl|MOF:LiTFSI|LPSCl|LiIn symmetric cells, showing excellent Li plating/stripping processes for over 2 months.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Confinement-Driven Heterogeneous Benzene Crystallization in Silica Nanopores

Nanoconfinement alters the thermodynamics, dynamics, and kinetics of fluids hosted in nanoscale solid nanopores to an extent that depends on the characteristics of the confining space and the chemistry of the confined fluids. Confinement-induced alterations in the phase behavior of confined energetic fluids under high pressure or low temperature are highly relevant to subsurface and subsea phenomena such as fluid flow in porous media, hydrate formation and dissociation, and gas storage capacity. Although extensive efforts have been directed toward understanding the phase behavior of confined fluids, the role of solid-liquid interfaces in the phase transitions of organic liquids has not been resolved yet. Here, we explore the onset and growth of benzene crystallization confined in 6 nm sized SBA-15 silica nanopores in the temperature range 300-200 K using in situ extended range small-angle and wide-angle X-ray scattering (SAXS/WAXS) measurements and atomistic classical molecular dynamics (MD) simulations. The crystallization onset of confined benzene depresses to 265 K compared to the freezing point of bulk benzene ( ~278 K), followed by the continuous growth of the emerged crystals in the pore space with complete crystallization at 200 K. The orientation of the emerged benzene crystals is dominated by parallel (π-π stacking) and perpendicular (T-shape stacking) orientations along the cylindrical pore radius and pore length, respectively. The onset of benzene crystals occurs heterogeneously on the pore surface and grows continuously toward the pore center. Further, confined benzene undergoes a dynamical crossover from fragile to strong dynamics behavior, inferred from the rotational and translational diffusion. The insights provided by this study have significant implications for the phase transitions of confined organic liquids that are relevant to a wide range of applications in the biological, geological, environmental, and chemical fields.

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

Decoding the Effect of Anion Identity on the Solubility of N -(2-(2-Methoxyethoxy)ethyl)phenothiazine (MEEPT)

Variations in the solubility of redox-active organic molecules (ROM) of interest for nonaqueous redox flow batteries (RFB), especially as the ROM state-of-charge changes during charge-discharge cycling, present significant molecular design challenges. The situation is further complicated as ROM solubility can be regulated by the choice of electrolyte salt and solvent that together with the ROM comprise the catholyte or anolyte (redox electrolyte) formulation, presenting materials design challenges. The ROM N-(2-(2-methoxyethoxy)ethyl)phenothiazine (MEEPT) is a viscous liquid at room temperature and is miscible in several organic solvents, including acetonitrile and propylene carbonate. The MEEPT radical cation (MEEPT +· ) paired with tetrafluoroborate (BF 4 - ) in acetonitrile presents a 0.5 M solubility, a dramatic decrease when compared to the viscous liquid of neutral MEEPT. Here, in this study, we present a joint experimental, regression modeling, and molecular dynamics (MD) simulations investigation to explore MEEPT-X (where X represents the counteranion) salt solubility variability as a function of concentration and counteranion chemistry in acetonitrile. We find a strong dependence of the salt solubility on the counteranion and relate these findings to explicit intermolecular interactions between MEEPT +· and the counteranion in the electrolyte solution.

Perera, Anton S. [Univ. of Kentucky, Lexington, KY↗
DOE OSTI2022

Lithium Solvation and Mobility in Ionic Liquid Electrolytes with Asymmetric Sulfonyl-Cyano Anion

The solvation structure and transport properties of Li + in ionic liquid (IL) electrolytes based on n-methyl-n-butylpyrrolidinium cyano(trifluoromethanesulfonyl)imide [PYR 14 ][CTFSI] and [Li][CTFSI] (0 ≤ x Li ≤ 0.7) were studied by Raman and Nuclear Magnetic Resonance (NMR) diffusometry, and molecular dynamics (MD) simulations. At x Li < 0.3, Li + coordination is dominated by the cyano group. As x Li is increased, free cyano-sites become limited, resulting in increased coordination via the sulfonyl group. Here, the 1:1 mixture of the symmetric anions bis(trifluoromethanesulfonyl)imide ([TFSI]) and dicyanamide ([DCA]) results in similar physical properties as the IL with [CTFSI]. However, anion asymmetry is shown to increase Li-salt solubility and promote Li+ transference. The lifetimes of Li + -cyano coordination for [CTFSI] are calculated to be shorter than those for [DCA], indicating that the competition from the sulfonyl group weakens its solvation with Li + . This resulted in higher Li + transference for the electrolyte with [CTFSI]. In relation to the utility of these electrolytes in energy storage, the Li–LiFePO 4 half cells assembled with IL electrolyte (x Li = 0.3, 0.5, and 0.7) demonstrated a nominal capacity of 140 mAh/g at 0.1C rate and 90 °C where the cell with x Li = 0.7 IL electrolyte demonstrated 61% capacity retention after 100 cycles and superior rate capability owing to increased electrochemical stability.

25 ENERGY STORAGE↗
DOE OSTI2020

Ligand-Binding-Site Refinement to Generate Reliable Holo Protein Structure Conformations from Apo Structures

The first important step in a structure-based virtual screening is the judicious selection of a receptor protein. In cases where the holo protein receptor structure is unavailable, significant reduction in virtual screening performance has been reported. In this work, we present a robust method to generate reliable holo protein structure conformations from apo structures using molecular dynamics (MD) simulation with restraints derived from holo structure binding-site templates. We perform benchmark tests on two different datasets: 40 structures from a directory of useful decoy-enhanced (DUD-E) and 84 structures from the Gunasekaran dataset. Our results show successful refinement of apo binding-site structures toward holo conformations in 82% of the test cases. In addition, virtual screening performance of 40 DUD-E structures is significantly improved using our MD-refined structures as receptors with an average enrichment factor (EF), an EF1% value of 6.2 compared to apo structures with 3.5. Docking of native ligands to the refined structures shows an average ligand root mean square deviation (RMSD) of 1.97 Å (DUD-E dataset and Gunasekaran dataset) relative to ligands in the holo crystal structures, which is comparable to the self-docking (i.e., docking of the native ligand back to its crystal structure receptor) average, 1.34 Å (DUD-E dataset) and 1.36 Å (Gunasekaran dataset). On the other hand, docking to the apo structures yields an average ligand RMSD of 3.65 Å (DUD-E) and 2.90 Å (Gunasekaran). These results indicate that our method is robust and can be useful to improve virtual screening performance of apo structures.

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

Neural Network for Principle of Least Action

The principle of least action is the cornerstone of classical mechanics, theory of relativity, quantum mechanics, and thermodynamics. Here, we describe how a neural network (NN) learns to find the trajectory for a Lennard-Jones (LJ) system that maintains balance in minimizing the Onsager–Machlup (OM) action and maintaining the energy conservation. The phase-space trajectory thus calculated is in excellent agreement with the corresponding results from the “ground-truth” molecular dynamics (MD) simulation. Furthermore, we show that the NN can easily find structural transformation pathways for LJ clusters, for example, the basin-hopping transformation of an LJ 38 from an incomplete Mackay icosahedron to a truncated face-centered cubic octahedron. Unlike MD, the NN computes atomic trajectories over the entire temporal domain in one fell swoop, and the NN time step is a factor of 20 larger than the MD time step. The NN approach to OM action is quite general and can be adapted to model morphometrics in a variety of applications.

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

Identification of Small-Molecule Inhibitors of Fibroblast Growth Factor 23 Signaling via In Silico Hot Spot Prediction and Molecular Docking to α-Klotho

Fibroblast growth factor 23 (FGF23) is a therapeutic target for treating hereditary and acquired hypophosphatemic disorders, such as X-linked hypophosphatemic (XLH) rickets and tumor-induced osteomalacia (TIO), respectively. FGF23-induced hypophosphatemia is mediated by signaling through a ternary complex formed by FGF23, the FGF receptor (FGFR), and α-Klotho. Currently, disorders of excess FGF23 are treated with an FGF23-blocking antibody, burosumab. Small-molecule drugs that disrupt protein/protein interactions necessary for the ternary complex formation offer an alternative to disrupting FGF23 signaling. In this study, the FGF23:α-Klotho interface was targeted to identify small-molecule protein/protein interaction inhibitors since it was computationally predicted to have a large fraction of hot spots and two druggable residues on α-Klotho. We further identified Tyr433 on the KL1 domain of α-Klotho as a promising hot spot and α-Klotho as an appropriate drug-binding target at this interface. Subsequently, we performed in silico docking of ~5.5 million compounds from the ZINC database to the interface region of α-Klotho from the ternary crystal structure. Following docking, 24 and 20 compounds were in the final list based on the lowest binding free energies to α-Klotho and the largest number of contacts with Tyr433, respectively. Herein, five compounds were assessed experimentally by their FGF23-mediated extracellular signal-regulated kinase (ERK) activities in vitro, and two of these reduced activities significantly. Both these compounds were predicted to have favorable binding affinities to α-Klotho but not have a large number of contacts with the hot spot Tyr433. ZINC12409120 was found experimentally to disrupt FGF23:α-Klotho interaction to reduce FGF23-mediated ERK activities by 70% and have a half maximal inhibitory concentration (IC50) of 5.0 ± 0.23 μM. Molecular dynamics (MD) simulations of the ZINC12409120:α-Klotho complex starting from in silico docking poses reveal that the ligand exhibits contacts with residues on the KL1 domain, the KL1–KL2 linker, and the KL2 domain of α-Klotho simultaneously, thereby possibly disrupting the regular function of α-Klotho and impeding FGF23:α-Klotho interaction. ZINC12409120 is a candidate for lead optimization.

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

Ion Solvation and Transport in Narrow Carbon Nanotubes: Effects of Polarizability, Cation-π Interaction, and Confinement

Understanding ion solvation and transport under confinement is critical for a wide range of emerging technologies, including water desalination and energy storage. While molecular dynamics (MD) simulations have been widely used to study the behavior of confined ions, considerable deviations between simulation results depending on the specific treatment of intermolecular interactions remain. In the following, we present a systematic investigation of the structure and dynamics of two representative solutions, that is, KCl and LiCl, confined in narrow carbon nanotubes (CNTs) with a diameter of 1.1 and 1.5 nm, using a combination of first-principles and classical MD simulations. Our simulations show that the inclusion of both polarization and cation-π interactions is essential for the description of ion solvation under confinement, particularly for large ions with weak hydration energies. Beyond the variation in ion solvation, we find that cation-π interactions can significantly influence the transport properties of ions in CNTs, particularly for KCl, where our simulations point to a strong correlation between ion dehydration and diffusion. Finally, our study highlights the complex interplay between nanoconfinement and specific intermolecular interactions that strongly control the solvation and transport properties of ions.

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

Integrated Quantum-Classical Protocol for the Realistic Description of Solvated Multinuclear Mixed-Valence Transition-Metal Complexes and Their Solvatochromic Properties

Linear cyanide-bridged polymetallic complexes, which undergo photoinduced metal-to-metal charge transfer, represent prototypical systems for studying long-range electron-transfer reactions and understanding the role played by specific solute–solvent interactions in modulating the excited-state dynamics. Here, to tackle this problem, while achieving a statistically meaningful description of the solvent and of its relaxation, one needs a computational approach capable of handling large polynuclear transition-metal complexes, both in their ground and excited states, as well as the ability to follow their dynamics in several environments up to nanosecond time scales. Here, we present a mixed quantum classical approach, which combines large-scale molecular dynamics (MD) simulations based on an accurate quantum mechanically derived force field (QMD-FF) and self-consistent QMD polarized point charges, with IR and UV–vis spectral calculations to model the solvation dynamics and optical properties of a cyano-bridged trinuclear mixed-valence compound (trans-[(NC) 5 Fe III (μ-CN)Ru II (pyridine) 4 (μ-NC)Fe III (CN) 5 ] 4– ). We demonstrate the reliability of the QMD-FF/MD approach in sampling the solute conformational space and capturing the local solute–solvent interactions by comparing the results with higher-level quantum mechanics/molecular mechanics (QM/MM) MD reference data. The IR spectra calculated along the classical MD trajectories in different solvents correctly predict the red shift of the CN stretching band in the aprotic medium (acetonitrile) and the subtle differences measured in water and methanol, respectively. By explicitly including the solvent molecules around the cyanide ligands and calculating the thermal averaged absorption spectra using time-dependent density functional theory calculations within the Tamm–Dancoff approximation, the experimental solvatochromic shift is quantitatively reproduced going from water to methanol, while it is overestimated for acetonitrile. This discrepancy can likely be traced back to the lack of important dispersion interactions between the solvent cyano groups and the pyridine substituents in our micro solvation model. The proposed protocol is applied to the ground state in water, methanol, and acetonitrile and can be flexibly generalized to study excited-state nonequilibrium solvation dynamics.

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

Generating Protein Structures for Pathway Discovery Using Deep Learning

Resolving the intricate details of biological phenomena at the molecular level is fundamentally limited by both length- and time scales that can be probed experimentally. Molecular dynamics (MD) simulations at various scales are powerful tools frequently employed to offer valuable biological insights beyond experimental resolution. However, while it is relatively simple to observe long-lived, stable configurations of, for example, proteins, at the required spatial resolution, simulating the more interesting rare transitions between such states often takes orders of magnitude longer than what is feasible even on the largest supercomputers available today. One common aspect of this challenge is pathway discovery, where the start and end states of a scientific phenomenon are known or can be approximated, but the mechanistic details in between are unknown. Here, we propose a representation-learning-based solution that uses interpolation and extrapolation in an abstract representation space to synthesize potential transition states, which are automatically validated using MD simulations. The new simulations of the synthesized transition states are subsequently incorporated into the representation learning, leading to an iterative framework for targeted path sampling. Our approach is demonstrated by recovering the transition of a RAS-RAF protein domain (CRD) from membrane-free to interacting with the membrane using coarse-grain MD simulations.

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

Computational Insights into Malononitrile-Based Carbanions for CO 2 Capture

Although anionic N and O sites have been widely used in chemisorption of CO 2 , carbanions are much less explored for CO 2 capture. Here, in this study, we employ ab initio calculations and quantum mechanical/molecular mechanical (QM/MM) molecular dynamics (MD) simulations to examine the interaction between CO 2 and the malononitrile carbanion, [CH(CN) 2 ] – . We have explored the potential energy surface of CO 2 binding by scanning the C–C distance between CO 2 and the central C site of the carbanion. We find that CO 2 prefers to bind to the nitrile group physically rather than to form a C–C bond via the carboxylation reaction at the sp 2 C site. Moreover, the two –CN groups can attract two CO 2 molecules at equal strength. The presence of an alkali metal ion enhances both physical and chemical interactions of CO 2 with the malononitrile carbanion. QM/MM MD simulations further confirm the preference of physical interaction in the condensed ionic liquid phase with a phosphonium cation. Our findings suggest that ionic liquids based on the malononitrile carbanion may have a high CO 2 solubility for carbon capture.

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

Hydrophilic and Apolar Hydration in Densely Grafted Cationic Brushes and Counterions with Large Mobilities

We employ all-atom molecular dynamics (MD) simulation framework to unravel water microstructure and ion properties for cationic [Poly(2-(methacryloyloxy)ethyl trimethylammonium chloride] (PMETAC) brushes with chloride ions as counterions. First, we identify locally separate water domains (or first hydration shells) each around the {N(CH 3 ) 3 } + and the C=O functional groups of the PMETAC chain and one around the Cl - ion. These first hydration shells around the respective moieties overlap and the extent of the overlap depends on the nature of the species triggering it. Second, despite the overlap, the water molecules in these domains demonstrate disparate properties dictated by the properties of the atoms/groups around which they are located. For example, the presence of the methyl groups make the {N(CH 3 ) 3 } + group trigger apolar hydration as evidenced by the corresponding orientation of the dipole of the water molecules around the {N(CH 3 ) 3 } + moiety. These water molecules around the N(CH 3 ) 3 } + group also have enhanced tetrahedrality as compared to the water molecules constituting the hydration layer around the C=O group and the Cl - counterion. Our simulations also identify that there is an intervening water layer between the Cl - ion and {N(CH 3 ) 3 } + group: this layer prevents the Cl - ion from coming very close to the {N(CH 3 ) 3 } + group. As a consequence, there is a significantly large mobility of the Cl - ions inside the PMETAC brush layer. Furthermore, the C=O group of the PE chain, due to the partial negative charge on the oxygen atom and the specific structure of the PMETAC brush system, demonstrates strongly hydrophilic behavior and enforces a specific dipole response of water molecules analogous to that experienced by water around anionic species of high charge density. As a result, our findings confirm that PMETAC brushes undergo hydrophilic hydration at one site, apolar hydration at another site, and ensures large mobility of the supported Cl- counterions.

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

Thermodynamic Analysis of Silk Fibroin–Graphite Hybrid Materials and Their Morphology

Silk fibroin (SF) is a β-sheet-rich protein that is responsible for the remarkable tensile strength of silk. In addition to its mechanical properties, SF is biocompatible and biodegradable, making it an attractive candidate for use in biotic/abiotic hybrid materials. A pairing of particular interest is the use of SF with graphene-based nanomaterials (GBNs). The properties of this interface drive the formation of well-ordered nanostructures and can improve the electronic properties of the resulting hybrid. It was previously demonstrated that SF can form lamellar nanostructures in the presence of graphite; however, the equilibrium morphology and associated driving interactions are not fully understood. Here, in this study, we characterize these interactions between SF and SF lamellar with graphite using molecular dynamics (MD) simulations and umbrella sampling (US). We find that SF lamellar nanostructures have strong orientational and spatial preferences on graphite that are driven by the hydrophobic effect, destabilizing solvent–protein interactions and stabilizing protein–protein and protein–graphite interactions. Finally, we show how careful consideration of these underlying interactions can be applied to rationally modify the nanostructure morphology.

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

Structural Properties of [N1888][TFSI] Ionic Liquid: A Small Angle Neutron Scattering and Polarizable Molecular Dynamics Study

In this study, we investigate the quaternary ammonium-based ionic liquid (QAIL), methyltrioctylammonium bis(trifluoromethylsulfonyl)imide, [N 1888 ][TFSI], utilizing small angle neutron scattering (SANS) measurements and polarizable molecular dynamics (MD) simulations to characterize the shortand long-range liquid structure. Scattering structure factors show signatures of three length scales in reciprocal space indicative of alternating polarity (k ~ 0.44 Å –1 ), charge (k ~ 0.75 Å –1 ), and neighboring or adjacent (k ~ 1.46 Å –1 ) domains. Excellent agreement between simulation and experimental scattering structure factors validates various simulation analyses that provide detailed atomistic characterization of the different length scale correlations. The first solvation shell structure is illustrated by obtaining radial, angular, dihedral, and combined distribution functions, where two dominant spatial motifs, N + ···N – and N + ···O – , compete for optimal packing around the polar head of the [N 1888 ] + cation. Intermediate and long-range structures are governed by the balance between local electroneutrality and octyl chain networking, respectively. By computing the charge-correlation structure factor, S ZZ , and the spatial extent of the octyl chain network using graph theory, the bulk-phase structure of [N 1888 ][TFSI] is characterized in terms of electrostatic screening and apolar domain formation length scales.

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

Probing the Solvation Shells of Lithium Ions in Glyme-Based Electrolytes

Glymes have been extensively studied as solvents for Li-battery electrolytes, most recently in equimolar mixtures with lithium bis(trifluoromethanesulfonyl)imide (LiTFSI), due to their ability to form stable solvates. However, directly quantifying free and coordinated glyme molecules in the liquid state has been challenging due to several experimental limitations. Here, in this work, new vibrational probes are demonstrated for studying the solvation structures of diglyme and triglyme in LiTFSI electrolytes. These IR probes make use of an amine group to report the solvation state of glymes at salt-to-solvent molar ratios ranging from 1:5 to 1:10. Characterization of the thermodynamic properties of the solvent exchange occurring in the first solvation shell of lithium ions (Li + ) showed an equilibrium constant for these probes close to unity at room temperature. This result demonstrates that the probes exhibit a similar solvation behavior to their glyme analogue. Concentration dependence studies also revealed a lack of significant amounts of contact ion pairs at the studied concentrations. Moreover, the first solvation shell of Li + appears to be formed by two partially chelating glyme molecules, establishing that even triglyme with multiple chelation sites does not fully coordinate the cation. Complementary molecular dynamics (MD) simulations agree with the experimental results and suggest that at these concentrations, TFSI – predominantly forms solvent-separated ion pairs. However, the simulations do not properly capture the partial solvation structure of the glyme molecules in the solvation shell of Li + as derived from the experiments.

electrolytes↗