Search NASA⌕ Search

SEARCH · Search NASA

Results for “consistent scheme”

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 19 records

Fundamental Studies of Bifunctional Catalysts for Tandem Reactions

During the past few years of our funding from DOE Basic Energy Sciences, we have focused on the catalytic conversion of biomass-derived 5-hydroxymethyl furfural (HMF) as a platform molecule for sustainable synthesis of value-added chemicals. The combination of the diene and aldehyde functionalities in HMF enables catalytic production of acetalized HMF derivatives with diol or epoxy reactants to allow reversible synthesis of norcantharimide derivatives upon Diels-Alder reaction with maleimides. Whereas the electron-withdrawing nature of aldehyde group inhibits the reactivity of diene group in HMF to undergo Diels-Alder reactions with dienophiles, by converting the aldehyde group in HMF to an electron-donating hydroxyl group, HMF can append maleimide-based chemicals through Diels-Alder coupling. Acetalization of HMF not only produces an electron-rich diene for Diels-Alder reaction, but it also alters the reversibility of acetal formation and thereby allows for the controlled release of chemicals by hydrolysis. Therefore, we synthesized various norcantharimide derivatives from biomass-derived HMF by acetalization over an acid catalyst, Amberlyst-15, followed by Diels-Alder reaction with maleimide. The norcantharimides release the starting materials by retro Diels-Alder reaction that is triggered by acetal hydrolysis under acidic (≤pH 3) conditions. In subsequent work, we explored the synthesis and properties of functional polyurethanes and polyesters with tunable properties from biomass-derived (HMF)-Acetone-HMF (HAH) monomers. HAH can be selectively hydrogenated over Cu and Ru catalysts to produce partially-hydrogenated HAH (PHAH) and fully-hydrogenated HAH (FHAH). HAH contains functional groups, including hydroxyl, furan, enone, and ketone functionalities, that can be exploited to further tune the polymer properties. The π-electron conjugation between the enone and furan groups in HAH monomer stabilizes the furan group to prevent the uncontrolled degradation of furans during reactions, such as etherification. Therefore, HAH is a renewable monomer that can be used to synthesize polymers, having high molecular weight, symmetric functionalities, an inexpensive production price, and being derived from renewable resources without competition with food resources. We demonstrated the synthesis of functional polyurethanes and a renewable polyester from HAH-derived monomers, and 4,4’-methylenebis(phenyl isocyanate) (MDI). We showed that it is essential to be able to selective hydrogenate C=C double bonds in HAH to be able to achieve effective Diels Alder reactions of the furan ring. Therefore, we studied the hydrogenation at temperatures from 313 – 393 K of HAH over Pd, Ru, and Cu based catalysts. HAH was selectively hydrogenated to produce partially-hydrogenated monomers (PHAH) over Cu and Ru catalysts and to fully-hydrogenated HAH monomers (FHAH) over the Ru catalyst. Pd based catalysts yielded a mixture of partially and fully hydrogenated monomers. Reaction kinetics models were employed to quantify the kinetic behavior for hydrogenation over Ru, Cu, and Pd catalysts. A 5-step pathway exhibited over Pd and Ru catalysts consists of both series and parallel reaction steps, where HAH is both converted to fully hydrogenated products sequentially via series reactions of partially hydrogenated intermediates, as well as converted directly in parallel reactions to form the fully hydrogenated products. In contrast, a 3-step pathway over the Cu catalyst consists only of the consecutive reaction steps, where the final product was formed via series reactions of intermediate products. Additionally, reaction over the Cu catalyst did not hydrogenate the furan rings of the HAH molecule and yielded a different final product than those hydrogenation over Pd and Ru catalysts. Using the results from our detailed reaction kinetics studies, we developed universal conditions for the maximum production of the various hydrogenated products for both batch and plug flow systems. In subsequent work, we employed reaction schemes to optimize the yields of products from selective hydrogenations of HAH in isopropanol solvent and in the presence of liquid water. Reaction schemes consisting of 7, 9, and 11 steps were examined to describe the rates of formation of the observed products and reaction intermediates for hydrogenation of HAH over Ru and Pd catalysts, and a 3-step scheme was studied over Cu catalysts. Rate constants and activation energies were calculated using these reaction schemes, and we then apply these reaction schemes to explore the effects of water addition on the hydrogenation pathways. The effects of water addition to isopropanol (IPA) solvents on the hydrogenation of HAH were markedly different over Pd, Ru, and Cu catalysts. Over the Pd catalyst, the addition of water to IPA increased hydrogenation rates and promoted hydrogenation of furan rings. The addition of water to IPA yielded significant carbon losses over the Ru catalyst, and slowed hydrogenation steps over Cu, while significantly inhibiting hydrogenation of the ketone group. This behavior opened routes toward increased production rates of PHAH=O (a partially hydrogenated form of HAH containing a C=O bond), a product in which the diene groups of the furan rings were not hydrogenated. Addition of water also allowed increased feed concentrations of HAH that were previously not possible in pure IPA solvents.

