SEARCH · Search NASA

Results for “First Principles”

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 595 records · Page 33

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Cluster-Graph Fingerprinting: A Framework for Quantitative Analysis of Machine-Learned Interatomic Model Training and Simulation Data

Machine-learned interatomic models represent a significant advancement in simulation methods, extending the predictive ability of first-principles methods to previously inaccessible length and time scales. However, the data-driven nature of these models can lead to difficult-to-detect errors that can compromise prediction accuracy. To address this challenge, we introduce a novel fingerprinting approach based on the Chebyshev Interaction Model for Efficient Simulation (ChIMES) ML-IAM graph-based descriptor. Our strategy enables efficient and statistically rigorous analysis of system configurations used in ML-IAM training and those generated by their application, e.g., in molecular dynamics simulations. We demonstrate that these fingerprints can effectively assess novelty of a configuration relative to an existing data set and determine dissimilarity among individual configurations, which are two key tasks in workflows for active learning-based ML-IAM training, data set curation, and on-the-fly uncertainty quantification.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Perspective on Many-Body Methods for Molecular Polaritonic Systems

Recent advances in strong light–matter interactions have revealed a wealth of new physical phenomena in molecules embedded in optical cavities, including modified chemical reactivity, altered excitation spectra, and novel quantum correlations. To describe these effects from first-principles, the field of ab initio quantum electrodynamics (QED) has emerged as a compelling extension of quantum chemistry that treats electronic and photonic degrees of freedom on equal footing. In this Perspective, we review the growing landscape of many-body QED methods, including Hartree–Fock, density functional theory (QEDFT), time-dependent DFT (QED-TDDFT), configuration interaction (QED-CI), complete active space (QED-CASSCF), coupled cluster (QED-CC), quantum Monte Carlo (QED-QMC), and density matrix renormalization group (QED-DMRG), highlighting recent developments and implementations. We further explore real-time methods, gradient and Hessian formalisms, and the integration of nonadiabatic nuclear dynamics. Applications range from benchmark simulations of polaritonic chemistry to quantum simulations on emerging quantum hardware. We conclude by outlining future directions for theory development and interdisciplinary efforts at the interface of quantum chemistry, condensed matter, and quantum optics.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Characterizing Defect Dynamics in Silicon Carbide Using Symmetry-Adapted Collective Variables and Machine Learning Interatomic Potentials

Silicon carbide (SiC) divacancies are attractive candidates for spin-defect qubits possessing long coherence times and optical addressability. The high activation barriers associated with SiC defect formation and motion pose challenges for their study by first-principles molecular dynamics. In this work, we develop and deploy machine learning interatomic potentials (MLIPs) to accelerate defect dynamics simulations while retaining ab initio accuracy. We employ an active learning strategy comprising symmetry-adapted collective variable discovery and enhanced sampling to compile configurationally diverse training data, calculation of energies and forces using density functional theory (DFT), and training of an E(3)-equivariant MLIP based on the Allegro model. Here, the trained MLIP reproduces DFT-level accuracy in defect transition activation free energy barriers, enables the efficient and stable simulation of multidefect 216-atom supercells, and permits an analysis of the temperature dependence of defect thermodynamic stability and formation/annihilation kinetics to propose an optimal annealing temperature to maximally stabilize VV divacancies.

Computer simulations↗

DOE OSTI2026

Charge Transport in Solvated Donor–Acceptor Functionalized Peptoids: Molecular Dynamics and Rate Theory

