Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “INTERPOLATION”

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 1,621 records · Page 90
All sourcesNASA NTRSDOE OSTI
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↗
DOE OSTI2025

Beam optics ramping in underconstrained lattice design: Application to Electron-Ion Collider hadron storage ring cooling section

This paper presents the lattice design and optics ramping strategy for the cooling section of the hadron storage ring at the electron-ion collider. The main challenge is that available tuning knobs exceed beam-optics constraints. Independently optimized injection and top-energy optics often yield disconnected solutions, making interpolation impossible. To address this, we propose two new methods. The first is an intermediate-penalty scheme that ensures ramping path continuity by penalizing constraint violations at intermediate states. The second is a continuation-based approach that adapts high-energy optics to low energy, guided by an adaptive weighting scheme to balance injection and ramping constraints. The solutions meet all beam dynamics and hardware limits. The two methods offer a general strategy for ramping in systems where the solution space is under-constrained and the starting and target configurations are far apart.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Real-time scattering in Ising field theory using matrix product states

We study scattering in Ising field theory (IFT) using matrix product states and the time-dependent variational principle. IFT is a one-parameter family of strongly coupled nonintegrable quantum field theories in 1+1 dimensions, interpolating between massive free fermion theory and Zamolodchikov's integrable massive 𝐸 8 theory. Particles in IFT may scatter either elastically or inelastically. In the postcollision wave function, particle tracks from all final-state channels occur in superposition; processes of interest can be isolated by projecting the wave function onto definite particle sectors, or by evaluating energy density correlation functions. Using numerical simulations we determine the time delay of elastic scattering and the probability of inelastic particle production as a function of collision energy. We also study the mass and width of the lightest resonance near the 𝐸 8 point in detail. Close to both the free fermion and 𝐸 8 theories, our results for both elastic and inelastic scattering are in good agreement with expectations from form-factor perturbation theory. Using numerical computations to go beyond the regime accessible by perturbation theory, we find that the high-energy behavior of the two-to-two particle scattering probability in IFT is consistent with a conjecture of Zamolodchikov. Our results demonstrate the efficacy of tensor-network methods for simulating the real-time dynamics of strongly coupled quantum field theories in 1+1 dimensions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

First Look at Quartic-in-Spin Binary Dynamics at Third Post-Minkowskian Order

We compute the conservative and radiation-reaction contributions to classical observables in the gravitational scattering between a spinning and a spinless black hole to the fourth order in spin and third order in the gravitational constant. The conservative results are obtained from two-loop amplitudes for the scattering process of a massive scalar with a massive spin-𝑠 field (𝑠=0, 1, 2) minimally coupled to gravity, employing the recently introduced spin interpolation method to resolve all spin-Casimir terms. The two-loop amplitude exhibits a spin-shift symmetry in both probe limits, which we conjecture to be a sign of yet unknown integrability of Kerr orbits through the quartic order in spin and to all orders in the gravitational constant. We obtain the radial action from the finite part of the amplitude and use it to compute classical observables, including the impulse and spin kick. This is done using the recently introduced covariant Dirac brackets, which allow for the computation of classical scattering observables for general (nonaligned) spin configurations. Finally, employing the radiation-reaction amplitude proposed by Alessio and Di Vecchia, together with the Dirac brackets, we obtain radiation-reaction contributions to observables at all orders in spin and beyond the aligned-spin limit. We find agreement with known results up to the quadratic order in spin for both conservative and radiation-reaction contributions. Our results advance the state of the art in the understanding of spinning binary dynamics in general relativity and demonstrate the power and simplicity of the Dirac bracket formalism for relating scattering amplitudes to classical observables.

general relativity↗
DOE OSTI2026

Di-nucleons do not form bound states at heavy pion mass