09 BIOMASS FUELS↗

A hybrid nodal-staggered pseudo-spectral electromagnetic particle-in-cell method with finite-order centering

Electromagnetic particle-in-cell (PIC) codes are widely used to perform computer simulations of a variety of physical systems, including fusion plasmas, astrophysical plasmas, plasma wakefield particle accelerators, and secondary photon sources driven by ultra-intense lasers. In a PIC code, Maxwell's equations are solved on a grid with a numerical method of choice. This article focuses on pseudo-spectral analytical time-domain (PSATD) algorithms and presents a novel hybrid PSATD PIC scheme that combines the respective advantages of standard nodal and staggered methods. The novelty of the hybrid scheme consists in using finite-order centering of grid quantities between nodal and staggered grids, in order to combine the solution of Maxwell's equations on a staggered grid with the deposition of charges and currents and the gathering of electromagnetic forces on a nodal grid. The correctness and performance of the novel hybrid scheme are assessed by means of numerical tests that employ different classes of PSATD equations in a variety of physical scenarios, ranging from the modeling of electron-positron pair creation in vacuum to the simulation of laser-driven and particle beam-driven plasma wakefield acceleration. It is shown that the novel hybrid scheme offers significant numerical and computational advantages, compared to purely nodal or staggered methods, for all the test cases presented.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Direct numerical simulation of a turbulent methane/air flame impinged by a sub-breakdown electric field

Here, the effects of an external electric field on a turbulent methane/air diffusion flame are analyzed in this work using direct numerical simulations. The analyzed configuration consists of a temporally evolving mixing layer of air and a mixture of methane and nitrogen at 1 atm that is impinged by an electric field in the direction normal to the mean mixing plane. The combustion and chemi-ionization reactions involved in the flow are modeled using finite-rate chemistry and a reduced reaction scheme consisting of 26 species and 134 reactions. The mass diffusion and ion-wind effects are modeled using a detailed description of the diffusion coefficients and electric mobilities based on kinetic theory. The presented calculations show that the turbulence generated within the mixing layer is mostly unaffected by the applied electric field for the configuration under exam. In fact, the electric body force is developed away from the mixing region where the flow is uniform. Conversely, the turbulence is able to introduce very high intermittency in the electrically charged species concentration and, consequently, in the electric body force. Such intermittency will constitute a challenge in the future formulation of combustion models that take into account ion-wind effects.

09 BIOMASS FUELS↗

DG-IMEX method for a two-moment model for radiation transport in the $\mathscr{O}$($v$/$c$) limit

Here, we consider neutral particle systems described by moments of a phase-space density and propose a realizability-preserving numerical method to evolve a spectral two-moment model for particles interacting with a background fluid moving with nonrelativistic velocities. The system of nonlinear moment equations, with special relativistic corrections to $\mathscr{O}$($v$/$c$), expresses a balance between phase-space advection and collisions and includes velocity-dependent terms that account for spatial advection, Doppler shift, and angular aberration. The model is conservative for the correct $\mathscr{O}$($v$/$c$) Eulerian-frame number density and is consistent, to $\mathscr{O}$($v$/$c$), with Eulerian-frame energy and momentum conservation. This model is closely related to the one promoted by Lowrie et al. and similar to models currently used to study transport phenomena in large-scale simulations of astrophysical environments. The proposed numerical method is designed to preserve moment realizability, which guarantees that the moments correspond to a nonnegative phase-space density. The realizability-preserving scheme consists of the following key components: (i) a strong stability-preserving implicit-explicit (IMEX) time-integration method; (ii) a discontinuous Galerkin (DG) phase-space discretization with carefully constructed numerical uxes; (iii) a realizability-preserving implicit collision update; and(iv) a realizability-enforcing limiter. In time integration, nonlinearity of the moment model necessitates solution of nonlinear equations, which we formulate as fixed-point problems and solve with tailored iterative solvers that preserve moment realizability with guaranteed global convergence. We also analyze the simultaneous Eulerian-frame number and energy conservation properties of the semi-discrete DG scheme and propose a "spectral redistribution" scheme that promotes Eulerian-frame energy conservation. Through numerical experiments, we demonstrate the accuracy and robustness of this DG-IMEX method and investigate its Eulerian-frame energy conservation properties.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Central Compact Finite‐Difference Scheme With High Spectral Resolution for KdV Equation

