Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “EIGENSTATES”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Quantum chaos, integrability, and late times in the Krylov basis

Quantum chaotic systems are conjectured to display a spectrum whose fine-grained features (gaps and correlations) are well described by random matrix theory (RMT). We propose and develop a complementary version of this conjecture: quantum chaotic systems display a Lanczos spectrum whose local means and covariances are well described by RMT. To support this proposal, we first demonstrate its validity in examples of chaotic and integrable systems. We then show that for Haar-random initial states in RMTs the mean and covariance of the Lanczos spectrum suffice to produce the full long-time behavior of general survival probabilities including the spectral form factor, as well as the spread complexity. In addition, for initial states with continuous overlap with energy eigenstates, we analytically find the long-time averages of the probabilities of Krylov basis elements in terms of the mean Lanczos spectrum. This analysis suggests a notion of eigenstate complexity, the statistics of which differentiate integrable systems and classes of quantum chaos. Lastly, we clarify the relation between spread complexity and the universality classes of RMT by exploring various values of the Dyson index and Poisson distributed spectra.

Combinatorics↗
DOE OSTI2021

Emergent Ergodicity at the Transition between Many-Body Localized Phases

Strongly disordered systems in the many-body localized (MBL) phase can exhibit ground state order in highly excited eigenstates. The interplay between localization, symmetry, and topology has led to the characterization of a broad landscape of MBL phases ranging from spin glasses and time crystals to symmetry protected topological phases. Understanding the nature of phase transitions between these different forms of eigenstate order remains an essential open question. In this work, we conjecture that no direct transition between distinct MBL orders can occur in one dimension; rather, an ergodic phase always intervenes. Motivated by recent advances in Rydberg-atom-based quantum simulation, we propose an experimental protocol where the intervening ergodic phase can be diagnosed via the dynamics of local observables.

1-dimensional spin chains↗
DOE OSTI2021

Rodeo Algorithm for Quantum Computing

We present a stochastic quantum computing algorithm that can prepare any eigenvector of a quantum Hamiltonian within a selected energy interval $\ [E-\epsilon, E+\epsilon]$. In order to reduce the spectral weight of all other eigenvectors by a suppression factor δ, the required computational effort scales as $\ O[|\log \delta|/(p \epsilon)]$, where p s the squared overlap of the initial state with the target eigenvector. The method, which we call the rodeo algorithm, uses auxiliary qubits to control the time evolution of the Hamiltonian minus some tunable parameter E . In this manner, we converge to the target eigenvector with exponential accuracy in the number of measurements. In addition to preparing eigenvectors, the method can also compute the full spectrum of the Hamiltonian. We illustrate the performance with several examples. For energy eigenvalue determination with error $\epsilon$, the computational scaling is $\ O[(\log \epsilon)^2/(p \epsilon)]$. For eigenstate preparation, the computational scaling is $\ O(\log \Delta/p)$, where $\Delta$ is the magnitude of the orthogonal component of the residual vector. The speed for eigenstate preparation is exponentially faster than that for phase estimation or adiabatic evolution.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Exotic edge states of C 3 high-fold fermions in honeycomb lattices

A generalization of the graphene honeycomb model to the case where each site in the honeycomb lattice contains a n -fold degenerate set of eigenstates of the C 3 symmetry has been recently proposed to describe several systems, including triangulene crystals and photonic lattices. These generalized honeycomb models are defined by ( n a , n b ) , the number of C 3 eigenstates in the a and b sites of the unit cell, resulting in n a + n b bands. Thus, the (1,1) case gives the coventional honeycomb model that describes the two low-energy bands in graphene. Generalizations, such as (2,1), (2,2), and (3,3) display several nontrivial features, such as coexisting graphenelike Dirac cones with flat bands, both at zero and finite energy, as well as robust degeneracy points where a flat band and a parabolic band meet at the Γ point. Here we explore the edge states of this class of crystals, using as reference triangulene crystals, and we find several types of edge states absent in the conventional (1,1) honeycomb case, associated to the nontrivial features of the two-dimensional bands of the high-fold case. First, we find dispersive edge states associated to the finite-energy flat bands, that occur both at the armchair and zigzag termination. Second, in the case of noncentrosymmetric triangulene crystals that lead to a S = 1 Dirac band, we have a bonding-antibonding pair of dispersive edge states, localized in the same edge so that their energy splitting is reduced as their localization increases, opposite to the conventional behavior of pairs of states localized in opposite edges. Third, for the (3,3) case, that hosts a gap separating a pair of flat conduction and valence bands, we find nondispersive edge states with E = 0 in all edge terminations. Published by the American Physical Society 2024

