Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Hamiltonian”

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,315 records · Page 73
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

GALIC: hybrid multi-qubitwise pauli grouping for quantum computing measurement

Abstract Observable estimation is a core primitive in NISQ-era algorithms targeting quantum chemistry applications. To reduce the state preparation overhead required for accurate estimation, recent works have proposed various simultaneous measurement schemes to lower estimator variance. Two primary grouping schemes have been proposed: full commutativity (FC) and qubit-wise commutativity (QWC), with no compelling means of interpolation. In this work we propose a generalized framework for designing and analyzing context-aware hybrid FC/QWC commutativity relations. We use our framework to propose a noise-and-connectivity aware grouping strategy: Generalized backend-Aware pauLI Commutation (GALIC). We demonstrate how GALIC interpolates between FC and QWC, maintaining estimator accuracy in Hamiltonian estimation while lowering variance by an average of 20% compared to QWC. We also explore the design space of near-term quantum devices using the GALIC framework, specifically comparing device noise levels and connectivity. We find that error suppression has a more than 10 × larger impact on device-aware estimator variance than qubit connectivity with even larger correlation differences in estimator biases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Characterization and thermometry of dissipatively stabilized steady states

In this work we study the properties of dissipatively stabilized steady states of noisy quantum algorithms, exploring the extent to which they can be well approximated as thermal distributions, and proposing methods to extract the effective temperature T. We study an algorithm called the relaxational quantum eigensolver (RQE), which is one of a family of algorithms that attempt to find ground states and balance error in noisy quantum devices. In RQE, we weakly couple a second register of auxiliary ‘shadow’ qubits to the primary system in Trotterized evolution, thus engineering an approximate zero-temperature bath by periodically resetting the auxiliary qubits during the algorithm’s runtime. Balancing the infinite temperature bath of random gate error, RQE returns states with an average energy equal to a constant fraction of the ground state. We probe the steady states of this algorithm for a range of base error rates, using several methods for estimating both T and deviations from thermal behavior. In particular, we both confirm that the steady states of these systems are often well-approximated by thermal distributions, and show that the same resources used for cooling can be adopted for thermometry, yielding a fairly reliable measure of the temperature. These methods could be readily implemented in near-term quantum hardware, and for stabilizing and probing Hamiltonians where simulating approximate thermal states is hard for classical computers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Lyapunov controlled counterdiabatic quantum optimization

We introduce a quantum algorithm that integrates counterdiabatic (CD) protocols with quantum Lyapunov control (QLC) to address combinatorial optimization problems. This approach offers versatility, allowing implementation as either a digital-analog or purely digital algorithm based on selected control strategies. By examining spin-glass Hamiltonians, we illustrate how the algorithm can explore alternative paths to enhance solution outcomes compared to conventional CD techniques. This method reduces dependence on extensive higher-order CD terms and on classical optimization techniques, making it more suitable for existing quantum computing platforms. The combination of digital compression via CD protocols and the adaptable nature of QLC methods positions this approach as a promising candidate for near-term quantum devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Spin pumping from antiferromagnetic insulator spin-orbit-proximitized by adjacent heavy metal: a first-principles Floquet-nonequilibrium Green function study

