Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nuclear Equation of State”

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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Constraining the Dense Matter Equation of State with Joint Analysis of NICER and LIGO/Virgo Measurements

The Neutron Star Interior Composition Explorer collaboration recently published a joint estimate of the mass and the radius of PSR J0030+0451, derived via X-ray pulse-profile modeling. Raaijmakers et al. explored the implications of this measurement for the dense matter equation of state (EOS) using two parameterizations of the high-density EOS: a piecewise-polytropic model, and a model based on the speed of sound in neutron stars (NSs). In this work we obtain further constraints on the EOS following this approach, but we also include information about the tidal deformability of NSs from the gravitational wave signal of the compact binary merger GW170817. Further, we compare the constraints on the EOS to those set by the recent measurement of a 2.14 Me pulsar, included as a likelihood function approximated by a Gaussian, and find a small increase in information gain. To show the flexibility of our method, we also explore the possibility that GW170817 was a NS–black hole merger, which yields weaker constraints on the EOS.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Non-Hermitian quantum mechanics approach for extracting and emulating continuum physics based on bound-state-like calculations: Detailed description

Here, this work applies a reduced basis method to study the continuum physics of a finite quantum system—either few or many-body. Specifically, I develop reduced-order models, or emulators, for the underlying inhomogeneous Schrödinger equation and train the emulators against the equation's bound-state-like solutions at complex energies. The emulators rapidly and accurately interpolate and extrapolate the matrix elements of the Hamiltonian resolvent operator (Green's function) across a parameter space that includes both complex energy and other real-valued physical inputs in the Schrödinger equation. The spectra, discretized and compressed as the result of emulation, and the associated resolvent matrix elements (or amplitudes), have the defining characteristics of non-Hermitian quantum mechanics calculations, featuring complex eigenenergies with negative imaginary parts and branch cuts moved below the real axis in the complex energy plane. Therefore, one now has a method that extracts continuum physics from bound-state-like calculations and emulates those extractions in the input parameter space. Building on a prior Letter [Zhang, Phys. Rev. Lett. 135, 242501 (2025)], this article provides the full theoretical details, a comprehensive analysis of the method's performance, and a brief discussion of how it can be coupled with existing continuum approaches to perform emulations in their input parameter spaces.

ab initio calculations↗
DOE OSTI2024

Irradiation Results of Commercial Neutron and Gamma Sensors at the Ohio State University Research Reactor

This report serves to present the evaluation results of commercial radiation detectors (SPNDs) with the potential to accomplish the data objectives—having sufficient gamma and fast neutron sensitivity—at temperatures near 650°C. The detectors chosen for evaluation are gamma ion chambers from Exosens (previously known as Photonis) models CRGA11 and CRGE32 and tantalum-based self-powered neutron detectors (Ta-SPND) from Mirion. The evaluation was performed in a series of heat irradiations from ambient to 850°C in the 9.5-inch dry tube furnace at the Ohio State University Research Reactor (OSURR). Ion chamber counting curves were measured to evaluate sensor operability at temperature. Detector sensitivity to reactor power and temperature were measured and presented in curve fit parameters. The curve fit equations were used to identify the suggested operational temperatures based on reactor power. Overall, it was evaluated that the CRGA11 was not significantly affected by temperatures up to 650°C and is operable—with higher temperature-contributed signals—up to 700°C. The CRGE32 was more affected by the high temperatures compared to the CRGA11. As a result of increasing temperature, the leakage current was a dominating factor. While the detector can operate up to 600°C and 700°C with lowered high voltage, it is not recommended unless a suitably strong gamma flux field is present. Finally, Ta SPND did not demonstrate good performance beyond 350°C due to the presence of an unknown phenomenon at changing temperatures. The study of the phenomenon is an active research topic outside the scope of this project.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2025

MELCOR Deposition Models Applied to Noble Metals

