SEARCH · Search NASA

Results for “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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 397 records · Page 22

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Urbanization exacerbates continental- to regional-scale warming

Urbanization is usually ignored when estimating past changes in large-scale climate and for future climate projections since cities historically covered a small fraction of the Earth’s surface. Here, by combining global land surface temperature observations with historical estimates of urban area, we demonstrate that the urban contribution to continental- to regional-scale warming has become non-negligible, especially for rapidly urbanizing regions and countries in Asia. Consequently, expected urban expansion over the next century suggests further increased urban influence on large-scale surface climate in the future (approximately 0.16 K for North America and Europe for high-emission scenario in 2100). Based on these results, also seen for air temperature, we argue that, in line with other forms of land use/land cover change, urbanization should be explicitly included in climate change assessments. This requires incorporation of dynamic urban extent and biophysics in current-generation Earth system models to quantify potential urban feedback on the climate system across scales.

Chakraborty, TC↗

DOE OSTI2025

Dynamics of McMillan mappings I. McMillan multipoles

In this article, we consider two dynamical systems: the McMillan sextupole and octupole integrable mappings, originally proposed by Edwin McMillan. Both represent the simplest symmetric McMillan maps, characterized by a single intrinsic parameter. While these systems find numerous applications across various domains of mathematics and physics, some of their dynamical properties remain unexplored. We aim to bridge this gap by providing a comprehensive description of all stable trajectories, including the parametrization of invariant curves, Poincaré rotation numbers, and canonical action–angle variables. In the second part, we establish connections between these maps and general chaotic maps in standard form. Our investigation reveals that the McMillan sextupole and octupole serve as first-order approximations of the dynamics around the fixed point, akin to the linear map and quadratic invariant (known as the Courant–Snyder invariant in accelerator physics), which represents zeroth-order approximations (referred to as linearization). Furthermore, we propose a novel formalism for nonlinear Twiss parameters, which accounts for the dependence of rotation number on amplitude. This stands in contrast to conventional betatron phase advance used in accelerator physics, which remains independent of amplitude. Notably, in the context of accelerator physics, this new formalism demonstrates its capability in predicting dynamical aperture around low-order resonances for flat beams, a critical aspect in beam injection/extraction scenarios.

43 PARTICLE ACCELERATORS↗

DOE OSTI2025

Fluorinated polythiophenes with ester side chains for boosting the V OC and efficiency in non-fullerene polymer solar cells

As the power conversion efficiency (PCE) of organic photovoltaics (OPVs) approximates the 19 % threshold, wide bandgap (WBG) polythiophenes (PTs) have gained increasing attention due to their simple chemical structures and ease of synthesis, making them promising candidates for large-scale production. However, the benchmark polymer poly(3-hexylthiophene) (P3HT) is limited by its high-lying HOMO energy level, which restricts the open-circuit voltage (V OC ) in solar cells. Here, in this study, we introduce a novel series of PT derivatives (PDC8-T, PDC8-3T, PDC16-3T, and PDC16-3T-2F) featuring ester side chains designed to fine-tune electronic properties through a streamlined three-step synthesis. Additionally, we incorporated a π-spacer to reduce steric hindrance and elongated alkyl side chains to improve solubility and processability. Compared to P3HT, these PT derivatives demonstrate a significant reduction in HOMO energy levels, lowering by approximately 0.3–0.4 eV. Among them, PDC16-3T-2F—with fluorine atom substitution—achieves the lowest HOMO energy level, induces a coplanar molecular conformation, and enhances polymer aggregation behavior. We evaluated these PT derivatives in inverted non-fullerene bulk-heterojunction (NFA BHJ) OPVs. The PDC8-3T device showed a relatively low PCE of 0.69 %, with a V OC of 0.76 V, a short-circuit current density (J SC ) of 3.32 mA/cm 2 , and a fill factor (FF) of 27.3 %. In contrast, the PDC16-3T-2F device achieved an impressive PCE of 7.21 %, with a V OC of 0.85 V, a J SC of 14.60 mA/cm 2 , and an FF of 58.4 %. This remarkable improvement is attributed to the fluorine substitution, which not only enhances molecular orientation but also downshifts the HOMO energy level and further boosts the V OC . Hence, these molecular design strategies have led to a fibrillar bicontinuous interpenetrating network with optimal nanoscale phase separation within the active layer.

