# Spin transport in a one-dimensional anisotropic Heisenberg model

###### Abstract

We analytically and numerically study spin transport in a one-dimensional Heisenberg model in linear-response regime at infinite temperature. It is shown that as the anisotropy parameter is varied spin transport changes from ballistic for to anomalous at the isotropic point , to diffusive for finite , ending up as a perfect isolator in the Ising limit of infinite . Using perturbation theory for large a quantitative prediction is made for the dependence of diffusion constant on .

###### pacs:

05.60.Gg, 75.10.Pq, 05.30.-d, 03.65.Yz, 05.70.LnThe one-dimensional spin Heisenberg model is one of the oldest quantum models Heisenberg:28 being also the simplest model of interacting quantum particles. Despite being exactly solvable by the Bethe ansatz Bethe:31 , calculating its dynamical properties, like transport, is by no means simple. Understanding transport in the Heisenberg model is important for several reasons. First, it is still not known what are the necessary requirements for a system to display phenomenological transport laws where a current is proportional to the gradient of the driving field. Second motivation comes from the condensed matter where one would like to understand transport in strongly-correlated electron systems. A paradigmatic model is the Hubbard model, thought to be related to the problem of high- superconductivity, whose low energy excitations can be described by the antiferromegnetic Heisenberg model. The Heisenberg model is with a very high accuracy realized also in the so-called spin-chain materials spin-chain . One-dimensional spin chains are receiving increased attention also from the cold-atoms community, where they can be experimentally realized cold . Particularly notorious to understand is the transport of magnetization. It has been actively studied for more than 20 years, yet the subject is still hotly debated. In particular, numerical calculations sometimes give conflicting results while there are almost no rigorous statements, exceptions are a finite Drude weight at Shastry:90 and recently also for Prosen:11 , both for , signaling ballistic transport. Various, mostly numerical approaches, range from the Mazur’s inequality Zotos:97 ; Prosen:11 , the Bethe ansatz calculation bethe , exact diagonalization exactdiag , quantum Monte Carlo Alvarez:02 , Lanczos method Samir:04 , wave-packet evolution Langer:09 , Luttinger liquid theory Affleck:09 , master equation JSTAT:09 , and correlation functions Steinigeweg:09 . Because the model is solvable one would be tempted to think that it is ballistic Castella:95 . However, recently a solvable diffusive quantum model has been found XXdephasing . There is also mounting evidence JSTAT:09 ; Langer:09 ; Steinigeweg:09 that it is diffusive for .

The results presented in the present paper, together with other recent works, enable us to give a complete picture of spin transport in the linear-response regime in the Heisenberg model. For the model is ballistic at infinite temperature Prosen:11 as well as at any finite or zero Shastry:90 temperature. At our results show anomalous transport at infinite temperature. It is also plausible to expect anomalous behavior at finite temperatures. For and infinite temperature we show diffusive transport. Recent work Jesenko:11 shows that as one decreases temperature, the diffusion constant increases, possibly exponentially fast in . At temperatures lower than the gap transport trivially stops. As one lowers towards the isotropic point the diffusion constant diverges at any temperature, see also Jesenko:11 .

To describe a nonequilibrium situation we couple boundary spins of the chain to magnetization reservoirs. Time evolution of the density matrix describing the chain evolves according to the Lindblad master equation,

(1) |

where the dissipative linear operator describing bath is expressed in terms of Lindblad operators , . Reservoirs are realized by two Lindblad operators at each end, at the left end and at the right end, . We fix . In the diffusive regime for the diffusion constant does not depend neither on the value of nor on the detailed Lindblad operators used to model the bath otherL . Our Lindblad operators are such that they always induce a nonequilibrium steady state (NESS) with almost zero average energy density, which therefore corresponds to an infinite temperature temp . The Hamiltonian is

(2) |

Small – In this section we are interested in small , so that the term acts as a perturbation of the XX model. Analytical perturbative results will serve as a spring-board for the discussion of transport for not so small values of . We split the Lindblad superoperator into two parts, , where is the part for , while is the perturbation, i.e., . For the NESS solution of this master equation, denoted by , is nondegenerate and ballistic and can be neatly written in a matrix product operator form with matrices of fixed dimension JPA10 . Perturbation series ansatz for the NESS for small is . We want to calculate the first two orders in . Also, in all terms we will be interested only in the lowest order terms in the driving , meaning that we study linear-response behavior. Plugging this ansatz into NESS equation , using , and then equating terms with the same order in , gets us two matrix equations, and . A known zeroth-order solution can be used to get inhomogeneous terms in the linear equations for , which can in turn be used to get . The only problem is that the number of linear equations is exponentially large in the length of the chain . Nevertheless, several general remarks can be made because we know (it contains only , spin current , and their products) and the action of and : the first order term does not contain any magnetization or spin current ; they appear only in the 2nd order term . Therefore, to the lowest order in perturbation the current and the magnetization profile do not change! They are the same as for the ballistic XX model at . Because the perturbation of a nondegenerate NESS is nonsingular, fixing , one can always find a sufficiently small such that the perturbative expansion will converge and the system is ballistic. However, for the transport in the thermodynamic limit the relevant order of limits is first fixing and and only then sending . This limit is more difficult to treat because higher order terms, for instance , can grow with faster than the first order term , causing the convergence radius to shrink as . As we will see, this is indeed what happens.