Noble metals are chemically inert elements that form particles that can circulate around the fuel loop of a liquid-fueled molten salt reactor (LF-MSR). Due to their health-hazardous abilities, it is important to track these elements and to understand their contributions to source term during a reaction, particularly for Te as it decays into iodine. This report describes three solid deposition models that were implemented into the MELCOR reactor modeling code. The models that were implemented were Dittus-Boelter for turbulent fluid deposition, Talbot thermophoretic transport equation and Richardson-Zaki gravitational deposition. These closures were used to calculate the mass transfer coefficient and was compared to those reported for the Molten Salt Reactor Experiment (MSRE). The impacts of salt, operation and particle parameters on how they affect solid deposition are described, and found that solid deposition is strongly dependent on salt flow and viscosity, as well as particle diffusion rate and size. The conclusions form the basis for activities to be carried out in FY26 and for pursuit of future activities.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Relativistic hybrid stars in light of the NICER PSR J0740 + 6620 radius measurement

We explore the implications of the recent radius determination of PSR J0740 + 6620 by the NICER experiment, combined with the neutron skin measurement by the PREX-II experiment and the associated inference of the slope of symmetry energy, for the structure of hybrid stars with a strong first-order phase transition from nucleonic to quark matter. Furthermore, we combine a covariant density-functional nucleonic equation of state (EOS) with a constant-speed-of-sound EOS for quark matter. We show that the radius and tidal deformability ranges obtained from GW170817 can be reconciled with the implication of the PREX-II experiment if there is a phase transition to quark matter in the low-mass compact star. In the high-mass segment, the EOS needs to be stiff to comply with the large-radius inference for PSR J0740 + 6620 and J0030 + 0451 with masses M ≃ 2 M ⊙ and M≃1.4 M⊙. We show that twin stars are not excluded, but the mass and radius ranges (with M ≥ M ⊙ ) are restricted to narrow domains ΔM twin ≲ 0.05 M ⊙ and ΔR twin ~ 1.0 km. We also show that the existence of twin configurations is compatible with the light companion in the GW190814 event being a hybrid star in the case of values of the sound-speed square s = 0.6 and s = 1/3.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

QCD–Gravity Double Copy in Regge Asymptotics: From \(2\rightarrow n\) Amplitudes to Radiation in Shockwave Collisions

This paper discusses multi-particle production in QCD and in gravity at ultrarelativistic energies, their double-copy relations, and strong parallels in emergent shockwave dynamics. Dispersive techniques are applied to derive the BFKL equation for multi-gluon production in Regge asymptotics. Identical methods apply in gravity and are captured by a gravitational Lipatov equation. The building blocks in both cases are Lipatov vertices and reggeized propagators satisfying double-copy relations; in gravity, Weinberg’s soft theorem is recovered as a limit of the Lipatov framework. BFKL evolution in QCD generates wee parton states of maximal occupancy characterized by an emergent semi-hard saturation scale. Renormalization group equations in the Color Glass Condensate (CGC) EFT describe wee parton correlations and their rapidity evolution. A shockwave picture of deeply inelastic scattering and hadron–hadron collisions follows, with multi-particle production described by Cutkosky’s rules in strong time-dependent fields. Gluon radiation in the CGC EFT has a double copy in gravitational shockwave collisions, with a similar correspondence applicable between gluon and graviton shockwave propagators. Possible extensions of this semi-classical double copy are outlined for computing multi-particle production in gravitational shockwave collisions, self-force and tidal contributions, and classical and quantum noise in the focusing of geodesics.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Unraveling Electronic and Vibrational Coherences Following a Charge Transfer Process in a Photosystem II Reaction Center

