Search NASA⌕ Search

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.

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

At least 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

First-principles investigation of the thermophysical properties of NaCl, PuCl 3 , and NaCl-PuCl 3 Molten salts

Molten salts have a variety of applications that span the nuclear and solar industries, and which involve thermal storage and heat transfer. There is a present knowledge gap in the thermophysical properties of molten salts, which limits the readiness level of molten salt applications. This is especially pertinent for molten salt reactors where the fissile material is dissolved within the molten salt. Ab initio molecular dynamics is a common method employed to investigate the structural and thermophysical properties at elevated temperatures. NaCl-PuCl 3 is a candidate fuel salt in molten salt reactors. The scope of this study is to investigate the seven unique compositions and a range of temperatures of NaCl-PuCl 3 and calculate the density, heat capacity, compressibility, enthalpy of mixing, and coefficient of thermal expansion. Within this work, the van der Waals (vdW) interactions are primarily handled with the vdW-DF2 functional but spot-checked with both the dDsC and DFT-D3 methods. The calculated densities are in good agreement with the NaCl and the eutectic compositions, and no experimental values exist for PuCl 3 at these temperatures. The densities are fit to a first-order Redlich-Kister expansion as a function of both composition and temperature. The heat capacity in the literature is scattered but the calculated values agree well with one of the experimental results. Here, the heat capacity increases as a linear function with respect to concentration at a rate of about 7.5 J/mol-K per 10 mol% PuCl 3 .

36 MATERIALS SCIENCE↗
DOE OSTI2024

First-principles investigation of lanthanides diffusion in HCP zirconium via vacancy-mediated transport

The diffusion of lanthanide fission products plays an important role in the growth of the fuel-cladding chemical interaction (FCCI) region in metallic fuels. The use of a Zr interdiffusion barrier may mitigate the transport of lanthanides from the fuel to the cladding, but the efficacy of such a liner is not yet known. In this paper, the stability and vacancy-mediated diffusion of La, Ce, Pr, and Nd in hexagonal close-packed (HCP) Zr is investigated via density functional theory (DFT) calculations and self-consistent mean field (SCMF) analysis. DFT is used to calculate the formation, binding, and migration energies of vacancies and vacancy-solute pairs. The DFT energetics are used in the KineCluE code to calculate the Onsager transport coefficients. La is found to be the fastest diffusing species in HCP Zr and experiences an almost isotropic diffusion behavior. The other three species (Ce, Pr, and Nd) demonstrate anisotropic diffusion where the diffusion in the basal planes is significantly faster than that along the c-axis. Further, the calculated lanthanide diffusivities in HCP Zr are fitted to an Arrhenius relation and the activation energies and prefactors are reported for the first time. Furthermore, the vacancy drag and the segregation tendencies were analyzed using the calculated off-diagonal transport coefficients. According to our vacancy-mediated diffusion model, lanthanides are expected to be enriched at vacancy sinks at low temperatures, while at high temperatures, lanthanides are depleted at sinks and will preferably diffuse into the bulk. The enrichment/depletion transition temperature depends on the diffusion direction (basal or axial) and hence will be controlled by the grain texture and orientation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

First-Principles Evaluation of Proton Hopping in Tetrahedral Oxide Motifs

Proton-conducting oxides (PCOs) are important materials used as ionic conductors for energy conversion technologies. Existing research efforts on PCO optimization and discovery generally focus on complex perovskite-based oxides that require doping and alloying to engineer oxygen deficiency and high proton conductivity. However, the variety of chemical compositions and coordination environments in oxides poses challenges for efficient materials design. In this computational study, we construct a database of simplified motifs to elucidate the relationship between fundamental materials chemistry and proton kinetics. Specifically, we focus on the zincblende crystal structure as a proxy for tetrahedral metal–oxide (M–O) coordination environments. We systematically quantified the effects of cation type, oxidation states, and M–O bond lengths on the proton hopping barrier, and found that strong M–O bonds and metal cations with large and variable oxidation states (e.g., Mo 6+ , V 5+ ) lead to smaller proton hopping barriers. By mapping the candidate cations and their preferred bond geometries onto materials databases such as the Inorganic Crystal Structure Database (ICSD) and Materials Project, we identified real materials containing the corresponding metal–oxide units. In general, we observed good agreement between the calculated proton hopping barriers obtained in real crystal structures and those predicted by our motif database. We also discuss the limitations of our model and possible future extensions to improve its predictive capabilities. Overall, our model provides a first step for the rational design and quick screening of energy-efficient PCOs.

