Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “GW method”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Life cycle assessment of utilizing freshly cut urban wood: A case study

Large amounts of urban wood resources are generated through tree removals in urban areas. Therefore, there is a growing interest to improve the environmental performance of the urban wood supply chain by quantifying the environmental impacts of creating high-value products from urban tree removals to enhance its utilization. By surveying existing urban wood utilization operations for primary data in two major cities: Baltimore, Maryland, and Milwaukee, Wisconsin, a new life cycle inventory (LCI) dataset per m 3 of dry urban hardwood lumber and live edge slabs was created. Incorporating the new LCI data into life cycle assessment method, the total global warming (GW) impact of converting urban trees to kiln-dried lumber was estimated to be 122 kg CO 2 eq and 336 kg CO 2 eq per m 3 of lumber produced for Baltimore and Milwaukee, respectively. In both cases, the total environmental impact of the product was dominated by the conversion of urban trees to kiln dried lumber at the sawmill processing stage. Using scenario analysis for the source of heat generation used in kiln (force) drying wood, substituting wood fuel from coproducts for natural gas in a boiler for Milwaukee resulted in a substantial reduction in fossil-derived GW impact. Urban and traditional hardwood lumber production follow similar GW impact trajectories over their life-cycle stages with the lumber processing stage having the highest environmental impacts for both. Finally, increasing the overall lumber yield and decreasing kiln (force) drying would substantially improve the environmental performance of urban hardwood lumber and help make it comparable to traditional lumber.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

The Challenges Ahead for Multimessenger Analyses of Gravitational Waves and Kilonova: A Case Study on GW190425

In recent years, there have been significant advances in multimessenger astronomy due to the discovery of the first, and so far only confirmed, gravitational wave event with a simultaneous electromagnetic (EM) counterpart, as well as improvements in numerical simulations, gravitational wave (GW) detectors, and transient astronomy. This has led to the exciting possibility of performing joint analyses of the GW and EM data, providing additional constraints on fundamental properties of the binary progenitor and merger remnant. Here, in this work, we present a new Bayesian framework that allows inference of these properties, while taking into account the systematic modeling uncertainties that arise when mapping from GW binary progenitor properties to photometric light curves. We extend the relative binning method presented in Zackay et al. to include extrinsic GW parameters for fast analysis of the GW signal. The focus of our EM framework is on light curves arising from r-process nucleosynthesis in the ejected material during and after merger, the so-called kilonova, and particularly on black hole–neutron star systems. As a case study, we examine the recent detection of GW190425, where the primary object is consistent with being either a black hole or a neutron star. We show quantitatively how improved mapping between binary progenitor and outflow properties, and/or an increase in EM data quantity and quality are required in order to break degeneracies in the fundamental source parameters.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Dynamic mode decomposition for extrapolating nonequilibrium Green's-function dynamics

The Hartree-Fock generalized Kadanoff-Baym ansatz (HF-GKBA) offers an approximate numerical procedure for propagating the two-time nonequilibrium Green's function (NEGF). Here, using the GW self-energy, we compare the HF-GKBA to exact results for a variety of systems with long- and short-range interactions, different two-body interaction strengths, and various nonequilibrium preparations. We find excellent agreement between the HF-GKBA and exact time evolution in models when more realistic long-range exponentially decaying interactions are considered. Furthermore, this agreement persists for long times and for intermediate to strong interaction strengths. In large systems, HF-GKBA becomes prohibitively expensive for long-time evolutions. For this reason, we look at the use of dynamical mode decomposition (DMD) to reconstruct long-time NEGF trajectories from a sample of the initial trajectory. Using no more than 16% of the total time evolution, we reconstruct the total trajectory with high fidelity. Our results show the potential for DMD to be used in conjunction with HF-GKBA to calculate long-time trajectories in large-scale systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Challenging conventional assumptions in PV: a high-throughput open-air approach to low-cost perovskite module production

Perovskite solar modules (PSMs) offer a promising pathway to low-cost photovoltaics, yet their commercialization is challenged by manufacturing scalability, device uniformity, additive costs, interlayer complexity, and module stability. This study introduces a comprehensive technoeconomic analysis of single junction PSM's and projections for tandem perovskite-Si modules that integrate all materials and manufacturing steps, module performances, projected lifetimes, and manufacturing costs across scales. Here, we highlight an open-air manufacturing approach to fabricate all active layers of serially interconnected PSMs, including electrodes and charge transport layers, enabling high-throughput production without inert or vacuum environments. The analysis reveals two orders of magnitude throughput enhancement and cost reductions of 24% in all-open-air production, escalating to over 60% at 1 GW factory capacity compared to conventional methods. Levelized cost of energy (LCOE) projections for utility-scale installations over 30 years, accounting for module replacement and recycling, demonstrate the potential to achieve the 2030 US target of $0.03 per kWh with realistic 7–11-year PSM lifetimes, outperforming incumbent silicon-based modules. Neither four terminal (4T) nor two terminal (2T) tandem-Si PSMs improve over single junction perovskite or silicon LCOE regardless of higher efficiencies at any modeled lifetime. Addressing PSM technical challenges with a cost-modeling framework guides commercialization efforts and provides a convincing pathway for challenging incumbent Si-based PV.

