Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hybrid density functional calculations”

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

Electronic Structure Calculations for Improving Catalyst Design for Use in the Hybrid Sulfur Cycle

The hybrid sulfur cycle is a promising means for efficient production of large quantities of hydrogen with a small carbon footprint. However, many challenges remain in designing catalysts that are stable at high temperatures and have high sulfur tolerance. Theoretical calculations can provide fundamental insights at the electronic scale to improve our understanding of reaction mechanisms and surface energetics, thereby facilitating improved catalyst design. Here, density functional theory was used to explore the reaction energetics of the hybrid sulfur cycle on Pt, Au, and Pt - Au surfaces. A database of information has been built that contains calculations from which adsorption energetics and density of states information can be extracted.

36 MATERIALS SCIENCE↗
DOE OSTI2022

DFT-FE 1.0: A massively parallel hybrid CPU-GPU density functional theory code using finite-element discretization

In this work, we present DFT-FE 1.0, building on DFT-FE 0.6 [Comput. Phys. Commun. 246, 106853 (2020)], to conduct fast and accurate large-scale density functional theory (DFT) calculations (reaching ~ 100,000 electrons) on both many-core CPU and hybrid CPU-GPU computing architectures. This work involves improvements in the real-space formulation—via an improved treatment of the electrostatic interactions that substantially enhances the computational efficiency—as well high-performance computing aspects, including the GPU acceleration of all the key compute kernels in DFT-FE. We demonstrate the accuracy by comparing the ground-state energies, ionic forces and cell stresses on a wide-range of benchmark systems against those obtained from widely used DFT codes. Further, we demonstrate the numerical efficiency of our implementation, which yields ~ 20× CPU-GPU speed-up by using GPU acceleration on hybrid CPU-GPU nodes. Notably, owing to the parallel-scaling of the GPU implementation, we obtain wall-times of 80–140 seconds for full ground-state calculations, with stringent accuracy, on benchmark systems containing ~ 6, 000 – 15,000 electrons.

pseudopotential↗
DOE OSTI2022

Accurate parameterization of the kinetic energy functional for calculations using exact-exchange

In this report, electronic structure calculations based on Kohn-Sham density functional theory (KSDFT) that incorporate exact exchange or hybrid functionals are associated with large computational expense, a consequence of the inherent cubic scaling bottleneck and large associated prefactor, which limits the length and time scales that can be accessed. Though orbital-free density functional theory (OFDFT) calculations scale linearly with system size and are associated with significantly smaller prefactor, they are limited by the absence of accurate density-dependent kinetic energy functionals. Therefore, the development of accurate density-dependent kinetic energy functionals is important for OFDFT calculations of large realistic systems. To this end, we propose a method to train kinetic energy functional models at the exact-exchange level of theory by using a dictionary of physically relevant terms that have been proposed in the literature in conjunction with linear or nonlinear regression methods to obtain the fitting coefficients. For our dictionary, we use gradient expansion of the kinetic energy, nonlocal models proposed in the literature and their nonlinear combinations, such as a model that incorporates spatial correlations between higher order derivatives of electron density at two points. The predictive capabilities of these models are assessed by using a variety of model one-dimensional systems that exhibit diverse bonding characteristics, such as a chain of eight hydrogens, LiF, LiH, C 4 H 2 , C 4 N 2 and C 3 O 2 . We show that by using data from model one-dimensional KSDFT calculations performed using the exact-exchange functional for only a few neutral structures, it is possible to generate models with high accuracy for charged systems as well as electron and kinetic energy densities during self-consistent field iterations. In addition, we show that it possible to learn both the orbital dependent terms, i.e., the kinetic energy and the exact-exchange energy, and models that incorporate additional nonlinearities in spatial correlations, such as a quadratic model, are needed to capture subtle features of the kinetic energy density that are present in exact-exchange-based KSDFT calculations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Spectroscopic signature of negative electronic compressibility from the Ti core-level of titanium carbonitride MXene