Motivated by the recent experiment on spin pumping from sub-THz radiation-driven uniaxial antiferromagnetic insulator (AFI) MnF 2 into heavy metal (HM) Pt hosting strong spin-orbit (SO) coupling, we compute and compare pumped spin currents in Cu/MnF 2 /Cu and Pt/MnF 2 /Cu heterostructures. Recent theories of spin pumping by AFI have relied on simplistic Hamiltonians (such as tight-binding) and the scattering approach to quantum transport yielding the so-called interfacial spin mixing conductance (SMC), but the concept of SMC ceases to be applicable when SO coupling is present directly at the interface. In contrast, we use a more general first-principles quantum transport approach which combines noncollinear density functional theory with Floquet-nonequilibrium Green's functions in order to take into account: SO-proximitized AFI as a new type of quantum material, different from isolated AFI and brought about by AFI hybridization with adjacent HM layer; strong SO coupling at the interface; and evanescent wavefunctions penetrating from Pt or Cu into AFI layer to make its interfacial region conducting rather than insulating as in the isolated AFI. The DC component of pumped spin current $I_\mathrm{DC}^{S_z}$ vs. precession cone angle $\theta_{\boldsymbol{l}}$ of the Néel vector l of AFI does not follow putative $I^{S_z}_\mathrm{DC} \propto \sin^2 \theta_{\boldsymbol{l}}$, except for very small angles $\theta_{\boldsymbol{l}} \lesssim 10^\circ$ for which we define an effective SMC from the prefactor and find that it doubles from MnF2/Cu to MnF2/Pt interface. In addition, the angular dependence $I^{S_z}_\mathrm{DC}(\theta_{\boldsymbol{l}})$ differs for opposite directions of precession of the Néel vector, leading to twice as large SMC for the right-handed than for the left-handed chirality of the precession mode.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Integrating subsystem embedding subalgebras and coupled cluster Green’s function: a theoretical foundation for quantum embedding in excitation manifold

Here, in this study, we introduce a novel approach to coupled-cluster Green's function (CCGF) embedding by seamlessly integrating conventional CCGF theory with the state-of-the-art sub-system embedding sub-algebras coupled cluster (SES-CC) formalism. This integration focuses primarily on delineating the characteristics of the sub-system and the corresponding segments of the Green's function, defined explicitly by active orbitals. Crucially, our work involves the adaptation of the SES-CC paradigm, addressing the left eigenvalue problem through a distinct form of Hamiltonian similarity transformation. This advancement not only facilitates a comprehensive representation of the interaction between the embedded sub-system and its surrounding environment but also paves the way for the quantum mechanical description of multiple embedded domains, particularly by employing the emergent quantum flow algorithms. Our theoretical underpinnings further set the stage for a generalization to multiple embedded sub-systems. This expansion holds significant promise for the exploration and application of non-equilibrium quantum systems, enhancing the understanding of system–environment interactions. In doing so, the research underscores the potential of SES-CC embedding within the realm of quantum computations and multi-scale simulations, promising a good balance between accuracy and computational efficiency.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Computation of the expectation value of the spin operator S^ 2 for the spin-flip Bethe–Salpeter equation

Spin-flip (SF) methods applied to excited-state approaches like the Bethe–Salpeter equation allow access to the excitation energies of open-shell systems, such as molecules and defects in solids. The eigenstates of these solutions, however, are generally not eigenstates of the spin operator S^ 2 . Even for simple cases where the excitation vector is expected to be, for example, a triplet state, the value of S^ 2 may be found to differ from 2.00; this difference is called 'spin contamination'. The expectation values S^ 2 must be computed for each excitation vector, to assist with the characterization of the particular excitation and to determine the amount of spin contamination of the state. Here, our aim is to provide for the first time in the SF methods literature a comprehensive resource on the derivation of the formulas for S^ 2 as well as its computational implementation. After a brief discussion of the theory of the SF Bethe–Salpeter equation (BSE) and some examples further illustrating the need for calculating S^ 2 , we present the derivation for the general equation for computing S^ 2 with the eigenvectors from an SF-BSE calculation, how it is implemented in a Python script, and timing information on how this calculation scales with the size of the SF-BSE Hamiltonian.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Exploring the nexus of many-body theories through neural network techniques: the tangent model

Abstract In this paper, we present a physically informed neural network (NN) representation of the effective interactions associated with coupled-cluster downfolding models to describe chemical systems and processes. The NN representation not only allows us to evaluate the effective interactions efficiently for various geometrical configurations of chemical systems corresponding to various levels of complexity of the underlying wave functions, but also reveals that the bare and effective interactions are related by a tangent function of some latent variables. We refer to this characterization of the effective interaction as a tangent model. We discuss the connection between this tangent model for the effective interaction with the previously developed theoretical analysis that examines the difference between the bare and effective Hamiltonians in the corresponding active spaces.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

