Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reactive molecular dynamics”

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 379 records · Page 21

Formation of carbon and oxygen rich surface layer on high purity magnesium by atmospheric carbon dioxide plasma

Carbon and oxygen-rich corrosion barrier layer formed on Mg by a simple and scalable CO 2 atmospheric plasma (CO 2 -AP) process. The reactive CO 2 -AP interacts with the Mg surface and forms a unique layered structure with the top MgCO 3 /MgO-intermixed particulates pillars and the bottom dense layer. The surface features were simultaneously formed on the nano-/micro-structured MgO layer by carbonate molecules, plasma-active CO 2 molecules, and/or other volatile organic compounds on the nano-/micro-structured MgO particle layer. The resulting surfaces after CO 2 -AP were either hydrophobic or hydrophilic and exhibited lower anodic current or high resistance for Mg corrosion. For the hydrophobic surfaces of CO 2 -AP treated Mg, molecular dynamic simulations were performed to understand the origin of hydrophobicity and identified that the amorphous carbon layers formed on the Mg surface are the source. The environmentally benign abundant-gas-based process enables the cost reduction associated with waste treatment, generation of by-product, and supply of raw material.

36 MATERIALS SCIENCE↗

Assessing the Limitations of Self-Interaction-Corrected Functionals for Describing the Hydrated Electron

Simulating the hydrated electron using density functional theory is challenging due to the prevalence of self-interaction error in standard functionals. Hybrid functionals like PBEh(40) can reasonably describe the chemistry of an excess electron in water and partially mitigate self-interaction error by incorporating exact Hartree–Fock exchange, but they are computationally expensive making them impractical for large-scale and long-time ab initio molecular dynamics simulations. Explicit self-interaction correction schemes that are applied on an orbital-by-orbital basis offer a potential alternative when the correction is limited to the singly occupied molecular orbital obtained with a generalized gradient approximation functional. Here, we examine whether the Perdew–Zunger self-interaction correction scheme applied to the revPBE functional can provide a computationally efficient and physically sensible alternative to PBEh(40) for the hydrated electron. We find that functionals incorporating a self-interaction correction scheme should be viewed with caution when applied to the hydrated electron and its reactivity. Furthermore, we show that it is critical to consider extensive sampling and diverse chemical environments when validating their performance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Allosteric control of olefin isomerization kinetics via remote metal binding and its mechanochemical analysis

Abstract Allosteric control of reaction thermodynamics is well understood, but the mechanisms by which changes in local geometries of receptor sites lower activation reaction barriers in electronically uncoupled, remote reaction moieties remain relatively unexplored. Here we report a molecular scaffold in which the rate of thermal E-to-Z isomerization of an alkene increases by a factor of as much as 10 4 in response to fast binding of a metal ion to a remote receptor site. A mechanochemical model of the olefin coupled to a compressive harmonic spring reproduces the observed acceleration quantitatively, adding the studied isomerization to the very few reactions demonstrated to be sensitive to extrinsic compressive force. The work validates experimentally the generalization of mechanochemical kinetics to compressive loads and demonstrates that the formalism of force-coupled reactivity offers a productive framework for the quantitative analysis of the molecular basis of allosteric control of reaction kinetics. Important differences in the effects of compressive vs. tensile force on the kinetic stabilities of molecules are discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Modeling Absolute Redox Potentials of Ferrocene in the Condensed Phase

Absolute thermodynamic quantities for critical chemical reactions are needed to determine the role of solvents and reactive environments in catalysis and electrocatalysis beyond the relative scales typically employed. In principle, theoretical methods can provide such quantification but are often hindered by the innate complexity of strong electron correlation and dynamic relaxation of solvent environments. Here, we present and validate a protocol for calculating the redox potentials of ferrocene/ferrocenium redox pair in the acetonitrile. Equation-of-motion ionization potential coupled-cluster single-double (EOM-IP-CCSD) and effective fragment potential (EFP) methods are used to characterize the adiabatic and vertical ionization potentials (IP) as well as the electron affinity processes. We benchmark molecular mechanics against the EFP model to show the differences in ferrocene electronic polarizability in two redox states. Our best estimate of the redox potential (4.94 eV) agrees well with the experimental value (4.93 eV). This demonstrated the ability of modern computational methods to predict absolute redox potentials quantitatively and, more critically, quantify the correlation of dynamic effects, which underlie their origin.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Potentials of mean force fail to describe chemical bond-breaking in solution

