Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “functioning”

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,495 records · Page 83
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Investigating the vertical distribution of the disk as a function of radial action

As heating processes can broaden the distributions of radial actions and the vertical distributions of the Galactic disks, we investigate the vertical distribution of the Galactic disks as a function of radial action based on Apache Point Observatory Galactic Evolution Experiment (APOGEE) and Gaia data in order to deepen our understanding of the formation and heating history of the Galactic disks. The vertical distributions of the thin and thick Galactic disks defined in the chemical plane were fitted with a simple exponential function with a free parameter of scale height in different radial action ranges. Therefore, we were able to analyze the scale height as a function of radial action for different disk populations. We find that the distributions of radial action for both the thin and thick disks can be approximately described by pseudo-isothermal distributions, which give a statistical measurement for the temperature of the disk as indicated by the mean radial action of the star sample. Estimations of the scale heights in different radial action ranges for these pseudo-isothermal distributions of the disks seem to show fixed relationships between radial action J R and scale height h. We describe these relationships with a two-parameter function of h = √(J R /a)+b, where a and b are free parameters. When testing with a three-parameter function of h = α √(J R /a)+b, we find that this two-parameter function describes the thin disk well, but we note the function should be used with care for the thick disk. When comparing the best-fit relationships between the inner and outer disk for both of the thin and thick disks, we find that the relationships are nearly the same for the thin disks but are different for the thick disks. The inner thick disk shows a nearly flattened relationship, while the outer thick disk presents a gradually increasing relationship. This work highlights an alternative way to unveil the heating history of the disks by investigating the relationship between scale height and radial action, as these relationships encode the formation and heating history of the Galactic disks.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Self-consistent implementation of meta-GGA functionals for the ONETEP linear-scaling electronic structure package

Accurate and computationally efficient exchange-correlation functionals are critical to the successful application of linear-scaling density functional theory (DFT). Local and semi-local functionals of the density are naturally compatible with linear-scaling approaches, having a general form which assumes the locality of electronic interactions and which can be efficiently evaluated by numerical quadrature. Presently, the most sophisticated and flexible semi-local functionals are members of the meta-generalized-gradient approximation (meta-GGA) family, and depend upon the kinetic energy density, τ, in addition to the charge density and its gradient. In order to extend the theoretical and computational advantages of τ-dependent meta-GGA functionals to large-scale DFT calculations on thousands of atoms, we have implemented support for τ-dependent meta-GGA functionals in the ONETEP program. In this paper we lay out the theoretical innovations necessary to implement τ-dependent meta-GGA functionals within ONETEP's linear-scaling formalism. We present expressions for the gradient of the τ-dependent exchange-correlation energy, necessary for direct energy minimization. We also derive the forms of the τ-dependent exchange-correlation potential and kinetic energy density in terms of the strictly localized, self-consistently optimized orbitals used by ONETEP. To validate the numerical accuracy of our self-consistent meta-GGA implementation, we performed calculations using the B97M-V and PKZB meta-GGAs on a variety of small molecules. Using only a minimal basis set of self-consistently optimized local orbitals, we obtain energies in excellent agreement with large basis set calculations performed using other codes. Lastly, to establish the linear-scaling computational cost and applicability of our approach to large-scale calculations, we present the outcome of self-consistent meta-GGA calculations on amyloid fibrils of increasing size, up to tens of thousands of atoms.

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

Orientation dependence of the work function for metal nanocrystals

We present that work function values measured at different surfaces of a metal are usually different. This raises an interesting question: What is the work function of a nano-size crystal, where differently oriented facets can be adjacent? Work functions of metallic nanocrystals are also of significant practical interest, especially in catalytic applications. Using real space pseudopotentials constructed within density functional theory, we compute the local work function of large aluminum and gold nanocrystals. Lastly, we investigate how the local work function follows the change of the surface plane orientation around multifaceted nanocrystals, and we establish the importance of the orbital character near the Fermi level in determining work function differences between facets.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Survival of the most transferable at the top of Jacob’s ladder: Defining and testing the ω B97M(2) double hybrid density functional

