Search NASA⌕ Search

SEARCH · Search NASA

Results for “Density functional theory simulations”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Accelerated Steam Methane Reforming by Dynamically Applied Charges

Catalyst design has traditionally focused on tuning active site properties to optimally bind reaction intermediates and balance the kinetic requirements of multiple competing chemical processes, as necessitated by the Sabatier principle. It has recently been proposed that for reactions following certain potential energy landscapes, the activity limit imposed by the Sabatier principle may be overcome by using programmed oscillations of surface electron density at the timescales of surface reactions (i.e., “catalytic resonance”). Here, we use a combination of density functional theory (DFT) simulations and transient kinetic models (TKMs) to simulate the kinetics of steam methane reforming (SMR) on Ru(211) surfaces under statically and dynamically applied charges. DFT-calculated binding energies of SMR intermediates and transition states exhibit strong sensitivity to positively applied charges and follow unique scaling relationships that deviate from linear periodic trends across transition metals. Our simulations demonstrate that applying a small positive charge to Ru dramatically enhances the steady-state turnover frequency (TOF) of SMR by up to 5 orders of magnitude above the TOF observed over neutral Ru. Thus, statically charging Ru catalysts may be an effective strategy to lower the temperature requirements for SMR. Dynamic square-wave oscillations in charge resulted in SMR catalytic resonance with an onset frequency f ∼ 106 Hz and the corresponding average TOFs exceeding the statically charged Ru surface by an additional 15%. Here, based on sensitivity analyses performed for the two end points of oscillation, we propose that dynamic TOF improvement beyond the Sabatier maximum can be expected when the system oscillates between two kinetic regimes that are uniquely controlled by distinct elementary steps.

Catalysts↗
DOE OSTI2025

Understanding Trigger Linkage Dynamics in Energetic Materials Using Mixed Picramide Nitrate Ester Explosives

The ability to predict the handling sensitivity of new organic energetic materials has been a longstanding goal. We report the synthesis and characterization of six new nitropicramide energetic materials with mixed functional groups that mimic known explosives such as nitroglycerin, erythritol tetranitrate (ETN), and pentaerythritol tetranitrate (PETN). The molecules have been studied theoretically using quantum molecular dynamics (QMD) simulations and density functional theory (DFT) calculations to identify the weakest bond in the reactants - the trigger-linkages - which control handling sensitivity, and to quantify their specific enthalpies of explosion. In good accord with the drop weight impact sensitivity data, our calculations predict that the sensitivities of the molecules are very similar owing to the small variations of the energy output and rates of trigger linkage rupture. In addition, both the QMD and DFT calculations point to the nitropicramide N–NO 2 bonds as the trigger linkages rather than the more typical O–NO 2 bonds. We propose that the switch of the trigger linkage from the nitrate esters to the nitramine groups arises from the strongly electron withdrawing character of the adjacent trinitrobenzene groups.

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

Theoretical Calculation Guided Design of Single-Atom Catalysts toward Fast Kinetic and Long-Life Li–S Batteries

Lithium–sulfur (Li–S) batteries are enticing next-generation energy storage technologies due to their high theoretical energy density, environmental friendliness, and low cost. Yet, low conductivity of sulfur species, dissolution of polysulfides, poor conversion from sulfur reduction, and lithium sulfide (Li 2 S) oxidation reactions during discharge–charge processes hinder their practical applications. Herein, under the guidance of density functional theory calculations, we have successfully synthesized large-scale single atom vanadium catalysts seeded on graphene to achieve high sulfur content (80 wt % sulfur), fast kinetic (a capacity of 645 mAh g –1 at 3 C rate), and long-life Li–S batteries. Both forward (sulfur reduction) and reverse reactions (Li 2 S oxidation) are significantly improved by the single atom catalysts. This finding is confirmed by experimental results and consistent with theoretical calculations. The ability of single metal atoms to effectively trap the dissolved lithium polysulfides (LiPSs) and catalytically convert the LiPSs/Li 2 S during cycling significantly improved sulfur utilization, rate capability, and cycling life. Our work demonstrates an efficient design pathway for single atom catalysts and offers solutions for the development of high energy/power density Li–S batteries.

