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 1,405 records · Page 78
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

First-Principles Calculations of the Structural, Electronic, Optical, and Mechanical Properties of 21 Pyrophosphate Crystals

Pyrophosphate crystals have a wide array of applications in industrial and biomedical fields. However, fundamental understanding of their electronic structure, optical, and mechanical properties is still scattered and incomplete. In the present research, we report a comprehensive theoretical investigation of 21 pyrophosphates A 2 M (H 2 P 2 O 7 ) 2 •2H 2 O with either triclinic or orthorhombic crystal structure. The molecule H 2 P 2 O 7 is the dominant molecular unit, whereas A = (K, Rb, NH 4 , Tl), M = (Zn, Cu, Mg, Ni, Co, Mn), and H 2 O stand for the cation elements, transition metals, and the water molecules, respectively. The electronic structure, interatomic bonding, partial charge distribution, optical properties, and mechanical properties are investigated by first-principles calculations based on density functional theory (DFT). Most of these 21 crystals are theoretically investigated for the first time. The calculated results show a complex interplay between A, M, H 2 P 2 O 7 , and H 2 O, resulting in either metallic, half-metallic, or semi-conducting characteristics. The novel concept of total bond order density (TBOD) is used as a single quantum mechanical metric to characterize the internal cohesion of these crystals to correlate with the calculated properties, especially the mechanical properties. This work provides a large database for pyrophosphate crystals and a road map for potential applications of a wider variety of phosphates.

36 MATERIALS SCIENCE↗
DOE OSTI2025

High-Temperature Gas Sensor Materials with Properties Predicted via First-Principles Calculations with Machine Learning Modeling and Experimental Corroboration

Understanding the temperature dependence of functional properties of sensing materials is vital for their applications in combustion environments. The electron-phonon coupling that derives the electronic structure change with temperatures is a key property of interest as it affects other sensing responses. Herein, we first assess the temperature dependence of band gap renormalization in sensing materials by employing Allen-Heine-Cardona (AHC) theory with density functional theory (DFT) simulations corroborated with experimental observation. As the AHC calculations are impractical for high-throughput screening of materials, we employ data-driven Gaussian process regression to predict the parameters employed in the O’Donnell empirical model from a set of physical features. To mitigate the reliability issues arising from the small size of the dataset, we apply a Bayesian technique to improve the generalizability of the data-driven models as well as to quantify the uncertainty associated with theoretical predictions. These models capture well the overall trend of the O’Donnell parameters with respect to a reduced feature set obtained by transforming the available physical features. Quantifying the associated uncertainty helps us understand the reliability of the predictions and, therefore, the variation of bandgap as a function of temperature for other novel materials. The predicted candidates from machine learning models are further validated by experiments and DFT calculations.

bandgap renormalization↗
DOE OSTI2022

Stochastic Vector Techniques in Ground-State Electronic Structure

Herein we review a suite of stochastic vector computational approaches for studying the electronic structure of extended condensed matter systems. These techniques help reduce algorithmic complexity, facilitate efficient parallelization, simplify computational tasks, accelerate calculations, and diminish memory requirements. While their scope is vast, we limit our study to ground-state and finite temperature density functional theory (DFT) and second-order many-body perturbation theory. More advanced topics, such as quasiparticle (charge) and optical (neutral) excitations and higher-order processes, are covered elsewhere. We start by explaining how to use stochastic vectors in computations, characterizing the associated statistical errors. Next, we show how to estimate the electron density in DFT and discuss effective techniques to reduce statistical errors. Finally, we review the use of stochastic vectors for calculating correlation energies within the second-order Møller-Plesset perturbation theory and its finite temperature variational form. Example calculation results are presented and used to demonstrate the efficacy of the methods.

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

First-principles calculations of metal surfaces. II. Properties of low-index platinum surfaces toward understanding electron emission