A meta-generalized gradient approximation, range-separated double hybrid (DH) density functional with VV10 non-local correlation is presented. The final 14-parameter functional form is determined by screening trillions of candidate fits through a combination of best subset selection, forward stepwise selection, and random sample consensus (RANSAC) outlier detection. The MGCDB84 database of 4986 data points is employed in this work, containing a training set of 870 data points, a validation set of 2964 data points, and a test set of 1152 data points. Following an xDH approach, orbitals from the ωB97M-V density functional are used to compute the second-order perturbation theory correction. The resu|ting functional, ωB97M(2) is benchmarked against a variety of leading double hybrid density functionals, including B2PLYP-D3(BJ), B2GPPLYP-D3(BJ), ωB97X-2(TQZ), XYG3, PTPSS-D3(0), XYGJ-OS, DSD-PBEP86-D3(BJ), and DSD-PBEPBE-D3(BJ). Encouragingly, the overall performance of ωB97M(2) on nearly 5000 datq points clearly surpasses that of all of the tested density funstionals. As a Rung 5 density functional, ωB97M(2) completes our family of combinatorially optimized functionals, complementing B97M-V on Rung 3, and ωB97X-V and ωB97M-V on Rung 4. The results suggest that ωB97M(2) has the potential to serve as a powerful predictive tool for accurate and efficient electronic structure calculations of main-group chemistry.

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

Effects of the guide field on electron distribution functions in the diffusion region of asymmetric reconnection

Crescent electron distribution functions in the electron diffusion region in asymmetric reconnection with a guide field (strength less than 50% of the upstream magnetic field) are investigated by theory and fully kinetic simulations. Electron motion in the electron diffusion region is a combination of meandering across the current sheet and gyration around the guide field. Combining the two canonical momentum conservation laws in a one-dimensional model with the energy conservation law, we derive the shape of crescent electron distribution functions in both velocity planes: perpendicular to the magnetic field and coplanner with the magnetic field. There are two major effects of the guide field on crescent distribution functions: (1) widening the opening angle of the crescent in a reduced distribution function and (2) the linear cut-off of a distribution function nonparallel to the magnetic field. We also discuss a translational mapping (perpendicular to the linear cut-off) of distribution functions between the X-line and another position. In addition, considering the gyration due to the normal magnetic field, we predict the shape of distribution functions in general positions in the electron diffusion region. Finally, the predictions will be useful to identify features of magnetic reconnection in a guide field by satellite observations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Shadow energy functionals and potentials in Born–Oppenheimer molecular dynamics

In Born–Oppenheimer molecular dynamics (BOMD) simulations based on the density functional theory (DFT), the potential energy and the interatomic forces are calculated from an electronic ground state density that is determined by an iterative self-consistent field optimization procedure, which, in practice, never is fully converged. The calculated energies and forces are, therefore, only approximate, which may lead to an unphysical energy drift and instabilities. Here, we discuss an alternative shadow BOMD approach that is based on backward error analysis. Instead of calculating approximate solutions for an underlying exact regular Born–Oppenheimer potential, we do the opposite. Instead, we calculate the exact electron density, energies, and forces, but for an underlying approximate shadow Born–Oppenheimer potential energy surface. In this way, the calculated forces are conservative with respect to the approximate shadow potential and generate accurate molecular trajectories with long-term energy stabilities. We show how such shadow Born–Oppenheimer potentials can be constructed at different levels of accuracy as a function of the integration time step, δt, from the constrained minimization of a sequence of systematically improvable, but approximate, shadow energy density functionals. For each energy functional, there is a corresponding ground state Born–Oppenheimer potential. These pairs of shadow energy functionals and potentials are higher-level generalizations of the original “zeroth-level” shadow energy functionals and potentials used in extended Lagrangian BOMD. The proposed shadow energy functionals and potentials are useful only within this extended dynamical framework, where also the electronic degrees of freedom are propagated as dynamical field variables together with the atomic positions and velocities. The theory is quite general and can be applied to MD simulations using approximate DFT, Hartree–Fock, or semi-empirical methods, as well as to coarse-grained flexible charge models.

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