25 ENERGY STORAGE↗
DOE OSTI2020

Magnetic Tunability in RE-DOBDC MOFs via NO x Acid Gas Adsorption

The magnetic susceptibility of NO x -loaded RE-DOBDC (rare earth (RE): Y, Eu, Tb, Yb; DOBDC: 2,5-dihydroxyterephthalic acid) metal–organic frameworks (MOFs) is unique to the MOF metal center. RE-DOBDC samples were synthesized, activated, and subsequently exposed to humid NO x . Each NO x -loaded MOF was characterized by powder X-ray diffraction, and the magnetic characteristics were probed by using a VersaLab vibrating sample magnetometer (VSM). Lanthanide-containing RE-DOBDC (Eu, Tb, Yb) are paramagnetic with a reduction in paramagnetism upon adsorption of NO x . Y-DOBDC has a diamagnetic moment with a slight reduction upon adsorption of NO x . The magnetic susceptibility of the MOF is determined by the magnetism imparted by the framework metal center. The electronic population of orbitals contributes to determining the extent of magnetism and change with NO x (electron acceptor) adsorption. Eu-DOBDC results in the largest mass magnetization change upon adsorption of NO x due to more available unpaired f electrons. Experimental changes in magnetic moment were supported by density functional theory (DFT) simulations of NO x adsorbed in lanthanide Eu-DOBDC and transition metal Y-DOBDC MOFs.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Task Sharing of Proton Incorporation in Vertically Aligned Nanocomposite Triple Conductors: Growth, Structure, and Surface Exchange Kinetics

As protonic ceramic electrolysis cells emerge for efficient H 2 production, there is a need to develop air electrode materials enabling fast, durable steam splitting and proton incorporation. Single-phase triple conductors may fail to satisfy the myriad performance/stability requirements, and their critical charge-carriers (holes, oxygen vacancies, and protons) are in competition, limiting their concentrations. Instead, we propose task-sharing, vertically aligned nanocomposites (VANs), comprising a proton conductor (BaZr 0.9 Y 0.1 O 3-δ ) and a redox-active mixed ionic electronic conductor (Ce 0.9 Pr 0.1 O 2-δ ), that may enable rapid proton surface exchange at the solid–gas interface and transport along the solid–solid heterointerfaces. We grew VANs by pulsed laser deposition and investigated the interplay between their processing conditions, structure, and proton and oxygen surface exchange kinetics. We varied the substrate temperature, laser repetition rate, laser fluence, and processing oxygen pressure. The crystallinity and phases were characterized by grazing-incidence X-ray diffraction, and the strain and structural order as a function of depth were evaluated by angle-dependent synchrotron X-ray pair distribution function analysis. To evaluate the potential for interdiffusion, the formation energies of substitutional defects were simulated with density functional theory. Corresponding structural analysis and elemental mapping were performed by scanning/transmission electron microscopy, energy-dispersive X-ray spectroscopy, and electron energy-loss spectroscopy, indicating distinct nanoscale compositional regions with a hierarchical structure embedded in individual VANs columns and minimal interdiffusion across a bilayer film. Proton and oxygen surface exchange coefficients (k H , k O ) and polarization resistances were evaluated by electrical and optical relaxations and impedance spectroscopy of VAN-incorporated protonic ceramic electrochemical cells, respectively, at 400–500 °C, demonstrating values comparable to some of the best-known triple and mixed conductors.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Dioxygen Activation Kinetics over Distinct Cu Site Types in Cu-Chabazite Zeolites