The stability of low-index platinum surfaces and their electronic properties is investigated with density functional theory, toward the goal of understanding the surface structure and electron emission, and identifying precursors to electrical breakdown, on nonideal platinum surfaces. Propensity for electron emission can be related to a local work function, which, in turn, is intimately dependent on the local surface structure. The (1×N) missing row reconstruction of the Pt(110) surface is systematically examined. The (1×3) missing row reconstruction is found to be the lowest in energy, with the (1×2) and (1×4) slightly less stable. In the limit of large (1×N) with wider (111) nanoterraces, the energy accurately approaches the asymptotic limit of the infinite Pt(111) surface. Furthermore, this suggests a local energetic stability of narrow (111) nanoterraces on free Pt surfaces that could be a common structural feature in the complex surface morphologies, leading to work functions consistent with those on thermally grown Pt substrates.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Ultrafast Optically Induced Perturbation of Oxygen Octahedral Rotations in Multiferroic BiFeO 3 Thin Films

The functional properties of complex oxides, including magnetism and ferroelectricity, are closely linked to subtle structural distortions. Ultrafast optical excitations provide the means to manipulate structural features and ultimately to affect the functional properties of complex oxides with picosecond-scale precision. Here, we report that the lattice expansion of multiferroic BiFeO 3 following above-bandgap optical excitation leads to distortion of the oxygen octahedral rotation (OOR) pattern. The continuous coupling between OOR and strain was probed using time-resolved X-ray free-electron laser diffraction with femtosecond time resolution. Density functional theory calculations predict a relationship between the OOR and the elastic strain consistent with the experiment, demonstrating a route to employing this approach in a wider range of systems. Ultrafast control of the functional properties of BiFeO 3 thin films is enabled by this approach because the OOR phenomena are related to ferroelectricity, and via the Fe-O-Fe bond angles, the superexchange interaction between Fe atoms.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Many‐body response of benzene at monolayer MoS 2 : Van der Waals interactions and spectral broadening

Abstract Models of surface enhancement of molecular electronic response properties are challenging for two reasons: (a) molecule‐surface interactions require a simultaneous solution of the molecular and the surface dynamic response (a daunting task), and (b) when solving for the electronic structure of the combined molecule + surface system, it is not trivial to single out the particular physical effects responsible for enhancement. To tackle this problem, in this work, we apply a formally exact decomposition of the system's response function into subsystem contributions by using subsystem density functional theory (DFT), which grants access to dynamic polarizabilities and optical spectra. In order to access information about the interactions between the subsystems, we extend a previously developed subsystem‐based adiabatic connection fluctuation‐dissipation theorem of DFT to separate the additive from the nonadditive correlation energy and identify the nonadditive correlation as the van der Waals interactions. As an example, we choose benzene adsorbed on monolayer MoS 2 . We isolate the contributions to benzene's dynamic response arising from the interaction with the surface, and for the first time, we evaluate the enhancements to the effectiveness of C 6 coefficients as a function of benzene‐MoS 2 distance and adsorption site. We also quantify the spectral broadening of the benzene's electronic excited states due to their interaction with the surface. We find that the broadening has a similar decay law with the molecule‐surface distance as the leading van der Waals interactions (ie, R −6 ) and that the surface enhancement of dispersion interactions between benzene molecules is less than 5% but is still large enough (0.5 kcal/mol) to likely play a role in the prediction of interface morphologies.

Umerbekova, Alina↗
DOE OSTI2020

Density-Matrix Based Extended Lagrangian Born–Oppenheimer Molecular Dynamics

Extended Lagrangian Born–Oppenheimer molecular dynamics [ Phys. Rev. Lett. 2008, 100, 123004] is presented for Hartree–Fock theory, where the extended electronic degrees of freedom are represented by a density matrix, including fractional occupation numbers at elevated electronic temperatures. In contrast to regular direct Born–Oppenheimer molecular dynamics simulations, no iterative self-consistent field optimization is required prior to the force evaluations. To sample regions of the potential energy landscape where the gap is small or vanishing, which leads to particular convergence problems in regular direct Born–Oppenheimer molecular dynamics simulations, an adaptive integration scheme for the extended electronic degrees of freedom is presented. The integration scheme is based on a tunable, low-rank approximation of a fourth-order kernel, K, that determines the metric tensor, T ≡ K T K, used in the extended harmonic oscillator of the Lagrangian that generates the dynamics of the electronic degrees of freedom. Here, the formulation and algorithms provide a general guide to implement extended Lagrangian Born–Oppenheimer molecular dynamics for quantum chemistry, density functional theory, and semiempirical methods using a density matrix formalism.

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

Correlated Anion Disorder in Heteroanionic Cubic TiOF 2