The information content of projected galaxy fields

ABSTRACT The power spectrum of the non-linearly evolved large-scale mass distribution recovers only a minority of the information available on the mass fluctuation amplitude. We investigate the recovery of this information in 2D ‘slabs’ of the mass distribution averaged over ≈100 h−1 Mpc along the line of sight, as might be obtained from photometric redshift surveys. We demonstrate a Hamiltonian Monte Carlo method to reconstruct the non-Gaussian mass distribution in slabs, under the assumption that the projected field is a point-transformed Gaussian random field, Poisson-sampled by galaxies. When applied to the Quijote N-body suite at z = 0.5 and at a transverse resolution of 2 h−1 Mpc, the method recovers ∼30 times more information than the 2D power spectrum in the well-sampled limit, recovering the Gaussian limit on information. At a more realistic galaxy sampling density of 0.01 h3 Mpc−3, shot noise reduces the information gain to a factor of 5 improvement over the power spectrum at resolutions of 4 h−1 Mpc or smaller.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Differentiable stochastic halo occupation distribution

ABSTRACT In this work, we demonstrate how differentiable stochastic sampling techniques developed in the context of deep reinforcement learning can be used to perform efficient parameter inference over stochastic, simulation-based, forward models. As a particular example, we focus on the problem of estimating parameters of halo occupation distribution (HOD) models that are used to connect galaxies with their dark matter haloes. Using a combination of continuous relaxation and gradient re-parametrization techniques, we can obtain well-defined gradients with respect to HOD parameters through discrete galaxy catalogue realizations. Having access to these gradients allows us to leverage efficient sampling schemes, such as Hamiltonian Monte Carlo, and greatly speed up parameter inference. We demonstrate our technique on a mock galaxy catalogue generated from the Bolshoi simulation using a standard HOD model and find near-identical posteriors as standard Markov chain Monte Carlo techniques with an increase of ∼8× in convergence efficiency. Our differentiable HOD model also has broad applications in full forward model approaches to cosmic structure and cosmological analysis.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

DAmodel: hierarchical Bayesian modelling of DA white dwarfs for spectrophotometric calibration

We use hierarchical Bayesian modelling to calibrate a network of 32 all-sky faint DA white dwarf (DA WD) spectrophotometric standards (⁠16.5 < V , 19.5⁠) alongside three CALSPEC standards, from 912 Å to 32 μm. The framework is the first of its kind to jointly infer photometric zero points and WD parameters (surface gravity log g⁠, effective temperature T eff ⁠, extinction A V ⁠, dust relation parameter R V ) by simultaneously modelling both photometric and spectroscopic data. We model panchromatic Hubble Space Telescope Wide Field Camera 3 (HST/WFC3) UVIS and IR photometry, HST/STIS UV spectroscopy, and ground-based optical spectroscopy to sub-per cent precision. Photometric residuals for the sample are the lowest yet yielding < 0.004 mag RMS on average from the UV to the NIR, achieved by jointly inferring time-dependent changes in system sensitivity and WFC3/IR count-rate nonlinearity. Our GPU-accelerated implementation enables efficient sampling via Hamiltonian Monte Carlo, critical for exploring the high-dimensional posterior space. The hierarchical nature of the model enables population analysis of intrinsic WD and dust parameters. Inferred spectral energy distributions from this model will be essential for calibrating the James Webb Space Telescope as well as next-generation surveys, including Vera Rubin Observatory’s Legacy Survey of Space and Time and the Nancy Grace Roman Space Telescope.

methods: statistical↗
DOE OSTI2022

Negative string tension of a higher-charge Schwinger model via digital quantum simulation

