Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Density functional”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

First principles density functional theory study of tritium species adsorption on Ni(111) surface and diffusion in nickel-sublayer for tritium storage

The nickel-plated zircaloy-4 is used as a tritium ( 3 H) getter in the tritium-producing burnable absorber rods (TPBARs) to capture 3 H produced in the 6 Li-riched annular γ-LiAlO 2 pellet under neutron irradiation. The experimental data and our previous theoretical results showed that the 3 H species produced from the γ-LiAlO 2 pellet were mainly 3 H 2 and 3 H 2 O. These 3 H species diffuse from the surface of the LiAlO 2 pellet across vacuum to the nickel-plated zircaloy-4 getter and then further diffuse into the getter to chemically form metal hydrides. While a number of studies show that oxygen binds strongly as compared to 3 H on the nickel (Ni) layer, the detailed mechanism of 3 H species absorption and diffusion across the Ni plate and Ni/Zr interface are still unclear. By employing density functional theory calculations, here we explored the 3 H 2 and 3 H 2 O species adsorption and dissociation on the Ni(111) surface and diffusion into the Ni sublayer. Our results indicated that the 3 H 2 and 3 H 2 O dissociate on the Ni(111) surface. The NiO x and Ni(O 3 H) x could be formed in the Ni layer due to the higher oxygen (O) diffusion energy barrier and formation of Ni vacancy defects. The oxygen was found to be retained in the Ni layer from diffusing across the Ni–Zr interface. This was revealed by comparing the diffusion barriers for 3 H with O. 3 H was found to have nearly three times smaller diffusion barrier than for O, making 3 H comparatively easier to diffuse through the Ni layer. In conclusion, the obtained results provide guidelines for experimental measurements on 3 H retention behavior in TPBARs and may open further avenues to explore the impurity effects on 3 H diffusion and storage at the Ni/zircaloy interfaces.