We have found the exact expressions for and for small ; the coefficients in front of all the terms are rational functions with a rather large denominators and as such not very transparent. We therefore do not give their precise form here but rather focus on their scaling with in order to infer the convergence radius of the perturbative expansion. Because we are interested in the spin transport, we need to know the behavior of the correction in that involves and spin current . It turns out that the correction in the current is the same for all sites while the correction to the magnetization depends on the position. For the spin current, the coefficient in front of the term is , for , for , while it is for and for . If one looks at the dependence of these coefficients on one notices that it is to a very good approximation linear. Fitting gives the dependence , with deviations being possibly exponentially small in . Similar corrections, all growing linearly with are also found for (the one for is in fact equal to times the one for the current, for other spins prefactors are larger). To sum up, the expectation value of the spin current is to order (and to linear order in ) equal to . Because the 2nd order correction grows with the system size, the perturbative expansion holds only for , and therefore breaks down in the thermodynamic limit. Unfortunately, from our analytical calculation one therefore can not decide about the nature of the spin transport at finite foot1 .

For larger we used time-dependent density matrix renormalization (tDMRG) simulations JSTAT:09 with to get expectations of magnetization and spin current in the NESS. Results are in Fig. 1.

For the current saturates for sufficiently large , with the saturation current monotonically decreasing with . The system is therefore ballistic. For small the characteristic at which converges to a constant value scales as , which is the same as the scaling of the analytical perturbative result. Even though perturbation theory fails as , the scaling apparently caries over beyond the perturbative result.

Isotropic– At the isotropic point the spin current at fixed driving scales as , nicely seen in tDMRG data in Fig. 2. The isotropic Heisenberg model at infinite temperature therefore display anomalous diffusion, with the diffusion constant diverging as . This is the first observation of an anomalous diffusion in a coherent (Hamiltonian) quantum system.

Furthermore, the magnetization profile along the chain has a nice scaling with and . As can be seen in Fig. 2, the dimensionless scaling function looks to be very close to , however, deviations seen in the figure seem to be larger than the finite-size or numerical accuracy effects.

Large – Perhaps the most interesting regime is for where numerical calculations point to a diffusive transport at infinite temperature JSTAT:09 ; Steinigeweg:09 . Because analytical treatment seems to be hard, we choose to study the case of large , where the Hamiltonian is close to the Ising one. The case of large can be equivalently reformulated with the Hamiltonian , with foot2 . The NESS for is exponentially degenerate. It is easy to see that the eigenvector of a dissipative bath part is at the left end and at the last spin. Because of , any operator of the form , with being an arbitrary combination of and , is a zero eigenstate of the Lindblad superoperator, i.e., the NESS. There are independent states of this form. Besides these, there are additional NESS states, namely, those for which holds and where now includes also . This increases the degeneracy even further. Because of this high degeneracy perturbation theory is more difficult than for small . An important thing to note is that the NESS for such Ising chain can support an arbitrary magnetization profile, while the spin current is always zero. The Ising spin chain is therefore a perfect insulator. Exponentially high degeneracy can now be understood also as being due to the isolation of the bulk from the boundaries, so that spins in the bulk “do not know” about the reservoirs at the boundaries. High degeneracy is therefore generic and can not be removed by a different choice of Lindblad operators. Perturbation breaks this high degeneracy making the NESS nondegenerate. The gap between two eigenvalues of the Lindblad superoperator with the largest real parts scales as for small perturbations. This means that if we want to calculate the lowest order corrections in the NESS exactly, we have to expand it upto order . Perturbative expansion can be for small written as . The resulting linear equations for unknown and are , and , where we have split the superoperator into an unperturbed part and the perturbation . Using the appropriate ansatz, we have obtained exact solutions for small , however, they are again complicated, involving many terms. We only point out features important for the spin transport. The first observation is that can contain only terms that are already present in the NESS for . This includes magnetization, but not the spin current. Therefore, operators will be present in , while the spin current will be present only in the first order term (because will result in the current). Spin current is therefore always proportional to (), while the magnetization scales as . Writing out the equation involving the coefficient in front of the spin current in , one gets , where is the coefficient in front of in and is the coefficient in front of in . Our exact analytical solutions for small ’s show that the term is always equal to for . Therefore, for large one has an exact relation . This states that if the magnetization profile is linear on average, then the spin current scales as and the transport is diffusive. Of course, showing that might be no easier than showing this for the current. Exact solutions give the expectation value of the current as , for . To access the limiting value we have looked at the convergence of with . If the transport is diffusive, this coefficient should converge to the diffusion constant . In the Fig. 3 we plot the values of these exact coefficients, together with numerically obtained ones for . The scaling seems to be linear in enabling us to obtain the limiting value of the coefficient as .