Many liquid phase studies assume that the potential energy surfaces of reacting molecules are the same as in the gas phase, neglecting complex solvent dynamics that can completely alter the nature of chemical reactivity. Even studies that include solvent effects typically only consider them in an average, equilibrium way as part of a potential of mean force (PMF). In this work, we use mixed quantum/classical simulations to compare how equilibrium and non-equilibrium solvent motions affect the photodissociation of a simple diatomic molecule, NaK + , in liquid tetrahydrofuran. A PMF analysis shows that as the excited-state molecule dissociates with the solvent at equilibrium, the bonding electron remains associated with K + at short bond distances but eventually localizes on Na + at the end of dissociation. When we examine non-equilibrium dynamical photodissociation trajectories, however, we find that they fall into three distinct categories: about a quarter of them have the bonding electron mainly associated with Na + , another quarter stay mainly associated with K + , and about half have the bonding electron shared roughly equally between the two ions. The results show that equilibrium PMFs cannot accurately describe the dynamics of bond-breaking chemical reactions in solution because there is insufficient time for the solvent to reach equilibrium on the time scale over which bond dissociation occurs. Furthermore, our analysis shows that the solvent coupling between the electronic energy surfaces is similar at and away from equilibrium, suggesting that other factors, such as solute velocity-driven solvent memory effects, play a more important role in explaining the failure of the equilibrium PMF to predict the non-equilibrium dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Methanol at Water–TiO 2 Interfaces: Free Energies of Water and Methanol Dissociation

Methanol adsorption on TiO 2 surfaces has long been studied due to its role in enhancing photocatalytic hydrogen evolution, yet how it modulates surface chemistry under aqueous conditions remains little understood. Using molecular dynamics with an ab initio-based deep neural network potential, we find that methanol adsorption induces markedly different effects on the aqueous surfaces of anatase and rutile, the two common phases of TiO 2 . In anatase, methanol adsorption significantly enhances water dissociation, which is otherwise rare at the neat water interface. This enhancement arises from an alternative dissociation pathway mediated by surface-bound methoxyl groups. In contrast, methanol adsorption tends to suppress water dissociation on rutile, replacing it with thermodynamically favored methanol dissociation. Overall, methanol adsorption in an aqueous environment alters not only the availability of key reactive intermediates involved in hydrogen evolution but also the hydrogen source, which turns out to be primarily methanol on rutile, whereas both water and methanol are consumed on anatase. These results provide mechanistic insights into the coupled roles of organic adsorbates and water at photocatalytic interfaces, with implications on how methanol enhances the activity of H 2 evolution.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Universal crossed beam imaging studies of polyatomic reaction dynamics

The marriage between high level quantum calculations and experimental advances in laser technology, quantum state control, and detection techniques have opened the door to the study of molecular collision dynamics at a new level of detail. However, one current challenge lies in adapting these powerful strategies to address questions beyond the scope of the small ground state systems that have largely been the focus of reaction dynamics investigations to-date. For molecules with intermediate or large size (more than 6 atoms), lack of spectroscopic information and spectral congestion limit quantum state preparation, control and detection for experiment, and the large number of degrees of freedom of the system makes accurate quantum dynamics calculations prohibitively expensive. Nevertheless, studies of the chemical dynamics of such systems can reveal novel aspects of reactivity not anticipated based upon the behavior of smaller model systems. This Perspective will highlight applications of soft vacuum ultraviolet photoionization at 157 nm as a universal probe in combination with crossed beams and DC slice velocity map ion imaging to study bimolecular reaction dynamics of molecules of intermediate or large size, illuminated with support of high-level ab initio calculations. In this work, we report on the chemical dynamics of atomic oxygen or chlorine reactions with organic compounds: propanol isomers, alkylamines (N(CH 3 ) 3 and NH(CH 3 ) 2 ), and isobutene ((CH 3 ) 2 CCH 2 ) studied using this approach. The polyatomic radical products from the hydrogen abstraction process have been detected by 157 nm photoionization and their slice ion images embody translational energy and angular information that directly reflect the underlying collision dynamics. Various reaction mechanisms (such as direct abstraction and addition–elimination) along with the involvement of roaming dynamics and novel intersystem crossing pathways are presented. These demonstrate the power of this technique to reveal fundamentally new aspects of reaction dynamics that arise in larger reaction systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Oxide– and Silicate–Water Interfaces and Their Roles in Technology and the Environment