Scalable solar-energy conversion requires photoactive materials that combine the efficiency of natural photosynthetic systems with the stability and processability needed for practical applications. Achieving reliable charge transport in soft, self-assembled organic materials remains challenging, as structural fluctuations and environmental effects strongly influence charge-transfer (CT) rates. Here, we present a broadly applicable computational framework for evaluating CT rates in the condensed phase, combining Fermi’s golden rule rate theory with inputs from all-atom molecular dynamics (MD) simulations and first-principles electronic-structure calculations. The approach does not rely on system-specific parametrization and is applicable to a wide range of soft and disordered materials. We demonstrate the applicability and usefulness of the framework on redox-active peptoids functionalized with iron–porphyrin (Fe–P) complexes, a bioinspired platform with programmable donor–acceptor units and tunable three-dimensional organization. The calculated CT rates exhibit strong sensitivity to molecular conformation, with variations spanning several orders of magnitude. This dependence is shown to arise from the pronounced variation in diabatic electronic coupling with the relative orientations and separations of the Fe–P complexes across the conformational ensemble. The framework provides a consistent route for connecting atomistic structure to CT kinetics in the condensed phase and enables analysis of structure–rate relations in organic semiconducting systems.

Charge transfer↗

DOE OSTI2025

Structure of Novel Phosphonium-Based Ionic Liquids with S and O Substitutions from Experiments and a Mixed Quantum-Classical Approach

This article presents experimental characterization information and synchrotron X-ray scattering measurements on a set of novel O- and S-substituted phosphonium-based ionic liquids (ILs) all coupled with the bis(fluorosulfonyl)imide (FSI-) anion. The ILs include the ethoxyethyltriethylphosphonium (P 222(2O2) + ) and triethyl[2-(ethylthio)ethyl]phosphonium (P 222(2S2) + ) cations, and we contrast results on these with those for unsubstituted triethylpentylphosphonium (P 2225 + ). The article also introduces a physics-based protocol that combines classical force field studies on larger simulation boxes with classical and first-principles studies on smaller boxes. The method produces significantly improved S(q) functions in the regime which in prior publications we have associated with inter- and intraionic adjacency correlations. By understanding which shorter-range structural changes improve S(q) in the q-regime of interest, we are also able to pinpoint specific deficiencies in the classical force field model. The approach we take should be quite general and could help study other complex liquids on different length scales.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Ionic Liquids under Radiation and the Dimer Radical Dicyanamide Anion

Ionic liquids (ILs) used in solar, electrochemical, nuclear, or even space exploration applications will necessarily be exposed to conditions in which electron detachment may occur and transient charge-depleted radical species may form. In a recent article (J. Am. Chem. Soc. 2025, 147, 23395–23398), we studied the different possible kinetic fates of an excess electron in the popular 1-butyl-1-methylpyrrolidinium dicyanamide IL but did not address the nature of electron-deficient transient radical species these leaving electrons produce, the so-called holes. Using pulse radiolysis measurements and first-principles calculations in the bulk and gas phases to interpret these, we show that the most likely hole species in this IL is a dimer radical anion that, somewhat surprisingly, has similar spectral features to the excess electron, including a broad near-infrared absorption band.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

THD-C Sheet: A Novel Nonbenzenoid Carbon Allotrope with Tetra-, Hexa-, and Dodeca-Membered Rings

Here, we propose a novel two-dimensional carbon-based structure with tetra-, hexa-, and dodeca-membered rings, which we refer to by the abbreviated name, THD-C. The structure presents a mixture of sp–sp 2 hybridization and can potentially be synthesized by the topological assembly of 4-ethynyldiphenylacetylene molecules. By employing first-principles calculations, the stability and ease of synthesis of the sheet are investigated and compared with various C-allotropes. We predict its metallic behavior and excellent kinetic and dynamic stability. Due to the crystal structure of the sheet, a strong mechanical anisotropy is observed. The effects of functionalization on the electronic properties of the material are also studied, and different semiconducting systems are obtained. The potential of THD-C for energy storage in metal-based batteries, hydrogen storage, and catalysis is also investigated, and we find a superior performance in comparison to graphite and other allotropes. The quantum confinement effect is investigated by constructing nanoribbons and nanotubes of various sizes. For ribbons, we find that tailor-made electronic and magnetic properties can be obtained and explored in potential spintronic devices. Additionally, we observe that nanotubes are conducting irrespective of their chirality and can potentially be used for capture, storage, and separation of industrially relevant small gas molecules.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Defect Properties, Anion Ordering, and Photochromic Mechanism in Yttrium Oxyhydride

