Search NASA⌕ Search

SEARCH · Search NASA

Results for “Density functional theory simulations”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Ab Initio Simulations of Tritium Diffusion in Al 2 O and Intermetallic Al 12 (TM) 2.34 Aluminide Coating Phases

Density functional theory simulations have been carried out to investigate the diffusion of interstitial tritium in Al 12 (TM) 2.34 and Al 2 O bulk phases. In Al 12 (TM) 2.34 the transition metal (TM) sites are occupied on average by 58.93 at.% Fe, 18.52 at.% Cr, and 22.54 at.% Ni. While the insertion of interstitial tritium in Al 2 O lead to a strong disordering of the structure, we only investigated interstitial tritium in Al 12 (TM) 2.34 and its iron end-member (i.e., Al 12 Fe 2.34 ). Nine diffusion pathways have been investigated along the a-, b-, and c -axis of Al 12 (TM) 2.34 . The direction of fastest diffusion for interstitial tritium is found to be along the b-axis, with a calculated diffusion coefficient of ≈10 -10 m 2 .s -1 at 600 K, which is at least one order of magnitude faster than those previously calculated for interstitial tritium in other Al-rich phases such as Fe2Alx, Fe4Al13, and FeNiAl5 for which D T =10 -11 m 2 .s -1 , ≈10 -12 m 2 .s -1 , and D T ≈10 -13 m 2 .s -1 respectively. The comparison of the energy landscape between Al 12 (TM) 2.34 and Al 12 Fe 2.34 for the fastest diffusion pathways in each axis direction, found that some pathways are not sensitive to transition metal mixing, while other are more affected, leading to higher energy barrier in Al 12 Fe 2.34 in part due to a more energetically favorable formation of Fe—T bond compared to Ni—T. However, we found that both materials have similar diffusion coefficients as the fastest diffusion pathways occurs along pathways that are not very sensitive to transition metal mixing.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Toward first principles-based simulations of dense hydrogen

Accurate knowledge of the properties of hydrogen at high compression is crucial for astrophysics (e.g., planetary and stellar interiors, brown dwarfs, atmosphere of compact stars) and laboratory experiments, including inertial confinement fusion. There exists experimental data for the equation of state, conductivity, and Thomson scattering spectra. However, the analysis of the measurements at extreme pressures and temperatures typically involves additional model assumptions, which makes it difficult to assess the accuracy of the experimental data rigorously. On the other hand, theory and modeling have produced extensive collections of data. They originate from a very large variety of models and simulations including path integral Monte Carlo (PIMC) simulations, density functional theory (DFT), chemical models, machine-learned models, and combinations thereof. At the same time, each of these methods has fundamental limitations (fermion sign problem in PIMC, approximate exchange–correlation functionals of DFT, inconsistent interaction energy contributions in chemical models, etc.), so for some parameter ranges accurate predictions are difficult. Recently, a number of breakthroughs in first principles PIMC as well as in DFT simulations were achieved which are discussed in this review. Here we use these results to benchmark different simulation methods. We present an update of the hydrogen phase diagram at high pressures, the expected phase transitions, and thermodynamic properties including the equation of state and momentum distribution. Furthermore, we discuss available dynamic results for warm dense hydrogen, including the conductivity, dynamic structure factor, plasmon dispersion, imaginary-time structure, and density response functions. We conclude by outlining strategies to combine different simulations to achieve accurate theoretical predictions that are based on first principles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Local ordering in Ge/Ge–Sn semiconductor alloy core/shell nanowires revealed by extended x-ray absorption fine structure (EXAFS)

Short-range atomic order in semiconductor alloys is a relatively unexplored topic that may promote design of new materials with unexpected properties. Here, local atomic ordering is investigated in Ge–Sn alloys, a group-IV system that is attractive for its enhanced optoelectronic properties achievable via a direct gap for Sn concentrations exceeding ≈10 at. %. The substantial misfit strain imposed on Ge–Sn thin films during growth on bulk Si or Ge substrates can induce defect formation; however, misfit strain can be accommodated by growing Ge–Sn alloy films on Ge nanowires, which effectively act as elastically compliant substrates. In this work, Ge core/Ge 1–x Sn x (x ≈ 0.1) shell nanowires were characterized with extended x-ray absorption fine structure (EXAFS) to elucidate their local atomic environment. Simultaneous fitting of high-quality EXAFS data collected at both the Ge K-edge and the Sn K-edge reveals a large (≈ 40%) deficiency of Sn in the first coordination shell around a Sn atom relative to a random alloy, thereby providing the first direct experimental evidence of significant short-range order in this semiconductor alloy system. Furthermore, comparison of path length data from the EXAFS measurements with density functional theory simulations provides alloy atomic structures consistent with this conclusion.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Understanding the fundamental driver of semiconductor radiation tolerance with experiment and theory