This work presents a combination of cell‐node and cell‐centered compact finite difference scheme for the approximation of third derivatives involved in Korteweg–de Vries (KdV) equations. This approach employs a half‐shifted derivative construction at cell centers, avoiding the need for compact interpolation, thereby removing transfer errors; hence, it improves spectral resolution and maintains high‐order accuracy. Fourier analysis is performed to show the spectral properties of the proposed formulation, which provides higher spectral resolutions as compared to node‐based compact schemes. A filtering strategy is incorporated to suppress high‐frequency oscillations without compromising the accuracy of the numerical scheme, and the total variation diminishing Runge Kutta (TVDRK3) method is applied for time integration. Numerical experiments on linear, nonlinear, and coupled KdV systems are conducted, and a comparative analysis with cell‐node compact schemes confirms that the proposed scheme consistently reduces errors by up to an order of magnitude and achieves high spectral resolution properties.

97 MATHEMATICS AND COMPUTING↗

Inferring colloidal interaction from scattering by machine learning

A machine learning solution for the potential inversion problem in elastic scattering is outlined. The inversion scheme consists of two major components, a generative network featuring a variational autoencoder which extracts the targeted static two-point correlation functions from experimentally measured scattering cross sections, and a Gaussian process framework which probabilistically infers the relevant structural parameters from the inverted correlation functions. Via a case study of charged colloidal suspensions, the feasibility of this approach for quantitative study of molecular interaction is critically benchmarked and its merit over existing deterministic approaches, in terms of numerical accuracy and computationally efficiency, is demonstrated.

36 MATERIALS SCIENCE↗

2D kinetic-ion simulations of inverted corona fusion targets

Laser-driven “inverted corona” fusion targets have attracted interest as a low-convergence neutron source and platform for studying kinetic physics. The scheme consists of a hollow or gas-filled spherical shell made of deuterated plastic. The shell has one or more laser entrance holes (LEH), resembling a spherical hohlraum. The laser passes through the LEH’s and illuminates the interior surface of the shell, ablating a plasma that travels inward towards the target center. Long ion mean free paths in the converging plasma can lead to significant interpenetration, atomic mix, and other kinetic effects. Here, in this work we report on numerical simulations of inverted corona targets using the kinetic-ion, fluid–electron hybrid particle-in-cell (PIC) approach in 2D RZ geometry. 2D simulations suggest that shape effects do not have a significant impact on plasma evolution and observed yield trends are primarily the result of 1D kinetic mix mechanisms. Simulations are also compared against available experimental data recorded at the OMEGA laser facility. In particular, synthetic x-ray emission images show good qualitative agreement with experimental results, albeit with an apparent timing discrepancy for the two-sided vacuum target. More generally, we demonstrate the potential of hybrid-PIC simulations for full-system modeling and experimental design, including collisional absorption of laser energy, plasma evolution, mix, and fusion burn.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Physics or Knob‐Tuning? Tropical Anvil Peak Is Captured by GCMs