We perform a high-statistics lattice QCD calculation of the low-energy two-nucleon scattering amplitudes. To address discrepancies in the literature, the calculation is performed at a heavy pion mass in the limit that the light quark masses are equal to the physical strange quark mass, 𝑚 𝜋 = 𝑚 𝐾 ≃ 714 MeV. Using a state-of-the-art momentum space method, we rule out the presence of a bound di-nucleon in both the isospin 0 (deuteron) and 1 (di-neutron) channels, in contrast with many previous results that made use of compact hexaquark creation operators. To diagnose the discrepancy, we add such hexaquark interpolating operators to our basis and find that they do not affect the determination of the two-nucleon finite-volume spectrum, and thus they do not couple to deeply bound di-nucleons that are missed by the momentum-space operators. Furthermore, we perform a high-statistics calculation of the HAL QCD potential on the same gauge ensembles and find qualitative agreement with our main results. We conclude that di-nucleons do not form bound states at heavy pion masses and that previous identification of deeply bound di-nucleons must have arisen from a misidentification of the spectrum from off-diagonal elements of a correlation function.

Physics - Physics of elementary particles and fiel↗
DOE OSTI2025

Block Lanczos algorithm for lattice QCD spectroscopy and matrix elements

Recent work introduced a new framework for analyzing correlation functions with improved convergence and signal-to-noise properties, as well as rigorous quantification of excited-state effects, based on the Lanczos algorithm and spurious eigenvalue filtering with the Cullum-Willoughby test. Here, we extend this framework to the analysis of correlation-function matrices built from multiple interpolating operators in lattice quantum chromodynamics (QCD) by constructing an oblique generalization of the block Lanczos algorithm, as well as a new physically motivated reformulation of the Cullum-Willoughby test that generalizes to block Lanczos straightforwardly. The resulting block Lanczos method directly extends generalized eigenvalue problem (GEVP) methods, which can be viewed as applying a single iteration of block Lanczos. Block Lanczos provides qualitative and quantitative advantages over GEVP methods analogous to the benefits of Lanczos over the standard effective mass, including faster convergence to ground- and excited-state energies, explicitly computable two-sided error bounds, straightforward extraction of matrix elements of external currents, and asymptotically constant signal-to-noise. No fits or statistical inference are required. Proof-of-principle calculations are performed for noiseless mock-data examples as well as two-by-two proton correlation-function matrices in lattice QCD.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Shear viscosity from perturbative quantum chromodynamics to the hadron resonance gas at finite baryon, strangeness, and electric charge densities

Through model-to-data comparisons from heavy-ion collisions, it has been shown that the quark gluon plasma has an extremely small shear viscosity at vanishing densities. At large baryon densities, significantly less is known about the nature of the shear viscosity from quantum chromodynamics (QCD). Within heavy-ion collisions, there are three conserved charges: baryon number (B), strangeness (S), and electric charge (Q). Here we calculate the shear viscosity in two limits using perturbative QCD (pQCD) and an excluded-volume hadron resonance gas at finite BSQ densities. We then develop a framework that interpolates between these two limits such that shear viscosity is possible to calculate across a wide range of finite BSQ densities. We find that the pQCD and hadron resonance gas calculations have different BSQ density dependencies such that a rather nontrivial shear viscosity appears at finite densities.

Phenomenology↗
DOE OSTI2025

Krylov complexity in mixed phase space

We investigate the Krylov complexity of thermofield double states in systems with mixed phase space, uncovering a direct correlation with the Brody distribution, which interpolates between Poisson and Wigner statistics. Our analysis spans two-dimensional random matrix models featuring (I) GOE-Poisson and (II) GUE-Poisson transitions and extends to higher-dimensional cases, including a stringy matrix model (GOE-Poisson) and the mass-deformed SYK model (GUE-Poisson). Krylov complexity consistently emerges as a reliable marker of quantum chaos, displaying a characteristic peak in the chaotic regime that gradually diminishes as the Brody parameter approaches zero, signaling a shift toward integrability. These results establish Krylov complexity as a powerful diagnostic of quantum chaos and highlight its interplay with eigenvalue statistics in mixed phase systems.

chaos & nonlinear dynamics↗
DOE OSTI2026

Back-to-back dijet production in DIS at arbitrary Bjorken x: TMD gluon distributions to twist-3 accuracy

