Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:1207.0862 · Journal of Physics Condensed Matter

Time evolution within a comoving window: scaling of signal fronts and magnetization plateaus after a local quench in quantum spin chains

Valentin Zauner-Stauber, Martin Ganahl, Hans Gerd Evertz, Tomotoshi Nishino

Published Oct 7, 2015• 65 citations
arXivHTMLPDF

Abstract

We present a modification of Matrix Product State time evolution to simulate the propagation of signal fronts on infinite one-dimensional systems. We restrict the calculation to a window moving along with a signal, which by the Lieb-Robinson bound is contained within a light cone. Signal fronts can be studied unperturbed and with high precision for much longer times than on finite systems. Entanglement inside the window is naturally small, greatly lowering computational effort. We investigate the time evolution of the transverse field Ising (TFI) model and of the S = 1/2 XXZ antiferromagnet in their symmetry broken phases after several different local quantum quenches. In both models, we observe distinct magnetisation plateaus at the signal front for very large times, resembling those previously observed for the particle density of tight binding (TB) fermions. We show that the normalised difference to the magnetisation of the ground state exhibits similar scaling behaviour as the density of TB fermions. In the XXZ model there is an additional internal structure of the signal front due to pairing, and wider plateaus with tight binding scaling exponents for the normalised excess magnetisation. We also observe parameter dependent interaction effects between individual plateaus, resulting in a slight spatial compression of the plateau widths. In the TFI model, we additionally find that for an initial Jordan-Wigner domain wall state, the complete time evolution of the normalised excess longitudinal magnetisation agrees exactly with the particle density of TB fermions.

DOI
10.1088/0953-8984/27/42/425602

Related papers

Inhomogeneous quenches in the transverse field Ising chain: scaling and front dynamics

Márton Kormos

We investigate the non-equilibrium dynamics of the transverse field quantum Ising chain evolving from an inhomogeneous initial state given by joining two macroscopically different semi-infinite chains. We obtain integral expressions for all two-point correlation functions of the Jordan--Wigner Majorana fermions at any time and for any value of the transverse field. Using this result, we analytically compute the profiles of various physical observables in the space-time scaling limit and show that they can be obtained from a hydrodynamic picture based on ballistically propagating quasiparticles. Going beyond the hydrodynamic limit, we analyze the approach to the non-equilibrium steady state and find that the leading late time corrections display a lattice effect. We also study the fine structure of the propagating fronts which are found to be described by the Airy kernel and its derivatives. Near the front we observe the phenomenon of energy back-flow where the energy locally flows from the colder to the hotter region.

arXiv:1704.03744•SciPost Physics•Sep 8, 2017•91 citations

Light-cone and local front dynamics of a single-particle extended quantum walk

Hemlata Bhandari, P. Durganandini

We study the light-cone and front dynamics of a single-particle continuous time extended quantum walk on a one-dimensional lattice with finite range hopping. We show that, in general, for an initially localized state, propagating wave fronts can be characterized as ordinary or extremal fronts with the latter exhibiting an anomalous subdiffusive scaling behavior in the front region. We investigate the dynamical global and local scaling properties of the cumulative probability distribution function for the extended walk with nearest- and next-nearest-neighbor hopping using analytical and numerical methods. The global scaling shows the existence of a ``causal light-cone'' corresponding to excitations traveling with a velocity smaller than a maximal ``light velocity.'' Maximal fronts moving with fixed light velocity bound the causal cone. The front regions spread with time subdiffusively exhibiting a local Airy scaling which leads to an internal staircase structure. At a certain critical next-nearest-neighbor hopping strength, there is a transition from a phase with one `causal cone' to a phase with two nested `causal cones' and the existence of an internal staircase structure in the corresponding cumulative distribution profiles. We also connect the study to that in spin chain systems and indicate that a single-particle quantum walk on the one-dimensional lattice already captures the many-body physics of a spin chain system. In particular, we suggest that the time evolution of a single-particle quantum walk on the one-dimensional lattice with an initially localized state is equivalent to the time evolution of a domain-wall initial state in a corresponding spin chain system.

