Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DENSITY OF STATES”

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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Modulating the Electron Affinity of Small Bipyridyl Molecules on Single Gold Nanoparticles for Plasmon-Driven Electron Transfer

Developing controlled platforms for plasmon-driven chemistry is of great importance in catalytic reactions at the nanoscale. We report anion radical formation for five bipyridyl complexes of varying degrees of electron affinity utilizing optically fo-cused intraband (594 nm) and interband (532 nm) pump excitation of single gold nanoparticles. The surface-enhanced Raman scattering (SERS) of anion radicals for the five non-resonant adsorbed molecules 2,2’-bipyridine (22BPY), 4,4’-bipyridine (44BPY), trans-1,2-bis(4-pyridyl)ethylene (BPE), 1,2-bis(4-pyridyl)acetylene (BPA), and 1,2-bis(4-pyridyl)ethane (BPEt) were detected using localized surface-plasmon resonance (LSPR) excitation with 785 nm. The electron affinity of the five bipyridyl complexes were determined using electrochemistry. Molecules with low electron affinity experienced high-er instances of radical anion formation under a plasmon-coupled intraband electron transfer excitation (594 nm) whereas molecules with high electron affinity showed a preference for anion radical formation under direct interband electron trans-fer excitation (532 nm). The lowest unoccupied molecular orbital (LUMO) energy levels for low electron affinity surface-bound molecules (22BPY, BPEt) are on average ~0.43 eV higher for high electron affinity surface-bound molecules (BPA, BPE, 44BPY) as calculated using time-dependent density functional theory, elucidating the importance of plasmon coupling to energy levels that facilitate charge transfer pathways. We also show the ability to ‘activate’ high vs low electron affinity single nanoparticles with the choice of pump excitation wavelength. The findings show the complex interplay between molecular electron affinity, orbital overlap with the density of states of the plasmonic metal, and excitation energetics of the pump laser wavelength. Potential applications of this work include enhanced control over molecular scale catalysis, biosensor design, and solar energy capture.

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

Ab-Initio Investigation of Finite Size Effects in Rutile Titania Nanoparticles with Semilocal and Nonlocal Density Functionals

In this work, we employ hybrid and generalized gradient approximation (GGA) level density functional theory (DFT) calculations to investigate the convergence of surface properties and electronic gap of rutile titania nanoparticles with particle size. The surface energies and electronic gaps are calculated for cuboidal particles with minimum dimension ranging from 3.7 Angstrom (24 atoms) to 10.3 Angstrom (384 atoms) using a highly-parallel real-space DFT code to enable hybrid level DFT calculations of larger nanoparticles than are typically practical. We deconvolute the geometric and electronic finite size effects in surface energy, and evaluate the influence of defects on electronic gap and density of states (DOS). The electronic finite size effects in surface energy vanish when the minimum length scale of the nanoparticles becomes greater than 10 Angstrom. We show that this length scale is consistent with a computationally efficient numerical analysis of the characteristic length scale of electronic interactions. The surface energy of nanoparticles having minimum dimension beyond this characteristic length can be approximated using slab calculations that account for the geometric defects. In contrast, the finite size effects on the electronic gap and DOS is highly dependent on the shape and size of these particles. Furthermore, the DOS for cuboidal particles and more realistic particles constructed using the Wulff algorithm reveal that defect states within the electronic gap play a key role in determining the eigen value distribution of nanoparticles and the electronic gap does not converge to the bulk limit for the particle sizes investigated.

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

Four-Terminal Electrochemistry: A Back-Gate Controls the Electrochemical Potential of a 2D Working Electrode