organic↗
DOE OSTI2021

First principles calculations in support of Pu aging: calculating the effects of lattice imperfections on thermodynamics

Plutonium (Pu) has, in theory, well defined crystal structures: its atoms are arranged in regular spatial patterns. But Pu is radioactive, and as its nuclei decay those regular spatial patterns are interrupted. The interruptions are lattice imperfections, which are known to affect how materials respond to their environment. To adequately model Pu, we need to know which lattice imperfections are present, how they interact with each other, and how they affect the material’s response to its environment. The work presented here aims to use density functional theory (DFT) calculations to begin to answer the latter, in particular, how individual lattice imperfections affect measurable effects including thermal expansion (the change in volume in response to a change in temperature), heat capacity (the amount of thermal energy needed to change a material’s temperature), and elastic moduli (a material’s resistance to applied stresses). Pu poses many computational challenges. The $\textit{f}$ electrons require special attention. Of all the elements, Pu has the largest number of electrons that must be included in the calculations. Thermal effects require calculating the phonons (the lattice vibrations), which for systems containing lattice imperfections demand large, complex computational cells - but computational resources limit the size and complexity. With careful restructuring of how the calculations are performed, all these challenges have been met to enable calculations that provide insight into how lattice imperfections affect Pu’s response to its environment. Reported here are the computational challenges and the advances developed to meet them, along with the first results showing the strong effect that one prototype lattice imperfection (an interstitial Pu atom in a delta-phase Pu lattice) has on Pu’s response to its environment.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Determining the Adsorption Energetics of 2,3-Butanediol on RuO 2 (110): Coupling First-Principles Calculations With Global Optimizers

As climate change continues to pose a threat to the Earth due to the disrupted carbon cycles and fossil fuel resources remain finite, new sources of sustainable hydrocarbons must be explored. 2,3-butanediol is a potential source to produce butene because of its sustainability as a biomass-derived sugar. Butene is an attractive product because it can be used as a precursor to jet fuel, categorizing this work in the alcohol-to-jet pathway. While studies have explored the conversion of 2,3-butanediol to butene, little is understood about the fundamental reaction itself. We quantify the energetics for three pathways that were reported in the literature in the absence of a catalyst. One of these pathways forms a 1,3-butadiene intermediate, which is a highly exothermic process and thus is unlikely to occur since 2,3-butanediol likely gets thermodynamically trapped at this intermediate. We further determined the corresponding energetics of 2,3-butanediol adsorption on an ensemble of predetermined binding sites when it interacts with a defect-free stoichiometric RuO 2 (110) surface. Within this ensemble of adsorption sites, the most favorable site has 2,3-butanediol covering a Ru 5–coordinated cation. This approach is compared to that obtained using the global optimization algorithm as implemented in the Northwest Potential Energy Surface Search Engine. When using such a global optimization algorithm, we determined a more favorable ground-state structure that was missed during the manual adsorption site testing, with an adsorption energy of -2.61 eV as compared to -2.34 eV when using the ensemble-based approach. We hypothesize that the dehydration reaction requires a stronger chemical bond, which could necessitate the formation of oxygen vacancies. As such, this study has taken the first step toward the utilization of a global optimization algorithm for the rational design of Ru-based catalysts toward the formation of butene from sustainable resources.

09 BIOMASS FUELS↗
DOE OSTI2025

First-Principles Insights into the Thermocatalytic Cracking of Ammonia-Hydrogen Blends on Fe(110). 2. Kinetics