Two-dimensional transition metal carbides, carbonitrides, and nitrides called MXenes exhibit high metallic conductivity, ion intercalation capability and reversible redox activity, prompting their applications in energy storage and conversion, electromagnetic interference (EMI) shielding, and electronics, among many other fields. It has been shown that replacement of about 50% of carbon atoms in the most popular MXene family member, titanium carbide (Ti 3 C 2 T x ), by nitrogen atoms, forming titanium carbonitride (Ti 3 CNT x ), leads to drastically different properties, such as very high negative charge in solution and extreme EMI shielding effectiveness, exceeding all known materials, even metals at comparable thicknesses. Here, by using ultraviolet photoemission spectroscopy (UPS), the electronic structures of Ti 3 CNT x and Ti 3 C 2 T x are systematically investigated and compared as a function of charge carrier density. We observe that, in contrast to Ti 3 C 2 T x , the Ti 3p core-level of Ti 3 CNT x exhibits a counterintuitive shift to a lower binding energy of up to approximately 250 meV upon increasing the electron density, which is a spectroscopic signature of negative electronic compressibility (NEC). These experimentally measured chemical potential shifts are well-captured by the density functional theory (DFT) calculation. The DFT results also further suggest that the hybridization of titanium-nitrogen bonding in Ti 3 CNT x helps promoting the available states of Ti atoms for receiving more electron above the Fermi level and leads to the observed NEC. Furthermore, our findings explain the differences in electronic properties between the two very important and widely studied MXenes and also suggest a new strategy to apply the NEC effect of Ti 3 CNT x in energy and charge storage applications.

2D materials↗
DOE OSTI2020

Passivating Nucleobases Bring Charge Transfer Character to Optically Active Transitions in Small Silver Nanoclusters

DNA-wrapped silver nanoclusters (DNA–AgNCs) are known for their efficient luminescence. However, their emission is highly sensitive to the DNA sequence, the cluster size, and its charge state. To get better insights into photophysics of these hybrid systems, simulations based on density functional theory (DFT) are performed. Our calculations elucidate the effect of the structural conformations, charges, solvent polarity, and passivating bases on optical spectra of DNA–AgNCs containing five and six Ag atoms. It is found that inclusion of water in calculations as a polar solvent media results in stabilization of nonplanar conformations of base-passivated clusters, while their planar conformations are more stable in vacuum, similar to the bare Ag5 and Ag6 clusters. Cytosines and guanines interact with the cluster twice stronger than thymines, due to their larger dipole moments. In addition to the base–cluster interactions, hydrogen bonds between bases notably contribute to the structure stabilization. While the relative intensity, line width, and the energy of absorption peaks are slightly changing depending on the cluster charge, conformations, and base types, the overall spectral shape with five well-resolved bands at 2.5–5.5 eV is consistent for all structures. Independent of the passivating bases and the cluster size and charge, the low energy optical transitions at 2.5–3.5 eV exhibit a metal to ligand charge transfer (MLCT) character with the main contribution emerging from Ag-core to the bases. Cytosines facilitate the MLCT character to a larger degree comparing to the other bases. However, the doublet transitions in clusters with the open shell electronic structure (Ag5 and Ag6 + ) result in appearance of additional red-shifted (<2.5 eV) and optically weak band with negligible MLCT character. The passivated clusters with the closed shell electronic structure (Ag5 + and Ag6) exhibit higher optical intensity of their lowest transitions with much higher MLCT contribution, thus having better potential for emission, than their open shell counterparts.

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

Phase Switching as the Origin of Large Piezoelectric Response in Organic-Inorganic Perovskites: A First-Principles Study

Piezoelectrics are critical functional components of many practical applications such as sensors, ultrasonic transducers, actuators, medical imaging, telecommunications. So far, the best performing piezoelectrics are ferroelectric ceramics, many of which are toxic, heavy, hard, and cost-ineffective. Recently, a groundbreaking discovery of extraordinarily large piezoelectric coefficients in the family of organic-inorganic perovskites gave a hope for a cheaper, environmentally friendly, inexpensive, light-weight and flexible alternative. However, the origin of such response in organic-inorganic ferroelectrics whose spontaneous polarization is an order of magnitude smaller than for inorganic counterparts remains unclear. In this study we employ first-principles simulations to predict that the mechanism associated with large piezoelectric constants is of extrinsic origin and associated with switching between the stable phase and previously overlooked energetically competitive metastable phase that can be stabilized by the external stress. Here, the phase switching changes the polarization direction and, therefore, gives origin to large piezoelectric response, similar to PbZr$_{1-x}$Ti$_x$O$_3$ near morphotropic phase boundary. Existence of such metastable phases is likely to manifest as the dynamical molecular disorder above the Curie temperature and, therefore, could be intrinsic to the entire family of organic-inorganic ferroelectrics with such disorder.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Molecular Fe-N 4 Moieties Coupled with Atomic Co-N 4 Sites Toward Improved Oxygen Reduction Performance