We study some properties of generalized global symmetry for the charge- q Schwinger model in the Hamiltonian formalism, which is the (1 + 1)D quantum electrodynamics with a charge- q Dirac fermion. This model has the |$\mathbb {Z}_q\, 1$| -form symmetry, which is a remnant of the electric |$U(1)\, 1$| -form symmetry in the pure Maxwell theory. It is known that, if we put the theory on closed space, then the Hilbert space is decomposed into q distinct sectors, called universes, and some states with higher energy density do not decay to the ground state due to the selection rule of the 1-form symmetry. Even with open boundaries, we can observe the stability of such states by seeing a negative string tension behavior, meaning that opposite charges repel each other. In order to see negative string tensions, the vacuum angle θ has to be large enough and the standard path-integral Monte Carlo method suffers from the sign problem. We develop a method based on the adiabatic state preparation to see this feature with digital quantum simulation and confirm it using a classical simulator of quantum devices. In particular, we measure the local energy density and see how it jumps between the inside and outside of the insertion of the probe charges. We explicitly see that the energy density inside is lower than that outside. This is a clear signature of the negative string tension.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

A detailed analysis of the dynamics of fast neutrino flavor conversions with scattering effects

We calculate the dynamics of fast neutrino flavor conversions with Boltzmann collisions of neutrino scatterings in a homogeneous system. We find the enhancement or suppression of the flavor conversions in various setups of the collision terms. We analyze the mechanism of fast flavor conversions based on the motion of polarization vectors in the cylindrical coordinate analogous to a pendulum motion. The phases of the all the polarization vectors synchronize in the linear evolution phase, and the phase deviation from the Hamiltonian governs the neutrino flavor conversions. In a non-linear regime of flavor conversions, the collision terms induce a spiral motion of the polarization vector and gradually make the phase space smaller. The collision terms align all of the polarization vectors, and the flavor conversions eventually settle into equilibrium when the distributions of neutrinos become isotropic. Though our current analysis does not fully clarify the non-linear phenomena of fast flavor conversions, the framework of the pendulum motion gives a new insight into this complicated phenomenon that will be helpful in further studies.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Bayesian differential programming for robust systems identification under uncertainty

This paper presents a machine learning framework for Bayesian systems identification from noisy, sparse and irregular observations of nonlinear dynamical systems. The proposed method takes advantage of recent developments in differentiable programming to propagate gradient information through ordinary differential equation solvers and perform Bayesian inference with respect to unknown model parameters using Hamiltonian Monte Carlo sampling. This allows an efficient inference of the posterior distributions over plausible models with quantified uncertainty, while the use of sparsity-promoting priors enables the discovery of interpretable and parsimonious representations for the underlying latent dynamics. A series of numerical studies is presented to demonstrate the effectiveness of the proposed methods, including nonlinear oscillators, predator–prey systems and examples from systems biology. Taken together, our findings put forth a flexible and robust workflow for data-driven model discovery under uncertainty. All codes and data accompanying this article are available at https://bit.ly/34FOJMj .

Science & Technology - Other Topics↗
DOE OSTI2022

Gaussian processes meet NeuralODEs: a Bayesian framework for learning the dynamics of partially observed systems from scarce and noisy data

We present a machine learning framework (GP-NODE) for Bayesian model discovery from partial, noisy and irregular observations of nonlinear dynamical systems. The proposed method takes advantage of differentiable programming to propagate gradient information through ordinary differential equation solvers and perform Bayesian inference with respect to unknown model parameters using Hamiltonian Monte Carlo sampling and Gaussian Process priors over the observed system states. This allows us to exploit temporal correlations in the observed data, and efficiently infer posterior distributions over plausible models with quantified uncertainty. The use of the Finnish Horseshoe as a sparsity-promoting prior for free model parameters also enables the discovery of parsimonious representations for the latent dynamics. A series of numerical studies is presented to demonstrate the effectiveness of the proposed GP-NODE method including predator–prey systems, systems biology and a 50-dimensional human motion dynamical system. This article is part of the theme issue ‘Data-driven prediction in dynamical systems’.

Science & Technology - Other Topics↗
DOE OSTI2025

Exploring isospin symmetry breaking in exotic nuclei: High-precision mass measurement of 23 Si and shell-model calculations of 𝑇 = 5/2 nuclei