Using this limit we can predict that the spin current goes as , resulting in the diffusion constant foot3

(3) |

The dependence of on has been discussed in Stein:10 . The 2nd order term (as well as the 3rd) does not contain any corrections to the current. From the analysis of the solutions for small it also does not appear that they would grow with . The convergence radius of the perturbative series is therefore finite and does not decrease as .

To verify the theoretical prediction for a diffusion constant (3) we have again performed tDMRG simulations for a range of as well as . For each a diffusive scaling of the current has been checked and the prefactor, being the diffusion constant, determined. The results are plotted in Fig.4, together with Eq.(3).

We can see a perfect agreement for large , with the perturbative result (3) holding upto quite small foot4 .

Conclusion– By using perturbation theory in the limit of small and large anisotropies as well as large-scale numerical simulations we have shown that for the Heisenberg model displays ballistic spin transport. At the isotropic point transport is anomalous, with the current scaling as . This is the first observation of an anomalous transport in a coherent quantum system and has strong implications for an unexplained high heat conductivity measured in spin-chain materials exper . For we show that the transport is diffusive and inversely proportional to for large anisotropies. Support by the Program P1-0044 and the Grant J1-2208 of the Slovenian Research Agency, the project 57334 by CONACyT, Mexico, and IN114310 by UNAM is acknowledged.

## References