Tropical anvil clouds peak near 200 hPa and significantly impact Earth's climate, yet its physical realism in coarse‐resolution General Circulation Models (GCMs) remains debated. We examine anvil cloud formation by performing simulations with a GCM with a hierarchy of cloud fraction schemes ranging from a complex prognostic Tiedtke scheme to a simple binary scheme. All schemes consistently reproduce the anvil peak. The robust anvil peak arises because extremely cold temperatures at the upper troposphere facilitate frequent saturation events, producing clouds that disproportionately influence mean cloud fraction. Sensitivity experiments with enhanced evaporation of cloud condensate unexpectedly show increased anvil coverage, highlighting how slight evaporative moistening reinforces local saturation in cold upper‐tropospheric conditions. These results demonstrate that the tropical anvil cloud peak emerges from fundamental thermodynamic constraints, rather than specific cloud fraction parameterization choices.

54 ENVIRONMENTAL SCIENCES↗

Optimized nuclear energy density functionals including long-range pion contributions

Nuclear energy density functionals successfully reproduce properties of nuclei across almost the entire nuclear chart. However, nearly all available functionals are phenomenological in nature and lack a rigorous connection to systematically improvable nuclear forces. This issue might be solved with an energy density functional obtained from first principles. As an intermediate step towards this goal we construct the GUDE family of functionals that is obtained from a hybrid scheme consisting of long-range pion-exchange contributions derived from chiral effective field theory at the Hartree-Fock level and a phenomenological Skyrme part. When including pion contributions beyond next-to-leading order in the chiral expansion, we find significant improvements over a reference Skyrme functional constructed following the same protocol. We analyze the importance of different pion contributions and identify which terms drive the observed improvements. Since pions are incorporated without adding further optimization parameters to the functionals, the improvements can be attributed to the functional form of these terms. Our work therefore suggests that the considered chiral contributions constitute useful ingredients for true energy density functionals. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Effect of thermochemical treatments on laser-induced luminescence spectra from strontium titanate: comparison with swift ion-beam irradiation experiments

Results recently reported on the effect of thermochemical treatments on the (He-Cd) laser-excited emission spectra of strontium titanate (STO) are re-analyzed here and compared with results obtained under ion-beam irradiation. Contributing bands centered at 2.4 eV and 2.8 eV, which appear under laser excitation, present intensities dependent upon previous thermal treatments in oxidizing (O 2 ) or reducing atmosphere (H 2 ). As a key result, the emission band centered at 2.8 eV is clearly enhanced in samples exposed to a reducing atmosphere. From a comparison with the ionoluminescence data, it is concluded that the laser-excited experiments can be rationalized within a framework developed from ion-beam excitation studies. In particular, the band at 2.8 eV, sometimes attributed to oxygen vacancies, behaves as expected for optical transitions from conduction-band (CB) states to the ground state level of the self-trapped exciton center. The band at 2.0 eV reported in ion-beam irradiated STO, and attributed to oxygen vacancies, is not observed in laser-excited crystals. As a consequence of our analysis, a consistent scheme of electronic energy levels and optical transitions can now be reliably offered for strontium titanate.

74 ATOMIC AND MOLECULAR PHYSICS↗

Robust Optimized Pulse Schemes for Atomic Fountain Interferometry

The robustness of an atomic fountain interferometer with respect to variations in the initial velocity of the atoms and deviations from the optimal pulse amplitude is examined. We numerically simulate the dynamics of an interferometer in momentum space with a maximum separation of and map out the expected signal contrast depending on the variance of the initial velocity distribution and the value of the laser field amplitude. We show that an excitation scheme based on rapid adiabatic passage significantly enhances the expected signal contrast, compared to the commonly used scheme consisting of a series of and pulses. We demonstrate further substantial increase of the robustness by using optimal control theory to identify splitting and swapping pulses that perform well on an ensemble average of pulse amplitudes and velocities. Our results demonstrate the ability of optimal control to significantly enhance future implementations of atomic fountain interferometry.

74 ATOMIC AND MOLECULAR PHYSICS↗

WRF-Comfort: simulating microscale variability in outdoor heat stress at the city scale with a mesoscale model

