Search NASA⌕ Search

SEARCH · Search NASA

Results for “density functional theory (DFT) calculations”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Feasibility of Prelithiation in LiFePO 4

Lithium iron phosphate (LiFePO 4 ) is widely applied as the cathode material for the energy storage Li-ion batteries due to its low cost and high cycling stability. However, the low theoretical specific capacity of LiFePO 4 makes its initial capacity loss more concerning. Therefore, lithium compensation by way of prelithiation and applications of sacrificial Li-rich additives in LiFePO 4 is imminent in elevating the energy density and/or prolonging the lifetime of the LiFePO 4 -based Li-ion batteries (LIBs). Prelithiation in LiFePO 4 is herein carried out by electrochemical and chemical methods and its feasibility is proved on the basis of the electrochemical evaluations such as the initial charge capacity and the cycling stability. In addition, the site of the pre-intercalated Li-ions is found via comprehensive physical characterizations and the density functional theory (DFT) calculations. Furthermore, these findings open a new avenue for elevating the energy density and/or prolonging the lifetime of the high-energy-density batteries.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Mapping the Complete Reaction Energy Landscape of a Metal–Organic Framework Phase Transformation

Crystalline materials undergo valuable phase transformations, and the energetic processes that underlie these transformations can be fully characterized through a combination of thermodynamic and kinetic studies. Here, we report the first complete reaction energy landscape of metal–organic framework (MOF) interpenetration, specifically in the phase transformation of NU-1200 to its doubly interpenetrated counterpart, STA-26. We characterized the thermodynamics of this phase transformation by pairing experiments with density functional theory (DFT) calculations. This analysis revealed that factors such as the increase in crystal density likely drive Zr- and Hf-NU-1200 to STA-26 interpenetration, while other chemical interactions such as steric repulsions prevent Th-NU-1200 from interpenetrating. Using time-resolved in situ X-ray diffraction, we monitored phase transformation reaction profiles and extracted quantitative kinetic information using the Avrami-Erofe’ev model. As a result, we obtained activation energies for the Zr- and Hf-NU-1200 transformations to Zr- and Hf-STA-26, respectively, revealing slower phase change kinetics for MOFs with stronger bonds. Finally, we paired the kinetic data with experimental observations to classify the mechanistic model of this phase transformation as partial dissolution. Here, we anticipate that this thermodynamic, kinetic, and mechanistic understanding will broadly inform further studies on the energetics of crystallization.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
NASA NTRS2013

Association Mechanisms of Unsaturated C2 Hydrocarbons with Their Cations: Acetylene and Ethylene

The ion-molecule association mechanism of acetylene and ethylene with their cations is investigated by ab initio quantum chemical methods to understand the structures, association energies, and the vibrational and electronic spectra of the products. Stable puckered cyclic isomers are found as the result of first forming less stable linear and bridge isomers. The puckered cyclic complexes are calculated to be strongly bound, by 87, 35 and 56 kcal/mol for acetylene-acetylene cation, ethylene-ethylene cation and acetylene-ethylene cation, respectively. These stable complexes may be intermediates that participate in further association reactions. There are no association barriers, and no significant inter-conversion barriers, so the initial linear and bridge encounter complexes are unlikely to be observable. However, the energy gap between the bridged and cyclic puckered isomers greatly differs from complex to complex: it is 44 kcal/mol in C4H4 +, but only 6 kcal/mol in C4H8 +. The accurate CCSD(T) calculations summarized above are also compared against less computationally expensive MP2 and density functional theory (DFT) calculations for structures, relative energies, and vibrational spectra. Calculated vibrational spectra are compared against available experiments for cyclobutadiene cation. Electronic spectra are also calculated using time-dependent DFT.

ethylene↗
DOE OSTI2020

Superconductivity enhancement in phase-engineered molybdenum carbide/disulfide vertical heterostructures

