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 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Ab Initio Electron-Phonon Interactions in Correlated Electron Systems

Electron-phonon ( e - ph ) interactions are pervasive in condensed matter, governing phenomena such as transport, superconductivity, charge-density waves, polarons, and metal-insulator transitions. First-principles approaches enable accurate calculations of e - ph interactions in a wide range of solids. However, they remain an open challenge in correlated electron systems (CES), where density functional theory often fails to describe the ground state. Therefore reliable e - ph calculations remain out of reach for many transition metal oxides, high-temperature superconductors, Mott insulators, planetary materials, and multiferroics. Here we show first-principles calculations of e - ph interactions in CES, using the framework of Hubbard-corrected density functional theory ( DFT + U ) and its linear response extension ( DFPT + U ), which can describe the electronic structure and lattice dynamics of many CES. We showcase the accuracy of this approach for a prototypical Mott system, CoO, carrying out a detailed investigation of its e - ph interactions and electron spectral functions. While standard DFPT gives unphysically divergent and short-ranged e - ph interactions, DFPT + U is shown to remove the divergences and properly account for the long-range Fröhlich interaction, allowing us to model polaron effects in a Mott insulator. Our work establishes a broadly applicable and affordable approach for quantitative studies of e - ph interactions in CES, a novel theoretical tool to interpret experiments in this broad class of materials.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Local lattice distortions and the structural instabilities in bcc Nb–Ta–Ti–Hf high-entropy alloys: An ab initio computational study

Local lattice distortions (LLD) and structural stability of body-centered cubic (bcc) Nb–Ta–Ti–Hf high-entropy alloys (HEAs) are studied as functions of composition employing ab initio density-functional theory calculations, with specific focus on the role of the relative concentrations of group IV (Ti and Hf) versus group V (Nb and Ta) elements. Calculated results are presented as a function of composition x in Nb x Ta 0.25 Ti (0.75-x)/2 Hf (0.75-x)/2 alloys, for elastic moduli, phonon spectral functions, LLD and structural energy differences for the bcc and competing hexagonal close-packed (hcp) and ω phases. The results highlight the important role of group V elements and LLD in stabilizing the bcc structure. They further reveal how composition x can be tuned to alter both the magnitude of the LLD and structural energy differences. Specifically, the magnitude of the structural energy differences, and elastic and dynamic stability of the bcc phase, are enhanced with increasing x, while the LLD increase in magnitude as this concentration is decreased. The results also show evidence of correlated LLD at lower values of x, reflecting local structural distortions towards the ω phase, but not hcp. The degree of ω-collapse is nevertheless partial i.e., transformation towards this phase is not observed to be complete due to the presence of Ta and Nb. At lower values of x we further find an energy landscape characterized by multiple, nearly degenerate local energy minima for different values of the LLD.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Phase Transformations among TiO2 Polymorphs

Polymorphs widely exist in nature and synthetic systems and are well known to determine material properties. Understanding phase transformation mechanisms among polymorphs enables the design of structures and tuning of phases to tailor material properties. However, current understanding is limited due to the lack of direct observations that track the processes and the structural evolution is not known at atomic scale. Here, integrating (semi) in situ transmission electron microscopy and density functional theory, we report atomic structural evolutions of phase transformation from anatase (A) to rutile (R), brookite (B), R-phase, and TiO. Besides the consistent paths with previous reports, we discover several unreported paths, including a [001] direction and (020) plane of anatase to [100]A and (011)R of rutile, respectively, ( , ) and , . Density functional theory analysis elucidates atomic structural evolution during the processes and over 16% of Ti-O bonds break and reform during the processes with energy barriers of ~0.7-1.0 eV per TiO2 formula unit. Under electron-beam irradiation, anatase particles transform into TiO2-R phase or TiO at high or room temperature, respectively. We also reveal the anisotropic nature of electron-beam effect, which is seldom discussed: dependence of crystallographic orientation with respect to electron-beam irradiation direction. Understanding the atomic structural evolution sheds light on interpreting and controlling TiO2 polymorphs and interface structures for various applications. The revealed electron-beam effects in our work provide guidance for in situ transmission electron microscopy studies.

Song, Miao↗
DOE OSTI2026

Ultrafast One-Dimensional Peierls-Distortion Dynamics in 1⁢T′−Re⁢S 2 Revealed by 4D Electron Microscopy