Madail, L. (ORCID:0000000208457748)↗
DOE OSTI2022

Improved measurement of the strong-phase difference $\delta _D^{K\pi }$ in quantum-correlated $D{\overline{D}}$ decays

The decay D → K - π + is studied in a sample of quantum-correlated D$\overline{D}$ pairs, based on a data set corresponding to an integrated luminosity of 2.93 fb -1 collected at the ψ(3770) resonance by the BESIII experiment. The asymmetry between CP-odd and CP-even eigenstate decays into K - π + is determined to be A Kπ = 0.132 ± 0.011 ± 0.007, where the first uncertainty is statistical and the second is systematic. This measurement is an update of an earlier study exploiting additional tagging modes, including several decay modes involving a $K^{0}_{L}$ meson. The branching fractions of the $K^{0}_{L}$ modes are determined as input to the analysis in a manner that is independent of any strong phase uncertainty. Using the predominantly CP-even tag D → π + π - π 0 and the ensemble of CP-odd eigenstate tags, the observable $A^{πππ^0}_{Kπ}$ is measured to be 0.130 ± 0.012 ± 0.008. The two asymmetries are sensitive to $r^{Kπ}_{D}$ cos $δ^{Kπ}_{D}$ where $r^{Kπ}_{D}$ and $δ^{Kπ}_{D}$ are the ratio of amplitudes and phase difference, respectively, between the doubly Cabibbo-suppressed and Cabibbo-favoured decays. In addition, events containing D → K - π + tagged by D → $K^{0}_{S,L}$π + π - are studied in bins of phase space of the three-body decays. This analysis has sensitivity to both $r^{Kπ}_{D}$ cos $δ^{Kπ}_{D}$ and $r^{Kπ}_{D}$, sin $δ^{Kπ}_{D}$. A fit to A Kπ , $A^{πππ^0}_{Kπ}$ and the phase-space distribution of the D → $K^{0}_{S,L}$π + π - tags yields $δ^{Kπ}_{D}$ = ($187.6^{+8.9+5.4}_{-9.7-6.4}$)°, where external constraints are applied for $r^{Kπ}_{D}$ and other relevant parameters. This is the most precise measurement of $δ^{Kπ}_{D}$ in quantum-correlated D$\overline{D}$ decays.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Enhancing scalability of a matrix-free eigensolver for studying many-body localization

We propose several techniques to enhance the parallel scalability of a matrix-free eigensolver designed for studying many-body localization (MBL) of quantum spin chain models with nearest-neighbor interactions and on-site disorder. This type of problem is computationally challenging because the dimension of the associated Hamiltonian matrix grows exponentially with respect to the number of spins L, and we need to average over different realizations of the random disorder to obtain relevant statistical behavior. For each disorder realization, we need to compute eigenvalues from different regions of the spectrum and their corresponding eigenvectors. In previous work, the interior eigenstates for a single eigenvalue problem are computed via the shift-and-invert Lanczos algorithm. Due to the extremely high memory footprint of the LU factorizations, this technique is not well suited for large L’s. For example, we need thousands of compute nodes on modern high performance computing infrastructures to go beyond L = 24. The matrix-free approach does not suffer from this memory bottleneck, however, its scalability is limited by a computation and communication load imbalance. To reduce this imbalance and to significantly enhance the scalability of the matrix-free eigensolver, we reorder the matrix and leverage the consistent space runtime, CSPACER. We also show its efficiency in managing irregular communication patterns at scale compared to optimized MPI non-blocking two-sided and one-sided RMA implementation variants. This effort enables us to study MBL for spin chains with a larger number of spins. The efficiency and effectiveness of the proposed algorithm is demonstrated by computing eigenstates on a massively parallel many-core high performance computer.

METIS↗
DOE OSTI2023

A time-dependent momentum-resolved scattering approach to core-level spectroscopies

