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 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
The Lieb-Robinson correlation function is one way to capture the propagation of quantum entanglement and correlations in many-body systems. We consider arrays of qubits described by the tranverse-field Ising model and examine correlations as the expanding front of entanglement first reaches a particular qubit. Rather than a new bound for the correlation function, we calculate its value, both numerically and analytically. A general analytical result is obtained that enables us to analyze very large arrays of qubits. The velocity of the entanglement front saturates to a constant value, for which an analytic expression is derived. At the leading edge of entanglement, the correlation function is well-described by an exponential reduced by the square root of the distance. This analysis is extended to arbitrary arrays with general coupling and topologies. For regular two and three dimensional qubit arrays with near-neighbor coupling we find the saturation values for the direction-dependent Lieb-Robinson velocity. The symmetry of the underlying 2D or 3D lattice is evident in the shape of surfaces of constant entanglement, even as the correlations front expands over hundreds of qubits.
The Lieb-Robinson bound (LRB) states that the range and strength of interactions between the constituents of a complex many-body system impose upper limits to how fast the signal can propagate. It manifests in a light conelike growth of correlation function connecting two distant subsystems. Here, we employ the techniques of quantum information to demonstrate that the LRB can be determined from local measurements performed on a single qubit that is connected to a many-body system. This formulation provides an operational recipe for estimating the LRB in complex systems, replacing the measurement of the correlation function with simple single-particle manipulations. We demonstrate the potency of this approach by deriving the upper limit to the speed of signal propagation in the XY spin chain.
The nonequilibrium response of a quantum many-body system defines its fundamental transport properties and how initially localized quantum information spreads. However, for long-range-interacting quantum systems little is known. We address this issue by analyzing a local quantum quench in the long-range Ising model in a transverse field, where interactions decay as a variable power law with distance ∝r(-α), α>0. Using complementary numerical and analytical techniques, we identify three dynamical regimes: short-range-like with an emerging light cone for α>2, weakly long range for 1<α<2 without a clear light cone but with a finite propagation speed of almost all excitations, and fully nonlocal for α<1 with instantaneous transmission of correlations. This last regime breaks generalized Lieb-Robinson bounds and thus locality. Numerical calculation of the entanglement spectrum demonstrates that the usual picture of propagating quasiparticles remains valid, allowing an intuitive interpretation of our findings via divergences of quasiparticle velocities. Our results may be tested in state-of-the-art trapped-ion experiments.
Dynamical correlations of the spin and the energy density are investigated in the critical region of the random transverse-field Ising chain by numerically exact calculations in large finite systems (L ≤ 128). The spin-spin autocorrelation function is found to decay proportional to (log t) −2xm and (log t) −2xs m in the bulk and on the surface, respectively, with xm and x s m the bulk and surface magnetization exponents, respectively. On the other hand the critical energy-energy autocorrelation functions have a power law decay, which are characterized by novel critical exponents ηe ≈ 2.2 in the bulk and η s e ≈ 2.5 at the surface, respectively. The numerical results are compared with the predictions of a scaling theory. 05.50.+q, 64.60.Ak, 68.35.Rh Typeset using REVTEX 1 The asymptotic behavior of the time-dependent correlation functions for interacting many-body systems turned out a very difficult subject of theoretical research. Exact results in this field are scarce, one can mention the one-dimensional spin 1/2 XY-model [1] and the Ising chain in a transverse-field [2]. Both models can be mapped onto a system of non-interacting fermions, where the equal-position correlation functions are calculated by the Pfaffian method utilizing the theory of Töplitz determinants. In this Letter we consider- at the first time- the critical dynamical correlations of an interacting quantum system in the presence of quenched (i.e. time-dependent) disorder. It has recently become clear that quenched disorder has rather different effects on phase transitions in quantum systems [3] than on those thermally driven phase transitions. For example, in the Griffiths phase, which is situated at the disordered side of the critical point, the susceptibility has an essential singularity in classical systems, whereas in a quantum system the corresponding singularity is stronger, it is in a power law form. Here we consider the prototype of random quantum systems the one-dimensional random transverse-field Ising model defined by the Hamiltonian:
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.