Rhenium disulfide (ReS 2 ), a prototypical 2D semiconductor with an anisotropic 1⁢T′ structure due to the pronounced Peierls distortion, has demonstrated great potential for polarization-dependent optoelectronics and photonics. Here, we report an ultrafast phase transition occurring within 1 ps, accompanied by a one-dimensional Peierls-distortion relaxation in 1⁢T′−ReS 2 revealed with combined 4D electron microscopy and time-dependent density functional theory calculations. Upon femtosecond laser pulse excitation, the Re-Re dimerization morphology rapidly transforms from a diamond cluster to zigzag chains and persists for several nanoseconds. Here, this ultrafast Peierls-distortion relaxation is further verified by transient changes in optical anisotropy via polarization-dependent transient absorption spectroscopy. Time-dependent density functional theory calculations attribute this novel phase transition to strong correlation between Peierls distortion and impulsive photoexcited carrier doping, predicting a transient band-gap collapse. This Letter opens up an exciting avenue for ultrafast control of Peierls-distortion-induced anisotropy and the metal-insulator transition in 1⁢T′−ReS 2 .

1T’-ReS2↗
DOE OSTI2022

Strongly correlated itinerant magnetism near superconductivity in NiTa 4 Se 8

Metallic ferromagnets with strongly interacting electrons often exhibit remarkable electronic phases such as ferromagnetic superconductivity, complex spin textures, and nontrivial topology. Here, in this report, we discuss the synthesis of a layered magnetic metal NiTa 4 Se 8 (or Ni 1/4 TaSe 2 ) with a Curie temperature of 58 Kelvin. Magnetization data and density functional theory calculations indicate that the nickel atoms host uniaxial ferromagnetic order of about 0.7µ B per atom, while an even smaller moment is generated in the itinerant tantalum conduction electrons. Strong correlations are evident in flat bands near the Fermi level, a high heat capacity coefficient, and a high Kadowaki-Woods ratio. Density functional theory calculations suggest that electron and hole Fermi surfaces in the ferromagnetic phase are associated with opposite spin polarization. When the system is diluted of magnetic ions, the samples become superconducting below about 2 Kelvin. We discuss possible mechanisms for superconductivity in this family.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2020

All-electron ab initio Bethe-Salpeter equation approach to neutral excitations in molecules with numeric atom-centered orbitals

The Bethe-Salpeter equation (BSE) based on GW quasiparticle levels is a successful approach for calculating the optical gaps and spectra of solids and also for predicting the neutral excitations of small molecules. We here present an all-electron implementation of the GW+BSE formalism for molecules, using numeric atom-centered orbital (NAO) basis sets. We present benchmarks for low-lying excitation energies for a set of small organic molecules, denoted in the literature as “Thiel’s set.” Literature reference data based on Gaussian-type orbitals are reproduced to about one millielectron-volt precision for the molecular benchmark set, when using the same GW quasiparticle energies and basis sets as the input to the BSE calculations. For valence correlation consistent NAO basis sets, as well as for standard NAO basis sets for ground state density-functional theory with extended augmentation functions, we demonstrate excellent convergence of the predicted low-lying excitations to the complete basis set limit. A simple and affordable augmented NAO basis set denoted “tier2+aug2” is recommended as a particularly efficient formulation for production calculations. We finally demonstrate that the same convergence properties also apply to linear-response time-dependent density functional theory within the NAO formalism.

Liu, Chi↗
DOE OSTI2022

First-principles wave-vector- and frequency-dependent exchange-correlation kernel for jellium at all densities

Here we propose a spatially and temporally nonlocal exchange correlation (XC) kernel for the spin-unpolarized fluid phase of ground-state jellium for use in time-dependent density functional and linear response calculations. The kernel is constructed to satisfy known properties of the exact XC kernel to accurately describe the correlation energies of bulk jellium and to satisfy frequency-moment sum rules at a wide range of bulk jellium densities, including those low densities that display strong correlation and symmetry breaking. These effects are easier to understand in the simple jellium model than in real systems. All exact constraints satisfied by the recent MCP07 kernel are maintained in the revised MCP07 (rMCP07) kernel, while others are added. The revision $f^{rMCP07}_{XC}$ (q, ω) differs from MCP07 only for nonzero frequencies ω. Only at densities much lower than those of real bulk metals is the frequency dependence of the kernel important for the correlation energy of jellium. As the wave vector q tends to zero, the kernel has a -4$πα(ω)/q^2$ divergence whose frequency-dependent ultranonlocality coefficient $α(ω)$ vanishes in jellium, and is predicted by rMCP07 to be extremely small for the real metals Al and Na.

36 MATERIALS SCIENCE↗
DOE OSTI2026