We derive the gluon transverse-momentum-dependent (TMD) operator structure of back-to-back\\\\r\\\\nquark–antiquark dijet production in deep inelastic scattering at arbitrary Bjorken-x to twist-3 ac\\\\r\\\\ncuracy. Working at leading order in the strong coupling and in the kinematic regime where the\\\\r\\\\ntransverse momentum imbalance of the jets is much smaller than their individual transverse mo\\\\r\\\\nmenta, we perform a systematic gradient expansion of the quark propagator in a background gluon\\\\r\\\\nfield. This expansion organizes multiple interactions with the target in terms of longitudinal Wilson\\\\r\\\\nlines and gauge-invariant field-strength insertions, yielding a TMD description valid beyond the\\\\r\\\\nstrict high-energy eikonal (x → 0) approximation. We obtain explicit cross sections for longitudi\\\\r\\\\nnally and transversely polarized virtual photons, identifying all contributing gluon TMD operators\\\\r\\\\nup to twist-3, including structures involving F+−, Fij, and three-gluon correlators. The full lon\\\\r\\\\ngitudinal phase eixP+z− associated with Bjorken-x is retained throughout. In the small-x limit,\\\\r\\\\nour results reproduce the known sub-eikonal expressions obtained in the Color Glass Condensate\\\\r\\\\nframework, establishing a direct connection between the general-x TMD expansion and high-energy\\\\r\\\\nfactorization. We further reduce the operator basis using equations of motion, minimizing the num\\\\r\\\\nber of independent nonperturbative matrix elements entering the cross section. This work provides\\\\r\\\\na systematic foundation for extending TMD analyses of dijet production beyond leading twist, es\\\\r\\\\ntablishing a unified operator framework valid at arbitrary Bjorken-x that smoothly interpolates\\\\r\\\\nbetween moderate- and small-x descriptions of gluon TMDs.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Thermodynamics and magnetism in the two-dimensional to three-dimensional crossover of the Hubbard model

The realization of antiferromagnetic (AF) correlations in ultracold fermionic atoms on an optical lattice is a significant achievement. Experiments have been carried out in one, two, and three dimensions, and have also studied anisotropic configurations with stronger tunneling in some lattice directions. Such anisotropy is relevant to the physics of cuprate superconductors and other strongly correlated materials. Moreover, this anisotropy might be harnessed to enhance AF order. Here we numerically investigate, using the determinant quantum Monte Carlo method, a simple realization of anisotropy in the three-dimensional (3D) Hubbard model in which the tunneling between planes, $t_⊥$, is unequal to the intraplane tunneling $\textit{t}$. This model interpolates between the three-dimensional isotropic ($t_⊥ = t$) and two-dimensional (2D;$t_⊥$ = 0) systems. In this work, we show that at fixed interaction strength to tunneling ratio ($\textit{U/t}$), anisotropy can enhance the magnetic structure factor relative to both 2D and 3D results. However, this enhancement occurs at interaction strengths below those for which the Néel temperature $T_{\text{Néel}}$ is largest, in such a way that the structure factor cannot be made to exceed its value in isotropic 3D systems at the optimal $\textit{U/t}$. We characterize the 2D-3D crossover in terms of the magnetic structure factor, real space spin correlations, number of doubly occupied sites, and thermodynamic observables. An interesting implication of our results stems from the entropy's dependence on anisotropy. As the system evolves from 3D to 2D, the entropy at a fixed temperature increases. Correspondingly, at fixed entropy, the temperature will decrease going from 3D to 2D. This suggests a cooling protocol in which the dimensionality is adiabatically changed from 3D to 2D.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Simulating non-native cubic interactions on noisy quantum machines

As a milestone for general-purpose computing machines, we demonstrate that quantum processors can be programed to efficiently simulate dynamics that are not native to the hardware. Moreover, on noisy devices without error correction, we show that simulation results are significantly improved when the quantum program is compiled using modular gates instead of a restricted set of standard gates. We demonstrate the general methodology by solving a cubic interaction problem, which appears in nonlinear optics, gauge theories, as well as plasma and fluid dynamics. To encode the non-native Hamiltonian evolution, we decompose the Hilbert space into a direct sum of invariant subspaces in which the nonlinear problem is mapped to a finite-dimensional Hamiltonian simulation problem. Furthermore, in a three-states example, the resultant unitary evolution is realized by a product of approximately 20 standard gates, using which approximately ten simulation steps can be carried out on state-of-the-art quantum hardware before results are corrupted by decoherence. In comparison, the simulation depth is improved by more than an order of magnitude when the unitary evolution is realized as a single cubic gate, which is compiled directly using optimal control. Alternatively, parametric gates may also be compiled by interpolating control pulses. Modular gates thus obtained provide high-fidelity building blocks for quantum Hamiltonian simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Computationally efficient zero-noise extrapolation for quantum-gate-error mitigation