Stacking layers of atomically thin transition-metal carbides and two-dimensional (2D) semiconducting transition-metal dichalcogenides, could lead to nontrivial superconductivity and other unprecedented phenomena yet to be studied. In this work, superconducting α-phase thin molybdenum carbide flakes were first synthesized, and a subsequent sulfurization treatment induced the formation of vertical heterolayer systems consisting of different phases of molybdenum carbide—ranging from α to γ' and γ phases—in conjunction with molybdenum sulfide layers. These transition-metal carbide/disulfide heterostructures exhibited critical superconducting temperatures as high as 6 K, higher than that of the starting single-phased α-Mo 2 C (4 K). We analyzed possible interface configurations to explain the observed moiré patterns resulting from the vertical heterostacks. Our density-functional theory (DFT) calculations indicate that epitaxial strain and moiré patterns lead to a higher interfacial density of states, which favors superconductivity. Such engineered heterostructures might allow the coupling of superconductivity to the topologically nontrivial surface states featured by transition-metal carbide phases composing these heterostructures potentially leading to unconventional superconductivity. Moreover, we envisage that our approach could also be generalized to other metal carbide and nitride systems that could exhibit high-temperature superconductivity.

heterostructure↗
DOE OSTI2024

2D-EFICACY: Control of Metastable 2D Carbide[1]Chalcogenide Heterolayers: Strain and Moire Engineering

The experimental isolation of graphene led to the discovery of an entirely new world of two-dimensional (2D) materials in which the 2D nature often leads to emergent behaviors not seen in bulk systems. 2D transition metal dichalcogenides (TMDs) exhibit physico-chemical properties that depend on the transition metal, polymorph, thickness, and presence and type of defects. Recently, a group of thin (10-100nm) transition metal carbides (TMCs), such as Mo2C, has been synthesized that exhibit a thickness-dependent superconducting critical temperature (Tc). These thin TMCs are different from MXenes, another class of 2D materials consisting of few layers of nitrides or carbides (<5nm) produced by chemical etching and delamination. The goal of this renewal proposal is to combine experiment and computation to synthesize and elucidate the guiding principles that control the growth, orientation and strain of heterostacks of thin TMCs and TMDs composed with Nb, Ti and W. We expect to stabilize metastable hybrid phases of TMCs sandwiched between TMDs (H-TMD/Cs) with unprecedented physico-chemical properties. As part of previous DOE-funded work by the Terrones/Sinnott groups, thin (10-100 nm thick) Mo2C flakes were successfully synthesized by chemical vapor deposition (CVD). By subsequently exposing Mo2C to H2S, partial chalcogenization was achieved, resulting in heterostacks of MoCx phases and MoS2. The formation of MoS2 led to a deficiency of Mo atoms in the underlying Mo2C, resulting in an inhomogeneous phase change from α-Mo2C to γ’-MoCx and then to γ-MoC. The γ’-MoCx is a strained metastable phase and the heterostack of all three phases demonstrated an increased Tc relative to that of α-Mo2C, from 4 to 6K; its interleaved layered structure consisting of superconducting and semiconducting phases is ideal for future studies of Josephson junction series arrays. Moiré patterns in these heterostacked systems could result in new phenomena, as moiré patterns in bilayer graphene showed unconventional superconductivity and moiré excitons have been observed in twisted TMD heterobilayers. The scientific hypothesis of the proposed synergistic computational and experimental research is that orientation and strain control within confined thin metastable TMCs, sandwiched by stable phases of TMCs and layered TMDs, will depend on kinetic and thermodynamic “knobs” that include fast temperature changes, chalcogen diffusion through preferred crystallographic planes, reaction times, pressure, reactive atmosphere, precursors, and surfactants, which will also tailor properties such as superconductivity, magnetism, ferroelectricity, piezoelectricity, and catalytic performance. We will develop the guiding principles for the synthesis and stabilization of metastable H-TMD/Cs based on Nb, Ti and W. In order to validate the hypothesis, four tasks are proposed: The first task will synthesize ultra-thin TMCs based on Nb, W and Ti, by: 1) adapting the CVD method used for Mo2C, 2) plasma assisted CVD, 3) defect-mediated CVD processes, and 4) cryo-milling of carbide powders. The second task will accomplish the synthesis and basic physico-chemical characterizations of H-TMD/Cs by chalcogenization of the materials synthesized in task one, and by carbonization of TMDs. H-TMD/Cs will also be investigated for their suitability in energy conversion applications such as supercapacitors, Li and multivalent ion batteries, and electrocatalysts, topics of interest to DoE. These tasks will be carried out in close conjunction with density functional theory (DFT) calculations with insights into energetics, lattice parameters, stability, phase diagrams, band structures, and density of states of H-TMD/Cs. The third task will characterize and evaluate strain and moiré patterns at the interfaces of different H-TMD/Cs by high-resolution scanning transmission electron microscopy (HR-STEM), scanning tunneling microscopy (STM), and conductive tip atomic force microscopy. Nudged elastic band calculations with DFT will be performed to understand the chalcogen diffusion process, which will provide insights into the interfaces between different phases of TMCs and TMDs. The fourth task aims at quantifying the stability and dynamics of H-TMD/Cs by in-situ TEM and Raman studies under heating, strain, and electrical biasing. Phonon calculations using DFT will provide a basis for interpreting Raman spectra. This coherent framework involving synthesis, characterization, and computation will result in a broad scientific impact for energy related applications. The ability to develop new H-TMD/Cs will enhance a range of applications that include batteries, catalysts, switches, sensors, quantum computing components and smart coatings.