Research on high-efficiency and cost-efficient catalysts for oxygen reduction reaction (ORR) is still a vital but challenging issue for commercializing metal–air batteries. Herein, a single-molecule/atom hybrid catalyst is developed to boost the ORR, in which iron phthalocyanine molecules containing molecular Fe-N 4 moieties couple with atomic Co-N 4 sites on the surface of polyhedral carbon. Density functional theory calculations reveal that face-to-face laminated construction of Fe-N 4 and Co-N 4 in the hybrid catalyst can effectively modulate the electronic structure of active iron atoms and reduce the energy barrier of the rate-determining step for ORR. As a result, this hybrid catalyst demonstrates excellent ORR performance, featuring a half-wave potential of 0.904 V, a peak power density of 238.3 mW cm -2 for zinc–air battery, and outstanding electrocatalytic stability. Here, this work offers a distinctive and robust molecular/atomic engineering approach to creating efficient electrocatalysts, advancing the fields of metal–air batteries.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Band-edge Engineering to Eliminate Radiation-Induced Defect States in Perovskite Scintillators

Under radiative environments such as extended hard X- or γ-rays, degradation of scintillation performance is often due to irradiation-induced defects. To overcome the effect of deleterious defects, novel design mitigation strategies are needed to identify and design more resilient materials. The potential for band-edge engineering to eliminate the effect of radiation-induced defect states in rare-earth-doped perovskite scintillators is explored, taking Ce 3+ -doped LuAlO 3 as a model material system, using density functional theory (DFT)-based DFT + U and hybrid Heyd–Scuseria–Ernzerhof (HSE) calculations. Furthermore, from spin-polarized hybrid HSE calculations, the Ce 3+ activator ground-state 4f position is determined to be 2.81 eV above the valence band maximum in LuAlO 3 . Except for the oxygen vacancies which have a deep level inside the band gap, all other radiation-induced defects in LuAlO 3 have shallow defect states or are outside the band gap, that is, relatively far away from either the 5d 1 or the 4f Ce 3+ levels. Finally, we examine the role of Ga doping at the Al site and found that LuGaO 3 has a band gap that is more than 2 eV smaller than that of LuAlO 3 . Specifically, the lowered conduction band edge envelopes the defect gap states, eliminating their potential impact on scintillation performance and providing direct theoretical evidence for how band-edge engineering could be applied to rare-earth-doped perovskite scintillators.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Conformal high-entropy oxide coatings enable fast and durable surface oxygen reactions

Developing active and durable air electrodes for efficient oxygen reactions is challenging for protonic ceramic cells (PCCs), especially at temperatures below 550°C. Here, in this study, we report a rationally designed conformal coating with a high-entropy PrNi 0.2 Mn 0.2 Co 0.2 Fe 0.2 Cu 0.2 O 3−δ (PNMCFC) perovskite structure on the surface of a state-of-the-art PrBaCo 2 O 5+δ (PBC) air electrode. The formed hybrid air electrode (PNMCFC-PBC) shows faster surface oxygen kinetics and a more stable phase structure in high-humidity air than the bare PBC electrode. Further density functional theory calculations suggest that the conformal coating mitigates Ba segregation at the interface and improves oxygen-related reactions, enhancing overall stability and electrocatalytic performance. The cells with the developed hybrid electrodes show encouraging electrochemical performance at 550°C: a polarization resistance of 0.72 Ω cm 2 , a peak power density of 1.30 W cm −2 , an electrolysis current density of −1.36 A cm −2 at 1.3 V, and reasonable operating stabilities (∼200 h at 550°C).

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2021

Octahedral rotations in Ruddlesden-Popper layered oxides under pressure from first principles

