# Thermal entanglement between non-nearest-neighbor spins on fractal lattices

## Abstract

We investigate thermal entanglement between two non-nearest-neighbor sites in ferromagnetic Heisenberg chain and on fractal lattices by means of the decimation renormalization-group (RG) method. It is found that the entanglement decreases with increasing temperature and it disappears beyond a critical value . Thermal entanglement at a certain temperature first increases with the increase of the anisotropy parameter and then decreases sharply to zero when is close to the isotropic point. We also show how the entanglement evolves as the size of the system becomes large via the RG method. As increases, for the spin chain and Koch curve the entanglement between two terminal spins is fragile and vanishes when , but for two kinds of diamond-type hierarchical (DH) lattices the entanglement is rather robust and can exist even when becomes very large. Our result indicates that the special fractal structure can affect the change of entanglement with system size.

###### pacs:

03.67.Mn, 73.43.Nq, 75.10.Pq, 64.60.ae^{1}

## I Introduction

An essential difference between quantum and classical physics is the possible existence of nonlocal correlation in quantum system which is called the entanglement 1 (). Recently, the quantum entanglement has been recognized as an crucial resource in various fields of quantum information such as quantum communication and computation 2 (); 3 (); 4 (); 5 (); 6 (); 7 (). Since the entanglement is fragile and sensitive to many environment factors, many efforts are devoted to studying stable entanglement for realistic system in finite temperature. Thus, thermal entanglement has naturally received much attention by its advantage of stability and requiring neither measurement nor controlled switching of interactions in the preparing process. For spin systems can be used for gate operation in quantum computer thermal entanglement on solid spin systems have been widely studied, for example, the Heisenberg spin chain both in the absence absence field () and presence presence field 1 (); presence field 2 () of an external magnetic field , spin rings spin rings () and spin clusters spin clusters 1 (); spin clusters 2 (). However, most of these works only focused on thermal entanglement between nearest, next-nearest or next-to-next-nearest neighbor spins 14 (); 15 (); 16 (); 17 (); 18 (); 19 (). This motivate us to propose two questions: (i) Can thermal entanglement exist between distant non-nearest-neighbor sites in spin system? (ii) How does thermal entanglement evolve as system size grows? But it is very difficult to obtain exact results on entanglement in spin systems on arbitrary lattices especially fractal lattices, since this usually requires the expression of the partition function which is too complicated to solve when the system size becomes very large.

In recent years, the entanglement at zero temperature in the large size system has been studied by adopting the renormalization-group (RG) method. In 2002, A. Osterloh et al first introduced the density matrix renormalization-group (DMRG) approach to study the entanglement close to the quantum phase transition (QPT) Nature () and reveal a profound difference between classical correlations and the non-local quantum correlation. Further,by applying the quantum renormalization-group (QRG) approach, M. Kargarian et al investigated the entanglement in the anisotropic Heisenberg model XXZ 1 (); XXZ 2 () and discussed the nonanalytic behaviors and the scaling close to the quantum critical point of the system. Recently We have calculated the block-block entanglement in the XY model without and with staggered Dzyaloshinskii-Moriya (DM) interaction by using this QRG method and have found the DM interaction can enhance the entanglement and influence the QPT of the system Fu-Wu Ma (); fu wu ma ().

Inspired by above idea, we apply the real-space renormalization-group (RSRG) approach which is developed in the Refs. 24 (); 25 (); 26 (); 27 (); 28 (); 29 (); 30 (); 31 (); 32 (); 33 (); 34 () to study the thermal entanglement between two end sites in the spin chain, Koch curve and on the diamond-type hierarchical (DH) lattices and analyze the influence of the temperature, the anisotropy parameter and the system size on the entanglement.

## Ii Model and method

The effective Hamiltonian of the spin- anisotropic ferromagnetic Heisenberg spin chain with sites is

(1) |

where denote the Pauli operators at site . , is the exchange coupling parameter, and is the Boltzmann constant. For simplicity, we assume that and . The sum is over all the nearest-neighbor spin pairs and is the anisotropy parameter. For and , the isotropic Heisenberg (XXX) and Ising model are obtained, respectively. The state of the above system in thermal equilibrium can be described by the density operator , where Tr is the partition function.