Abstract. Urban overheating and its ongoing exacerbation due to global warming and urban development lead to increased exposure to urban heat and increased thermal discomfort and heat stress. To quantify thermal stress, specific indices have been proposed that depend on air temperature, mean radiant temperature (MRT), wind speed, and relative humidity. While temperature and humidity vary on scales of hundreds of meters, MRT and wind speed are strongly affected by individual buildings and trees and vary on the meter scale. Therefore, most numerical thermal comfort studies apply microscale models to limited spatial domains (commonly representing urban neighborhoods with building blocks) with resolutions on the order of 1 m and a few hours of simulation. This prevents the analysis of the impact of city-scale adaptation and/or mitigation strategies on thermal stress and comfort. To solve this problem, we develop a methodology to estimate thermal stress indicators and their subgrid variability in mesoscale models – here applied to the multilayer urban canopy parameterization BEP-BEM within the Weather Research and Forecasting (WRF) model. The new scheme (consisting of three main steps) can readily assess intra-neighborhood-scale heat stress distributions across whole cities and for timescales of minutes to years. The first key component of the approach is the estimation of MRT in several locations within streets for different street orientations. Second, mean wind speed and its subgrid variability are downscaled as a function of the local urban morphology based on relations derived from a set of microscale LES and RANS simulations across a wide range of realistic and idealized urban morphologies. Lastly, we compute the distributions of two thermal stress indices for each grid square, combining all the subgrid values of MRT, wind speed, air temperature, and absolute humidity. From these distributions, we quantify the high and low tails of the heat stress distribution in each grid square across the city, representing the thermal diversity experienced in street canyons. In this contribution, we present the core methodology as well as simulation results for Madrid (Spain), which illustrate strong differences between heat stress indices and common heat metrics like air or surface temperature both across the city and over the diurnal cycle.

Geology↗

Hyperspectral remote sensing-based plant community map for region around NGEE-Arctic intensive research watersheds at Seward Peninsula, Alaska, 2017-2019

Using airborne hyperspectral remote sensing data from NASA Airborne Visible-Infrared Imaging Spectrometer- Next Generation (AVIRIS-NG) platforms in a region near NGEE-Arctic intensive watersheds at Seward peninsula of Alaska, high resolution (5m) maps of plant community distribution were developed and included in this data collected. AVIRIS-NG data collected over 2017-2019 period were used to develop deep neural networks, trained using vegetation plot observations collected at NGEE-Arctic watersheds at Kougarok, Council and Teller. A hierarchical vegetation classification scheme consisting of six classes at Level I, and 16 classes at Level II contained in two .txt files were used to developed the plant community maps for the region. Two geospatial raster data files (.tif) at both thematic levels are shared in this data collection. Data files in this collection use Alaska Albers Equal Area projection. Readme files available in three formats (*.html, *.md, *.pdf) and one *.png visualization map.The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research.The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska.Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗

Study of Seamless Microgrid Transition Operation Using Grid-Forming Inverters

This paper investigates operational techniques to achieve seamless (smooth) microgrid (MG) transitions by dispatching a grid-forming (GFM) inverter. In traditional approaches, the GFM inverter must switch between grid-following (GFL) and GFM control modes during MG transition operation. Today's inverter technology allows GFM inverters to always operate in GFM control mode, so it is worth exploring how to use them to achieve smooth MG transition operation. This paper proposes three operational techniques: a traditional scheme of switching between GFL and GFM control; a new scheme of consistent GFM control and shifting the droop intercept up before islanding operation; and a new scheme of consistent GFM control and shifting the droop intercept up before synchronization operation. A full hardware setup is established to compare the three techniques and showcase their implementations in real-world applications. The results show that the third technique outperforms the others and exhibits the best transition performance because the GFM inverter maintains the same operating points during the transition operation. Therefore, we conclude that ensuring smooth MG transition operation requires that the GFM inverter(s) maintain the same operating points (v, f, P, Q, and phase angle) during the transition operation in addition to minimizes the point of common coupling power flow.

grid-forming control↗

Study of Seamless Microgrid Transition Operation Using Grid-Forming Inverters

This paper investigates operational techniques to achieve seamless (smooth) microgrid (MG) transitions by dispatching a grid-forming (GFM) inverter. In traditional approaches, the GFM inverter must switch between grid-following (GFL) and GFM control modes during MG transition operation. Today's inverter technology allows GFM inverters to always operate in GFM control mode, so it is worth exploring how to use them to achieve smooth MG transition operation. This paper proposes three operational techniques: a traditional scheme of switching between GFL and GFM control; a new scheme of consistent GFM control and shifting the droop intercept up before islanding operation; and a new scheme of consistent GFM control and shifting the droop intercept up before synchronization operation. A full hardware setup is established to compare the three techniques and showcase their implementations in real-world applications. The results show that the third technique outperforms the others and exhibits the best transition performance because the GFM inverter maintains the same operating points during the transition operation. Therefore, we conclude that ensuring smooth MG transition operation requires that the GFM inverter(s) maintain the same operating points (v, f, P, Q, and phase angle) during the transition operation in addition to minimizes the point of common coupling power flow.

