Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density functional theory development”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
NASA NTRS1977

An altitude-dependent spacecraft charging model

A model for the altitude dependence of the hot plasma parameters responsible for the electrostatic charging of spacecraft was developed. Based upon plasma orbit theory, the directed velocity is a function of the ambient magnetic field flux density. A consequence of this approach is that while the thermal velocity distributions (assumed to be Maxwellian) of the plasma particles are independent of the magnetic field strength (and hence altitude), the particle densities increase with magnetic field strength. Thus, according to this model, while the equilibrium voltage is independent of altitude, the charging current density increases with decreasing altitude. However, the probability of such spacecraft charging decreases with decreasing altitude.

Haffner, J. W.↗
DOE OSTI2021

CHARGED POINT DEFECTS AND DEFECT COMPLEXES IN PHOSPHORENE AND STRUCTURAL PREDICTION OF TAU-11, A DENSITY FUNCTIONAL THEORY APPROACH

With the advent of powerful computers, first principle calculations have become a robust and reliable tool to complement experiments in the study of materials. First principle calculations such as density functional theory can be used to guide experiments to save time and effort when it comes to new and advanced materials discovery. Study of charge defects in materials and crystal structure prediction are such areas in which density functional theory calculations facilitate the prediction of defect formation energy, charge transition levels and elemental site occupations. Experimentally, studying charge defects requires inferring results from various techniques. Likewise, determining site occupancies for similar elements using experimental techniques is challenging and sometimes impossible. In this work, we study the charged intrinsic and extrinsic defects and defect complexes in phosphorene and determine the site occupancies of Al, Fe and Si in the τ11 phase using density functional theory. In the defect work, we calculate the defect formation energy, charge transition level, binding energies and Stokes shift for vacancy, dopant substitution and dopant-vacancy defect complexes in phosphorene. We found that vacancy defect in phosphorene becomes negatively charged in n-type doping and may passivate the dopants and reduce carrier concentration and mobility. For non-metal dopants in phosphorene, we predict that O, S and Si prefer to form dopant-vacancy complexes removing the vacancy defect states from the band gap. Mn dopant-vacancy defect complex exhibits possibility of switching between two magnetic spin states. Lastly, using density functional theory, we complement the results from neutron powder diffraction to determine the site occupancies of Al and Si in τ11 phase which was used to determine the alloying element to develop a high temperature and low density alloy.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Impact of fission product inclusion on phase development in U 3 Si 2 fuel

Due to its high thermal conductivity and uranium density, U 3 S i2 has been considered as a candidate for use as an accident tolerant fuel (ATF). In order to fully assess its suitability and performance as a fuel, the impact of fission products (FPs) on the stability and performance of U 3 Si 2 must be investigated. The interactions of FPs and U 3 Si 2 have had relatively little study until now and require experimental and computational examination. U 3 Si 2 was doped with individual FPs to explore U-Si-FP interactions and phase equilibria that may impact the performance of the ATF during irradiation. Elemental Ce, Mo, Y, or Zr were used to individually dope U 3 Si 2 at a concentration of 5 wt% FP. A diffusion couple of a 1:1 Mo:Zr alloy and U 3 Si 2 was heated to 1200 °C in order to consider the impacts of multiple FPs on the stability and structure of the fuel. Samples were characterized for FP solubility and secondary phase formation using electron microscopy, energy dispersive spectroscopy, and x-ray diffraction. First principles density functional theory calculations complemented the experimental effort to understand FP behavior. Experimental and computational findings were used in the development of a thermodynamic database containing 8 major FPs and their associated silicide phases. Finally, fuel compositions generated from depletion calculations were used to thermodynamically model the equilibrium phases of the fuel undergoing burnup.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

SPARC-X: Quantum simulations at extreme scale - reactive dynamics from first principles

We have developed the massively parallel electronic structure code SPARC-X: a computational framework for performing Kohn-Sham Density Functional Theory (DFT) calculations that can scale linearly with the number of atoms in the system, while being able to leverage petascale and emerging exascale parallel computers to study chemical phenomena at unprecedented length and time scales. SPARC-X exploits a recent breakthrough in electronic structure methodologies: systematically improvable, strictly local, orthonormal, discontinuous real-space bases that efficiently and systematically capture the local chemistry of the system. With further adaptation using new machine-learning techniques and the use of the massively parallel Spectral Quadrature (SQ) electronic structure method, the algorithmic complexity and prefactor associated with DFT calculations involving semilocal as well as hybrid functionals are dramatically reduced. Using petascale computational resources, SPARC-X enables quantum mechanical simulations at length and time scales previously accessible only by empirical approaches, e.g., 1,000,000 atoms for a few picoseconds using semilocal functionals or 1,000 atoms for a few picoseconds using hybrid functionals. Using exascale resources, the sizes and times targeted are two orders of magnitude larger. Such a capability has applications in a wide variety of chemical sciences, including reactive interfaces where large length- and/or long time-scales are needed and traditional force fields fail. This is particularly important in dynamic catalysis, where bond breaking and formation must be understood in detail. We developed, tested, and employed the SPARC-X framework to understand the photocatalytic properties of TiO 2 nanoparticles, revealing finite size effects that cannot be captured with standard model systems or functionals. This integrated development and application strategy ensures that SPARC-X remains a robust, efficient, and scalable software package for quantum simulations on current petascale and emerging exascale computing resources.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1994

