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 73 records · Page 4

iSPECTRON: a simulation interface for linear and nonlinear spectra with ab-initio quantum chemistry software

We introduce iSPECTRON, an open source (under the Educational Community License version 2.0) program that parses data from common quantum chemistry software (NWChem, OpenMolcas, Gaussian, Cobramm, etc.), produces the input files for the simulation of linear and nonlinear spectroscopy of molecules with the Spectron code, and analyzes the spectra with a broad range of tools. Vibronic spectra are expressed in term of the electronic eigenstates, obtained through quantum chemistry computations, and vibrational/bath effects are incorporated in the framework of the displaced harmonic oscillator model, where all required quantities are computed at the Franck-Condon point. The code capabilities are illustrated by simulating linear absorption, transient absorption and two dimensional electronic spectra of the pyrene molecule. Two levels of electronic structure theory, TDDFT (with NWChem) and RASSCF/RASPT2 (with OpenMolcas), are compared where possible. Acknowledgements: F.S., A.N., D.R.N., N.G., S.M, M.G. acknowledge support from the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences, Chemical Sciences, Geosciences, and Biosciences Division under Award Nos. DE-SC0019484, KC-030103172684. The Spectron code was developed with support from the National Science Foundation (Grant CHE- 1953045). This research benefited from computational resources provided by EMSL, a DOE Office of Science User Facility sponsored by the Office of Biological and Environmental Research and located at PNNL. PNNL is operated by Battelle Memorial Institute for the United States Department of Energy under DOE Contract No. DE-AC05-76RL1830.

Segatta, Francesco↗

The QCD axion sum rule

Abstract We demonstrate that the true QCD axion that solves the strong CP problem can be found in all generality outside the customary standard QCD band, with QCD being the sole source of Peccei-Quinn breaking. The essential reason is that the basis of axion-gluon interactions does not need to coincide with the mass basis. Specifically, we consider the case in which the QCD axion field is not the only singlet scalar in Nature but it mixes with other singlet scalars (besides theη′). We determine the exact mathematical condition for an arbitraryN-scalar potential to be Peccei-Quinn invariant. Such potentials provide extra sources of mass for the customary axion without enlarging the Standard Model gauge symmetry. The contribution to the axion mass stemming from the QCD topological susceptibility is shown to be shared then among theNaxion eigenstates through a precise sum rule. Their location can only be displaced to the right of the standard QCD band. We demonstrate that the axion closest to this band can be displaced from it by a factor of$$ \sqrt{N} $$ N at most, and this corresponds to the case in which all axion signals are maximally deviated. Conversely, if one axion is found on the standard QCD band, the other eigenstates will be out of experimental reach. Our results imply that any ALP experiment which finds a signal to the right of the standard QCD axion band can be solving the strong CP problem within QCD, with the associatedN− 1 excitations to be found in an area of parameter space that we determine. We illustrate the results and phenomenology in some particular cases.

Physics↗

Quantum chaos and unitary black hole evaporation

The formation and evaporation of small AdS black holes in a theory with a holographic dual is governed by the usual rules of quantum mechanics. The eigenstate thermalization hypothesis explains the validity of semiclassical gravity for local bulk observables and can be used to quantify the magnitude of quantum corrections to the semi-classical approximation. The holographic dual produces a basis of black hole states with finite energy width, and observables that are smooth functions on the classical phase space will self-average over a large number of energy eigenstates, exponential in the Bekenstein-Hawking entropyS, leading to results that are consistent with semiclassical gravity up to small corrections of order e –S/2 . As expected, the semiclassical description breaks down for transition amplitudes that reflect the unitary evolution of the holographic theory.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Celestial amplitudes from UV to IR

