SEARCH · Engineering Papers

Results for “first principles density functional theory”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 829 records · Page 46

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Unraveling the intrinsic atomic physics behind x-ray absorption line shifts in warm dense silicon plasmas

Here, we present a novel free-energy density functional theory (DFT)-based methodology for optical property calculations of warm dense matter to cover a wide range of thermodynamic conditions and photon energies including the entire x-ray range. It uses Mermin–Kohn–Sham density functional theory with exchange-correlation (XC) thermal effects taken into account via a fully temperature dependent generalized gradient approximation XC functional. The methodology incorporates a combination of the ab initio molecular dynamics (AIMD) snapshotted Kubo–Greenwood optic data with a single-atom in simulation cell calculations to close the photon energy gap between the L and K edges and extend the K-edge tail toward many-keV photon energies. This gap arises in the standard scheme due to a prohibitively large number of bands required for the Kubo–Greenwood calculations with AIMD snapshots. Kubo–Greenwood data on snapshots provide an accurate description of optic properties at low photon frequencies slightly beyond the L-edge and x-ray absorption near edges structure (XANES) spectra, while data from periodic calculations with single atoms cover the tail regions beyond the edges. To demonstrate its applicability to mid-Z materials where the standard DFT-based approach is not computationally feasible, we have applied it to opacity calculations of warm dense silicon plasmas. These first-principles calculations revealed a very interesting phenomon of red-to-blue shift in K–L (1s → 2p) and K-edge absorptions along both isotherm and isochore, which are absent in most continuum-lowering models of traditional plasma physics. This new physics phenomenon can be attributed to the underlying competition between the screening of deeply bound core electrons and the screening of outer-shell electrons caused by warm-dense-plasma conditions. We further demonstrate that the ratio of 1s → 2p to the K-edge x-ray absorptions can be used to characterize warm-dense-plasma conditions. Eventually, based on our absorption calculations, we have established a first-principles opacity table (FPOT) for silicon in a wide range of material densities and temperatures.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Reactivity of Carbonate Solvent Electrolytes on Lithium Silicon Anodes

Silicon (Si) is promising for lithium-ion battery (LIB) anodes due to their high theoretical capacity and low electrochemical potential. However, significant challenges remain, including severe volumetric expansion during cycling and the electrochemical instability of electrolytes, which leads to the formation of a nonuniform solid electrolyte interphase (SEI). To investigate SEI formation mechanisms, computational molecular dynamics simulations offer valuable insights. In this work, we examine the trajectories and charge transfer behavior of lithium hexafluorophosphate (LiPF 6 ) salt with various solvent compositions using density functional theory (DFT) and ab initio molecular dynamics (AIMD). Among the tested electrolyte systems, LiPF 6 with vinylene carbonate (VC) added to ethyl methyl carbonate (EMC) exhibits the lowest reactivity with the Si anode. In contrast, the effects of fluoroethylene carbonate (FEC) and VC depend on whether the primary solvent is EMC alone or a mixture of ethylene carbonate (EC) and EMC. Moreover, we show that electrolyte reactivity varies with the degree of lithiation of the Si anode (LiSi vs Li 15 Si 4 ) and under different charge states. To decouple electrolyte reactivity from surface effects, we analyze the dissociation and formation energies of individual species from solvated configurations. Overall, these first-principles-based findings provide a strategic foundation for electrolyte design to improve cycling stability and extend calendar life in LIBs using Si anodes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Electrochemical Sodiation Mechanism in Magnetite Nanoparticle-Based Anodes: Understanding of Nanoionics-Based Sodium Ion Storage Behavior of Fe 3 O 4