We demonstrate that ultrathin semiconductor working electrodes integrated into metal–insulator–semiconductor (MIS) stacks are an enabling platform for understanding non-Faradaic semiconductor electrochemistry. Furthermore, 5 nm thick ZnO electrodes were deposited on 30 nm HfO 2 dielectric on a Pd “gate” electrode. Application of a bias V G between the Pd gate and the ZnO electrode causes electrons to accumulate in the ZnO layer as measured by recording the in-plane sheet conductance. By contacting the top surface of the ZnO layer with the electrolyte in a conventional three-electrode electrochemical cell, we show that the gate voltage V G modulates the electrochemical potential V ZnO of the ZnO film with respect to a reference electrode. Electrochemical potential changes ΔV ZnO up to –1 V vs Ag/Ag + are achieved for V G = +7 V. Furthermore, by measuring V ZnO vs V G , we extract the quantum capacitance CQ of the ZnO film as a function of the Fermi-level position, which provides a direct measure of the ZnO electronic density of states (DOS). Finally, we demonstrate that the gated ZnO working electrodes can disentangle the two principal components of electrochemical potential, namely, the Fermi-level shift Δδ and the double-layer charging energy eΔΦ EDL . This disentanglement hinges on a fundamental difference between back-gating and normal electrochemical control, namely, that electrochemical control requires double-layer charging, while back-gate control does not. Collectively, the results show that the backside gate electrode is an effective fourth terminal that enables measurements that are difficult to achieve in conventional three-terminal electrochemical setups.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Strain-Induced Shifts in Defective Graphite Phonon Modes Predicted by Density Functional Theory

Carbon fiber composites have gained attention as a structural material because of their high strength-to-weight ratio, and understanding the effect of defects on reactivity and mechanical properties is important for the longevity and safety of the composite. Although it is known that strain causes the underlying graphitic vibrational modes to redshift, it is not clear how strain may alter reactivity and defect-induced vibrational changes. To investigate the strain-induced phonon changes of defective carbon fiber composites, density functional theory calculations of graphite are used, including intercalated hydrogen and fluorine defects. By comparing changes in the bond lengths, formation energies, and phonon density of states for uniaxially and biaxially strained graphite, strain was found to generally make defect formation more favorable and the specific behavior changes are dependent on the strain direction and defect identity. Specifically, intercalated fluorine phonons are more sensitive to strain than hydrogen intercalation phonons, and strain applied along the zigzag direction alters the calculated properties more than strain along the armchair direction. Finally, these results highlight the importance of understanding the microstructural effect of deviations from the ideal material because small changes in strain or defect type can significantly alter the behavior of the carbon fiber composite core.

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

Fundamental Insights into Cathode Stability: Linking Compositional Tuning and Local Coordination in Complex Metal Oxides under Aqueous Transformations

Compositional tuning of complex metal oxides in Li-ion battery materials influences their performance as well as their end-of-life behavior, in particular, the tendency to release toxic metal cations in aqueous solution. We modeled ternary variants of a parent LiCoO 2 delafossite structure by varying the metal identity and relative amounts. This yielded ten model formulations of Li(A 4/6 B 1/6 C 1/6 )O 2 , where the material is enriched with the A metal and doped with B and C, with Ni, Mn, Co, Fe, Al, V, and Ti as constituent metals. To assess their stability in aqueous conditions, metal release energetics were calculated using a combination of Density Functional Theory calculations and thermodynamics. Metal release in ternary oxides is dictated by subtle variations in the coordination environment of the leaving group. To identify governing chemical features across diverse compositions with varying local coordination environments, we leverage random forest regression and descriptor importance analysis. A key result is that metal–oxygen orbital hybridization, quantified using a projected density-of-states-derived descriptor, H d/p , provides a physically grounded measure of interaction strength that governs metal release energetics. This refined perspective goes beyond conventional oxidation state considerations and offers more robust insights for materials science. Finally, we model defect surface-bound O 2 dimer formation as a proxy for reactive oxygen species (ROS) generation. The results show that Ni-rich compositions more readily stabilize spin-polarized O 2 dimers, corroborating experimental reports of an increased ROS-driven biological response. In conclusion, our results establish a compositional and electronic basis for metal release and surface oxygen reactivity that form a rationale for complex metal oxide design principles.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Fundamental Interactions of Bimetallic Cu x Pd y ( x + y = 4) Clusters Supported on the α-WC(0001) Surface and Their Performance for CO 2 Adsorption and Dissociation