14 SOLAR ENERGY↗
DOE OSTI2023

Bending effects and optical properties of WSe 2 nanoribbons of topological phase

A WSe 2 monolayer of 1T' phase is a large band gap quantum spin Hall insulator, supporting dissipationless charge and spin transports through the topologically protected edge states. Here, in this work, we explore the nanoribbon forms of 1T' phase WSe 2 by first-principles density functional calculations and the many-body perturbation GW and Bethe-Salpeter equation method. We found that the 1T' WSe 2 nanoribbon can show topological edge states with a ribbon width of ~4-6 nm. Those edge bands show crossing through the Fermi level an odd number of times, with one kind of spin-polarization connecting the valence band continuum and conduction band continuum. The topological features of the edge bands hold even under small and medium bending in the nanoribbon, while large bending induces large band splitting, resulting in a topological switch-off in the edge bands. The semiconducting 1T' WSe2 nanoribbon shows a large tunability with bending in optical absorption spectra and exciton states. The lowest-energy exciton is changed from optically dark in the flat nanoribbon to bright in the bent nanoribbons. These properties in the 1T' WSe 2 nanoribbons suggest potential applications in controllable quantum electronics and exciton-based quantum information processes.

36 MATERIALS SCIENCE↗
DOE OSTI2022

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↗
DOE OSTI2023

New insights into the electronic structure of α-U and δ-Pu

Here, this work presents the results of a theoretical study of the electronic structure of two actinide metals, α -U and δ -Pu. We compare our ab-initio results obtained with the recently developed self-consistent Vertex corrected GW approach with previously published experimental measurements such as photo-electron spectroscopy, for the occupied density of states, and bremsstralung isochromat spectroscopy (BIS) and inverse photo-electron spectroscopy (IPES), for the unoccupied density of states. Our ab-initio approach includes all important relativistic effects (it is based on Dirac’s equation) and it represents the first application of the Vertex corrected GW approach in the physics of actinides. Overall, our theoretical results are in good agreement with the experimental data, which supports the level of approximations which our theoretical method is based upon. By comparing our vertex corrected GW results with our results obtained with less sophisticated approaches (local density approximation and self-consistent GW) we differentiate the strength of correlation effects in Uranium and Plutonium. Also, our theoretical results allow us to elucidate the subtle differences between the previously published experimental BIS and IPES data on the unoccupied density of states in α -U.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Fabrication of THz corrugated wakefield structure and its high power test

We present overall process for developing terahertz (THz) corrugated structure and its beam-based measurement results. 0.2-THz corrugated structures were fabricated by die stamping method as the first step demonstration towards GW THz radiation source and GV/m THz wakefield accelerator. 150-μm thick disks were produced from an OFHC (C10100) foil by stamping. Two types of disks were stacked alternately to form 46 mm structure with ~ 170 corrugations. Custom assembly was designed to provide diffusion bonding with a high precision alignment of disks. The compliance of the fabricated structure have been verified through beam-based wakefield measurement at Argonne Wakefield Accelerator Facility. Both measured longitudinal and transverse wakefield showed good agreement with simulated wakefields. Measured peak gradients, 9.4 MV/m/nC for a long single bunch and 35.4 MV/m/nC for a four bunch trains, showed good agreement with the simulation.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Spin-polarization anisotropy controlled by bending in tungsten diselenide nanoribbons and tunable excitonic states

A WSe 2 monolayer shows many interesting properties due to its spin-orbit coupling induced spin splitting in bands around the Fermi level and the spin-valley configuration. The orientation of the spin polarization in the relevant bands is crucial for the nature of exciton states and the optical valley selectivity. In this work, we studied the WSe 2 nanoribbons under different mechanical bending curvatures and electron/hole doping with density functional theory and their optical absorption and excitonic states with many-body perturbation GW and BSE (Bethe-Salpeter equation) methods. We found that the WSe 2 nanoribbons can exhibit an enhanced SOC effect and a spatially varying spin polarization in bands around the Fermi level under bending. The spin-polarization can show an anisotropy (or asymmetry) in those nearly degenerate bands, leading to a controllable magnetism via bending and electron/hole doping to the nanoribbons, suggesting a potential application in compact and controllable magnetic nanodevices and spintronics. The optical absorption spectrum of the nanoribbon presents a large tunability with bending within the near infrared region of about 0.4 to 1.5 eV, showing an enhanced absorption at a large bending condition. Finally, the exciton states generally show mixed or various spin configuration in the electron and hole pairs that are controlled by bending, potentially useful for applications in spin-based quantum information processes.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Optically controlled single-valley exciton doublet states with tunable internal spin structures and spin magnetization generation