A reaction center is a unique biological system that performs the initial charge separation within a Photosystem II (PSII) multiunit enzyme, which eventually drives the catalytic water-splitting in plants and algae. The possible role of quantum coherences coinciding with the energy and charge transfer processes in PSII reaction center is one of the active areas of research. Here, we study these quantum coherences by using a numerically exact method on an excitonic dimer model, including linear vibronic coupling and employing optimal parameters from experimental two-dimensional coherent spectroscopic measurements. This enables us to precisely capture the excitonic interaction between pigments and the dissipation of the energy from electronic and charge-transfer (CT) states to the protein environment. We employ the time nonlocal (TNL) quantum master equation to calculate the population dynamics, which yields numerically reliable results. The calculated results show that, due to the strong dissipation, the lifetime of electronic coherence is too short to have direct participation in the charge transfer processes. However, there are long-lived vibrational coherences present in the system at frequencies close to the excitionic energy gap. These are strongly coupled with the electronic coherences, which makes the detection of the electronic coherences with conventional techniques very challenging. Additionally, we unravel the strong excitonic interaction of radical pair (PD1 and PD2) in the reaction center, which results in a long-lived electronic coherence of >100 fs, even at room temperature. Our work provide important physical insight to the charge separation process in PSII reaction center, which may be helpful for better understanding of photophysical processes in other natural and artificial light-harvesting systems.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Greedy emulators for nuclear two-body scattering

Applications of reduced basis method emulators are increasing in low-energy nuclear physics because they enable fast and accurate sampling of high-fidelity calculations, enabling robust uncertainty quantification. Here, in this paper, we develop, implement, and test two model-driven emulators based on the (Petrov-)Galerkin projection using the prototypical test case of two-body scattering with the Minnesota potential and a more realistic local chiral potential. The high-fidelity scattering equations are solved with the matrix Numerov method, a reformulation of the popular Numerov recurrence relation for solving special second-order differential equations as a linear system of coupled equations. A novel error estimator based on reduced-space residuals is applied to an active learning approach (a greedy algorithm) to choosing training samples (“snapshots”) for the emulator and contrasted with a proper orthogonal decomposition (POD) approach. Both approaches allow for computationally efficient offline-online decompositions, but the greedy approach requires many fewer snapshot calculations. These developments set the groundwork for emulating scattering observables based on chiral nucleon-nucleon and three-nucleon interactions and optical models, where computational speed-ups are necessary for Bayesian uncertainty quantification. Our emulators and error estimators are widely applicable to linear systems.

Bayesian methods↗
DOE OSTI2025

Effects of Quark Core Sizes of Baryons in Neutron Star Matter

We describe the quark substructure of hadrons and the equation of state of high-density neutron star matter by using the Nambu–Jona-Lasinio (NJL) model, which is an effective quark theory based on QCD. The interaction between quarks fully respects the chiral and flavor symmetries. Guided by the success of various low-energy theorems, we assume that the explicit breaking of these symmetries occurs only via the current quark masses, and all other symmetry breakings are of dynamical nature. In order to take into account the effects of the finite quark core sizes of the baryons on the equation of state, we make use of an excluded volume framework that respects thermodynamic consistency. The effects generated by the swelling quark cores generally act repulsively and lead to an increase in the pressure with increasing baryon density. On the other hand, in neutron star matter, these effects also lead to a decrease in the density window where hyperons appear because it becomes energetically more favorable to convert the faster moving nucleons into hyperons. Our quantitative analysis shows that the net effect of the excluded volume is too small to solve the long-standing “hyperon puzzle”, which is posed by the large observed masses of neutron stars. Thus, the puzzle persists in a relativistic effective quark theory which takes into account the short-range repulsion between baryons caused by their finite and swelling quark core sizes in a phenomenological way.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Investigation into Na and Cs activity coefficients in high salt solutions to support Cs removal in Hanford tank waste

The treatment of Hanford tank waste is one of the most technically challenging environmental cleanup activities for the U.S. Department of Energy to date. To expedite the processing of liquid waste stored in underground tanks in southeastern Washington state, it is necessary to remove the significant dose contributor, 137 Cs. Toward this effort, ion exchange with crystalline silicotitanate (CST) has been employed as part of the Tank Side Cesium Removal system. The model used to predict Cs exchange onto CST was developed using activity coefficients calculated from the Bromley equation. A series of batch contact tests that varied in [Na] were conducted to look at the impact of Na concentration on Cs distribution. Experimental distribution ratios ( K d ) were compared to the distribution ratios predicted using three different activity coefficient models: (1) commercially available HSC software, (2) the Bromley equation, and (3) a simplified approach adapted from Marcos-Arroyo et al. Ultimately, the Bromley method underpredicted the effect of ionic strength on the Na activity coefficient ( γ Na+ ), HSC overestimated the impact of ionic strength on the expected performance due to the Cs activity coefficient ( γ Cs+ ), but the simplified approach predicted the experimental K d values quite well in a binary matrix. In conclusion, expansion of this approach in complex matrices is necessary for application to Hanford tank waste.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2020

