Entropy of entanglement and correlations induced by a quench: Dynamics of a quantum phase transition in the quantum Ising model

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 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 xi proportional to the square root of the quench time tau(Q): xi= root 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 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 1/6 log(2) xi, as might have been expected from the Kibble-Zurek mechanism. Close to the critical point, the entropy saturates when the block size L approximate to xi, but-in the subsequent evolution in the ferromagnetic phase-a somewhat larger length scale l=root tau(Q) ln 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 approximate to 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 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 xi. We also derive probability distribution for the number of kinks in a finite spin chain after the transition.