We report electrochemical ion storage behaviors of Fe 3 O 4 nanoparticles, as a representative transition metal oxide for an environmentally benign and low-cost anode for a sodium-ion battery, are thoroughly investigated through a combination of electrochemical analysis and diagnostics of Fe 3 O 4 electrode cells, X-ray-based and spectroscopic analysis of material structure evolution as functions of depth of discharge (DoD) and state of charge (SoC), and first principle modeling. The gravimetric capacity is found to be 50 mA h/g for bulk Fe 3 O 4 (50 nm average crystallite size) and 100 mA h/g--about a tenth of the theoretical prediction for complete conversion--for Fe 3 O 4 nanoparticles (8.7 nm average particle size), respectively. A fundamental and mechanistic study of material evolution as functions of DoD and SoC shows that Fe 3 O 4 does not allow electrochemical incorporation of Na + ions into the empty cation positions of the inverse spinel structure, leading to our assertion that electrochemical intercalation of Na + ions to conversion of the Fe 3 O 4 anode in sodium-ion batteries is nonviable. A density functional theory investigation points to the impracticality of the intercalation of Na + ions into Fe 3 O 4 and further validates our experimental findings. We propose several possible mechanisms corresponding to the observed low capacity, including formation of solid electrolyte interphases with unfavorable properties and adsorption of Na + ions onto surfaces of nanoparticles and/or at heterointerfaces in Fe 3 O 4 composite electrodes in a NaPF 6 -based electrolyte system.

25 ENERGY STORAGE↗

DOE OSTI2024

A First-Principles Approach to Modeling Surface Site Stabilities on Multimetallic Catalysts

The study of multimetallic alloys and the multitude of possible surface compositions have sparked a tremendous interest in engineering low-cost materials with high activity and selectivity in heterogeneous catalysis. Multimetallic systems provide complementary functionalities and an unprecedented tunability when designing catalyst formulations. However, due to their immense structural and compositional complexity, the investigation and identification of an optimal catalyst is a tedious and time-consuming process, both experimentally and theoretically. Therefore, theoretical design principles are highly desirable to accelerate the screening of catalyst structures across the vast compositional space. In this paper, we introduce a simple and general model for predicting the site stability of multimetallic surfaces and nanoparticles, which is based on physical principles. The model requires only a small set of density functional theory (DFT) calculations of metal atom binding energies on monometallic and dilute alloy surface slabs to optimize the parameters in the simple model. The resulting model allows for the quantification of the stability of any particular atom site in any conceivable chemical environment across a wide range of morphologies, sizes, and arrangements by interpolating the derived parameters from a monometallic system to a completely diluted alloyed system. Herein, we demonstrate the robustness of the model across an extensive data set of transition metal alloy surfaces and 147-atoms cuboctahedral nanoparticles (NPs) composed of IrRhRu and PtPdRu. In conclusion, our approach yields mean absolute errors of ≈0.15 (IrRhRu), 0.20 (PtPdRu), 0.19 (IrRhRu NP), and 0.26 (PtPdRu NP) eV relative to site binding energies calculated using DFT.

36 MATERIALS SCIENCE↗

DOE OSTI2022

First-principles study of the T center in silicon

The T center in silicon is a well-known carbon-based color center that has been recently considered for quantum technology applications. Using first-principles computations, we show that the excited state is formed by a defect-bound exciton made of a localized defect state occupied by an electron to which a hole is bound. The localized state is of strong carbon p character and reminiscent of the localization of the unpaired electron in the ethyl radical molecule. The radiative lifetime for the defect-bound exciton is calculated to be on the order of microseconds, much longer than for other quantum defects such as the nitrogen vacancy center in diamond and in agreement with experiments. The longer lifetime is associated with the small transition dipole moment as a result of the very different nature of the localized and delocalized states forming the defect-bound exciton. Finally, we use first-principles calculations to assess the stability of the T center. Here, we find the T center to be stable against decomposition into simpler defects when keeping the stoichiometry fixed. However, we identify that the T center is easily prone to (de)hydrogenation and so requires very precise annealing conditions (temperature and atmosphere) to be efficiently formed.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Disorder-induced magnetoelastic behaviors of MnTexSbyBi1-x-y alloys