Ammonia (NH 3 ) is an energy-rich molecule that is routinely synthesized from nitrogen (N 2 ) and hydrogen (H 2 ). NH 3 ’s more favorable physical properties compared to H 2 suggests it may offer a way to more conveniently store, transport, and, when needed, extract H 2 via thermal decomposition. However, the high kinetic barrier and endoergicity to decompose to H 2 and N 2 require high temperatures. The standard reaction free energy indicates nearly 100% thermodynamic conversion to the diatomic molecules only at ~673 K and higher. However, even at these temperatures, a catalyst, e.g., iron (Fe), is needed for favorable kinetic conversion. Here, in this study, we explore via density functional theory the kinetics of NH 3 decomposition on the most stable facet of body-centered cubic Fe, namely, (110), under typical high-temperature and finite-pressure operando conditions. We predict coverage-dependent energetics of elementary surface reactions, often neglected in atomic-scale modeling. From these models, we find the recombinative desorption of adsorbed N as N 2 is rate-determining at 573.15–773.15 K and even at an extreme case of 1173.15 K. From microkinetic modeling, we find that the steady-state turnover frequencies (TOFs) for N 2 and H 2 generation rates (r$_{H_2}$) depend exponentially on temperature. The catalyst achieves a steady-state TOF of 36.4 s –1 and an r$_{H_2}$ of 0.107 μmol cm –2 s –1 for a feed of 1.8 bar NH 3 with 0.2 bar H 2 at 1173.15 K. However, at 773.15 K, with the same feed composition and velocity, the steady-state TOF and r$_{H_2}$ decrease to 0.14 s –1 and 4.10 × 10 –4 μmol cm –2 s –1 , respectively, as the process is significantly hindered by slow N 2 desorption. Although at first glance counterintuitive, our simulations suggest that surface modifications that reduce Fe’s reactivity toward NH x species should enhance its overall NH 3 decomposition activity.

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

Nuclear Theory from First Principles to Forefront Experiments (Final Report for DE-SC0018638)

The Lee research group is a part of the Nuclear Lattice Effective Field Theory Collaboration, which has developed and performed ab initio lattice simulations of nuclear structure, scattering, reactions, and many-body systems. The other senior members of the collaboration include Ulf-G. Meißner at Bonn/Julich, Evgeny Epelbaum and Hermann Krebs at Bochum, Timo Lahde and Thomas Luu at Julich, and Gautam Rupak at Mississippi State. Our letter “Ab initio alpha-alpha scattering” was featured in a Nature News and Views article. Another letter “Nuclear binding near a quantum phase transition” was highlighted in a Viewpoint article in the online APS journal Physics as well as a news article in the IOP publication Physics World (September 21, 2016). Our letter “Eigenvector continuation with subspace learning” was also highlighted a Synopsis article in Physics. There are many promising ab initio approaches being used to calculate the properties of few-and many-nucleon systems. This includes the no-core shell model, symmetry-adapted no-core shell model quantum Monte Carlo, auxiliary-field diffusion Monte Carlo, self-consistent Green’s functions, many-body perturbation theory, in-medium similarity renormalization group, and coupled cluster methods. Lattice effective field theory is another ab initio approach which combines the framework of effective field theory with lattice Monte Carlo methods to allow favorable scaling from few- to many-body systems. Perhaps the most important aspect of lattice effective field theory is that its strengths and weaknesses are orthogonal to that of other ab initio methods. For example, lattice effective field theory has little difficulty in probing cluster structures and collectivity, while such features are much more difficult using other methods. Furthermore it can be used to compute superfluid condensate fractions as well as the phase diagram of strongly matter and the density and temperature dependence of clustering. Lattice effective field theory was first used in simulations of infinite nuclear matter and infinite neutron matter with pions and without pions. In addition to the efforts by our collaboration, there have been recent efforts by other groups as well.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

First principles investigation of dopants and defect complexes in CdSe$_x$Te$_{1-x}$

Se alloying is a common approach to improve the performance of CdTe solar cells by tuning the bandgap, defect levels, and carrier density. A fundamental understanding of these improvements, specifically the effect of Se alloying on the behavior of defects and dopants in CdTe, remains unclear. Here, in this work, we present a density functional theory (DFT) study of point defect energetics in CdTe and CdSe x Te 1-x with x = 0.25, leading to a comparison of how native defects, dopants (As and Cu), impurities (Cl and O), and related defect complexes behave in CdTe vs CdSe x Te 1-x . Our calculations, performed by combining semi-local and nonlocal hybrid functionals, show a general lowering of the formation energies of native defects as well as substitutional defects formed by As and Cl upon Se addition. For successful p-type doping with As, destabilizing Cl-based defects in the CdSeTe lattice would be essential. We find evidence for some low-energy defect complexes of As, Cl, and O in CdSe 0.25 Te 0.75 . The computed defect formation energies further enable estimates of temperature-dependent defect concentrations and self-consistent Fermi levels. A comparison of defect energetics with the energies of impurity phases reveals that As, Cu, Cl, and O overwhelmingly prefer being segregated to unwanted As 2 O 5 , AsCl 3 , Cd 2 AsCl 2 , and CuO x phases rather than remain at defect sites, but such segregation is less likely to happen in CdSe 0.25 Te 0.75 than in CdTe. Overall, our work presents a list of likely defects and complexes in CdTe and Se-incorporated CdTe, paving the way to explain and mitigate limited dopant activation in experimental observations.