The entanglement of two-qubit system in the thermal state can be calculated by the negativity Negativity () which is based on the partial transpose method Partial transpose (). The negativity is defined as

(2) |

where is the negative eigenvalue of , and denotes the partial transpose with respect to the first subsystem. According this definition one can easily obtain thermal entanglement of the two-spin chain with the Hamiltonian (such as Fig 1. (a) shown). However, when the size of system becomes large, the density matrix is difficult or impossible to gain. The entanglement of two terminal spins on this system can not be directly worked out.

We apply the decimation RSRG method to solve the above problem. This decimation RSRG method 26 (); 28 (); 37 () has proved to be successful in spin chain and especially the fractal lattices. For the spin chain, this decimation procedure is illustrated in Fig. 1 (a). Simply, the generator is taken out in the infinite system. The generator with the Hamiltonian is renormalized into the new two-site chain with the Hamiltonian by integrating the internal site with the partition function being preserved. This transformation can be described as

(3) |

We can obtain the recurrence relation between the original parameters and the new parameters by solving the trace Tr with the method developed in Refs 28 (); 30 (). Combining this relation and the negativity of two-site system, the entanglement between two terminal spins in three-qubit chain can be obtained as follow

(4) |

The entanglement between two distant terminal spins in the Heisenberg chain can be calculated after many iterations of the recurrence relation. We take use of the same method to study thermal entanglement between two terminal sites in Koch curve and on the DH lattices as shown in Fig. 1 (b) (c) and (d). The analytical expression about the negativity in Eq. (4) is difficult to obtain, we will show some numerical results.

## Iii Heisenberg spin chain

We first study how the entanglement between terminal sites in spin chain with different number of sites varies with temperature at (shown in Fig. 2). For different cases of , the results have the similar feature that the entanglement decreases monotonically with increasing temperature and it vanishes beyond the critical temperature . At , the system is in the entangled ground state. As temperature increases, the entanglement decreases due to the mixture of the unentangled excited state with the ground state. At , the system is governed by the unentangled excited state completely, and therefore the entanglement vanishes. Comparing the entanglement of terminal sites in the different system, it is found that the entanglement of the ground state decreases sharply with increasing . Different from the maximally entangled Bell ground state for system, the ground state for system becomes a degenerate and related to state which cause the decrease of entanglement. The energy gap between the ground and the unentangled excited state increases with little range but the thermal energy of system increases with large range when increases. For , the system can easily overcome the gap and enter the unentangled state. That leads to the decrease of . This phenomenon reflects that thermal fluctuation of internal sites may suppress quantum effect.

The influence of the anisotropic parameter on the entanglement between two terminal sites at a fixed temperature is plotted in Fig. 3. As can be seen, for case, the system firstly keeps in the ground state with the entanglement . Then, as approaches the isotropic point , the energy gap between the ground and the unentangled excited states becomes so small that the system can jump to the unentangled excited state. Therefore there exists a sharp decease for entanglement and the entanglement vanishes when is close to zero. This result also accords with that the entanglement can not exist in an isotropic Heisenberg ferromagnetic chain in Ref Nielsen (). For cases, it is found that the entanglement increases firstly with increasing because the ground state is related to and will change with . All entanglement jump down to zero when reaches to zero. From above results, we can see that the entanglement is fragile and when the entanglement does not exist whatever the temperature and the anisotropic parameter are.

## Iv Fractal lattices

The properties of phase transition on different fractal lattices have been studied by the RSRG method, and the entanglement on these self-similar lattices remains to be explored. We turn to the study of the entanglement between end sites in Koch curves with non-integer fractal dimension and plot the numerical results of negativity versus and for different in Fig. 4. Compared with the entanglement in the spin chain (), it has similar properties that the entanglement decreases with and the maximal of entanglement and are approximately equal. But the entanglement variation versus is very different from that in spin chain, i.e., the range that the entanglement can exist is smaller, the maximal of entanglement is lager. The entanglement decreases quickly as becomes large and there exists no entanglement any longer when . We can deduce this result from the similar Hamiltonian and open boundary conditions of these two systems.

