Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “extrapolation”

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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
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 OSTI2025

Out-of-Distribution Generalization for Learning Quantum Channels with Low-Energy Coherent States

When experimentally learning the action of a continuous-variable quantum process by probing it with inputs, there will often be some restriction on the input states used. One experimentally simple way to probe a quantum channel is to use low-energy coherent states. Learning a quantum channel in this way presents difficulties, due to the fact that two channels may act similarly on low-energy inputs but very differently for high-energy inputs. They may also act similarly on coherent-state inputs but differently on nonclassical inputs. Extrapolating the behavior of a channel for more general input states from its action on the far more limited set of low-energy coherent states is a case of out-of-distribution generalization. To be sure that such generalization gives meaningful results, one needs to relate error bounds for the training set to bounds that are valid for all inputs. We show that for any pair of channels that act sufficiently similarly on low-energy coherent-state inputs, one can bound how different the input-output relations are for any (high-energy or highly nonclassical) input. This proves that out-of-distribution generalization is always possible for learning quantum channels using low-energy coherent states, as long as enough samples are used.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Emerging tension between isospin symmetry and nuclear deformation at 𝐴 = 78

Since the nuclear interaction is, to first order, isospin independent, simple relationships can be derived for 𝐸⁢2 strengths across isobaric analog states within an isobaric triplet. At the same time, there is a well-known relationship between the 2 + excitation energy and the 𝐵(𝐸⁢2;2 + → 0 + ) transition strength in deformed even-even nuclei. Here, we highlight a tension between these two trends in the 𝐴 = 78 isobaric triplet. Recent level-energy determinations indicate that the proton-rich member of the triplet, 78 Zr, is strongly deformed while recent lifetime measurements in 78 Sr and 78 Y, as extrapolated into 78 Zr, would imply that 78 Zr is weakly deformed. Furthermore, we explore the implications of this apparent contradiction and explore potential solutions.

59 ≤ A ≤ 89↗
DOE OSTI2026

Towards a Quintic Ginzburg-Landau Description of the (2,7) Minimal Model

We discuss dimensional continuation of the massless scalar field theory with the 𝑖⁢𝜙 5 interaction term. It preserves the so-called 𝒫⁢𝒯 symmetry, which acts by 𝜙 →−𝜙 accompanied by 𝑖 →−𝑖. Below its upper critical dimension 10/3, this theory has interacting infrared fixed points. We argue that the fixed point in 𝑑 = 2 describes the nonunitary minimal conformal model 𝑀⁡(2,7). We identify the operators 𝜙 and 𝜙 2 with the Virasoro primaries 𝜙 1,2 and 𝜙 1,3 , respectively, and 𝑖⁢𝜙 3 with a quasiprimary operator, which is a Virasoro descendant of 𝜙 1,3 . Our identifications appear to be consistent with the operator product expansions and with considerations based on integrability. Using constrained Padé extrapolations, we provide estimates of the critical exponents in 𝑑 = 3. We also comment on possible lattice descriptions of 𝑀⁡(2,7) and discuss RG flows to and from this conformal field theory (CFT). Finally, we conjecture that the minimal models 𝑀⁡(2,2⁢𝑛 +1) are described by the massless scalar field theories with the 𝑖⁢𝜙 2⁢𝑛−1 interaction terms.

Conformal field theory↗
DOE OSTI2025

Advancing quantum simulations of the nuclear shell model with Gray-code–based resource-efficient protocols

