Time-evolved State Research Articles

Distribution of quantum gravity induced entanglement in many-body systems

Abstract Recently, it was shown that if two distant test masses, each in a spatially superposed quantum state, become entangled due to their mutual gravitational interaction, then this entanglement could serve as evidence of the quantum nature of gravity We extend this treatment to a many-body system in a general setup and study the entanglement properties of the time-evolved state. We exactly compute the time-dependent I-concurrence for every bipartition and obtain the necessary and sufficient condition for the creation of genuine many-body entanglement. We further show that this entanglement is of generalised GHZ type when certain conditions are met. We also evaluate the amount of multipartite entanglement in the system using a set of generalised Meyer-Wallach measures.

Measurable Krylov spaces and eigenenergy count in quantum state dynamics

In this work, we propose a quantum-mechanically measurable basis for the computation of spread complexity. Current literature focuses on computing different powers of the Hamiltonian to construct a basis for the Krylov state space and the computation of the spread complexity. We show, through a series of proofs, that time-evolved states with different evolution times can be used to construct an equivalent space to the Krylov state space used in the computation of the spread complexity. Afterwards, we introduce the effective dimension, which is upper-bounded by the number of pairwise distinct eigenvalues of the Hamiltonian. The computation of the spread complexity requires knowledge of the Hamiltonian and a classical computation of the different powers of the Hamiltonian. The computation of large powers of the Hamiltonian becomes increasingly difficult for large systems. The first part of our work addresses these issues by defining an equivalent space, where the original basis consists of quantum-mechanically measurable states. We demonstrate that a set of different time-evolved states can be used to construct a basis. We subsequently verify the results through numerical analysis, demonstrating that every time-evolved state can be reconstructed using the defined vector space. Based on this new space, we define an upper-bounded effective dimension and analyze its influence on finite-dimensional systems. We further show that the Krylov space dimension is equal to the number of pairwise distinct eigenvalues of the Hamiltonian, enabling a method to determine the number of eigenenergies the system has experimentally. Lastly, we compute the spread complexities of both basis representations and observe almost identical behavior, thus enabling the computation of spread complexities through measurements.

Relaxation dynamics in the alternating XY chain following a quantum quench

We investigate the relaxation dynamics of the fermion two-point correlation function Cmn(t)=〈ψ(t)∣cm†cn∣ψ(t)〉 in the XY chain with alternating nearest-neighbor hopping interaction after a quench. We find that the deviation δ C mn (t) = C mn (t) − C mn (∞) decays with time following the power law behavior t −μ , where the exponent μ depends on whether the quench is to the commensurate phase (μ = 1) or incommensurate phase ( μ=12 ). This decay of δ C mn (t) arises from the transient behavior of the double-excited quasiparticle occupations and the transitions between different excitation spectra. Furthermore, we find that the steady value C mn (∞) only involves the average fermion occupation numbers (i.e. the average excited single particle) over the time-evolved state, which are different from the ground state expectation values. We also observe nonanalytic singularities in the steady value C mn (∞) for the quench to the critical points of the quantum phase transitions (QPTs), suggesting its potential use as a signature of QPTs.

Canonical Typicality for Other Ensembles than Micro-canonical

AbstractWe generalize Lévy’s lemma, a concentration-of-measure result for the uniform probability distribution on high-dimensional spheres, to a much more general class of measures, so-called GAP measures. For any given density matrix $$\rho $$ ρ on a separable Hilbert space $${\mathcal {H}}$$ H , $${\textrm{GAP}}(\rho )$$ GAP ( ρ ) is the most spread-out probability measure on the unit sphere of $${\mathcal {H}}$$ H that has density matrix $$\rho $$ ρ and thus forms the natural generalization of the uniform distribution. We prove concentration-of-measure whenever the largest eigenvalue $$\Vert \rho \Vert $$ ‖ ρ ‖ of $$\rho $$ ρ is small. We use this fact to generalize and improve well-known and important typicality results of quantum statistical mechanics to GAP measures, namely canonical typicality and dynamical typicality. Canonical typicality is the statement that for “most” pure states $$\psi $$ ψ of a given ensemble, the reduced density matrix of a sufficiently small subsystem is very close to a $$\psi $$ ψ -independent matrix. Dynamical typicality is the statement that for any observable and any unitary time evolution, for “most” pure states $$\psi $$ ψ from a given ensemble the (coarse-grained) Born distribution of that observable in the time-evolved state $$\psi _t$$ ψ t is very close to a $$\psi $$ ψ -independent distribution. So far, canonical typicality and dynamical typicality were known for the uniform distribution on finite-dimensional spheres, corresponding to the micro-canonical ensemble, and for rather special mean-value ensembles. Our result shows that these typicality results hold also for $${\textrm{GAP}}(\rho )$$ GAP ( ρ ) , provided the density matrix $$\rho $$ ρ has small eigenvalues. Since certain GAP measures are quantum analogs of the canonical ensemble of classical mechanics, our results can also be regarded as a version of equivalence of ensembles.