While new light sources allow for unprecedented resolution in experiments with X-rays, a theoretical understanding of the scattering cross-section is lacking. In the particular case of strongly correlated electron systems, numerical techniques are quite limited, since conventional approaches rely on calculating a response function (Kramers-Heisenberg formula) that is obtained from a perturbative analysis of scattering processes in the frequency domain. This requires a knowledge of a full set of eigenstates in order to account for all intermediate processes away from equilibrium, limiting the applicability to small tractable systems. In this work, we present an alternative paradigm, recasting the problem in the time domain and explicitly solving the time-dependent Schrödinger equation without the limitations of perturbation theory: a faithful simulation of the scattering processes taking place in actual experiments, including photons and core electrons. We show how this approach can yield the full time and momentum resolved Resonant Inelastic X-Ray Scattering (RIXS) spectrum of strongly interacting many-body systems. We demonstrate the formalism with an application to Mott insulating Hubbard chains using the time-dependent density matrix renormalization group method, which does not require a priory knowledge of the eigenstates and can solve very large systems with dozens of orbitals. This approach can readily be applied to systems out of equilibrium without modification and generalized to other spectroscopies.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Real-Time Krylov Theory for Quantum Computing Algorithms

Quantum computers provide new avenues to access ground and excited state properties of systems otherwise difficult to simulate on classical hardware. New approaches using subspaces generated by real-time evolution have shown efficiency in extracting eigenstate information, but the full capabilities of such approaches are still not understood. In recent work, we developed the variational quantum phase estimation (VQPE) method, a compact and efficient real-time algorithm to extract eigenvalues on quantum hardware. Here we build on that work by theoretically and numerically exploring a generalized Krylov scheme where the Krylov subspace is constructed through a parametrized real-time evolution, which applies to the VQPE algorithm as well as others. We establish an error bound that justifies the fast convergence of our spectral approximation. We also derive how the overlap with high energy eigenstates becomes suppressed from real-time subspace diagonalization and we visualize the process that shows the signature phase cancellations at specific eigenenergies. We investigate various algorithm implementations and consider performance when stochasticity is added to the target Hamiltonian in the form of spectral statistics. To demonstrate the practicality of such real-time evolution, we discuss its application to fundamental problems in quantum computation such as electronic structure predictions for strongly correlated systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Search for Heavy Neutral Leptons at the MINER$\nu$A detector

Heavy Neutral Leptons (HNL) are particles hypothesised to provide a mass generation mechanism for the active (observed) neutrino species, which are known to have nonzero mass from the definitive observation of neutrino oscillations. HNL are eigenstates of mass of the order $\mathcal{O}(0.1 − 1\,\,\mathrm{GeV}/c^{2})$, which mix into the active flavour eigenstates through the extended leptonic mixing matrix. Apart from neutrino mass, they could provide a natural dark matter candidate and a mechanism for matter-antimatter asymmetry in the early Universe, giving rise to today’s matter-dominated cosmos. Searches for HNL typically attempt to either confirm the existence of HNL through an excess in data that is most compatible with an HNL hypothesis, or by setting limits on the HNL parameter space $\left(M_{N4}, \left|U_{\alpha 4}\right|^{2}\right)$ in the case of statistically insignificant excess. Such searches are intensifying as part of a global research programme at both colliders and accelerator / atmospheric neutrino experiments. One such setting is the MINER$\nu$A experiment, located in the NuMI beamline at Fermilab. MINER$\nu$A has collected a large amount of data over seven years of operation to measure the cross-sections of neutrino-nucleus interactions necessary to drive systematic uncertainties down, in order for neutrino oscillation experiments to achieve sensitivity to CP violation in the neutrino sector. With the high-energy, high-intensity NuMI beam, and with good timing, position, and angular resolution leading to sensitivity to HNL decays, MINER$\nu$A is fertile ground for an HNL search, which is performed in this thesis. Specifically, a novel, experiment-agnostic and general simulation of HNL production and decay is presented and deployed; this simulation has been incorporated in the ubiquitous GENIE neutrino event generator for use with neutrino experiments in the future. The event selection and background characterisation is discussed in detail, including the primary background coming from charged-current coherent and diffractive pion production from neutrino-nucleus interactions, and constraints on the background using control regions from MINER$\nu$A data are derived. Finally, by means of a fake-data study given a background of $\mathcal{O}(500)$ events, the discovery potential and limit-setting capacity of MINER$\nu$A is demonstrated, and ways to improve this capacity are expounded upon.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Hilbert-space fragmentation, multifractality, and many-body localization