2-Dimensional materials↗
DOE OSTI2024

2D-EFICACY: Control of Metastable 2D Carbide Chalcogenide Heterolayers: Strain and Moire Engineering

The experimental isolation of graphene led to the discovery of an entirely new world of two-dimensional (2D) materials in which the 2D nature often leads to emergent behaviors not seen in bulk systems. 2D transition metal dichalcogenides (TMDs) exhibit physico-chemical properties that depend on the transition metal, polymorph, thickness, and presence and type of defects. Recently, a group of thin (10-100nm) transition metal carbides (TMCs), such as Mo 2 C, has been synthesized that exhibit a thickness-dependent superconducting critical temperature (Tc). These thin TMCs are different from MXenes, another class of 2D materials consisting of few layers of nitrides or carbides (<5nm) produced by chemical etching and delamination. The goal of this renewal proposal is to combine experiment and computation to synthesize and elucidate the guiding principles that control the growth, orientation and strain of heterostacks of thin TMCs and TMDs composed with Nb, Ti and W. We expect to stabilize metastable hybrid phases of TMCs sandwiched between TMDs (H-TMD/Cs) with unprecedented physico-chemical properties. As part of previous DOE-funded work by the Terrones/Sinnott groups, thin (10-100 nm thick) Mo 2 C flakes were successfully synthesized by chemical vapor deposition (CVD). By subsequently exposing Mo 2 C to H 2 S, partial chalcogenization was achieved, resulting in heterostacks of MoCx phases and MoS 2 . The formation of MoS 2 led to a deficiency of Mo atoms in the underlying Mo 2 C, resulting in an inhomogeneous phase change from α-Mo 2 C to γ’-MoCx and then to γ-MoC. The γ’-MoCx is a strained metastable phase and the heterostack of all three phases demonstrated an increased Tc relative to that of α-Mo 2 C, from 4 to 6K; its interleaved layered structure consisting of superconducting and semiconducting phases is ideal for future studies of Josephson junction series arrays. Moiré patterns in these heterostacked systems could result in new phenomena, as moiré patterns in bilayer graphene showed unconventional superconductivity and moiré excitons have been observed in twisted TMD heterobilayers. The scientific hypothesis of the proposed synergistic computational and experimental research is that orientation and strain control within confined thin metastable TMCs, sandwiched by stable phases of TMCs and layered TMDs, will depend on kinetic and thermodynamic “knobs” that include fast temperature changes, chalcogen diffusion through preferred crystallographic planes, reaction times, pressure, reactive atmosphere, precursors, and surfactants, which will also tailor properties such as superconductivity, magnetism, ferroelectricity, piezoelectricity, and catalytic performance. We will develop the guiding principles for the synthesis and stabilization of metastable H-TMD/Cs based on Nb, Ti and W. In order to validate the hypothesis, four tasks are proposed: The first task will synthesize ultra-thin TMCs based on Nb, W and Ti, by: 1) adapting the CVD method used for Mo 2 C, 2) plasma assisted CVD, 3) defect-mediated CVD processes, and 4) cryo-milling of carbide powders. The second task will accomplish the synthesis and basic physico-chemical characterizations of H-TMD/Cs by chalcogenization of the materials synthesized in task one, and by carbonization of TMDs. H-TMD/Cs will also be investigated for their suitability in energy conversion applications such as supercapacitors, Li and multivalent ion batteries, and electrocatalysts, topics of interest to DoE. These tasks will be carried out in close conjunction with density functional theory (DFT) calculations with insights into energetics, lattice parameters, stability, phase diagrams, band structures, and density of states of H-TMD/Cs. The third task will characterize and evaluate strain and moiré patterns at the interfaces of different H-TMD/Cs by high-resolution scanning transmission electron microscopy (HR-STEM), scanning tunneling microscopy (STM), and conductive tip atomic force microscopy. Nudged elastic band calculations with DFT will be performed to understand the chalcogen diffusion process, which will provide insights into the interfaces between different phases of TMCs and TMDs. The fourth task aims at quantifying the stability and dynamics of H-TMD/Cs by in-situ TEM and Raman studies under heating, strain, and electrical biasing. Phonon calculations using DFT will provide a basis for interpreting Raman spectra. This coherent framework involving synthesis, characterization, and computation will result in a broad scientific impact for energy related applications. The ability to develop new H-TMD/Cs will enhance a range of applications that include batteries, catalysts, switches, sensors, quantum computing components and smart coatings.