Exact time-evolving scattering states in open quantum-dot systems with an interaction: discovery of time-evolving resonant states

We study exact time-evolving many-electron states of an open double quantum-dot system with an interdot Coulomb interaction. A systematic construction of the time-evolving states for arbitrary initial conditions is proposed. For any initial states of one- and two-electron plane waves on the electrical leads, we obtain exact solutions of the time-evolving scattering states, which converge to known stationary scattering eigenstates in the long-time limit. For any initial states of localized electrons on the quantum dots, we find exact time-evolving states of a new type, which we refer to as time-evolving resonant states. In contrast to stationary resonant states, whose wave functions spatially diverge and not normalizable, the time-evolving resonant states are normalizable since their wave functions are restricted to a finite space interval due to causality. The exact time-evolving resonant states enable us to calculate the time-dependence of the survival probability of electrons on the quantum dots for the system with the linearized dispersions. It decays exponentially in time on one side of an exponential point of resonance energies while, on the other side, it oscillates during the decay as a result of the interference of the two resonance energies.

Notes on time entanglement and pseudo-entropy

Following arXiv:2210.12963 [hep-th], we investigate aspects of the time evolution operator regarded as a density operator and associated entanglement-like structures in various quantum systems. These involve timelike separations and generically lead to complex-valued entropy, although there are interesting real subfamilies. There are many parallels and close relations with reduced transition matrices and pseudo-entropy, which we discuss and clarify. For instance, a related quantity involves the time evolution operator along with a projection onto some initial state, which amounts to analysing pseudo-entropy for the initial state and its time-evolved final state.

Adiabatic theorem for classical stochastic processes

We apply adiabatic theorems developed for quantum mechanics to stochastic annealing processes described by the classical master equation with a time-dependent generator. When the instantaneous stationary state is unique and the minimum decay rate g is nonzero, the time-evolved state is basically relaxed to the instantaneous stationary state. By formulating an asymptotic expansion rigorously, we derive conditions for the annealing time τ that the state is close to the instantaneous stationary state. Depending on the time dependence of the generator, typical conditions are written as τ>const.×g−α with 1⩽α⩽2 . We also find that a rigorous treatment gives the scaling τ>const.×g−2|ln⁡g| .

Converting Long-Range Entanglement into Mixture: Tensor-Network Approach to Local Equilibration.

In the out-of-equilibrium evolution induced by a quench, fast degrees of freedom generate long-range entanglement that is hard to encode with standard tensor networks. However, local observables only sense such long-range correlations through their contribution to the reduced local state as a mixture. We present a tensor network method that identifies such long-range entanglement and efficiently transforms it into mixture, much easier to represent. In this way, we obtain an effective description of the time-evolved state as a density matrix that captures the long-time behavior of local operators with finite computational resources.

Adiabatic Evolution of Low-Temperature Many-Body Systems