Zero noise extrapolation (ZNE) is a widely used technique for gate error mitigation on near term quantum computers because it can be implemented in software and does not require knowledge of the quantum computer noise parameters. Traditional ZNE requires a significant resource overhead in terms of quantum operations. A recent proposal using a targeted (or random) instead of fixed identity insertion method (riim versus fiim) requires significantly fewer quantum gates for the same formal precision. We start by showing that riim can allow for ZNE to be deployed on deeper circuits than fiim but requires many more measurements to maintain the same statistical uncertainty. We develop two extensions to fiim and riim. The List Identity Insertion Method (liim) allows to mitigate the error from certain cnot gates, typically those with the largest error. Set Identity Insertion Method (siim) naturally interpolates between the measurement-efficient fiim and the gate-efficient riim allowing to trade off fewer cnot gates for more measurements. Finally, we investigate a way to boost the number of measurements, namely to run ZNE in parallel, utilizing as many quantum devices as are available. We explore the performance of riim in a parallel setting where there is a non-trivial spread in noise across sets of qubits within or across quantum computers.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Quantum error mitigation by layerwise Richardson extrapolation

A widely used method for mitigating errors in noisy quantum computers is Richardson extrapolation, a technique in which the overall effect of noise on the estimation of quantum expectation values is captured by a single parameter that, after being scaled to larger values, is eventually extrapolated to the zero-noise limit. We generalize this approach by introducing layerwise Richardson extrapolation (LRE), an error mitigation protocol in which the noise of different individual layers (or larger chunks of the circuit) is amplified and the associated expectation values are linearly combined to estimate the zero-noise limit. The coefficients of the linear combination are analytically obtained from the theory of multivariate Lagrange interpolation. LRE leverages the flexible configurational space of layerwise unitary folding, allowing for a more nuanced mitigation of errors by treating the noise level of each layer of the quantum circuit as an independent variable. Furthermore, we provide numerical simulations demonstrating scenarios where LRE achieves superior performance compared to traditional (single-variable) Richardson extrapolation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Deterministic generation of remote entanglement with active quantum feedback

We develop and study protocols for deterministic remote entanglement generation using quantum feedback, without relying on an entangling Hamiltonian. In order to formulate the most effective experimentally feasible protocol, we introduce the notion of average-sense locally optimal feedback protocols, which do not require real-time quantum state estimation, a difficult component of real-time quantum feedback control. We use this notion of optimality to construct two protocols that can deterministically create maximal entanglement: a semiclassical feedback protocol for low-efficiency measurements and a quantum feedback protocol for high-efficiency measurements. The latter reduces to direct feedback in the continuous-time limit, whose dynamics can be modeled by a Wiseman-Milburn feedback master equation, which yields an analytic solution in the limit of unit measurement efficiency. Our formalism can smoothly interpolate between continuous-time and discrete-time descriptions of feedback dynamics and we exploit this feature to derive a superior hybrid protocol for arbitrary nonunit measurement efficiency that switches between quantum and semiclassical protocols. Lastly, we show using simulations incorporating experimental imperfections that deterministic entanglement of remote superconducting qubits may be achieved with current technology using the continuous-time feedback protocol alone.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Effective Landau-type model of a Hf x Zr 1 - x O 2 -graphene nanostructure