Cu-exchanged zeolites activate dioxygen to form active sites for partial methane oxidation (PMO), nitrogen oxide decomposition, and carbon monoxide oxidation. Apparent rates of O 2 activation depend both on the intrinsic kinetics of distinct Cu site types and the distributions of such sites within a given zeolite, which depend on the density and arrangement of the framework Al atoms. In this work, we use hydrothermal synthesis methods to control the arrangement of framework Al sites in chabazite (CHA) zeolites and, in turn, the distinct Cu site types formed. Time-resolved in situ Raman spectroscopy reveals the kinetics of O 2 adsorption and activation within these well-defined Cu-CHA materials and the concomitant structural evolution of copper–oxygen (Cu x O y ) complexes, which are interpreted alongside Cu(I) oxidation kinetics extracted from in situ X-ray absorption spectroscopy (XAS). Raman spectra of several plausible Cu x O y species simulated using density functional theory suggest that experimental spectra (λ ex = 532 nm) capture the formation of mono(μ-oxo)dicopper species (ZCuOCuZ). Transient experiments show that the timescales required to form Cu x O y structures that no longer change in Raman spectra correspond to the durations of oxidative treatments that maximize CH 3 OH yields in stoichiometric PMO cycles (approximately 2 h). Yet, these periods extend well beyond the timescales for the complete conversion of the initial Cu(I) intermediates to their Cu(II) states (<0.3 h, reflected in X-ray absorption near edge spectroscopy spectra), which demonstrates that Cu x O y complexes continue to evolve structurally following rapid Cu(I) oxidation. The dependence of ZCuOCuZ formation rates on O 2 pressure, H 2 O pressure, and temperature is consistent with a mechanism in which ZCuOH reduces to form ZCu + sites that bind molecular oxygen and form ZCu-O 2 intermediates. Subsequent reaction with proximate ZCu + forms bridging peroxo dicopper complexes that cleave O–O bonds to form ZCuOCuZ in steps facilitated by water. These data and interpretations provide evidence for the chemical processes that link rapid and kinetically irrelevant Cu oxidation steps (frequently probed by XAS and UV–vis spectroscopy) to the relatively slow genesis of reactive Cu complexes that form CH 3 OH during PMO. In doing so, we reveal previously unrecognized complexities in the processes by which Cu ions in zeolites activate O 2 to form active Cu x O y complexes, which underscore the insight afforded by judicious combinations of experimental and theoretical techniques.

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

Modulating Operational Conditions to Mitigate Deactivation in Formate Dehydrogenation on Pd Phases

In heterogeneous catalysis, poisoning by surface-bound intermediates poses a major barrier to sustained catalyst performance in (de)hydrogenation reactions. Formate/bicarbonate systems, as liquid organic hydrogen carriers (LOHCs), offer a CO 2 -integrated, low-temperature pathway for hydrogen storage and release, making them attractive for circular energy applications. However, their lower hydrogen density and susceptibility to catalyst deactivation limit their competitiveness compared to conventional LOHCs like methylcyclohexane. Here, this study investigates the mechanistic origins of formate (HCOO – ) dehydrogenation and associated deactivation on Pd interfaces. Using density functional theory (DFT) simulations, we show that under thermocatalytic conditions, strongly bound formate accumulates on the catalyst surface (Pd(111)), blocking active sites, raising activation barriers, and leading to progressive performance loss. Because formate adsorption involves charge transfer, we exploit its sensitivity to electronic structure by modulating the electrochemical potential of the catalyst. Our results reveal that hydrogen transfer from water and formate exhibits opposing potential dependencies, providing insights into the opposing driving forces behind both catalytic activity and poisoning. To further probe this phenomenon, we examine electrochemically induced phase transitions in Pd, focusing on PdO(100) and PdH(110), which are stable under oxidizing and reducing potentials, respectively, and demonstrate enhanced dehydrogenation activity between −0.4 and 0.2 V vs standard hydrogen electrode (SHE). Complementary thermal treatments help decouple kinetic and thermodynamic contributions to intermediate binding. These findings underscore the critical role of the catalyst phase and external stimuli in dictating poison-active site interactions and highlight phase engineering as a promising strategy to mitigate deactivation. This work offers mechanistic insights and design principles for developing more resilient and efficient catalysts for LOHC applications under realistic operating conditions.

Pd phase↗
DOE OSTI2024

Optimizing Nonaqueous Sodium–Polysulfide Redox-Flow Batteries: The Role of Solvation Effects with Glyme Solvents