Improvements to CTF Code Verification and Unit Testing (FY2020)

In 2010, the U.S. Department of Energy created its first Energy Innovation Hub, which focuses on improving Light Water Reactors (LWRs) through Modeling and Simulation. This hub, named the Consortium for the Advanced Simulation of LWRs (CASL), attempts to characterize and understand LWR behavior under normal operating conditions and use any gained insights to improve their efficiency. In collaboration with North Carolina State University (NCSU), CASL has worked extensively on the thermal-hydraulic subchannel code Coolant Boiling in Rod Arrays—Three Field (COBRA-TF). The NCSU/CASL version of COBRA-TF has been rebranded as CTF. This document focuses on code verification test problems that ensure CTF converges to the correct answer for the intended application. The suite of code verification tests are mapped to the underlying conservation equations of CTF, and significant gaps are addressed. Convergence behavior and numerical errors are quantified for each of the tests. Tests that converge at the correct rate to the corresponding analytic solution are incorporated into the CTF automated regression suite. A new verification utility is created for this purpose, which enables code verification by generalizing the process. For problems that do not behave correctly, the results are reported but the problem is not included in the regression suite. In addition to verification studies, this document also quantifies the existing tests of constitutive models. A few existing gaps are addressed by adding new unit tests.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

FY20 Improvements to CTF Code Verification and Unit Testing

In 2010, the U.S. Department of Energy created its first Energy Innovation Hub, which focuses on improving Light Water Reactors (LWRs) through Modeling and Simulation. This hub, named the Consortium for the Advanced Simulation of LWRs (CASL), attempts to characterize and understand LWR behavior under normal operating conditions and use any gained insights to improve their efficiency. In collaboration with North Carolina State University (NCSU), CASL has worked extensively on the thermal-hydraulic subchannel code Coolant Boiling in Rod Arrays–Three Field (COBRA-TF). The NCSU/CASL version of COBRA-TF has been rebranded as CTF. This document focuses on code verification test problems that ensure CTF converges to the correct answer for the intended application. The suite of code verification tests are mapped to the underlying conservation equations of CTF, and significant gaps are addressed. Convergence behavior and numerical errors are quantified for each of the tests. Tests that converge at the correct rate to the corresponding analytic solution are incorporated into the CTF automated regression suite. A new verification utility is created for this purpose, which enables code verification by generalizing the process. For problems that do not behave correctly, the results are reported but the problem is not included in the regression suite. In addition to verification studies, this document also quantifies the existing tests of constitutive models. A few existing gaps are addressed by adding new unit tests.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Systematic study of the validity of the eikonal model including uncertainties

Nuclear reactions at intermediate beam energies are often interpreted using the eikonal model. In the analysis of complex reaction probes, where few-body reaction methods are needed, the eikonal method may be used as an efficient way for describing the fragment-target reaction process. In this work, we perform a systematic study to test the validity of the eikonal approximation for nucleon-nucleus reactions. We also quantify uncertainties due to the nucleon optical potential on reaction observables. We inspect the validity of the eikonal model and its semiclassical correction by comparing it to exact solutions (obtained from solving the optical-model equation with a finite-differences method) for a wide range of reactions. We also study the effect of relativistic corrections, both kinematic and dynamic, by effectively incorporating the relativistic effects at intermediate energies. The uncertainties from a Bayesian global optical potential (KDUQ) are propagated to the observables of interest. Our study includes neutron and proton reactions on 27 Al , 40 Ca , 90 Zr , and 208 Pb , for a wide range of energies 𝐸 lab = 0–400 MeV. We calculate neutron-total cross sections (elastic and reactions) as well as proton-absorption cross sections as a function of beam energy, using the eikonal model, the eikonal model with a semiclassical correction, and the exact solution. Here, we also compute angular distributions for the methods above. Our results show that for the proton-absorption cross section, the eikonal model can be used down to around 60 MeV and the semiclassical correction extends its use to 30 MeV. However, the validity of the eikonal model for the neutron-total cross section only goes down to ≈120 MeV, a range extended to ≈ 50 MeV when using the semiclassical correction. We find the semiclassical correction to the eikonal model to be less effective in describing the angular distributions. The 1⁢𝜎 uncertainty intervals on the observables we studied is less than 5% for most of the energies considered, but increases rapidly for higher energies, namely energies outside the range of KDUQ (𝐸 lab > 200MeV).