2-Dimensional Materials↗
DOE OSTI2023

Conventional High-Temperature Superconductivity in Metallic, Covalently Bonded, Binary-Guest C–B Clathrates

Inspired by the synthesis of XB 3 C 3 (X = Sr, La) compounds in the bipartite sodalite clathrate structure, density functional theory (DFT) calculations are performed on members of this family containing up to two different metal atoms. A DFT chemical pressure analysis on systems with X = Mg, Ca, Sr, Ba reveals that the size of the metal cation, which can be tuned to stabilize the B–C framework, is key for their ambient-pressure dynamic stability. High-throughput density functional theory calculations on 105 Pm$\bar{3}$ symmetry XYB 6 C 6 binary-guest compounds (where X, Y are electropositive metal atoms) find 22 that are dynamically stable at 1 atm, expanding the number of potentially synthesizable phases by 19 (18 metals and 1 insulator). Here, the density of states at the Fermi level and superconducting critical temperature, T c , can be tuned by changing the average oxidation state of the metal atoms, with T c being highest for an average valence of +1.5. KPbB 6 C 6 , with an ambient-pressure Eliashberg T c of 88 K, is predicted to possess the highest T c among the studied Pm$\bar{3}$n XB 3 C 3 or Pm$\bar{3}$XYB 6 C 6 phases, and calculations suggest it may be synthesized using high-pressure high-temperature techniques and then quenched to ambient conditions.

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

Mechanistic Study of 1,2-Dichloroethane Hydrodechlorination on Cu-Rich Pt–Cu Alloys: Combining Reaction Kinetics Experiments with DFT Calculations and Microkinetic Modeling

Cu-rich Pt–Cu bimetallic catalysts are among the most promising candidates for actively catalyzing the hydrodechlorination of 1,2-dichloroethane (1,2-DCA) toward ethylene production. Combining reaction kinetics experiments with density functional theory (DFT) calculations and mean-field microkinetic modeling, we present a systematic mechanistic study for 1,2-DCA hydrodechlorination on Cu-rich Pt–Cu alloy catalysts. Our DFT (PBE+(TS+SCS)) results suggest that increasing Cu content in the Pt–Cu alloy destabilizes C 2 -species adsorption while stabilizing the binding of atomic chlorine. The reaction energetics of all the elementary steps in the 1,2-DCA reaction network were calculated on a Pt 1 Cu 3 (111) model surface. The DFT results were then used to construct a microkinetic model, and the model-predicted reaction rates were compared with our reaction kinetics experimental results on a Cu-rich SiO 2 -supported Pt–Cu alloy catalyst through a parameter estimation procedure. Both the reaction kinetics experiments and the microkinetic model after parameter adjustments yielded 100% selectivity to ethylene. The microkinetic model pointed to a reaction pathway involving two sequential chlorine-removal steps on the Pt–Cu alloy catalyst, a mechanism distinct from the one previously identified on pure Pt/SiO 2 catalysts, which involved an initial hydrogen-removal step. Adjustments to the DFT-derived parameters indicate the possible formation of chlorine-induced Cu-enriched surface sites during 1,2-DCA hydrodechlorination conditions, sites that are more active than those encountered in the bulk Pt 1 Cu 3 (111) alloy surface. Furthermore, our study offers valuable initial insights on the 1,2-DCA hydrodechlorination reaction mechanism and the nature of the active sites on PtCu bimetallic catalysts.

