Search NASA⌕ Search

SEARCH · Search NASA

Results for “Hamiltonian evolution”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Robust negativity in the quantum-to-classical transition of Kerr dynamics

Here, we quantify the quantum-to-classical transition of the single-mode Kerr nonlinear dynamics in the presence of loss. We establish three timescales that govern the dynamics, each with distinct characteristics. For times short compared with the Ehrenfest time, the evolution is classical, characterized by Gaussian dynamics. For sufficiently long times, as we increase the initial photon number, unitary Kerr evolution would generate macroscopic superpositions of coherent states (so-called kitten states). However, this is severely restricted in the presence of small photon loss, and the expectation values of observables coincide with their classical values. The intermediate timescale, however, shows resilient quantum behavior in the macroscopic limit. We show that in the mean-field non-Gaussian regime, the Kerr Hamiltonian (with small photon loss) generates a significant amount of Wigner-negativity, and classical flow is recovered only if the loss rate grows with system size. Our results broaden the usual understanding of quantum-to-classical transitions and demonstrate the potential for creating robust nonclassical resources for continuous-variable quantum information processing in the presence of loss.

Raza, Mohsin [University of New Mexico, Albuquerqu↗
DOE OSTI2022

Dynamical preparation of an atomic condensate in a Hofstadter band

The creation of a Hamiltonian in the quantum regime which has nontrivial topological features is a central goal of the cold-atom community, enabling widespread exploration of novel phases of quantum matter. A general scheme to synthesize such Hamiltonians is based on dynamical modulation of optical lattices which thereby generate vector potentials. At the same time, the modulation can lead to heating and serious difficulties with equilibration. Here we show that these challenges can be overcome by demonstrating how a Hofstadter Bose-Einstein condensate (BEC) can be dynamically realized, using experimental protocols. From Gross-Pitaevskii simulations our study reveals a complex, multistage evolution; this includes a chaotic intermediate "heating" stage followed by a spontaneous reentrance to the BEC. The observed behavior is reminiscent of evolution in cosmological models.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS1974

Possible histories of the obliquity of Mercury

Possible evolutionary paths of the spin axis of Mercury are studied theoretically in terms of generalized Cassini's laws using the Hamiltonian approach. Three stable positions (Cassini states) of the spin axis are discerned. The Cassini states in which the spin axis remains fixed in the frame rotating with the orbit precession are extremes in the tide-free Hamiltonian under the condition of constant angular momentum. The stability of the 3/2 spin resonance for Mercury as a function of obliquity is discussed, and it is shown how control of the resonance is passed from one to another term in the averaged potential as the obliquity is increased with a corresponding loss of stability at each transfer. A sampling technique is used to infer the evolution due to tidal interaction. Two tidal models are used - one with a constant value of the specific dissipation function Q and the other with Q inversely proportional to frequency. An approximate analytic solution for precession trajectories very close to a Cassini state is developed.

Peale, S. J.↗
DOE OSTI2024

Thermal evolution of spin excitations in honeycomb Ising antiferromagnetic FePSe 3

We use elastic and inelastic neutron scattering (INS) to study the antiferromagnetic (AF) phase transitions and spin excitations in the two-dimensional (2D) zig-zag antiferromagnet FePSe 3 . By determining the magnetic order parameter across the AF phase transition, we conclude that the AF phase transition in FePSe 3 is first-order in nature. In addition, our INS measurements reveal that the spin waves in the AF ordered state have a large easy-axis magnetic anisotropy gap, consistent with an Ising Hamiltonian, and possible biquadratic magnetic exchange interactions. On warming across T N , we find that dispersive spin excitations associated with three-fold rotational symmetric AF fluctuations change into FM spin fluctuations above T N . These results suggest that the first-order AF phase transition in FePSe 3 may arise from the competition between C 3 symmetric AF and C 1 symmetric FM spin fluctuations around T N , in place of a conventional second-order AF phase transition.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Alternative approach to quantum imaginary time evolution

There is increasing interest in quantum algorithms (QAs) that are based on the imaginary time evolution (ITE), a successful classical numerical approach to obtain ground states. However, most of the proposals so far require heavy postprocessing computational steps on a classical computer, such as solving linear equations. Here we provide an alternative approach to implement ITE. A key feature in our approach is the use of an orthogonal basis set: the propagated state is efficiently expressed in terms of orthogonal basis states at every step of the evolution. We argue that the number of basis states needed at those steps to achieve an accurate solution can be kept on the order of n , the number of qubits, by controlling the precision (number of significant digits) and the imaginary time increment. The number of quantum gates per imaginary time step is estimated to be polynomial in n . Additionally, while in many QAs the locality of the Hamiltonian is a key assumption, in our algorithm this restriction is not required. This characteristic of our algorithm renders it useful for studying highly nonlocal systems, such as the occupation-representation nuclear shell model. Here, we illustrate our algorithm through numerical implementation on an IBM quantum simulator.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Foundations of variational discrete action theory

Variational wave functions and Green's functions are two important paradigms for solving quantum Hamiltonians, each having their own advantages. Here we detail the variational discrete action theory (VDAT), which exploits the advantages of both paradigms in order to approximately solve the ground state of quantum Hamiltonians. VDAT consists of two central components: the sequential product density matrix (SPD) ansatz and a discrete action associated with the SPD. The SPD is a variational ansatz inspired by the Trotter decomposition and characterized by an integer $\mathscr{N}$, recovering many well-known variational wave functions, in addition to the exact solution for $\mathscr{N}$ = ∞. The discrete action describes all dynamical information of an effective integer time evolution with respect to the SPD. We generalize the path integral to our integer time formalism, which converts a dynamic correlation function in integer time to a static correlation function in a compound space. We also generalize the usual many-body Green's function formalism to integer time, which results in analogous but distinct mathematical structures, yielding integer time versions of the generating functional, Dyson equation, and Bethe-Salpeter equation. We prove that the SPD can be exactly evaluated in the multiband Anderson impurity model (AIM) by summing a finite number of diagrams. For the multiband Hubbard model, we prove that the self-consistent canonical discrete action approximation (SCDA), which is the integer time analog of the dynamical mean-field theory, exactly evaluates the SPD for d = ∞. VDAT within the SCDA provides an efficient yet reliable method for capturing the local physics of quantum lattice models, which will have broad applications for strongly correlated electron materials. More generally, VDAT should find applications in various many-body problems in physics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Preparations for quantum simulations of quantum chromodynamics in 1+1 dimensions. II. Single-baryon β-decay in real time

Here, a framework for quantum simulations of real-time weak decays of hadrons and nuclei in a two-flavor lattice theory in one spatial dimension is presented. A single generation of the Standard Model is found to require 16 qubits per spatial lattice site after mapping to spin operators via the Jordan-Wigner transformation. Both quantum chromodynamics and flavor-changing weak interactions are included in the dynamics, the latter through four-Fermi effective operators. Quantum circuits that implement time evolution in this lattice theory are developed and run on Quantinuum’s H1-1 20-qubit trapped ion system to simulate the β-decay of a single baryon on one lattice site. These simulations include the initial state preparation and are performed for both one and two Trotter time steps. The potential intrinsic error-correction properties of this type of lattice theory are discussed, and the leading lattice Hamiltonian required to simulate 0νββ-decay of nuclei induced by a neutrino Majorana mass term is provided.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1993

Alternative descriptions of wave and particle aspects of the harmonic oscillator

The dynamical properties of the wave and particle aspects of the harmonic oscillator can be studied with the help of the time-dependent Schroedinger equation (SE). Especially the time-dependence of maximum and width of Gaussian wave packet solutions allow to show the evolution and connections of those two complementary aspects. The investigation of the relations between the equations describing wave and particle aspects leads to an alternative description of the considered systems. This can be achieved by means of a Newtonian equation for a complex variable in connection with a conservation law for a nonclassical angular momentum-type quantity. With the help of this complex variable, it is also possible to develop a Hamiltonian formalism for the wave aspect contained in the SE, which allows to describe the dynamics of the position and momentum uncertainties. In this case the Hamiltonian function is equivalent to the difference between the mean value of the Hamiltonian operator and the classical Hamiltonian function.

Schuch, Dieter↗
DOE OSTI2021

Strong-field ionization of the triplet ground state of O 2

Using strong-field ionization as a probe, we observe highly nonperiodic evolution of the spin-rotation wave packet launched by a nonionizing femtosecond pulse in oxygen. The nonperiodicity is readily apparent only in rotationally cold molecules that are pumped with a weak alignment pulse. We show that this behavior is a consequence of the spin-rotation and the spin-spin couplings in the triplet ground state of the neutral molecule. A model that includes these couplings in the field-free Hamiltonian but in neither the alignment nor the ionization step explains most of the observed dynamics, suggesting that neither process depends explicitly on the electronic spin. We also show that the angle dependence of strong-field ionization can be retrieved from the delay-dependent signal even when coupling to spin complicates the rotational dynamics.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Testing isospin symmetry breaking in ab initio nuclear theory

In this work we present the first steps towards benchmarking isospin symmetry breaking in ab initio nuclear theory for calculations of superallowed Fermi β decay. Using the valence-space in-medium similarity renormalization group, we calculate b and c coefficients of the isobaric multiplet mass equation, starting from two different Hamiltonians constructed from chiral effective field theory. We compare results to experimental measurements for all T = 1 isobaric analog triplets of relevance to superallowed β decay for masses A = 10 to A = 74 and find an overall agreement within approximately 250 keV of experimental data for both b and c coefficients. A greater level of accuracy, however, is obtained by a phenomenological Skyrme interaction or a classical charged-sphere estimate. Lastly, we show that evolution of the valence-space operator does not meaningfully improve the quality of the coefficients with respect to experimental data, which indicates that higher-order many-body effects are likely not responsible for the observed discrepancies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Control optimization for parametric Hamiltonians by pulse reconstruction

Optimal control techniques provide a means to tailor the control pulses required to generate customized quantum gates, which helps to improve the resilience of quantum simulations to gate errors and device noise. However, the significant amount of (classical) computation required to generate customized gates can quickly undermine the effectiveness of this approach, especially when pulse optimization needs to be iterated. We propose a method to reduce the computational time required to generate the control pulse for a Hamiltonian that is parametrically dependent on a time-varying quantity. We use simple interpolation schemes to accurately reconstruct the control pulses from a set of pulses obtained in advance for a discrete set of predetermined parameter values. We obtain a reconstruction with very high fidelity and a significant reduction in computational effort. We report the results of the application of the proposed method to device-level quantum simulations of the unitary (real) time evolution of two interacting neutrons based on superconducting qubits.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Basis light-front quantization for the Λ b and Σ b baryons

Within the basis light-front quantization framework, we compute the masses and light-front wave functions of the Λ b baryon and its isospin triplet counterparts Σ b + , Σ b 0 , and Σ b − using a light-front effective Hamiltonian in the leading Fock sector. These wave functions are obtained as eigenstates of the effective Hamiltonian, which incorporates the one-gluon exchange interaction with fixed coupling and a three-dimensional confinement potential. With the quark masses and the couplings as adjustable parameters, the computed masses are set within the experimental range. The resulting predictions for their electromagnetic properties align well with other theoretical calculations. Additionally, the parton distribution functions (PDFs) of these baryons are obtained for the first time, with gluon and sea quark distributions dynamically generated through QCD evolution of the valence quark PDFs. Published by the American Physical Society 2025

Meng, Lingdi (ORCID:0009000096439904)↗
DOE OSTI2024

Hamiltonian switching control of noisy bipartite qubit systems

Abstract We develop a Hamiltonian switching ansatz for bipartite control that is inspired by the quantum approximate optimization algorithm, to mitigate environmental noise on qubits. We demonstrate the control for a central spin coupled to bath spins via isotropic Heisenberg interactions, and then make physical applications to the protection of quantum gates performed on superconducting transmon qubits coupling to environmental two-level-systems (TLSs) through dipole-dipole interactions, as well as on such qubits coupled to both TLSs and a Lindblad bath. The control field is classical and acts only on the system qubits. We use reinforcement learning with policy gradient to optimize the Hamiltonian switching control protocols, using a fidelity objective for specific target quantum gates. We use this approach to demonstrate effective suppression of both coherent and dissipative noise, with numerical studies achieving target gate implementations with fidelities over 0.9999 (four nines) in the majority of our test cases and showing improvement beyond this to values of 0.999 999 999 (nine nines) upon a subsequent optimization by GRadient Ascent Pulse Engineering (GRAPE). We analyze how the control depth, total evolution time, number of environmental TLS, and choice of optimization method affect the fidelity achieved by the optimal protocols and reveal some critical behaviors of bipartite control of quantum gates.

Physics↗
DOE OSTI2021

Theory of Trotter Error with Commutator Scaling

The Lie-Trotter formula, together with its higher-order generalizations, provides a simple approach to decomposing the exponential of a sum of operators. Despite significant effort, the error scaling of such product formulas remains poorly understood. We develop a theory of Trotter error that overcomes the limitations of truncating the Baker-Campbell-Hausdorff expansion. Our analysis directly exploits the commutativity of operator summands, producing tighter error bounds for both real- and imaginary-time evolutions. Whereas previous work achieves similar goals for systems with geometric locality or Lie-algebraic structure, our approach holds in general. We give a host of improved algorithms for digital quantum simulation and quantum Monte Carlo methods, nearly matching or even outperforming the best previous results. Our applications include: (i) a simulation of second-quantized plane-wave electronic structure, nearly matching the interaction-picture algorithm of Low and Wiebe; (ii) a simulation of $k$-local Hamiltonians almost with induced one-norm scaling, faster than the qubitization algorithm of Low and Chuang; (iii) a simulation of rapidly decaying power-law interactions, outperforming the Lieb-Robinson-based approach of Tran et al.; (iv) a hybrid simulation of clustered Hamiltonians, dramatically improving the result of Peng, Harrow, Ozols, and Wu; and (v) quantum Monte Carlo simulations of the transverse field Ising model and quantum ferromagnets, tightening previous analyses of Bravyi and Gosset. We obtain further speedups using the fact that product formulas can preserve the locality of the simulated system. Specifically, we show that local observables can be simulated with complexity independent of the system size for power-law interacting systems, which implies a Lieb-Robinson bound nearly matching a recent result of Tran et al. Our analysis reproduces known tight bounds for first- and second-order formulas. We further investigate the tightness of our bounds for higher-order formulas. For quantum simulation of a one-dimensional Heisenberg model with an even-odd ordering of terms, our result overestimates the complexity by only a factor of $5$. Our bound is also close to tight for power-law interactions and other orderings of terms. This suggests that our theory can accurately characterize Trotter error in terms of both the asymptotic scaling and the constant prefactor.

quantum computing, numerical analysis↗
DOE OSTI2025

Triaxiality and shape dynamics in 70 Ge

The electromagnetic properties of low-lying states in 70 Ge were investigated via multi-step Coulomb excitation of a 70 Ge beam impinging on a 208 Pb target at the ATLAS facility of the Argonne National Laboratory. A total of 27 transitional elements and six diagonal matrix elements coupling 11 low-lying states, were extracted from the measured cross sections. These were used to calculate reduced transition probabilities, spectroscopic quadrupole moments, and rotational invariant shape parameters, providing enhanced precision and expanding on previous studies. The experimental data were compared within several theoretical frameworks, including the generalized triaxial rotor model, configuration interaction shell-model calculations, and computations within the combined frameworks of relativistic density functional theory and the five-dimensional collective Hamiltonian. The results demonstrate a good agreement with the experimental data and, in conjunction with calculations using a two-state mixing model, support significant triaxiality and strong mixing between the 0$^{+}_{1}$ and 0$^{+}_{2}$ states. This results in the magnitudes of their respective quadrupole deformations [β rms (0$^{+}_{1}$) = 0.228 (3), β rms (0$^{+}_{2}$) = 0.273 (1)] being more similar than previously observed. Furthermore, the implications of these results for understanding the complex shape coexistence phenomena, the role of triaxiality, and shape evolution along the Ge isotopic chain are discussed.

59 ≤ A ≤ 89↗
DOE OSTI2023

Quantum computing for fusion energy science applications

This is a review of recent research exploring and extending present-day quantum computing capabilities for fusion energy science applications. We begin with a brief tutorial on both ideal and open quantum dynamics, universal quantum computation, and quantum algorithms. Then, we explore the topic of using quantum computers to simulate both linear and nonlinear dynamics in greater detail. Because quantum computers can only efficiently perform linear operations on the quantum state, it is challenging to perform nonlinear operations that are generically required to describe the nonlinear differential equations of interest. In this work, we extend previous results on embedding nonlinear systems within linear systems by explicitly deriving the connection between the Koopman evolution operator, the Perron–Frobenius evolution operator, and the Koopman–von Neumann evolution (KvN) operator. We also explicitly derive the connection between the Koopman and Carleman approaches to embedding. Extension of the KvN framework to the complex-analytic setting relevant to Carleman embedding, and the proof that different choices of complex analytic reproducing kernel Hilbert spaces depend on the choice of Hilbert space metric are covered in the appendixes. Finally, we conclude with a review of recent quantum hardware implementations of algorithms on present-day quantum hardware platforms that may one day be accelerated through Hamiltonian simulation. We discuss the simulation of toy models of wave–particle interactions through the simulation of quantum maps and of wave–wave interactions important in nonlinear plasma dynamics.

Joseph, I. (ORCID:0000000255400840)↗
DOE OSTI2022

On the S-matrix of Ising field theory in two dimensions

We explore the analytic structure of the non-perturbative S-matrix in arguably the simplest family of massive non-integrable quantum field theories: the Ising field theory (IFT) in two dimensions, which may be viewed as the Ising CFT deformed by its two relevant operators, or equivalently, the scaling limit of the Ising model in a magnetic field. Our strategy is that of collider physics: we employ Hamiltonian truncation method (TFFSA) to extract the scattering phase of the lightest particles in the elastic regime, and combine it with S-matrix bootstrap methods based on unitarity and analyticity assumptions to determine the analytic continuation of the 2 → 2 S-matrix element to the complex s-plane. Focusing primarily on the “high temperature” regime in which the IFT interpolates between that of a weakly coupled massive fermion and the E 8 affine Toda theory, we will numerically determine 3-particle amplitudes, follow the evolution of poles and certain resonances of the S-matrix, and exclude the possibility of unknown wide resonances up to reasonably high energies.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