Entanglement in a fermionic spin chain containing a single mobile boson under decoherence

Entanglement in a fermionic spin chain containing a single mobile boson under decoherence

Hamid Arian Zad Young Researchers and Elite Club, Mashhad Branch, Islamic Azad University, Mashhad, Iran e-mail: arianzad.hamid@mshdiau.ac.ir    Majid Moradi (MKMajid) Physics Department, Shahrood University of Technology, Shahrood, Iran e-mail: majidofficial@gmail.com
Abstract

The concurrence between first and the last sites of a fermionic spin chain containing a single boson is rigorously investigated at finite low temperature in the vicinity of a weak homogeneous magnetic field. We consider the boson as a mobile spin-1 particle through the chain and study concurrence without/under decoherence and express some interesting phase flip and bit flip reactions of the pairwise entanglement between first and the last half-spins in the chain. Our investigations show that the concurrence between two considered half-spins has different behavior for various positions of the single boson along the chain. Indeed, we realize that the single boson mobility has an essential role to probe the pairwise entanglement intensity between two spins located at the opposite ends of a fermionic chain. Interestingly, the entanglement remains alive for higher temperatures when the boson is the nearest neighbor of the first fermion. When the boson is at the middle of chain, it is demonstrated that the threshold temperature (at which the concurrence vanishes) versus decoherence rate can be considered as a threshold temperature boundary. These results pave the way to set and interpret the numerical and analytical expressions for utilizing quantum information in realistic scenarios such as quantum state transmission, quantum communication science and quantum information processing, where the both fermion-fermion and fermion-boson correlations should be taken in to account.

1 Introduction

Up to now, there was a tremendous interest for investigating entanglement of various Heisenberg models ABGARYAN20155 (); Sengupta2009 (); Horodecki2009 (); Ananikian2012 (); Wootters1998 (); Langari2008 (); TWelang ().There have been a great number of studies focusing on the applications of entanglement in quantum teleportation, spintronics and quantum information processing Mcmahon2007 (); Barmettler2010 (); Dolde2013 (); Furman2012 (). Spin models are known as simple popular quantum systems in solid state physics and quantum comuting which exhibits remarkable entanglement features Ollivier2001 (); Canosa2004 (); Zhang2005 (); GHAHARI2007 (); ZHANG2011 (), specially quantum antiferromagnetic spin chains have been the subject of numerous experimental as well as theoretical studies for a long period of time Tonegawa1998 (). A wide variety of the analytical and numerical techniques have been developed in order to study Ising-Heisenberg spin models ABGARYAN20155 (); Torrico2014 (); Chuang2008 (); Zad2017 (); Rojas2012 (); ABGARYAN201510 () due to enhance the applications of quantum information theory in solid state physics and condensed matter such as investigating thermal entanglement for the mixed spin chains Tonegawa1998 (); Durga2013 (); Movahedian2017 (); Arian2018 (); Arian2018S (); Zhe2009 ().

Motivated by the real materials it has already been proposed the class of one-dimensional plaquette chains to investigate the magnetic properties as well as quantum state transferring from theoretical and analytical points of view Miyahara1999 (); Souza2016 (). Besides the frustrated spin- orthogonal-dimer chain, Koga et al. investigated in detail the first-order quantum phase transition of the orthogonal-dimer spin chain Koga2000 (); Koga2002 (). A wide variety of the orthogonal dimer plaquette investigations have been carried out from theoretical and experimental view points. Driven by the comments given in Ref. Paulinelli2013 (), we here consider the orthogonal dimer plaquette chain with antiferromagnetic coupling as schematically described in Fig. 1 with an analytical Hamiltonian. On the other hand, the spin-impurity atoms were also widely discussed by some authors, for instance, T. Fukuhara et al. Fukuhara2013 () studied the quantum dynamics created by the spin-impurity atom, as it propagates in a one-dimensional spin chain. A specific probability distribution of the impurity at different times using single-site-resolved imaging of bosonic atoms have been probed as well. Their experimental results open new prospects to studying mobile impurities in quantum spin systems specially quantum fluids. Accordingly, we consider the dimer plaquette as a single mobile boson in a one-dimensional spin chain which dramatically affects on the thermal entanglement between a pair of fermions.