Exploring the exact limits of the real-time equation-of-motion coupled cluster cumulant Green’s functions

In this paper, we analyze the properties of the recently proposed real-time equation-of-motion coupled-cluster (RT-EOM-CC) cumulant Green’s function approach [Rehr et al., J. Chem. Phys. 152, 174113 (2020)]. We specifically focus on identifying the limitations of the original time-dependent coupled cluster (TDCC) ansatz and propose an enhanced double TDCC ansatz, ensuring the exactness in the expansion limit. In addition, we introduce a practical cluster-analysis-based approach for characterizing the peaks in the computed spectral function from the RT-EOM-CC cumulant Green’s function approach, which is particularly useful for the assignments of satellite peaks when many-body effects dominate the spectra. Our preliminary numerical tests focus on reproducing, approximating, and characterizing the exact impurity Green’s function of the three-site and four-site single impurity Anderson models using the RT-EOM-CC cumulant Green’s function approach. The numerical tests allow us to have a direct comparison between the RT-EOM-CC cumulant Green’s function approach and other Green’s function approaches in the numerical exact limit.

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

Chemical synthesis rewriting of a bacterial genome to achieve design flexibility and biological functionality

Understanding how to program biological functions into artificial DNA sequences remains a key challenge in synthetic genomics. Here, we report the chemical synthesis and testing of Caulobacter ethensis-2.0 ( C. eth-2.0 ), a rewritten bacterial genome composed of the most fundamental functions of a bacterial cell. We rebuilt the essential genome of Caulobacter crescentus through the process of chemical synthesis rewriting and studied the genetic information content at the level of its essential genes. Within the 785,701-bp genome, we used sequence rewriting to reduce the number of encoded genetic features from 6,290 to 799. Overall, we introduced 133,313 base substitutions, resulting in the rewriting of 123,562 codons. We tested the biological functionality of the genome design in C. crescentus by transposon mutagenesis. Our analysis revealed that 432 essential genes of C. eth-2.0 , corresponding to 81.5% of the design, are equal in functionality to natural genes. These findings suggest that neither changing mRNA structure nor changing the codon context have significant influence on biological functionality of synthetic genomes. Discovery of 98 genes that lost their function identified essential genes with incorrect annotation, including a limited set of 27 genes where we uncovered noncoding control features embedded within protein-coding sequences. In sum, our results highlight the promise of chemical synthesis rewriting to decode fundamental genome functions and its utility toward the design of improved organisms for industrial purposes and health benefits.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

The Impulse Approximation Scattering Function and Its Use in Monte Carlo Photon Transport Simulations

The viability of using the impulse approximation scattering function in Monte Carlo photon transport simulations is explored. This scattering function can be constructed from the double differential incoherent scattering cross section developed by Ribberfors and Berggren. A commonly used method for modeling photon Doppler broadening, which is referred to as the hybrid Doppler broadening method, can also be derived from this cross section. A new photon Doppler broadening method, called the consistent Doppler broadening method, is derived and discussed. This method eliminates some of the commonly employed approximations in the hybrid Doppler broadening method, in part, by using the impulse approximation scattering function. Integrated incoherent cross sections generated using the impulse approximation scattering function and the widely used Waller-Hartree scattering function are in good agreement above 20 keV. Below 20 keV, differences as high as 70% are observed, which differs from the roughly 5% differences observed by Ribberfors for some of the materials. Integral and spectral quantities for two problems are also generated using the Monte Carlo photon transport capabilities of the Framework for REsearch in Nuclear ScIence and Engineering (FRENSIE). Due to the small, relative result differences observed when using the impulse approximation scattering function, it is considered a viable alternative to the Waller-Hartree scattering function. In addition, some small, but expected, differences in spectral fluxes at low energies can be avoided by adopting the consistent Doppler broadening method.

42 ENGINEERING↗
DOE OSTI2011

Partial pair correlation functions and viscosity of liquid Al-Si hypoeutectic alloys via high-energy X-ray diffraction experiments