Nonaqueous redox-flow batteries (NARFBs) that use economical alkali metals and the corresponding metal polysulfides are highly attractive for grid-scale energy storage. Although sodium–sulfur systems have been recognized as promising candidates and have been the focus of many studies due to their high earth abundance and energy density, an understanding of the role of the solvation chemistry of commonly used glyme solvents is missing. Herein, we report a systematic investigation into the solvation effects of glyme-based Na-S electrolytes through comprehensive physiochemical experiments and Density Functional Theory (DFT) simulations. Our findings revealed, on one hand, that an optimal coordination strength between glymes and Na + could maintain a relatively smooth Na + diffusion. On the other hand, glyme solvents with extended chain lengths shift the reduction potential of S 8 2– negatively to elevate the formation barrier of undesirable short-chain polysulfides (S n 2– , n ≤ 4) that have high membrane permeability. This solvation phenomenon not only mitigates capacity fading but also extends the operational longevity of the Na-S NARFBs. In conclusion, the results underscore the critical roles of balanced solvent–cation interactions and controlled redox potentials in improving the stability and efficiency of Na-S NARFB systems, marking a significant advancement in the development of sustainable energy storage solutions.

25 ENERGY STORAGE↗
DOE OSTI2024

Epitaxial Metal Electrodeposition Controlled by Graphene Layer Thickness

Control over material structure and morphology during electrodeposition is necessary for material synthesis and energy applications. One approach to guide crystallite formation is to take advantage of epitaxy on a current collector to facilitate crystallographic control. Single-layer graphene on metal foils can promote “remote epitaxy” during Cu and Zn electrodeposition, resulting in growth of metal that is crystallographically aligned to the substrate beneath graphene. However, the substrate–graphene–deposit interactions that allow for epitaxial electrodeposition are not well understood. Here, we investigate how different graphene layer thicknesses (monolayer, bilayer, trilayer, and graphite) influence the electrodeposition of Zn and Cu. Scanning transmission electron microscopy and electron backscatter diffraction are leveraged to understand metal morphology and structure, demonstrating that remote epitaxy occurs on mono- and bilayer graphene but not trilayer or thicker. Density functional theory (DFT) simulations reveal the spatial electronic interactions through thin graphene that promote remote epitaxy. This work advances our understanding of electrochemical remote epitaxy and provides strategies for improving control over electrodeposition.

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

Single-Phase Spinel NiCo 2 O 4 as Highly Active and Stable Electrocatalysts for Urea Oxidation Reaction in Urea Electrolysis

Exploring and designing a stable and active catalyst for the urea electro-oxidation reaction (UOR, CO(NH 2 ) 2 + 6OH – → CO 2 + N 2 + 5H 2 O + 6e – ) is crucial for the long-term sustainability of ecological systems and clean energy production. We found that spinel NiCo 2 O 4 is a stable and active electrocatalyst for UOR at a relatively low anodic potential without triggering the competing oxygen evolution reaction (OER). A urea electrolysis cell (CO(NH 2 )2 + H 2 O → CO 2 + N 2 + 2H 2 ) utilizing a spinel NiCo 2 O 4 anode and a commercial Pt cathode was further characterized through galvanostatic polarization tests, demonstrating excellent structural stability at various current densities. Post-mortem analysis of long-term urea electrolysis measurements suggested that NiCo 2 O 4 electrocatalysts maintained a stable spinel structure. However, redistribution of Ni 3+ to Ni 2+ valence on the catalyst surface was observed, in contrast to the intact Co valence, indicating that (i) Ni sites are active toward urea adsorption and sequential electro-oxidation; (ii) while urea oxidation proceeds primarily through the direct electro-oxidation mechanism, chemical reactions between the Ni 3+ site and urea occur during long-term electrochemical UOR operation. Density functional theory (DFT) simulations were used to calculate the adsorption energies of urea molecules on NiO, Co 3 O 4 , and NiCo 2 O 4 , revealing the importance of regulating the configuration of adsorbed urea molecules on the NiCo 2 O 4 surface.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Impact of Hydrogen Bonds on CO 2 Binding in Eutectic Solvents: An Experimental and Computational Study toward Sorbent Design for CO 2 Capture