Here, we present a high-precision mass measurement of the proton-rich nucleus 23 Si, performed with the LEBIT Penning trap at the Facility for Rare Isotope Beams (FRIB) utilizing the time-of-flight ion cyclotron resonance (TOF-ICR) technique. We determined a mass excess of 23362.9(5.8) keV, which agrees with a recent storage-ring measurement from the experimental Cooler-Storage Ring (CSRe) in Lanzhou but has a factor of 20 improved precision 23 Si is hence the nucleus with the most precisely known mass among all nuclei with an isospin projection of 𝑇 𝑧 = −5/2. We performed shell-model calculations with the USDC and USDCm Hamiltonians to study binding energy differences and Thomas-Ehrmann shifts in mirror systems with an isospin up to 𝑇 = 5/2. Our experimental result and other recently reported masses of neutron-deficient sd-shell nuclei agree well with the theoretical predictions, demonstrating that isospin symmetry breaking in sd-shell nuclei—even at high isospin values—is well described by modern shell-model calculations.

20 ≤ A ≤ 38↗
DOE OSTI2025

Geometric Phase in Anisotropic Kepler Problem: Perspective for Realization in Rydberg Atoms

We predict a gyroscopic effect that can be demonstrated with Rydberg atoms following the dynamics of a Kepler Hamiltonian with an additional uniaxial anisotropy induced by optical ponderomotive force. This effect is analogous to the rotation of the Foucault pendulum in response to the Earth’s rotation. Furthermore, we argue that in Rydberg states with a large principal quantum number a similar geometric angle can be generated by mechanical rotations of an atomic-optical setup on timescales between 1 μ⁢s and 1 ms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Engineering micromotion in Floquet prethermalization via space-time symmetries

We present a systematic framework for Floquet prethermalization under strong resonant driving, emphasizing the pivotal role of dynamical space-time symmetries. Our approach demonstrates how dynamical space-time symmetries map onto the projective static symmetry group of the prethermal Hamiltonian governing the prethermal regime. We introduce techniques for detecting dynamical symmetries through the time evolution of local observables, facilitating a detailed analysis of micromotion within each period and surpassing the limitations of conventional stroboscopic Floquet prethermal dynamics. To implement this framework, we present a prethermal protocol that preserves order-2 dynamical symmetry in a spin-ladder model, confirming the predicted relationships between the expectation values of local observables at distinct temporal points in the Floquet cycle, linked by this symmetry.

floquet systems↗
DOE OSTI2026

Scattering Processes from Quantum Simulation Algorithms for Scalar Field Theories

We provide practical simulation methods for scalar field theories on a quantum computer that yield improved asymptotics as well as concrete gate estimates for the simulation and physical qubit estimates using the surface code. We achieve these improvements through two optimizations. First, we consider a finite volume approach for estimating the elements of the S-matrix. This approach is appropriate in general for 1+1D and for certain low-energy elastic collisions in higher dimensions. Second, we implement our approach using a series of different fault-tolerant simulation algorithms for Hamiltonians formulated both in the field occupation basis and field amplitude basis. Our algorithms are based on either second-order Trotterization or qubitization. The cost of Trotterization in occupation basis scales as O ( λ N 7 | Ω | 3 / ( M 5 / 2 ϵ 3 / 2 ) ) where λ is the coupling strength, N is the occupation cutoff, | Ω | is the volume of the spatial lattice, M is the mass of the particles and ϵ is the uncertainty in the energy calculation used for the S -matrix determination. Qubitization in the field basis scales as O ( | Ω | 2 ( k 2 Λ + k M 2 ) / ϵ ) , where k is the cutoff in the field and Λ is a scaled coupling constant. We find in both cases that the bounds suggest physically meaningful simulations can be performed using on the order of 4 × 10 6 physical qubits and 10 12 T -gates which corresponds to roughly one day on a superconducting quantum computer with surface code and a cycle time of 100 ns. This places the simulation of scalar field theory within striking distance of the gate counts for the best available chemistry simulation results.

Hardy, Andrew [Toronto U.] (ORCID:0000000235817382↗