Background: Some of the computational limitations in solving the nuclear many-body problem could be overcome by utilizing quantum computers. The nuclear shell-model calculations providing deeper insights into the properties of atomic nuclei are one such case with high demand for resources, as the size of the Hilbert space grows exponentially with the number of particles involved. Quantum algorithms are being developed to overcome these challenges and advance such calculations. Purpose: To develop quantum circuits for the nuclear shell-model, leveraging the capabilities of noisy intermediate-scale quantum (NISQ) devices. Here, we aim to minimize resource requirements (specifically in terms of qubits and gates) and strive to reduce the impact of noise by employing relevant mitigation techniques. Methods: We achieve noise resilience by designing an optimized Ansatz for the variational quantum eigensolver (VQE) based on Givens rotations and incorporating qubit-ADAPT-VQE in combination with variational quantum deflation (VQD) to compute ground and excited states, incorporating the zero-noise extrapolation mitigation technique. Furthermore, the qubit requirements are significantly reduced by mapping the basis states to qubits using Gray-code encoding and generalizing transformations of fermionic operators to efficiently represent many-body states. Results: By employing the resource-efficient protocols, we achieve the ground and excited state energy levels of 38 Ar and 6 Li with better accuracy. These energy levels are presented for noiseless simulations, noisy conditions, and after applying noise mitigation techniques. Results are compared for Jordan-Wigner and Gray-code encoding using VQE, qubit-ADAPT-VQE, and VQD. Conclusions: Our work highlights the potential of resource-efficient protocols to leverage the full potential of NISQ devices in scaling the nuclear shell model calculations, offering a pathway toward more complex quantum simulations in nuclear physics. This approach establishes a framework for studying other nuclear systems with improved quantum resource efficiency, marking a significant advancement in applying quantum computing to realistic nuclear physics applications.

Physics - Nuclear physics and radiation physics↗
DOE OSTI2026

Coupled-channel approach to isotensor π π π scattering from lattice QCD

The quest to understand three-body dynamics from first-principle QCD includes the study of nonresonant and resonant systems. The isospin I = 2 system is of particular interest having no three-body resonance but featuring a resonance in a subchannel, while also being a coupled-channel problem. In this study, we calculate the finite-volume spectrum from lattice QCD at two different pion masses, map the amplitude to the infinite volume through a generalized Finite-Volume Unitarity three-body quantization condition, investigate the limit of a narrow ρ , and compare with an effective Lagrangian prediction at leading order. Chiral extrapolations between different pion masses are performed.

