# Power-law relaxation behavior of an initially localized state in the spin-1/2 Heisenberg chain

## Abstract

We present power-law relaxation behavior of the local magnetizations in the equilibration dynamics of the spin-1/2 Heisenberg spin chain as an isolated integrable quantum system. We perform the exact time evolution of the expectation values of the local spin operators by evaluating them with the determinant formula of the form factors. We construct such an initial quantum state that has a localized profile of the local magnetizations, and perform the exact time evolution over a very long period of time. We show that the local magnetization relaxes as some power of the time variable with no definite time scale, while the fidelity relaxes very fast with its relaxation time being proportional to the inverse of the energy width, i.e. the Boltzmann time.

###### pacs:

75.10.Pq, 02.30.Ik, 05.70.Ln## I Introduction

The non-equilibrium dynamics of isolated integrable quantum systems such as quantum quench has attracted much attention, recently (1); (2); (3); (4); (5); (6); (7); (8); (9); (10); (11); (12); (13); (14); (15). Due to novel experimental demonstrations of out-of-equilibrium dynamics of closed quantum systems made of cold atoms (16); (17); (18), it has become an extensive research subject how thermalization occurs in isolated quantum systems (3). For integrable quantum systems, which have a large number of conserved quantities, equilibration dynamics is not trivial. It has been conjectured that integrable quantum systems after long time the expectation values of local quantities approach those of the Generalized Gibbs Ensemble (GGE) (2). It has been investigated for particular integrable models and initial states whether the GGE conjecture holds or not (10); (11); (12); (13); (14); (15).

The studies of non-equilibrium quantum dynamics are closely related to the recent renewed interest in the foundation of quantum statistical mechanics (19); (25); (20); (21); (22); (23); (24); (26). One of the fundamental questions is whether and how equilibrium distribution functions are realized only through the unitary time evolution of a quantum state in an isolated system (27). It has been demonstrated that in some general or concrete settings a pure initial state of an isolated quantum system equilibrates or thermalizes in a certain mathematical sense (19); (25); (20); (21); (28); (29); (30); (31). It is related to the idea that a typical pure state of a macroscopic quantum system can fully describe thermal equilibrium (22); (23); (24); (26). The question is often discussed in terms of the eigenstate thermalization hypothesis (32); (33).

Although the above mentioned studies show that a wide class of isolated quantum systems thermalize or equilibrates to GGE, the time scale for thermalization has been studied only very recently (28); (29). Furthermore, the quantum dynamics of a pair of spins have been performed numerically (34); (35) in an experimental setting (18). However, it has not been explicitly shown how the expectation value of a certain local operator of an interacting integrable quantum system such as the Heisenberg chain of finite size equilibrates in time. It has not been clear even whether it approaches some constant value or not. Here we do not take the time average of the expectation value. It is natural to conjecture that as the system size increases, the expectation value of a local operator at a given time becomes close to the limiting value of the infinite system. However, we expect that there are many aspects in the way how it approaches the equilibrium value in time and they are related to interesting topics of non-equilibrium dynamics close to integrable systems such as prethermalization (36) and pre-relaxation (37), which should occur between the integrability and its breaking under non-integrable perturbations. Moreover, it is important to give a concrete example of equilibration of a local operator in a long period of time for interacting quantum systems from the viewpoint of stability and equilibrium properties of KMS states (38).

Let us introduce the Hamiltonian of the anti-ferromagnetic Heisenberg spin chain under the periodic boundary conditions (PBC): for .

(1) |

Here, denote the Pauli matrices. We also call the system the XXZ chain. In the present paper we show the exact relaxation dynamics of local magnetizations of the spin-1/2 Heisenberg chain for an initially localized quantum state. We also study it for other states constructed from spinon states. Although it shows a strong oscillating behavior in time, we argue that the amplitudes of oscillation show a power-law decay, by evaluating the square deviations of the local magnetizations. We suggest the power-law decay may be universal for the expectation values of other local operators while the decaying exponent depends on the initial state.

The paper consists of the following. We briefly review the Bethe-ansatz equations and the scheme to evaluate local magnetizations, which is due to the recent development of the algebraic Bethe-ansatz (40); (39). We first show that the relaxation time of the fidelity is given by the Boltmann time and is consistent with recent rigorous results. Here, the fidelity is given by the square amplitude of correlation between the initial state and state at time : . We then show the time evolution of the profile of local magnetizations , and observe how initially localized profile collapses in time. It seems that the resulting oscillations of local magnetizations continue forever. However, we shall show that the square deviations of the local magnetization exhibits a power-law decay.