Here, to describe the charge-polarization coupling in the nanostructure formed by a thin Hf x Zr 1-x O 2 film with a single-layer graphene as a top electrode, we develop the “effective” Landau-Ginzburg-Devonshire model. This approach is based on the parametrization of the Landau expansion coefficients for the polar (FE) and antipolar (AFE) orderings in thin Hf 1-x Zr x O 2 films from a limited number of polarization-field curves and hysteresis loops. The Landau expansion coefficients are nonlinearly dependent on the film thickness h and Zr/[Hf+Zr] ratio x, in contrast to h-independent and linearly x-dependent expansion coefficients of a classical Landau energy. We explain the dependence of the Landau expansion coefficients by the strong nonmonotonic dependence of the Hf 1-x Zr x O 2 film polar properties on the film thickness, grain size and surface energy. The proposed Landau free energy with five “effective” expansion coefficients, which are interpolation functions of x and h, describes the continuous transformation of polarization dependences on applied electric field and hysteresis loop shapes induced by the changes of x and h in the range 0 < x < 1 and 5 nm < h < 35 nm. Using the effective free energy, we demonstrated that the polarization of Hf 1-x Zr x O 2 film influences strongly on the graphene conductivity, and the full correlation between the distribution of polarization and charge carriers in graphene is revealed. In accordance with our modeling, the polarization of the (5 – 25) nm thick Hf 1-x Zr x O 2 films, which are in the ferroelectric-like or antiferroelectric-like states for the chemical compositions 0.35 ≤ x ≤ 0.95, determine the concentration of carriers in graphene and can control its field dependence. The result can be promising for creation of next generation Si-compatible nonvolatile memories and graphene-ferroelectric FETs, because the working voltages applied to the Hf 1-x Zr x O 2 film (which acts as a gate) can be relatively low (less than 2 V). These low voltages are sufficient to induce the pronounced hysteresis of ferroelectric polarization in the Hf 1-x Zr x O 2 gate, which, due to the strong electric coupling, induces the hysteresis of the graphene charge.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Continuous momentum dependence in the dynamical cluster approximation

The dynamical cluster approximation (DCA) is a quantum cluster extension to the single-site dynamical mean-field theory that incorporates spatially nonlocal dynamic correlations systematically and nonperturbatively. The DCA + algorithm addresses the cluster shape dependence of the DCA and improves the convergence with cluster size by introducing a lattice self-energy with continuous momentum dependence. However, we show that the DCA + algorithm is plagued by a fundamental problem when its self-consistency equations are formulated using the bare Green's function of the cluster. This problem is most severe in the strongly correlated regime at low doping, where the DCA + self-energy becomes overly metallic and local, and persists to cluster sizes where the standard DCA has long converged. In view of the failure of the DCA + algorithm, we propose to complement DCA simulations with a post-interpolation procedure for single-particle and two-particle correlation functions to preserve continuous momentum dependence and the associated benefits in the DCA. We demonstrate the effectiveness of this practical approach with results for the half-filled and hole-doped two-dimensional Hubbard model.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Nonmonotonic crossover in electronic phase separated manganite superlattices driven by the superlattice period

Studying manganite superlattices [(LCMO) 2n / (PCMO) n ] t made of La 0.625 Ca 0.375 MnO 3 (LCMO) and Pr 0.625 Ca 0.375 MnO 3 (PCMO), we found an unexpected behavior varying the period n. At small n, the ensemble is a three-dimensional ferromagnetic metal due to interfacial charge transfer. At large n , the LCMO layers dominate transport. However, rather than a smooth interpolation between these limits a sharp transport and magnetic anomaly is found at an intermediate critical PCMO thickness n*. Magnetic force microscopy reveals that the phase-separation length scale also maximizes at n* where, unexpectedly, it becomes comparable to that of the (La 1-y Pr y ) 0.625 Ca 0.375 MnO 3 (LPCMO) alloy. Finally, we conjecture the phenomenon originates in a disorder-related length scale: Large charge-ordered clusters as in LPCMO can only nucleate when Pr-rich regions reach a critical size related to n*.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Temperature dependence of London penetration depth anisotropy in superconductors with anisotropic order parameters

Here, we study the effects of anisotropic order parameters on the temperature dependence of London penetration depth anisotropy γ λ ( T ) . After MgB 2 , this dependence is commonly attributed to distinct gaps on multiband Fermi surfaces in superconductors. We have found, however, that the anisotropy parameter may depend on temperature also in one-band materials with anisotropic order parameters Δ ( T , k F ) ; a few such examples are given. We have also found that for different order parameters, the temperature dependence of Δ ( T ) / Δ ( 0 ) can be represented with good accuracy by the interpolation suggested by Einzel, which simplifies considerably the evaluation of γ λ ( T ) . Of particular interest are mixed order parameters of two symmetries for which γ λ ( T ) may go through a maximum for a certain relative weight of two phases. Also, for this case we find that the ratio Δ max ( 0 ) / T c may exceed substantially the weak-coupling limit of 1.76. It, however, does not imply strong coupling; rather, it is due to significantly anisotropic angular variation of Δ .

36 MATERIALS SCIENCE↗