Yttrium oxyhydride (YHO) undergoes a reversible photochromic transition when exposed to ultraviolet light. However, the mechanism for this transformation is not fully understood, and the structure and precise chemical composition of YHO remain under debate. Here, we use first-principles density functional theory calculations with a hybrid functional to study the structure, chemical stability, and point defect properties of YHO. As experiments have shown, we find that YHO prefers a cubic structure, with H and O anions present in equal concentrations and located on tetrahedral sites. Stoichiometric and ordered YHO is chemically stable, but it has a wide band gap of 5.01 eV, considerably larger than that measured in experiments (2.4–3.8 eV). On the other hand, Y4H10O has a smaller band gap of 2.97 eV and also has a region of chemical stability; thus, the actual material may include some fraction of this H-rich structure. The defect chemistry of YHO is dominated by anionic antisite species (H O and O H ), with hydrogen interstitials (H i ) and vacancies (V H ) also present in reasonably high concentrations. We show that antisite disorder lowers the band gap relative to the perfectly ordered structure, bringing the magnitude of the gap into closer agreement with experiment. Based on our calculations of defect migration and the positions of defect states relative to the band edges, we link the onset of photochromic behavior to the reaction H O – → V O 0 + H i – , which follows photoexcitation of a H O + defect. H i – can subsequently migrate away and be trapped by additional H O + defects, contributing to the persistence of the reaction, while the resultant oxygen vacancy, V O 0 , introduces an occupied defect state that leads to optical absorption at visible wavelengths. Our results can explain reported discrepancies between experimental and computational results for YHO, and they allow us to propose specific atomic-scale processes that can lead to photochromism. In conclusion, understanding these mechanisms is key for unlocking YHO’s application in devices ranging from smart windows and optoelectronics to electrochemical synapses for neural networks.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Large Electrically and Chemically Tunable Rashba–Dresselhaus Effects in Ferroelectric CsGeX 3 (X = Cl, Br, I) Perovskites

Rashba–Dresselhaus effects, which originate from spin–orbit coupling and allow for spin manipulations, are actively explored in materials, following the pursuit of spintronics and quantum computing. However, materials that possess practically significant Rashba–Dresselhaus effects often contain toxic elements and offer little opportunity for the tunability of the effects. We used first-principles simulations to reveal that the recently discovered halide ferroelectrics in the CsGeX 3 (X = Cl, Br, I) family possess large and tunable Rashba-Dresselhaus effects. In particular, they give origin to the spin splitting of up to 171 meV in the valence band of CsGeI 3 . The value is chemically tunable and can decrease by 25% and 70% for CsGeBr 3 and CsGeCl 3 , respectively. Such chemical tunability could result in the engineering of desired values through a solid solution technique. Application of an electric field was found to result in structural changes that could decrease and increase spin splitting, leading to electrical tunability of the effect. In the vicinity of conduction and valence band extrema, the spin textures are mostly of the Rashba type, which is promising for spin-to-charge conversion applications. The spin directions are coupled with the polarization direction, leading to Rashba-ferroelectricity cofunctionality. Furthermore, our work identifies lead-free perovskite halides as excellent candidates for spin-based applications and is likely to stimulate further research in this direction.

Electric fields↗

DOE OSTI2025

Efficient Modeling of Structural, Electronic, and Optical Properties of Silver and Gold Metal Nanoclusters and Alloys Using Optimized SCC-DFTB Parameters

