Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “GW approximation”

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
All sourcesNASA NTRSDOE OSTI
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 OSTI2026

Tunable Electronic Energy Level Alignment and Exciton Diversity in Organic–Inorganic van der Waals Heterostructures

van der Waals stacking of two-dimensional (2D) materials offers a powerful platform for engineering material interfaces with tailored electronic and optical properties. While most van der Waals multilayers have featured inorganic monolayers, incorporating molecular monolayers introduces additional degrees of tunability and functionality. Here, in this study, we investigate hybrid bilayers composed of atomically thin perylene-based molecular crystals interfaced with monolayer transition metal dichalcogenides (TMDs), specifically MoS 2 and WS 2 . Using the ab initio many-body perturbation theory within the GW approximation and the Bethe-Salpeter equation approach, we predict emergent properties beyond those of the isolated constituent systems. Notably, we find substantial renormalization of monolayer molecular crystal band gap due to TMD-induced polarization. Furthermore, by varying the TMD monolayer, we demonstrate tuning of the energy level alignment of the bilayer and subsequent control over a diversity of lowest-energy excitons, which include strongly bound hybrid excitons and long-lived charge-transfer excitons. These findings establish organic-inorganic van der Waals heterostructures as a promising class of materials for tunable optoelectronic devices and quantum excitonic phenomena, expanding the design space for low-dimensional systems.

GW-BSE calculations↗
DOE OSTI2024

Electronic Properties of Ultra‐Wide Bandgap B x Al 1− x N Computed from First‐Principles Simulations

Abstract Ultra‐wide bandgap (UWBG) materials such as AlN and BN hold great promise for future power electronics due to their exceptional properties. They exhibit large bandgaps, high breakdown fields, high thermal conductivity, and high mechanical strengths. AlN and BN have been extensively researched, however, their alloys, B x Al 1− x N, are much less studied despite their ability to offer tunable properties by adjusting x . In this article, the electronic properties of 17 recently predicted ground states of B x Al 1− x N in the x = 0 − 1 range are predicted using first‐principles density functional theory and many‐body perturbation theory within GW approximation. All the B x Al 1− x N structures are found to be UWBG materials and have bandgaps that vary linearly from that of wurtzite‐phase ( w ) AlN (6.19 eV) to that of w ‐BN (7.47 eV). The bandstructures of B x Al 1− x N show that a direct‐to‐indirect bandgap crossover occurs near x = 0.25. Furthermore, it is found that B x Al 1− x N alloys have much larger dielectric constants than the constituent bulk materials (AlN = 9.3 ɛ 0 or BN = 7.3 ɛ 0 ), with values reaching as high as 12.1 ɛ 0 . These alloys are found to exhibit large dielectric breakdown fields in the range 9–35 MV cm −1 with a linear dependence on x . This work provides the much needed advancement in the understanding of the properties of B x Al 1− x N to aid their application in next‐generation devices.

Milne, Cody L.↗
DOE OSTI2026

Stochastic GW -GPU: Rapid Quasi-Particle Energies for Molecules beyond 10,000 Atoms

StochasticGW is a code for computing accurate quasi-particle (QP) energies of molecules and material systems in the GW approximation. StochasticGW utilizes the stochastic Resolution of the Identity (sROI) technique to enable a massively parallel implementation with computational costs that scale semilinearly with system size, allowing the method to access systems with tens of thousands of electrons. Here, we introduce a new implementation, StochasticGW-GPU, for which the main bottleneck steps have been ported to GPUs and give substantial performance improvements over previous versions of the code. We showcase the new code by computing band gaps of hydrogenated silicon clusters (Si x H y ) containing up to 10,001 atoms and 35,144 electrons, and we obtain individual QP energies with a statistical precision of better than ±0.03 eV with times-to-solution of less than 1 h.