Investigating many-body localization (MBL) using exact numerical methods is limited by the exponential growth of the Hilbert space. However, localized eigenstates display multifractality and only extend over a vanishing fraction of the Hilbert space. Here, building on this remarkable property, we develop a simple yet efficient decimation scheme to discard the irrelevant parts of the Hilbert space of the random-field Heisenberg chain. This leads to a Hilbert space fragmentation in small clusters, allowing to access larger systems at strong disorder. The MBL transition is quantitatively predicted, together with a geometrical interpretation of MBL multifractality as a shattering of the Hilbert space.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Anisotropy of localized states in an anisotropic disordered medium

Highlights: • Mobility edges of the Anderson transition are unaffected by anisotropy alone. • The critical exponent of the Anderson localization transition is also unaffected. • Anderson localized eigenstates in an anisotropic medium retain substantial anisotropy. • Their anisotropy is weaker than expected from purely geometric considerations. • Modes with the longest lifetimes are found to be the most anisotropic. We study Anderson localization of a scalar wave in an ensemble of resonant point scatterers embedded in an anisotropic background medium. For uniaxial anisotropy of moderate strength, the mobility edges and the critical exponent of the localization transition are found to be unaffected by the anisotropy provided that the determinant of the anisotropy tensor is kept equal to one upon introducing the anisotropy. Localized modes have anisotropic spatial shapes although their anisotropy is weaker than the one expected from purely geometric considerations. The modes with the longest lifetimes are found to be the most anisotropic and their anisotropy increases with the size of the disordered medium.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Delocalization of a non-Hermitian quantum walk on random media in one dimension

Highlights: • We study the localization-delocalization transition of a non-Hermitian quantum walk. • We find that the phase transition is similar to the one in the Hatano-Nelson model. • All eigenvectors get extended and all eigenvalues become complex at the transition. • This implies that the localization lengths of all eigenvectors are the same. We first review the localization–delocalization transition of a non-Hermitian random tight-binding Anderson model, called the Hatano–Nelson model. We then report a new result for a non-Hermitian extension of a discrete-time quantum walk on a one-dimensional random medium; we numerically find a delocalization transition similar to one of the Hatano–Nelson model. As a common feature to both models, at the transition point, an eigenvector gets delocalized and at the same time the corresponding energy eigenvalue (for the latter quantum-walk model, the imaginary unit times the phase of the eigenvalue of the time-evolution operator) becomes complex. One of the unique properties of the present non-Hermitian quantum walk is that the localization length of all eigenvectors is the same, and thereby all eigenstates simultaneously undergo the delocalization transition and all energy eigenvalues become complex at the same time when we turn up a non-Hermitian parameter.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Adaptive variational quantum eigensolvers for highly excited states

Highly excited states of quantum many-body systems are central objects in the study of quantum dynamics and thermalization that challenge classical computational methods due to their volume-law entanglement content. In this work, we explore the potential of variational quantum algorithms to approximate such states. We propose an adaptive variational quantum eigensolver (VQE) for excited states (X) that self-generates a variational ansatz for arbitrary eigenstates of a many-body Hamiltonian H by attempting to minimize the energy variance with respect to H. We benchmark the method by applying it to an Ising spin chain with integrable and nonintegrable regimes, where we calculate various quantities of interest, including the total energy, magnetization density, and entanglement entropy. We also compare the performance of adaptive VQE-X to an adaptive variant of the folded-spectrum method. For both methods, we find a strong dependence of the algorithm's performance on the choice of operator pool used for the adaptive construction of the ansatz. In particular, an operator pool including long-range two-body gates accelerates the convergence of both algorithms in the nonintegrable regime. Here, we also study the scaling of the number of variational parameters with system size, finding that an exponentially large number of parameters may be necessary to approximate individual highly excited states. Nevertheless, we argue that these methods lay a foundation for the use of quantum algorithms to study finite-energy-density properties of many-body systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Quantum many-body scars from Einstein-Podolsky-Rosen states in bilayer systems

Quantum many-body scar states are special eigenstates of nonintegrable models with distinctive entanglement features that give rise to infinitely long-lived coherent dynamics under quantum quenches from certain initial states. Here, we elaborate on a construction of quantum many-body scar states in which they emerge from Einstein-Podolsky-Rosen states in systems with two layers, wherein the two layers are maximally entangled. We apply this construction to spin systems as well as systems of itinerant fermions and bosons and demonstrate how symmetries can be harnessed to enhance its versatility. We show that several well-known examples of quantum many-body scars, including the tower of states in the spin-1 XY model and the η-pairing states in the Fermi-Hubbard model, can be understood within this formalism. We also demonstrate how an infinite tower of many-body scar states can emerge in bilayer Bose-Hubbard models with charge conservation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Lindblad many-body scars