Choline-based amino acid ionic liquids with anions glycinate, β-alaninate, phenylalaninate, and prolinate were synthesized and mixed with ethylene glycol to form lower-viscosity benign eutectic solvents for CO 2 capture. The highest capacity measured was 0.7 moles of CO 2 per mole of ionic liquid (2 moles CO 2 per kg solvent) for a 1 to 2 mole ratio mixture of choline prolinate to ethylene glycol at 1 bar of CO 2 and 25 °C. Under 5000 ppm of CO 2 , half of this capacity was realized. Here, through a combined study of quantitative 13 C NMR spectroscopy, molecular dynamics simulations and density functional theory calculations, we show that hydrogen bonding in the eutectic solvent prevents proton-transfer between prolinate anions upon CO 2 absorption, which occurs in the absence of ethylene glycol and deactivates binding sites. Blocking this proton transfer leads to a higher binding capacity compared to neat choline prolinate. This work demonstrates the impact of hydrogen bonding on the CO 2 binding mechanism and energetics, as well as physical and thermal properties in eutectic solvents, thus addressing an unmet need and informing future studies on the development of benign sorbents for capturing CO 2 from dilute streams.

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

Pathways for Electron Transfer at MgO–Water Interfaces from Ab Initio Molecular Dynamics

The nature of electron transfer across metal oxide–water interfaces depends significantly on the band gap of the oxide and its band edge energies relative to the potentials of relevant aqueous redox couples. Here we focus on the water interface with MgO, a prototypical wide band gap oxide whose conduction band edge is close in energy to that of water. We investigate the behavior of an excess electron at and out of equilibrium near the interface using ab initio molecular dynamics based on hybrid density functional theory. Our simulations show that under equilibrium conditions the excess electron (donated by an Al impurity in MgO) localizes to a midgap defect state comparable in energy and shape to a hydrated electron in bulk water. To characterize the electron transfer from the conduction band of MgO to interfacial product states, we dope near-equilibrium configurations of the pristine MgO–water system with Al and run short trajectories of these instantaneously out-of-equilibrium systems. Here we observe two distinct products associated with the excess electron: a surface-localized electron (e surf –) and an aqueous hydrogen radical (H • ). The H • pathway exhibits a much higher activation barrier despite being more exoergic, making e surf – the kinetic product. Our characterization of the pathways on the basis of Marcus theory is consistent with the poor observed utility of MgO for water radiolysis. Moreover, we anticipate that the computational framework employed here will be broadly applicable to assessing electron transfer mechanisms at aqueous, photocatalytic interfaces.

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

Atomic-Scale Imaging Reveals Polar-π Interactions in Two-Dimensional Molecular Superlattices

Controlling coassembly of synthetic oligomers into binary superlattices at the atomic level is challenging. Here, we report a strategy for programming polar-π interactions in oligomeric peptoids, a class of sequence-defined peptidomimetics, facilitating the formation of homogeneous two-dimensional (2D) superlattices. N-2-phenylethyl and N-(2-perfluorophenyl)ethyl side chains, similar in size, but with contrasting electrostatic characteristics, were introduced at defined sequence positions to generate favorable dipolar aromatic interactions. The resulting nanosheets exhibit different crystal motifs depending on the side chain interactions: systems containing only one type of aromatic side chain form a parallel V-shaped motif driven by π-π interactions, whereas a combination of both types of aromatic side chains, either within one backbone or through the coassembly of two distinct peptoids, adopt an antiparallel V-shaped superlattice with higher thermal stability, driven by polar-π interactions. Cryogenic transmission electron microscopy directly resolved the packing arrangement of perfluorophenyl and phenyl rings in individual nanosheet superlattices, confirming that intermolecular polar-π interaction dominates the superlattice motifs and increases lattice stability. Molecular dynamics simulations and density functional theory calculations further substantiate the energetic favorability of polar-π interactions over π-π interactions, rationalizing the formation of homogeneous superlattices with enhanced thermal stability. Our discoveries establish a design principle for binary coassembly using sequence-defined oligomers, which enables control over unit cell geometry, lattice stability, and molecular registration through aromatic side chain polarization and sequence control. This ability to program atomic-scale binary superlattices opens new avenues for designing functional 2D soft materials.