HfTi_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys hafnium-titanium (Hf-Ti). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Hf and Ti. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions. The calculations have been collected on Air Force HPC11 cluster using the VASP 6.5.1. Additional methodology and file structure information is available in the dataset README.txt file.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Structural Evidence of Interanionic Hydrogen Bonding in Phosphoric Acid Solutions

Interanionic hydrogen bonding (IAHB) is a noncovalent interaction between like-charged ions that challenges conventional electrostatic understanding. This study provides direct structural evidence of IAHB in concentrated aqueous phosphoric acid (PA) solutions, which exhibit >60% dissociation under these conditions. Oxygen K-edge X-ray absorption fine structure spectroscopy, combined with electron affinity time-dependent density functional theory calculations, reveals the formation of stable, cyclic phosphate-phosphate IAHB dimers at PA concentrations ≥7 M. Extended X-ray absorption fine structure data show distinct long-range ordering consistent with these dimers, and near-edge X-ray absorption fine structure spectra confirm a concentration-dependent transition from monomeric to dimeric species. Energy decomposition analysis through density functional theory shows that the formation of solution-phase IAHB is energetically favored and is attributed to polarization of and the charge transfer between the two fragments driven by the surrounding solvent molecules, in addition to the permanent electrostatics. These findings offers crucial structural insights into the H-bonded networks in concentrated PA, highlighting the critical role of solvent in facilitating anion–anion association.

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

Crystal structure of alectinib hydrochloride Type I, C 30 H 35 N 4 O 2 Cl