This dataset contains input and output files from density functional theory (DFT) simulations used to study the disorder-induced magnetoelastic behaviors of MnTexSbyBi1-x-y (0 ≤ x + y ≤ 1) alloys and their binary end members MnTe, MnSb, and MnBi. The alloys adopt the hexagonal NiAs-type (nickeline) structure and span ternary (MnTexSb1-x, MnTexBi1-x, MnBixSb1-x), and quaternary compositions across the full MnTe–MnSb–MnBi composition triangle. For each alloy composition, the dataset provides DFT calculations in three magnetic configurations: A-type antiferromagnetic (AFM), C-type AFM, and ferromagnetic (FM). Every magnetic configuration folder contains the fully relaxed crystal structure (CONTCAR), VASP input parameters (INCAR), and the main VASP output file (OUTCAR), from which total electronic energies, Mn magnetic moments, lattice parameters, and percent volume changes between magnetic states are extracted. These data are used to construct compositional phase diagrams, evaluate thermodynamic stability (formability), and map magnetoelastic responses across the alloy space. For A-type AFM and FM configurations, additional data are provided as follows: (i) FORCE_CONSTANTS and thermal_properties.yaml files at the top level of A-type_AFM/ and FM/ folders — present only for compositions marked with an asterisk (*) in Table I of the main text. These are derived from Phonopy finite-displacement calculations on full disordered 128-atom supercells and provide vibrational free energy, entropy (Svib)contribution from explicit disorder calculations. (Table I of the associated main manuscript) (ii) A VCA/ subfolder within A-type_AFM/ and FM/, containing FORCE_CONSTANTS and thermal_properties.yaml from Virtual Crystal Approximation phonon calculations (without spin-orbit coupling). VCA data are available for all compositions and are used to estimate vibrational contributions to the Gibbs free energy across the full composition space. (iii) A SOC/ subfolder containing CONTCAR, INCAR, and OUTCAR from spin-orbit coupling calculations, providing relativistic corrections to electronic energies and lattice parameters (Tables S2–S3 of the SM, and Table I of the main manuscript). (iv) A SOC/VCA/ subfolder containing FORCE_CONSTANTS and thermal_properties.yaml from VCA phonon calculations performed within the SOC framework, combining relativistic and vibrational thermodynamic corrections. The computed properties are used to map the AFM–FM magnetic crossover near MnTe0.75Sb0.25, demonstrate disorder- and spin-induced phonon broadening, identify a semiconductor-to-metal crossover, and quantify the pronounced magnetoelastic volume response near the magnetic phase boundary.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Factors influencing quantum evaporation of helium from polar semiconductors from first principles

While there is much indirect evidence for the existence of dark matter (DM), to date it has evaded detection. Current efforts focus on DM masses over ~GeV—to push the sensitivity of DM searches to lower masses, new DM targets and detection schemes are needed. In this work, we focus on the latter—a novel detection scheme recently proposed to detect 10–100 meV phonons in polar target materials. Previous work showed that well-motivated models of DM can interact with polar semiconductors to produce an athermal population of phonons. This new sensing scheme proposes that these phonons then facilitate quantum evaporation of 3 He from a van der Waals film deposited on the target material. However, a fundamental understanding of the underlying process is still unclear, with several uncertainties related to the precise rate of evaporation and how it can be controlled. In this work, we use ab initio density functional theory calculations to compare the adsorption energies of helium atoms on a polar target material, sodium iodide, to understand the underlying evaporation physics. We explore the role of surface termination, monolayer coverage, and elemental species on the rate of He evaporation from the target material. Using this, we discuss the optimal target features for He-evaporation experiments and their range of tunability through chemical and physical modifications such as applied field and surface termination.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Relationship between A-site cation and magnetic structure in 3d–5d–4f double perovskite iridates Ln 2 NiIrO 6 ( Ln = La, Pr, Nd)

