Search NASA⌕ Search

SEARCH · Search NASA

Results for “Face-centered cubic crystals”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Iron isotope fractionation between solid and liquid metal in the Fe-P±Ni system: Experimental constraints and implications for meteorites

Iron meteorites record a range of Fe isotope compositions that hold valuable information regarding the evolution of their parent bodies. Interpreting this isotopic variability, however, requires experimental constraints on the equilibrium isotope fractionation between phases. It is thought that the cores of many iron meteorite parent bodies experienced fractional crystallization, during which crystallization of solid iron-nickel occurs from an increasingly non-metal-rich liquid alloy. Phosphorus is one component of this alloy, and this study provides the first constraints on Fe-isotope fractionation between solid and liquid alloys in the Fe-Ni-P system. Experiments comprising Fe and P show a clear enrichment in the light isotopes of Fe in the liquid phase, which increases with the amount of phosphorus. Nickel-bearing samples are offset from the trend defined by Ni-free experiments, which is accounted for by the change in the solid alloy phase from a body-centered cubic to face-centered cubic structure upon the addition of Ni. The increasing light isotope enrichment of the liquid with increasing P content suggests interstitial solution of P, which is known to lengthen Fe-Fe bonds in Fe-P liquids (Waseda and Shiraishi 1977). Results suggest a negligible effect of P on Fe isotope fractionation during planetesimal core crystallization. Iron isotopes may, however, prove useful for identifying the petrogenesis of schreibersite in pallasites and iron meteorites.

58 GEOSCIENCES↗
NASA NTRS1975

Aluminum under high pressure. I - Equation of state

A curve of applied pressure versus lattice constant is calculated for single-crystal aluminum. It results from an application of the method of structural expansions for deriving the energies of simple metals, a method known to give reasonable results for the elastic constants even at second order in the effective electron-ion interaction. The latter is taken from Fermi-surface analysis, and it is verified that the extant face-centered cubic structure remains the preferred crystalline phase up to the highest pressures considered. Arguments are given to suggest that the curve should have reasonable a priori accuracy and can admit possible improvement if experimental data in the intermediate-pressure region can be provided to refine the (in principle) energy-dependent pseudopotential. At three megabars, the lattice constant is reduced by only 22 per cent; the ion cores at this pressure are still very well separated.

Friedli, C.↗
NASA NTRS1975

Aluminum under high pressure: 1 equation of state

A curve of applied pressure P versus lattice constant a is calculated for single crystal aluminum. It results from an application of the method of structural expansions for deriving the energies of simple metals, a method known to give reasonable results for the elastic constants even at second order in the effective electron-ion interaction. The method is taken from Fermi surface analysis . This essentially experimental information is used to verify that the extant face-centered cubic structure remains the preferred crystalline phase up to the highest pressures considered. Arguments are given to suggest that the P versus a curve should have reasonable a priori accuracy, which may be improved if experimental data in the intermediate pressure region can be provided to refine the energy dependent pseudo-potential. At 3 megabars the lattice constant is reduced by only 22%; the ion cores at this pressure are still very well separated.

Friedli, C.↗
DOE OSTI2022

A high-dimensional neural network potential for molecular dynamics simulations of condensed phase nickel and phase transitions

A high-dimensional neural network interatomic potential was developed and used in molecular dynamics simulations of condensed phase Ni and Ni systems with liquid–solid phase coexistence. The reference data set was generated by sampling the potential energy surface over a broad temperature-pressure domain using ab initio MD simulations to train a unified potential. Excellent agreement was achieved between bulk face-centered cubic nickel thermal expansion simulations and relevant experimental data. The same potential also yields accurate structures and diffusivities in the liquid state. The phase transition between liquid and solid phases was simulated using the two-phase interface method. The predicted melting point temperature is within a few kelvins of the literature value. Here, the general methodology could be applied to describe crystals with much more complex phase behaviors.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Experimental evidence of disordered crystalline premixing in sputter-deposited Ni(V)/Al multilayers