The liquid structure of Al–Si hypoeutectic binary alloys was characterized by diffraction experiments using a high-energy X-ray (synchrotron) beam source. The diffraction experiments were carried out for liquid pure Al, Al–3 wt% Si, Al–7 wt% Si, Al–10 wt% Si and Al–12.5 wt% Si alloys at several temperatures. The salient structure information such as structure factor (SF), pair distribution function (PDF), radial distribution function (RDF), coordination number (CN) and atomic packing densities (PD) were quantified as a function of Si concentration and melt temperatures. Reverse Monte Carlo (RMC) analysis was carried out using the diffraction experimental data to quantify the partial pair correlation functions, such as partial structure factor, partial pair distribution function (PPDF) and partial radial distribution function. Furthermore, the partial pair distribution function and the liquid atomic structure information were used in a semi-empirical model to evaluate the viscosity of these liquid alloys at various melt temperatures. The results show that the viscosity determined by semi-empirical methods using the atomic structure information is in good agreement with the experimentally determined viscosity values.

Srirangam, Prakash↗
DOE OSTI2011

COSMOLOGICAL EVOLUTION OF SUPERMASSIVE BLACK HOLES. I. MASS FUNCTION AT 0 < z {approx}< 2

We present the mass function of supermassive black holes (SMBHs) over the redshift range z = 0-2, using the latest deep luminosity and mass functions of field galaxies to constrain the masses of their spheroids, which we relate to SMBH mass through the empirical correlation between SMBH and spheroid mass (the M{sub .}-M{sub sph} relation). In addition to luminosity fading of the stellar content of the spheroids, we carefully consider the variation of the bulge-to-total luminosity ratio of the galaxy populations and the M{sub .}/M{sub sph} ratio, which, according to numerous recent studies, evolves rapidly with redshift. The SMBH mass functions derived from the galaxy luminosity and mass functions show very good agreement, both in shape and in normalization. The resultant SMBH mass function and integrated mass density for the local epoch (z Almost-Equal-To 0) match well those derived independently by other studies. Consistent with other evidence for cosmic downsizing, the upper end of the mass function remains roughly constant since z Almost-Equal-To 2, while the space density of lower mass black holes undergoes strong evolution. We carefully assess the impact of various sources of uncertainties on our calculations. A companion paper uses the mass function derived in this work to determine the radiative efficiency of black hole accretion and constraints that can be imposed on the cosmological evolution of black hole spin.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

MASS GROWTH AND MERGERS: DIRECT OBSERVATIONS OF THE LUMINOSITY FUNCTION OF LRG SATELLITE GALAXIES OUT TO z = 0.7 FROM SDSS AND BOSS IMAGES

We present a statistical study of the luminosity functions of galaxies surrounding luminous red galaxies (LRGs) at average redshifts (z) = 0.34 and (z) = 0.65. The luminosity functions are derived by extracting source photometry around more than 40,000 LRGs and subtracting foreground and background contamination using randomly selected control fields. We show that at both studied redshifts the average luminosity functions of the LRGs and their satellite galaxies are poorly fitted by a Schechter function due to a luminosity gap between the centrals and their most luminous satellites. We utilize a two-component fit of a Schechter function plus a log-normal distribution to demonstrate that LRGs are typically brighter than their most luminous satellite by roughly 1.3 mag. This luminosity gap implies that interactions within LRG environments are typically restricted to minor mergers with mass ratios of 1:4 or lower. The luminosity functions further imply that roughly 35% of the mass in the environment is locked in the LRG itself, supporting the idea that mass growth through major mergers within the environment is unlikely. Lastly, we show that the luminosity gap may be at least partially explained by the selection of LRGs as the gap can be reproduced by sparsely sampling a Schechter function. In that case LRGs may represent only a small fraction of central galaxies in similar mass halos.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

USING KAPPA FUNCTIONS TO CHARACTERIZE OUTER HELIOSPHERE PROTON DISTRIBUTIONS IN THE PRESENCE OF CHARGE-EXCHANGE