Tafen, De Nyago [National Energy Technology Lab. (↗
NASA NTRS1998

EUPDF: Eulerian Monte Carlo Probability Density Function Solver for Applications With Parallel Computing, Unstructured Grids, and Sprays

The success of any solution methodology used in the study of gas-turbine combustor flows depends a great deal on how well it can model the various complex and rate controlling processes associated with the spray's turbulent transport, mixing, chemical kinetics, evaporation, and spreading rates, as well as convective and radiative heat transfer and other phenomena. The phenomena to be modeled, which are controlled by these processes, often strongly interact with each other at different times and locations. In particular, turbulence plays an important role in determining the rates of mass and heat transfer, chemical reactions, and evaporation in many practical combustion devices. The influence of turbulence in a diffusion flame manifests itself in several forms, ranging from the so-called wrinkled, or stretched, flamelets regime to the distributed combustion regime, depending upon how turbulence interacts with various flame scales. Conventional turbulence models have difficulty treating highly nonlinear reaction rates. A solution procedure based on the composition joint probability density function (PDF) approach holds the promise of modeling various important combustion phenomena relevant to practical combustion devices (such as extinction, blowoff limits, and emissions predictions) because it can account for nonlinear chemical reaction rates without making approximations. In an attempt to advance the state-of-the-art in multidimensional numerical methods, we at the NASA Lewis Research Center extended our previous work on the PDF method to unstructured grids, parallel computing, and sprays. EUPDF, which was developed by M.S. Raju of Nyma, Inc., was designed to be massively parallel and could easily be coupled with any existing gas-phase and/or spray solvers. EUPDF can use an unstructured mesh with mixed triangular, quadrilateral, and/or tetrahedral elements. The application of the PDF method showed favorable results when applied to several supersonic-diffusion flames and spray flames. The EUPDF source code will be available with the National Combustion Code (NCC) as a complete package.

Raju, M. S.↗
DOE OSTI2022

The exact exchange–correlation potential in time-dependent density functional theory: Choreographing electrons with steps and peaks

The time-dependent exchange–correlation potential has the unusual task of directing fictitious non-interacting electrons to move with exactly the same probability density as true interacting electrons. This has intriguing implications for its structure, especially in the non-perturbative regime, leading to step and peak features that cannot be captured by bootstrapping any ground-state functional approximation. Here, we review what has been learned about these features in the exact exchange–correlation potential of time-dependent density functional theory in the past decade or so and implications for the performance of simulations when electrons are driven far from any ground state.

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

Electronic structure of chromium trihalides beyond density functional theory

In this work, we explore the electronic band structure of freestanding monolayers of chromium trihalides Cr X 3 , X = Cl, Br, I, within an advanced ab initio theoretical approach based on the use of Green's function functionals. We compare the local density approximation with the quasiparticle self-consistent GW (QS GW ) approximation and its self-consistent extension ( QS G W ^ ) by solving the particle-hole ladder Bethe-Salpeter equations to improve the effective interaction W . We show that, at all levels of theory, the valence band consistently changes shape in the sequence Cl → Br → I , and the valence band maximum shifts from the M point to the Γ point. By analyzing the dynamic and momentum-dependent self-energy, we show that QS G W ^ adds to the localization of the systems in comparison with QS GW , thereby leading to a narrower band and reduced amount of halogens in the valence band manifold. Further analysis shows that X = Cl is most strongly correlated, and X = I is least correlated (most bandlike) as the hybridization between Cr d and X p enhances in the direction Cl → Br → I . For CrBr 3 and CrI 3 , we observe remarkable differences between the QS GW and QS G W ^ valence band structures, while their eigenfunctions are very similar. We show that weak perturbations, like moderate strain, weak changes to the d - p hybridization, and adding small U , can flip the valence band structures between these two solutions in these materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Evidence for Ga clusters in β-Ga 2 O 3 from Raman spectroscopy and density functional theory

Monoclinic gallium oxide (β-Ga 2 O 3 ) single crystals have a Raman mode at ∼250 cm −1 that is strongly correlated with free-electron density. Prior work attributed this peak to an electronic excitation of a shallow donor impurity band. However, heavily n-type thin films grown by metalorganic chemical vapor deposition or molecular beam epitaxy do not have the peak. In the present work, an alternate model is proposed: the 250 cm −1 Raman peak arises from Ga clusters, defined as two or more Ga atoms that form Ga–Ga bonds. Raman mapping reveals variations in the frequency that are consistent with a distribution of cluster sizes. The intensity of the peak decreases as the temperature is raised, attributed to melting of the Ga clusters. First-principles calculations indicate that the 250 cm −1 mode is due to Ga–Ga bond-stretching vibrations. As the Fermi energy is raised, the formation of Ga–Ga dimers becomes energetically favorable, explaining the correlation between n-type conductivity and the appearance of the Raman peak.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Simulations of x-ray absorption spectra for CO desorbing from Ru(0001) with transition-potential and time-dependent density functional theory approaches

The desorption of a carbon monoxide molecule from a Ru(0001) surface was studied by means of X-ray Absorption Spectra (XAS) computed with Transition Potential (TP-DFT) and Time Dependent (TD-DFT) DFT methods. By unraveling the evolution of the CO electronic structure upon desorption, we observed that at 2.3 Å from the surface, the CO molecule has already predominantly gas-phase character. While C 1s XAS is quite insensitive to changes in the C–O bond length, the O 1s excitation is very sensitive with the π* coming down in energy upon CO bond stretching, which competes with the increase in orbital energy due to the repulsive interaction with the metallic surface. We show in a systematic way that the TP-DFT method can describe the XAS rather well at the endpoints (chemisorbed and gas phase) but is affected by artificial charge transfer and/or incorrect spin treatment in the transition region in cases like CO, where there are low-lying π* orbitals and large exchange interactions between the core 1s and valence-acceptor π* orbitals. As an alternative, we demonstrate by comparing with experimental data that a linear response approach using TD-DFT employing common exchange-correlation functionals and finite-size clusters can yield a good description of the spectral evolution of the 1s → π* transition with correct spin and gas-to-chemisorbed chemical shifts in good agreement with experiment.

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

Towards accurate orbital-free simulations: A generalized gradient approximation for the noninteracting free energy density functional

For orbital-free ab initio molecular dynamics, especially on systems in extreme thermodynamic conditions, we provide the first pseudo-potential-adapted generalized gradient approximation (GGA) functional for the non-interacting free energy. This is achieved by systematic finite-temperature extension of our recent LKT ground state non-interacting kinetic energy GGA functional (Phys. Rev. B 98, 041111(R) (2018)). We test the performance of the new functional first via static lattice calculations on crystalline aluminum and silicon. Then we compare deuterium equation of state results against both path-integral Monte Carlo and conventional (orbital-dependent) Kohn-Sham results. The new functional, denoted LKTF, outperforms the previous best semi-local free energy functional, VT84F (Phys. Rev. B 88, 161108(R) (2013)), and provides modestly faster simulations. Finally, we discuss subtleties of identification of kinetic and entropic contributions to non-interacting free-energy functionals obtained by extension from ground state orbital-free kinetic energy functionals.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Interactions involved in the adsorption of ethylene glycol and 2-hydroxyethoxide on the Au(111) surface: a Density Functional Theory study

Abstract Context:The monolayers of ethylene glycol and 2-hydroxyethoxide on gold surfaces have been used in hybrid materials as biosensors. In this article, the adsorption of ethylene glycol and 2-hydroxyethoxide on the Au(111) surface was analyzed. For the first system, ethylene glycol on Au(111), there are Au$$\cdot \cdot \cdot $$ · · · O and Au$$\cdot \cdot \cdot $$ · · · H interactions. To the best of our knowledge, the Au$$\cdot \cdot \cdot $$ · · · H interaction has been overlooked until now. However, in this work, there is strong evidence that this interaction is important to stabilize the system. For the second system, the atomic interactions mentioned previously are also predicted, although there is an additional interaction between 2-hydroxyethoxide molecules. Such an interaction induces the link -O-H-O-, with high values of the electron density at the critical points of the corresponding bond path of the O-H interaction. These links suggest the forming of ethylene glycol chains.Methods:The calculations were performed using two exchange-correlation functionals: BEEF-vdW and C09$$_{x}$$ x -vdW; both functionals incorporate dispersion effects within the Kohn-Sham approach in Density Functional Theory as implemented in GPAW code and ASE computational packages. The contacts between the molecules considered in this article and the Au(111) surface were analyzed through the Quantum Theory of Atoms in Molecules implemented in GPUAM code.

Biochemistry & Molecular Biology↗
DOE OSTI2021

Exact time-dependent density-functional theory for nonperturbative dynamics of the helium atom

By inverting the time-dependent Kohn-Sham equation for a numerically exact dynamics of the helium atom, we show that the dynamical step and peak features of the exact correlation potential found previously in one-dimensional models persist for real three-dimensional systems. We demonstrate that the Kohn-Sham and true current-densities differ by a rotational component. The results have direct implications for approximate TDDFT calculations of atoms and molecules in strong fields, emphasizing the need to go beyond the adiabatic approximation, and highlighting caution in quantitative use of the Kohn-Sham current.

74 ATOMIC AND MOLECULAR PHYSICS↗