The sputter deposition of alternating layers of Ni(V) and Al forms a reactive multilayer known to undergo self-propagating formation reactions when ignited. The sequential deposition process leads to nanometer-scale premixing of reactants at each included interface, which ultimately affects multilayer exothermicity. This work performs the direct measurement of a disordered face-centered cubic (FCC) solid solution premixed phase at the interfaces of Ni(V)/Al multilayers via scanning transmission electron microscopy. The crystallinity of the observed phase differs from previously reported a priori predictions of an amorphous interlayer. The disordered FCC phase retains its symmetry after annealing for 16 h at 135 ± 5 °C, but the lattice parameter shifts consistently with an Al-rich composition. The existence of a crystalline premix in Ni(V)/Al is attributed to the electronic contribution to the entropy of crystallization. The importance of electronic entropy to the phase formation of energetic materials motivates its inclusion when constructing digital twins for atomistic kinetics and ignition sensitivity.

Crystal structure↗
DOE OSTI2024

In situ measurement of three-dimensional intergranular stress localizations and grain yielding under elastoplastic axial-torsional loading

The three-dimensional grain-averaged response of solid bar samples under non-proportional (NP) elastoplastic axial-torsional loading was investigated using in situ high energy diffraction microscopy (HEDM) and companion crystal plasticity finite element (CPFE) modeling. Important stress metrics including applied shear (σ θZ ) and axial (σ ZZ ) stress tensor components, stress and stress deviator tensor invariants (I 1 , J 2 , and J 3 ), von Mises equivalent stress (σ$^{grain}_{VM}$), maximum resolved shear stress (mRSS), stress triaxiality (η), and lode angle parameter ($\barθ$) values were tracked for ~300 grains under two different loading conditions: (1) Torsion-dominated loading (low NP) and (2) Tension-torsion loading (high NP) in equiatomic NiCoCr, a representative multicomponent face-centered cubic (FCC) superalloy. Overall, significant stress localizations existed within both samples as evidenced by the radial dependence of grain-resolved σ θZ , σ$^{grain}_{VM}$, and J 2 ; by comparison, I 1 , J 3 , η, and $\barθ$ metrics did not show discernible trends within the volume. These stress localizations reveal a complex interplay between axial and shear stress components (e.g., stress coupling) resulting in grain yielding near the sample surface largely driven by shear stress, whereas internal grain yielding was largely accommodated by axial stress. Grain-resolved stress localization trends were described well by the CPFE model, although some discrepancies in magnitude occurred, particularly for volumetric stress metrics (I 1 and η) due to initial type II residual stress distributions. The superposition of initial residual stress states onto CPFE grain-resolved data significantly improved model accuracy for η. This suggests that residual stresses more strongly influence the simulation of volumetric rather than deviatoric (yield) stress metrics.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Dataset of simulated vibrational density of states and X-ray diffraction profiles of mechanically deformed and disordered atomic structures in Gold, Iron, Magnesium, and Silicon

This dataset is comprised of a library of atomistic structure files and corresponding X-ray diffraction (XRD) profiles and vibrational density of states (VDoS) profiles for bulk single crystal silicon (Si), gold (Au), magnesium (Mg), and iron (Fe) with and without disorder introduced into the atomic structure and with and without mechanical loading. Included with the atomistic structure files are descriptor files that measure the stress state, phase fractions, and dislocation content of the microstructures. All data was generated via molecular dynamics or molecular statics simulations using the Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) code. This dataset can inform the understanding of how local or global changes to a materials microstructure can alter their spectroscopic and diffraction behavior across a variety of initial structure types (cubic diamond, face-centered cubic (FCC), hexagonal close-packed (HCP), and body-centered cubic (BCC) for Si, Au, Mg, and Fe, respectively) and overlapping changes to the microstructure (i.e., both disorder insertion and mechanical loading).

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗
DOE OSTI2026

Quantifying the dislocation content of atomically resolved grain boundary line defects using the Nye tensor

The Nye tensor, which quantifies the density of Burgers vector for a given dislocation line direction, can be effectively used to characterize dislocation content in bulk crystals from atomic-resolution transmission electron microscopy images. The Nye tensor can be calculated from these images, in part because the reference state is simply defined by the lattice of the perfect crystal. The application of the Nye tensor to interfacial line defects, for which the natural reference state is the dichromatic pattern of the two grains in their reference orientation, poses additional challenges. In this work, we present a method that employs the Nye tensor to characterize the edge dislocation content of line defects at grain boundaries from atomic-resolution images. This approach enables us to rapidly characterize all edge dislocation content along a grain boundary. Additionally, the Nye tensor provides information about line defect core structure. Finally, we demonstrate this method on two exemplar defects: a twin boundary disconnection and a facet junction in face-centered cubic Au.