Heretofore, the dynamic behavior of spin chains under decoherence has been particularly studied Li2011 (); Maziero2009 (); Zad2016 (). Maziero et.al. investigated rigorously the dynamics of quantum and classical correlations under decoherence and carried out some interesting outcomes: classical correlation remains constant over time and quantum correlation decays monotonically over time Maziero2009 (). During last decade, similar tasks has been made to understand defect production due to the dynamics of a spin model through a critical point Sengupta2009 (); Zurek2015 (); Sengupta2008 ().

Because quantum communications allows to have unconditionally stable communication, and on the other hand, a quantum processor could dramatically speed up computations, there is a great deal of attentions in techniques to protect quantum information from noise Smith1812 (). Hence, in this paper we are going to study and compare the entanglement of three different types of antiferromagnetic fermionic chains with defect caused by a boson under decoherence. In this regard, we consider the entanglement between first and the last half-spins in a finite size fermionic chain consisting of a single mobile boson. Three types of the suggested model with different spin arrangement are determined by changing the position of the single boson along the chain.

The paper is organized as follows. In Sec. 2 we first introduce the three sub-models with analytical Hamiltonians and describe the concurrence as a measure of entanglement by characterizing pairwise reduced density matrix. In Sec. 3 we show the numerical calculations and simulations of the pairwise entanglement under decoherence. Section 4 is devoted to conclusions.

2 The Model

In the present work, we assume an antiferromagnetic spin-1/2 chain with the total number of sites containing a single boson and four fermions with three different types of the spin arrangement, namely, we consider various locations for the boson along the chain and create three dissimilar sub-models. We witness that different positions of the boson along the chain introduce different properties for the entanglement between two fermions located at the opposite ends of the chain. To study the behavior of such systems, their Hamiltonians should be introduced. Firstly, we start our investigations by contemplating the boson at the middle of chain. So, the Hamiltonian can be expressed as

(1)

where is the total number of particles (here ), denotes exchange interaction between half-spins, is the exchange interaction between the boson and its nearest fermions. and are anisotropy properties between, respectively, fermion-fermion and fermion-boson. is the homogeneous magnetic field applied along the z-axis. Here, are the Pauli operators and

(2)

After driving eigenvalues and eigenvectors of the Hamiltonian (1), in the standard basis, the total density matrix in equilibrium can be given by . Regarding this, we can extract reduced density matrix of the half-spins at two opposite ends of the chain as

(3)

where is thermal partition function of the total system. (we set ) and is the temperature.

In order to study thermal pairwise entanglement, we use the concurrence which can be defined by
, where and are square roots of the eigenvalues of the inner product with ( is the complex conjugate of matrix ).

3 decoherency

Recently, it has been revealed that the concurrence of bipartite subsystem (1/2,1/2) in a tripartite mixed spin chain (1/2,1,1/2) undergoing only the phase flip channel by means of the Kraus operators, is different in behavior Zad2016 (). In this work, we investigate the thermal entanglement for two half-spins at two opposite ends of the spin chain including a single boson under decoherence by means of the concurrence. Here, the state of bipartite system introduced by (3) is exactly investigated when it undergoes both phase flip and bit flip.

3.1 Phase flip channel.

Kraus operators for the phase flip channel are given by , , and that describe the noise channels and , where . The evolved state under the phase flip channel can be characterized by a completely positive trace preserving (CPTP) map
, where is the first site of the chain under consideration and is its last half-spin.

First we turn our attention to investigate the thermal entanglement between two spins located at the ends of the fermionic chain when the single boson is set at the middle of chain (Fig. 1) with Hamiltonian (1).

Figure 1: Schematic representation of the spin chain containing a single boson located in the middle with spin arrangement (1/2,1/2,1,1/2,1/2)(sub-model 1).

For the mentioned spin formation, we investigate the concurrence at low temperature and weak magnetic field under decoherence for two half-spins at two ends of the chain and compare it with the case where the phase damping rate is zero (Fig. 2). By inspecting Fig. 2(a), one can see that the concurrence gradually decreases upon rising the temperature. As increases, the maximum amount of concurrence decreases at low temperature. Also, it is obvious that by increasing until , slope of the curve is reduced and the concurrence vanishing goes further on the temperature axis, this actually means that the entanglement can stay alive for higher temperatures. When increases more, the slope of the concurrence curve increases and the bipartite entanglement vanishes at lower temperatures.

Figure 2: The concurrence for two half-spins located at the opposite ends of the chain for various fixed values of . (a) When the single boson is fixed at the middle of chain, and (b) when the boson is fixed at the left hand side of the chain. denotes the concurrence under the phase flip. Hamiltonian Parameters have been taken as , , , and .