Computation of optical properties using conventional time-dependent density functional theory (TD-DFT) is time-consuming and memory-intensive. In this study, we investigate the accuracy and efficiency of the density functional tight binding (DFTB) framework with newly optimized Slater–Koster (SK) parameters for modeling the structural, electronic properties, and absorption spectra of silver and gold nanoclusters and their alloys. Our investigation of the ground state (GS) properties demonstrates that the newly developed GS-SK parameters enable DFTB to closely approximate DFT-calculated bond lengths for octahedron, tetrahedron, icosahedra, and truncated octahedron with sizes Ag n /Au n (n = 19, 20, 38, 55), nanoclusters and Ag 20 /Au 20 nanoalloys, with a maximum deviation of approximately 0.15 Å. Formation energy results indicate that the GS-SK parameters can closely estimate changes in formation energies with alloy composition, and the comparison of electronic structures for Ag 20 , Au 20 , and AgAu alloy nanoclusters using the DFTB approximation reveals good agreement in the projected density of states (DOS) profiles and energy levels. A second set of SK parameters, ES-SK, has been developed to describe excited state (ES) properties, including the absorption spectra of silver octahedron Ag 19 , tetrahedral Ag n (n = 20, 56, 84), truncated octahedron Ag 38 , and icosahedra Ag 55 closed-shell clusters and their gold and alloy counterparts over a broad range of alloy compositions. This parametrization uses TD-DFTB calculations and fine-tunes the d and p eigenvalues by comparing them to reference absorption spectra from first-principles TD-DFT. This enables the generation of absorption spectra that closely match the reference spectra when plasmon excitation is dominant, as demonstrated by studying the plasmonic properties of icosahedral Ag n and Au n (n = 309 and 561) nanoparticles. This includes the rapid loss in plasmon quality when Au partially replaces Ag in alloy clusters. Furthermore, these results provide a foundation for addressing computational bottlenecks in plasmonics and with new prospects for applications in the quantum plasmonics for bimetallic alloys.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Morphology, Deformations, and Photocatalytic Activity of Thermally Treated Brookite Titanium Dioxide Thin Films

Metastable states, in which the coupling between long-range lattice deformations and electronic properties can be controlled, provide a pathway to tailoring the behavior of photocatalytic materials by directing the flow of photoinduced charge carriers. Brookite is a metastable polymorph of earth-abundant TiO 2 that exhibits photocatalytic function and, due to its high energy relative to the anatase and rutile polymorphs, may serve as a precursor for the formation of transitional metastable structures. In this work, facile thermal annealing is employed to promote the formation of predominantly brookite-phase films, regulate the brookite lattice distortions, and determine the effect of these distortions on charge separation, ultimately directed at enhancing photocatalytic activity. Profile fitting of X-ray diffraction patterns and peak shifts in Raman spectra revealed structural distortions of the brookite lattice. Structural defects, including lattice gliding, dislocations, stacking faults, and twin boundaries, were observed using scanning transmission electron microscopy. First-principles simulations reveal how the lattice distortions associated with stacking faults induce band bending, thus increasing the photocatalytic activity of brookite. In conclusion, this study provides insight into the microstructural tuning of metastable phases to enhance their unique functionalities.

band bending↗

DOE OSTI2025

Computational Study of Indium Oxide Photoelectrodes

Using a combination of first-principles molecular dynamics simulations and electronic structure calculations, we characterize the atomistic structure and vibrational properties of a photocatalytic surface of In 2 O 3 , a promising photoelectrode for the production of hydrogen peroxide. We then investigate the surface in contact with water and show that the electronic states of In 2 O 3 are appropriately positioned in energy to facilitate the two-electron water oxidation reaction (WOR) over the competing four-electron oxygen evolution reaction. In conclusion, we further propose that the use of strained thin films interfaced with water is beneficial in decreasing the optical gap of In 2 O 3 and thus in utilizing a wider portion of the solar spectrum for the WOR.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Theoretical Investigations of Anharmonic Effects and Phonon Transport in the Cubic Phase of Crystalline Perovskite CsPbCl 3