14 SOLAR ENERGY↗

DOE OSTI2025

Assessment of flamelet/progress variable methods for supersonic combustion

Tabulated chemistry models, including the flamelet/progress variable approach, have been successfully used for a variety of turbulent flame simulations. The progress variable describes the progress of reactions in a system and parameterizes a lookup table of thermochemical variables. This approach reduces the cost of simulations, transporting only one scalar (progress variable) instead of the many species mass fractions required for detailed chemistry. Originally developed for low Mach number flame simulations, recent works have focused on extensions of this approach to compressible flames, supersonic combustion, and detonations, with applications such as scramjet combustors and rotating detonation engines. Unlike low Mach simulations, compressible flow simulations require solving the energy transport equation, which is coupled to the equation of state. This leads to additional modeling challenges regarding the thermodynamics and its impact on the chemistry. The validity of modeling assumptions, for example the relationship between energy and temperature, also varies with the combustion regime. The present work provides a detailed assessment of the existing strategies for chemistry tabulation for compressible/supersonic combustion, including detonations. A priori analysis indicates that approximations which are reasonable for weakly compressible flames may break down for shock-induced combustion. Furthermore, the analysis identifies specific assumptions and approximations that do not hold for detonations, emphasizing that care must be taken when applying tabulated chemistry models outside their intended combustion regimes.

Detonations↗

DOE OSTI2026

Directional dependence of equilibrated TLD-400 chips in multiple radiation fields

Fielding of thermoluminescent dosimeters (TLDs) for measurement of photon radiation dose in experiments is the standard practice for γ-irradiation facilities, pulsed power x-ray facilities, and reactor facilities at Sandia National Laboratories (SNL). Due to the high-dose experimental conditions and the mixed ( 1 n, γ) fields in these facilities, SNL radiation metrologists have historically used CaF 2 :Mn TLDs (also known as TLD-400). Recent inquiries to the radiation metrology staff have raised concerns that the aluminum-equilibrated TLD-400 chips may exhibit a directionally dependent response. The metrologists were asked whether the dose measured by the chip may be impacted by the angle of incidence on the equilibrated TLD. To provide a thorough answer to this query, a set of adjoint Monte Carlo radiation transport calculations was performed for three different equilibrated TLDs as well as bare TLD chips using the Integrated Tiger Series (ITS) code. The key feature of the adjoint calculations performed in this study is that each photon escaping the modeled geometry was tallied into angular bins to provide information on both the energy and angular dependence of the equilibrated chip. After the adjoint Monte Carlo calculations were completed, the resulting energy-dependent response function for each angular bin was convolved with multiple photon energy spectra representing various radiation facilities at SNL. The radiation facilities selected for analysis span a range of photon energies from approximately 1 keV up to approximately 20 MeV. Thus, the presented results are applicable to a wide variety of radiation facilities around the world. Although the bare TLD-400 chip was expected to display the largest variation due to the photon angle of impact on the dosimeter, the dosimeter with the thinnest aluminum equilibrator (SNL thin equilibrated TLD-400) was determined to have the biggest differential between the impact angle with the maximum dose ( D max ) and the impact angle with the minimum dose ( D min ). However, the SNL normal equilibrated TLD-400 chip demonstrated a dramatic reduction in that differential between maximum and minimum dose angles. The reduction in this differential is one of the dominant factors in experimenters’ choice to field these dosimeters at SNL radiation facilities. The results from the PNNL (Hanford) energy-flattening field capsule are consistent with the previous publications. The differential dose responses due to the photon impact angle indicate that experimenters should strive to field their TLD-400 dosimeters in a consistent manner to avoid additional uncertainty in the measurements based on dosimeter orientation.