Space is the operating environment of a multitude of systems that our society is heavily reliant on today, however, maintaining operability there necessitates special consideration of the electronic systems' tolerance of space radiation. Electronic systems are critically dependent on the electronic properties of their semiconductor components, which are modified by space radiation with an adverse impact on the space system performance. What innate property allows some semiconductors to sustain little damage while others accumulate defects rapidly with dose is poorly understood, which limits the extent to which radiation tolerance can be implemented as a design criterion. Here, to gain insight into what properties are drivers of semiconductor radiation tolerance, the first step is to generate a dataset of the relative radiation tolerance of a broad sampling of semiconductors. To accomplish this, Rutherford backscatter channeling experiments are used to compare the displaced lattice atom buildup in InAs, InP, GaP, GaN, ZnO, MgO, and Si as a function of stepwise alpha particle dose. With this experimental information on radiation-induced incorporation of interstitial defects in hand, hybrid density functional theory electron densities (and their derived quantities) are calculated and their gradient and Laplacian are evaluated to obtain key fundamental information about the interactions in each material. It is shown that simple, undifferentiated values (which are typically used to describe bond strength) are insufficient to predict radiation tolerance. Instead, the curvature of the electron density at bond critical points provides a measure of radiation tolerance consistent with the experimental results obtained. This curvature and associated forces surrounding bond critical points have the potential to disfavor the localization of displaced lattice atoms at these points, favoring their diffusion toward perfect lattice positions. With this criterion to predict radiation tolerance, simple density functional theory simulations can be conducted on potential new materials to gain insight into how they may operate in demanding high radiation environments.

36 MATERIALS SCIENCE↗
DOE OSTI2022

First-principles simulation of light-ion microscopy of graphene

The extreme sensitivity of 2D materials to defects and nanostructure requires precise imaging techniques to verify presence of desirable and absence of undesirable features in the atomic geometry. Helium-ion beams have emerged as a promising materials imaging tool, achieving up to 20 times higher resolution and 10 times larger depth-of-field than conventional or environmental scanning electron microscopes. Here, we offer first-principles theoretical insights to advance ion-beam imaging of atomically thin materials by performing real-time time-dependent density functional theory simulations of single impacts of 10–200 keV light ions in free-standing graphene. Here we predict that detecting electrons emitted from the back of the material (the side from which the ion exits) would result in up to three times higher signal and up to five times higher contrast images, making 2D materials especially compelling targets for ion-beam microscopy. This predicted superiority of exit-side emission likely arises from anisotropic kinetic emission. The charge induced in the graphene equilibrates on a sub-fs time scale, leading to only slight disturbances in the carbon lattice that are unlikely to damage the atomic structure for any of the beam parameters investigated here.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Sulfonated polybenzimidazole membrane with graphene oxide additive for 2,3-butanediol/water separation: A molecular simulation

Membrane separation for 2,3-butanediol (2,3-BDO) recovery from fermentation broth is highly valued for sustainable and renewable processes, but it requires efficient membrane materials. Here, this work evaluates the sulfonated polybenzimidazole (sPBI) and its graphene oxide (GO) doped composite membrane for separating 2,3-BDO and water via atomistic simulations. Density functional theory calculations are applied to identify various forms of sPBI structures and quantify their binding interactions with 2,3-BDO and water. Classical molecular dynamic simulations are used to evaluate the structural changes, diffusivity, and selectivity of 2,3-BDO and water in different sPBI models, GO surfaces, and GO-doped sPBI composite models. Our results suggest that sPBI slightly increases the crystallinity of the membrane structures, enhances the adsorption strength for both 2,3-BDO and water, and improves the water/2,3-BDO selectivity by 2–3 times. The GO surfaces display a maximum selectivity at a surface coverage of 0.1–0.15 for both hydroxyl and epoxy surface groups. The addition of GO flakes to sPBI creates new interaction sites for 2,3-BDO and water at the interface of sPBI and GO, and the water/2,3-BDO selectivity of GO-doped sPBI models is further increased up to 3 times. This work illustrates how the integrated addition of sPBI and GO flakes offers a promising approach to selective separation of 2,3-BDO and water, providing theoretical guidance for polybenzimidazole-based membranes in the potential application of 2,3-BDO recovery.