Resolving anion configurations in heteroanionic materials is crucial for understanding and controlling their properties. For anion-disordered oxyfluorides, conventional Bragg diffraction cannot fully resolve the anionic structure, necessitating alternative structure determination methods. We have investigated the anionic structure of anion-disordered cubic (ReO 3 -type) TiOF 2 using X-ray pair distribution function (PDF), 19 F MAS NMR analysis, density functional theory (DFT), cluster expansion modeling, and genetic-algorithm structure prediction. Our computational data predict short-range anion ordering in TiOF 2 , characterized by predominant cis-[O 2 F 4 ] titanium coordination, resulting in correlated anion disorder at longer ranges. To validate our predictions, we generated partially disordered supercells using genetic-algorithm structure prediction and computed simulated X-ray PDF data and 19 F MAS NMR spectra, which we compared directly to experimental data. To construct our simulated 19 F NMR spectra, we derived new transformation functions for mapping calculated magnetic shieldings to predicted magnetic chemical shifts in titanium (oxy)fluorides, obtained by fitting DFT-calculated magnetic shieldings to previously published experimental chemical shift data for TiF 4 . We find good agreement between our simulated and experimental data, which supports our computationally predicted structural model and demonstrates the effectiveness of complementary experimental and computational techniques in resolving anionic structure in anion-disordered oxyfluorides. From additional DFT calculations, we predict that increasing anion disorder makes lithium intercalation more favorable by, on average, up to 2 eV, highlighting the significant effect of variations in short-range order on the intercalation properties of anion-disordered materials.

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

When van der Waals Met Kagome: A 2D Antimonide with a Vanadium-Kagome Network

2D materials showcase unconventional properties emerging from quantum confinement effects. In this work, a “soft chemical” route allows for the deintercalation of K + from the layered antimonide KV 6 Sb 6 , resulting in the discovery of a new metastable 2D-Kagome antimonide K 0.1(1) V 6 Sb 6 with a van der Waals gap of 3.2 Å. The structure of K 0.1(1) V 6 Sb 6 was determined via the synergistic techniques, including X-ray pair distribution function analysis, advanced transmission electron microscopy, and density functional theory calculations. The K 0.1(1) V 6 Sb 6 compound crystallizes in the monoclinic space group C2/m (a = 9.57(2) Å, b = 5.502(8) Å, c = 10.23(2) Å, β = 97.6(2)°, Z = 2). The [V 6 Sb 6 ] layers in K 0.1(1) V 6 Sb 6 are retained upon deintercalation and closely resemble the layers in the parent compound, yet deintercalation results in a relative shift of the adjacent [V 6 Sb 6 ] layers. The magnetic properties of the K 0.1(1) V 6 Sb 6 phase in the 2–300 K range are comparable to those of KV 6 Sb 6 and another Kagome antimonide KV 3 Sb 5 , consistent with nearly temperature-independent paramagnetism. Electronic band structure calculation suggests a nontrivial band topology with flat bands and opening of band crossing afforded by deintercalation. Transport property measurements reveal a metallic nature for K 0.1(1) V 6 Sb 6 and a low thermal conductivity of 0.6 W K –1 m –1 at 300 K. Additionally, ion exchange in KV 6 Sb 6 via a solvothermal route leads to a successful partial exchange of K + with A + (A = Na, Rb, and Cs). Here, this study highlights the tunability of the layered structure of the KV 6 Sb 6 compound, providing a rich playground for the realization of new 2D materials.

2D↗
DOE OSTI2023

Pressure-induced charge orders and their postulated coupling to magnetism in hexagonal multiferroic LuFe 2 O 4