arXiv:1903.07838•Physical Review A•Mar 11, 2019•5 citations

The Lieb-Robinson correlation function for long disordered transverse-field Ising chains

Brendan J. Mahoney, Craig S. Lent

The Lieb-Robinson correlation function can be used to characterize the propagation of quantum information in Ising chains. Considerable work has been done to establish bounds on this correlation function in various circumstances. To actually calculate the value of the correlation function directly typically requires a state space which grows exponentially with system size, and so is intractable for all but relatively small systems. We employ a recently-developed method that enables direct calculation of the Lieb--Robinson correlation function and which scales linearly with system size. This enables the computation for systems with more than a thousand qubits, revealing the propagation of quantum information down the chain. We extend this technique to the problem of Ising chains with randomly disordered coupling strengths. Increasing disorder causes localization of the quantum correlations and halts propagation of quantum information down the chain. We see light-cone isocontours therefore becoming vertical. The decay of correlations down the chain is characterized by a correlation length which in the paramagnetic phase increases quadratically with mean coupling strength and inverse quadratically with disorder level. The correlation length peaks at the phase transition. We map our general expression for the correlation function in the Ising system onto the probability density in the single particle tight-binding Su-Schrieffer-Heeger model for a specific initial condition.

quant-pharXiv:2512.22395•Physical Review Research•Dec 30, 2025•0 citations

Lieb-Robinson correlation function for the quantum transverse-field Ising model

Brendan J. Mahoney, Craig S. Lent

The Lieb-Robinson correlation function is the norm of a commutator between local operators acting on separate subsystems at different times. This provides a useful state-independent measure for characterizing the specifically quantum interaction between spatially separated qubits. The finite propagation velocity for this correlator defines a “light cone” of quantum influence. We calculate the Lieb-Robinson correlation function for one-dimensional qubit arrays described by the transverse-field Ising model. Direct calculations of this correlation function have been limited by the exponential increase in the size of the state space with the number of qubits. We introduce a technique that avoids this barrier by transforming the calculation to a sum over Pauli walks which results in linear scaling with system size. We can then explore propagation in arrays of hundreds of qubits and observe the effects of the quantum phase transition in the system. We observe the emergence of two distinct velocities of propagation: a correlation front velocity, which is affected by the phase transition, and the Lieb-Robinson velocity which is not. The correlation front velocity is equal to the maximum group velocity of single quasiparticle excitations. The Lieb-Robinson velocity describes the extreme leading edge of correlations when the value of the correlation function itself is still very small. For the semi-infinite chain of qubits at the quantum critical point, we derive an analytical result for the correlation function. Published by the American Physical Society 2024

arXiv:2402.11080•Physical Review Research•Jun 17, 2024•3 citations

Dynamics of the entanglement spectrum in spin chains

Giacomo Torlai, Luca Tagliacozzo, Gabriele De Chiara

We study the dynamics of the entanglement spectrum, that is the time evolution of the eigenvalues of the reduced density matrices after a bipartition of a one-dimensional spin chain. Starting from the ground state of an initial Hamiltonian, the state of the system is evolved in time with a new Hamiltonian. We consider both instantaneous and quasi adiabatic quenches of the system Hamiltonian across a quantum phase transition. We analyse the Ising model that can be exactly solved and the XXZ for which we employ the time-dependent density matrix renormalisation group algorithm. Our results show once more a connection between the Schmidt gap, i.e. the difference of the two largest eigenvalues of the reduced density matrix and order parameters, in this case the spontaneous magnetisation.

arXiv:1311.5509•Journal of Statistical Mechanics Theory and Experiment•Jun 9, 2014•55 citations

Scaling and front dynamics in Ising quantum chains

Thierry Platini, Dragi Karevski