1,2-Dichloroethane↗
DOE OSTI2021

Surfactant inhibition mechanisms of carbonate mineral dissolution in shale

Surfactants are common additives to hydraulic fracturing and enhanced oil recovery (EOR) fluids, and are under consideration for amendment to supercritical carbon dioxide for geological carbon sequestration (GCS). The effect of a common anionic surfactant, internal olefin sulfonate (IOS), on mineral dissolution from shale into brine was evaluated. When added to brine at concentrations exceeding the critical micelle concentration (94 mg/L), IOS inhibited carbonate mineral dissolution in an Eagle Ford shale, as well as dissolution of optical quality calcite (the dominant carbonate in the shale). Laser profilometry images provide spatial resolution across > 3 orders of magnitude, and indicate that IOS addition to brine both enhances the formation of new etch pits in calcite, and impedes their further growth. Time-of-flight secondary ion mass spectrometry surface profiles show for the first time that IOS preferentially adsorbs at calcite pit edges versus flat calcite surfaces (i.e., terraces). Surface pressure calculations, sulfur K-edge near edge X-ray absorption fine structure (NEXAFS) spectroscopy results, and density functional theory (DFT) calculations support this observation; the DFT results indicate that the sulfonate head group of the IOS molecule binds strongly to the calcite step site as compared to the terrace site. The S K-edge NEXAFS results indicate that IOS adsorbed more to etched calcite surfaces compared to smooth calcite surfaces. Overall, the results indicate that weak adsorption on flat calcite surfaces (i.e., terraces) disrupts water structure and enhances mass transfer of dissolution, while strong adsorption on calcite pit edges displaces adsorbed water and inhibits further etch pit growth. This work provides the first direct evidence of preferential adsorption of IOS to etched calcite surfaces and links it to macroscopic dissolution kinetics. Finally, this work has implications for surfactant-containing fluids used in hydraulic fracturing, EOR and potentially GCS for subsurface injection into carbonate rich reservoirs.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Coexistence of Dirac fermion and charge density wave in the square-net-based semimetal LaAuSb 2

Here, we report a comprehensive study of magnetotransport properties, angle-resolved photoemission spectroscopy (ARPES), and density functional theory (DFT) calculations on self-flux grown LaAuSb 2 single crystals. Resistivity and Hall measurements reveal a charge density wave (CDW) transition at 77 K. MR and de Haas-Van Alphen (dHvA) measurements indicate that the transport properties of LaAuSb 2 are dominated by Dirac fermions that arise from Sb square nets. ARPES measurements and DFT calculations reveal an electronic structure with a common feature of the square-net-based topological semimetals, which is in good agreement with the magnetotransport properties. Our results indicate the coexistence of CDW and Dirac fermion in LaAuSb 2 , both of which are linked to the bands arising from the Sb-square net, suggesting that the square net could serve as a structural motif to explore various electronic orders.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Atomic Manipulation of Single Erbium on GaAs(110): From Adatom Motion to Creation and Control of Bistable Substitutions