Manipulating quantum states through light–matter interactions has been actively pursued in two-dimensional materials research. Significant progress has been made toward the optical control of the valley degrees of freedom in semiconducting monolayer transition-metal dichalcogenides, based on doubly degenerate excitons from their two distinct valleys in reciprocal space. Here, we introduce a type of optically controllable doubly degenerate exciton states that come from a single valley, dubbed as single-valley exciton doublet (SVXD) states. They are unique in that their constituent holes originate from the same valence band, making possible the direct optical control of the spin structure of the excited constituent electrons. Combining ab initio GW plus Bethe–Salpeter equation ( GW -BSE) calculations and a theoretical analysis method, we demonstrate such SVXD in substrate-supported monolayer bismuthene—which has been successfully grown using molecular beam epitaxy. In each of the two distinct valleys in the Brillouin zone, strong spin–orbit coupling and C 3 v symmetry lead to a pair of degenerate 1s exciton states (the SVXD states) with opposite spin configurations. Any coherent linear combinations of the SVXD in a single valley can be excited by light with a specific polarization, enabling full manipulation of their internal spin configurations. In particular, a controllable net spin magnetization can be generated through light excitation. Our findings open routes to control quantum degrees of freedom, paving the way for applications in spintronics and quantum information science.

2D materials↗
DOE OSTI2023

Characterizing a supernova’s standing accretion shock instability with neutrinos and gravitational waves

Here, we perform a novel multi-messenger analysis for the identification and parameter estimation of the Standing Accretion Shock Instability (SASI) in a core collapse supernova with neutrino and gravitational wave (GW) signals. In the neutrino channel, this method performs a likelihood ratio test for the presence of SASI in the frequency domain. For gravitational wave signals we process an event with a modified constrained likelihood method. Using simulated supernova signals, the properties of the Hyper-Kamiokande neutrino detector, and O3 LIGO Interferometric data, we produce the two-dimensional probability density distribution (PDF) of the SASI activity indicator and calculate the probability of detection P D as well as the false identification probability P FI . We discuss the probability to establish the presence of the SASI as a function of the source distance in each observational channel, as well as jointly. Compared to a single-messenger approach, the joint analysis results in P D (at P FI = 0.1) of SASI activities that is larger by up to ≈ 40% for a distance to the supernova of 5 kpc. We also discuss how accurately the frequency and duration of the SASI activity can be estimated in each channel separately. Our methodology is suitable for implementation in a realistic data analysis and a multi-messenger setting.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Exploiting a derivative discontinuity estimate for accurate G0W0 ionization potentials and electron affinities

Abstract The GW approximation has become an important tool for predicting charged excitations of isolated molecules and condensed systems. Its popularity can be attributed to many factors, including a favorable scaling and relatively good accuracy. In practical applications, the GW is often performed as a one-shot perturbation known as G 0 W 0 . Unfortunately, G 0 W 0 suffers from a strong starting point dependence and is often not as accurate as one would need. Self-consistent GW methodologies alleviate these problems but come with a marked increase in computational cost. In this manuscript, we propose the use of an estimate of the exchange-correlation derivative discontinuity to provide a remarkably good starting point for G 0 W 0 calculations, yielding ionization potentials and electron affinities with eigenvalue self-consistent GW quality at no additional cost. We assess the quality of the resulting methodology with the GW 100 benchmark set and compare its advantages over other similar methods.

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

Importance of nonuniform Brillouin zone sampling for ab initio Bethe-Salpeter equation calculations of exciton binding energies in crystalline solids

Excitons are prevalent in semiconductors and insulators, and their binding energies are critical for optoelectronic applications. The state-of-the-art method for first-principles calculations of excitons in extended systems is the ab initio GW-Bethe-Salpeter equation (BSE) approach, which can require a fine sampling of reciprocal space to accurately resolve solid-state exciton properties. Here, in this study, we show, for a range of semiconductors and insulators, that the commonly employed approach of uniformly sampling the Brillouin zone can lead to underconverged exciton binding energies, as impractical grid sizes are required to achieve adequate convergence. We further show that nonuniform sampling of the Brillouin zone, focused on the region of reciprocal space where the exciton wave function resides, enables efficient rapid numerical convergence of exciton binding energies at a given level of theory. We propose a well-defined convergence procedure, which can be carried out at relatively low computational cost and which in some cases leads to a correction of previous best theoretical estimates by almost a factor of 2, qualitatively changing the predicted exciton physics. These results call for the adoption of nonuniform sampling methods for ab initio GW-BSE calculations and for revisiting previously computed values for exciton binding energies of many systems.