Quantum many-body scars have received much recent attention for being both intriguing nonergodic states in otherwise quantum chaotic systems and promising candidates to encode quantum information efficiently. So far, these studies have mostly been restricted to Hermitian systems. Here, we study many-body scars in many-body quantum chaotic systems coupled to a Markovian bath, which we term Lindblad many-body scars. They are defined as simultaneous eigenvectors of the Hamiltonian and dissipative parts of the vectorized Liouvillian. Importantly, because their eigenvalues are purely real, they are not related to revivals. The number and nature of the scars depend on both the symmetry of the Hamiltonian and the choice of jump operators. For a dissipative four-body Sachdev-Ye-Kitaev (SYK) model with 𝑁 fermions, either Majorana or complex, we construct analytically some of these Lindblad scars while others could only be obtained numerically. As an example of the former, we identify 𝑁/2+1 scars for complex fermions due to the 𝑈⁡(1) symmetry of the model and two scars for Majorana fermions as a consequence of the parity symmetry. Similar results are obtained for a dissipative XXZ spin chain. We also characterize the physical properties of Lindblad scars. First, the operator size is independent of the disorder realization and has a vanishing variance. By contrast, the operator size for nonscarred states, believed to be quantum chaotic, is well described by a distribution centered around a specific size and a finite variance, which could be relevant for a precise definition of the eigenstate thermalization hypothesis in dissipative quantum chaos. Moreover, the entanglement entropy of these scars has distinct features such as a strong dependence on the partition choice and, in certain cases, a large entanglement.

Eigenstate thermalization↗
DOE OSTI2022

Three-point functions in $\mathrm{ABJM}$ and Bethe Ansatz

We develop an integrability-based framework to compute structure constants of two sub-determinant operators and a single-trace non-BPS operator in ABJM theory in the planar limit. In this first paper, we study them at weak coupling using a relation to an integrable spin chain. We first develop a nested Bethe ansatz for an alternating SU(4) spin chain that describes single-trace operators made out of scalar fields. We then apply it to the computation of the structure constants and show that they are given by overlaps between a Bethe eigenstate and a matrix product state. We conjecture that the determinant operator corresponds to an integrable matrix product state and present a closed-form expression for the overlap, which resembles the so-called Gaudin determinant. We also provide evidence for the integrability of general sub-determinant operators. The techniques developed in this paper can be applied to other quantities in ABJM theory including three-point functions of single-trace operators.

1/N Expansion↗
DOE OSTI2024

A principle of maximum ignorance for semiclassical gravity

The principle of maximum ignorance posits that the coarse-grained description of a system is maximally agnostic about its underlying microscopic structure. We briefly review this principle for random matrix theory and for the eigenstate thermalization hypothesis. We then apply this principle in holography to construct ensembles of random mixed states. This leads to an ensemble of microstates which models our microscopic ignorance, and which on average reproduces the effective semiclassical physics of a given bulk state. We call this ensemble the state-averaging ansatz. The output of our model is a prediction for semiclassical contributions to variances and higher statistical moments over the ensemble of microstates. The statistical moments provide coarse-grained — yet gravitationally non-perturbative — information about the microstructure of the individual states of the ensemble. We show that these contributions exactly match the on-shell action of known wormhole configurations of the gravitational path integral. These results strengthen the view that wormholes simply parametrize the ignorance of the microstructure of a fundamental state, given a fixed semiclassical bulk description.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Constraints on a generalization of geometric quantum mechanics from neutrino and B0-$$ \overline{B^0} $$ oscillations

Abstract Nambu Quantum Mechanics, proposed in Phys. Lett. B536, 305 (2002), is a deformation of canonical Quantum Mechanics in which the manifold over which the “phase” of an energy eigenstate time evolves is modified. This generalization affects oscillation and interference phenomena through the introduction of two deformation parameters that quantify the extent of deviation from canonical Quantum Mechanics. In this paper, we constrain these parameters utilizing atmospheric neutrino oscillation data, andB 0 -$$ \overline{B^0} $$ B 0 ¯ oscillation data from Belle. Surprisingly, the bound from atmospheric neutrinos is stronger than the bound from Belle. Various features of Nambu Quantum Mechanics are also discussed.

Physics↗