We demonstrate atomic-scale control over the position and configuration of single erbium (Er) atoms at the GaAs(110) surface using scanning tunneling microscopy (STM). By applying voltage pulses from the STM tip, we drive Er adatoms to move sequentially among three distinct adsorption sites and controllably substitute into Ga lattice sites. The resulting substitutional Er atoms exhibit two isoelectronic configurations that differ by less than one angstrom in lateral position and can be reversibly switched by additional voltage pulses. Differential conductance spectroscopy shows differences in the unoccupied local electronic states for the two substitutional configurations, and an absence of mid-gap states. Preliminary density functional theory (DFT) calculations support these findings, identifying the most stable adsorption and substitutional sites, their relative binding energies, charge states, and local density of states (LDOS) characteristics. Together, the experimental and theoretical results establish a method for precise, reversible manipulation of rare-earth dopants at semiconductor surfaces, opening a pathway to engineer atomic-scale quantum and optoelectronic functionalities.

Poche, Steven [Ohio State University]↗
DOE OSTI2026

An Algorithm for Atom-Centered Lossy Compression of the Atomic Orbital Basis in Density Functional Theory Calculations

Large atomic-orbital (AO) basis sets of at least triple and preferably quadruple-ζ (QZ) size are required to adequately converge Kohn–Sham density functional theory (DFT) calculations toward the complete basis set limit. However, incrementing the cardinal number by one nearly doubles the AO basis dimension, and the computational cost scales as the cube of the AO dimension, so this is very computationally demanding. Here, in this work, we develop and test a threshold-based natural atomic orbital (NAO) scheme in which ϵ-NAOs are obtained as eigenfunctions of atomic blocks of the density matrix in a one-center orthogonalized representation. This enables compression of the AO basis that is optimal for a given threshold, 10 –ϵ , by discarding NAOs with occupation numbers below that threshold. Extensive pilot test calculations using the Hartree–Fock functional and taking the converged density matrix as input suggest that a threshold of 10 –5 can yield a compression factor (ratio of AO to compressed ϵ-NAO dimension) between 2.5 and 4.5 for the QZ pc-3 basis. The errors in relative energies are typically less than 0.1 kcal/mol when the compressed basis is used instead of the uncompressed basis. Between 10 and 100 times smaller errors (i.e., usually less than 0.01 kcal/mol) can be obtained with a threshold 10 –7 , while the compression factor is typically between 2 and 2.5.

basis sets↗
DOE OSTI2020

Transition Metal Nitrides as Promising Catalyst Supports for Tuning CO/H 2 Syngas Production from Electrochemical CO 2 Reduction

The electrochemical carbon dioxide reduction reaction (CO 2 RR) to produce synthesis gas (syngas) with tunable CO/H 2 ratios has been studied by supporting Pd catalysts on transition metal nitride (TMN) substrates. Combining experimental measurements and density functional theory (DFT) calculations, Pd-modified niobium nitride (Pd/NbN) is found to generate much higher CO and H 2 partial current densities and greater CO Faradaic efficiency than Pd-modified vanadium nitride (Pd/VN) and commercial Pd/C catalysts. In-situ X-ray diffraction identifies the formation of PdH in Pd/NbN and Pd/C under CO 2 RR conditions, whereas the Pd in Pd/VN is not fully transformed into the active PdH phase. DFT calculations show that the stabilized *HOCO and weakened *CO intermediates on PdH/NbN are critical to achieving higher CO 2 RR activity. This work suggests that NbN is a promising substrate to modify Pd, resulting in an enhanced electrochemical conversion of CO 2 to syngas with a potential reduction in precious metal loading.

25 ENERGY STORAGE↗
DOE OSTI2022

Real-Time Modulation of Hydrogen Evolution Activity of Graphene Electrodes Using Mechanical Strain

This paper reports the effect of mechanically applied elastic strain on the hydrogen evolution reaction (HER) activity of graphene under acidic conditions. Here, applied tensile strain of 0.2% on a graphene electrode is shown to lead to a 1–3% increase of HER current. Tensile strain increases HER activity while compressive strain decreases it. Density functional theory (DFT) calculations using a periodic graphene slab model show increasing in the adsorption energy of H atom with growing tensile strain, leading to a corresponding enhancement of current density in HER, similar to that observed experimentally.

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

All-electron APW+ lo calculation of magnetic molecules with the SIRIUS domain-specific package