The role of anharmonic effects on lattice dynamics and thermal transport has been investigated in the cubic phase of the CsPbCl 3 perovskite. Limitations of the harmonic approximation, which lead to phonon instabilities in the Brillouin zone, were addressed based on self-consistent phonon theory in combination with first-principles calculations and by incorporation of frequency renormalization effects from bubble and loop diagrams. This theoretical approach demonstrates significant improvements compared to results obtained using a self-consistent phonon dynamic matrix. Both the cubic-to-tetragonal phase transition temperature as well as the predicted lattice thermal conductivity values show good agreement with selected sets of available experimental data. In conclusion, the accuracy of the predicted thermal conductivity data is also tested against systems significantly larger than those allowed by quantum calculations, by using molecular dynamic simulations with machine-learning force fields.

Lattices↗

DOE OSTI2025

An Atomistic Study of Reactivity in Solid-State Electrolyte Interphase Formation for Li/Li7P3S11

Lithium metal batteries offer superior volumetric and gravimetric specific capacities compared to those based on traditional graphite anodes. Although advancements in solid-state electrolytes address safety concerns, challenges remain, particularly regarding interphase formation in lithium metal anodes. This work presents a computational framework based on high-throughput first-principles density functional theory and machine-learning interatomic potentials (MLIPs) including automated iterative, active learning to enable robust computational exploration of interphase formation between lithium metal anodes and an inorganic solid-state electrolyte. As a demonstration, we apply the framework to a Li/Li7P3S11 interface and find that it accurately identifies the experimentally observed, thermodynamically stable interphase products as well as their overall spatial arrangement within a heterogeneous, amorphous layered structure, with Li2S domains of nanocrystallinity. Our simulations show two stages, a fast and slow diffusion reaction regime, that corroborate the relative phase formation rate of Li x P, Li2S, and Li3P. Using the Onsager transport theory, we capture time-dependent ionic diffusion within the reacting interface, including cross-correlation effects. We found that cross-correlation effects between Li-P and P-S ionic motion significantly influence P-ion diffusion, making it highly sensitive to the local environment and potentially leading to "kinetic trapping" of Li-P phases. The passivation of the interface is shown as the ionic fluxes all approach zero, effectively halting interphase growth.

Diffusion↗

DOE OSTI2025

Stabilization of Preternatural Barium Oxidation States as an Unexpected Byproduct of β-Decay: Discovery of a New Halide Semiconductor Alloy

137 Cs has a wide range of roles in the nuclear industry. The solid material, safely encapsulated in CsCl as 137 CsCl, is stored as fission product waste from nuclear power production and legacy waste from nuclear weapons production; it has also served as a radiation source in food and sewage irradiators as well as medical devices. However, because of the solubility of the chloride salt and the relatively high specific activity of 137 Cs, damaged or broken capsules can lead to severe radiological accidents. Safe capsule design and material recycling are complicated by the unclear structural evolution during β-decay, which remains ambiguous due to the differing oxidation states of Cs (1+) and Ba (2+). Here, in this study, we use first-principles calculations to investigate the evolving structure–property relationships of Cs 1–x Ba x Cl during β-decay. Despite the well-established 2+ formal oxidation state of alkali-earth metals, we find that Ba 1+ can be stabilized in the form of a mixed-valence alloy at low concentrations. Specifically, we identify three regimes for the β-decay of 137 Cs into CsCl: Ba-doped CsCl (Ba ≤ 14%), wherein Ba has the expected 2+ oxidation state; Cs–Ba–Cl alloys, where Ba has a mix of the usual Ba 2+ and highly unusual Ba 1+ oxidation state in the form of a quasi-disordered mixed-valence alloy (Ba = 25%); and phase separation into a CsCl + BaCl 2 + Ba (m) mechanical mixture, where Ba reverts to its expected 2+ oxidation state (Ba > 25%). Surprisingly, the Cs 0.75 Ba 0.25 Cl mixed-valence alloy is a narrow indirect band gap semiconductor (1.05 eV) despite the insulating nature of both CsCl and BaCl 2 . It also exhibits strongly excitonic polarized optical properties, has glass-like ultralow thermal conductivity (directional average of 0.21 W/mK at 300 K), and shows greater resistance to deformation under both tensile and volumetric strengths compared with the original CsCl structure (e.g., shear and Young’s modulus of 9.04 and 31.62 GPa, respectively). These findings imply that transmutation of 137 Cs leads to highly unusual chemical bonding that stabilizes Ba 1+ in local regions of the quasi-disordered Cs 0.25 Ba 0.75 Cl, resulting in anomalous physical properties. Moreover, this discovery provides valuable insight for safe nuclear waste capsule design, which can aid in preventing environmental or human exposure to radioactive materials.