Celestial amplitudes represent 4D scattering of particles in boost, rather than the usual energy-momentum, eigenstates and hence are sensitive to both UV and IR physics. We show that known UV and IR properties of quantum gravity translate into powerful constraints on the analytic structure of celestial amplitudes. For example the soft UV behavior of quantum gravity is shown to imply that the exact four-particle scattering amplitude is meromorphic in the complex boost weight plane with poles confined to even integers on the negative real axis. Would-be poles on the positive real axis from UV asymptotics are shown to be erased by a flat space analog of the AdS resolution of the bulk point singularity. The residues of the poles on the negative axis are identified with operator coefficients in the IR effective action. Far along the real positive axis, the scattering is argued to grow exponentially according to the black hole area law. Exclusive amplitudes are shown to simply factorize into conformally hard and conformally soft factors. The soft factor contains all IR divergences and is given by a celestial current algebra correlator of Goldstone bosons from spontaneously broken asymptotic symmetries. The hard factor describes the scattering of hard particles together with the boost-eigenstate clouds of soft photons or gravitons required by asymptotic symmetries. These provide an IR safe S-matrix for the scattering of hard particles.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Non-Boolean quantum amplitude amplification and quantum mean estimation

This paper generalizes the quantum amplitude amplification and amplitude estimation algorithms to work with non-Boolean oracles. The action of a non-Boolean oracle $U_\varphi $ on an eigenstate $\mathinner {|{x}\rangle }$ is to apply a state-dependent phase-shift $\varphi (x)$. Unlike Boolean oracles, the eigenvalues $\exp (i\varphi (x))$ of a non-Boolean oracle are not restricted to be $\pm 1$. Two new oracular algorithms based on such non-Boolean oracles are introduced. The first is the non-Boolean amplitude amplification algorithm, which preferentially amplifies the amplitudes of the eigenstates based on the value of $\varphi (x)$. Starting from a given initial superposition state $\mathinner {|{\psi _0}\rangle }$, the basis states with lower values of $\cos (\varphi )$ are amplified at the expense of the basis states with higher values of $\cos (\varphi )$. The second algorithm is the quantum mean estimation algorithm, which uses quantum phase estimation to estimate the expectation $\mathinner {\langle {\psi _0|U_\varphi |\psi _0}\rangle }$, i.e., the expected value of $\exp (i\varphi (x))$ for a random x sampled by making a measurement on $\mathinner {|{\psi _0}\rangle }$. It is shown that the quantum mean estimation algorithm offers a quadratic speedup over the corresponding classical algorithm. Both algorithms are demonstrated using simulations for a toy example. Potential applications of the algorithms are briefly discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Efimov-like states and quantum funneling effects on synthetic hyperbolic surfaces

Engineering lattice models with tailored inter-site tunnelings and onsite energies could synthesize essentially arbitrary Riemannian surfaces with highly tunable local curvatures. Here, we point out that discrete synthetic Poincaré half-planes and Poincaré disks, which are created by lattices in flat planes, support infinitely degenerate eigenstates for any nonzero eigenenergies. Such Efimov-like states exhibit a discrete scaling symmetry and imply an unprecedented apparatus for studying quantum anomaly using hyperbolic surfaces. Furthermore, all eigenstates are exponentially localized in the hyperbolic coordinates, signifying the first example of quantum funneling effects in Hermitian systems. As such, any initial wave packet travels towards the edge of the Poincaré half-plane or its equivalent on the Poincaré disk, delivering an efficient scheme to harvest light and atoms in two dimensions. Furthermore, our findings unfold the intriguing properties of hyperbolic spaces and suggest that Efimov states may be regarded as a projection from a curved space with an extra dimension.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Precise determination of the $B^0_S - \bar {B}^0_S$ oscillation frequency

Mesons comprising a beauty quark and strange quark can oscillate between particle ($B^0_s$) and antiparticle ($\bar{B}^0_s$) flavour eigenstates, with a frequency given by the mass difference between heavy and light mass eigenstates, Δm s . Here we present a measurement of Δm s using $B^0_s$ → $D^–_sπ^+$ decays produced in proton–proton collisions collected with the LHCb detector at the Large Hadron Collider. The oscillation frequency is found to be Δm s = 17.7683 ± 0.0051 ± 0.0032 ps –1 , where the first uncertainty is statistical and the second is systematic. This measurement improves on the current Δm s precision by a factor of two. We combine this result with previous LHCb measurements to determine Δm s = 17.7656 ± 0.0057 ps –1 , which is the legacy measurement of the original LHCb detector.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Anomalous localization in a kicked quasicrystal