The combination of reduced dimensionality and tunable structural distortions in layered perovskite oxides makes these materials ideal platforms for designing novel properties and functionalities. One example is hybrid improper ferroelectricity in n = 2 Ruddlesden-Popper oxides, where the combination of a layered crystal structure and rotations of the metal-oxide octahedra break symmetry and induce a polarization. Precisely controlling the octahedral rotation distortions, for example by the application of hydrostatic pressure, provides a pathway to tune and optimize the properties of these materials. Here we combine group theoretic methods, density functional theory calculations, and Landau theory analysis to investigate how octahedral rotations respond to pressure in the hybrid improper ferroelectrics Sr 3 Zr 2 O 7 , Ca 3 Ti 2 O 7 , and Sr 3 Sn 2 O 7 . We find that factors that are known to control the pressure response of ABO 3 perovskites—the formal charge of the A- and B-site cations, tolerance factor, and B-site chemistry—also impact the pressure response of these layered perovskites. We also show that coupling between the octahedral rotation and strain order parameters plays a key role in determining the overall pressure response. Despite some similarities, we find that these layered perovskites display a distinct pressure response compared to their ABO 3 perovskite analogs. By identifying trends and underlying mechanisms that control octahedral rotations in Ruddlesden-Popper oxides under pressure, this work lays the foundation for tailoring the structure and properties of these materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Charged point defect benchmark of Hybrid and GGA-PBE

Data used by the publication "High-throughput calculations of charged point defect properties with semi-local density functional theory - performance benchmarks for materials screening applications." This work presented an in-depth benchmark analysis of automated, semi-local point defect calculations with a-posteriori corrections, compared to 245 “gold standard” hybrid calculations previously published. We considered three different a-posteriori correction sets for semi-local calculations, implemented in a fully automated workflow, and consider the qualitative and quantitative differences for four different categories of defect information: thermodynamic transition levels, formation energies, fermi levels, and dopability limits. We highlighted the type of qualitative information about point defect properties that can be extracted from high-throughput calculations based on semi-local DFT methods, while also demonstrating the limits of quantitative accuracy that can be achieved by these approaches.

Broberg, Danny↗
DOE OSTI2024

Uncovering an Interfacial Band Resulting from Orbital Hybridization in Nickelate Heterostructures

The interaction of atomic orbitals at the interface of perovskite oxide heterostructures has been investigated for its profound impact on the band structures and electronic properties, giving rise to unique electronic states and a variety of tunable functionalities. In this study, we conducted an extensive investigation of the optical and electronic properties of epitaxial NdNiO 3 thin films grown on a series of single crystal substrates. Unlike films synthesized on other substrates, NdNiO 3 on SrTiO 3 (NNO/STO) gives rise to a unique band structure which features an additional unoccupied band situated above the Fermi level. Our comprehensive investigation, which incorporated a wide array of experimental techniques and density functional theory calculations, revealed that the emergence of the interfacial band structure is primarily driven by the orbital hybridization between Ti 3d orbitals of the STO substrate and O 2p orbitals of the NNO thin film. Furthermore, exciton peaks have been detected in the optical spectra of the NNO/STO film, attributable to the pronounced electron-electron (e-e) and electron-hole (e-h) interactions propagating from the STO substrate into the NNO film. Further, these findings underscore the substantial influence of interfacial orbital hybridization on the electronic structure of oxide thin-films, thereby offering key insights into tuning their interfacial properties.

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

Stability, electronic quantum states, and magnetic interactions of Er 3+ ions in Ga 2 ⁢O 3

Here, we report an ab initio study of phase stability, defect formation, electronic structure, and multiple magnetic, Dzyaloshinskii-Moriya, optical, hyperfine, and crystal field interactions in erbium (Er)-doped wide band gap 𝛼- and 𝛽-gallium oxides (Ga 2 ⁢O 3 ), critically important to make a foundation for both optoelectronic and quantum information applications. The chemical, structural, mechanical, and dynamical stabilities of the pristine phases are confirmed from respective negative formation energies, negative cohesive energies, favorable elastic constants, and positive phonon frequencies. The phonon dispersions indicate that the Ga-O bonds are uniform in the 𝛼-phase, while they vary in the 𝛽-phase due to the anisotropic polyhedral movement. The defect formation energy analysis confirms that both Er-doped 𝛼- and 𝛽−Ga 2 ⁢O 3 prefer Er 3+ (neutral) state. The underestimated band gaps of the pristine phases from standard density functional theory (DFT) calculations as compared to experimental values are corrected by employing the hybrid functional calculations, resulting in the indirect band gaps of 5.21 eV in 𝛼−Ga 2 ⁢O 3 and 4.94 eV in 𝛽−Ga 2 ⁢O 3 . The site preference energy analysis indicates partial occupation of Er in the octahedral site of Ga. The anisotropic nature of hyperfine tensor coefficients of Er is similar in both phases, which may be due to the occupation of Er in the same octahedral Ga site. On the other hand, the calculated magnetic exchange interaction between two Er dopants is negative for 𝛼 and positive for 𝛽, indicating an antiferromagnetic (AFM) ground state in the former and a ferromagnetic (FM) ground state in the latter. Large values of Dzyaloshinskii-Moriya interactions (DMIs) are obtained along the 𝑥 direction in the 𝛼 and along the 𝑦 direction in the 𝛽. The large DMI may support exotic magnetic textures, a promising direction for spintronic applications. The analysis of dielectric constants and refractive indices of both pristine and Er-doped phases shows a good agreement with available experimental values. The calculated optical anisotropy is slightly higher in 𝛽 than those in 𝛼, which is due to the involvement of lower symmetry in 𝛽. The crystal field coefficients (CFCs) calculated from DFT are used to analyze 4⁢𝑓 multiplets and 4⁢𝑓 −4⁢𝑓 transitions. Thus calculated lowest energy level of the first excited state to the lowest energy level of the ground state is about 1.53 µ⁢m, which is in a good agreement with available experiments, and it falls within the quantum telecommunication wavelength range.