According to our investigations, the concurrence shows different behaviors with respect to the various positions of the boson. Namely, the single boson movement through the chain can alter the thermal entanglement property between two considered half-spins. For instance, suppose that the boson to be the nearest neighbor of first site of the chain (sub-model 2).

Figure 3: The threshold temperature against decoherence rate when the boson is situated at the middle of chain (blue solid line) in comparison with that of the boson is placed in the left hand side of the chain (red dashed-dot line) under phase flip.

The concurrence of this model is presented in Fig. 2(b), one can see that the maximum value of the concurrence has increased in comparison with Fig. 2(a) for all fixed values of . On the other hand, the concurrence vanishing occurs at higher or equal temperatures. We call these temperatures at which the concurrence becomes zero as threshold temperatures. The threshold temperature versus is illustrated in Fig. 3, revealing that the optimum point for which the entanglement survives the longest is about . We understood that the threshold temperature values experience a special sweep point under the considered circumstances. It also shows that changing the position of boson to left hand side of the chain (sub-model 2) causes the threshold temperature slightly increase (red dashed-dot line). As a matter of fact, the threshold temperature diagram when the boson is situated at the middle of chain (blue-solid line) can be considered as the threshold temperature boundary.

Figure 4: Top model shows the spin chain containing a single boson located in the left side with spin configuration (1/2,1,1/2,1/2,1/2) (sub-model 2). Bottom model shows the spin chain consisting of a single boson localized in the right side with spin configuration (1/2,1/2,1/2,1,1/2) (sub-model 3).

Furthermore, we examine the concurrence for the chain under consideration in which the single boson is the nearest neighbor of the last half-spin (sub-model 3 presented in Fig. 4). In this case, we realize that there is no quantum entanglement between two fermions located at the opposite ends of the chain under the considered circumstances, namely, the concurrence is zero. Moreover, we witness that the pairwise entanglement remains zero for a such spin configuration under the phase flip. Finally, we conclude that the position of the single boson in the chain plays an important role to describe the pairwise entanglement intensity between the first and the last half-spins.

3.2 Bit flip channel

Now, we investigate three different types of the suggested model under bit flip. The Kraus operators for the bit flip channel are , , and that describe another kind of noise channels between sites and . With similar calculations, we obtain the concurrence of state (3) under the bit flip as well as the phase flip case. For the case where the single boson is localized at the middle of chain, in Fig. 5 we display the temperature dependence of the thermal entanglement between two fermions at the opposite ends of the chain by means of the concurrence.

Figure 5: The concurrence of half-spins at opposite ends of the chain versus temperature for various fixed values of under bit flip. The boson is localized at the middle of chain for the same set of the coupling constants as used in Fig. 2.

The curves indicate the concurrence behavior with respect to the various fixed values of . We see that the thermal pairwise entanglement is very sensitive upon increasing the decoherence rate . When increases from zero, the maximum value of the concurrence decreases, until for the range the concurrence disappears at low temperature and weak magnetic field.

Figure 6: The concurrence of two half-spins at opposite ends of the chain when the boson is localized at the left side under bit flip (a) in the interval , and (b) for .

In Fig. 6 we display the concurrence as a function of temperature for various fixed values of for the case where the single boson is boosted to the left side of the chain. To more clarity, we investigate the behavior of concurrence versus two intervals of , i.e., and separately. The maximum amount of concurrence remarkably decreases upon increasing for the considered interval (Fig. 6 (a)), further increase of (Fig. 6 (b)) leads to sharply increasing the concurrence. Also it is obvious that for , the concurrence continues to exist for higher temperatures. So the entanglement vanishes at higher temperatures. To study this case more accurately, we plot the threshold temperature against decoherence rate in Fig. 7.

Figure 7: The threshold temperature against decoherence rate when the boson is the nearest neighbor of the first fermion under bit flip.

In the interval the maximum value of concurrence decreases monotonically until reaches a recursive point () at which the maximum value starts its increasing with further increase of . Moreover, in this interval of the decoherence rate, with increase of the threshold temperature tends to the lower values until reaches a recursive point at . By more increasing the decoherence rate (), the maximum amount of concurrence shifts toward higher values, on the other hand, the threshold temperature shifts toward higher temperatures as well. There is a steep increase of the temperature at . With further increase of the threshold temperature begins to decrease again. As a result, the threshold temperature is strongly dependent on the decoherence rate .