Abstract Quantum transport can distinguish between dynamical phases of matter. For instance, ballistic propagation characterizes the absence of disorder, whereas in many-body localized phases, particles do not propagate for exponentially long times. Additional possibilities include states of matter exhibiting anomalous transport in which particles propagate with a non-trivial exponent. Here we report the experimental observation of anomalous transport across a broad range of the phase diagram of a kicked quasicrystal. The Hamiltonian of our system has been predicted to exhibit a rich phase diagram, including not only fully localized and fully delocalized phases but also an extended region comprising a nested pattern of localized, delocalized and multifractal states, which gives rise to anomalous transport. Our cold-atom realization is enabled by new Floquet engineering techniques, which expand the accessible phase diagram by five orders of magnitude. Mapping transport properties throughout the phase diagram, we observe disorder-driven re-entrant delocalization and sub-ballistic transport, and we present a theoretical explanation of these phenomena based on eigenstate multifractality.

Physics↗

Quantum scars of bosons with correlated hopping

Abstract Recent experiments on Rydberg atom arrays have found evidence of anomalously slow thermalization and persistent density oscillations, which have been interpreted as a many-body analog of the phenomenon of quantum scars. Periodic dynamics and atypical scarred eigenstates originate from a “hard” kinetic constraint: the neighboring Rydberg atoms cannot be simultaneously excited. Here we propose a realization of quantum many-body scars in a 1D bosonic lattice model with a “soft” constraint in the form of density-assisted hopping. We discuss the relation of this model to the standard Bose-Hubbard model and possible experimental realizations using ultracold atoms. We find that this model exhibits similar phenomenology to the Rydberg atom chain, including weakly entangled eigenstates at high energy densities and the presence of a large number of exact zero energy states, with distinct algebraic structure.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Evaluation of the excitation spectra with diffusion Monte Carlo on an auxiliary bosonic ground state

We aim to improve upon the variational Monte Carlo (VMC) approach for excitations replacing the Jastrow factor by an auxiliary bosonic (AB) ground state and multiplying it by a fermionic component factor. The instantaneous change in imaginary time of an arbitrary excitation in the original interacting fermionic system is obtained by measuring observables via the ground-state distribution of walkers of an AB system that is subject to an auxiliary effective potential. The effective potential is used to (i) drive the AB system’s ground-state configuration space toward the configuration space of the excitations of the original fermionic system and (ii) subtract from a diffusion Monte Carlo (DMC) calculation contributions that can be included in conventional approximations, such as mean-field and configuration interaction (CI) methods. In this novel approach, the AB ground state is treated statistically in DMC, whereas the fermionic component of the original system is expanded in a basis. The excitation energies of the fermionic eigenstates are obtained by sampling a fermion–boson coupling term on the AB ground state. We show that this approach can take advantage of and correct for approximate eigenstates obtained via mean-field calculations or truncated interactions. We demonstrate that the AB ground-state factor incorporates the correlations missed by standard Jastrow factors, further reducing basis truncation errors. Relevant parts of the theory have been tested in soluble model systems and exhibit excellent agreement with exact analytical data and CI and VMC approaches. In particular, for limited basis set expansions and sufficient statistics, AB approaches outperform CI and VMC in terms of basis size for the same systems. The implementation of this method in current codes, despite being demanding, will be facilitated by reusing procedures already developed for calculating ground-state properties with DMC and excitations with VMC.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Interpretations of ground-state symmetry breaking and strong correlation in wavefunction and density functional theories