Interfacial reactions drive all elemental cycling on Earth and play pivotal roles in human activities such as agriculture, water purification, energy production and storage, environmental contaminant remediation, and nuclear waste repository management. The onset of the 21st century marked the beginning of a more detailed understanding of mineral aqueous interfaces enabled by advances in techniques that use tunable high-flux focused ultrafast laser and X-ray sources to provide near-atomic measurement resolution, as well as by nanofabrication approaches that enable transmission electron microscopy in a liquid cell. This leap into atomic- and nanometer-scale measurements has uncovered scale-dependent phenomena whose reaction thermodynamics, kinetics, and pathways deviate from previous observations made on larger systems. A second key advance is new experimental evidence for what scientists hypothesized but could not test previously, namely, interfacial chemical reactions are frequently driven by “anomalies” or “non-idealities” such as defects, nanoconfinement, and other nontypical chemical structures. Third, progress in computational chemistry has yielded new insights that allow a move beyond simple schematics, leading to a molecular model of these complex interfaces. In combination with surface-sensitive measurements, we have gained knowledge of the interfacial structure and dynamics, including the underlying solid surface and the immediately adjacent water and aqueous ions, enabling a better definition of what constitutes the oxide– and silicate–water interfaces. This critical review discusses how science progresses from understanding ideal solid–water interfaces to more realistic systems, focusing on accomplishments in the last 20 years and identifying challenges and future opportunities for the community to address. Further, we anticipate that the next 20 years will focus on understanding and predicting dynamic transient and reactive structures over greater spatial and temporal ranges as well as systems of greater structural and chemical complexity. Closer collaborations of theoretical and experimental experts across disciplines will continue to be critical to achieving this great aspiration.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

First principles simulation of reacting hypersonic flow over a blunt wedge

This article presents molecular-level analysis of a reactive, near-continuum, Mach 21 nitrogen flow over a blunt wedge using the direct molecular simulation (DMS) method. The flow conditions lead to internal energy excitation and dissociation in the flow field, resulting in thermal and chemical nonequilibrium in the flow. Thermal nonequilibrium in the vibrational mode is observed to extend to the molecular level, where the vibrational energy distributions at various points in the flow field are observed to be non-Boltzmann. Furthermore, this is the first reactive DMS calculation where the wall is assumed to be isothermal and full momentum accommodation of the particles is enforced, hence incorporating viscous wall effects. Since the DMS method uses a quantum mechanically generated interaction potential as its only modeling input, all thermochemical and transport properties of the flow field can directly be attributed to the ab initio potential energy surface. Using the DMS solution as a benchmark, this article assesses the performance of Navier–Stokes computational fluid dynamics solutions using lower fidelity two-temperature models. Two models are chosen as points of comparison: the well-known Park two-temperature model and the recently developed modified Marrone and Treanor model.

Mechanics↗

Coupled Electronic and Nuclear Motions during Azobenzene Photoisomerization Monitored by Ultrafast Electron Diffraction

Ultrafast electron diffraction is a powerful technique that can resolve molecular structures with femtosecond and angstrom resolutions. We demonstrate theoretically how it can be used to monitor conical intersection dynamics in molecules. Specific contributions to the signal are identified which vanish in the absence of vibronic coherence and offer a direct window into conical intersection paths. A special focus is on hybrid scattering from nuclei and electrons, a process that is unique to electron (rather than X-ray) diffraction and monitors the strongly coupled nuclear and electronic motions in the vicinity of conical intersections. Furthermore, an application is made to the cis to trans isomerization of azobenzene, computed with exact quantum dynamics wavepacket propagation in a reactive two-dimensional nuclear space.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Hidden features in the OH-stretching spectra of amino acid decorated air–water interfaces