Hexagonal LuFe 2 O 4 is a promising charge order (CO) driven multiferroic material with high charge and spin-ordering temperatures. The coexisting charge and spin orders on Fe 3+ /Fe 2+ sites result in magnetoelectric behaviors, but the coupling mechanism between the charge and spin orders remains elusive. Here, by tuning external pressure, we reveal three charge-ordered phases with suggested correlation to magnetic orders in LuFe 2 O 4 : (i) a centrosymmetric incommensurate three-dimensional CO with ferrimagnetism, (ii) a non-centrosymmetric incommensurate quasi-two-dimensional CO with ferrimagnetism, and (iii) a centrosymmetric commensurate CO with antiferromagnetism. Experimental in situ single-crystal X-ray diffraction and X-ray magnetic circular dichroism measurements combined with density functional theory calculations suggest that the charge density redistribution caused by pressure-induced compression in the frustrated double-layer [Fe 2 O 4 ] cluster is responsible for the correlated spin-charge phase transitions. The pressure-enhanced effective Coulomb interactions among Fe-Fe bonds drive the frustrated (1/3, 1/3) CO to a less frustrated (1/4, 1/4) CO, which induces the ferrimagnetic to antiferromagnetic transition. Our results not only elucidate the coupling mechanism among charge, spin, and lattice degrees of freedom in LuFe 2 O 4 , but also provide a new way to tune the spin-charge orders in a highly controlled manner.

36 MATERIALS SCIENCE↗
DOE OSTI2024

3D Heisenberg universality in the van der Waals antiferromagnet NiPS 3

Van der Waals (vdW) magnetic materials are comprised of layers of atomically thin sheets, making them ideal platforms for studying magnetism at the two-dimensional (2D) limit. These materials are at the center of a host of novel types of experiments, however, there are notably few pathways to directly probe their magnetic structure. We confirm the magnetic order within a single crystal of NiPS 3 and show it can be accessed with resonant elastic X-ray diffraction along the edge of the vdW planes in a carefully grown crystal by detecting structurally forbidden resonant magnetic X-ray scattering. We find the magnetic order parameter has a critical exponent of β ~ 0.36, indicating that the magnetism of these vdW crystals is more adequately characterized by the three-dimensional (3D) Heisenberg universality class. We verify these findings with first-principles density functional theory, Monte-Carlo simulations, and density matrix renormalization group calculations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Emergence of the temperature–density relation in the low-density intergalactic medium

We examine the evolution of the phase diagram of the low-density intergalactic medium during the Epoch of Reionization in simulation boxes with varying reionization histories from the Cosmic Reionization on Computers project. The probability density function (PDF) of gas temperature at fixed density exhibits two clear modes: a warm and a cold temperature mode, corresponding to the gas inside and outside of ionized bubbles. We find that the transition between the two modes is ‘universal’ in the sense that its timing is accurately parametrized by the value of the volume-weighted neutral fraction for any reionization history. This ‘universality’ is more complex than just a reflection of the fact that ionized gas is warm and neutral gas is cold: it holds for the transition at a fixed value of gas density, and gas at different densities transitions from the cold to the warm mode at different values of the neutral fraction, reflecting a non-trivial relationship between the ionization history and the evolving gas density PDF. Furthermore, the ‘emergence’ of the tight temperature–density relation in the warm mode is also approximately ‘universally’ controlled by the volume-weighted neutral fraction for any reionization history. In particular, the ‘emergence’ of the temperature–density relation (as quantified by the rapid decrease in its width) occurs when the neutral fraction is 10 –4 ≲ X HI ≲ 10 –3 for any reionization history. Our results indicate that the neutral fraction is a primary quantity controlling the various properties of the temperature–density relation, regardless of reionization history.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Collisional stopping power of ions in warm dense matter

Here, a model for the collisional stopping of ions on free electrons in warm dense matter is developed and explored. It is based on plasma kinetic theory, but with modifications to address the warm dense matter regime. Specifically, it uses the Boltzmann-Uehling-Uhlenbeck kinetic equation to incorporate effects of Fermi degeneracy of electrons. The cross section is computed from quantum scattering of electrons and ions occurring via the potential of mean force derived from an average atom model, which incorporates effects of strong Coulomb correlations. Predictions from this model show comparable accuracy to results from time-dependent density functional theory calculations for deuterium near solid density and a temperature of several electronvolts, at a fraction of the computational cost. Further, the model captures the transition of a plasma from the classical limit to the degenerate limit, including qualitative behaviors of solid state theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Carrier- and strain-tunable intrinsic magnetism in two-dimensional $MAX_3$ transition metal chalcogenides