It is worthy to note that the concurrence is zero under both phase flip and bit flip for the case where the boson is localized in the right side of the chain at low temperature and weak magnetic field.

4 Conclusions

We considered three variants of the antiferromagnetic half-spin chains including a single boson in three different positions (nearest neighbor of the first site, middle of the chain and nearest neighbor of the last site) under decoherence at low temperature and weak magnetic field. We rigorously investigated the thermal entanglement between two fermions at the opposite ends of the chain under the both phase flip and bit flip by means of the concurrence. We realized that the localization of the boson plays an essential role in determining bipartite entanglement intensity in the chain, i.e., by moving the boson along the chain the maximum concurrence remarkably changes. We compared pairwise entanglement of two fermions at the opposite ends of the chain when the single boson is localized at the middle of the chain with the case where the boson is the nearest neighbor of the first fermion and concluded that the maximum value of the concurrence will increase at low temperature. Interestingly in this manner, the entanglement survives for higher temperatures. Furthermore, when the boson is nearest neighbor of the last fermion, the thermal entanglement vanishes under both bit flip and phase flip.

Also, the pairwise thermal entanglement is strongly dependent on the decoherence rate for the various spin arrangements under the phase flip. Indeed, the entanglement death occurs at higher temperatures upon increasing decoherence rate until a threshold point at which the entanglement death turns toward lower temperatures. We investigated the threshold temperature versus decoherence rate and revealed that the threshold temperature increases by increasing the decoherence rate. The maximum amount of the threshold temperature occurs at a specific decoherence rate. With further increasing decoherence rate, the value of the threshold temperature decreases. Here, we proved that the threshold temperature diagram versus the decoherence rate for the case when the boson is localized at the middle of the chain can be considered as the minimum boundary limit of the threshold temperature for all spin configurations of the chain.

Under the bit flip, we concluded that the concurrence decreases upon increasing both temperature and decoherence rate at weak magnetic field. The same as phase flip, the concurrence is dependent on the decohrence rate. We observed a monotonically decreasing behavior of the threshold temperature until a jumping point at which the threshold temperature suddenly increases to its maximum value. After that, the threshold temperature decreases with a gentle slope.

The method introduced in this paper can be straightforwardly generalized to larger quantum spin systems containing different spin arrangements and seeking new interesting outcomes. However, all results should be compatible with the remarks and outcomes expressed in this paper.

References