Chemical reactivity at the air–water interface is governed by the interfacial solvation of reactive species. For instance, during aqueous amino acid-based CO 2 absorption, water reorganizes around the reactive sites and couples dynamically with reaction pathways, facilitating the reaction. In this context, surface-sensitive vibrational sum-frequency generation (vSFG) spectroscopy can probe the OH stretch vibrations of interfacial water and determine the solvation structures around reactants and products, thereby furthering our understanding of the role of interfacial solvation. However, vSFG spectra of the air–water interface in the presence of charged species can be remarkably complex; key species-bound local water structures with distinct orientations may be hidden beneath prominent vSFG peaks arising from water–water hydrogen bonds and remain difficult to resolve. Here, we measure and compute vSFG spectra of the water OH stretch at air–water interfaces decorated with amino acids in their zwitterionic and anionic forms, as well as equimolar mixtures of these forms with bicarbonate. The latter represents post-CO 2 -absorption conditions. We find that computing depth- and frequency-dependent spectral densities—decomposed into contributions from water molecules hydrogen-bonded exclusively to other water molecules, exclusively to amines, exclusively to carboxylates, or shared between these polar/charged groups—is indispensable for accurate interpretation of the vSFG spectra. Key findings include orientational flip-flop in water sub-layers, strong carboxylate-water H-bonding, and water orientational ordering extending into the bulk aqueous phase induced by anionic amino acids. Here, this study provides a computational spectroscopic platform for improved understanding of interfacial solvation relevant to interfacial reactivity.

Air-water interface↗

Mechanistic Insights into Defect-Mediated Crystallization Revealed by Lattice Strain Evolution

Structural defects and lattice strain are intrinsic to many crystalline materials, yet their roles in controlling chemical reaction mechanisms and directing crystallization pathways remain poorly understood. Here, in this study, we revealed the three-dimensional evolution of strain and dislocation defects at the nanoscale during the growth of heterogeneously nucleated barite (BaSO 4 ) and calcite (CaCO 3 ) crystals by using coherent X-ray scattering, electron microscopy, and molecular simulations. Unlike barite, which formed with minimal internal strain, calcite developed dislocation defects and exhibited spatially varying strain that increased during growth. During growth in Sr-rich solutions, calcite likely incorporates Sr 2+ into the defects, which further modulates the local lattice structure and increases both the compressive and tensile strain. These findings suggest that calcite crystallization was likely dominated by attachment of precursor phases, which gave rise to defect-enriched domain structures not predicted by classical growth models. By linking defect formation to ion incorporation and growth dynamics, this work provides fundamental insight into how lattice-level strain heterogeneity governs the chemical reactivity of ionic crystals.

Bragg coherent diffractive imaging↗

New Perspectives on Vibrational Energy Transfer in Energetic Materials: Insights from Pressure-Tuned Ultrafast Spectroscopy

Understanding the manner in which vibrational energy flows between molecular and lattice vibrations is of great interest in physical chemistry due to its central role in reactivity and energy dissipation in molecular materials. Here, in this feature article, we highlight our recent efforts employing ultrafast broadband infrared spectroscopy toward understanding the interplay between molecular and lattice vibrations in energetic materials, motivated by the open questions surrounding the role of vibrational energy transfer (VET) in reaction initiation in these materials. Our work addresses the ongoing debate on the participation of doorway modes in VET. We further present new results from high-pressure ultrafast experiments on RDX, a hydrogen-bonded material, and BNFF, a hydrogen-free material, to explore how intermolecular interaction strength governs VET pathways and time scales. Collectively, our findings reveal that vibrational dynamics in these systems occurs across three distinct time regimes, with VET being incomplete out to hundreds of picoseconds, suggesting the importance of considering nonstatistical reactions in the modeling of these materials. These time scales vary as intermolecular interaction strength is indirectly modified by application of static pressure, indicating dramatic changes to the vibrational structure of these materials under shock-relevant conditions. Our results thus shed light on how intermolecular interactions shape vibrational energy redistribution in molecular materials, and highlight the need for further theoretical and experimental investigation.

