We study the nonequilibrium dynamics of an isolated bipartite quantum system, the sunburst quantum Ising model, under interaction quench. The prequench limit of this model is two noninteracting integrable systems, namely a transverse Ising chain and finite number of isolated qubits. As a function of interaction strength, the spectral fluctuation property goes from Poisson to Wigner-Dyson statistics. We chose entanglement entropy as a probe to study the approach to thermalization or lack of it in postquench dynamics. In the near-integrable limit, as expected, the linear entropy displays oscillatory behavior, while in the chaotic limit it saturates. Along with the chaotic nature of the time evolution generator, we show the importance of the role played by the coherence of the initial state in deciding the nature of thermalization. We further show that these findings are general by replacing the Ising ring with a disordered XXZ model with disorder strength putting it in the many-body localized phase.
We study the nonequilibrium dynamics of an isolated bipartite quantum system, the sunburst quantum Ising model, under interaction quench. The prequench limit of this model is two noninteracting integrable systems, namely a transverse Ising chain and finite number of isolated qubits. As a function of interaction strength, the spectral fluctuation property goes from Poisson to Wigner-Dyson statistics. We chose entanglement entropy as a probe to study the approach to thermalization or lack of it in postquench dynamics. In the near-integrable limit, as expected, the linear entropy displays oscillatory behavior, while in the chaotic limit it saturates. Along with the chaotic nature of the time evolution generator, we show the importance of the role played by the coherence of the initial state in deciding the nature of thermalization. We further show that these findings are general by replacing the Ising ring with a disordered XXZ model with disorder strength putting it in the many-body localized phase.
Abstract We study the post-quench unitary dynamics of a quantum sunburst spin model, composed of a transverse-field quantum Ising ring which is suddenly coupled to a set of independent external qubits along the longitudinal direction, in a way to respect a residual translation invariance and the Ising Z 2 symmetry. Starting from the different equilibrium quantum phases of the system, we characterize the decoherence and the energy storage in the external qubits, which may be interpreted as a probing apparatus for the inner Ising ring. Our results show that, in proximity of the quantum transitions of the Ising ring, either first-order or continuous, it is possible to put forward dynamic finite-size scaling frameworks which unveil peculiar scaling regimes, depending on the way in which the large-size limit is taken: either by fixing the number n of probing qubits, or their interspace distance b . In any case, the dependence of the various observables on n can be reabsorbed into a redefinition of the quench parameter by a n prefactor. We also address the role of a nearest-neighbor coupling between the external qubits.
Quantum quench is a non-equilibrium process where the Hamiltonian is suddenly changed during the quantum evolution. The change can be made by spatially local perturbations (local quench) or globally switching to a completely different Hamiltonian (global quench). This thesis investigates the post-quench non-equilibrium dynamics with an emphasis on the time dependence of the quantum entanglement. We inspect the scaling of entanglement entropy (EE) to learn how correlation and entanglement built up in a quench. We begin with two local quench examples. In Chap. 2, we apply a local operator to the groundstate of the quantum Lifshitz model and monitor the change of the EE. We find that the entanglement grows according to the dynamical exponent z = 2 and then saturates to the scaling dimension of the perturbing operator -- a value representing its strength. In Chap. 3, we study the evolution after connecting two different one-dimensional critical chains at their ends. The Loschmidt echo which measures the similarity between the evolved state and the initial one decays with a power law, whose exponent is the scaling dimension of the defect (junction). Among other conclusions, we see that the local quench dynamics contain universal information of the (critical) theory. In the global quench scenario, the change of the Hamiltonian affects all parts of the system. In this thesis, we focus on the global chaotic quench driven by generic non-integrable Hamiltonians. In Chap. 4, we propose to use the operator entanglement entropy of the unitary operator as a probe. Its fast linear entanglement production is sharply contrasted to the slow logarithmic spreading of the many-body localized system. The entanglement saturation suggests that the evolution operator in the long time can be modeled by a random unitary matrix. In Chap. 5, we construct a random tensor network which consists of random unitary matrices connected locally to model chaotic evolution with local interactions. We find that the entanglement dynamics is mapped to the statistical mechanics of interacting random walks. This appealing emergent picture allows us to understand the universal linear growth as well as the fluctuations of entanglement in a chaotic quench.
We study the ground-state properties of a quantum sunburst model, composed of a quantum Ising spin ring in a transverse field, symmetrically coupled to a set of ancillary isolated qubits, to maintain a residual translation invariance and also a Z_{2} symmetry. The large-size limit is taken in two different ways: either by keeping the distance between any two neighboring ancillary qubits fixed or by fixing their number while increasing the ring size. Substantially different regimes emerge, depending on the various Hamiltonian parameters: For small energy scale δ of the ancillary subsystem and small ring-qubit interaction κ, we observe rapid and nonanalytic changes in proximity to the quantum transitions of the Ising ring, both first order and continuous, which can be carefully controlled by exploiting renormalization-group and finite-size scaling frameworks. Smoother behaviors are instead observed when keeping δ>0 fixed and in the Ising disordered phase. The effect of an increasing number n of ancillary spins turns out to scale proportionally to sqrt[n] for sufficiently large values of n.
Investigating the scrambling of quantum information (also known as quantum chaos) in many-body quantum systems provides valuable insights into the postquench evolution of a system. While various methods are available to understand the chaotic dynamics of these quantum systems, scrambling of mutual information following a quench has recently also emerged as a diagnostic of chaotic dynamics. In this work, we consider a chain of spin-1/2 particles of a finite length L evolving with the mixed-field Ising Hamiltonian and impose open boundary conditions. We simulate the time evolution of entanglement entropy and mutual information following a quench from the N\'eel state in this system using tensor networks and find that the entanglement entropy for nonintegrable systems saturates to a constant value at late times, while it continues to oscillate for integrable systems. We also find that mutual information peaks as a function of distance between intervals decay faster for nonintegrable systems compared to integrable systems, in agreement with the results in literature for XXZ chains. We compare the oscillations in entanglement entropy evolution obtained from simulations in the integrable case with analytic results from the quasiparticle picture and find agreement.
In one-dimensional low-density Jaynes-Cummings Hubbard (JCH) models [Phys. Rev. E 106, 064107 (2022)2470-004510.1103/PhysRevE.106.064107], we proved that the eigenstate thermalization hypothesis (ETH) is valid when the tunneling strength and coupling strength are of the same order. Surprisingly, at the weak tunneling limit, we observed that the entanglement entropy and scaling law of kinetic energy operators also exhibit obvious quantum chaotic properties, this is an unexpected result. To substantiate these findings, we further discuss their nonequilibrium dynamics in this paper. Our analysis reveals that when the model is a weak tunneling limit after the quench and the initial state is an equilibrium state of chaos, the system reaches an equilibrium state. This observation supports the conclusion that the low-density JCH model at the weak tunneling limit is nonintegrable, corroborating our previous results [Phys. Rev. E 106, 064107 (2022)2470-004510.1103/PhysRevE.106.064107]. Additionally, by discussing the validity of the fluctuation-dissipation theorem (FDT) and the evolution behavior of entanglement entropy and fidelity, we numerically demonstrate the differences between the one-dimensional low-density JCH model and general nonintegrable systems. Specifically, in the low-density JCH model, when the Hamiltonian after the quench is integrable, the validity of FDT depends on the thermal behavior of the initial Hamiltonian, and a metastable state is observed during the evolution of entanglement entropy. Our research presents an an intriguing and unique nonintegrable model, enriching the current understanding of nonintegrable systems.
The transverse field Ising model of N fully connected spins provides a mathematical accessible, yet non-trivial, system to investigate the non-thermal dynamics after a sudden quantum quench. The focus in this thesis is three-fold. First, in a 1/N expansion the dynamics of the mean magnetization within its variance is obtained analytically in the large N limit. The variance constitutes a leading order correction to the mean field limit at infinite N. While mean field theories are one of the most accessible approximations to deal with the quantum complexity, its validity in time out of equilibrium has to be checked on a case by case basis. By studying the dynamics of the variance, we see that the mean field approximation can break down at surprisingly early times scaling with the square root of N. The underlying mechanism in the fully connected Ising model is identified as a dephasing effect between nearby effective orbits. The dynamics of the magnetization variance discriminates four qualitatively different regions in the dynamical phase diagram that are indistinguishable by just looking at the mean magnetization. Second, bipartite entanglement between disjoint groups of the spins in the fully connected Ising model is studied. In the large N limit the entanglement Hamiltonian is shown to be a harmonic oscillator. A quantitative relation between the angular frequency of the oscillator and a squeezing measure of the spin state is derived. Explicit expressions for the Renyi entanglement entropies are given. In contrast to the variance, which vanishes in the thermodynamic limit, the entanglement entropy saturate to an N-independent positive value. The dynamics of the entanglement Hamiltonian and the entanglement entropies is discussed. Linear growth as a function of time is linked to the unstable fixed point and the homoclinic orbit of a pitchfork bifurcation, while logarithmic increase is explained due to a dephasing mechanism of nearby effective orbits. A generalization to other models that permit a semiclassical effective description (on a possibly higher dimensional phase space) is given. As a byproduct the connection between entangled Gaussian states and the symplectic eigenvalues of its reduced covariance is derived. An upper bound on the von Neumann entanglement entropy of any bipartite pure state is given in terms of the symplectic properties of its reduced covariance. Third, the existence of dynamical phase transitions (DPTs) in the fully connected Ising model is argued by two different approaches. DPTs are defined as non-analyticities in the time-dependence of the return probability rate function in the thermodynamic limit. Like entanglement entropy, the behavior of the return rate function is an observable-independent quantity used to characterize non-equilibrium states. A numerically based connection between the DPT and the analytic property of the infinite time averaged order parameter as a function of quench strength is discussed. The concept of Feynman history states is used to devise the notion of history entanglement. An information theoretic interpretation of the history entanglement entropy as the precision of a quantum clock follows from Holevo’s bound. The dependence of this entanglement entropy on the observation time constitutes an operator-independent characterization of the quantum dynamics in closed systems. For observation times comparable to the Heisenberg time, the history entanglement becomes sensitive to the statistics of the energy gaps, and thus distinguishes integrable and non-integrable dynamics. This is demonstrated for two examples, a random matrix ensemble, and a one-dimensional hardcore boson model with integrability breaking next nearest neighbor interaction.