We consider finite-range, many-body fermionic lattice models and we study the evolution of their thermal equilibrium state after introducing a weak and slowly varying time-dependent perturbation. Under suitable assumptions on the external driving, we derive a representation for the average of the evolution of local observables via a convergent expansion in the perturbation, for small enough temperatures. Convergence holds for a range of parameters that is uniform in the size of the system. Under a spectral gap assumption on the unperturbed Hamiltonian, convergence is also uniform in temperature. As an application, our expansion allows us to prove closeness of the time-evolved state to the instantaneous Gibbs state of the perturbed system, in the sense of expectation of local observables, at zero and at small temperatures. As a corollary, we also establish the validity of linear response. Our strategy is based on a rigorous version of the Wick rotation, which allows us to represent the Duhamel expansion for the real-time dynamics in terms of Euclidean correlation functions, for which precise decay estimates are proved using fermionic cluster expansion.

Complete Hilbert-Space Ergodicity in Quantum Dynamics of Generalized Fibonacci Drives.

Ergodicity of quantum dynamics is often defined through statistical properties of energy eigenstates, as exemplified by Berry's conjecture in single-particle quantum chaos and the eigenstate thermalization hypothesis in many-body settings. In this work, we investigate whether quantum systems can exhibit a stronger form of ergodicity, wherein any time-evolved state uniformly visits the entire Hilbert space over time. We call such a phenomenon complete Hilbert-space ergodicity (CHSE), which is more akin to the intuitive notion of ergodicity as an inherently dynamical concept. CHSE cannot hold for time-independent or even time-periodic Hamiltonian dynamics, owing to the existence of (quasi)energy eigenstates which precludes exploration of the full Hilbert space. However, we find that there exists a family of aperiodic, yet deterministic drives with minimal symbolic complexity-generated by the Fibonacci word and its generalizations-for which CHSE can be proven to occur. Our results provide a basis for understanding thermalization in general time-dependent quantum systems.

Tensor network variational optimizations for real-time dynamics: Application to the time-evolution of spin liquids

Within the Projected Entangled Pair State (PEPS) tensor network formalism, a simple update (SU) method has been used to investigate the time evolution of a two-dimensional U(1) critical spin-1/2 spin liquid under Hamiltonian quench [Phys. Rev. B 106, 195132 (2022)]. Here we introduce two different variational frameworks to describe the time dynamics of SU(2)-symmetric translationally-invariant PEPS, aiming to improve the accuracy. In one approach, after using a Trotter-Suzuki decomposition of the time evolution operator in term of two-site elementary gates, one considers a single bond embedded in an environment approximated by a Corner Transfer Matrix Renormalization Group (CTMRG). A variational update of the two tensors on the bond is performed under the application of the elementary gate and then, after symmetrization of the site tensors, the environment is updated. In the second approach, a cluster optimization is performed on a finite (periodic) cluster, maximizing the overlap of the exact time-evolved state with a symmetric finite-size PEPS ansatz. Observables are then computed on the infinite lattice contracting the infinite-PEPS (iPEPS) by CTMRG. We show that the variational schemes outperform the SU method and remain accurate over a significant time interval before hitting the entanglement barrier. Studying the spectrum of the transfer matrix, we find that the asymptotic correlations are very well preserved under time evolution, including the critical nature of the singlet correlations, as expected from the Lieb-Robinson (LR) bound theorem. Consistently, the system (asymptotic) boundary is found to be described by the same Conformal Field Theory of central charge c = 1c=1 during time evolution. We also compute the time-evolution of the short distance spin-spin correlations and estimate the LR velocity.

Time evolution of spread complexity in quenched Lipkin–Meshkov–Glick model

We use the spread complexity (SC) of a time-evolved state after a sudden quantum quench in the Lipkin–Meshkov–Glick (LMG) model prepared in the ground state as a probe of the quantum phase transition when the system is quenched toward the critical point. By studying the growth of the effective number of elements of the Krylov basis that contributes to the SC more than a preassigned cutoff, we show how the two phases of the LMG model can be distinguished. We also explore the time evolution of spread entropy after both non-critical and critical quenches. We show that the sum contributing to the spread entropy converges slowly in the symmetric phase of the LMG model compared to that in the broken phase, and for a critical quench, the spread entropy diverges logarithmically at late times.

Non-Markovian quantum state diffusion for spin environments