Kappa functions have long been used in the analysis and modeling of suprathermal particles in various space plasmas. In situ observations of the supersonic solar wind show its distribution contains a cold ion core and power-law tail, which is well-represented by a kappa function. In situ plasma observations by Voyager, as well as observations of energetic neutral atom (ENA) spectra by the Interstellar Boundary Explorer (IBEX), showed that the compressed and heated inner heliosheath (IHS) plasma beyond the termination shock can also be represented by a kappa function. IBEX exposes the IHS plasma properties through the detection of ENAs generated by charge-exchange in the IHS. However, charge-exchange modifies the plasma as it flows through the IHS, and makes it difficult to ascertain the parent proton distribution. In this paper we investigate the evolution of proton distributions, initially represented by a kappa function, that experience losses due to charge-exchange in the IHS. In the absence of other processes, it is no longer representable by a single kappa function due to the energy-dependent, charge-exchange process. While one can still fit a kappa function to the evolving proton distribution over limited energy ranges, this yields fitting parameters (pseudo-density, pseudo-temperature, pseudo-kappa index) that depend on the energy range of the fit. We discuss the effects of fitting a kappa function to the IHS proton distribution over limited energy ranges, its dependence on the initial proton distribution properties at the termination shock, and implications for understanding the observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Combining protein sequences and structures with transformers and equivariant graph neural networks to predict protein function

Abstract Motivation Millions of protein sequences have been generated by numerous genome and transcriptome sequencing projects. However, experimentally determining the function of the proteins is still a time consuming, low-throughput, and expensive process, leading to a large protein sequence-function gap. Therefore, it is important to develop computational methods to accurately predict protein function to fill the gap. Even though many methods have been developed to use protein sequences as input to predict function, much fewer methods leverage protein structures in protein function prediction because there was lack of accurate protein structures for most proteins until recently. Results We developed TransFun—a method using a transformer-based protein language model and 3D-equivariant graph neural networks to distill information from both protein sequences and structures to predict protein function. It extracts feature embeddings from protein sequences using a pre-trained protein language model (ESM) via transfer learning and combines them with 3D structures of proteins predicted by AlphaFold2 through equivariant graph neural networks. Benchmarked on the CAFA3 test dataset and a new test dataset, TransFun outperforms several state-of-the-art methods, indicating that the language model and 3D-equivariant graph neural networks are effective methods to leverage protein sequences and structures to improve protein function prediction. Combining TransFun predictions and sequence similarity-based predictions can further increase prediction accuracy. Availability and implementation The source code of TransFun is available at https://github.com/jianlin-cheng/TransFun.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Taxonomic and functional restoration of tallgrass prairie soil microbial communities in comparison to remnant and agricultural soils

Abstract Restoring ecosystems requires the re-establishment of diverse soil microbial communities that drive critical ecosystem functions. In grasslands, restoration and management require the application of disturbances like fire and grazing. Disturbances can shape microbial taxonomic composition and potentially functional composition as well. We characterized taxonomic and functional gene composition of soil communities using whole genome shotgun metagenomic sequencing to determine how restored soil communities differed from pre-restoration agricultural soils and original remnant soils, how management affects soil microbes, and whether restoration and management affect the number of microbial genes associated with carbohydrate degradation. We found distinct differences in both taxonomic and functional diversity and composition among restored, remnant, and agricultural soils. Remnant soils had low taxonomic and functional richness and diversity, as well as distinct composition, indicating that restoration of agricultural soils does not re-create soil microbial communities that match remnants. Prescribed fire management increased functional diversity, which also was higher in more recently planted restorations. Finally, restored and post-fire soils included high abundances of genes encoding cellulose-degrading enzymes, so restorations and their ongoing management can potentially support functions important in carbon cycling.

Mason, Cayla N. (ORCID:0000000345983824)↗
DOE OSTI2021

Information content of higher order galaxy correlation functions