Crystallographic defects↗
DOE OSTI2022

ORNL_AISD_NiNb

This dataset describes the nickel-niobium solid solution binary alloy, where the two constituent elements nickel (Ni) and niobium (Nb) are randomly placed on an underlying crystal lattice. This dataset for nickel-niobium (Ni-Nb) alloys available includes the formation energy and bulk modulus for each crystal structure. Each atomic sample has a disordered phase which is obtained starting from an initial regular crystal structure of type body-centered cubic (BCC), face-centered cubic (FCC), or hexagonal compact packed (HCP). The geometry optimization ensures that all the alloy samples reached the equilibrium with negative formation energy. We perform geometry optimizations using the LAMMPS simulation package [1], a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. We utilized the embedded atom model (EAM) potential for Ni and Nb developed in a previous study [2]. The potential could describe behaviors of the liquid and solid phases of Ni-Nb alloy. The structural factors and angular distributions of three atoms are well-matched with X-ray and ab initio-based molecular dynamics data. We prepared the three different crystals with different initial lattice parameters (3.52 Ã… for FCC, 3.32 Ã… for BCC, and 3.5 Ã… for HCP). We performed energy minimization in two steps. Firstly, we minimized the structures with an isotropic unit cell to minimize the side effects from our arbitrary lattice parameters for all other compositions. Then, we applied geometry optimization with a triclinic (non-orthogonal) unit cell to fully minimize the stress components to calculate the elastic constants. In this procedure, we chose 10,000 as the maximum number of allowable steps aimed at obtaining fully relaxed atomic geometries. The dataset consists of three sets of crystal structures. The first set contains 46,086 irregular crystal structures, each of them with 54 atoms, obtained through optimization starting from a regular BCC crystal structure. The second set contains 24,543 irregular crystal structures, each of them with 32 atoms, obtained through optimization starting from a regular FCC crystal structure. The third set contains 39,303 irregular crystal structures, each of them with 48 atoms, obtained through optimization starting from a regular HCP crystal structure. The atomic configurations within each set span the possible compositional range. The three sets have been unified in a global dataset, which is extremely heterogeneous in terms of crystal structures, lattice volumes, and atomic configurations. Organization of files inside the dataset: the dataset contains three subdirectories called • BCC_opt • FCC_opt • HCP_opt based on the type of initial regular structure used to start the geometry optimization. Inside each of these folders, every atomic structure is identified by a string “A_B_Câ€, where A denotes the number of Nb in the system, B denotes index of structure with a given Nb number, and C denotes the total number of structures generated with a given Nb number. For each optimized crystal structure identified by the unique string of characters “A_B_Câ€, three files are provided: • A_B_C_opt.xyz: The optimized geometries in xyz format • A_B_C_opt.cfg: The optimized geometries in cfg format. It includes cell information and atomic energy, and forces calculated from LAMMPS. • A_B_C.elastic: Raw data of 21 elastic constants from LAMMPS output. • A_B_C.bulk: Calculated upper and lower bounds of bulk modulus and averaged one based on Voigt-Reuss-Hill approach from *.elastic. References: [1] A. P. Thompson, H. M. Aktulga, R. Berger, D. S. Bolintineanu, W. M. Brown, P. S. Crozier, P. J. in 't Veld, A. Kohlmeyer, S. G. Moore, T. D. Nguyen, R. Shan, M. J. Stevens, J. Tranchida, C. Trott, and S. J. Plimpton. LAMMPS - a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. Comp. Phys. Comm., 271:108171, 2022. [2] Y Zhang, R Ashcraft, MI Mendelev, CZ Wang, and KF Kelton. Experimental and molecular dynamics simulation study of structure of liquid and amorphous ni62nb38 alloy. The Journal of chemical physics, 145(20):204505, 2016.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Evaluation of Equiatomic CrMnFeCoNiCu System and Subsequent Derivation of a Non-Equiatomic MnFeCoNiCu Alloy