Cluster models↗
DOE OSTI2020

Influence of 73 Rb on the ashes of accreting neutron stars

We find that the proton separation energy, S(p), of 73 Rb is –640(40) keV, deduced from the observation of β-delayed ground-state protons following the decay of 73 Sr. This lower-limit determination of the proton separation energy of 73 Rb coupled with previous upper limits from nonobservation, provides a full constraint on the mass excess with ΔM ( 73 Rb) = –46.01 ± 0.04 MeV. With this new mass excess and the excitation energy of the J π = 5/2 – isobaric-analog state (T = 3/2) in 73 Rb, an improved constraint can be put on the mass excess of 73Sr using the isobaric-multiplet mass equation (IMME), and we find ΔM( 73 Sr) = –31.98 ± 0.37 MeV. Furthermore, these new data were then used to study the composition of ashes on accreting neutron stars following Type I x-ray bursts. Counterintuitively, we find that there should be an enhanced fraction of A > 102 nuclei with more negative proton separation energies at the 72 Kr rp-process waiting point. Larger impurities of heavier nuclei in the ashes of accreting neutron stars will impact the cooling models for such astrophysical scenarios.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Shape of atomic nuclei in heavy ion collisions

In the hydrodynamic model description of heavy ion collisions, the final-state anisotropic flows v n are linearly related to the strengths of the multipole shape of the distribution of nucleons in the transverse plane, ε n : v n ∝ε n . The ε n , for n=1, 2, 3, 4, are sensitive to the shapes of the colliding ions, characterized by the quadrupole β 2 , octupole β 3 , and hexadecapole β 4 deformations. This sensitivity is investigated analytically and also in a Monte Carlo Glauber model. One observes a robust linear relation, $\langleε^{2}_{n}\rangle$ = $a^{'}_{n} + b^{'}_{n}β^{2}_{n}$, for events in a fixed centrality. The $\langleε^{2}_{1}\rangle$ has a contribution from β 3 and β 4 , and $\langleε^{2}_{3}\rangle$ from β 4 . In ultracentral collisions, there are little cross contributions between β 2 and ε 3 and between β 3 and ε 2 , but clear cross contributions are present in noncentral collisions. Additionally, $\langleε^{2}_{n}\rangle$ are insensitive to nonaxial shape parameters such as the triaxiality. This is good news because the measurements of v 2 , v 3 , and v 4 can be used to constrain simultaneously the β 2 , β 3 , and β 4 values. This is best done by comparing two colliding ions with similar mass numbers and therefore nearly identical $a^{'}_{n}$, to obtain a simple equation that relates the β n of the two species. Finally, this opens up the possibility to map the shape of the atomic nuclei at a timescale (<10 –24 s) much shorter than probed by low-energy nuclear structure physics (<10 –21 s), which ultimately may provide information complementary to that obtained in the nuclear structure experiments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Effects of Accident Thermal Loading on In-Plane Shear Behavior of Steel-Plate Composite Walls

Structural walls in safety-related nuclear facilities are required to be designed for seismic and accident thermal (due to postulated high-energy pipe break events) loading combination. Current U.S. and international codes provide limited guidance for analysis and design of walls for this loading combination. This paper describes the experimental results and observations from tests conducted on a laboratory-scale (1:4 to 1:5) test unit representing steel-plate composite (SC) walls subjected to combined in-plane (seismic) and accident thermal loading. The test unit was subjected to surface temperatures of up to 450°F in combination with cyclic in-plane loading. Results of similar experiments recently conducted in Japan are also summarized (with surface temperatures up to 570°F). Surface heating combined with the low thermal conductivity and high specific heat of concrete resulted in nonlinear thermal gradients through the thickness of the specimens. These nonlinear thermal gradients and the associated self- or internal restraint led to extensive concrete cracking. This concrete cracking reduced the initial and secant stiffness of the specimens. The initial stiffness of the heated specimens reduced to 30 to 40% of the initial stiffness of the control (unheated) specimen. The secant stiffness of the heated specimens reduced up to 50% of the secant stiffness of the control (unheated) specimen. However, the in-plane shear strength of the heated SC specimens was still approximately 10 to 30% greater than the nominal in-plane shear strength calculated, for the limit state of steel plate von Mises yielding, using AISC N690 equations and measured material properties. Evaluation of the experimental results and observations suggests that the in-plane shear strength of SC walls subjected to typical accident surface temperatures (up to 570°F) can be estimated conservatively using the current provisions of AISC N690. The stiffness for accident thermal loading combinations can be considered to reduce from cracked composite stiffness at ambient temperature to fully cracked—that is, steel-only stiffness as the surface temperature increases up to typical accident value.

Construction & Building Technology↗
DOE OSTI2020

The Multiphysics on Advanced Platforms Project

In 2015, the Lawrence Livermore National Laboratory started development of next-generation multiphysics simulation capabilities for the National Nuclear Security Administration under the Advanced Technologies Development and Mitigation (ATDM) element of the Advanced Simulation and Computing program in collaboration with the Exascale Computing Project (ECP). A key driver for this effort across the NNSA tri-lab was the emergence of advanced high performance computing (HPC) architectures based on heterogeneous compute capabilities, including GPU based systems, as part of the national drive toward exascale computing platforms at multiple Department of Energy (DOE) facilities. Developing a multiphysics code capable of meeting the various simulation needs of the NNSA as defined by the current generation of integrated codes (or ICs), initially developed as part of the Accelerated Strategic Computing Initiative (ASCI) program beginning in 1996, and able to scale to the current 100 petaflop class pre-exascale systems, as well the forthcoming exaflop class computers, is a daunting challenge. To accomplish this ambitious goal, LLNL has embraced two key themes: use of high-order numerical methods and a modular approach to code development. The LLNL next generation effort is organized under the Multi-Physics on Advanced Platforms Project (MAPP). A foundational component of MAPP is the Axom computer science (CS) toolkit which provides infrastructure for the development of modular, performance portable, multi-physics application codes. MARBL is a next-generation application code built on the Axom base to address the modeling needs of the high energy density physics (HEDP) community for simulating high-explosive, magnetic or laser driven experiments such as inertial confinement fusion (ICF), pulsed-power magneto-hydrodynamics (MHD), equation of state (EOS) and material strength studies as part of the NNSA’s stockpile stewardship program (SSP).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Non-Hermitian Quantum Mechanics Approach for Extracting and Emulating Continuum Physics Based on Bound-State-like Calculations

Here, this Letter introduces a unified emulation framework for studying continuum physics in finite quantum systems. Using a reduced basis method, we construct powerful emulators for the inhomogeneous Schrödinger equation that operate in a combined parameter space of complex energy (𝐸) and other inputs (𝜽). Within the space, the emulators simultaneously perform analytical continuation in 𝐸—extracting continuum physics from numerically simpler bound-state-like calculations—and interpolate this entire process across 𝜽. This yields a small, non-Hermitian system whose properties (e.g., resonances and scattering observables) can be rapidly predicted for any 𝜽. Crucially, the complex-𝐸 emulation provides a pathway to compute continuum observables for complex systems where advanced bound-state methods exist but direct continuum calculations are yet to be developed, while the 𝜽 emulation enables rapid parameter-space exploration and can be adapted to accelerate other existing continuum calculations. Demonstrations with two- and three-body systems highlight the method’s effectiveness and suggest its connection to (near-)optimal rational approximation. This Letter presents the key results, with further details reserved for a companion paper.

ab initio calculations↗