We study the relaxation behaviour of the quantum Ising chain, focusing our attention onto the non-equilibrium dynamics of the transverse magnetization. The initial states, from which the magnetization relaxes, are product states such as those generated by setting in contact several systems, each initially equilibrated at a given temperature. Due to the free fermionic structure of the chain, the dynamics of the transverse magnetization is easily expressed in a compact form. In the completely factorized initial state, corresponding to a situation where all the spins are thermalized independently, we obtain in the scaling limit the Green function associated to the transverse magnetization. The relaxation behaviour is also considered in the system-bath case, where part of the chain called the system is thermalized at a temperature Ts and the remaining part is at a temperature Tb. The magnetization profiles show a scaling behaviour. Moreover, in the extreme case Tb=∞ and Ts=0, it is shown that the magnetization relaxes in quantized steps in the strong transverse field region.

•Jan 1, 2005•34 citations
We present a modification of Matrix Product State time evolution to simulate the propagation of signal fronts on infinite one-dimensional systems. We restrict the calculation to a window moving along with a signal, which by the Lieb-Robinson bound is contained within a light cone. Signal fronts can be studied unperturbed and with high precision for much longer times than on finite systems. Entanglement inside the window is naturally small, greatly lowering computational effort. We investigate the time evolution of the transverse field Ising (TFI) model and of the S = 1/2 XXZ antiferromagnet in their symmetry broken phases after several different local quantum quenches. In both models, we observe distinct magnetisation plateaus at the signal front for very large times, resembling those previously observed for the particle density of tight binding (TB) fermions. We show that the normalised difference to the magnetisation of the ground state exhibits similar scaling behaviour as the density of TB fermions. In the XXZ model there is an additional internal structure of the signal front due to pairing, and wider plateaus with tight binding scaling exponents for the normalised excess magnetisation. We also observe parameter dependent interaction effects between individual plateaus, resulting in a slight spatial compression of the plateau widths. In the TFI model, we additionally find that for an initial Jordan-Wigner domain wall state, the complete time evolution of the normalised excess longitudinal magnetisation agrees exactly with the particle density of TB fermions.
We investigate the non-equilibrium dynamics of the transverse field quantum Ising chain evolving from an inhomogeneous initial state given by joining two macroscopically different semi-infinite chains. We obtain integral expressions for all two-point correlation functions of the Jordan--Wigner Majorana fermions at any time and for any value of the transverse field. Using this result, we analytically compute the profiles of various physical observables in the space-time scaling limit and show that they can be obtained from a hydrodynamic picture based on ballistically propagating quasiparticles. Going beyond the hydrodynamic limit, we analyze the approach to the non-equilibrium steady state and find that the leading late time corrections display a lattice effect. We also study the fine structure of the propagating fronts which are found to be described by the Airy kernel and its derivatives. Near the front we observe the phenomenon of energy back-flow where the energy locally flows from the colder to the hotter region.
We study the light-cone and front dynamics of a single-particle continuous time extended quantum walk on a one-dimensional lattice with finite range hopping. We show that, in general, for an initially localized state, propagating wave fronts can be characterized as ordinary or extremal fronts with the latter exhibiting an anomalous subdiffusive scaling behavior in the front region. We investigate the dynamical global and local scaling properties of the cumulative probability distribution function for the extended walk with nearest- and next-nearest-neighbor hopping using analytical and numerical methods. The global scaling shows the existence of a ``causal light-cone'' corresponding to excitations traveling with a velocity smaller than a maximal ``light velocity.'' Maximal fronts moving with fixed light velocity bound the causal cone. The front regions spread with time subdiffusively exhibiting a local Airy scaling which leads to an internal staircase structure. At a certain critical next-nearest-neighbor hopping strength, there is a transition from a phase with one `causal cone' to a phase with two nested `causal cones' and the existence of an internal staircase structure in the corresponding cumulative distribution profiles. We also connect the study to that in spin chain systems and indicate that a single-particle quantum walk on the one-dimensional lattice already captures the many-body physics of a spin chain system. In particular, we suggest that the time evolution of a single-particle quantum walk on the one-dimensional lattice with an initially localized state is equivalent to the time evolution of a domain-wall initial state in a corresponding spin chain system.
The Lieb-Robinson correlation function can be used to characterize the propagation of quantum information in Ising chains. Considerable work has been done to establish bounds on this correlation function in various circumstances. To actually calculate the value of the correlation function directly typically requires a state space which grows exponentially with system size, and so is intractable for all but relatively small systems. We employ a recently-developed method that enables direct calculation of the Lieb--Robinson correlation function and which scales linearly with system size. This enables the computation for systems with more than a thousand qubits, revealing the propagation of quantum information down the chain. We extend this technique to the problem of Ising chains with randomly disordered coupling strengths. Increasing disorder causes localization of the quantum correlations and halts propagation of quantum information down the chain. We see light-cone isocontours therefore becoming vertical. The decay of correlations down the chain is characterized by a correlation length which in the paramagnetic phase increases quadratically with mean coupling strength and inverse quadratically with disorder level. The correlation length peaks at the phase transition. We map our general expression for the correlation function in the Ising system onto the probability density in the single particle tight-binding Su-Schrieffer-Heeger model for a specific initial condition.
The Lieb-Robinson correlation function is the norm of a commutator between local operators acting on separate subsystems at different times. This provides a useful state-independent measure for characterizing the specifically quantum interaction between spatially separated qubits. The finite propagation velocity for this correlator defines a “light cone” of quantum influence. We calculate the Lieb-Robinson correlation function for one-dimensional qubit arrays described by the transverse-field Ising model. Direct calculations of this correlation function have been limited by the exponential increase in the size of the state space with the number of qubits. We introduce a technique that avoids this barrier by transforming the calculation to a sum over Pauli walks which results in linear scaling with system size. We can then explore propagation in arrays of hundreds of qubits and observe the effects of the quantum phase transition in the system. We observe the emergence of two distinct velocities of propagation: a correlation front velocity, which is affected by the phase transition, and the Lieb-Robinson velocity which is not. The correlation front velocity is equal to the maximum group velocity of single quasiparticle excitations. The Lieb-Robinson velocity describes the extreme leading edge of correlations when the value of the correlation function itself is still very small. For the semi-infinite chain of qubits at the quantum critical point, we derive an analytical result for the correlation function. Published by the American Physical Society 2024
We study the dynamics of the entanglement spectrum, that is the time evolution of the eigenvalues of the reduced density matrices after a bipartition of a one-dimensional spin chain. Starting from the ground state of an initial Hamiltonian, the state of the system is evolved in time with a new Hamiltonian. We consider both instantaneous and quasi adiabatic quenches of the system Hamiltonian across a quantum phase transition. We analyse the Ising model that can be exactly solved and the XXZ for which we employ the time-dependent density matrix renormalisation group algorithm. Our results show once more a connection between the Schmidt gap, i.e. the difference of the two largest eigenvalues of the reduced density matrix and order parameters, in this case the spontaneous magnetisation.
We study the relaxation behaviour of the quantum Ising chain, focusing our attention onto the non-equilibrium dynamics of the transverse magnetization. The initial states, from which the magnetization relaxes, are product states such as those generated by setting in contact several systems, each initially equilibrated at a given temperature. Due to the free fermionic structure of the chain, the dynamics of the transverse magnetization is easily expressed in a compact form. In the completely factorized initial state, corresponding to a situation where all the spins are thermalized independently, we obtain in the scaling limit the Green function associated to the transverse magnetization. The relaxation behaviour is also considered in the system-bath case, where part of the chain called the system is thermalized at a temperature Ts and the remaining part is at a temperature Tb. The magnetization profiles show a scaling behaviour. Moreover, in the extreme case Tb=∞ and Ts=0, it is shown that the magnetization relaxes in quantized steps in the strong transverse field region.