Here, we present a density functional theory study of the carrier-density and strain dependence of magnetic order in two-dimensional (2D) $MAX_3$ ($\textit{M}$ = V, Cr, Mn, Fe, Co, Ni; $\textit{A}$ = Si, Ge, Sn; and $\textit{X}$ = S, Se, Te) transition metal trichalcogenides. Our ab initio calculations show that this class of compounds includes wide and narrow gap semiconductors, metals, and half-metals, and that most of these compounds are magnetic. Although antiferromagnetic order is most common, ferromagnetism is predicted in $\textit{M}$SiSe 3 for $\textit{M}$ = Mn and Ni; in $\textit{M}$SiTe 3 for $\textit{M}$ = V and Ni; in MnGeSe 3 ; $\textit{M}$GeTe 3 for $\textit{M}$ = Cr, Mn, and Ni; in FeSnS 3 ; and in MSnTe 3 for $\textit{M}$ = V, Mn, and Fe. Among these compounds CrGeTe 3 , VSnTe 3 , and CrSnTe 3 are ferromagnetic semiconductors. Our calculations suggest that the competition between antiferromagnetic and ferromagnetic order can be substantially altered by strain engineering, and in the semiconductor case also by gating. The associated critical temperatures can be enhanced by means of carrier doping and strains.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Modeling dark photon oscillations in our inhomogeneous Universe

Not provided.A dark photon may kinetically mix with the Standard Model photon, leading to observable cosmological signatures. The mixing is resonantly enhanced when the dark photon mass matches the primordial plasma frequency, which depends sensitively on the underlying spatial distribution of electrons. Crucially, inhomogeneities in this distribution can have a significant impact on the nature of resonant conversions. We develop and describe, for the first time, a general analytic formalism to treat resonant oscillations in the presence of inhomogeneities. Our formalism follows from the theory of level crossings of random fields and only requires knowledge of the one-point probability density function (PDF) of the underlying electron number density fluctuations. We validate our formalism using simulations and illustrate the photon-to-dark photon conversion probability for several different choices of PDFs that are used to characterize the low-redshift Universe.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Coupling of organic cation and inorganic lattice in methylammonium lead halide perovskites: Insights into a pressure-induced isostructural phase transition

Organic-inorganic hybrid metal-halide perovskite (OIHP) materials provide a tunable platform for engineering their optoelectronic properties. Although several high-pressure studies have been conducted from the OIHP family of single crystals and films, the exact nature of the dynamic coupling of the CH 3 NH 3 (MA) cation with the octahedral lattice framework and the mechanisms responsible for the structural phase transformation under pressure are not well captured. By combined photoluminescence (PL), synchrotron-based x-ray diffraction, and Raman-scattering studies as a function of pressure from methylammonium lead bromide (MAPbBr 3 ), we shed light on an isostructural phase transition due to the coupling of the MA cation and the PbBr 6 lattice through hydrogen bonding. The sharp discontinuities at ~ 1 GPa and ~ 3 GPa in the PL peak positions correlate with the structural changes observed in high-pressure XRD and Raman-scattering studies. The electronic band edge as a function of pressure is calculated within density-functional theory. The PL peak position, intensity and width of the excitonic peak show significant changes at 2 GPa, which corroborate the changes observed in high-pressure Raman-scattering studies. The frequencies of the lattice modes and the C-H/N-H bending and stretching modes of the MA cation show anomalous changes and other nuances at 2 GPa. The suppression of rotational and orientational disorder of the organic moiety is initiated at 2 GPa and the ordering is completed by 3.0 GPa, leading to an order-disorder type cubic II (Im$\overline{3}$) to orthorhombic ( Pnma ) phase transition. Along with the revelation of an isostructural transformation at 2 GPa, this paper highlights the impact of molecular vibrations on the electronic properties of MAPbBr 3 under pressure.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Accurate point defect energy levels from non-empirical screened range-separated hybrid functionals: The case of native vacancies in ZnO

We use density functional theory (DFT) with non-empirically tuned screened range-separated hybrid (SRSH) functionals to calculate the electronic properties of native zinc and oxygen vacancy point defects in ZnO, and we predict their defect levels for thermal and optical transitions in excellent agreement with available experiments and prior calculations that use empirical hybrid functionals. Furthermore, the ability of this non-empirical first-principles framework to accurately predict quantities of relevance to both bulk- and defect-level spectroscopy enables high-accuracy DFT calculations with non-empirical hybrid functionals for defect physics, at a reduced computational cost.

Defects↗