The tungsten carbide α-WC(0001) surface, an active system for the activation of H 2 and important hydrogenation processes involving unsaturated hydrocarbons, can serve as a support of bimetallic clusters to produce materials with unique catalytic properties, opening routes for a wide range of technical applications. In particular, Cu x Pd y clusters are of particular interest because they combine metals with different properties. A stochastic method was applied to obtain the geometry of Cu x Pd y (x + y = 4) bare clusters, evaluating thousands of possibilities to obtain stable structures, yielding one isomer for Cu 4 , Cu 2 Pd 2 , Cu 1 Pd 3 , and Pd 4 and two isomers for Cu 3 Pd 1 . These clusters were supported on C and W terminations of the tungsten carbide (0001) surface, exploring all of the binding possibilities. The adsorption energies on the C and W terminations are in the ranges from −2.51 to −3.02 eV and from −2.26 to −3.30 eV, respectively. The strongest and weakest binding was seen for monometallic Cu 4 and Pd 4 clusters on both C and W terminations, while the Cu-Pd bimetallics have intermediate adsorption energies but lack a clear trend in terms of composition. The location of Cu x Pd y clusters over the (0001) surface induces a decrease in the work function relative to the pristine surface, while the cluster-surface Bader charge transfer and variations in the partial density of states point to changes in the electronic structure of the carbide atoms upon binding of the metallic clusters. The d-band center of the Cu x Pd y deposited on WC(0001) indicates an intermediate reactivity among Cu(111) and Pd(111) surfaces, modulating the reactivity with small numbers of Cu and Pd atoms, i.e., atom economy in catalyst design. The likelihood of existence of the most stable Cu x Pd y (x + y = 4) clusters in the temperature range of 298-400 K is 100%. The composite Cu x Pd y /α-WC(0001) (x + y = 4), is a nontrivial system since 22 isomers are needed to completely describe its structural properties. Among the isomers, seven structures are necessary to represent Cu 3 Pd 1 /α-WC(0001), five for Pd 4 /α-WC(0001), two for Cu 4 /α-WC(0001), and four for Cu 2 Pd 2 /α-WC(0001) and Cu 1 Pd 3 /α-WC(0001). The large number of cluster isomers supported on the tungsten carbide surface opens the door for several applications in the heterogeneous catalysis of the Cu x Pd y /α-WC(0001) composite, with the possibility of modulating the geometric, electronic, and chemical properties according to a desired application. Test studies for the adsorption of CO 2 indicate that the Cu x Pd y /α-WC(0001) composites are highly active for the adsorption and decomposition of the molecule, with bimetallic and admetal-carbide interactions playing a key role in the binding performance. In conclusion, this high activity indicates that these systems should be useful as catalysts for the conversion of CO 2 to oxygenates or light alkanes.

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

Electronic Structure of the Au Nanoparticle-TiO 2 Heterojunction: Influence of Nanoparticle Size, Shape, Oxygen Vacancies, and Temperature

The electronic structure of the gold nanoparticle-titanium dioxide (AuNP-TiO 2 ) heterojunction plays a critical role for charge transfer and recombination dynamics that underpin its photocatalytic function. However, building a representative model to capture the key physics remains a significant challenge. Here, we investigate the influence of the AuNP size and shape, as well as oxygen vacancy (V O ) defects at the anatase-phase TiO 2 (101) surface and the temperature of the heterojunction, on its interfacial electronic properties. Using density functional theory (DFT), we compare the closed-shell Au 20 and open-shell Au 19 clusters interfaced with pristine and V O defect TiO 2 surfaces. We find that the presence of a V O defect transforms pure TiO 2 from a p-type to an n-type semiconductor, reversing the interfacial band bending from downward to upward. For the heterosystem, density of states (DOS) analysis shows that V O minimally affects the open-shell Au 19 -TiO 2 system, but it significantly alters the closed-shell Au 20 -TiO 2 system, converting it from p-type to n-type at the Γ-point at 0 K. Furthermore, ab initio molecular dynamics (AIMD) simulations at 300 K reveal significant thermal fluctuations in AuNP positions relative to the TiO 2 surface. These fluctuations result in dynamic variations in the gap between the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) in all systems studied: Au 19 NP-Pristine TiO 2 , Au 20 NP-Pristine TiO 2 , Au 19 NP-V O defect TiO 2 , and Au 20 NP-V O defect TiO 2 . Complementary Bader charge analysis performed at both 0 K and finite-temperature AIMD snapshots supports the emergence of an upward band bending and the formation of a Schottky barrier at the V O containing heterojunctions. Notably, we find that at finite temperature, an Au atom can dynamically passivate the V O , leading to a pronounced widening of the HOMO-LUMO gap in the AuNP-V O defect TiO 2 heterojunctions. Our computational findings underscore the pivotal role of the V O defect and thermal effects in modulating interfacial band alignment, electronic states, and HOMO-LUMO gaps, providing insights for designing Au-TiO 2 heterojunctions with tailored electronic properties.