Center for modeling of turbulence and transition: Research briefs, 1993

This research brief contains the progress reports of the research staff of the Center for Modeling of Turbulence and Transition (CMOTT) from June 1992 to July 1993. It is also an annual report to the Institute for Computational Mechanics in Propulsion located at Ohio Aerospace Institute and NASA Lewis Research Center. The main objectives of the research activities at CMOTT are to develop, validate, and implement turbulence and transition models for flows of interest in propulsion systems. Currently, our research covers eddy viscosity one- and two-equation models, Reynolds-stress algebraic equation models, Reynolds-stress transport equation models, nonequilibrium multiple-scale models, bypass transition models, joint scalar probability density function models, and Renormalization Group Theory and Direct Interaction Approximation methods. Some numerical simulations (LES and DNS) have also been carried out to support the development of turbulence modeling. Last year was CMOTT's third year in operation. During this period, in addition to the above mentioned research, CMOTT has also hosted the following programs: an eighteen-hour short course on 'Turbulence--Fundamentals and Computational Modeling (Part I)' given by CMOTT at the NASA Lewis Research Center; a productive summer visitor research program that has generated many encouraging results; collaborative programs with industry customers to help improve their turbulent flow calculations for propulsion system designs; a biweekly CMOTT seminar series with speakers from within and without the NASA Lewis Research Center including foreign speakers. In addition, CMOTT members have been actively involved in the national and international turbulence research activities. The current CMOTT roster and organization are listed in Appendix A. Listed in Appendix B are the abstracts of the biweekly CMOTT seminar. Appendix C lists the papers contributed by CMOTT members.

Liou, William W.↗
DOE OSTI2026

Ab initio thermodynamics of Ni and Co incorporation in Mg hydroxide, carbonate, and hydroxycarbonate minerals

Ni and Co are critical elements needed for modern technologies, and a better understanding of the ability of Mg-based minerals to incorporate these elements would benefit strategy development for Ni and Co recovery from mafic and ultramafic deposits. Here, in this work, we performed density functional theory (DFT) calculations of Ni and Co incorporation in six potential products of the carbonation of mafic and ultramafic silicates: brucite (Mg(OH) 2 ), magnesite (MgCO 3 ), nesquehonite (MgCO 3 ⸱3H 2 O), lansfordite (MgCO 3 ⸱5H 2 O), artinite (Mg 2 CO 3 (OH) 2 ⸱3H 2 O), and hydromagnesite (Mg 5 (CO 3 ) 4 (OH) 2 ⸱4H 2 O). The DFT results were used in an ab initio thermodynamics framework to explore the pH 2 O–pCO 2 conditions at which the Mg-based minerals were predicted to be thermodynamically stable and to quantify the Gibbs free energy of Ni and Co substitution at Mg sites. Among the six Mg-based minerals, brucite and magnesite were predicted to have the lowest Ni and Co substitution free energy. An analysis of the effect of temperature indicated that, at low temperature (<100 K), brucite more readily accommodated Ni and Co, while, at higher temperature (>335 K), magnesite more favorably incorporated Ni and Co. Between 100 K and 335 K, Ni was predicted to preferentially substitute for Mg in brucite and Co for Mg in magnesite, thus leading to a driving force for separating Ni and Co in conditions where brucite and magnesite both form. Insights gained in this work could therefore help select experimental conditions that either promote or inhibit incorporation of these critical elements into Mg-based mineral phases.

Critical elements↗
DOE OSTI2025

Separation of Linear and Cyclic Siloxanes in Pure Silica Zeolites