## Ii Method

The eigenvectors and eigenvalues with down spins of the Hamiltonian are constructed from solutions of the Bethe-ansatz equations (BAE) (41); (42).

(2) | |||||

Here, quantum numbers are given by integers or half-integers according to the condition: . Here we take the branch: for any . We denote by the eigenvector assoiated with a solution of the BAE, . The energy eigenvalue of the eigenstate is given by

(3) |

We consider a state consisting of a superposition of the Bethe-ansatz eigenvectors

(4) |

Here we have abbreviated eigenstate by and denotes a certain set of solutions of BAE.

We denote by the local magnetization on the th site at time , and define it by the expectation value of the local spin operator with respect to the quantum state at time as follows.

(5) |

Here the time evolution of the state is given by where is given by . We evaluate the local magnetization by the form factor expansion:

(6) |

Here we recall that denotes the energy of the Bethe state .

In order to evaluate the matrix elements of the local operator between two Bethe eigenstates and , we make use of the determinant formula of the matrix elements of derived by Kitanine, Maillet and Terras (39). For numerical evaluation, we further divide it out by the Cauchy determinant similarly as that of the 1D Bose gas (40).

We remark that in order to show equilibration of local observables systematically, we may consider the Bloch vector, which is given by the expectation values of the traceless operators in the Hilbert space of the XXX spin chain with respect to a given quantum state (23). We consider the local magnetizations at sites , for are elements of the Bloch vector.

## Iii Fidelity

Let us evaluate the fidelity for pure states which are given the sums over spinon states with equal weight in the energy width for , 0.05 and 0.1, etc., specifically. Each spinon eigenstate corresponds to a real solution of BAE (2) specified by putting two holes in the set of quantum integers s (41); (42). Here we assume that the system size is given by , and the number of down-spins by half the system size minus 1: .

(7) |

where symbol denotes a subset of the set of all spinon eigenstates, and in the case of equal weight we set , where denotes the number of elements in the set . We plot in Fig. 1 the time evolution of the fidelity for the three states after initial time with for , 0.05 and 0.1: .

We observe that the three graphs of fidelity versus time shown in Fig. 1 are well approximated by the following expression (43)

(8) |

There is only one fitting parameter , practically. Here we assume that the inverse temperature is given by the energy width: . The initial Gaussian behavior should be consistent with other approaches (44); (45). Non-exponential behavior in the fidelity of interacting many-body systems is recently argued (46).

We estimate the relaxation time of fidelity through the numerical plots in Fig. 1. In each graph we determine it by the time when the fidelity takes the minimum value first after initial time. The red dots in Fig. 1 show the points of time when we determine the relaxation time .

The estimates of relaxation time of fidelity are plotted in Fig. 2. They are consistent with the line: . Here we remark that in the units of the present paper we take . Thus, it is given by the Boltzmann time assuming . Furthermore, it is interesting to observe that it coincides with a rigorous theoretical estimate (29) of thermalization time in a macroscopic quantum system for a typical non-equilibrium subspace: .

For the XXX chain of we evaluated the fidelity for considering all the Bethe-ansatz eigenstates associated with not only real solutions but also complex solutions as in Ref. (47) (see also (48)). We compared the time evolution of fidelity for the state consisting of only real solutions and that for the state of all solutions including complex solutions. However, there is no clear difference in the time evolution of fidelity between the states of only real solutions and those of all solutions.

## Iv Local magnetization

Let us now consider mainly the all-spinon state , which is given by the sum over all the spinon states for the XXX spin chain with an even system size in the sector of down-spins, i.e. in the case of an almost half-filled lattice. We evaluate the local magnetization at time by eq. (6).

We observe in Fig. 3 that the profile of local magnetizations for at initial time is localized around at the middle site of and almost vanishes on other sites. We may consider it as the “positive” analogue of a “quantum dark soliton” constructed in the one-dimensional Bose gas with delta-function interactions (31), which are given by the sum over the type-II excitations in a branch. Here we remark that the average of is given by for .