Investigation into non-equiatomic high-entropy alloys has grown in recent years due to questions about the role of entropy stabilization in forming single-phase solid solutions. Non-equiatomic alloys have been shown to retain the outstanding mechanical properties exhibited by their equiatomic counterparts and even improve electrical, thermal, and magnetic properties, albeit with relaxed composition bounds. However, much remains to understand the processing–structure–property relationships in all classes of so-called high-entropy alloys (HEAs). Here, we are motivated by the natural phenomena of crystal growth and equilibrium conditions to introduce a method of HEA development where controlled processing conditions determine the most probable and stable composition. This is demonstrated by cooling an equiatomic CrMnFeCoNiCu alloy from the melt steadily over 3 days (cooling rate ~4 °C/h). The result is an alloy containing large Cr-rich precipitates and an almost Cr-free matrix exhibiting compositions within the MnFeCoNiCu system (with trace amounts of Cr). From this juncture, it is argued that the most stable composition is within the CrMnFeCoNiCu system rather than the CrMnFeCoNi system. With further optimization and evaluation, a unique non-equiatomic alloy, Mn 17 Fe 21 Co 24 Ni 24 Cu 14 , is derived. The alloy solidifies and recrystallizes into a single-phase face-centered cubic (FCC) polycrystal. In addition to possible applications where Invar is currently utilized, this alloy can be used in fundamental studies that contrast its behavior with its equiatomic counterpart and shed light on the development of HEAs.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Determination of the vacancy distribution over Al cation sites in γ-Al 2 O 3

Although gamma-alumina (γ-Al 2 O 3 ) is an extensively used material with wide-ranging applications due to its inherently high surface area and acidity, its atomic structure is still not fully understood. γ-Al 2 O 3 is described as having a spinel-like structure, where the O sublattice has a face-centered cubic (FCC) arrangement and Al cations are placed in the spinel tetrahedral and octahedral interstitial sites. Achieving the correct stoichiometry of Al 2 O 3 , however, requires the introduction of Al vacancies into some of the interstitial sites. Despite the importance of accurately describing the structure of γ-Al 2 O 3 , the distribution of vacancies between tetrahedral and octahedral sites remains unclear, in part because of the usually poor crystalline quality of γ-Al 2 O 3 that has often been used in previous studies. To determine the actual cation distribution in γ-Al 2 O 3 , single-crystalline γ-Al 2 O 3 was investigated using a correlative approach of experimental and simulated selected-area electron diffraction (SAED) and high-resolution electron energy-loss spectroscopy (EELS). Comparison of the reflection intensities in single-crystal SAED to simulated SAED from models with varied vacancy distributions revealed that vacancies exist primarily on tetrahedral sites, contrary to the placement of vacancies on octahedral sites proposed in several common models. Comparison of EELS spectra—acquired with the highest energy resolution reported so far for γ-Al 2 O 3 —with ab initio multiple scattering EELS simulations confirmed the distribution of vacancies on tetrahedral sites. These results enable more accurate modeling of γ-Al 2 O 3 to better predict its properties in existing and future applications.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Experimental and computational studies on high-entropy carbide MoNbTaVWC 5 under high pressures

High-entropy carbide, MoNbTaVWC 5 , was synthesized from oxide precursors of the constituent metals, mixed with graphite powder in a microwave-generated hydrogen plasma at 26.66 kPa and 2100 °C. Ambient x-ray diffraction analysis confirms the full conversion of oxide precursors into a single-phase, face-centered cubic structure with a lattice parameter a = 4.3309 Å. Nanoindentation measured a hardness of 24.5 ± 1.3 GPa and an elastic modulus of 386 ± 22 GPa. The synthesized sample, mixed with a copper pressure marker, was studied by the radial x-ray diffraction technique with beryllium gasketing in a diamond anvil cell up to 70 GPa. The experimentally measured pressure–volume curve and shear strength were compared with theoretical predictions using the special quasi-random structure technique and density functional theory. MoNbTaVWC 5 achieved a 12% volume compression at 70 GPa and exhibited a high shear strength of 6.6 GPa. The present study demonstrates that the high-entropy carbide MoNbTaVWC 5 exhibits exceptional incompressibility and high strength under extreme conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Chemical Insights into the Formation of Colloidal Iridium Nanoparticles from In Situ X-ray Total Scattering: Influence of Precursors and Cations on the Reaction Pathway

