We investigate the dynamical spreading of spatial correlations after a quantum quench starting from a magnetically disordered state in the transverse-field Ising model at one (1D) and two spatial dimensions (2D). We analyze specifically the longitudinal and transverse spin-spin correlation functions at equal time with use of several methods. From the comparison of the results in 1D obtained by the linear spin-wave approximation (LSWA) and those obtained by the rigorous analytical approach, we show that the LSWA can asymptotically reproduce the exact group velocity in the limit of strong transverse fields while it fails to capture the detailed time dependence of the correlation functions. By applying the LSWA to the 2D case, in which the rigorous analytical approach is unavailable, we estimate the propagation velocity to be $Ja/(2\ensuremath{\hbar})$ at the strong-field limit, where
J is the Ising interaction and
a is the lattice spacing. We also utilize the tensor-network method based on the projected-entangled pair states for 2D and quantitatively compute the time evolution of the correlation functions for a relatively short time. Our findings provide useful benchmarks for quantum simulation experiments of correlation spreading and theoretical refinement of the Lieb-Robinson bound in the future.
In this master thesis, 3 methods of solving of Quantum Ising model have been outlined and compared, they include ABC. We investigate the entanglement dynamics and find energy and spin values for 1-dimensional Quantum Ising chain and 2 and 3-dimensional Quantum Ising models. System starts in a paramagnetic regime and switches to ferromagnetic spin-spin coupling constant J. In this work, we consider instantaneous (step function) switch, which demonstrates the instantaneous switch from paramagnetic to ferromagnetic state. Solution of the system is made by quantum simulation using Microsofts’ Corporations computational simulator called The Language-Integrated Quantum Operations (LIQUi|>) Simulator, 1/Z expansion method and Exact Diagonalization method. Accuracy and precision of results of first two methods are compared with Exact Diagonalization that provides the exact solution of the system. With that it can be seen to what extent 1/Z expansion method is applicable for solving of Quantum Ising model.
We study the dynamics of quantum correlations in a class of exactly solvable Ising-type models. We analyze in particular the time evolution of initial Bell states created in a fully polarized background and on the ground state. We find that the pairwise entanglement propagates with a velocity proportional to the reduced interaction for all the four Bell states. Singletlike states are favored during the propagation, in the sense that tripletlike states change their character during the propagation under certain circumstances. Characteristic for the anisotropic models is the instantaneous creation of pairwise entanglement from a fully polarized state; furthermore, the propagation of pairwise entanglement is suppressed in favor of a creation of different types of entanglement. The ``entanglement wave'' evolving from a Bell state on the ground state turns out to be very localized in space time. Our findings agree with a recently formulated conjecture on entanglement sharing; some results are interpreted in terms of this conjecture.
We study dynamics of entanglement in the thermalization process of an isolated quantum many-body system. We propose a simple setup for measuring the propagation speed of entanglement entropy (EE) in numerical simulations and apply it to the integrable/non-integrable spin models in 1D - the transverse Ising (TI) model, the chaotic Ising (CI) model, and the extended chaotic Ising (ECI) model. We find that two distinct time-scales
t∗ and
tdiff arise in the dynamics of EE in the thermalization process: the former represents the time-scale for the saturation of EE and the latter characterizes spreading of EE over the entire system. Evaluating the propagation speed of entanglement from
tdiff, we find that entanglement propagates ballistically with a constant velocity irrespective of the integrability of the model. The propagation speed of entanglement is found to coincide with the maximum group velocity of quasi-particle excitations in the TI model. We also evaluate the propagation speed of entanglement by mutual information and find the characteristic time-scale
tMI. We show that the propagation speeds of entanglement evaluated by
tMI and
tdiff agree well. We discuss the condition for thermalization based on the numerical results and propose that scrambling of the entire system has to take place before saturation of EE for thermalization.
Abstract We study dynamics of quantum entanglement in smooth global quenches with a finite rate, by computing the time evolution of entanglement entropy in 1 + 1 dimensional free scalar theory with time-dependent masses which start from a nonzero value at early time and either crosses or approaches zero. The time-dependence is chosen so that the quantum dynamics is exactly solvable. If the quenches asymptotically approach a critical point at late time, the early-time and late-time entropies are proportional to the time and subsystem size respectively. Their proportionality coefficients are determined by scales: in a fast limit, an initial correlation length; in a slow limit, an effective scale defined when adiabaticity breaks down. If the quenches cross a critical point, the time evolution of entropy is characterized by the scales: the initial correlation length in the fast limit and the effective correlation length in the slow limit. The entropy oscillates, and the entanglement oscillation comes from a coherence between right-moving and left-moving waves if we measure the entropy after time characterized by the quench rate. The periodicity of the late-time oscillation is consistent with the periodicity of the oscillation of zero modes which are zero-momentum value of two point functions of a fundamental field and its conjugate momentum.
Quantum Ising model in one dimension is an exactly solvable example of a quantum phase transition. We investigate its behavior during a quench caused by a gradual turning off of the transverse bias field. The system is then driven at a fixed rate characterized by the quench time ${\ensuremath{\tau}}_{Q}$ across the critical point from a paramagnetic to ferromagnetic phase. In agreement with Kibble-Zurek mechanism (which recognizes that evolution is approximately adiabatic far away, but becomes approximately impulse sufficiently near the critical point), quantum state of the system after the transition exhibits a characteristic correlation length $\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{\ensuremath{\xi}}$ proportional to the square root of the quench time ${\ensuremath{\tau}}_{Q}$: $\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{\ensuremath{\xi}}=\sqrt{{\ensuremath{\tau}}_{Q}}$. The inverse of this correlation length is known to determine average density of defects (e.g., kinks) after the transition. In this paper, we show that this same $\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{\ensuremath{\xi}}$ controls the entropy of entanglement, e.g., entropy of a block of
L spins that are entangled with the rest of the system after the transition from the paramagnetic ground state induced by the quench. For large
L, this entropy saturates at $\frac{1}{6}\phantom{\rule{0.2em}{0ex}}{\mathrm{log}}_{2}\phantom{\rule{0.2em}{0ex}}\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{\ensuremath{\xi}}$, as might have been expected from the Kibble-Zurek mechanism. Close to the critical point, the entropy saturates when the block size $L\ensuremath{\approx}\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{\ensuremath{\xi}}$, but---in the subsequent evolution in the ferromagnetic phase---a somewhat larger length scale $l=\sqrt{{\ensuremath{\tau}}_{Q}}\phantom{\rule{0.2em}{0ex}}\mathrm{ln}\phantom{\rule{0.2em}{0ex}}{\ensuremath{\tau}}_{Q}$ develops as a result of a dephasing process that can be regarded as a quantum analog of phase ordering, and the entropy saturates when $L\ensuremath{\approx}l$. We also study the spin-spin correlation using both analytic methods and real time simulations with the Vidal algorithm. We find that at an instant when quench is crossing the critical point, ferromagnetic correlations decay exponentially with the dynamical correlation length $\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{\ensuremath{\xi}}$, but (as for entropy of entanglement) in the following evolution length scale
l gradually develops. The correlation function becomes oscillatory at distances less than this scale. However, both the wavelength and the correlation length of these oscillations are still determined by $\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{\ensuremath{\xi}}$. We also derive probability distribution for the number of kinks in a finite spin chain after the transition.