Thomas, Phillip S. [Lawrence Berkeley National Lab↗
DOE OSTI2025

Near-half-metallic state in the half-Heusler PtMnSb film on a III-V substrate

The interplay between half-metallic ferromagnetism and spin-orbit coupling within the inversion-symmetry-broken structure of half Heuslers provides an ideal platform for various spintronics functionalities. Taking advantage of good lattice matching, it is highly desired to epitaxially integrate promising Heuslers into III-V semiconductor-based devices. PtMnSb is one of the first half Heuslers predicted to be an above-room-temperature half-metal with large spin orbit coupling, however, its half-metallicity and potential as a spintronics material has remained elusive due to lack of high quality samples. Here, in this study, we demonstrate epitaxial growth of single crystal PtMnSb(001) film on GaSb(001) substrates using molecular beam epitaxy. Direct observation of the band structure via angle-resolved photoemission spectroscopy and many-body perturbation theory within the quasi-particle self-consistent GW approximation (QPGW) reveal that PtMnSb hosts rather a near-half-metallic state with both spin bands crossing the Fermi level and with high spin polarization over 90%. Temperature dependence of magnetization also shows an anomalous enhancement below 60 K, which can be associated with the development of such a near-half-metallic state at low temperatures. Epitaxial growth of high crystalline PtMnSb on a III-V paves the way for systematic clarification of its spin transport properties with fine-tuning of strain in heterostructure devices.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Queued Up: 2025 Edition – Characteristics of Power Plants Seeking Transmission Interconnection As of the End of 2024 [Slides]

Electric transmission system operators (ISOs, RTOs, or utilities) require proposed power plants seeking to connect to the transmission grid to undergo a series of impact studies before they can be built. This process establishes what new transmission equipment or upgrades may be needed before a project can connect to the system and assigns the costs of that equipment. The lists of projects in this process are known as “interconnection queues”. In collaboration with interconnection.fyi, Berkeley Lab compiled, aggregated, and cleaned interconnection queue data from >50 transmission grid operators (7 ISO/RTOs and 49 non-ISO balancing areas), which collectively represent ~97% of currently installed U.S. electric generating capacity. The dataset includes requests submitted to queues through the end of 2024, and only includes requests seeking to connect to the transmission grid (not distribution-connected or behind-the-meter projects). The files below include both a PDF report and an Excel data file. The PDF report analyzes interconnection data and metrics through the end of 2024. The Excel data file includes (a) the full project-level interconnection queue dataset through 2024, (b) a codebook (data dictionary) describing each data field, and (c) 35 additional tabs featuring tables summarizing a range of interconnection metrics. Key highlights from the Queued Up: 2025 Edition (featuring data through 2024) include: • As of the end of 2024, there were ~10,300 projects actively seeking grid interconnection in the U.S., representing 1,400 GW of generation and approximately 890 GW of storage. • Historic withdrawal rates alongside relatively fewer new requests resulted in a 12% decrease in total active queue volume compared to the prior year. • Active natural gas capacity (136 GW, +72% year-over-year) increased in 2024, while solar (956 GW, -12%), storage (890 GW, -13%), and wind (271 GW, -26%) capacity decreased. • 408 GW of capacity already has a draft or executed interconnection agreement (IA) but has not yet reached commercial operations. • The time projects spend in queues before reaching COD is increasing. For the regions with available data, the median duration from IR to COD has doubled from <2 years for projects built in 2000-2007 to over 4 years for those built in 2018-2024. • Ultimately, most of this proposed capacity will not be built. Only 13% of capacity that submitted interconnection requests from 2000-2019 had reached commercial operations by the end of 2024; 77% of that capacity had been withdrawn and 10% was still active. • FERC Order 2023 and various other reforms are being implemented. These are important measures to reduce interconnection bottlenecks and enhance grid system reliability, but it is too early to measure and assess their full impact. • New additions for the 2025 edition include: (a) additional detail on data processing and gaps; (b) updates on interconnection reforms; (c) new analysis on interconnection agreements, and more.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

AquaPV: Regulatory and Environmental Considerations for Floating Photovoltaic Projects Located on Federally Controlled Reservoirs in the United States

To meet the nation's decarbonization goals, the U.S. Department of Energy's Solar Futures study forecasts that installed solar photovoltaic (PV) capacity must increase nearly tenfold, from 80 gigawatts (GW) in 2020 to approximately 760 GW cumulative installed capacity by 2035. Ground-mounted PV is expected to dominate future solar deployment and will require more than 3.5 million acres of land to meet annual demand projections (of nearly 45 GW) by 2030. However, various competing demands for land (e.g., agricultural production, conservation) and high land acquisition costs in specific locations could be challenges to meeting future PV demand solely with ground-mounted PV deployment. Floating photovoltaics (FPV) may be an alternative in locations where ground-mounted PV is not feasible and aid in reaching the nation's PV deployment and decarbonization goals. FPV is a newer siting approach in which a PV array is affixed to a floating apparatus and sited on a water body like a reservoir behind a dam. FPV systems may be stand-alone or co-located at new or existing hydroelectric facilities or pumped storage hydropower (PSH) facility reservoirs. Co-located FPV systems may or may not be operationally paired and work in tandem with the hydroelectric or PSH facility. This report provides novel analysis to understand the opportunities and challenges associated with developing stand-alone and co-located FPV projects on reservoirs in the United States. Specifically, the report explores potential environmental and energy benefits and environmental impacts associated with the siting, construction, and operation of FPV projects. The report also identifies and analyzes U.S. federal- and state-issued permits and authorizations required by federal laws to understand the licensing pathways and regulatory requirements for FPV projects sited on reservoirs licensed by the Federal Energy Regulatory Commission and on powered and non-powered reservoirs owned by the Bureau of Reclamation or U.S. Army Corps of Engineers.

ENERGY PLANNING, POLICY, AND ECONOMY,SOLAR ENERGY↗
DOE OSTI2025

Land of Opportunity: Potential for Renewable Energy on Federal Lands

Renewable energy (RE) in the United States has historically been deployed primarily on private lands, but the growing interest in RE development, generally across the country and specifically on federal lands, raises questions about the potential for RE on public lands. This study seeks to estimate RE technical potential on federal lands and project the amount of RE to be developed on federal lands under decarbonization scenarios for the contiguous United States. The study applied a combination of high-resolution geospatial analysis and power sector modeling and relied on multiple partner federal agencies and land administrators from the Bureau of Land Management, the U.S. Fish and Wildlife Service, the U.S. Forest Service, the U.S. Department of Defense, and the U.S. Department of Energy. In our reference siting access case, we estimate 44 million acres of federal land across the contiguous United States is potentially suitable for UPV development, which corresponds to a capacity technical potential of 5,750 gigawatts (GW) (Figure ES-1). Federal land area available for wind development is similar to UPV but wind’s generating capacity technical potential is lower (875 GW). The technical potential is estimated for two geothermal technologies, hydrothermal and enhanced geothermal systems (EGS), both of which have a smaller amount of federal land available for development (12 and 27 million acres, respectively). However, in terms of capacity, the technical potential for EGS (975 GW) is approximately equal to wind’s technical potential and there is an estimated 130 GW of hydrothermal potential. We also developed cases with more-limited land available for UPV and wind (for federal and non-federal lands) resulting in 96% reduction of wind capacity potential and 70% reduction for UPV. In a case with additional constraints applied to non-federal lands only (Limited Private), the overall (federal and non-federal) technical potential declines but the share of that technical potential on federal lands is higher than in the other siting cases. The technical potential is the maximum amount that could be developed, but only a small fraction would be developed or needed in the future. Across seven scenarios that achieve 100% carbon-free electricity by 2035, we estimate 26 GW to 270 GW of RE capacity could be deployed on federal lands by 2035. The three central scenarios have 51–84 GW of RE deployed by 2035 and requiring about 500,000 to 1,000,000 acres of total land area. Direct land consumption is less than total land use required, thus enabling co-use opportunities. The large technical potential estimates and the increasing deployment projections from the collection of scenarios show the opportunities for RE development on federal lands. Capturing these opportunities-while minimizing conflicts with other land uses, federal department or agency missions, and public interest, and simultaneously maximizing the economic, grid, and social value of the projects-would require collaborative planning among federal land administrators, grid planners, project developers, the public, and other stakeholders.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
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 OSTI2026

The past, present, and future of peaking thermal power plants in the United States

Today's power systems rely on "peaker plants" to reliably serve load during peak demand periods. In this study we consider the present competitiveness of different peaking options, how much and what kinds of plants have provided U.S. peaking capacity, and potential future peaking fleet compositions. We explore how capital intensity impacts the breakeven capacity factor between two potential resources: combustion turbines (CTs) and combined cycle plants (CCs). CTs outcompete CCs below 12%-17% annual capacity factor at today's prices, but can shift with changes to fuel or start costs. Historically, gas CTs and petroleum steam plants most closely fit the role of peakers. Peaking capacity could grow from approximately 280 GW today to 460-770 GW in 2050 composed of a wider range of resources. We conclude by discussing implications of this shift, with a focus on the potential planning considerations for shifting to a greater utilization of CCs for peaking needs.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2026

Impacts of Solar Land Availability on the Evolution of the U.S. Power System

Long-term planning scenarios of the United States electricity system have not identified solar resource availability to be a limiting factor in how much solar capacity might be deployed. However, there are some smaller regions where siting challenges indicate solar deployment could be limited by solar resource availability. In this work, we use geospatial modeling tools to develop a restrictive siting scenario that significantly reduces the amount and location of utility-scale solar resources. We use that restrictive siting regime, along with other previously defined siting scenarios, to examine how solar resource availability impacts the buildout of the electricity system using business-as-usual and decarbonization futures in a national-scale capacity expansion model. The restrictive siting scenario uses observed siting pressures to exclude many types of land on which development is allowed today and results in a greater than 90% reduction in solar resource potential-to approximately 3000 GW. Even with this drastic reduction, utility-scale solar capacity grows by more than 3 times through 2050 in the business-as-usual scenarios and by more than 6 times in the decarbonization scenarios, though less solar is deployed in scenarios with lower available solar resource. We quantify how system costs and generation mixes change as solar resource availability changes.

14 SOLAR ENERGY↗
DOE OSTI2025

Multiple Slater Determinants and Strong Spin-Fluctuations as Key Ingredients of the Electronic Structure of Electron- and Hole-Doped Pb 10− x Cu x (PO 4 ) 6 O

LK-99, with chemical formula Pb 10−x Cu x (PO 4 ) 6 O, was recently reported to be a room-temperature superconductor. While this claim has met with little support in a flurry of ensuing work, a variety of calculations (mostly based on density-functional theory) have demonstrated that the system possesses some unusual characteristics in the electronic structure, in particular flat bands. We have established previously that within DFT, the system is insulating with many characteristics resembling the classic cuprates, provided the structure is not constrained to the P3(143) symmetry nominally assigned to it. Here we describe the basic electronic structure of LK-99 within self-consistent many-body perturbative approach, quasiparticle self-consistent GW (QSGW) approximation and their diagrammatic extensions. QSGW predicts that pristine LK-99 is indeed a Mott/charge transfer insulator, with a bandgap gap in excess of 3 eV, whether or not constrained to the P3(143) symmetry. When Pb 9 Cu(PO 4 ) 6 O is hole-doped, the valence bands modify only slightly, and a hole pocket appears. However, two solutions emerge: a high-moment solution with the Cu local moment aligned parallel to neighbors, and a low-moment solution with Cu aligned antiparallel to its environment. In the electron-doped case the conduction band structure changes significantly: states of mostly Pb character merge with the formerly dispersionless Cu d state, and high-spin and low spin solutions once again appear. Thus we conclude that with suitable doping, the ground state of the system is not adequately described by a band picture, and that strong correlations are likely. Irrespective of whether this system class hosts superconductivity or not, the transition of Pb 10 (PO 4 ) 6 O from being a band insulator to Pb 9 Cu(PO 4 ) 6 O, a Mott insulator, and multi-determinantal nature of doped Mott physics make this an extremely interesting case-study for strongly correlated many-body physics.

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

Assessing the Costs of Hydropower at Non-Powered Dams Using a Reference Site Model

Hydropower capacity in the United States currently stands at approximately 103 GW, and there are remaining water resources that could help meet the rapidly increasing demand for electricity and ancillary grid services. Existing dams that do not generate power, known as non-powered dams, are a near-term solution to enhance the contribution of hydropower to the US power grid. However, there are thousands of such sites, and the lack of detailed information about their specific characteristics and associated costs presents significant challenges for stakeholders. This study addresses these challenges by developing a reference site model to evaluate the potential costs of hydropower at non-powered dams using currently available technologies. An application of the model reveals a wide range of estimates for capacity, capital costs, levelized cost of electricity, and cost components. While many sites are competitive with current technologies, the majority would need cost-reducing innovations to be viable. Despite the limited available information, the model offers valuable insights into the relative competitiveness of hydropower projects at non-powered dams. The simulation results highlight the need for continued technological advancements in hydropower and provide a basis for evaluating the benefits of new innovations.

cost drivers↗
DOE OSTI2025

E-Area Low-Level Waste Facility Inadvertent Human Intruder Limits and Doses in Support of the PA2022

This report documents the inadvertent human intruder (IHI) analysis for the E-Area Low-Level Waste Facility (ELLWF) at the Savannah River Site (SRS), near Aiken, South Carolina. This analysis supports the revised ELLWF Performance Assessment (PA), complying with the Department of Energy standard for operation of low-level waste disposal facilities (USDOE, 2017). The ELLWF is an operating waste disposal facility and is scheduled to continue accepting waste to 2065. One task of the revised PA is to establish waste inventory limits for the various disposal units at ELLWF. This is done by modeling future contaminant release and transport through applicable pathways to human receptors, comparing predicted doses per disposed curie with applicable performance measures, to obtain inventory limits which will assure that doses to receptors do not exceed performance measures. This report documents results of modeling future doses to one class of receptor, the inadvertent human intruder. It is assumed that after site closure, public knowledge of the site is lost, and IHIs will engage in activities on the ELLWF that will disrupt the closure cap, causing dose to the IHI. Following USDOE (2017), six different stylized exposure scenarios are considered, simulating activities by an IHI which could result in a radiological dose. The six scenarios are: • Acute – Basement Construction: IHI constructs a basement and encounters waste during excavation which is inadvertently mixed with clean soil and diluted. • Acute – Well Drilling: IHI drills a water well through waste and is exposed to drill cuttings mixed with clean soil that are brought to the surface. • Acute – Discovery: IHI begins constructing a basement but stops when encountering the riprap in the final closure cap and is exposed to photon radiation from unexcavated material residing in the undisturbed waste zone. • Chronic – Agriculture: Resident IHI is exposed to waste that was excavated for basement construction and mixed with native soil in the intruder’s vegetable garden. • Chronic – Post-Drilling: Resident IHI is exposed to waste from drill cuttings mixed with native soil and scattered in the garden area. • Chronic – Residential: Resident IHI is exposed to external radiation while in home located above waste with shielding provided by the concrete basement floor and any soil or engineered material remaining between the basement and waste. Dose calculations are performed using the SRNL Dose Toolkit (Aleman, 2023), following the approach of Smith et al (2019). Calculations are performed separately for 27 of the 33 disposal units (DUs) at ELLWF and are radionuclide specific. The results of the IHI analysis include: • Dose Factors: mrem per disposed curie (acute) and mrem/yr per disposed curie (chronic) for each parent radionuclide, for each DU. • Inventory Limits: in curies, for each parent radionuclide, for each DU. • Estimated Dose to IHI: mrem (acute) and mrem/yr (chronic), for each DU, given its projected closure inventory without inventory biases applied. Most DU-specific IHI inventory limits are in the range of 10 3 to 10 7 curies per nuclide. The lowest inventory limits are associated with gamma-emitters such as Sn-126, Ra-226, Th-232, and Cm-248. Radionuclides with short half-lives such as Pu-241, and nuclides which are pure beta emitters or which decay by electron capture, such as Ni-59 and Ni-63, have the highest limits. For the 27 evaluated DUs, predicted IHI doses are shown in Table ES-1. The maximum acute dose is 1.18 mrem, at ST23, much less than the DOE performance measure of 500 mrem (USDOE, 2017). The highest chronic dose is 37.2 mrem/yr at ST02, below the DOE performance measure of 100 mrem/yr. Also shown are estimated inventory sums of fractions (SOFs) at closure in 2065, for groundwater (GW) and IHI pathways. For each DU, the inventory is constrained by the GW pathway. For most DUs, the IHI SOFs are approximately 1000 times lower than the GW SOF values, and the IHI pathway does not drive risk for any disposal unit.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

Energy storage in combined gas-electric energy transitions models: The case of California

California’s vision for a net-zero future by 2045 relies heavily on variable renewable energy systems. Thus, energy storage - particularly long-duration storage - could play a fundamental role in reliably supplying low-carbon electricity. We study energy storage using the BRIDGES model, a combined gas-electric capacity expansion model for California across multiple investment periods (2025-2045), modeled with progressively decreasing carbon emission targets to a zero emissions by 2045. This least-cost optimization model includes renewable gas production via power-to-gas, long-term storage of energy in gaseous form, electric energy storage such as through batteries and hydrogen storage, and renewable energy generation, all with capacity tracking and investment. Multiple scenarios are evaluated to examine the sensitivity of the optimal storage portfolio to system-level and sector-level parameters. The scenario results show that all electric energy storage systems - which vary in storage duration - are deployed and required in a net-zero California in 2045, amounting to around 75 GW of storage capacity. Lithium ion systems make up approximately 80% of this power capacity and supply most short-run storage needs. Hydrogen storage - in the form of a power-to-gas-to-power system - emerges as a replacement to conventional natural gas storage, comprising most of the total energy storage capacity (~ 4 TWh). This capacity is less than 5% of the current natural gas storage capacity (94 TWh), indicating sufficient room for repurposing part of the gas infrastructure. A demand-side sensitivity analysis proves that higher electricity demand correlates with more builds of Li-ion batteries, while higher industrial heat demand leads to more builds of long-duration storage systems in a net-zero economy. Furthermore, power-to-gas systems satisfy part of the industrial heat demand by locally supplying renewable gas, which overtakes the traditional centralized gas storage and transfers through pipelines, casting significant doubts on the future of the large-scale gas infrastructure.

03 NATURAL GAS↗
DOE OSTI2025

The discriminant power of bubble wall velocities: gravitational waves and electroweak baryogenesis

A precise determination of the bubble wall velocity v$_{w}$ is crucial for making accurate predictions of the baryon asymmetry and gravitational wave (GW) signals in models of electroweak baryogenesis (EWBG). Working in the local thermal equilibrium approximation, we exploit entropy conservation to present efficient algorithms for computing v$_{w}$, significantly streamlining the calculation. We then explore the parameter dependencies of v$_{w}$, focusing on two sample models capable of enabling a strong first-order electroweak phase transition: a ℤ$_{2}$-symmetric singlet extension of the SM, and a model for baryogenesis with CP violation in the dark sector. We study correlations among v$_{w}$ and the two common measures of phase transition strength, α$_{n}$ and v$_{n}$/T$_{n}$. Interestingly, we find a relatively model-insensitive relationship between v$_{n}$/T$_{n}$ and α$_{n}$. We also observe an upper bound on α$_{n}$ for the deflagration/hybrid wall profiles naturally compatible with EWBG, the exact value for which varies between models, significantly impacting the strength of the GW signals. In summary, our work provides a framework for exploring the feasibility of EWBG models in light of future GW signals.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Range-separated thermal hybrid exchange-correlation density functional for accurate band-gap calculations of warm dense matter

Here, in this study, we present a range-separated thermal hybrid exchange-correlation (XC) functional, termed RS-KDT0, which provides highly accurate band-gap predictions at elevated temperature while reducing the computational cost relative to thermal global hybrid functionals. At zero temperature, RS-KDT0 recovers the widely used ground-state range-separated hybrid functional HSE06. To assess its performance across temperature regimes, we conducted static calculations of electronic pressure, band gaps, and band structures. RS-KDT0 exhibits the same level of accuracy as HSE06 in the low-temperature regime and demonstrates improved agreement with finite-temperature GW calculations at elevated temperatures when compared to other thermal XC functionals such as the generalized gradient approximation (GGA)-level KDT16 and the global hybrid KDT0.

density functional theory↗
DOE OSTI2024

GPU-Accelerated Solution of the Bethe–Salpeter Equation for Large and Heterogeneous Systems

We present a massively parallel GPU-accelerated implementation of the Bethe–Salpeter equation (BSE) for the calculation of the vertical excitation energies (VEEs) and optical absorption spectra of condensed and molecular systems, starting from single-particle eigenvalues and eigenvectors obtained with density functional theory. The algorithms adopted here circumvent the slowly converging sums over empty and occupied states and the inversion of large dielectric matrices through a density matrix perturbation theory approach and a low-rank decomposition of the screened Coulomb interaction, respectively. Further computational savings are achieved by exploiting the nearsightedness of the density matrix of semiconductors and insulators to reduce the number of screened Coulomb integrals. We scale our calculations to thousands of GPUs with a hierarchical loop and data distribution strategy. The efficacy of our method is demonstrated by computing the VEEs of several spin defects in wide-band-gap materials, showing that supercells with up to 1000 atoms are necessary to obtain converged results. We discuss the validity of the common approximation that solves the BSE with truncated sums over empty and occupied states. In conclusion, we then apply our GW-BSE implementation to a diamond lattice with 1727 atoms to study the symmetry breaking of triplet states caused by the interaction of a point defect with an extended line defect.

Absorption spectra↗