There are several procedures at different time steps as shown in Fig. 3 for the whole time evolution in which the initially localized profile of local magnetizations collapses to a uniformly oscillating profile in the final stage. First, the initially-localized profile separates to two localized waves moving in different directions. Then, the two localized waves propagate and collide each other at around at the 1st site. They further propagate with smaller height of localized waves, and finally they merge into such a profile that is roughly uniform in space but strongly oscillating in time.

Thus, the initially localized profile in Fig. 3 dynamically collapses to the roughly uniform profile oscillating in time. Here, the exact time evolution of the local magnetization at a site shows strong oscillating behavior in Fig. 4. It seems that the fluctuations of the local magnetizations do not easily relax to the order of the inverse of the number of the Bethe-ansatz eigenstates consisting of the quantum state , i.e. , which is the order of magnitude expected statistically.

However, by taking the site-average of the squared deviations of the local magnetizations over all sites we observe a power law decay over a long period of time as in Fig.5. We define the site-average of squared deviations of local magnetizations on the th sites by

(9) |

Here we remark that due to the symmetry of the Hamiltonian the sum of all the local magnetization is given by a constant: .

For the all-spinon state of , the squared deviations of local magnetizations decay almost as an inverse of time initially (), and as an inverse power of time with a smaller exponent (). For the all-spinon state of the squared deviations of local magnetizations are shown in Fig. 5. For various other states which are given by the sums over sets of spinons, we observed similar power-law decay behavior of squared deviations (9). Here, the estimates of the exponent are different among the states. For instance, the yrast state, which is given by the sum over lowest excitations among spinon states, the estimate of the exponent is close to zero. We suggest that power-law decay may be universal in the time evolution of the expectation values of local operators for interacting integrable quantum systems, while the exponent depends on the states.

We now argue that the time scale for equilibration of the local magnetizations is very long. We suggest that the time scale of the collapse of the initially localized magnetization profile in the first stage is given by the system size divided by the spinon velocity. It is proportional to the system size. Here we recall that two localized waves travel in time in opposite directions in Fig. 3. However, the power-law decay of the squared deviations of the local magnetizations shown in Fig. 5 shows that there is no definite relaxation time in the whole long-time equilibration dynamics of the local magnetizations.

In summary we have performed exact relaxation dynamics of the local magnetization profile for the initially localized state in the XXX chain. Through the square deviations of local magnetizations we presented the power-law decay of the local magnetization in a long time. We suggest that the expectation values of other local operators may also decay as a power of the time variable for interacting integrable quantum systems.

## Acknowledgement

The authors would like to thank E. Kaminishi, S. Moriya and J. Sato for useful comments. The present study is partially supported by Grant-in-Aid for Scientific Research No. 15K05204.

### References