Significance The ground state of a quantum mechanical system is the lowest-energy eigenstate of the Hamiltonian. In isolation, it persists unchanged forever, with symmetries dictated by those of the Hamiltonian. But near-eigenstates of broken symmetry can persist for long times, even on the scale of human measurement. The appearance of broken symmetries of the electron density or spin density in a density functional calculation can reveal strong correlations among the electrons that are present in a symmetry-unbroken wavefunction. Symmetry breaking can arise when a wave-like fluctuation drops to zero frequency. The presented examples are the stretched hydrogen molecule, antiferromagnetism in solids, and the static charge-density wave in a low-density jellium, which is shown quantitatively to be a zero-frequency plasma wave.

42 ENGINEERING↗

Surrogate models for plasma displacement and current in 3D perturbed magnetohydrodynamic equilibria in tokamaks

Abstract A numerical database of over one thousand perturbed three-dimensional (3D) equilibria has been generated, constructed based on the MARS-F (Liu et al 2000 Phys. Plasmas 7 3681) computed plasma response to the externally applied 3D field sources in multiple tokamak devices. Perturbed 3D equilibria with the n = 1–4 ( n is the toroidal mode number) toroidal periodicity are computed. Surrogate models are created for the computed perturbed 3D equilibrium utilizing model order reduction (MOR) techniques. In particular, retaining the first few eigenstates from the singular value decomposition (SVD) of the data is found to produce reasonably accurate MOR-representations for the key perturbed quantities, such as the perturbed parallel plasma current density and the plasma radial displacement. SVD also helps to reveal the core versus edge plasma response to the applied 3D field. For the database covering the conventional aspect ratio devices, about 95% of data can be represented by the truncated SVD-series with inclusion of only the first five eigenstates, achieving a relative error (RE) below 20%. The MOR-data is further utilized to train neural networks (NNs) to enable fast reconstruction of perturbed 3D equilibria, based on the two-dimensional equilibrium input and the 3D source field. The best NN-training is achieved for the MOR-data obtained with a global SVD approach, where the full set of samples used for NN training and testing are stretched and form a large matrix which is then subject to SVD. The fully connected multi-layer perceptron, with one or two hidden layers, can be trained to predict the MOR-data with less than 10% RE. As a key insight, a better strategy is to train separate NNs for the plasma response fields with different toroidal mode numbers. It is also better to apply MOR and to subsequently train NNs separately for conventional and low aspect ratio devices, due to enhanced toroidal coupling of Fourier spectra in the plasma response in the latter case.

3D equilibrium↗

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↗

Filtering states with total spin on a quantum computer

Starting from a general wave-function described on a set of spins or qubits, we propose several quantum algorithms to extract the components of this state on eigenstates of total spin S 2 and its z component S z . The method plays the role of total spin projection and gives access to the amplitudes of the initial state on a total spin basis. Here, the different algorithms have various degrees of sophistication depending on the requested tasks. They can either solely project onto the subspace with good total spin or completely uplift the degeneracy in this subspace. In the former case, when the projection on the total spin j is made, we show that the number of operations for the projection can be reduced from a quadratic to a linear dependence in j. After each measurement, the state collapses to one of the spin eigenstates that could be used for postprocessing. For this reason, we call the method total quantum spin filtering.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Asymptotic errors in adiabatic evolution

The adiabatic theorem in quantum mechanics implies that if a system is in a discrete eigenstate of a Hamiltonian and the Hamiltonian evolves in time arbitrarily slowly, the system will remain in the corresponding eigenstate of the evolved Hamiltonian. Understanding corrections to the adiabatic result that arise when the evolution of the Hamiltonian is slow—but not arbitrarily slow—has become increasingly important, especially since adiabatic evolution has been proposed as a method of state preparation in quantum computing. Here, this paper identifies two regimes, an adiabatic regime in which corrections are generically small and can depend on details of the evolution throughout the path, and a hyperadiabatic regime in which the error is given by a form similar to an asymptotic expansion in the inverse of the evolution time with the coefficients depending principally on the behavior at the endpoints. However, the error in this hyperadiabatic regime is neither given by a true asymptotic series nor solely dependent on the endpoints: the coefficients combine the contributions from both endpoints, with relative phase factors that depend on the average spectral gaps along the trajectory, multiplied by the evolution time. The central result of this paper is to identify a quantity, referred to as the typical error, which is obtained by appropriately averaging the error over evolution times that are small compared to the evolution time itself. This typical error is characterized by an asymptotic series and depends solely on the endpoints of the evolution, remaining independent of the details of the intermediate evolution.