We report a comprehensive investigation of Ln 2 NiIrO 6 (Ln=La, Pr, Nd) using thermodynamic and transport properties, neutron powder diffraction, resonant inelastic x-ray scattering, and density-functional theory (DFT) calculations to investigate the role of A-site cations on the magnetic interactions in this family of hybrid 3d–5d–4f compositions. Magnetic structure determination using neutron diffraction reveals antiferromagnetism for La 2 NiIrO 6 , a collinear ferrimagnetic Ni and Ir state that is driven to long-range antiferromagnetism upon the onset of Nd ordering in Nd 2 NiIrO 6 , and a noncollinear ferrimagnetic Ni and Ir sublattice interpenetrated by a ferromagnetic Pr lattice for Pr 2 NiIrO 6 . For Pr 2 NiIrO 6 , heat-capacity results reveal the presence of two independent magnetic sublattices, and transport resistivity indicates insulating behavior and a conduction pathway that is thermally mediated. Here, a first principles DFT calculation elucidates the existence of the two independent magnetic sublattices within Pr 2 NiIrO 6 and offers insight into the behavior in La 2 NiIrO 6 and Nd 2 NiIrO 6 . Resonant inelastic x-ray scattering is consistent with spin-orbit coupling splitting the t 2 g manifold of octahedral Ir 4+ into a J eff = 1/2 and J eff = 3/2 state for all members of the series considered.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Majorana modes with side features in magnet-superconductor hybrid systems

Magnet-superconductor hybrid (MSH) systems represent promising platforms to host Majorana zero modes (MZMs), the elemental building blocks for fault-tolerant quantum computers. Theoretical description of such MSH structures is mostly based on simplified models, not accounting for the complexity of real materials. Here, based on density functional theory, we derive a superconducting 80-band model to study an MSH system consisting of a magnetic manganese chain on the s wave superconductor niobium. For a wide range of values of the superconducting order parameter, the system is a topological superconductor, with MZMs exhibiting non-universal spatial patterns and a drastic accumulation of spectral weight on both sides along the magnetic chain. These side feature states can be explained by an effective model which is guided by the ab initio results. Performing scanning tunneling spectroscopy experiments on the same system, we observe a spatial structure in the low-energy local density of states that is consistent with the theoretical findings. Our results open a first-principle approach to the discovery of topological superconductors.

36 MATERIALS SCIENCE↗

DOE OSTI2026

From Spin to Pseudospin Symmetry: The Origin of Magic Numbers in Nuclear Structure

Magic numbers lie at the heart of nuclear structure, reflecting enhanced stability in nuclei with closed shells. While the emergence of magic numbers beyond 20 is commonly attributed to strong spin-orbit coupling, the microscopic origin of the spin-orbit potential remains elusive, owing to its dependence on the resolution scale and renormalization scheme of nuclear forces. Here, we investigate the evolution of nuclear shell structure with varying momentum resolution in nuclear interactions derived from chiral effective field theory, using the similarity renormalization group to link different scales. We uncover a novel transition from spin symmetry to pseudospin symmetry as the resolution scale decreases, during which magic numbers emerge naturally. A similar pattern is found in calculations using relativistic one-boson-exchange potentials, underscoring the robustness of the phenomenon. This establishes a direct connection between realistic nuclear forces with a high resolution scale and effective nuclear forces at coarse-grained scales, offering a first-principles explanation for the origin of magic numbers and pseudospin symmetry in nuclear shell structure and new insights into the structure of exotic nuclei far from stability

Energy levels↗

DOE OSTI2023

First-Principles Understanding of Optical Excitations within Low-Dimensional Materials

The objective of this project is to utilize first-principles computational approaches to understand optical excitations within two-dimensional monolayers and heterostructures, with the ultimate goal of designing new materials by modifying the chemical and physical structure on the nanoscale. In 2D heterostructures, the nature and movement of the exciton is determined to a large extent by the electron-electron and electron-phonon strengths within individual monolayers, as well as the coupling between layers, all of which in turn depend on chemistry, solid-state screening properties, and quantum confinement. Controlling the nature and migration properties of excitons requires an understanding of the complex relationship between electrons, phonons, defects, disorder, and dynamics. The aim of this project is to determine, by analysis of highly accurate density functional theory (DFT) and many-body perturbation theory (MBPT) calculations, how these complexities can be decomposed into simple, tunable parameters. The specific goals of this project are to 1) understand electron conductivity, absorption and transparency of monolayer semiconductors and metals; 2) develop a theory of excitons near defects in low dimensions; and 3) understand the influence of inter-layer interactions in heterostructures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