Mora Perez, Carlos [Lawrence Berkeley National Lab↗
DOE OSTI2020

A Metastable CaSH 3 Phase Composed of HS Honeycomb Sheets that is Superconducting Under Pressure

Evolutionary searches have predicted a number of ternary Ca–S–H phases that could be synthesized at pressures of 100–300 GPa. P6 3 /mmc CaSH 2 , Pnma CaSH 2 , Cmc2 1 CaSH 6 , and I$\bar{4}$CaSH 20 were composed of a Ca–S lattice along with H 2 molecules coordinated in a “side-on” fashion to Ca. The H–H bond lengths in these semiconducting phases were elongated because of H 2 σ → Ca d donation, and Ca d → H 2 σ* back-donation, via a Kubas-like mechanism. P$\bar{6}$m2 CaSH 3 , consisting of two-dimensional HS and CaH 2 sheets, was metastable and metallic above 128 GPa. The presence of van Hove singularities increased its density of states at the Fermi level and concomitantly the superconducting critical temperature, which was estimated to be as high as ~100 K at 128 GPa. Furthermore, this work will inspire the search for superconductivity in materials based upon honeycomb HX (X = S, Se, Te) and MH 2 (M = Mg, Ca, Sr, Ba) layers under pressure.

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

Spatially Resolved Investigation of Mixed Valence and Insulator-to-Metal Transition in an Organic Salt

Using scanning tunneling microscopy/spectroscopy (STM/STS), we investigate the evolution of electronic structures across the boundaries of 7,7,8,8-tetracyanoquinodimethane (TCNQ) and K-TCNQ assemblies on a weakly interacting substrate. Despite the semiconducting/insulating nature of TCNQ (TCNQ 0 ) and K-TCNQ (TCNQ -1 ), a continuum metallic-like density of states extending deep (~1.5 nm) into the TCNQ assembly is observed near domain boundary. We attribute the formation of these states to the abrupt change of molecular valence, which perturbs the electrostatics of the junction and creates local electric fields as evidenced by the band bending near domain boundary. To the best of our knowledge, this study provides the first microscopic understanding of the crucial physics occurring near domain boundaries of mixed valence in K-TCNQ, or broadly speaking charge transfer complexes, which highlights these boundaries as potential “weak” points to initiate the electric field induced insulator-to-metal transition.

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

Hydrogen Dynamics in Supercritical Water Probed by Neutron Scattering and Computer Simulations

In this work, an investigation of supercritical water is presented combining inelastic and deep inelastic neutron scattering experiments and molecular dynamics simulations based on a machine-learned potential of ab initio quality. The local hydrogen dynamics is investigated at 250 bar and in the temperature range of 553–823 K, covering the evolution from subcritical liquid to supercritical gas-like water. The evolution of libration, bending, and stretching motions in the vibrational density of states is studied, analyzing the spectral features by a mode decomposition. Moreover, the hydrogen nuclear momentum distribution is measured, and its anisotropy is probed experimentally. It is shown that hydrogen bonds survive up to the higher temperatures investigated, and we discuss our results in the framework of the coupling between intramolecular modes and intermolecular librations. Results show that the local potential affecting hydrogen becomes less anisotropic within the molecular plane in the supercritical phase, and we attribute this result to the presence of more distorted hydrogen bonds.

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

Isolated [B 2 (CN) 6 ] 2– : Small Yet Exceptionally Stable Nonmetal Dianion

Here, we report the observation of a small, yet remarkably stable, metal-free hexacyanodiborate dianion [B 2 (CN) 6 ] 2– in the gas phase. Negative ion photoelectron spectroscopy (NIPES) was employed to measure its spectra at multiple laser wavelengths, yielding a 1.9 eV electron binding energy (EBE) —a remarkably high value of electronic stability and a ~2.60 eV repulsive Coulomb barrier (RCB) for electron detachment. This rationalizes the observation of this dianion, although homolytic charge-separation dissociation into two [B(CN) 3 ] •– is energetically favorable. Quantum chemical calculations demonstrate a $D_{3d}$ staggered conformation for both the dianion and radical monoanion, and the calculated EBE and RCB match the experimental values well. The simulated density of states spectrum reproduces all measured electronic transitions, while the simulated vibrational progressions for the ground state transition cover a much narrower EBE range compared to the experimental band, indicating appreciable auto-photodetachment via electronically excited dianion resonances

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

Spectral Signatures of a Negative Polaron in a Doped Polymer Semiconductor: Energy Levels and Hubbard U Interactions

The modern picture of negative charge carriers on conjugated polymers invokes the formation of a singly occupied (spin-up/spin-down) level within the polymer gap and a corresponding unoccupied level above the polymer conduction band edge. The energy splitting between these sublevels is related to on-site Coulomb interactions between electrons, commonly termed Hubbard U. However, spectral evidence for both sublevels and experimental access to the U value is still missing. Here, in this work, we provide evidence by n-doping the polymer P(NDI 2 OD-T 2 ) with [RhCp*Cp] 2 , [N-DMBI] 2 , and cesium. Changes in the electronic structure after doping are studied with ultraviolet photoelectron and low-energy inverse photoemission spectroscopies (UPS, LEIPES). UPS data show an additional density of states (DOS) in the former empty polymer gap while LEIPES data show an additional DOS above the conduction band edge. These DOS are assigned to the singly occupied and unoccupied sublevels, allowing determination of a U value of ~1 eV.

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

Transferable Water Potentials Using Equivariant Neural Networks

Machine learning interatomic potentials (MLIPs) have emerged as a technique that promises quantum theory accuracy for reduced cost. It has been proposed [J. Chem. Phys. 2023, 158, 084111] that MLIPs trained on solely liquid water data cannot accurately transfer to the vapor–liquid equilibrium while recovering the many-body decomposition (MBD) analysis of gas-phase water clusters. This suggests that MLIPs do not directly learn the physically correct interactions of water molecules, limiting transferability. In this work, we show that MLIPs using equivariant architecture and trained on 3200 liquid water structures reproduces liquid-phase water properties (e.g., density within 0.003 g/cm 3 between 230 and 365 K), vapor–liquid equilibrium properties up to 550 K, the MBD analysis of gas-phase water cluster up to six-body interactions, and the relative energy and the vibrational density of states of ice phases. We show that potentials developed using equivariant MLIPs allow transferability for arbitrary phases of water that remain stable in nanosecond long simulations.

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

Electron–Phonon and Spin–Lattice Coupling in Atomically Thin Layers of MnBi 2 Te 4

MnBi 2 Te 4 represents a new class of magnetic topological insulators in which novel quantum phases emerge at temperatures higher than those found in magnetically doped thin films. Here in this paper, we investigate how couplings between electron, spin, and lattice are manifested in the phonon spectra of few-septuple-layer thick MnBi 2 Te 4 . After categorizing phonon modes by their symmetries, we study the systematic changes in frequency, line width, and line shape of a spectrally isolated A1g mode. The electron–phonon coupling increases in thinner flakes as manifested in a broader phonon line width, which is likely due to changes of the electron density of states. In 4- and 5-septuple thick samples, the onset of magnetic order below the Néel temperature is concurrent with a transition to an insulating state. We observe signatures of a reduced electron–phonon scattering across this transition as reflected in the reduced Fano parameter. Finally, spin–lattice coupling is measured and modeled from temperature-dependent phonon frequency.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Hole Relaxation Bottlenecks in CdSe/CdTe/CdSe Lateral Heterostructures Lead to Bicolor Emission

Concentric lateral CdSe/CdTe/CdSe heterostructures show bicolor photoluminescence from both a red charge transfer band of the CdSe/CdTe interface and a green fluorescence from CdSe. Further, this work uses visible and near-infrared transient spectroscopy measurements to demonstrate that the deviation from Kasha's rule arises from a hole relaxation bottleneck from CdSe to CdTe. Hole transfer can take up to 1 ns, which permits radiative relaxation of excitons remaining in CdSe. Simulations indicate that the hole relaxation bottleneck arises due to the sparse density of states and poor spatial overlap of hole states at energies near the CdSe band edge. The divergent kinetics of transfer for band edge and hot holes is exploited to vary the ratio of green and red photoluminescence with excitation wavelength, providing another knob to control emission color. These findings support the use of lateral heterojunctions as a method for slowing carrier relaxation in two-dimensional materials.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2024

Microscopic Origin of Twist-Dependent Electron Transfer Rate in Bilayer Graphene

Using molecular simulation and continuum dielectric theory, we consider how electrochemical kinetics are modulated by the twist angle in bilayer graphene electrodes. By establishing a connection between the twist angle and the screening length of charge carriers within the electrode, we investigate how tunable metallicity modifies the statistics of the electron transfer energy gap. Constant potential molecular simulations show that the activation free energy for electron transfer increases with screening length, leading to a non-monotonic dependence on the twist angle. Here, the twist angle alters the density of states, tuning the number of thermally accessible channels for electron transfer and the reorganization energy by affecting the stability of the vertically excited state through attenuated image charge interactions. Understanding these effects allows us to express the Marcus rate of interfacial electron transfer as a function of the twist angle in a manner consistent with experimental observations.

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

Nano-photocurrent Mapping of Local Electronic Structure in Twisted Bilayer Graphene

In this work, we report a combined nano-photocurrent and infrared nanoscopy study of twisted bilayer graphene (TBG) enabling access to the local electronic phenomena at length scales as short as 20 nm. We show that the photocurrent changes sign at carrier densities tracking the local superlattice density of states of TBG. We use this property to identify domains of varying local twist angle by local photothermoelectric effect. Consistent with the photocurrent study, infrared nanoimaging experiments reveal optical conductivity features dominated by twist-angle-dependent interband transitions. Our results provide a fast and robust method for mapping the electronic structure of TBG and suggest that similar methods can be broadly applied to probe electronic inhomogeneities of Moiré superlattices in other van der Waals heterostructures.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Stability of Calcium Ion Battery Electrolytes: Predictions from Ab Initio Molecular Dynamics Simulations

Multivalent batteries, such as magnesium-ion, calcium-ion, and zinc-ion batteries, have attracted significant attention as next-generation electrochemical energy storage devices to complement conventional lithium-ion batteries (LIBs). Among them, calcium-ion batteries (CIBs) are the least explored due to the difficult reversible Ca deposition-dissolution. In this work, we examined the stability of four different Ca salts with weakly coordinating anions and three different solvents commonly employed in existing battery technologies to identify suitable candidates for CIBs. By employing Born-Oppenheimer molecular dynamics (BOMD) simulations on salt-Ca and solvent-Ca interfaces, we find that the tetraglyme solvent and carborane salt are promising candidates for CIBs. Due to the strong reducing nature of the calcium surface, the other salts and solvents readily decompose. Further, we explain the microscopic mechanisms of salt/solvent decomposition on the Ca surface using time-dependent projected density of states, time-dependent charge-transfer plots, and climbing-image nudged elastic band calculations. Collectively, this work presents the first mechanistic assessment of the dynamical stability of candidate salts and solvents on a Ca surface using BOMD simulations, and provides a predictive path toward designing stable electrolytes for CIBs.

Born-Oppenheimer molecular dynamics↗