Feng, Yuchuan [The George Washington University] (↗
DOE OSTI2025

Robustness of Vacancy-Bound Non-Abelian Anyons in the Kitaev Model in a Magnetic Field

Non-Abelian anyons in quantum spin liquids (QSLs) provide a promising route to fault-tolerant topological quantum computation. In the exactly solvable Kitaev honeycomb model, such anyons of the QSL state can be bound to nonmagnetic spin vacancies and endowed with non-Abelian statistics by an infinitesimal magnetic field. Here, we investigate how this approach for stabilizing non-Abelian anyons extends to a finite magnetic field represented by a proper Zeeman term. Through large-scale density-matrix renormalization group simulations, we compute the vacancy-anyon binding energy as a function of magnetic field for both the ferromagnetic and antiferromagnetic Kitaev models. Here, we find that anyon binding remains robust within the entire QSL phase for the ferromagnetic Kitaev model but breaks down already inside this phase for the antiferromagnetic Kitaev model. To compute a binding energy several orders of magnitude below the magnetic energy scale, we introduce both a refined definition and an extrapolation scheme based on carefully tailored perturbations.

Xiao, Bo [Oak Ridge National Laboratory (ORNL), Oa↗
DOE OSTI2025

Mechanically and electrically tunable Rashba-Edelstein effect in ferroelectric semiconductors, CsGe⁢𝑋 3 (𝑋 = I, Br, Cl)

The ability of materials to convert charge current into spin current is fundamental to many spintronics applications. One means of realizing this conversion is via Rashba-Edelstein effect (REE). Using density functional theory simulations, we predict that REE can be induced in the recently discovered family of semiconducting ferroelectrics, CsGe⁢𝑋 3 (𝑋 = I, Br, Cl). The effect is quantified through Rashba-Edelstein coefficients, 𝜒 𝑥⁢𝑦 =−𝜒 𝑦⁢𝑥 , which are nonzero in valence and conduction bands. The largest values, obtained for CsGeI 3 , are 3.45 × 10 10 $\frac{ℏ}{Λ cm}$ and 0.97 × 10 10 $\frac{ℏ}{Λ cm}$ in the conduction and valence bands, respectively. The values are comparable to, and sometimes exceed, those of other inorganic materials, although the maximal values occur away from the band edges. The coefficients' sign couples to the direction of spontaneous polarization, offering opportunities for nonvolatile spin current manipulation via external electric field. Furthermore, these coefficients are highly tunable through strain engineering owing to strain-induced variations in spin textures and energy isosurfaces. Specifically, 𝜒 in the valence band of CsGeI 3 is enhanced to 3.61 × 10 10 $\frac{ℏ}{Λ cm}$ under 5% biaxial strain. Given this potential, we synthesized CsGeI 3 to validate our simulation structure and found excellent agreement between experiment and simulation, thereby allowing extrapolation of our findings to practically significant temperatures. Furthermore, our study identifies promising materials for nonvolatile, multifunctional spintronic applications.

Electronic structure↗
DOE OSTI2025

New state of matter between the hadronic phase and the quark-gluon plasma?

Lattice-quantum chromodynamics (QCD) simulations and theoretical arguments hint at the existence of an intermediate phase of strongly interacting matter between a confined hadron gas and a deconfined quark-gluon plasma (QGP). We qualitatively and semiquantitatively explore and differentiate the phase structures in the temperature window from the QCD pseudocritical temperature 𝑇 c ≃ 160 MeV to the pure gluonic deconfinement temperature 𝑇 d ≃ 285 MeV. We propose a three-regime picture using a hadron resonance gas description augmented with the exponential spectrum of strings, corresponding to highly excited mesons and glueballs, based on the analysis of a large number 𝑁 c of colors. We estimate the entropy density from our model to confirm that the lattice-QCD data are bracketed with three regimes, i.e., a hadron gas, a QGP, and a new phase for 𝑇 c ≲ 𝑇 ≲ 𝑇 d . In this new phase, which we name a spaghetti of quarks with glueballs (SQGBs), thermal degrees of freedom of quarks are liberated, yet gluons remain confined in glueballs. Since the Hagedorn temperature 𝑇 H ∼ 285 MeV is universal in the meson and the glueball sectors, in the infinite-𝑁 c limit, the phase diagram in the plane of the baryon chemical potential and the temperature is reduced to one with the confined and deconfined phases and quarkyonic matter at high density. At large but finite 𝑁 c , an SQGB window may open between these phases. We point out that the SQGB has interesting similarities with quarkyonic matter and that this matter in the large-𝑁 c limit is confined as measured by the interaction between heavy quarks, but behaves in other respects like a quasifree gas of quarks. As a result of the extrapolation to 𝑁 c = 3, we present a revised phase diagram with the SQGB phase bounded by thermal crossovers. Finally, we give a quantitative analysis of chiral-symmetry restoration in the SQGB phase.

Color confinement↗
DOE OSTI2026

Proton isovector helicity PDF at NNLO and the twist-3 moment $\tilde{d}$ 2 from lattice QCD at physical quark masses

We present a lattice quantum chromodynamics calculation of the 𝑥-dependent isovector quark helicity parton distribution function (PDF) of the proton in the large momentum effective theory (LaMET) framework. Through operator product expansion (OPE) we also extract the $\tilde{d}$ 2 moment of the twist-3 PDF 𝑔 𝑇 ⁡(𝑥) for the first time in the $\overline{MS}$ scheme, which is proportional to the average color Lorentz force experienced by the quark in the proton. This calculation is performed on a lattice of spacing 𝑎 =0.076 fm at physical quark masses. The quasi-PDF matrix elements are measured in proton states boosted to momenta 𝑃 𝑧 ={0,0.25,1.02,1.53} GeV. We first extract the lowest few helicity PDF moments from the renormalization-group (RG) invariant ratios of the matrix elements with OPE. Combined with the matrix elements relevant for 𝑔 𝑇 ⁡(𝑥), we obtain $\tilde{d}$$^{u-d}_2$⁡(2 GeV) =0.0024⁢(46) at next-to-leading order in $\overline{MS}$. Then, the helicity quasi-PDF matrix elements are renormalized in the hybrid scheme with linear renormalon resummation and Fourier transformed to the 𝑥-space after an asymptotic extrapolation. The quasi-PDF is perturbatively matched to the $\overline{MS}$ PDF with RG and threshold resummations at next-to-leading power and next-to-next-to-leading logarithmic accuracies. After resummations, we determine the PDF in the region 𝑥 ∈[0.25,0.75]. The end-point regions are then parametrized, combined with the LaMET prediction at moderate 𝑥, and fitted to the short-distance matrix elements in coordinate space.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Ab initio calculations of monopole sum rules: From finite nuclei to infinite nuclear matter

We compute moments of the isoscalar monopole response of 𝑁 = 𝑍 closed-shell nuclei based on chiral nucleon-nucleon plus three-nucleon interactions. We employ the random-phase approximation (RPA) and two ab initio many-body approaches, the in-medium similarity renormalization group (IMSRG) and coupled-cluster theory (CC). In the IMSRG framework, the moments are obtained as ground-state expectation values, whereas in the CC approach, they are evaluated through excited-state calculations. We find good agreement between the IMSRG and CC results across all nuclei studied. RPA provides a reasonable approximation to the correlated methods if the interaction is soft. From the calculated moments, we extract average energies of the monopole response, compute finite-nucleus incompressibilities, and estimate the incompressibility of symmetric nuclear matter by a fit to a leptodermous expansion. Our extrapolated values are lower than those obtained in nuclear-matter calculations with the same interactions, but the values are consistent with phenomenological ranges.

Bonaiti, Francesca [Michigan State Univ., East Lan↗
DOE OSTI2021

Resolving Discrepancies in Spin-Torque Ferromagnetic Resonance Measurements: Lineshape versus Linewidth Analyses

When spin-orbit torques are measured using spin-torque ferromagnetic resonance, two alternative ways of analyzing the results to extract the torque efficiencies—lineshape analysis and analysis of the change in linewidth versus direct current—often give inconsistent results. We identify a source for these inconsistencies. We show that fits of spin-torque ferromagnetic resonance data to the standard analysis framework leave significant residuals that we identify as due to (i) current-induced excitations of a small volume of magnetic material with magnetic damping much larger than that of the bulk of the magnetic layer, which we speculate is associated with the heavy-metal–magnet interface and (ii) oscillations of the sample magnetization at the modulation frequency due to heating. Here, the dependence of the residual signals on direct current can interfere with an accurate extraction of spin-torque efficiencies by the linewidth method. We show that the discrepancies between the two types of analysis can be largely eliminated by extrapolation of the window of magnetic fields used in the linewidth fits to small values so as to minimize the influence of the residual signals.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Neutrino-detector design for safeguarding small modular reactors

Nuclear reactors have long been a favored source for antineutrino measurements for estimates of power and burnup. With appropriate detector parameters and background rejection, an estimate of the reactor power can be derived from the measured antineutrino event rate. Antineutrino detectors are potentially attractive as a safeguards technology that can monitor reactor operations and thermal power from a distance. Advanced reactors have diverse features that may present challenges for current safeguards methods. By comparison, neutrino detectors offer complementary features, including a remote, continuous, unattended, and near-real-time monitoring capability, that may make them useful for safeguarding certain classes of advanced reactors. This study investigates the minimum depth and size of an antineutrino detector for a small modular reactor to meet safeguards needs for advanced reactors. Furthermore, extrapolating performance from several prior reactor antineutrino experiments, this study uses an analytical approach to develop a possible design for a remote antineutrino-based monitoring device.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Pressure tuning of structural and magnetic transitions in EuAg 4 As 2

We report temperature-dependent measurements of ambient-pressure specific heat, magnetic susceptibility, anisotropic resistivity, and thermal expansion as well as in-plane resistivity under pressure up to 20.8 kbar on single crystals of EuAg 4 As 2 . Based on thermal expansion and in-plane electrical transport measurements at ambient pressure this compound has two, first-order, structural transitions in the 80–120 K temperature range. Ambient-pressure specific heat, magnetization, and thermal expansion measurements show a cascade of up to seven transitions between 8 and 16 K associated with the ordering of the $Eu^{2+}$ moments. In-plane electrical transport is able to detect the more prominent of these transitions, at 15.5, 9.9, and 8.7 K, as well as a weak feature at 11.8 K at ambient pressure. Pressure-dependent electrical transport data show that the magnetic transitions shift to higher temperatures under pressure, as does the upper structural transition, whereas the lower structural transition is suppressed and ultimately vanishes. A jump in resistivity, associated with the upper structural transition, decreases under pressure with an extrapolated disappearance (or a change of sign) by 30–35 kbar. In the 10–15 kbar range a kink in the pressure dependency of the upper structural transition temperature as well as the high- and low-temperature in-plane resistivities suggests that a change in the electronic structure may occur in this pressure range. The results are compared with the literature data for EuAg 4 As 2 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

First-principles coupled cluster theory of the electronic spectrum of transition metal dichalcogenides

The electronic properties of two-dimensional transition metal dichalcogenides (2D TMDs) have attracted much attention during the last decade. Here we show how a diagrammatic ab initio coupled cluster singles and doubles (CCSD) treatment paired with a careful thermodynamic limit extrapolation in two dimensions can be used to obtain converged band gaps for monolayer materials in the MoS 2 family. We find CCSD gaps to lie in the upper range of the spread of $\textit{GW}$ approximation based on density functional theory (DFT) simulations, and also find slightly higher effective hole masses compared to previous reports. We also investigate the ability of CCSD to describe trion states, finding a reasonable qualitative structure, but poor excitation energies due to the lack of screening of three-particle excitations in the effective Hamiltonian. Our study provides an independent high-level benchmark of the role of many-body effects in 2D TMDs and showcases the potential strengths and weaknesses of diagrammatic coupled cluster approaches for realistic materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Tensor network study of the m = 1/2 magnetization plateau in the Shastry-Sutherland model at finite temperature

The two-dimensional infinite projected entangled pair state tensor network is evolved in imaginary time with the full update (FU) algorithm to simulate the Shastry-Sutherland model in a magnetic field at finite temperature directly in the thermodynamic limit. We focus on the phase transition into the m = 1/2 magnetization plateau, which was observed in experiments on SrCu 2 (BO 3 ) 2 . For the largest simulated bond dimension, the early evolution in the high-temperature regime is simulated with the simple update (SU) scheme and then, as the correlation length increases, continued with the FU scheme towards the critical regime. We apply a small symmetry-breaking bias field and then extrapolate towards zero bias using a simple scaling theory in the bias field. The combined SU + FU scheme provides an accurate estimate of the critical temperature, even though the results could not be fully converged in the bond dimension in the vicinity of the transition. Here, the critical temperature estimate is improved with a generalized scaling theory that combines two divergent length scales: One due to the bias, and the other due to the finite bond dimension. The obtained results are consistent with the transition being in the universality class of the two-dimensional classical Ising model. The estimated critical temperature is 3.5(2) K, which is well above the temperature 2.1 K used in the experiments.

2-dimensional systems↗
DOE OSTI2021

Cohesion and excitations of diamond-structure silicon by quantum Monte Carlo: Benchmarks and control of systematic biases

We have carried out quantum Monte Carlo (QMC) calculations of silicon crystal focusing on the accuracy and systematic biases that affect the electronic structure characteristics. The results show that 64 and 216 atom supercells provide an excellent consistency for extrapolated energies per atom in the thermodynamic limit for ground, excited, and ionized states. We have calculated the ground state cohesion energy with both systematic and statistical errors below ≈ 0.05 eV. The ground state exhibits a fixed-node error of only 1.3(2)% of the correlation energy, suggesting an unusually high accuracy of the corresponding single-reference trial wave function. We obtain a very good agreement between optical and quasiparticle gaps that affirms the marginal impact of excitonic effects. Our most accurate results for band gaps differ from the experiments by about 0.2 eV. This difference is assigned to a combination of residual finite-size and fixed-node errors. In conclusion, we have estimated the crystal Fermi level referenced to vacuum that enabled us to calculate the edges of valence and conduction bands in agreement with experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Plasmonically assisted channels of photoemission from metals

We analyze recently measured nonlinear photoemission spectra from Ag surfaces that reveal resonances whose energies do not scale with the applied photon energy but stay pinned to multiples of bulk plasmon energy $\hbar$ω p above the Fermi level. To elucidate these unexpected and peculiar features we investigate the spectra of plasmons generated in a solid by the optically pumped electronic polarization and their effect on photoemission. By combining quadratic response formalism for calculations of photoemission yield, a nonperturbative approach to inelastic electron scattering, and first-principles calculations for the electronic structure, we demonstrate the dependence of probability amplitude for single- and multiplasmon excitations on the basic parameters characterizing the photon pulse and the system. The resulting multiexcitation spectrum evolves towards a truncated plasmonic coherent state. Analogous concept is extrapolated to interpret plasmon generation by multiphoton excited electronic polarization. Based on this we elaborate a scenario that the thus created real plasmons act as supplementary frequency-locked pump field for non-Einsteinian plasmonically assisted channels of photoemission from metals. The established paradigm enables assignment and assessment of the observed linear $\hbar$ω p and nonlinear 2$\hbar$ω p electron yields from Ag. Such effects may be exploited for selective filtering of optical energy conversion in electronic systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