Iridium nanoparticles are important catalysts for several chemical and energy conversion reactions. Studies of iridium nanoparticles have also been a key for the development of kinetic models of nanomaterial formation. However, compared to other metals such as gold or platinum, knowledge on the nature of prenucleation species and structural insights into the resultant nanoparticles are missing, especially for nanoparticles obtained from Ir$_x$Cl$_y$ precursors investigated here. We use in situ X-ray total scattering (TS) experiments with pair distribution function (PDF) analysis to study a simple, surfactant-free synthesis of colloidal iridium nanoparticles. The reaction is performed in methanol at 50 °C with only a base and an iridium salt as precursor. From different precursor salts–IrCl 3 , IrCl 4 , H 2 IrCl 6 , or Na 2 IrCl 6 –colloidal nanoparticles as small as Ir ~55 are obtained as the final product. The nanoparticles do not show the bulk iridium face-centered cubic ($fcc$) structure but show decahedral and icosahedral structures. The formation route is highly dependent on the precursor salt used. Using IrCl 3 or IrCl 4 , metallic iridium nanoparticles form rapidly from Ir$_x$Cl$_y$$^{n–}$ complexes, whereas using H 2 IrCl 6 or Na 2 IrCl 6 , the iridium nanoparticle formation follows a sudden growth after an induction period and the brief appearance of a crystalline phase. With H 2 IrCl 6 , the formation of different Ir$_n$ ($n$ = 55, 55, 85, and 116) nanoparticles depends on the nature of the cation in the base (LiOH, NaOH, KOH, or CsOH, respectively) and larger particles are obtained with larger cations. As the particles grow, the nanoparticle structure changes from partly icosahedral to decahedral. In conclusion, the results show that the synthesis of iridium nanoparticles from Ir$_x$Cl$_y$ is a valuable iridium nanoparticle model system, which can provide new compositional and structural insights into iridium nanoparticle formation and growth.

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

First-principles theory for cerium predicts three distinct face-centered cubic phases

We show results from first-principles calculations for cerium at very high compressions. These reveal a most remarkable behavior in a material; depending on atomic volume, cerium adopts three distinct face-centered cubic (fcc) phases driven by different physical mechanisms. The two well-known a and phases are vigorously debated in the literature, but we focus on the a phase as a metal with delocalized character of the 4f electron. The ultimate high compression fcc phase, here named ω, is driven partly by electrostatics. Our density-functional theory (DFT) study excellently reproduces the experimentally known compression behavior of cerium up to a few Mbar but goes beyond those pressures with structural transitions to tetragonal, hexagonal, and cubic (fcc) phases occurring before 100 Mbar (10000 GPa or 10 TPa). The 4f-electron contribution to the chemical bonding is shown to rule phase transitions and compressibility. The change of 4f occupation nicely explains the pressure dependence of the structural axial ratio in the tetragonal phase. At very high pressure, structures known at low pressures return because of band broadening, electrostatic ion repulsion, and an increase in hybridization between states that under normal conditions can be considered core (atomic like) states and the valence-band states.

36 MATERIALS SCIENCE↗
DOE OSTI2024

High-pressure phase transitions in a laser directed energy deposited Fe-33Cu Alloy

Here, additively manufactured Fe-Cu alloys contain both equilibrium face-centered cubic (FCC) and metastable body-centered cubic (BCC) crystal structure Cu precipitates depending on their size. However, the stability of these nanoscale precipitates under extreme conditions such as high pressures has not been reported. This study investigates the phase transformations and microstructural stability of laser directed energy deposition (DED-LB) made Fe67Cu33 alloy (nominal composition in at.%) under high static pressure deformation using in-situ synchrotron X-ray diffraction under pressure, postmortem high-resolution scanning transmission electron microscopy (HR-STEM), and molecular dynamics (MD) simulations. In-situ XRD results reveal a reversible phase transformation from BCC to the hexagonal close-packed (HCP) structure in the Fe grains at an onset pressure of 16.4 GPa, significantly higher than reported in the literature for pure Fe. Although no high-pressure phase transition was observed in the FCC Cu grains through XRD, HR-STEM analysis uncovers a phase transition to the HCP structure in nanoscale metastable BCC Cu precipitates within the BCC Fe matrix. After decompression, the Fe matrix reverted back to the BCC structure with periodic lath martensite, while regions of the nanoscale BCC Cu precipitates retained the metastable HCP structure. MD simulations support the BCC → HCP transition in the embedded nanoscale coherent Cu precipitates, consistent with the classical Burgers mechanism. Thus, by leveraging in-situ XRD observation, postmortem high-resolution S/TEM microscopy, and MD simulations, this study offers profound insights into the distinctive phase transformations induced by high pressure in DED-LB Fe67Cu33 alloy, which is distinguished by its hierarchical microstructure.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Crystal Nucleation Kinetics and Mechanism: Influence of Interaction Potential