Now we consider two kinds of DH lattices with fractal dimensions (lattice A, for simplicity) and (lattice B, for simplicity). The RG transformation on these DH lattices respectively have been shown in Fig. 1 (c) and Fig. 1 (d). We first discuss the dependence of the entanglement between terminal sites on with In Fig. 5, one can find some similar behaviors that the entanglement is the maximal value at , and the entanglement decreases with increasing temperature and vanishes beyond the critical temperature. However, some different phenomenons are also observed that the entanglement between end sites on these two lattices decrease more slowly with increasing and it still exists even though the system size becomes very large (). For the case of on lattice B, as increases, the entanglement at zero temperature decreases but the corresponding increases. It is obvious that the entanglement for different crosses at . This result indicates that these two fractal lattices have special energy level structure. The energy gap between the entangled ground state and the unentangled excited state is so large that the system can jump to the unentangled excited state only at higher .

The variation of the entanglement between end sites on DH lattices versus at is also discussed. As can be seen in Fig. 6 (a) for the lattice A, the entanglement firstly increases as increases, and then it quickly decays to zero when reaches the isotropic point. The entanglement decrease very slowly when becomes very large and there exist no cross point when reaches zero. For the case of lattice B, Fig. 6 (b) shows that the entanglement also exhibits stable and it changes very little when is not very close to zero. The entanglement mainly remains robust with the increase of . In this graph, we also observe that a ”entanglement crossing” occurs at . At a fixed temperature, the thermal excited energy of the system is determined. Only when is very close to zero, the energy gap between the entangled ground state and the unentangled excited state can become so small that the system can enters the unentangled state. It also indicates that the different fractal structure can influence the entanglement by changing the energy level structure of system.

## V Conclusion

We have investigated thermal entanglement between two end spins in Heisenberg chain, Koch curve and on two kinds of DH lattices with and by the decimation RG method. The effect of the temperature and the anisotropy parameter on thermal entanglement is discussed. It is found that the symmetry of system and the thermal fluctuation can suppress or promote the quantum effect at different conditions. We also have noticed that the entanglement on some special lattices may exhibit different property when the system size becomes large. The entanglement on two kinds of DH lattices is quite robust and it can survive even though becomes very large in contrast to that in spin chain. The phenomenon of the ”entanglement crossing” indicates that the special fractal structure does influence on the entanglement.

###### Acknowledgements.

This work is supported by the National Natural Science foundation of China under Grant No. 10775088, the Shandong Natural Science foundation under Grant No. Y2006A05, and the Science foundation of Qufu Normal University. The authors would like to thank Hong-Bing Li, Yin-Fang Li, and Cong-Fei Du for fruitful discussions and useful comments.Figure captions:

Fig. 1. The procedure of the RG transformation. From (a) to (d), it shows the transformation of one-dimensional spin chain, Koch curve, the diamond-type hierarchical lattice with fractal dimension and .

Fig. 2. The entanglement between two end sites in Heisenberg ferromagnetic chain versus temperature at for different number of sites (from top to bottom, , , and ).

Fig. 3. The entanglement between two end sites in Heisenberg ferromagnetic chain as anisotropy parameter at for different number of sites (from top to bottom, , , and ).

Fig. 4. The variation of the entanglement between end sites in Koch curve: (a) the entanglement versus temperature at . (b) the entanglement versus at .

Fig. 5. Upper panel: the negativity of two end sites on the DH lattice with versus at . Lower panel: the negativity of two end sites on the DH lattice with versus at . The negativity for different value of has a cross point at .

Fig. 6. Upper panel: the evolution of the entanglement between terminal sites on DH lattice with as increases at for different value of Lower panel: the entanglement between terminal sites on DH lattice with versus at for different value of . The entanglement for different has a cross point at .

### Footnotes

- thanks: Corresponding author

### References