Adjoint radiation transport↗

DOE OSTI2024

Minimizing grid energy consumption in wastewater treatment plants: Towards green energy solutions, water sustainability, and cleaner environment

Wastewater treatment plants (WWTPs) consume significant amount of energy to sustain their operation. From this point, the current study aims to enhance the capacity of these facilities to meet their energy needs by integrating renewable energy sources. The study focused on the investigation of two primary solar energy systems in As Samra WWTP in Jordan. The first system combines parabolic trough collectors (PTCs) with thermal energy storage (TES). This system primarily serves to fulfill the thermal energy demands of the plant by reducing the demands from boiler units, which allows more biogas for electricity generation. The second system is a photovoltaic (PV) system with Lithium-Ion batteries, which directly produces electricity that will be used to cover part of the electrical energy demands of plant. To assess the optimal configuration, two distinct scenarios have been formulated and compared to the current case scenario (SC#1). The first scenario focuses on maximizing the net present value (NPV) and minimizing the levelized cost of electricity (LCOE). The second scenario is centred on minimizing the levelized cost of heat (LCOH). The findings indicate that both scenarios succeeded in reducing the reliance on the grid to a value that reach 1 %. Moreover, they both reduced biogas percentage in energy production from 88 % to approximately 65 % through the integration of the PV system. In terms of thermal demand, SC#2 reduced the reliance on biogas boiler units from 100 % to 25 %, while SC#3 achieved an even more impressive reduction to just 8 %. The best LCOE value was attained in SC#2, at 0.0895 USD/kWh, with an NPV of 10.54 million USD. Conversely, SC# 3 yielded an LCOH value of 0.0432 USD/kWh th compared to 0.0534 USD/kWh th USD for SC#2. In conclusion, despite their relatively high capital and operating costs, SC#2 and SC#3 managed to substantially decrease the annual electricity expenditure from approximately 2 million USD to 86,000 USD and 0 USD, respectively.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2025

A predictive analytical model of electrical transport in multi-principal-element alloys

A predictive analytical model is presented for the electrical conductivity of multi-principal-element alloys (MPEAs), including those containing aluminum, transition metals, and refractory metals. Given that the lattice parameter of the Wigner-Seitz cell of an MPEA is similarly variable to a bulk metallic glass, it is postulated that electron scattering can be approximated by a series of two-level systems. Here, the resulting reduced-order model enabled an accurate determination of electrical resistivity and electron thermal conductivity based on the scattering of electrons in a two-level system across a Bloch-potential-based virtual crystal approximation. Model results are compared to experimental four-point probe electrical resistivity measurements between 300 K and 700 K for Al 0.3 CoCrCuFeNi, CoCrFeMnNi, (CoCrFeMnNi) 0.98 W 0.02 , (CoCrFeMnNi) 0.95 W 0.05 , and Nb 4 Ta 4 V 3 Ti, for model validation.

Analytical model↗

DOE OSTI2026

Suppression and enhancement of electromagnetic pulses from laser–target interactions by strong magnetic fields

Laser–target interactions generate intense electromagnetic pulses (EMPs) that can interfere with measurements and damage equipment. In this paper we show that applying a magnetic field to nanosecond pulse laser–target interactions decreases the magnitude of EMPs. We demonstrate this effect in two experiments with different geometries (spherical versus planar), laser intensities ( ∼10 13 versus ∼10 15 W/cm 2 ) and applied field strength (12 versus 0.1 T) that both observed suppression of EMPs in approximately the 1 GHz band (by factors of 0.65× and 0.32×, respectively). We then observe the opposite effect at high intensities with a picosecond pulse: for planar experiments with laser intensities of approximately 10 19 W/cm 2 and magnetic fields of 6–10 T, the magnitude of EMPs is increased by a factor of 1.75×. These results provide a benchmark for models of EMP generation, but suggest that magnetic fields are not a viable solution for mitigating EMPs in the high-intensity laser facilities where they are most damaging.

electromagnetic pulse↗

DOE OSTI2025

A Colebrook equation for impinging radial wall jets

Here, in this study, we evaluate the skin coefficient of friction for steady turbulent radial wall jets across smooth and rough surfaces. Although the Colebrook equation has been used successfully for many decades to evaluate friction factors for flows through smooth and rough pipes, how roughness affects the skin friction coefficient for steady turbulent radial wall jets remains unclear. Here we explore a Colebrook-type equation for skin friction coefficients associated with single-phase turbulent radial wall jets arising from orthogonally impinging circular jets. The fully iterative solution, based on well-established concepts of turbulent wall-bounded flow, is presented along with a power-law approximation and a non-iterative approximation for the friction coefficient derived therefrom. We find the skin coefficient of friction defined on the peak radial velocity to be a function of position over rough but not smooth surfaces in contrast to pipe friction factors that remain independent of axial position. These results follow expected trends, explain prior heterogeneity in power-law expressions for the skin friction coefficient and have significant implications for the industrial use of jets in mixing vessels.

friction losses↗

DOE OSTI2024

Latest Pleistocene glacial chronology and paleoclimate reconstruction for the East River watershed, Colorado, USA

Reconstructing Pleistocene glaciation timing and extent is vital for understanding paleoclimate. Whereas late Pleistocene glaciation has been studied extensively in western North American mountain ranges, the glacial history of the western Elk Range in Colorado remains understudied, particularly in the East River watershed, a site of intense scientific focus. Here we use cosmogenic nuclide exposure and depth–profile dating methods to determine the timing of glaciation in the East River watershed. We use glacier modeling to reconstruct paleoglacier extents and quantify past climate conditions. Our findings indicate that the East River glacier retreated from its maximum position approximately 17–18 ka, moving to recessional positions between 13 and 15 ka, before experiencing more substantial retreat to high-elevation cirques around 13 ka. Glacier modeling suggests that the maximum ice extents at 17–18 ka could have been sustained by temperature depressions of approximately −6.5°C compared to modern conditions, assuming consistent precipitation. Additionally, the ice position at 13–15 ka could have been supported by temperature depressions of around −4.0°C. These results offer insights into the deglaciation timeline in the East River watershed and broader western Elk Range as well as paleoclimate conditions during the late Pleistocene, which may aid future research on critical zone evolution in the East River watershed.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Coil optimization for quasi-helically symmetric stellarator configurations

Filament-based coil optimizations are performed for several quasi-helical stellarator configurations, beginning with the one from Landreman & Paul ( Phys. Rev. Lett. , vol. 128, 2022, 035001), demonstrating that precise quasi-helical symmetry can be achieved with realistic coils. Several constraints are placed on the shape and spacing of the coils, such as low curvature and sufficient plasma–coil distance for neutron shielding. The coils resulting from this optimization have a maximum curvature 0.8 times that of the coils of the Helically Symmetric eXperiment (HSX) and a mean squared curvature 0.4 times that of the HSX coils when scaled to the same plasma minor radius. When scaled up to reactor size and magnetic field strength, no fast particle losses were found in the free-boundary configuration when simulating 5000 alpha particles launched at $3.5\,\mathrm {MeV}$ on the flux surface with a normalized toroidal flux of $s=0.5$ . An analysis of the tolerance of the coils to manufacturing errors is performed using a Gaussian process model, and the coils are found to maintain low particle losses for smooth, large-scale errors up to amplitudes of approximately $0.15\,\mathrm {m}$ . Another coil optimization is performed for the Landreman–Paul configuration with the additional constraint that the coils are purely planar. Visual inspection of the Poincaré plot of the resulting magnetic field-lines reveal that the planar modular coils alone do a poor job of reproducing the target equilibrium. Additional non-planar coil optimizations are performed for the quasi-helical configuration with $5\,\%$ volume-averaged plasma beta from Landreman et al. ( Phys. Plasma , vol. 29, issue 8, 2022, 082501), and a similar configuration also optimized to satisfy the Mercier criterion. The finite beta configurations had larger fast-particle losses, with the free-boundary Mercier-optimized configuration performing the worst, losing approximately $5.5\,\%$ of alpha particles launched at $s=0.5$ .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Enhanced Collisional Losses from a Magnetic Mirror Using the Lenard-Bernstein Collision Operator

Collisions are crucial in governing particle and energy transport in plasmas confined in a magnetic mirror trap. Modern gyrokinetic codes model transport in magnetic mirrors, but some utilize approximate model collision operators. This study focuses on a Pastukhov-style method of images calculation of particle and energy confinement times using a Lenard-Bernstein model collision operator. Prior work on parallel particle and energy balances used a different Fokker-Planck plasma collision operator. The method must be extended in non-trivial ways to study the Lenard-Bernstein operator. To assess the effectiveness of our approach, we compare our results with a modern finite element solver. Our findings reveal that the particle confinement time scales like a exp( a 2 ) using the Lenard-Bernstein operator, in contrast to the more accurate scaling that the Coulomb collision operator would yield a 2 exp( a 2 ), where a 2 is approximately proportional to the ambipolar potential. We propose that codes solving for collisional losses in magnetic mirrors utilizing the Lenard-Bernstein or Dougherty collision operator scale their collision frequency of any electrostatically confined species. This study illuminates the collision operator’s intricate role in the Pastukhov-style method of images calculation of collisional confinement.

fusion plasma↗

DOE OSTI2026

Oil in Alaska North Slope Gas Hydrate Reservoir: Micro-CT and Flow Simulation Insights into Permeability

Gas hydrate-bearing sands on the Alaska North Slope (ANS) host minor volumes of crude oil whose impact on formation permeability has never been quantified relative to a gas hydrate reservoir system. Here, we combine in situ pressure-core microcomputed-tomography (μ-CT), thermogravimetric analysis (TGA), gas-chromatography–mass-spectrometry (GC–MS), scanning-electron microscopy coupled with energy-dispersive spectroscopy (SEM-EDS), and Stokes-flow simulation to (i) quantify oil saturation, (ii) infer its source, and (iii) evaluate its influence on permeability after gas-hydrate dissociation. μ-CT and SEM-EDS imaging identify the host sediment as silty, with some grain-coating clay present. Eight representative 3003-voxel subvolumes extracted from a preserved pressure core (Core 15P-3c, from the B1 sand (unit B), HYDRATE 02 Geo Data Well (GDW), 886.96–887.29 m measured depth, MD) exhibit porosities of 0.34–0.37 and an average oil saturation of 0.09 ± 0.03, in comparison to an independent TGA-based estimate of 0.16. GC-MS fingerprints obtained from an adjacent depressurized core (Core 17P-1, 891.37–891.44 m MD), together with oil pore habit revealed by μ-CT, suggest that the oil is partly native but also includes contributions from mineral oil-based drilling fluid contamination. Further analysis of oil–porewater interactions and flushing tests confirms that native oil saturation in the analyzed samples ranges from 0.04 to 0.08 with an average of approximately 0.06, and the oil is highly immobile. Flow simulations demonstrate that the native oil saturation of ≈0.06 reduces permeability by approximately 50% due to both pore blockage and increased flow-path tortuosity. These findings reveal the potential presence of native oil within the ANS gas hydrate reservoir and suggest that future simulation models may need to account for its impact on permeability to improve long-term performance predictions of gas and water production.

02 PETROLEUM↗

DOE OSTI2024

How Does HF-DFT Achieve Chemical Accuracy for Water Clusters?

Bolstered by recent calculations of exact functional-driven errors (FEs) and density-driven errors (DEs) of semilocal density functionals in the water dimer binding energy , we investigate approximate FEs and DEs in neutral water clusters containing up to 20 monomers, charged water clusters, and alkali- and halide-water clusters. Our proxy for the exact density is r 2 SCAN 50, a 50% global hybrid of exact exchange with r 2 SCAN, which may be less correct than r 2 SCAN for the compact water monomer but importantly more correct for long-range electron transfers in the noncompact water clusters. We show that SCAN makes substantially larger FEs for neutral water clusters than r 2 SCAN, while both make essentially the same DEs. Unlike the case for barrier heights, these FEs are small in a relative sense and become large in an absolute sense only due to an increase in cluster size. SCAN@HF, short for SCAN evaluated on the Hartree–Fock (HF) density, produces a cancellation of errors that makes it chemically accurate for predicting the absolute binding energies of water clusters. Likewise, adding a long-range dispersion correction to r 2 SCAN@HF, as in the composite method HF-r 2 SCAN-DC4, makes its FE more negative than in r 2 SCAN@HF, permitting a near-perfect cancellation of FE and DE. r 2 SCAN by itself (and even more so, r 2 SCAN evaluated on the r 2 SCAN 50 density), is almost perfect for the energy differences between water hexamers, and thus probably also for liquid water away from the boiling point. Thus, the accuracy of composite methods like SCAN@HF and HF-r 2 SCAN-DC4 is not due to the HF density being closer to the exact density, but to a compensation of errors from its greater degree of localization. We also give an argument for the approximate reliability of this unconventional error cancellation for diverse molecular properties. Lastly, we confirm this unconventional error cancellation for the SCAN description of the water trimer via Kohn–Sham inversion of the CCSD(T) density.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Liquid Phase Modeling in Porous Media: Adsorption of Methanol and Ethanol in H-MFI in Condensed Water

Zeolites are used in the chemical and separation industries for their exceptional selectivity, adsorption capacity, regenerability, and stability in gas and liquid phase processing. Here, we developed an explicit solvation method for predicting solvent/condensed phase effects on adsorption free energies in microporous media such as zeolites based on the hybrid quantum mechanical/molecular mechanical free energy perturbation (QM/MM-FEP) technique. Our explicit solvation method for zeolite systems, called eSZS, aims to capture site-specific interactions during the adsorption process at the Brønsted acid sites of H-MFI zeolite while still considering the diverse configuration space of the solvent molecules. This strategy is ideal for chemical reactions or adsorbates that interact with the microporous medium in few distinct adsorbate/transition state configurations, i.e., the harmonic or similar approximations are acceptable for the adsorbate/transition state while such approximations break down for the solvent molecules that require extensive configuration space sampling. In this way, our approach effectively overcomes the limitations of implicit solvation models and classical force field methods for describing solvation effects on chemical reactions within porous materials such as zeolites. Specifically, in this study, we investigated various aspects of our hybrid QM/MM approach, including QM cluster size dependencies in a periodic electrostatically embedded cluster model (PEECM), rules for link atoms at the QM/MM boundary, and functional and basis set considerations for converged and reasonably accurate gas and aqueous phase methanol and ethanol adsorption free energy predictions in H-MFI. For gas phase adsorption of methanol and ethanol in H-MFI at a Brønsted acid site in T12 position, we compute adsorption free energies at 298 K of −0.61 and −0.75 eV, respectively, using a PEECM containing 50 Si and 1 Al atom with ωB97x-D/def2-TZVP level of theory. For solvent effect calculations, we sample the aqueous phase using grand canonical Monte Carlo (GCMC) simulations to (1) obtain a mean field of electrostatic interactions in the reaction system and (2) perform a rigorous free energy perturbation calculation. Similar to the experimentally and computationally observed endergonic solvation effects observed for hydrocarbon adsorption on metal surfaces, we also observe that a condensed aqueous environment destabilizes methanol and ethanol at these acid sites in H-MFI at 298 K. Specifically, the computed solvation free energies of adsorption (ΔΔG solv ) for methanol and ethanol are +0.44 and +0.54 eV, respectively. From this study, it is evident that adsorbates (methanol and ethanol) are competing with water for adsorption space inside the H-MFI zeolite, leading to an endergonic solvation effect. Here, we expect that the endergonic, aqueous solvent effect during adsorption in microporous zeolites is highly tunable by changing the pore size and hydrophobicity of the microporous material as this will affect the water density inside the pore structure.

Adsorption↗

DOE OSTI2025

A Perspective on Quantum Computing Applications in Quantum Chemistry Using 25-100 Logical Qubits

The intersection of quantum computing and quantum chemistry represents a promising frontier for achieving quantum utility in domains of both scientific and societal relevance. Owing to the exponential growth of classical resource requirements for simulating quantum systems, quantum chemistry has long been recognized as a natural candidate for quantum computation. This perspective focuses on identifying scientifically meaningful use cases where early fault-tolerant quantum computers, which are considered to be equipped with approximately 25-100 logical qubits, could deliver tangible impact. While recent advances in classical computing have pushed the boundaries of tractable simulations to unprecedented scales, this logical-qubit regime represents the first window where quantum devices can pursue qualitatively distinct strategies, such as polynomial-scaling phase estimation, direct simulation of quantum dynamics, and active-space embedding, that remain challenging for classical solvers, such as multireference charge-transfer and conical-intersection states central to photochemistry and materials design. We highlight near-term opportunities in algorithm and software design, discuss representative chemical problems suited for quantum acceleration, and propose strategic roadmaps and collaborative pathways for advancing practical quantum utility in quantum chemistry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Reducing the Cost of Energy Differences in Variational Monte Carlo with Spotlight Sampling

Here, we investigate an approximate sampling scheme that can significantly reduce the cost scaling of variational Monte Carlo when it is employed to predict the energy differences associated with local chemical changes. Inspired by side-chaining and embedding methods, this spotlight sampling approach adopts an approximate fragmented Hamiltonian and correlated sampling to reduce cost scaling to the point that it is essentially linear with system size, with the potential to go sublinear if certain conditions are met. In tests on bond stretching energies in alcohols, hydrogen dimer chains, and molecules with various degrees of π-system delocalization, we observe the anticipated linear scaling and an explicit cost crossover with standard variational Monte Carlo.

Bumann, Sonja [University of California, Berkeley,↗

DOE OSTI2025

Overset-Grid Method with Smooth Orbital Partitioning for Molecular Scattering Calculations

To solve molecular photoionization and electron scattering problems, we use an overset-grid representation of electronic continuum functions, which has an extended central spherical grid that overlaps small spherical grids (subgrids) centered on each atom of a polyatomic molecule. Here, in this work, we present an improved algorithm that smoothly partitions the total wave function between the central grid and the atomic subgrids. The smooth partitioning allows one to use approximately one-fourth the number of partial waves on the central grid compared to our previous implementation with switching functions. The resulting numerical method for treating electron scattering and photoionization of polyatomic molecules combines the accuracy and flexibility of pure numerical grid representations with the rapid convergence of hybrid combinations of atom-centered basis-set expansions and grid methods. The overset-grid representation is implemented using the complex Kohn variational principle for scattering and photoionization amplitudes. The faster convergence with respect to the number of central grid partial waves is demonstrated and accuracy is verified by comparisons with the previous implementation and with far more computationally demanding single-center numerical expansions in electron-molecule scattering and photoionization calculations on the neon dimer (Ne 2 ) system, carbon tetrafluoride (CF 4 ) molecule, and the pyridine (C 5 H 5 N) molecule in the static-exchange approximation.

Molecules↗