2,3-butanediol↗
DOE OSTI2026

Discovery of magnetic-field-tunable density modulations and spin tilting in a layered altermagnet

Altermagnets recently emerged as a new class of magnetic materials, arising from specific spin crystal symmetries. They exhibit a spin-polarized electronic band structure similar to ferromagnets, yet possess zero net magnetization, promising exotic properties. Here we study a layered triangular lattice altermagnet, cobalt-intercalated NbSe 2 using scanning tunneling microscopy and spectroscopy (STM/S). Spectroscopic-imaging STM and spin-polarized STM reveals emergent 2 a 0 tri-directional charge and spin density modulations on the selenium surface. Density functional theory simulations suggest these modulations reflect the underlying cobalt superstructure. We discover that an out-of-plane magnetic field tunes the modulation amplitudes and the electronic density-of-states in a manner dependent on the field direction and strength. This behavior is attributed to the field-induced tilting of cobalt spins, which can have profound implications on the electronic properties of the altermagnet. Our results provide atomic-scale insights to uncover a magnetic-field tunable altermagnetic band structure, highlight the importance of understanding spin canting in altermagnets.

condensed-matter physics↗
DOE OSTI2024

Influence of temperature, doping, and amorphization on the electronic structure and magnetic damping of iron

Hybrid magnonic quantum systems have drawn increased attention in recent years for coherent quantum information processing, but too large magnetic damping is a persistent concern when metallic magnets are used. Their intrinsic damping is largely determined by electron-magnon scattering induced by spin-orbit interactions. In the low scattering limit, damping is dominated by intra-band electronic transitions, which has been theoretically shown to be proportional to the electronic density of states at the Fermi level. In this work, we focus on body-centered-cubic iron as a paradigmatic ferromagnetic material. We comprehensively study its electronic structure using first-principles density functional theory simulations and account for finite lattice temperature, boron (B) doping, and structure amorphization. Our results indicate that temperature induced atomic disorder and amorphous atomic geometries only have a minor influence. Instead, boron doping noticeably decreases the density of states near the Fermi level with an optimal doping level of 6.25%. In addition, we show that this reduction varies significantly for different atomic geometries and report that the highest reduction correlates with a large magnetization of the material. Furthermore, this may suggest materials growth under external magnetic fields as a route to explore in experiment.

36 MATERIALS SCIENCE↗
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 OSTI2022

On the origin of multihole oxygen evolution in haematite photoanodes

The oxygen evolution reaction (OER) plays a crucial role in (photo)electrochemical devices that use renewable energy to produce synthetic fuels. Recent measurements on semiconducting oxides have found a power law dependence of the OER rate on surface hole density, suggesting a multihole mechanism. In this study, using transient photocurrent measurements, density functional theory simulations and microkinetic modelling, we have uncovered the origin of this behaviour in haematite. We show here that the OER rate has a third-order dependence on the surface hole density. We propose a mechanism wherein the reaction proceeds by accumulating oxidizing equivalents through a sequence of one-electron oxidations of surface hydroxy groups. The key O–O bond formation step occurs by the dissociative chemisorption of a hydroxide ion involving three oxyl sites. At variance with the case of metallic oxides, the activation energy of this step is weakly dependent on the surface hole coverage, leading to the observed power law.

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

The impact of nickel concentration and stacking fault energy on deformation mechanisms in high-purity austenitic Fe-Cr-Ni alloys