36 MATERIALS SCIENCE↗
DOE OSTI2025

PAH101: A GW+BSE Dataset of 101 Polycyclic Aromatic Hydrocarbon (PAH) Molecular Crystals

Abstract The excited-state properties of molecular crystals are important for applications in organic electronic devices. TheGWapproximation and Bethe-Salpeter equation (GW+BSE) is the state-of-the-art method for calculating the excited-state properties of crystalline solids with periodic boundary conditions. We present the PAH101 dataset ofGW+BSE calculations for 101 molecular crystals of polycyclic aromatic hydrocarbons (PAHs) with up to ~500 atoms in the unit cell. To the best of our knowledge, this is the firstGW+BSE dataset for molecular crystals. The data records include theGWquasiparticle band structure, the fundamental band gap, the static dielectric constant, the first singlet exciton energy (optical gap), the first triplet exciton energy, the dielectric function, and optical absorption spectra for light polarized along the three lattice vectors. The dataset can be used to (i) discover materials with desired electronic/optical properties, (ii) identify correlations between DFT andGW+BSE quantities, and (iii) train machine learned models to help in materials discovery efforts.

Science & Technology - Other Topics↗
DOE OSTI2024

ComDMFT v.2.0: Fully self-consistent ab initio GW+EDMFT for the electronic structure of correlated quantum materials

ComDMFT is a parallel computational package designed to study the electronic structure of correlated quantum materials from first principles. Our approach is based on the combination of first-principles methods and dynamical mean field theories. In version 2.0, we implemented fully-diagrammatic GW+EDMFT from first-principles self-consistently. In this approach, correlated electrons are treated within full GW+EDMFT and the rest are treated within full-GW, seamlessly. Further, this implementation enables the electronic structure calculation of quantum materials with weak, intermediate, and strong electron correlation without prior knowledge of the degree of electron correlation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Static Subspace Approximation for Random Phase Approximation Correlation Energies: Applications to Materials for Catalysis and Electrochemistry

Modeling complex materials using high-fidelity, ab initio methods at low cost is a fundamental goal for quantum chemical software packages. The GW approximation and random phase approximation (RPA) provide a unified description of both electronic structure and total energies using the same physics in a many-body perturbative approach that can be more accurate than generalized-gradient density functional theory (DFT) methods. However, GW/RPA implementations have historically been limited to either specific materials classes or application toward small chemical systems. Here, the static subspace approximation allows for reduced cost full-frequency GW/RPA calculations and has previously been benchmarked thoroughly for GW calculations. Here, we describe our approach to including partial occupations of electronic orbitals in full-frequency GW and RPA calculations for the study of electrocatalysts. We benchmarked RPA total energy calculations using the subspace approximation across a diverse test suite of materials for a variety of computational parameters. The benchmarking quantifies the impact of different extrapolation procedures for representing the static polarizability at infinite screened cutoff, and shows that using screened cutoffs above 20-25 Ryd result in diminishing accuracy returns for predicting RPA total energies. Additionally, for moderately sized electrocatalytic models, 2-3 times fewer computational resources are used to compute RPA total energies by representing the static polarizability with 20-30% of the static subspace basis, with an error of approximately 0.01 eV or better in RPA adsorption energy calculations. Finally, we show that for these electrochemical models RPA can shift DFT adsorption energy shifts by up to 0.5 eV and that GW can frequently shift DFT eigenvalues of surface and adsorbate states by approximately 0.5-1 eV.

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

Time-dependent density functional theory with the orthogonal projector augmented wave method

The projector augmented wave (PAW) method of Blöchl linearly maps smooth pseudo wavefunctions to the highly oscillatory all-electron DFT orbitals. Compared to norm-conserving pseudopotentials (NCPP), PAW has the advantage of lower kinetic energy cutoffs and larger grid spacing at the cost of having to solve for non-orthogonal wavefunctions. We earlier developed orthogonal PAW (OPAW) to allow the use of PAW when orthogonal wavefunctions are required. In OPAW, the pseudo wavefunctions are transformed through the efficient application of powers of the PAW overlap operator with essentially no extra cost compared to NCPP methods. Previously, we applied OPAW to DFT. Here, we take the first step to make OPAW viable for post-DFT methods by implementing it in real-time time-dependent (TD) DFT. Using fourth-order Runge–Kutta for the time-propagation, we compare calculations of absorption spectra for various organic and biological molecules and show that very large grid spacings are sufficient, 0.6–0.7 bohr in OPAW-TDDFT rather than the 0.4–0.5 bohr used in traditional NCPP-TDDFT calculations. This reduces the memory and propagation costs by around a factor of 3. Our method would be directly applicable to any post-DFT methods that require time-dependent propagations such as the GW approximation and the Bethe–Salpeter equation.

Chemistry↗