We present a computational assessment of pure-silica zeolites for separating linear and cyclic siloxanes. We developed a force field (FF) for pure silica zeolites, when combined with our previously developed FF for siloxanes using standard Lorentz–Berthelot combining rules, shows good agreement with dispersion-corrected density functional theory calculations. We used molecular dynamics simulations to investigate diffusion of siloxanes in pure silica zeolites and identified a pure silica zeolite with the structure code FAU that enables kinetic separation of linear and cyclic siloxanes. FAU allows the diffusion of linear siloxanes (L2–L6) while excluding cyclic siloxanes. D4 siloxane does not diffuse in any of the investigated zeolites, eliminating the potential of pure silica zeolites to achieve equilibrium-based separations of linear and cyclic siloxanes.

adsorption↗
DOE OSTI2024

Electromagnetic properties of indium isotopes illuminate the doubly magic character of 100 Sn

Understanding the nuclear properties in the vicinity of 100 Sn, which has been suggested to be the heaviest doubly magic nucleus with proton number Z equal to neutron number N, has been a long-standing challenge for experimental and theoretical nuclear physics. In particular, contradictory experimental evidence exists regarding the role of nuclear collectivity in this region of the nuclear chart. Here, we provide further evidence for the doubly magic character of 100 Sn by measuring the ground-state electromagnetic moments and nuclear charge radii of indium (Z = 49) isotopes as N approaches 50 from above using precision laser spectroscopy. Our results span almost the complete range between the two major closed neutron shells at N = 50 and N = 82 and reveal parabolic trends as a function of the neutron number, with a clear reduction towards these two closed neutron shells. Furthermore, a detailed comparison between our experimental results and numerical results from two complementary nuclear many-body frameworks (density functional theory and ab initio methods) exposes deficiencies in nuclear models and establishes a benchmark for future theoretical developments.

100Sn↗
DOE OSTI2025

Bridging material models across scales: An integrated approach to equation of state and molecular dynamics modeling of copper

New uncertainty-aware equation of state (EOS) and electrical conductivity (EC) models for copper have been developed. The multiphase EOS/EC models are fit to experimental solid/liquid EC isobar measurements as well as density-functional theory molecular dynamics (DFT-MD) EC calculations in both expanded and compressed regimes (0.1–16 g/ cm 3 ⁠). The liquid and solid EOS phases were fit to available experimental data along with additional DFT-MD data over the same range as the EC. Leveraging the DFT-MD data, a corresponding machine-learned interatomic potential (MLIAP) for copper was trained using genetic-algorithm optimization. The copper MLIAP was constrained by EOS shock points at high compressions. The final EOS bounded MLIAP proves to be stable over a large density range (approximately 0.1–20 g/ cm 3 ) with good agreement to an isothermal compression curve, shock Hugoniot, and liquid speed of sound measurements at high pressures (100s of GPa).

Acoustic measurements and instrumentation↗
DOE OSTI2020

Temperature-dependent lattice dynamics in iridium

The characterization of simple elemental systems is key to benchmarking first-principles modeling of electronic and vibrational behaviors of materials. A large body of literature has been built for most elemental systems in the periodic table. However, surprisingly little neutron work has been performed to probe the vibrational properties of iridium, likely due to its large neutron absorption cross section. Nonetheless, iridium is of significant scientific and technological interest due to large relativistic electron effects and electron-phonon coupling, particularly in strongly correlated iridate compounds. In this paper, we report temperature-dependent inelastic neutron scattering measurements of the vibrational properties of iridium, from which we extract key thermodynamic properties. To overcome the challenge of the large neutron absorption of iridium, we developed a simple post-processing correction procedure. The measured densities of phonon states compare well with quasiharmonic density functional theory calculations, although the obtained experimental phonon Grüneisen parameters are much larger than expected, reaching as high as γ=4.5, indicating substantial anharmonicity.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Ground electronic state description of thiourea coordination in homoleptic Zn 2+ , Ni 2+ and Co 2+ complexes using sulfur K -edge X-ray absorption spectroscopy

Sulfur K-edge X-ray absorption spectroscopy (XAS) was employed to experimentally characterize the coordinative bond between the thiourea (TU) or thiocarbamide ligand and transition metal (TM) ions Zn 2+ , Co 2+ and Ni 2+ in distorted tetrahedral and octahedral homoleptic coordination environments. Comparisons of XAS spectra of the free TU ligand and [Zn(TU) 4 ] 2+ , [Co(TU) 4 ] 2+ and [Ni(TU) 6 ] 2+ complexes clearly identify spectral features unique to TM 2+ –S(TU) bonding. Quantitative analysis of pre-edge intensities describes the covalency of Ni 2+ —S(TU) and Co 2+ —S(TU) bonding to be at most 21% and 9% as expressed by the S 3p contributions per TM 3d electron hole. Using relevant Ni 2+ complexes with dithiocarbamate and thioether ligands, we evaluated the empirical S 1s → 3p transition dipole integrals developed for S-donor ligands and their dependence on heteroatom substitutions. With the aid of density functional theory-based ground electronic state calculations, we found evidence for the need of using a transition dipole that is dependent on the presence of conjugated heteroatom (N) substitution in these S-donor ligands.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
NASA NTRS1985