Understanding how composition affects deformation mechanisms in austenitic stainless steels is essential for developing accurate predictive models of stress-induced failures and stress corrosion cracking. Nickel (Ni), an element classified as a critical element, plays a crucial role in these processes. It is important to examine how Ni concentration influences stacking fault energy (SFE) and, consequently, the deformation mechanisms of austenitic stainless steels. However, in commercial stainless steels, the effects of other alloying elements and impurities can obscure Ni's role, complicating efforts to isolate its impact. Here, in this study, we use two high-purity Fe-Cr-Ni alloys to investigate how Ni concentration and SFE interact to alter deformation mechanisms and induce martensitic transformation. By combining in situ synchrotron X-ray diffraction (XRD) tensile testing and post-mortem electron microscopy with density functional theory simulations, we gain precise insights into these phenomena. We find that the Fe18Cr10Ni (wt%) alloy, with its low SFE, exhibits higher stacking fault probability, deformation-induced martensitic transformation, and a lesser increase in dislocation density with plastic strain. In contrast, the Fe18Cr14Ni (wt%) alloy, with its higher SFE, shows enhanced deformation twinning and greater dislocation density with increasing strain. These findings from high-purity ternary alloys provide valuable insights that can guide the search for alternative elements to replace Ni while achieving similar effects on phase stability and deformation behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Two-Dimensional Conductive Ni-HAB as a Catalyst for the Electrochemical Oxygen Reduction Reaction

Catalytic systems whose properties can be systematically tuned via changes in synthesis conditions are highly desirable for the next-generation catalyst design and optimization. Herein, we present a two-dimensional (2D) conductive metal–organic framework consisting of M-N 4 units (M = Ni, Cu) and a hexaaminobenzene (HAB) linker as a catalyst for the oxygen reduction reaction. By varying synthetic conditions, we prepared two Ni-HAB catalysts with different crystallinities, resulting in catalytic systems with different electric conductivities, electrochemical activity, and stability. We show that crystallinity has a positive impact on conductivity and demonstrate that this improved crystallinity/conductivity improves the catalytic performance of our model system. Additionally, density functional theory simulations were performed to probe the origin of M-HAB’s catalytic activity, and they suggest that M-HAB’s organic linker acts as the active site with the role of the metal being to modulate the linker sites’ binding strength.

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

Nanoprecipitates to Enhance Radiation Tolerance in High-Entropy Alloys

The growth of advanced energy technologies for power generation is enabled by the design, development, and integration of structural materials that can withstand extreme environments, such as high temperatures, radiation damage, and corrosion. High-entropy alloys (HEAs) are a class of structural materials in which suitable chemical elements in four or more numbers are mixed to typically produce single-phase concentrated solid solution alloys (CSAs). Many of these alloys exhibit good radiation tolerance like limited void swelling and hardening up to relatively medium radiation doses (tens of displacements per atom (dpa)); however, at higher radiation damage levels (>50 dpa), some HEAs suffer from considerable void swelling limiting their near-term acceptance for advanced nuclear reactor concepts. In this study, we developed a HEA containing a high density of Cu-rich nanoprecipitates distributed in the HEA matrix. The Cu-added HEA, NiCoFeCrCu 0.12 , shows excellent void swelling resistance and negligible radiation-induced hardening upon irradiation up to high radiation doses (i.e., higher than 100 dpa). The void swelling resistance of the alloy is measured to be significantly better than NiCoFeCr CSA and austenitic stainless steels. Density functional theory simulations predict lower vacancy and interstitial formation energies at the coherent interfaces between Cu-rich nanoprecipitates and the HEA matrix. The alloy maintained a high sink strength achieved via nanoprecipitates and the coherent interface with the matrix at a high radiation dose (~50 dpa). From our experiments and simulations, the effective recombination of radiation-produced vacancies and interstitials at the coherent interfaces of the nanoprecipitates is suggested to be the critical mechanism responsible for the radiation tolerance of the alloy. Furthermore, the materials design strategy based on incorporating a high density of interfaces can be applied to high-entropy alloy systems to improve their radiation tolerance.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Tuning spin–orbit coupling in (6,5) single-walled carbon nanotube doped with sp 3 defects