grid-forming↗

Study of Seamless Microgrid Transition Operation Using Grid-Forming Inverters: Preprint

This paper investigates operational techniques to achieve smooth microgrid (MG) transitions by dispatching the grid-forming (GFM) inverter. In traditional approaches, the GFM inverter must switch between grid-following (GFL) and GFM control mode during MG transition operation. Today's inverter technology allows GFM inverters to always operate in GFM control mode, so it is worthy to explore how to use them to achieve smooth MG transition operation. This paper proposes three operational techniques: a traditional scheme of switching between GFL and GFM control; a new scheme of consistent GFM control and shifting the droop intercept up before islanding operation; and a new scheme of consistent GFM control and shifting the droop intercept up before synchronization operation. A full hardware setup is established to compare the three techniques and showcase their implementations in real-world applications. As expected, the third technique outperforms the others and exhibits the best transition performance because the GFM inverter maintains the same operating points during the transition operation. Therefore, we conclude that ensuring smooth MG transition operation requires the GFM inverter(s) to maintain the same operating points (v, f, P, Q, and phase angle) during the transition operation in addition to minimizes the point of common coupling power flow.

grid-following control↗

Influence of Markovianity and self-consistency on time-resolved spectral functions of driven quantum systems

We present a systematic comparison of the real-time Dyson expansion (RTDE) with established nonequilibrium Green's function (GF) approaches for simulating driven, interacting quantum systems. Focusing on density matrix dynamics, time-off-diagonal GFs, and time-resolved photoemission spectra, we benchmark RTDE against fully self-consistent Kadanoff-Baym equation (KBE) calculations, the generalized Kadanoff-Baym ansatz, and exact diagonalization for small systems using second-order many-body perturbation theory. Using a driven two-band Hubbard model, we show that mean-field single-particle density matrix trajectories provide a reliable baseline for RTDE across a broad range of interaction strengths and excited-carrier populations. Further, RTDE accurately captures correlation effects in the GFs, including long-lived oscillations and revivals that are strongly suppressed by the overdamping inherent to self-consistent KBE schemes. As a consequence, RTDE resolves rich nonequilibrium spectral structure in time-resolved photoemission, such as interaction- and population-dependent quasiparticle splittings and band gap renormalization, which are largely washed out in self-consistent approaches yet are present in exact solutions. Furthermore, our results demonstrate that RTDE bridges the gap between mean-field propagation and full two-time KBE simulations, retaining favorable linear scaling while capturing essential dynamical correlations relevant for ultrafast spectroscopy.

Electronic structure↗

Optimally tuned starting point for single-shot GW calculations of solids

The dependence of ab initio many-body perturbation theory within the $\textit{GW}$ approximation on the eigensystem used in calculating quasiparticle corrections limits this method's predictive power. Here, we investigate the accuracy of the recently developed Wannier-localized optimally tuned screened range-separated hybrid (WOT-SRSH) functional as a generalized Kohn-Sham starting point for single-shot $\textit{GW}$ ($G_0W_0$) calculations for a range of semiconductors and insulators. Comparison to calculations based on well-established functionals, namely, PBE, PBE0, and HSE, as well as to self-consistent GW schemes and to experiment, shows that band gaps computed via $G_0W_0$@WOT-SRSH have a level of precision and accuracy that is comparable to that of more advanced methods such as quasiparticle self-consistent $\textit{GW}$ and eigenvalue self-consistent $\textit{GW}$. In this work, we also find that $G_0W_0$@WOT-SRSH improves the description of states deeper in the valence band manifold. Finally, we show that $G_0W_0$@WOT-SRSH significantly reduces the sensitivity of computed band gaps to ambiguities in the underlying WOT-SRSH tuning procedure.

36 MATERIALS SCIENCE↗