Applications of Nickelate perovskites for neuromorphic computing from electronic structure and Machine Learning

While the limit of Moore's law is presently being reached with current microelectronic technologies, we need to develop new paradigms that overcome this limitation. In that respect, neuromorphic computing is a concept that emulates the neural behavior and response of the human brain, and it has been recognized as a promising alternative approach. In this research project, we will perform multi-fidelity scale bridging to explore the potential use of materials with metal to insulator transition for neuromorphic applications. In particular, rare earth nickelates are promising for such purposes, as the transition in these materials is quite sensitive to a broad set of different external stimuli. Our multi-fidelity approach will bridge the high-fidelity electronic structure calculations with classical potentials. We will bridge dynamical mean field theory with a classical atomistic representation via a deep learning force field. The neural network is trained with energies, charges, and forces obtained by accurate electronic structure theories based on Dynamical Mean Field Theory. The configurational space is generated from known crystal phases, ab initio molecular dynamics with exchange-correlation functionals corrected with the Hubbard model, disordered phases with different concentrations of oxygen vacancies, and nonsymmetrical positions and induced strain by grain interfaces or contact with a substrate. Strategies to train the model with a reduced number of training examples are obtained from active learning methods, and new structures for improving the learning process are generated by using machine learning autoencoders. This classical potential will be validated through a diversity of electronic structure methods and represents an important step to combine the flexibility and accuracy of first-principles with the speed of classical potentials. The generated multi-fidelity surrogate model will be used to understand the role of strain, oxygen vacancies, proton doping, the variation of the crystal phase, substrate effects, vibrational effects as the octahedral rotation, grain boundaries and defect effects on the response of a Metal to Insulator Transition (MIT) in correlated materials. Long time and large-scale simulations will help understand the role of different stimuli to control the hysteresis of the MIT, as it has been experimentally suggested. Selected configurations will be analyzed with higher-level theories to provide an accurate electronic description and to study how the orbitals and charges are rearranged under different conditions.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Hugoniot, sound speed, and phase transitions of single-crystal sapphire for pressures 0.2–2.1 TPa

Sapphire (Al 2 O 3 ) is a major constituent of the Earth's mantle and has significant contributions to the field of high-pressure physics. Constraining its Hugoniot over a wide pressure range and identifying the location of shock-driven phase transitions allows for development of a multiphase equation of state and enables its use as an impedance-matching standard in shock physics experiments. In this paper we present measurements of the principal Hugoniot and sound velocity from direct impact experiments using magnetically launched flyers on the Z machine at Sandia National Laboratories. The Hugoniot was constrained for pressures from 0.2–2.1 TPa and a four-segment piecewise linear shock-velocity–particle-velocity fit was determined. First-principles molecular dynamics simulations were conducted and agree well with the experimental Hugoniot. Sound-speed measurements identified the onset of melt between 450 and 530 GPa, and the Hugoniot fit refined the onset to 525 ± 13 GPa. A phase diagram which incorporates literature diamond-anvil cell data and melting measurements is presented.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2026

Transfer Learning Meets Embedded Correlated Wavefunction Theory for Chemically Accurate Molecular Simulations: Application to Calcium Carbonate Ion Pairing

Achieving chemical accuracy for molecular simulations remains a central challenge in computational chemistry. Here, we present an embedded correlated wavefunction transfer learning (ECW-TL) framework for accurately simulating molecular dynamics in the condensed phase. ECW-TL incorporates high-level electron exchange and correlation effects in ECW theory while preserving the training and computational efficiency of machine-learned interatomic potentials. We demonstrate the framework on Ca 2+ –CO 3 2– ion pairing in aqueous solution, a key process underlying CO 2 mineralization in seawater. As proof of principle, we first show that fine-tuning a DFT-revPBE-D3(BJ) baseline model with embedded-DFT-SCAN data reproduces the DFT-SCAN free-energy surface within 1 kcal/mol across all solvation states. Extending the framework to embedded MP2 and localized natural-orbital CCSD(T) further refines the free-energy profile, revealing the crucial role of exact electron exchange and correlation in determining ion-pair stability and structure. The computed ion-pair association free energy is in quantitative agreement with experimental measurements, further validating the accuracy of the ECW-TL framework. ECW-TL thus provides a general, data-efficient route for transferring CW accuracy to efficient simulations of complex aqueous and interfacial chemical processes.