The shapes of galaxy N-point correlation functions can be used as standard rulers to constrain the distance–redshift relationship. The cosmological density fields traced by late-time galaxy formation are initially nearly Gaussian, and hence, all the cosmological information can be extracted from their two-point correlation function. Subsequent non-linear evolution under gravity, as well as halo and then galaxy formation, generates higher order correlation functions. Since the mapping of the initial to the final density field is, on large scales, invertible, it is often claimed that the information content of the initial field’s power spectrum is equal to that of all the higher order functions of the final, non-linear field. This claim implies that reconstruction of the initial density field from the non-linear field renders analysis of higher order correlation functions of the latter superfluous. We show that this claim is false when the N-point functions are used as standard rulers. Constraints available from joint analysis of the two and three-point correlation functions can, in some cases, exceed those offered by the initial power spectrum. We provide a mathematical justification for this claim and demonstrate it using a large suite of N-body simulations. In particular, we show that for the z = 0 real-space matter field in the limit of vanishing shot-noise, taking modes up to k max = 0.2 h Mpc −1 , using the bispectrum alone offers a factor of 2 reduction in the variance on the cosmic distance scale relative to that available from the linear power spectrum.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Estradiol associations with brain functional connectivity in postmenopausal women

Abstract Objective Previous studies have found that estrogens play a role in functional connectivity in the brain; however, little research has been done regarding how estradiol is associated with functional connectivity in postmenopausal women. The purpose of this study was to examine the relationship between estradiol and functional connectivity in postmenopausal women. Methods Structural and blood oxygenation level–dependent resting-state magnetic resonance imaging scans of 88 cognitively healthy postmenopausal individuals were obtained along with blood samples collected the same day as the magnetic resonance imaging to assess hormone levels. We generated connectivity values in CONN toolbox version 20.b, an SPM-based software. Results A regression analysis was run using estradiol level and regions of interest (ROI), including the hippocampus, parahippocampus, dorsolateral prefrontal cortex, and precuneus. Estradiol level was found to enhance parahippocampal gyrus anterior division left functional connectivity during ROI-to-ROI regression analysis. Estradiol enhanced functional connectivity between the parahippocampal gyrus anterior division left and the precuneus as well as the parahippocampal gyrus anterior division left and parahippocampal gyrus posterior division right. An exploratory analysis showed that years since the final menstrual period was related to enhanced connectivity between regions within the frontoparietal network. Conclusions These results illustrated the relationship between estradiol level and functional connectivity in postmenopausal women. They have implications for understanding how the functioning of the brain changes for individuals after menopause that may eventually lead to changes in cognition and behavior in older ages.

Obstetrics & Gynecology↗
DOE OSTI2023

Assessment of exchange-correlation functionals on oxygen vacancy formation energies of metal oxides

Oxygen vacancies are a common defect in oxides and play vital roles in many technological applications such as oxygen separation, catalytic reactors, solid oxide fuel cells, and solar thermochemical water splitting. The oxygen vacancy formation energy is directly related to the oxide reduction enthalpy and the vacancy concentration and, therefore, is a key quantity underlying these applications. The complexity of measuring oxygen vacancy formation energies experimentally suggests the utility of calculations based on density functional theory (DFT). However, the calculated results are strongly dependent on the exchange-correlation functionals and other parameters such as the Hubbard U . In this work, we compare the performance of the meta-GGA strongly constrained and appropriately normed (SCAN) functional and the commonly used semilocal generalized gradient approximation (GGA) functionals with experimental values of structural parameters, band gaps, magnetic structures, and oxygen vacancy formation energies for six representative oxides, i.e., CaMnO 3 , SrMnO 3 , LaMnO 3 , YMnO 3 , LaFeO 3 , and CeO 2 . Further, our results show that SCAN usually has better agreement with the experimental lattice constants and band gaps and larger magnetic moments than the commonly used GGA functionals. Although SCAN overestimates the oxygen vacancy formation energies of transition metal oxides and therefore requires unusually large Hubbard U values to reproduce the experimental reduction enthalpies, it does predict the correct oxygen vacancy formation energy of CeO 2 , which challenges the commonly used GGA and hybrid functionals. Our results underscore the challenges that exist in describing these complex oxides by DFT and may shed light on developing more accurate exchange-correlation functionals for oxygen vacancy formation energy calculations.

36 MATERIALS SCIENCE↗