CdTe↗
DOE OSTI2024

First Principles Study of Aluminum Doped Polycrystalline Silicon as a Potential Anode Candidate in Li‐ion Batteries

Addressing sustainable energy storage remains crucial for transitioning to renewable sources. While Li‐ion batteries have made significant contributions, enhancing their capacity through alternative materials remains a key challenge. Micro‐sized silicon is a promising anode material due to its tenfold higher theoretical capacity compared to conventional graphite. However, its substantial volumetric expansion during cycling impedes practical application due to mechanical failure and rapid capacity fading. A novel approach is proposed to mitigate this issue by incorporating trace amounts of aluminum into the micro‐sized silicon electrode using ball milling. Density functional theory (DFT) is employed to establish a theoretical framework elucidating how grain boundary sliding, a key mechanism involved in preventing mechanical failure is facilitated by the presence of trace aluminum at grain boundaries. This, in turn, reduces stress accumulation within the material, reducing the likelihood of failure. To validate the theoretical predictions, capacity retention experiments are conducted on undoped and Al‐doped micro‐sized silicon samples. In conclusion, the results demonstrate significantly reduced capacity fading in the doped sample, corroborating the theoretical framework and showcasing the potential of aluminum doping for improved Li‐ion battery performance.

25 ENERGY STORAGE↗
DOE OSTI2020

Understanding initial zeolite oligomerization steps with first principles calculations

Zeolites are porous aluminosilicate materials that find various applications in the chemical industry in separations, catalysis, ion exchange, and so forth. However, despite their widespread use, the reaction mechanisms occurring during zeolite growth are still unclear. In the study, we use density functional theory calculations to gain insights into the thermodynamics of oligomerization, which constitute the initial steps of zeolite growth. By taking into consideration solvent and temperature effects, our results demonstrate that the growth of aluminosilicate systems is significantly more exothermic than their pure silicate counterparts. Under pH neutral conditions, water prefers to dissociate on the early-growth-stage aluminosilicate complexes rather than desorb, thus generating potential Brønsted acid sites on the oligomers. Additionally, (alumino)silicate growth pathways are evaluated in the presence of Na + cation, as well as the Ca 2+ cation for the pure silicate pathway. The presence of cations increases the exothermicity of growth, with Ca 2+ exhibiting the most energetically favorable growth environment for the silicate systems. Importantly, we demonstrate through reaction extent analysis that the presence of cations modulates the speciation of the formed oligomers, with Na + favoring linear species in addition to the generally preferred cyclic ones. Overall, this work provides a fundamental understanding of the thermodynamics of complex reaction paths that occur during early stages of zeolite growth and suggests that the initial growth steps can have significant impact on the final zeolite structure.

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

Rhodium Single–Atom Catalysts on Titania for Reverse Water Gas Shift Reaction Explored by First Principles Mechanistic Analysis and Compared to Nanoclusters

The thermocatalytic reduction of CO2 by H2 often proceeds via two competing reaction mechanisms—the reverse water gas shift reaction (rWGSR, CO2 + H2 ⇌ CO + H2O) and methanation (CO2 + 4H2 ⇌ CH4 + 2H2O). Atomically dispersed Rh1 catalysts on TiO2 show high selectivity toward the rWGSR compared with larger Rh nanoclusters, but the origin of this size-dependent selectivity has not been fully explained. Here we report density functional theory (DFT) calculations and microkinetic simulations that clarify the Rh1 active sites and rWGSR pathway on anatase TiO2(101), as well as the high rWGSR selectivity of Rh1 compared with supported Rhx (x = 2–8 atoms) nanoclusters. DFT-computed formation energies, vibrational frequency analysis, and microkinetic modeling suggest three plausible active sites: Rh1 on titania (Rh1/TiO2(101)), Rh1 with a nearby hydroxyl group (Rh1OH/TiO2(101)), and Rh1 near an oxygen vacancy at a three-fold coordinated site (Rh1 near O3cvac). Predicted turnover frequencies and apparent activation barriers for Rh1 indicate a faster reaction involving CO2 dissociation assisted by a support oxygen vacancy via Rh1 near O3cvac, as well as slower reactions involving Rh1OH/TiO2(101) or Rh1/TiO2(101) through a COOH intermediate. Furthermore, these Rh1 sites are selective toward CO rather than CH4 because of the weak adsorption of CO, large barrier for C-O bond dissociation, and the lack of nearby metal sites for H2 dissociation, in contrast to Rhx nanoclusters, including Rh2 dimers.

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