References

  • (1) V. S. Abgaryan and N. S. Ananikian and L. N. Ananikyan and V. Hovhannisyan, Solid State Communications 203, 0038-1098 (2015)
  • (2) K. Sengupta and D. Sen, Phys. Rev. A 80, 032304 (2009)
  • (3) R. Horodecki, P. Horodecki, M. Horodecki and K. Horodecki, Rev. Mod. Phys., 81, 865-942 (2009),
  • (4) N. S. Ananikian, L. N. Ananikyan, L. A . Chakhmakhchyan and O. Rojas, J Phys.: Cond. Matt. 24, 256001 (2012)
  • (5) W. K., Wootters, Phys. Rev. Lett. 80, 2245–2248 (1998); W. K., Wootters, Quantum Information and Computation 1, 27-44 (2001)
  • (6) R. Jafari, M. Kargarian, A. Langari, A. and M. Siahatgar, Phys. Rev. B 78, 214414 (2008)
  • (7) Werlang, T. and Trippe, C. and Ribeiro, G. A. P. and Rigolin, Gustavo, Phys. Rev. Lett. 105, 095702 (2010)
  • (8) D. McMahon, Quantum Computing Explained (John Wiley, New York, 2008)
  • (9) P. Barmettler, M. Punk, V. Gritsev, E. Demler and E. Altman, New J. Phys. 12, 055017 (2010)
  • (10) F. Dolde and I. Jakobi and B. Naydenov and N. Zhao and S. Pezzagna and C. Trautmann and J. Meijer and P. Neumann and F. Jelezko and J. Wrachtrup, Nat. Phys. 9, 139 (2013)
  • (11) G. B. Furman, V. M. Meerovich and V. L. Sokolovsky, Quant. Inf. Pro. 11, 1603 (2012)
  • (12) H. Ollivier and . Zurek, Phys. Rev. Lett. 88, 017901, (2001)
  • (13) N. Canosa and R. Rossignoli, Phys. Rev. A 69, 052306 (2004)
  • (14) G. F. Zhang and S. S. Li Phys. Rev. A 72, 034302 (2005)
  • (15) F. Ghahari and V. Karimipour and R. Shahrokhshahi Phys. Lett. A 363, 271 - 276 (2007)
  • (16) G.-F. Zhang, Z.-T. Jiang, A. Abliz, Ann. of Phys. 326, 867 - 875 (2011)
  • (17) T. Tonegawa, T. Hikihara, M. Kaburagi, T. Nishino, S. Miyashita and H. J. Mikeska, J. Phys. Soc. Jap. 67, 3 (1998)
  • (18) J. Torrico and M. Rojas and S. M. de Souza and O. Rojas and N. S. Ananikian, EPL 108, 50007 (2014)
  • (19) D. C. Li, X. P. Wang and Z. L. Cao, J. of Phys.: Cond. Matt. 20, 325229 (2008)
  • (20) H. Arian Zad and N. Ananikian, J. Phys. Condens. Matt. 29, 455402 (2017)
  • (21) O. Rojas, M. Rojas, N. S. Ananikian and S. M. D. Souza, Phys. Rev. A 86, 042330 (2012)
  • (22) V. S. Abgaryan, N. S. Ananikian, L. N. Ananikyan and V. Hovhannisyan, Solid State Comm. 224, 15 (2015)
  • (23) V. M. L. Durga Prasad Goli, S. Sahoo, S. Ramasesha and D. Sen, J. of Phys.: Cond. Matt. 25, 125603 (2013)
  • (24) H. Arian Zad and H. Movahhedian,Int. J. Mod. Phys. B 31,1750094 (2017)
  • (25) H. Arian Zad and N. Ananikian, J. Phys. Condens. Matt. 30, 165403 (2018)
  • (26) H. Arian Zad and N. Ananikian, Solid State Commun. 276, 24-27 (2018)
  • (27) Z. Sun, X. M.Lu and H. N. Xiong and J. Ma, New J. Phys. 11, 113005 (2009)
  • (28) S. Miyahara and K. Ueda, Phys. Rev. Lett. 82, 3701 (1999); H. Kageyama, K. Yoshimura, R. Stern, N. V. Mushnikov, K. Onizuka, M. Kato, K. Kosuge, C. P. Slichter, T. Goto, and Y. Ueda, Phys. Rev. Lett. 82, 3168 (1999)
  • (29) M. Rojas, S. M. de Souza and O. Rojas, Ann. of Phys. 377, 506 (2017)
  • (30) A. Koga, K. Okunishi, and N. Kawakami, Phys. Rev. B 62, 5558 (2000)
  • (31) A. Koga and N. Kawakami, Phys. Rev. B 65, 214415 (2002)
  • (32) H. G. Paulinelli, S. de Souza and Onofre Rojas, J. Phys.: Condens. Matter 25 306003 (2013)
  • (33) T. Fukuhara et al., Nat. Phys. 9, 235 (2013)
  • (34) B. Li, Z. X. Wang and S. M. Fei, Phys. Rev. A 83, 022321 (2011)
  • (35) J. Maziero, L. C. Celeri, R. M. Serra and V. Vedral, Phys. Rev. A 80, 044102 (2009)
  • (36) H. Arian Zad, Chinese Physics Letters 33, 090302 (2016)
  • (37) W. H. Zurek, U. Dorner and P. Zoller, Phys. Rev. Lett. 95, 105701 (2005)
  • (38) K. Sengupta, D. Sen and S. Mondal, Phys. Rev. Lett. 100, 077204 (2008)
  • (39) G. Smith and J. Yard, Science 321, 5897 (2008)
  • Comments 0
    Request Comment
    You are adding the first comment!
    How to quickly get a good reply:
    • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
    • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
    • Your comment should inspire ideas to flow and help the author improves the paper.

    The better we are at sharing our knowledge with each other, the faster we move forward.
    ""
    The feedback must be of minimum 40 characters and the title a minimum of 5 characters
       
    Add comment
    Cancel
    Loading ...
    214323
    This is a comment super asjknd jkasnjk adsnkj
    Upvote
    Downvote
    ""
    The feedback must be of minumum 40 characters
    The feedback must be of minumum 40 characters
    Submit
    Cancel

    You are asking your first question!
    How to quickly get a good answer:
    • Keep your question short and to the point
    • Check for grammar or spelling errors.
    • Phrase it like a question
    Test
    Test description