- E. Barouch, B.M. McCoy and M. Dresden, Phys. Rev. 2, 1075 (1970).
- M. Rigol, V. Dunjko, V. Yurovsky and M. Olshanii, Phys. Rev. Lett 98, 050405 (2007);
- M. Rigol, V. Dunjko and M. Olshanii, Nature 452, 854 (2008).
- D. Rossini, S. Suzuki, G. Mussardo, G.E. Santoro and A. Silva, Phys. Rev. B 82, 144302 (2010).
- J. Mossel and J.-S. Caux, New J. Phys. 12, 055028 (2010).
- F. Iglói and H. Rieger, Phys. Rev. Lett. 106, 035701 (2011).
- A.C. Cassidy, C.W. Clark, and M. Rigol, Phys. Rev. Lett. 106, 140405 (2011).
- P. Calabrese, F.H.L. Essler and M. Fagotti, Phys. Rev. Lett. 106, 227203 (2011); J. Stat. Mech. (2012) P07016; J. Stat. Mech. (2012) P07022.
- J.-S. Caux and F.H.L. Essler, Phys. Rev. Lett. 110, 257203 (2013).
- M. Brockmann, B. Wouters, D. Fioretto, J. De Nardis, R. Vlijm and J.-S. Caux, J. Stat. Mech., P12009 (2014).
- B. Pozsgay, M. Mestyán, M.A. Werner, M. Kormos, G. Zaránd and G. Takács, Phys. Rev. Lett. 113, 117203 (2014).
- B. Wouters, J. De Nardis, M. Brockmann, D. Fioretto, M. Rigol, and J.-S. Caux, Phys. Rev. Lett. 113, 117202 (2014).
- M. Mierzejewski, P. Prelov̌sek and T. Prozen, Phys. Rev. Lett. 114, 140601 (2015).
- E. Ilievski, M. Medenjak and T. Prosen, arXiv:1506.05049 .
- E. Ilievski, J. De Nardis, B. Wouters, J.-S. Caux, F.H.L. Essler, and T. Prosen, arXiv:1507.02993 .
- T. Kinoshita, T. Wenger and D.S. Weiss, Science 305, 1125 (2004); Phys. Rev. Lett. 95, 190406 (2005); Nature 440, 900 (2006).
- M. Gring, M. Kuhnert, T. Langen, T. Kitagawa, B. Rauer, M. Schreitl, I. Mazets, D. Adu Smith, E. Demler, J. Schmiedmayer, Science 337, 1318 (2012).
- T. Fukuhara, P. Schauß, M. Endres, S. Hild, M. Cheneau, I. Bloch and C. Gross, Nature 502, 76 (2013).
- H. Tasaki, Phys. Rev. Lett. 80, 1373 (1998).
- S. Goldstein, J.L. Lebowitz, C. Mastrodonato, R. Tumulka, and N. Zanghi, Phys. Rev. E 81, 011109 (2010).
- P. Reimann and M. Kastner, New J. Phys. 14, 043020 (2012)
- S. Goldstein, J.L. Lebowitz, R. Tumulka, and N. Zanghi, Phys. Rev. Lett. 96, 050403 (2006)
- A. Sugita, Nonlinear Phenom. Complex Syst. 10, 192 (2007).
- P. Reimann, Phys. Rev. Lett. 99, 160404 (2007).
- P. Reimann, Phys. Rev. Lett. 101, 190403 (2008).
- S. Sugiura and A. Shimizu, Phys. Rev. Lett. 108, 240401 (2012); Phys. Rev. Lett. 111, 010401 (2013).
- R. Tumulka, Eur. Phys. J. H 35, 201 (2010) [J. von Neumann, Z. Phys. 57, 30 (1929)]
- S. Goldstein, T. Hara and H. Tasaki, Phys. Rev. Lett. 111, 140401 (2013).
- S. Goldstein, T. Hara and H. Tasaki, New J. Phys. 17, 045002 (2015)
- R.V. Jensen and R. Shanker, , Phys. Rev. Lett. 54, 1879 (1985).
- J. Sato, R. Kanamoto, E. Kaminishi, and T. Deguchi, Phys. Rev. Lett. 108, 110401 (2012).
- J. M. Deutsch, Phys. Rev. A 43, 2046 (1991).
- M. Srednicki, Phys. Rev. E 50, 888 (1994).
- M. Ganahl, E. Rabel, F. H. L. Essler, and H. G. Evertz, Phys. Rev. Lett. 108, 077206 (2012).
- W. Liu and N. Andrei, Phys. Rev. Lett. 112, 257204 (2014).
- J. Berges, S. Borsanyi, C. Wetterrich, Phys Rev. Lett. 93, 142002 (2004)
- B. Bertini and M. Fagotti, arXiv:1501.07260
- O. Bratteli and D.W. Robinson, Operator Algebras and Quantum Statistical Mechanics 2: Equilibrium States (Berlin Heidelberg New York: Springer-Verlag) (1997).
- N. Kitanine, J.M. Maillet and V. Terras, Nucl. Phys. B 554 [FS] (1999) 647–678
- N. A. Slavnov, Teor. Mat. Fiz. 79, 232 (1989) ; 82, 389 (1990);
- V. Korepin, N. Bogoliubov and A. Izergin, “Quantum Inverse Scattering Method and Correlation Functions” (Cambridge: Cambridge University Press) (1993).
- M. Takahashi, “Thermodynamics of One-Dimensional Solvable Models”, (Cambridge: Cambridge University Press) (1999).
- T. Monnai, J. Phys. Soc. Jpn. 83, 064001 (2014).
- T. Prosen, Phys. Rev. E 65, 036208 (2002).
- E. J. Torres-Herrera and Lea F. Santos, Phys. Rev. A 90, 033623 (2014)
- A. del Campo, arXiv:1504.01620 .
- R. Hagemans and J.-S. Caux, J. Phys. A: Math. Theor. 40 14605-14647 (2007).
- P. R. Giri and T. Deguchi, J. Stat. Mech., P07007 (2015)