Theoretical Studies of the Surface Tension of Liquid Metal System

A major goal of this project is to understand the surface tension and other thermophysical properties of liquid metals and alloys from a fundamental viewpoint. The approach is to calculate these quantities by a first principles technique which combines the statistical-mechanical theory of the liquid state with an electronic pseudopotential theory of electrons in metals. The inhomogeneity of the surface is treated using an ionic-density-functional formalism developed with the support of NASA. Of particular interest are the variation of surface tension with temperature and impurity concentration: such variations strongly influence the types of convection which make take place in a low-gravity environment. Some progress has already been achieved in computing the reduction of surface tension due to the presence of low-surface-tension impurities, and the corresponding surface segregation of such impurities. In the coming year, it is planned to concentrate on the surface properties of materials of particular interest to the MSA program: Si, Ga and GaSn alloys. An additional goal is to gain some theoretical understanding of the high temperature thermophysical properties of liquid metals, particularly high melting point materials which have not been studied extensively from a theoretical viewpoint.

Stroud, D. G.↗
NASA NTRS2021

Utilizing local phase transformation strengthening for nickel-based superalloys

Almost 75 years of research has been devoted to producing superalloys capable of higher operating temperatures in jet turbine engines, and there is an ongoing need to increase operating temperature further. Here, a new disk Nickel-base superalloy is designed to take advantage of strengthening atomic-scale dynamic complexions. This local phase transformation strengthening provides the alloy with a three times improvement in creep strength over similar disk superalloys and comparable strength to a single crystal blade alloy at 760 °C. Ultra-high-resolution chemical mapping reveals that the improvement in creep strength is a result of atomic-scale η (D024) and χ (D019) formation along superlattice stacking faults. To understand these results, the energy differences between the L12 and competing D024 and D019 stacking fault structures and their dependence on composition are computed by density functional theory. This study can help guide researchers to further optimize local phase transformation strengthening mechanisms for alloy development.

Timothy M. Smith↗
DOE OSTI2023

Efficient computational design of two-dimensional van der Waals heterostructures: Band alignment, lattice mismatch, and machine learning

Here, we develop a computational database, website applications (web-apps), and machine-learning (ML) models to accelerate the design and discovery of two-dimensional (2D) heterostructures. Using density functional theory (DFT) based lattice parameters and electronic band energies for 674 nonmetallic exfoliable 2D materials, we generate 226 779 possible bilayer heterostructures. We classify these heterostructures into type-I, -II, and -III systems according to Anderson’s rule, which is based on the relative band alignments of the noninteracting monolayers. We find that type II is the most common and type III the least common heterostructure type. We subsequently analyze the chemical trends for each heterostructure type in terms of the Periodic Table of constituent elements. The band alignment data can also be used for identifying photocatalysts and high-work-function 2D metals for contacts. We validate our results by comparing them to experimental data as well as hybrid-functional predictions. Additionally, we carry out DFT calculations of a few selected systems (MoS 2 /WSe 2 , MoS 2 /h-BN, and MoSe 2 /CrI 3 ), to compare the band-alignment description with the predictions from Anderson’s rule. We develop web-apps to enable users to virtually create combinations of 2D materials and predict their properties. Additionally, we use ML tools to predict band-alignment information for 2D materials. The web-apps, tools, and associated data will be distributed through the JARVIS-HETEROSTRUCTURE website. Our analysis, results, and the developed web-apps can be applied to the screening and design applications, such as finding alternative photocatalysts, photodetectors, and high-work-function (WF) 2D-metal contacts.

2-dimensional systems↗
DOE OSTI2020

Wannier–Koopmans method calculations for transition metal oxide band gaps

The widely used density functional theory (DFT) has a major drawback of underestimating the band gaps of materials. Wannier–Koopmans method (WKM) was recently developed for band gap calculations with accuracy on a par with more complicated methods. WKM has been tested for main group covalent semiconductors, alkali halides, 2D materials, and organic crystals. Here we apply the WKM to another interesting type of material system: the transition metal (TM) oxides. TM oxides can be classified as either with d 0 or d 10 closed shell occupancy or partially occupied open shell configuration, and the latter is known to be strongly correlated Mott insulators. We found that, while WKM provides adequate band gaps for the d0 and d10 TM oxides, it fails to provide correct band gaps for the group with partially occupied d states. This issue is also found in other mean-field approaches like the GW calculations. We believe that the problem comes from a strong interaction between the occupied and unoccupied d-state Wannier functions in a partially occupied d-state system. We also found that, for pseudopotential calculations including deep core levels, it is necessary to remove the electron densities of these deep core levels in the Hartree and exchange–correlation energy functional when calculating the WKM correction parameters for the d-state Wannier functions.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Tackling Disorder in γ‐Ga 2 O 3