Carlson, Daniel Ryan [Sandia National Laboratories↗

Consequences of 131 I Transmutation in Gas Phase Radioiodine Molecules and Adsorbed on Graphite Surface

In the case of nuclear power plant accidents, several radionuclides are being dispersed into the environment. In particular, the release of radioiodine 131 I represents a major health concern because of its absorption by the thyroid gland. To reduce the release of 131 I into the environment, capture materials, such as activated carbon, provide a mean to efficiently retain radioiodine-bound molecular species. Here, we use ab initio molecular dynamic calculations to investigate the effects of 131 I radioactive decay on the molecular stability of I 2 , acid iodine (HI, HIO, HIO 2 , and HIO 3 ), and methyl-iodine species in both gas phase and adsorbed on graphite. In gas phase, we found that I 2 is the least resistant to 131 I transmutation, followed by the HIO x family of molecules, then HI and CH 3 I which are the most resistant. However, the graphite surface is found to greatly favorize the fragmentation of adsorbed molecular species. While the identification of gas phase fragments can offer new chemical route to consider in chemical kinetics models, the formation reactive fragment products at the graphite surface can favorize their reaction with surface carbon atoms which overtime could reduce the number of adsorption sites of capture materials available, and hence, impacting their capture capability.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Analysis of leading edge and trailing edge cover glass samples before and after treatment with advanced satellite contamination removal techniques

Two samples from Long Duration Exposure Facility (LDEF) experiment M0003-4 were analyzed for molecular and particulate contamination prior to and following treatment with advanced satellite contamination removal techniques (CO2 gas/solid jet spray and oxygen ion beam). The pre- and post-cleaning measurements and analyses are presented. The jet spray removed particulates in seconds. The low energy reactive oxygen ion beam removed 5,000 A of photo polymerized organic hydrocarbon contamination in less than 1 hour. Spectroscopic analytical techniques were applied to the analysis of cleaning efficiency including: Fourier transform infrared, Auger, x ray photoemissions, energy dispersive x ray, and ultraviolet/visible. The results of this work suggest that the contamination studied here was due to spacecraft self-contamination enhanced by atomic oxygen plasma dynamics and solar UV radiation. These results also suggest the efficacy for the jet spray and ion beam contamination control technologies for spacecraft optical surfaces.

Hotaling, S. P.↗

Modeling Functional Organic Chemistry in Arctic Rivers: An Idealized Siberian System

Rivers of the Arctic will become ever more important for the global climate, since they carry a majority of continental dissolved organic carbon flux into the rapidly changing polar ocean. Aqueous organics comprise a wide array of functional groups, several of which are likely to impact coastal and open water biophysical properties. Light attenuation, interfacial films, aerosol formation, gas release and momentum exchange can all be cited. We performed Lagrangian kinetic modeling for the evolution of riverine organic chemistry as the molecules in question make their way from the highlands to Arctic outlets. Classes as diverse as the proteins, sugars, lipids, re-condensates, humics, bio-tracers and small volatiles are all included. Our reduced framework constitutes an idealized northward flow driving a major hydrological discharge rate and primarily representing the Russian Lena. Mountainous, high solute and tundra sources are all simulated, and they meet up at several points between soil and delta process reactors. Turnover rates are parameterized beginning with extrapolated coastal values imposed along a limited tributary network, with connections between different terrestrial sub-ecologies. Temporal variation of our total dissolved matter most closely resembles the observations when we focus on the restricted removal and low initial carbon loads, suggesting relatively slow transformation along the water course. Thus, channel combinations and mixing must play a dominant role. Nevertheless, microbial and photochemical losses help determine the final concentrations for most species. Chemical evolution is distinct for the various functionalities, with special contributions from pre- and post-reactivity in soil and delta waters. Several functions are combined linearly to represent the collective chromophoric dissolved matter, characterized here by its absorption. Tributaries carry the signature of lignin phenols to segregate tundra versus taiga sources, and special attention is paid to the early then marine behaviors of low molecular weight volatiles. Heteropolycondensates comprise the largest percentage of reactive carbon in our simulations due to recombination/accumulation, and they tend to be preeminent at the mouth. Outlet concentrations of individual structures such as amino acids and absorbers lie above threshold values for biophysical influence, on the monolayer and light attenuation. The extent of coastal spreading is examined through targeted regional box modeling, relying on salinity and color for calibration. In some cases, plumes reach the scale of peripheral arctic seas, and amplification is expected during upcoming decades. Conclusions are mapped from the Lena to other boreal discharges, and future research questions are outlined regarding the bonding type versus mass release as permafrost degrades. Dynamic aqueous organic coupling is recommended for polar system models, from headwaters to coastal diluent.

54 ENVIRONMENTAL SCIENCES↗

Gaseous Surrogate Hydrocarbons for a Hifire Scramjet that Mimic Opposed Jet Extinction Limits for Cracked JP Fuels

This paper describes, first, the top-down methodology used to define simple gaseous surrogate hydrocarbon (HC) fuel mixtures for a hypersonic scramjet combustion subtask of the HiFIRE program. It then presents new and updated Opposed Jet Burner (OJB) extinction-limit Flame Strength (FS) data obtained from laminar non-premixed HC vs. air counterflow diffusion flames at 1-atm, which follow from earlier investigations. FS represents a strain-induced extinction limit based on cross-section-average air jet velocity, U(sub air), that sustains combustion of a counter jet of gaseous fuel just before extinction. FS uniquely characterizes a kinetically limited fuel combustion rate. More generally, Applied Stress Rates (ASRs) at extinction (U(sub air) normalized by nozzle or tube diameter, D(sub n or t) can directly be compared with extinction limits determined numerically using either a 1-D or (preferably) a 2-D Navier Stokes simulation with detailed transport and finite rate chemistry. The FS results help to characterize and define three candidate surrogate HC fuel mixtures that exhibit a common FS 70% greater than for vaporized JP-7 fuel. These include a binary fuel mixture of 64% ethylene + 36% methane, which is our primary recommendation. It is intended to mimic the critical flameholding limit of a thermally- or catalytically-cracked JP-7 like fuel in HiFIRE scramjet combustion tests. Our supporting experimental results include: (1) An idealized kinetically-limited ASR reactivity scale, which represents maximum strength non-premixed flames for several gaseous and vaporized liquid HCs; (2) FS characterizations of Colket and Spadaccini s suggested ternary surrogate, of 60% ethylene + 30% methane + 10% n-heptane, which matches the ignition delay of a typical cracked JP fuel; (3) Data showing how our recommended binary surrogate, of 64% ethylene + 36% methane, has an identical FS; (4) Data that characterize an alternate surrogate of 44% ethylene + 56% ethane with identical FS and nearly equal molecular weights; this could be useful when systematically varying the fuel composition. However, the mixture liquefies at much lower pressure, which limits on-board storage of gaseous fuel; (5) Dynamic Flame Weakening results that show how oscillations in OJB input flow (and composition) can weaken (extinguish) surrogate flames up to 200 Hz, but the weakening is 2.5x smaller compared to pure methane; and finally, (6) FS limits at 1-atm that compare with three published 1-D numerical OJB extinction results using four chemical kinetic models. The methane kinetics generally agree closely at 1-atm, whereas, the various ethylene models predict extinction limits that average ~ 45% high, which represents a significant problem for numerical simulation of surrogate-based flameholding in a scramjet cavity. Finally, we continue advocating the FS approach as more direct and fundamental for assessing idealized scramjet flameholding potentials than measurements of "unstrained" premixed laminar burning velocity or blowout in a Perfectly Stirred Reactor.

Pellett, Gerald L.↗

Understanding Solid-Electrolyte Interphase Formation at the Lithium Metal Anode of Lithium-Sulfur Batteries

Lithium-sulfur (Li-S) batteries are one of the most promising chemistries for the next generation of energy storage devices. In a standard Li-S cell, the metallic lithium anode is an essential component of the battery due to its low density, extremely high theoretical specific capacity, and very low negative electrochemical potential. However, several challenges related to the use of Li metal have prevented this battery technology from becoming commercially available. For instance, Li-anodes are highly reactive, which results in the continuous decomposition of the electrolyte and the formation of the solid-electrolyte interphase (SEI) layer. This, combined with the non-uniform deposition of Liions during plating at the anode surface, can seriously affect the performance, cycling, and safety of the battery. In addition, it has been suggested that a controlled SEI formation at the metallic Li anode can yield enhanced battery performance. Hence, a comprehensive understanding of how the SEI layer is formed can help elucidate improvements in this battery technology. In this work, the formation of the SEI layer in Li-S batteries is investigated using density functional theory (DFT) and ab initio molecular dynamics (AIMD) simulations. The stability of liquid organic-based electrolyte components, salt concentration, electron-rich environments, and the use of inorganic solid-state electrolyte (SSE) materials are explored in order to provide molecular-level fundamental insights into how the nature and composition of the electrolyte can alter the initial stages of the SEI formation. The products of electrolyte decomposition and electro-/chemical reaction pathways of liquid electrolyte component such as solvent and salts are predicted from AIMD simulations of electrolyte-anode interfaces. Energetics of reactions from DFT ground-state optimizations are also presented to confirm the decomposition mechanisms. Some significant differences are then drawn regarding the use of low and high salt concentrations. Simulations under electron-rich environments show additional multielectron electrochemical reactions of solvent and salt decomposition taking place due to the excess of electrons and the presence of radical anions in the solution. Finally, the initial stages of formation, stability, and main constituents of interfaces between Li-metal and Sbased electrodes with sulfide-based SSEs are characterized in detail.

25 ENERGY STORAGE↗