cluster chemistry↗

DOE OSTI2025

Data-Driven Approach for Controlled Icosahedral Boron- Rich Compound Growth

This final technical report summarizes the research accomplishments and research highlights at the end of the funding period. This project aimed to leverage existing and new computational data produced from first-principles and molecular dynamics simulations to understand the thermodynamic, mechanical, and electronic properties of icosahedral boron compounds. The goal is to achieve targeted material properties by controlling the synthesis routes of these boron-rich compounds.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Theoretical framework for new magnetic materials for quantum computing and information storage. Final report for the Award No. DE-SC0018910

The focus of this grant was on molecular magnetic materials for information storage and quantum computing. We have been developing robust, first-principle methods for computing relevant electronic and magnetic properties of molecular building blocks (SMMs) of novel magnetic materials and quantum computers. These tools enable theoretical modeling of SMMs’ behavior, facilitating the interpretation of experimental studies and aiding the design of novel magnetic materials. Our strategy is based on the spin-flip (SF) approach, which extends the hierarchy of black-box single-reference methods to strongly correlated systems. Specifically, we developed general scalable algorithms and computer codes for calculating molecular properties, with an emphasis on spin-related properties, such as zero-field splittings, hyperfine couplings, and g-tensors. While our primary focus was on SF wave functions and SF-TDDFT, the underlying theory and computer codes were formulated using reduced density matrices, such that these tools are applicable to a broader class of methods. To extend the scope of applicability of wave-function-based SF methods to larger systems, we developed reduced-scaling approaches for the equation-of-motion coupled-cluster (EOM-CC) methods and continue developing libtensor (our open-source general tensor contraction library for many-body methods). We carried out extensive benchmarks and also carried out several applications.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Diffusion energy barrier of Au on Bi 2 Se 3 : theory and experiment

The stability and diffusion of ultra-thin thermally deposited Au films on Bi 2 Se 3 was studied using scanning tunneling microscopy and density functional theory calculations. The Au/Bi 2 Se 3 interface is of interest as gold is predicted to provide excellent electrical contact while maintaining the spin-polarized characteristics of the electronic states in Bi 2 Se 3 that make the material attractive for spintronic applications. When deposited at room temperature, Au 10 covers the surface with tightly packed islands of nanometer scale dimension. The surface morphology is stable up to 400K. At this annealing temperature, Au atoms have sufficient energy to diffuse across the surface and aggregate into larger nanostructures. At 550K, the Bi 2 Se 3 surface is only sparsely covered, and the Au has formed clusters with length scales 5-10 times larger than the original islands formed at room temperature. Comparison of the experiment and first principle calculation lead to the conclusion that the diffusion energy barrier for Au on Bi 2 Se 3 is as high as 0.47 eV, 15 which is much larger than diffusion barriers on other van der Waals materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

Hybrid physics-based and machine learning tools for materials assessment

In this work, we develop novel physics-based and machine learning computational techniques to predict fundamental properties of metallic systems that affect radiation damage behavior in reactor structural materials. Oftentimes, atomistic predictions of engineering alloys simplify chemical compositions to a single element to reduce computational cost and complexity, introducing large sources of uncertainty and potentially missing important behavior. In addition, engineering alloys such as SS 316 are particularly challenging to simulate with first principles methods because of the additional degrees of freedom introduced by magnetic natures of the constituent elements, and very little data of this type exists within the literature. These novel methods aim to improve the qualitative prediction of radiation damage in engineering alloys by more accurately simulating their compositions, both by accelerating the computations and by developing novel analyses.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