- (1) W. Heisenberg, Zur Theorie des Ferromagnetismus, Z. Phys. 49, 619 (1928).
- (2) H. Bethe, Zur Theorie der Metalle I. Eigenwerte und Eigenfunktionen der linearen Atomkette, Z. Phys. A 71, 205 (1931).
- (3) A. V. Sologubenko, T. Lorenz, H. R. Ott, and A. Freimuth, Thermal conductivity via magnetic excitations in spin-chain materials, J. Low Temp. Phys. 147, 387 (2007); F. Heidrich-Meisner, A. Honecker, and W. Brenig, Transport in quasi one-dimensional spin-1/2 systems, Eur. Phys. J. Special Topics 151, 135 (2007).
- (4) S. Trotzky, P. Cheinet, S. Fölling, M. Feld, U. Schnorrberger, A. M. Rey, A. Polkovnikov, E. A. Demler, M. D. Lukin, and I. Bloch, Time-resolved observation and control of superexchange interactions with ultracold atoms in optical lattices, Science 319, 295 (2008); J. T. Barreiro, M. Müller, P. Schindler, D. Nigg, T. Monz, M. Chwalla, M. Hennrich, C. F. Roos, P. Zoller, and R. Blatt, An open-system quantum simulator with trapped ions, Nature 470, 486 (2011); J. Simon, W. S. Bakr, R. Ma, M. E. Tai, P. M. Preiss, and M. Greiner, Quantum simulation of antiferromagnetic spin chains in an optical lattice, Nature 472, 307 (2011).
- (5) B. S. Shastry and B. Sutherland, Twisted boundary conditions and effective mass in Heisenberg-Ising and Hubbard rings, Phys. Rev. Lett. 65, 243 (1990).
- (6) T. Prosen, Open XXZ spin chain: Nonequilibrium steady state and a strict bound on ballistic transport, Phys. Rev. Lett. 106, 217206 (2011).
- (7) X. Zotos, F. Naef, and P. Prelovšek, Transport and conservation laws, Phys. Rev. B 55, 11029 (1997); P. Jung and A. Rosch, Spin conductivity in almost integrable spin chains, Phys. Rev. B 76, 245108 (2007).
- (8) X. Zotos, Finite temperature Drude weight of the one-dimensional spin-1/2 Heisenberg model, Phys. Rev. Lett. 82, 1764 (1999); J. Benz, T. Fukui, A. Klümper, and C. Scheeren, On the finite temperature Drude weight of the anisotropic Heisenberg chain, J. Phys. Soc. Jpn. Supp. 74, 181 (2005).
- (9) X. Zotos and P. Prelovšek, Evidence for ideal insulating or conducting state in a one-dimensional integrable system, Phys. Rev. B 53, 983 (1996); B. N. Narozhny, A. J. Millis, and N. Andrei, Transport in the XXZ model, Phys. Rev. B 58, R2921 (1998); F. Heidrich-Meisner, A. Honecker, D. C. Cabra, and W. Brenig, Zero-frequency transport properties of one-dimensional spin-1/2 systems, Phys. Rev. B 68, 134436 (2003).
- (10) J. V. Alvarez and C. Gros, Low-temperature transport in Heisenberg chains, Phys. Rev. Lett. 88, 077203 (2002); S. Grossjohann and W. Brenig, Hydrodynamic limit for the spin dynamics of the Heisenberg chain from quantum Monte Carlo calculations, Phys. Rev. B 81, 012404 (2010).
- (11) P. Prelovšek, S. El Shawish, X. Zotos, and M. Long, Anomalous scaling of conductivity in integrable fermion systems, Phys. Rev. B 70, 205129 (2004).
- (12) S. Langer, F. Heidrich-Meisner, J. Gemmer, I. P. McCulloch, and U. Schollwöck, Real-time study of diffusive and ballistic transport in spin-1/2 chains using the adaptive time-dependent density matrix renormalization group method, Phys. Rev. B 79, 214409 (2009).
- (13) J. Sirker, R. G. Pereira, and I. Affleck, Diffusion and ballistic transport in one-dimensional quantum systems, Phys. Rev. Lett. 103, 216602 (2009).
- (14) T. Prosen and M. Žnidarič, Matrix product simulations of non-equilibrium steady states of quantum spin chains, J. Stat. Mech. P02035 (2009); M. Žnidarič, Dephasing-induced diffusive transport in the anisotropic Heisenberg model, New J. Phys. 12, 043001 (2010).
- (15) R. Steinigeweg and J. Gemmer, Density dynamics in translationally invariant spin-1/2 chains at high temperatures: A current-autocorrelation approach to finite time and length scales, Phys. Rev. B 80, 184402 (2009); R. Steinigeweg, H. Wichterich, and J. Gemmer, Density dynamics from current auto-correlations at finite time- and length-scales, Europhys. Lett. 88, 10004 (2009).
- (16) H. Castella, X. Zotos, and P. Prelovšek, Integrability and ideal conductance at finite temperatures, Phys. Rev. Lett. 74, 972 (1995).
- (17) M. Žnidarič, Exact solution for a diffusive nonequilibrium steady state of an open quantum chain, J. Stat. Mech., L05002 (2010); M. Žnidarič, Solvable quantum nonequilibrium model exhibiting a phase transition and a matrix product representation, Phys. Rev. E 83, 011108 (2011).
- (18) S. Jesenko and M. Žnidarič, Finite-temperature magnetization transport of the Heisenberg model, eprint arXiv:1105.6340.
- (19) We have checked that a 2-spin bath involving 16 Lindblad operators on the boundary two spins, and used before in JSTAT:09 , leads to essentially the same results at and . For the influence of Lindblad operators is more delicate and we do not discuss it here.
- (20) Energy density is in linear-response tDMRG simulations found to be close to zero (perturbative arguments show that is proportional to ). Temperature can be determined from the equilibrium state, which is for proportional to , or from expectation values of energy density and magnetization, as described in M. Žnidarič, T. Prosen, G. Benenti, G. Casati, and D. Rossini, Thermalization and ergodicity in one-dimensional many-body open quantum systems, Phys. Rev. E 81, 051135 (2010). In the linear-response regime the two methods agree.
- (21) M. Žnidarič, A matrix product solution for a nonequilibrium steady state of an XX chain, J. Phys. A 43, 415004 (2010).
- (22) In order to make sure that there are no faster growing terms in 3rd and higher orders we verified the perturbative prediction by tDMRG simulations, by which we can treat larger systems. It is found that the perturbative result indeed breaks down when .
- (23) The NESS in such formulation is exactly the same as for the Hamiltonian with the coupling strength (the coupling to reservoars should scale the same as the largest parameter in ).
- (24) Expression is for .
- (25) R. Steinigeweg and R. Schnalle, Projection operator approach to spin diffusion in the anisotropic Heisenberg chain at high temperatures, Phys. Rev. E 82, 040103(R) (2010).
- (26) A small discrepancy at largest can be attributed to finite-size effects as the numerical simulations get more difficult there and we could simulate only chains of order , whereas we could go upto at .
- (27) N. Hlubek, P. Ribeiro, R. Saint-Martin, A. Revcolevschi, G. Roth, G. Behr, G. Büchner, and C. Hess, Ballistic heat transport of quantum spin excitations as seen in , Phys. Rev. B 81, 020405(R) (2010).