The crystal structure of alectinib hydrochloride has been solved and refined using synchrotron X-ray powder diffraction data and optimized using density functional theory techniques. Alectinib hydrochloride crystallizes in space groupP2 1 /n(#14) with the following parameters:a= 12.67477(7),b= 10.44076(5),c= 20.38501(12) Å,β= 93.1438(7)°,V= 2693.574(18) Å 3 , andZ= 4 at 295 K. The crystal structure consists of stacks of molecules along theb-axis, and the stacks contain chains of strong N–H⋯Cl hydrogen bonds. One density functional theory calculation moved a proton from an N atom to the Cl, but another calculation yielded a more chemically reasonable result. The powder pattern has been submitted to ICDD for inclusion in the Powder Diffraction File™ (PDF®)

Materials Science↗
DOE OSTI2021

The FMO-DFTB Method

Although the fragment molecular orbital (FMO) method enables electronic structure calculations with near-linear scaling behavior with respect to system size, the computational cost of ab initio methods typically employed in conjunction with FMO is still prohibitive for routine calculations of very large systems or long timescale molecular dynamics simulations. We, therefore, combined the FMO and density-functional tight-binding (DFTB) method, which is one of the emerging semi-empirical quantum chemical methods, and have demonstrated that FMO-DFTB is capable of performing geometry optimizations for systems containing up to one million atoms using limited computational resources. In this chapter, we will review the basics of the DFTB method first before introducing FMO-DFTB, focusing on the relationship with density functional theory and other FMO methodologies. We also demonstrate the latest scalings of FMO-DFTB2 and DFTB3 using three-dimensional water clusters, showing that the most favorable scaling is O(N1.16). Applications of FMO-DFTB to various systems are briefly summarized, and an outlook to future applications is provided.

Nishimoto, Yoshio↗
DOE OSTI2023

The role of electron correlations in the electronic structure of putative Chern magnet TbMn6Sn6

Abstract A member of the RMn 6 Sn 6 rare-earth family materials, TbMn 6 Sn 6 , recently showed experimental signatures of the realization of a quantum-limit Chern magnet. In this work, we use quantum Monte Carlo (QMC) and density functional theory with Hubbard U (DFT + U ) calculations to examine the electronic structure of TbMn 6 Sn 6 . To do so, we optimize accurate, correlation-consistent pseudopotentials for Tb and Sn using coupled-cluster and configuration–interaction (CI) methods. We find that DFT + U and single-reference QMC calculations suffer from the same overestimation of the magnetic moments as meta-GGA and hybrid density functional approximations. Our findings point to the need for improved orbitals/wavefunctions for this class of materials, such as natural orbitals from CI, or for the inclusion of multi-reference effects that capture the static correlations for an accurate prediction of magnetic properties. DFT + U with Mn magnetic moments adjusted to the experiment predict the Dirac crossing in bulk to be close to the Fermi level, within ~120 meV, in agreement with the experiments. Our non-stoichiometric slab calculations show that the Dirac crossing approaches even closer to the Fermi level, suggesting the possible realization of Chern magnetism in this limit.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Mass enhancement in 3 d and s - p perovskites from symmetry breaking

In some d-electron oxides, the measured effective mass m* exptl has long been known to be significantly larger than the model effective mass m* model deduced from mean-field band theory, i.e., m* exptl =βm* model , where β>1 is the “mass-enhancement” or “mass-renormalization” factor. Previous applications of density functional theory (DFT), based on a symmetry-restricted structure with the smallest number of possible magnetic, orbital, and structural degrees of freedom, missed such mass enhancement. This fact has been taken as evidence of strong electronic correlation, often described via the symmetry-restricted dynamic mean-field approach of the many-body theory, being the exclusive enabling physics. This paper uses instead a static density functional approach that does not restrict positional or spin degrees of freedom (symmetry-broken structures). This approach analyzes the contributions of different symmetry-broken modalities to mass enhancement for a few nominally highly correlated d-electron perovskites as well as the nominally uncorrelated, closed-shell s-p bonding perovskites. It shows that the energy-lowering symmetry-broken spin effects (e.g., nonzero local moment in the paramagnetic phase) and structural effects (e.g., atomic displacement) as described in mean-field DFT already manifest mass enhancement for both electrons and holes in a range of d-electron perovskites SrVO 3 , SrTiO 3 , BaTiO 3 , and LaMnO 3 , as well as p-electron perovskites CsPbI 3 and SrBiO 3 , including both metals (SrVO 3 ) and insulators (the rest). This is revealed only when enlarged unit cells of the same parent global symmetry, which are large enough to allow for symmetry-breaking distortions and concomitant variations in spin order, are explored for their ability to lower the total energy. Positional symmetry breaking that leads to mass enhancement includes octahedral rotation in halide perovskites such as CsPbI 3 , Jahn-Teller-like Q${+}\atop{2}$ distortion in LaMnO 3 , and bond disproportionation in SrBiO 3 , while magnetic symmetry breaking resulting in mass enhancement includes the formation of a distribution of local moments in SrVO3 that averages to zero in the paramagnetic phase. Not all symmetry breaking leads to significant mass enhancement, e.g., the rather small octahedral rotations in the nearly perfectly cubic SrTiO 3 cause negligible mass enhancement, as do the paraelectric displacements in BaTiO 3 . In principle, by ergodicity, the two descriptions, i.e., the symmetry-restricted dynamic approach with a single time-fluctuating magnetic moment and the symmetry-broken mean-field approach with a static distribution of spatially fluctuated local moments, are not mutually exclusive but are a choice of representation and consequently, a choice of computational efficiency. In approximate implementations, the symmetry-broken mean-field approach appears to remove much of what was strong correlation in dynamically correlated symmetry-restricted solutions, leaving smaller (“weak”) residual correlation with respect to the exact solution.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Revealing quasi-excitations in the low-density homogeneous electron gas with model exchange–correlation kernels

Time-dependent density functional theory within the linear response regime provides a solid mathematical framework to capture excitations. The accuracy of the theory, however, largely depends on the approximations for the exchange–correlation (xc) kernels. Away from the long-wavelength (or q = 0 short wave-vector) and zero-frequency (ω = 0) limit, the correlation contribution to the kernel becomes more relevant and dominant over exchange. The dielectric function, in principle, can encompass xc effects relevant to describe low-density physics. Furthermore, besides collective plasmon excitations, the dielectric function can reveal collective electron–hole excitations, often dubbed “ghost excitons.” Besides collective excitons, the physics of the low-density regime is rich, as exemplified by a static charge-density wave that was recently found for r s > 69, and was shown to be associated with softening of the plasmon mode. These excitations are seen to be present in much higher density 2D homogeneous electron gases of r s ≳ 4. Here, in this work, we perform a thorough analysis with xc model kernels for excitations of various nature. The uniform electron gas, as a useful model of real metallic systems, is used as a platform for our analysis. We highlight the relevance of exact constraints as we display and explain screening and excitations in the low-density region.

36 MATERIALS SCIENCE↗
DOE OSTI2026

6⁢𝐻-perovskite dimer lattice with antiferromagnetic interactions: Ba 3⁢ Zn 1−𝑥⁢ Ca 𝑥 ⁢Ru 2 ⁢O 9

Here, we investigate the magnetic behavior of the 6⁢𝐻-perovskite dimer lattice Ba 3 ⁢Zn 1−𝑥 ⁢Ca 𝑥 ⁢Ru 2 ⁢O 9 using analytical theory, density functional theory, inelastic neutron scattering, and modeling of historical magnetization and neutron-scattering data. A dimer mean-field theory built upon classical Luttinger–Tisza analysis generates a phase diagram revealing a transition from a nonmagnetic singlet to a finite-moment ground state as interdimer couplings increase. A (generalized) linear spin-wave theory captures multiplet mixing, excitation gap closing, and fluctuation-induced moment suppression. Density-functional-theory calculations on selected compounds, together with neutron spectroscopy of dilute Ba 3 ⁢Zn⁢(Ru 1−𝑥 ⁢Sb 𝑥 ) 2 ⁢O 9 , confirm the exchange hierarchy, enabling quantification of previously published experiments within this framework. Our results identify three mechanisms for magnetic moment suppression—quantum fluctuations, ligand hybridization, and nonmagnetic-singlet/magnetic-multiplet mixing.

Pajerowski, Daniel M. [Oak Ridge National Laborato↗
DOE OSTI2025

Computational Investigation of a CO 2 Conversion Strategy via Diels–Alder Reaction in a Carbon Capture Solvent

Molecular-level insights into reactive separations are crucial for the design of new conversion pathways of carbon dioxide (CO 2 ). This work explores a postulated pathway that directs CO 2 to undergo inverse-electron-demand Diels–Alder reactions to produce heterocycles using the CO 2 chemically fixed on water-lean solvent molecules. Density functional theory calculations are applied to evaluate the lowest unoccupied molecular orbital (LUMO) energies of three types of reactants (1,3-butadiene, 1,3-cyclohexadiene, and 1,2,4,5-tetrazine) with various functional substituents. These calculations also provide a data set (5.8k data) for developing a machine learning model to efficiently predict LUMO energies. A computational screening of LUMO energies for an additional 47k diene and tetrazine candidates is performed, and a list of candidates with lowered LUMO energies by electron-withdrawing substituents is provided. These candidates are further examined by their reaction energy barriers computed from the interatomic potential or density functional theory. Two major energy barriers are identified, one for the proton transfer within the water-lean solvent and the other for the CO 2 transfer from the solvent molecule to the reactant candidate (diene or tetrazine). The functional substituents have a more significant impact on the second barrier but a very slight one on the first barrier. This exploratory work demonstrates a new possibility for guiding experimental efforts toward the chemical conversion of fixated CO 2 to value-added compounds.

Chemical reactions↗
DOE OSTI2026

Structural Propensities in Cs2MBiX6 (M=Na, Ag; X=Cl, Br) Bismuth Halide Double Perovskites

A previously unreported low-temperature phase transition in the bismuth halide double perovskite Cs2AgBiCl6 is reported, thereby establishing trends in the structural ground state across Cs2NaBiCl6, Cs2AgBiCl6, and Cs2AgBiBr6. Using the combined toolkit of variable-temperature synchrotron X-ray and neutron powder diffraction, Raman spectroscopy, and density-functional theory–based electronic structure modeling, we demonstrate a cubic Fm¯3m → tetragonal I4/m transition upon cooling with distinct onset temperatures. Neutron powder diffraction refinements permit the unambiguously assignment of the low-temperature phase of Cs2NaBiCl6 to I4/m, correcting prior reports of an I4/mmm ground state. Cs2AgBiCl6 is also found to transforms to a structure crystallizing in the I4/m space group at low temperatures. Temperaturedependent Raman data and density-functional theory-based modeling capture the softening and freezing of out-of-phase octahedral-tilt modes and quantify relative instabilities. Solid-state nuclear magnetic resonance spectroscopy at room temperature completes the characterization and helps underpin the subtle differences in covalency across the compounds. Trends in the phase transition temperature Ts and tilt magnitudes emerge from coupled effects of halide identity, M(I)–site bonding character, and a mismatch between interatomic distances. These results establish the structure– dynamics–bonding framework for tuning tilt-driven instabilities in halide double perovskites.

Tian, Haowen↗
DOE OSTI2020

Comparison of GW band structure to semiempirical approach for an FeSe monolayer

We present the G 0 W 0 band structure, core levels, and deformation potential of monolayer FeSe in the paramagnetic phase based on a starting mean field of the Kohn-Sham density functional theory (DFT) with the Perdew, Burke, and Ernzerhof functional. We find the GW correction increases the bandwidth of the states forming the M pocket near the Fermi energy, while leaving the Γ pocket roughly unchanged. We then compare the G 0 W 0 quasiparticle band energies with the band structure from a simple empirical + A approach, which was recently proposed to capture the renormalization of the electron-phonon interaction going beyond DFT in FeSe, when used as a starting point in density functional perturbation theory. We show that this empirical correction succeeds in approximating the GW nonlocal and dynamical self-energy in monolayer FeSe and reproduces the GW band structure near the Fermi surface, the core energy levels, and the deformation potential (electron-phonon coupling).

36 MATERIALS SCIENCE↗