Lee, Yen Jea [Lawrence Berkeley National Laborator↗
DOE OSTI2024

Evaluating the Role of Metastable Surfaces in Mechanochemical Reduction of Molybdenum Oxide

Mechanochemistry and mechanocatalysis are gaining increasing attention as environmentally friendly chemical processes because of their solvent-free nature and scalability. Significant effort has been devoted for studying continuum-scale phenomena in mechanochemistry, such as temperature and pressure gradients, but the atomic-scale mechanisms remain relatively unexplored. In this work, we focus on the mechanochemical reduction of MoO3 as a case study. We use experimental techniques to determine the mechanochemical reduction conditions and density functional theory (DFT) simulations to establish an atomistic framework for identifying the metastable surfaces that are most likely to enable this process. Our results show that metastable surfaces can significantly lower or remove thermodynamic barriers for surface reduction and that kinetic energy from milling can facilitate the formation of metastable surfaces that have high surface fracture energies and are not thermally accessible. These findings indicate that metastable surfaces are an important aspect of mechanochemistry along with hot spots and other continuum-scale phenomena.

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

99 TcO 4 - removal from legacy defense nuclear waste by an alkaline-stable 2D cationic metal organic framework

Removal of 99 TcO 4 - from legacy defense nuclear tank waste at Savannah River Site is highly desirable for the purpose of nuclear safety and environmental protection, but currently not achievable given the extreme conditions including high alkalinity, high ionic strength, and strong radiation field. Herein, we present a potential solution to this long-term issue by developing a two-dimensional cationic metal organic framework SCU-103, showing ultrahigh stability in alkaline aqueous media and great resistance to both β and γ radiation. More importantly, it is very effective for 99 TcO 4 - separation from aqueous media as demonstrated by fast exchange kinetics, high sorption capacity, and superior selectivity, leading to the successful removal of 99 TcO 4 - from actual Savannah River Site high level tank waste for the first time, to the best of our knowledge. In addition, the uptake mechanism is comprehensively elucidated by molecular dynamics simulation and density functional theory calculation, showing a unique chemical recognition of anions with low charge density.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2021

Stability and molecular pathways to the formation of spin defects in silicon carbide

Abstract Spin defects in wide-bandgap semiconductors provide a promising platform to create qubits for quantum technologies. Their synthesis, however, presents considerable challenges, and the mechanisms responsible for their generation or annihilation are poorly understood. Here, we elucidate spin defect formation processes in a binary crystal for a key qubit candidate—the divacancy complex (VV) in silicon carbide (SiC). Using atomistic models, enhanced sampling simulations, and density functional theory calculations, we find that VV formation is a thermally activated process that competes with the conversion of silicon (V Si ) to carbon monovacancies (V C ), and that VV reorientation can occur without dissociation. We also find that increasing the concentration of V Si relative to V C favors the formation of divacancies. Moreover, we identify pathways to create spin defects consisting of antisite-double vacancy complexes and determine their electronic properties. The detailed view of the mechanisms that underpin the formation and dynamics of spin defects presented here may facilitate the realization of qubits in an industrially relevant material.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Solvent oligomerization pathways facilitated by electrolyte additives during solid-electrolyte interphase formation

The solid–electrolyte interphase (SEI) layer formation is known to play an important role in determining the lifetime of lithium-ion batteries. A thin, stable SEI layer is linked to overall improved battery performance and longevity, however, the factors and mechanisms that lead to optimal SEI morphology and composition are not well understood. Inclusion of electrolyte additives (fluoroethylene carbonate, FEC; and vinylene carbonate, VC) is often necessary for improving SEI characteristics. To understand how these electrolyte additives impact SEI formation, we employed molecular dynamics (MD) and density functional theory (DFT) simulations to study the reaction networks and oligomerization pathways, respectively, for three systems containing ethylene carbonate (EC), a lithium ion, and FEC or VC. MD simulations suggest radical oligomerization pathways analogous to traditional oligomerization with nucleophilic alkoxide species via SN1 reaction mechanisms. Both S N 1 and S N 2 mechanisms were studied for all three systems using DFT. Oligomerization reactions were studied with both a standard alkoxide species and a ring-opened EC radical as the nucleophiles and EC, FEC, and VC as the electrophiles. For all cases, FEC and VC exhibited lower free energy barriers and more stable adducts when compared with EC. We conclude that one of the role of additives is to modify the oligomerization process of EC by introducing branching points (FEC) or termination points (VC).

25 ENERGY STORAGE↗