- J. S. Bell, Physics (Lang Island city), 1, 195 (1964).
- Artur K. Ekert, Phys. Rev. Lett. 67, 661 (1991).
- C.H. Bennett and S.J. Wiesner, Phys. Rev. Lett. 69, 2881 (1992).
- C.H. Bennett et al, Phys. Rev. Lett. 70, 1895 (1993).
- B. E. Kane, Nature (Londond) 393, 133 (1998).
- C. H. Bennett and D. P. Divincenzo, Nature (Londond) 404, 247 (2000d).
- M. A. Nielson and I. L. Chuang, Quantum Computation and Quantum Information (Cambridge University Press, Cambridge, U.K., 2000).
- W. K. Wootters, Phys. Rev. A 63, 052302 (2001).
- M. C. Arnesen, S. Bose, V. Vedral, and S. Bose, Phys. Rev. Lett. 87, 017901 (2001).
- X. Wang, Phys. Rev. A 64, 012313 (2001).
- X. Wang, Phys. Rev. A 66, 034302 (2002).
- Indrani Bose and Amit Tribedi, Phys. Rev. A 72, 022314 (2005).
- A. M. Souza, M. S. Reis, D. O. Soares-Pinto, I. S. Oliveira and R. S. Sarthour, Phys. Rev. B 77, 104402 (2008).
- K. Audenaert, J. Eisert, and M. B. Plenio R. F. Werner Phys. Rev. A 66, 042327 (2002).
- L. Zhou, H. S. Song, Y. Q. Guo, and C. Li, Phys. Rev. A 68, 024301 (2003).
- Guo-Feng Zhang, and Shu-Shen Li, Phys. Rev. A 72, 034302 (2005).
- Zhao-Yu Sun, Kai-Lun Yao, Wei Yao, De-Hua Zhang, and Zu-Li Liu, Phys. Rev. B 77, 014416 (2008).
- Fardin Kheirandish, S. Javad Akhtarshenas, and Hamidreza Mohammadi Phys. Rev. A 77, 042309 (2008).
- A. M. Souza, M. S. Reis, D. O. Soares-Pinto, I. S. Oliveira, and R. S. Sarthour, Phys. Rev. B 77, 104402 (2008).
- A. Osterloh, L. Amico, G. Falci, and R. Fazio, Nature 416, 608 (2002).
- M. Kargarian, R. Jafari, and A. Langari, Phys. Rev. A 76, 060304 (2007).
- M. Kargarian, R. Jafari and A. Langari, Phys. Rev. A 77, 032346 (2008)
- Fu-Wu Ma, Sheng-Xin Liu and Xiang-Mu Kong, Phys. Rev. A 83, 062309 (2011)
- Fu-Wu Ma, Sheng-Xin Liu and Xiang-Mu Kong, Phys. Rev. A 84, 042302 (2011).
- A. P. Young and R. B. Stinchcombe, J. Phys. C: Solid State Phys 9, 4419 (1976).
- R. Jullien, J. Fields and S. Doniach, Phys. Rev. Lett. 38, 1500 (1977).
- M. Suzuki and H. Takano, Phys. Lett. A 69, 426 (1979).
- H. Takano and M. Suzuki, J. Stat. Phys. 26, 635 (1981).
- Anibal O. Caride, Constantino Tsallis, and Susana I. Zanette, Phys. Rev. Lett. 51, 145 (1983).
- A. M. Mariz, R. M. Zorzenon Dos Santos, C. Tsallis and R. R. Dos Santos, Phys. Lett. A 108, 95 (1985).
- A. M. Mariz, C. Tsallis and A. O. Caride, J. Phys. C 18 4198 (1985).
- S. R. White and R. M. Noack, Phys. Rev. Lett. 68, 3487 (1992).
- J. Ricardo de Sousa, Phys. Lett. A 216, 321 (1996).
- N. S. Branco, J. Ricardo de Sousa, Phys. Rev. B 62, 5742 (2000).
- J. Ricardo de Sousa, N. S. Branco, B. Boechat and Claudette Cordeiro, Physica A 328, 167 (2003).
- G. Vidal, and R. F. Werner, Phys. Rev. A 65, 032314 (2002).
- A. Peres, Phys. Rev. Lett. 77, 1413 (1996); M. Horodecki, P. Horodecki and R. Horodecki, Phys. Lett. A 223, 1 (1996).
- Andre M. C. de Souza, Phys. Rev. B 48, 3744 (1993).
- M. A. Nielsen, Ph. D. thesis, University of Mexico, 1198, e-print arXiv: quant-ph/0011036.