3-dimensional systems↗
DOE OSTI2021

Kohn-Sham Solver (KSSOLV) v2.0

KSSOLV is a MATLAB toolbox for solving Kohn-Sham density functional theory based electronic structure eigenvalue problems. It uses an object oriented features of MATLAB to represent atom, molecules, wavefunctions and Hamiltonians and their operations. It is designed to make it easier for users to prototype and test new algorithms for solving the Kohn-Sham problem. KSSOLV2.0 contains significant improvement over the original KSSOLV described in a paper published in ACM Transaction on Mathematical Software (attached). In addition to performing ground state calculation for small molecules, it can also perform geometry optimization for both molecules and solids. It uses standard pseudopotentials and implements local density approximation, generalized gradient approximation and hybrid functionals. Future releases will also include time-dependent DFT and post DFT calculations such as the GW quasi-particle energy calculation and Bethe-Salpeter equation solver for optical absorption.

Yang, Chao↗
DOE OSTI2021

Beyond the Anderson rule: importance of interfacial dipole and hybridization in van der Waals heterostructures

Vertical stacking of two-dimensional materials with weak van der Waals (vdW) interactions has laid the ground for breakthroughs in physics as well as in technological applications. Although vdW interactions dominate interlayer binding, interlayer electronic coupling may not be negligible and can lead to properties beyond the superposition of constituent monolayers. Here, studying heterobilayers of transition-metal dichalcogenides (MQ 2 ; M = Mo, Ni, Pt; Q = S, Se) by means of density functional theory calculations, we show two mechanisms that influence the band gaps of vdW heterostructures beyond the Anderson rule: (1) interfacial hybridization (mainly involving out-of-plane states, such as chalcogen p z -states), which leads to an upshift in the valence band maxima and accordingly a decrease in the band gap. (2) Formation of an interfacial electric dipole, resulting in an effective gap increase in type-II junctions. While the former is material specific, depending on the proximity of p z -states to each other and the valence band maxima, the latter can be generally described using a model based on the charge density decay outside the monolayers and the pristine band edge positions with respect to the vacuum level, irrespective of junction type.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Emergent and robust ferromagnetic-insulating state in highly strained ferroelastic LaCoO 3 thin films

Transition metal oxides are promising candidates for the next generation of spintronic devices due to their fascinating properties that can be effectively engineered by strain, defects, and microstructure. An excellent example can be found in ferroelastic LaCoO 3 with paramagnetism in bulk. In contrast, unexpected ferromagnetism is observed in tensile-strained LaCoO 3 films, however, its origin remains controversial. Here we simultaneously reveal the formation of ordered oxygen vacancies and previously unreported long-range suppression of CoO 6 octahedral rotations throughout LaCoO 3 films. Supported by density functional theory calculations, we find that the strong modification of Co 3d-O 2p hybridization associated with the increase of both Co-O-Co bond angle and Co-O bond length weakens the crystal-field splitting and facilitates an ordered high-spin state of Co ions, inducing an emergent ferromagnetic-insulating state. Our work provides unique insights into underlying mechanisms driving the ferromagnetic-insulating state in tensile-strained ferroelastic LaCoO 3 films while suggesting potential applications toward low-power spintronic devices.

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↗