Single-walled carbon nanotubes (SWCNTs) containing sp 3 defects are a promising class of optoelectronic materials with bright photoluminescence and demonstrated single-photon emission. Using density functional theory simulations, complemented by measurements, we investigate the electronic structure of a series of quantum defects attached to (6,5) SWCNT with the goal of tuning the spin–orbit coupling by introduction of a heavy atom in the defect structure. We characterize the ground state electronic and spin properties of four synthesized and three potential defects on the tube and find that all of the synthesized defects considered introduce a localized midgap defect-centered state containing a single electron, ≈0.2–0.3 eV above the valence band. The spin density is located at the sp 3 defect site with negligible spin–orbit coupling even with the presence of a Pd atom. Three additional functional groups were tested via computation to increase spin localization near the metal, thereby increasing spin–orbit coupling. We predict that only the chlorodiphosphanepalladium(II)– [Cl(PH 3 ) 2 Pd(II)–] defect results in increased spin–orbit splitting of the defect state and the conduction band associated with the pristine-like SWCNT, a measure of the spin–orbit coupling of excited state transitions. This study suggests that for unpassivated sp 3 defects in (6,5) SWCNT, forming a direct bond between a heavy atom and the sp 3 carbon allows for tuning of spin–orbit coupling.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Solvothermal Synthesis of [Cr 7 S 8 (en) 8 Cl 2 ]Cl 3 ⋅ 2H 2 O with Magnetically Frustrated [Cr 7 S 8 ] 5+ Double‐Cubes**

Abstract A novel transition metal chalcohalide [Cr 7 S 8 (en) 8 Cl 2 ]Cl 3 ⋅ 2H 2 O, with [Cr 7 S 8 ] 5+ dicubane cationic clusters, has been synthesized by a low temperature solvothermal method, using dimethyl sulfoxide (DMSO) and ethylenediamine ( en ) solvents. Ethylenediamine ligand exhibits bi‐ and monodentate coordination modes; in the latter case ethylenediamine coordinates to Cr atoms of adjacent clusters, giving rise to a 2D polymeric structure. Although magnetic susceptibility shows no magnetic ordering down to 1.8 K, a highly negative Weiss constant, θ =−224(2) K, obtained from Curie‐Weiss fit of inverse susceptibility, suggests strong antiferromagnetic (AFM) interactions between S =3/2 Cr(III) centers. Due to the complexity of the system with (2 S +1) 7 =16384 microstates from seven Cr 3+ centers, a simplified model with only two exchange constants was used for simulations. Density‐functional theory (DFT) calculations yielded the two exchange constants to be J 1 =−21.4 cm −1 and J 2 =−30.2 cm −1 , confirming competing AFM coupling between the shared Cr 3+ center and the peripheral Cr 3+ ions of the dicubane cluster. The best simulation of the experimental data was obtained with J 1 =−20.0 cm −1 and J 2 =−21.0 cm −1 , in agreement with the slightly stronger AFM exchange within the triangles of the peripheral Cr 3+ ions as compared to the AFM exchange between the central and peripheral Cr 3+ ions. This compound is proposed as a synthon towards magnetically frustrated systems assembled by linking dicubane transition metal‐chalcogenide clusters into polymeric networks.

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

Multiple incommensurate magnetic states in the kagome antiferromagnet Na 2 Mn 3 Cl 8

The kagome lattice can host exotic magnetic phases arising from frustrated and competing magnetic interactions. However, relatively few insulating kagome materials exhibit incommensurate magnetic ordering. Here, we present a study of the magnetic structures and interactions of antiferromagnetic Na 2 Mn 3 Cl 8 with an undistorted Mn 2+ kagome network. Using neutron-diffraction and bulk magnetic measurements, we show that Na 2 Mn 3 Cl 8 hosts two different incommensurate magnetic states, which develop at T N1 = 1.6 K and T N2 = 0.6 K. Magnetic Rietveld refinements indicate magnetic propagation vectors of the form q=(q x , q y , $\frac{3}{2}$), and our neutron-diffraction data can be well described by cycloidal magnetic structures. By optimizing exchange parameters against magnetic diffuse-scattering data, we show that the spin Hamiltonian contains ferromagnetic nearest-neighbor and antiferromagnetic third-neighbor Heisenberg interactions, with a significant contribution from long-ranged dipolar coupling. This experimentally determined interaction model is compared with density-functional-theory simulations. Using classical Monte Carlo simulations, we show that these competing interactions explain the experimental observation of multiple incommensurate magnetic phases and may stabilize multi-q states. Here, our results expand the known range of magnetic behavior on the kagome lattice.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