We introduce an exact open system method to describe the dynamics of quantum systems that are strongly coupled to specific types of environments comprising of spins, such as central spin systems. Our theory is similar to the established non-Markovian quantum state diffusion theory, but for a spin bath instead of a Gaussian bath. The method allows us to represent the time-evolved reduced state of the system as an ensemble average of stochastically evolving pure states. We present a comprehensive theory for arbitrary linear spin environments at both zero and finite temperatures. Furthermore, we introduce a hierarchical expansion method that enables the numerical computation of the time evolution of the stochastic pure states, facilitating a numerical solution of the open system problem in relevant strong coupling regimes.

Generalized speed limits for classical stochastic systems and their applications to relaxation, annealing, and pumping processes

We extend the speed limit of a distance between two states evolving by different generators for quantum systems [K. Suzuki and K. Takahashi, Phys. Rev. Res. 2, 032016(R) (2020)] to the classical stochastic processes described by the master equation. We demonstrate that the trace distance between arbitrary evolving states is bounded from above by using a geometrical metric. The geometrical bound reduces to the Fisher information metric for the distance between the time-evolved state and the initial state. We compare the bound in relaxation and annealing processes with a different type of bound known for nonequilibrium thermodynamical systems. For dynamical processes such as annealing and pumping processes, the distance between the time-evolved state and the instantaneous stationary state becomes a proper choice and the bound is represented by the Fisher information metric of the stationary state. The metric is related to the counterdiabatic term defined from the time dependence of the stationary state.

Optimal Hamiltonian simulation for time-periodic systems

The implementation of time-evolution operators U(t), called Hamiltonian simulation, is one of the most promising usage of quantum computers. For time-independent Hamiltonians, qubitization has recently established efficient realization of time-evolution U(t)=e−iHt, with achieving the optimal computational resource both in time t and an allowable error ε. In contrast, those for time-dependent systems require larger cost due to the difficulty of handling time-dependency. In this paper, we establish optimal/nearly-optimal Hamiltonian simulation for generic time-dependent systems with time-periodicity, known as Floquet systems. By using a so-called Floquet-Hilbert space equipped with auxiliary states labeling Fourier indices, we develop a way to certainly obtain the target time-evolved state without relying on either time-ordered product or Dyson-series expansion. Consequently, the query complexity, which measures the cost for implementing the time-evolution, has optimal and nearly-optimal dependency respectively in time t and inverse error ε, and becomes sufficiently close to that of qubitization. Thus, our protocol tells us that, among generic time-dependent systems, time-periodic systems provides a class accessible as efficiently as time-independent systems despite the existence of time-dependency. As we also provide applications to simulation of nonequilibrium phenomena and adiabatic state preparation, our results will shed light on nonequilibrium phenomena in condensed matter physics and quantum chemistry, and quantum tasks yielding time-dependency in quantum computation.

Effects of solvent polarity on the novel excited-state intramolecular thiol proton transfer and photophysical property compared with the oxygen proton transfer

Recently, the dual-fluorescent phenomena of excited state intramolecular thiol proton transfer (ESIPT) for 3-thiolflavone derivative (3NTF) were reported by Chou and coworkers for the first time [J. Am. Chem. Soc. 143 (2021) 12715–12724], which opened a new chapter in the field of ESIPT. Based on density functional theory (DFT) and time-dependent density functional theory (TDDFT), the proton transfer processes of 3NTF in toluene, dichloromethane and acetonitrile were studied. By optimizing the structure of the ground (S0) state and first excited (S1) state of 3NTF in different solvents, the hydrogen-bond parameters and proton-transfer potential energy curves were calculated. It was shown that although photo-excitation enhanced the intramolecular hydrogen bonding strength and thus promoted the occurrence of ESIPT, the solvent polarities inhibited the enhancement of the hydrogen bond of S1 state, which was not conducive to ESIPT. The electron spectra analyses were consistent with experimental data, which confirmed the rationality of molecular configurations. The time-evolved excited state dynamics simulation was performed based on the optimized structure of 3NTF, indicating that the ESIPT was an ultrafast photochemical reaction less than 180 fs. Moreover, we compared the potential energy surfaces of ESIPT, electronic structures based on natural transition orbitals (NTOs) method and electron-hole isosurfaces for the 3NTF and the traditional flavone molecule (3NHF), concluded that the unusually large Stokes shift fluorescence of 3NTF was mainly caused by the coupling of ESIPT and twisting intramolecular charge transfer (TICT), and the 3NTF isomer had the more nπ* character in the electron transition process. The nπ* ICT significantly increased with the decrease of solvent polarities, affecting the molecular photophysical properties, this made it more widely used in biomedical, photochemical, materials science and other fields.