Abstract Ga 2 O 3 and its polymorphs are attracting increasing attention. The rich structural space of polymorphic oxide systems such as Ga 2 O 3 offers potential for electronic structure engineering, which is of particular interest for a range of applications, such as power electronics. γ‐Ga 2 O 3 presents a particular challenge across synthesis, characterization, and theory due to its inherent disorder and resulting complex structure–electronic‐structure relationship. Here, density functional theory is used in combination with a machine‐learning approach to screen nearly one million potential structures, thereby developing a robust atomistic model of the γ‐phase. Theoretical results are compared with surface and bulk sensitive soft and hard X‐ray photoelectron spectroscopy, X‐ray absorption spectroscopy, spectroscopic ellipsometry, and photoluminescence excitation spectroscopy experiments representative of the occupied and unoccupied states of γ‐Ga 2 O 3 . The first onset of strong absorption at room temperature is found at 5.1 eV from spectroscopic ellipsometry, which agrees well with the excitation maximum at 5.17 eV obtained by photoluminescence excitation spectroscopy, where the latter shifts to 5.33 eV at 5 K. This work presents a leap forward in the treatment of complex, disordered oxides and is a crucial step toward exploring how their electronic structure can be understood in terms of local coordination and overall structure.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Influence of elevated temperature and oxygen on the capture of radioactive iodine by silver functionalized silica aerogel

Reprocessing is considered a competent strategy for spent nuclear fuel management, yet radioiodine ( 129 I) is emitted in reprocessing off-gas as a hazardous byproduct. Silver functionalized silica aerogel (Ag 0 -aerogel), a promising iodine capture material, experiences a reduction in its capacity after prolonged exposure to off-gas components at elevated temperatures, a phenomenon termed as aging. To fully understand this process, we isolated the contribution of each aging factor, exposing Ag 0 -aerogel samples to N 2 and dry air gas streams, respectively, at 150 °C for different time periods. Aged samples were loaded with I 2 to examine the capacity change and comprehensively characterized to investigate the evolution of their properties. Results show that temperature alone did not alter Ag 0 -aerogel's capacity but triggered Ag 0 nanoparticles sintering and generated organic sulfur species. The presence of O 2 reduced the capacity by ~20 %, causing (i) formation of silver sulfide (Ag 2 S) crystals and (ii) oxidation of Ag-thiolate (Ag-S-r) to Ag sulfonate (Ag-SO 3 -r). Given that Ag 2 S readily adsorbs I 2 , the formation of Ag-SO 3 -r is the major inhibitor for iodine adsorption. This hypothesis was supported by density functional theory (DFT) simulations. These findings unraveled key mechanisms of Ag 0 -aerogel aging, which are useful in the development of materials that withstand realistic spent-nuclear-fuel-reprocessing off-gas conditions.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Accurate Prediction of pKb in Amines: Validation of the CAM-B3LYP/6-311+G(d,p)/SMD Model

Amines play several key roles in chemistry and biology and are involved in numerous industrial processes, often with significant economic impacts. Recently, amines are also garnering interest as catalysts for polymer synthesis and for CO 2 fixation, incentivizing the need to rapidly design and screen new amino compounds. Hence, developing reliable methods to predict their physicochemical properties, e.g., the base dissociation constant (pKb), is pivotal. Here, a density functional theory (DFT)-based approach was employed to compute the pKb of substituted amines, exploring the impact of several key parameters, including (i) the number of explicit water molecules at the reaction center, (ii) the van der Waals (vdW) surface, and (iii) solvent polarizability. In previous work, it was determined that including two explicit water molecules at the reaction center resulted in highly accurate pKb estimates for primary amines. Here, we find that including a third water molecule at the reaction center is essential for accurate pKb for secondary and tertiary amines. The revised methodology was then applied to a wider selection of amines, obtaining a minimum average error (MAE) < 0.4. In conclusion, this result represents an extension of our “easy-to-use method,” a simple and direct DFT approach exploiting CAM-B3LYP/SMD/6-311G+(d,p) to compute pKb without post facto modifications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