Fuhr, Addis S. [Oak Ridge National Laboratory (ORN↗

DOE OSTI2026

Toward Chemical Accuracy for Chemi- and Physisorption with an Efficient Density Functional

Understanding molecular adsorption on surfaces underpins many problems in chemistry and materials science. Accurately and efficiently describing the adsorption has been a challenging task for first-principles methods as the process can involve both short-range chemical bond formations and long-range physical interactions, e.g., van der Waals (vdW) interaction. Density functional theory presents an appealing choice for modeling adsorption reactions, although calculations with many exchange-correlation density functional approximations struggle to accurately describe both chemical and physical molecular adsorptions. Here, we propose an efficient density functional approximation that is accurate for both chemical and physical adsorption by concurrently optimizing its semilocal component and the long-range vdW correction against the prototypical adsorption CO/Pt(111) and Ar 2 binding energy curve. The resulting function opens the door to accurate and efficient modeling of general molecular adsorption.

adsorption↗

DOE OSTI2026

Machine Learning Interatomic Potentials for Modeling Framework Flexibility and Water Uptake in NbOFFIVE-1-Ni Metal–Organic Framework

Metal–organic frameworks (MOFs), with their distinctive porous structures and tunable chemical properties, have shown immense promise in the separation and storage of gases. Currently, the accurate simulation of their adsorptive properties remains challenging, especially for systems where the molecules fit very tightly into the pores. Traditional simulation methods often approximate the frameworks as rigid and do not account for the framework flexibility seen in materials such as NbOFFIVE-1-Ni. First-principles molecular dynamics (FPMD) simulations offer the desired accuracy in modeling this flexibility but are limited by their extensive computational demands, rendering them impractical for long simulations. Conversely, classical force field-based simulations offer computational efficiency but lack the necessary accuracy. Here, to break this accuracy-efficiency trade-off, we have developed machine learning interatomic potentials trained on energies and forces from FPMD to model the framework flexibility of NbOFFIVE-1-Ni in the presence of water over nanosecond time scales. Furthermore, by integrating MLIP-driven molecular dynamics (MLIP-MD) with grand canonical Monte Carlo (GCMC) simulations, we further incorporated framework flexibility into adsorption predictions, yielding water adsorption isotherms that better align with experimental data compared to those of conventional GCMC simulations. These advances offer new opportunities for the design and optimization of MOFs in gas storage and separation applications.

adsorption↗

DOE OSTI2024

The Gas Phase Protonation Sites of Six Naturally Occurring Nicotinoids

We have determined the gas phase protonation sites of six naturally occurring nicotinoids, namely nicotine (NIC), nornicotine (NOR), anabasine (ANB), anatabine (ANT), cotinine (COT) and myosmine (MYO) from experimental cryogenic ion trap infrared spectroscopy and first principles electronic structure calculations. The structural motif of these compounds consists of a common Pyridine and a differing non-Pyridine rings. The protonation site on either of these two rings correlates with the nicotinoid’s biological activity (addiction). At room temperature NIC is a mixture of Pyridine and Pyrrolidine (non-Pyridene) protomers, NOR, ANB, ANT and COT are pure Pyridine protomers and finally MYO is mostly a Pyroline (non-Pyridine) protomer. Their biological activity (addiction) correlates with the protonation preference in the nonPyridine protomer, with the exception of MYO, which, although is protonated in the Pyroline (non-Pyridine) protomer, it is not addictive because of a conjugated p system that induces a flat structure that is not optimally interacting with the environment in the in the binding pocket of the human brain.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