Entanglement between uncoupled modes with time-dependent complex frequencies

In this work we present the general unified description for the unitary time evolution generated by time-dependent non-Hermitian Hamiltonians embedding the bosonic representations of $\mathfrak{su}(1,1)$ and $\mathfrak{su}(2)$ Lie algebras. We take into account a time-dependent Hermitian Dyson maps written in terms of the elements of those algebras with the relation between non-Hermitian and its Hermitian counterpart being independent of the algebra realization. As a direct consequence, we verify that a time-evolved state of uncoupled modes modulated by a time-dependent complex frequency may exhibits a nonzero entanglement even when the cross operators, typical of the interaction between modes, are absent. This is due the nonlocal nature of the nontrivial dynamical Hilbert space metric encoded in the time-dependent parameters of the general Hermitian Dyson map, which depend on the imaginary part of the complex frequency. We illustrate our approach by setting the $\mathcal{PT}$-symmetric case where the imaginary part of frequency is linear on time for the two-mode bosonic realization of Lie algebras.

Transient phenomena in the axion assisted Schwinger effect

Particle production induced by a time-dependent background is well understood as the projection of the time-evolved initial state onto a set of final states. While the asymptotic initial and final states are well defined in the usual way, the definition of particles and antiparticles at intermediate times in the presence of external fields is ambiguous. These external fields moreover induce divergences which require regularization. In this paper we clarify some subtleties in the computation of transient effects in physical quantities for fermions in a homogeneous axion background, including Schwinger production in background electromagnetic fields. The presence of the axion requires particular care as well as knowledge of the UV theory when regulating the theory and computing the vacuum contribution to the fermion energy.

Performance evaluation of invariant-based inverse engineering by quantum speed limits

Quantum speed limits for two time-evolved states are introduced and applied to overlap between true dynamics and approximate dynamics. In particular, we point out that the present idea is suitable for invariant-based inverse engineering, i.e., the worst case performance of invariant-based inverse engineering can be evaluated by using a designed time-evolved state. As demonstrations, we apply the present method to stimulated Raman adiabatic passage and quantum annealing, and then we find that the present idea brings valuable insight to control scheduling.

Energy-filtered random-phase states as microcanonical thermal pure quantum states

We propose a method to calculate finite-temperature properties of a quantum many-body system for a microcanonical ensemble by introducing a pure quantum state named here an energy-filtered random-phase state, which is also a potentially promising application of near-term quantum computers. In our formalism, a microcanonical ensemble is specified by two parameters, i.e., the energy of the system and its associated energy window. Accordingly, the density of states is expressed as a sum of Gaussians centered at the target energy with its spread corresponding to the width of the energy window. We then show that the thermodynamic quantities such as entropy and temperature are calculated by evaluating the trace of the time-evolution operator and the trace of the time-evolution operator multiplied by the Hamiltonian of the system. We also describe how these traces can be evaluated using random diagonal-unitary circuits appropriate to quantum computation. The pure quantum state representing our microcanonical ensemble is related to a state of the form introduced by Wall and Neuhauser for the filter diagonalization method [M. R. Wall and D. Neuhauser, J. Chem. Phys. 102, 8011 (1995)], and therefore we refer to it as an energy-filtered random-phase state. The energy-filtered random-phase state is essentially a Fourier transform of a time-evolved state whose initial state is prepared as a random-phase state, and the cut-off time in the time-integral for the Fourier transform sets the inverse of the width of the energy window. The proposed method is demonstrated numerically by calculating thermodynamic quantities for the one-dimensional spin-1/2 Heisenberg model on small clusters up to 28 qubits, showing that the method is most efficient for the target energy around which the dense distribution of energy eigenstates is found.