The α-WC(0001) Surface as a Hydrogen Sponge: A First Principle Study of H 2 Dissociation and Formation of Low and High Coverages

Tungsten carbide (WC) displays a Pt-like behavior in catalysis, applied in hydrogenation processes. Numerous theoretical studies have modeled the behavior and use of adsorbed hydrogen without obtaining a general picture, missing basic links between H 2 dissociation and generation of high surface coverage (Θ H >0.5 ML). Here, in this study, the capability of C- and W-terminations of the α-WC(0001) surface is analyzed to dissociate several H 2 molecules to produce coverages, Θ H , ranging from low to very high values (0.13<Θ H <2.00 ML). Density functional theory and an ab initio atomistic thermodynamic were used to achieve the conditions for H 2 dissociation. The WC–C surface has higher capacity to dissociate H 2 molecules than WC–W. However, both surfaces can reach full surface coverage, Θ H =1 ML, at mild ambient conditions, T=300 K and P=1 atm, and even up to 500 K at low and high pressures. The H-adatoms on WC–W are more labile than on WC–C. The binding of adsorbates is hindered at high Θ H , implying a need to modulate Θ H according to the application. The results give the basis to understand the capabilities of WC-based catalysts in hydrogenation-related reactions, with the advantage of WC being a hydrogen reservoir at mild practical catalytic conditions.

03 NATURAL GAS↗
DOE OSTI2021

Algorithm advances and applications of time‐dependent first‐principles simulations for ultrafast dynamics

Abstract Far from equilibrium phenomenon is a central theme of contemporary material research. Such phenomenon can exhibit itself in atomic structure and dynamics, but very often it also happens as non‐equilibrium phenomenon in the electronic structure. In ab initio material simulation, density functional theory (DFT) has played an essential role in studying electronic ground state problems. For excited states, besides many‐body perturbation theory, another powerful tool is the time dependent DFT (TDDFT) method. In particular, the real‐time TDDFT (rt‐TDDFT) method can be used to simulate many non‐equilibrium phenomena directly. Here we introduce our works on some algorithm advances based on our recently rt‐TDDFT method. This method uses the plane‐wave basis set, and significantly accelerates its efficiency by increasing the time step from 0.1–1 as in traditional methods to 0.2–0.5 fs. The noncollinear magnetic moments and spin–orbit coupling have also been included in our rt‐TDDFT method. Furthermore, a Boltzmann‐TDDFT algorithm has been developed to solve the hot carrier overheating problem in Ehrenfest dynamics, and a natural orbital branching algorithm has been developed to overcome the mean‐field approximation in Ehrenfest dynamics nuclear trajectory, thus allows stochastic multiple paths in chemical reactions. Utilizing these methods, we have studied the photoinduced ultrafast demagnetization, ultrafast phase transition, energy transfer between plasmon and hot carriers, as well as the high‐energy ion implantation and low‐energy atomic diffusion in semiconductors. We believe the tools as the ones introduced here can enable us to study a wide range of phenomena which are of great interest in modern day material research. This article is categorized under: Structure and Mechanism > Computational Materials Science Electronic Structure Theory > Ab Initio Electronic Structure Methods Electronic Structure Theory > Density Functional Theory

Liu, Wen‐Hao↗
DOE OSTI2021

Black hole induced false vacuum decay from first principles

We provide a method to calculate the rate of false vacuum decay induced by a black hole. The method uses complex tunneling solutions and consistently takes into account the structure of different quantum vacua in the black hole metric via boundary conditions. The latter are connected to the asymptotic behavior of the time-ordered Green’s function in the corresponding vacua. We illustrate the technique on a two-dimensional toy model of a scalar field with inverted Liouville potential in an external background of a dilaton black hole. We analytically derive the exponential suppression of tunneling from the Boulware, Hartle-Hawking and Unruh vacua and show that they are parametrically different. The Unruh vacuum decay rate is exponentially smaller than the decay rate of the Hartle-Hawking state, though both rates become unsuppressed at high enough black hole temperature. We interpret the vanishing suppression of the Unruh vacuum decay at high temperature as an artifact of the two-dimensional model and discuss why this result can be modified in the realistic case of black holes in four dimensions.

2D Gravity↗