We report APW+lo (augmented plane wave plus local orbital) density functional theory (DFT) calculations of large molecular systems using the domain specific SIRIUS multi-functional DFT package. The APW and FLAPW (full potential linearized APW) task and data parallelism options and the advanced eigen-system solver provided by SIRIUS can be exploited for performance gains in ground state Kohn–Sham calculations on large systems. This approach is distinct from our prior use of SIRIUS as a library backend to another APW+lo or FLAPW code. We benchmark the code and demonstrate performance on several magnetic molecule and metal organic framework systems. We show that the SIRIUS package in itself is capable of handling systems as large as a several hundred atoms in the unit cell without having to make technical choices that result in the loss of accuracy with respect to that needed for the study of magnetic systems.

Chemistry↗
DOE OSTI2026

Thorium Monosilicide, ThSi: An Experimental and Theoretical Study

The present theoretical and experimental combination study investigates the ThSi molecule in detail. Computationally, we utilized high-level multireference and coupled-cluster levels of theory conjoined with large correlation consistent basis sets to study a series of electronic and spin–orbit states of ThSi. Here, we report potential energy curves (PECs), electron configurations at equilibrium distances, spectroscopic constants, energetics, and spin–orbit coupling effects for 16 electronic states of ThSi. The studied 16 electronic states are arranged tightly within 0.9 eV, highlighting the complexity of the electronic spectrum of ThSi. The ground electronic state of ThSi is a single-reference 1 1 Σ + state that derives from the 1σ 2 2σ 2 1π 4 electronic configuration. The Ω = 0 + spin–orbit ground state of ThSi is composed of 1 1 Σ + (47%) and 13Π (44%) electronic states. Our measured bond energy (D0) of ThSi, obtained using resonant two-photon ionization (R2PI) spectroscopy is 3.146(4) eV, where the assigned error limit is given in parentheses in units of the last quoted digits. The computed D0 of ThSi (Ω = 0 + ) at the CBS-C-CCSD(T)-δT(Q)-δDK-δSO level (3.181 eV) is in good agreement with the experimental value. Our derived enthalpy of formation for ThSi, Δ f H 0K o (ThSi(g)), is 971.8(6.0) kJ/mol. Finally, we have performed density functional theory (DFT) calculations for ThSi(1 1 Σ + ) using 16 exchange correlation functionals that span multiple rungs of “Jacob’s ladder” of density functional approximation (DFA) to assess the DFT errors on D 0 , r e , and ω e of ThSi with respect to experimental and ab initio coupled-cluster values.

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

Spectroscopic Insight on Neodymium Solvation in Lithium Borohydride-Supported Electrolyte

Borohydride-based electrolytes have recently emerged as promising media for the electrodeposition of electropositive metals, including rare earth (RE) elements. While the presence of supporting alkali metal cations and RE counteranions provides essential electrochemical conductivity for achieving fast metal electrodeposition, interactions between the host ligand and solvated neodymium (Nd) complexes remain unclear. This study provides insights into the coordination structure of a concentrated and directly solvated Nd salt in a lithium borohydride-supported electrolyte. Our spectroscopic results indicate that the RE coordination environment is significantly influenced by the solvation mechanism, which can vary between metathesis and complexation pathways, primarily dictated by stoichiometric factors. Under dilute conditions, nearly complete metathesis of anions leads to a high coordination number for the host ligand (borohydride), consistent with the previously reported solvated Nd speciation in chlorine-free electrolytes. In contrast, concentrated dissolution of the Nd salt in the supported electrolyte is dominated by a complexation pathway featuring a Li-ion-paired complex with a low coordination number of the host ligand. Density functional theory (DFT) calculations indicated that the observed blue shift in the borohydride vibration was the result of an increase in electron density drawn into the terminal B–H interbond region from the hydride as the coordination changed from Li to Nd. In conjunction with DFT results, vibrational analyses allowed correlation of the experimental shifts associated with changes in Nd ligation and coordination spheres, further consolidating the prevalence of highly chloride-coordinated species under concentrated conditions. In conclusion, the outcomes of this work illuminate the distinctive and heterogeneous coordination structures that the electroactive RE species can adopt at high concentrations in lithium borohydride-supported electrolytes, as a key step to comprehend the reported metal electrodeposition performance in these media.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