Modulating liquid-to-solid transitions and the resulting crystalline structure for tailored properties is much desired. Colloidal systems are exemplary to this end, but the fundamental knowledge gaps in relating the influence of intermolecular interactions to crystallization behavior continue to hinder progress. In this study, we address this knowledge gap by studying nucleation and growth in systems with modified Lennard-Jones potential. Specifically, we study the commonly used 12-6 potential and a softer 7-6 potential. The thermodynamic state point for the study is chosen such that both systems are investigated at the same level of supercooling and pressure. Under these conditions, we find that the nucleation rate for both systems is comparable. Interestingly, the nucleation pathways and resulting crystal structures are different. In the 12-6 system, nucleation and growth occur predominantly through the FCC structure. Softening the potential alters the critical nucleus composition and introduces two distinct nucleation pathways. One pathway predominantly leads to the nucleus with a body-centered cubic (BCC) structure, while the other favors the face-centered cubic (FCC) arrangement. Our study illustrates that polymorph selection can be achieved through modifications to intermolecular interactions without impacting nucleation kinetics. The results have significant implications in designing approaches for polymorph selection and modulating self-assembly mechanisms.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Nanolamellar phase transition in an additively manufactured eutectic high-entropy alloy under high pressures

Much is unknown about how phase transitions link to micro-/nano-structures in high-entropy systems, especially under extreme pressure and temperature conditions. This work studies the evolution of dual-phase nanolamellar eutectic high-entropy alloy phases of AlCoCrFeNi 2.1 generated by laser powder-bed fusion (L-PBF) for pressures up to 42 GPa. We compare quasi-hydrostatic high pressure synchrotron x-ray diffraction studies on L-PBF printed cylindrical samples up to 5.5 GPa (large-volume Paris–Edinburgh cell) to those carried out on an L-PBF printed foil in a diamond anvil cell where the pressure reached 42 GPa. Our results show that the initially alternating face-centered cubic (FCC) and body-centered cubic (BCC) nanolamellar structure of AlCoCrFeNi 2.1 transformed into single-phase FCC nanolamellae under high pressure with BCC–FCC phase transformation completion at 21 ± 3 GPa. Our results indicate a diffusionless BCC–FCC transformation in this additively manufactured far-from-equilibrium microstructure and demonstrate that the FCC phase is stable up to very high pressures. The measured equation of state for the FCC phase of AlCoCrFeNi 2.1 is presented up to 42 GPa and shows excellent agreement between the data obtained in large-volume press and diamond anvil cell experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Imposing correct jellium response is key to predict the density response by orbital-free DFT

Orbital-free density functional theory constitutes a computationally highly effective tool for modeling electronic structures of systems ranging from room-temperature materials to warm dense matter. Its accuracy critically depends on the employed kinetic energy (KE) density functional, which has to be supplied as an external input. In this work we consider several nonlocal and Laplacian-level KE functionals and use an external harmonic perturbation to compute the static density response at T = 0 K in the linear and beyond-linear response regimes. We test for the satisfaction of exact conditions in the limit of uniform densities and for how approximate KE functionals reproduce the density response of realistic materials (e.g., Al and Si) against the Kohn-Sham DFT reference, which employs the exact KE. The results illustrate that several functionals violate exact conditions in the uniform electron gas (UEG) limit. We find a strong correlation between the accuracy of the KE functionals in the UEG limit and in the strongly inhomogeneous case. This empirically demonstrates the importance of imposing the limit of UEG response for uniform densities and validates the use of the Lindhard function in the formulation of kernels for nonlocal functionals. This conclusion is substantiated by additional calculations for bulk aluminum (Al) with a face-centered cubic (fcc) lattice and silicon (Si) with an fcc lattice, body-centered cubic (bcc) lattice, and semiconducting crystal diamond state. The analysis of fcc Al, and fcc as well as bcc Si data follows closely the conclusions drawn for the UEG, allowing us to extend our conclusions to realistic systems that are subject to density inhomogeneities induced by ions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