adiabatic approximation↗

Magic-angle semimetals with chiral symmetry

Here, we construct and solve a two-dimensional, chirally symmetric model of Dirac cones subjected to a quasiperiodic modulation. In real space, this is realized with a quasiperiodic hopping term. This hopping model, as we show, at the Dirac node energy has a rich phase diagram with a semimetal-to-metal phase transition at intermediate amplitude of the quasiperiodic modulation, and a transition to a phase with a diverging density of states (DOS) and subdiffusive transport when the quasiperiodic hopping is strongest. We further demonstrate that the semimetal-to-metal phase transition can be characterized by the multifractal structure of eigenstates in momentum space and can be considered as a unique “unfreezing” transition. This unfreezing transition in momentum space generates flat bands with a dramatically renormalized bandwidth in the metallic phase similar to the phenomena of the band structure of twisted bilayer graphene at the magic angle. We characterize the nature of this transition numerically as well as analytically in terms of the formation of a band of topological zero modes. For pure quasiperiodic hopping, we provide strong numerical evidence that the low-energy DOS develops a divergence and the eigenstates exhibit Chalker (quantum-critical) scaling despite the model not being random. At particular commensurate limits the model realizes higher-order topological insulating phases. We discuss how these systems can be realized in experiments on ultracold atoms and metamaterials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Extended nonergodic regime and spin subdiffusion in disordered SU(2)-symmetric Floquet systems

In this work, we explore thermalization and quantum dynamics in a one-dimensional disordered SU(2)-symmetric Floquet model, where a many-body localized phase is prohibited by the non-Abelian symmetry. Despite the absence of localization, we find an extended nonergodic regime at strong disorder where the system exhibits nonthermal behaviors. In the strong disorder regime, the level spacing statistics exhibit neither a Wigner-Dyson nor a Poisson distribution, and the spectral form factor does not show a linear-in-time growth at early times characteristic of random matrix theory. The average entanglement entropy of the Floquet eigenstates is subthermal, although violating an area-law scaling with system sizes. We further compute the expectation value of local observables and find strong deviations from the eigenstate thermalization hypothesis. The infinite-temperature spin autocorrelation function decays at long times as $t^{–β}$ with $\textit{β}$ < 0.5, indicating subdiffusive transport at strong disorders.

36 MATERIALS SCIENCE↗

Exact quantum scars in the chiral nonlinear Luttinger liquid

While the chiral linear Luttinger liquid is integrable via bosonization, its nonlinear counterpart does not admit for an analytic solution. In this work, we find a subextensive number of exact eigenstates for a large family of density-density interaction terms. These states are embedded in a continuum of strongly correlated excited states. The real-space entanglement entropy of some exact states scales logarithmically with system size while that of others has volume-law scaling. We introduce momentum-space entanglement as an unambiguous differentiator between these exact states and the remaining excited states. With regard to momentum space, the exact states behave as bona fide quantum many-body scars: they exhibit identically zero momentum-space entanglement, while typical eigenstates behave thermally. We corroborate this finding by a level statistics analysis. Here we detail the general formalism for systematically finding all interaction terms and associated exact states, and present a number of infinite exact state sequences extending to arbitrarily high energies. Unlike many previous examples of quantum many-body scars, the exact states uncovered here do not lie at equidistant energies and do not follow from a special operator algebra. Instead, they are uniquely enabled by the interplay of Fermi statistics and chirality.

36 MATERIALS SCIENCE↗