# Entanglement from density measurements: analytical density-functional for the entanglement of strongly correlated fermions

###### Abstract

We derive an analytical density functional for the single-site entanglement of the one-dimensional homogeneous Hubbard model, by means of an approximation to the linear entropy. We show that this very simple density functional reproduces quantitatively the exact results. We then use this functional as input for a local density approximation to the single-site entanglement of inhomogeneous systems. We illustrate the power of this approach in a harmonically confined system, which could simulate recent experiments with ultracold atoms in optical lattices as well as in a superlattice and in an impurity system. The impressive quantitative agreement with numerical calculations – which includes reproducing subtle signatures of the particle density stages – shows that our density-functional can provide entanglement calculations for actual experiments via density measurements. Next we use our functional to calculate the entanglement in disordered systems. We find that, at contrast with the expectation that disorder destroys the entanglement, there exist regimes for which the entanglement remains almost unaffected by the presence of disordered impurities.

###### pacs:

03.67.Mn, 31.15.es, 37.10.Jk, 71.10.Fd## I Introduction

Entanglement, one of the most intriguing features of quantum mechanics, is an essential ingredient in quantum information theory horodecki (), and it is also receiving increasing attention in the study of fundamental phenomena such as quantum phase transitions (QPT) rev.solids (). Entanglement has been investigated in many different physical systems, such as photons in optical cavities davi (), ultracold atoms interacting with light vvf2007 (); porras (), and in solid state systems rev.solids (), where it has been associated to particle-particle correlations plenio ().

Recently, special attention has been devoted to solids, since they are considered one of the most promising systems for the development of quantum information devices stoneham (). In particular, the Hubbard model hubb (), which captures most important physics of interacting many-body systems tasaki (); paiva (); machida (); derz (), has been intensively used to understand entanglement in solids gu1 (); larsson1 (); vvf2006 (); vvf2008 (); carr (); coe2010 (); epl2010 (). In this context, creating, distributing, quantifying, and controlling entanglement represent great challenges. Much progress has been achieved and some well defined entanglement measures, such as entropies and concurrence, have been proposed rev.solids (). The von Neumann entropy has been used for measuring the bipartite entanglement of pure states rev.solids (). For systems with very many degrees of freedom, the linear entropy – which gives an indication of the number and spread of terms in the Schmidt decomposition of the state – has been proposed as a viable alternative buscemi2007 (); irene2008 (); shirwan2009 ().

For the homogeneous Hubbard model, the single-site entanglement zanardi () and its properties have been investigated vvf2006 () and associated to QPT’s oster (); gu1 (); larsson1 (). A local-density approximation (LDA) for the entanglement entropy of any generic inhomogeneous system has been proposed vvf2008 () and applied to the Hubbard model for different types of inhomogeneities. Also block-block entanglement has been investigated in connection to QPT gu2 () and from the view point of its universal and non-universal contributions vvf2008pra ().

Another extremely important – but less addressed – question concerning entanglement quantifiers, is whether and how theoretical predictions can be connected to quantities that are measured experimentally. It has been shown that magnetic susceptibility can be an entanglement witness gosh (); vedral (); almeida (), which is used to point out the presence of entanglement. Also momentum-momentum correlation lukin (); flo (), energies irene2008 (); shirwan2009 (), and probability density at highly symmetric points shirwan2009 () have been proposed as entanglement indicators, which reproduce the general trend of entanglement. However, a quantity able to predict entanglement quantitatively and, in addition, with the potential of being experimentally measurable, is still missing. Here we show that a very promising candidate for this is the particle density.

As a direct consequence of the Hohenberg-Kohn theorem hk (), the pillar of Density Functional Theory (DFT), ground-state entanglement must be a functional of the density. For the one-dimensional homogeneous Hubbard model this functional may be numerically accessed ba (), but an explicit expression cannot be obtained due to the implicit dependence of the ground-state energy on the density. Different analytical and numerical techniques have been proposed for this calculation gu1 (); larsson1 (); vvf2006 (). Nevertheless, obtaining numerically the ground-state energy has been so far the main limiting factor for a complete analytical understanding of the entanglement in the homogeneous Hubbard model. For inhomogeneous systems the scenario is even more complicated, despite the application of powerful DFT methods vvf2008 ().

From the experimental point of view, it has been recently shown that the measured density in systems of ultracold atoms in optical lattices with single-site resolution greiner3 (); bloch2 () can be very accurately reproduced by the single-site density of the Hubbard model. Several theoretically predicted phenomena have already been verified in such experiments, e.g., the superfluid to Mott-insulator transition greiner2 () and BEC-BCS crossovers bart (); via DFT, the density may now become a powerful tool for measuring entanglement in experiments.

Motivated by this, in the present work we provide an analytical expression for the density-functional for the entanglement of homogeneous strongly correlated systems. Our approach is based on an approximation of the linear entropy for the one-dimensional Hubbard model valid, in principle, for a specific range of interactions and densities. The functional is then extended to a vast range of parameters by means of well established transformations. We test the performance of our density functional in homogeneous and, via an LDA, in three different inhomogeneous systems. We find that the entanglement in many-body systems can be quantitatively obtained directly from our very simple density functional. Finally we use our method for investigating the entanglement in disordered systems. To achieve a good statistics, this analysis requires many realizations and, therefore, could not be explored so far, while our functional greatly reduces its computational cost. We find that in general disorder destroys entanglement. Surprisingly though, there exist regimes for which the degree of entanglement remains almost unaffected.

## Ii Analytic density functional for the homogeneous system

We start from the homogeneous Hubbard model,

(1) |

in which all sites are equivalent. Here , the on-site interaction, and , the external potential, are rescaled by the hopping parameter between neighboring sites; and are fermionic creation and annihilation operators, and is the number operator at site with spin .

The single-site entanglement is quantified using the linear entropy

(2) |

where the reduced density matrix of site , , is the trace of the total density matrix with respect to the remaining sites. For the homogeneous Hubbard model , with the single spin occupation probability, the double occupation, and the zero occupation probabilities zanardi ().

This allows us to express the single-site entanglement of the homogeneous Hubbard model as an explicit function of the occupation probabilities

(3) |

For non-magnetic systems the probabilities are given by

(4) | |||||

(5) | |||||

(6) |

where is the per-site ground-state energy.

In order to express the entanglement as an explicit function of the density, we will approximate Eq.(4). Initially, we focus on repulsive interaction and , but as shown later, the results will be extended to and to by means of particle-hole transformations. An exact analytical expression for the energy, and consequently for , is unknown in this regime, except at , for which the exact Lieb-Wu solution ba () gives

(7) |

where are Bessel functions of order .

For low-density systems the double occupancy is extremely reduced due to the repulsive interactions w2 (); for higher densities, increases and can not be discarded. We then consider when , and assume Eq.(7) as an approximation for otherwise foot1 ().

Then the linear entropy can be written as the explicit density functional

where is a step function with for and for . We extend our functional to the regime of using the particle-hole symmetry, by replacing with in Eq.(II). Such extension is necessary to describe a larger variety of systems, and also in view of extending the results via DFT methods to more realistic inhomogeneous systems.

For attractive interactions the double occupancy can neither be discarded nor approximated by the Lieb-Wu result. The extension of our analysis to can be though obtained by means of the mapping between the energies of repulsive and attractive systems campo (),

(9) |

which is valid for any density. From (4) and (9) we obtain

(10) |

and the density-functional for the entanglement becomes

(11) | |||||

In Figure 1 we compare our density-functionals for the entanglement, Eq.(II) and Eq.(11), to the exact entanglement Eq.(3). The agreement is very good, considering the simplicity of our density-functional expressions. For the exact results are very accurately reproduced because the double occupation is fully taken into account in Eq.(11). For the results reveal that the relevant dependency on for the entanglement is taken into account by the other occupation probabilities, i.e., the main contribution of the double occupancy to the entanglement is almost independent on the density, therefore the simplified form of adopted in Eq.(II) is enough for predicting the entanglement with good precision. In particular Figure 1 shows that there are ranges of densities, and , for which many-body interactions do not play a relevant role: in a homogeneous system, up to (and for when considering particle-hole transformation), the double occupation is in fact strongly reduced as here (). In turn in the Hubbard model, where only on-site interactions are considered, this implies a lesser role of the particle-particle interactions note_many_body ().

However there is a small range of combined parameters – namely for (see inset of Fig. 1) – for which our approximation to is not as good: here also because on-site repulsion is small, may become significant. The maximum error in this specific region decays exponentially from at to at . Considering that typical interactions in solids correspond to u.solids () and in experiments with optical lattices the interaction can be easily tuned in a wide range (up to ) u.opt (), for most systems of practical interest our approach can provide the entanglement with a precision better than for any density.

The very good agreement of our density functionals with exact calculations shows that the entanglement of actual experiments described by the Hubbard model greiner3 (); bloch2 () could be obtained accurately and directly from density measurements. Also our results imply much simpler calculations which opens possibilities for the entanglement investigation in more complex systems, such as the disordered media we will explore in the last section.

## Iii Density Functional for Inhomogeneous Systems

In order to investigate the performance of our density functional in inhomogeneous systems, we present in what follows numerical calculations for the Hubbard model with an inhomogeneous external potential . We adapt to this purpose the LDA formalism developed in Ref.vvf2008 (). The average single-site linear entropy of inhomogeneous systems is then given by

(12) |

In Eq.(12) the inputs are our density-functionals, Eqs.(II) or (11), calculated at the local -site density , and averaged over the sites foot2 ().

We start with the inhomogeneous harmonic potential, which may model, for example, the trap in systems of ultracold atoms in optical lattices greiner3 () or quantum dots quantumdot (); coe2010 (); epl2010 (). This inhomogeneous system contains very interesting physics: it undergoes subtle transitions in the density profile as the potential increases opticaltrap () and the single-site entanglement is able to point out each and every one of these phases vvf2008 (). We have chosen , which represents a typical density reached in the experiments morsch (), and , as experimentally the interaction can be tuned up to u.opt (). In Figure 2–a we compare our density-functional for the entanglement to accurate numerical results. The agreement is indeed very good and all the stages of the particle density vvf2008 () are correctly pointed out by the derivative of our functional (see the entanglement derivative in Fig.2–a).

We have also checked the performance of our approach for different inhomogeneities and different values of , such as a superlattice defined by a periodic modulation in the external potential (), and the presence of a single impurity with intensity () as shown in Fig. 2–b and Fig. 2–c, respectively. The agreement is quantitatively similar to that found for the harmonic confinement: for the superlattice the deviations are smaller than , while for the impurity system we found deviations of about for all impurity intensities (note the scale in the y-axis).

Our results for three different physical systems confirms that the entanglement can be calculated directly from density measurements by using our density-functional even for inhomogeneous systems. As all the above tests have confirmed the reliability of our functional, we next use our density-functional for investigating disordered systems.

## Iv Entanglement in Disordered Systems

The investigation of any physical property in disordered systems usually requires an average over a large number of samples in order to improve statistics and avoid results peculiar to specific realizations. Thus, for example, for the generation of a single point in a system of 100 sites considering 100 samples, it would be necessary to perform calculations of the desired average single-site property. This represents a serious limiting factor for estimating the entanglement via methods involving the computationally-expensive ground-state energy calculation. In this sense our density-functional prescription represents a great advantage.

Here we consider the average single-site entanglement over 100 realizations of a disordered system defined by a concentration of randomly distributed impurities with intensity . In Figure 3 we present the entanglement obtained with our density functional as a function of . The entanglement is maximum for non disordered chains () and for a plateau is observed for . This plateau in the entanglement reflects the competition between the repulsive interaction and the impurities’ attractive potential at the Mott-insulator transition induced by at the impurities’ sites. A plateau is also observed in systems with a single impurity (see Fig. 2–c), where the percentage drop of entanglement with the impurity strength is however much smaller. Interestingly, for , we also find a plateau for (see Fig..3, upper panel). The physics behind both plateaus is similar, but for repulsive and one needs a higher impurity concentration to reach the Mott-insulator transition, as in this case it is induced at the sites with no impurities. The reversed role between sites with and without impurities manifests itself in the reflection symmetry with respect to for systems with and , as shown in Fig.3 – lower panel. The width and height of the plateaus increase towards smaller impurity intensities with increasing (decreasing) impurity concentrations for (); the plateaus’ width also increases with as , as a result of the competition between the two potentials (see inset of Fig.3).

In disordered chains, similarly to what has been observed for other inhomogeneous systems, inhomogeneities are in general destructive to the entanglement, and the higher and , the more the entanglement is destroyed. However our calculations show that there are regimes for which the degree of entanglement is almost unaffected. The entanglement in this region is almost as higher as the homogeneous case, . This property could be useful for indicating the regime of parameters for which real-life devices may present a high degree of entanglement despite the presence of impurities.

## V Conclusions

In conclusion we have derived an analytical density functional for the entanglement in the Hubbard model represented by a very simple mathematical expression. We have shown that it reproduces the entanglement quantitatively in homogeneous and inhomogeneous systems. We tested our functional on very different inhomogeneous systems, i.e. an harmonic confinement, which simulates experiments of trapped atoms in optical lattices, as well as on a superlattice and a single-impurity system.

We applied our density-functional prescription to the investigation of the entanglement in disordered systems: this analysis would be computationally very expensive – and as such difficult to achieve – with other methods. We found that, although in general the disorder destroys the entanglement, there exist a regime of parameters, both for repulsive and attractive impurity potential, for which the degree of entanglement remains almost unaffected. This result is promising as it suggests potential robustness for future device production.

In addition of allowing simpler numerical calculations and paving the road to analytical analysis of many-body and quantum information properties, our findings are of direct interest for measuring the entanglement from experimental setups.

## Vi Acknowledgements

We thank E. A. L. Henn for fruitful discussions. VVF acknowledges funding from Brazilian agency CAPES(4101/09-0) and IDA from EPSRC grant EP/F016719/1.

## References

- (1) R. Horodecki, P. Horodecki, M. Horodecki, and K. Horodecki, Rev. Mod. Phys. 81, 865 (2009).
- (2) L. Amico, R. Fazio , A. Osterloh, and V. Vedral, Rev. Mod. Phys. 80, 517 (2008).
- (3) S. P. Walborn, P. H. Souto Ribeiro, L. Davidovich, F. Mintert, and A. Buchleitner, Phys. Rev. A 75, 032338 (2007).
- (4) V. V. França and G. A. Prataviera, Phys. Rev. A 75, 043604 (2007).
- (5) D. Porras and J. I. Cirac, Phys. Rev. A 78, 053816 (2008).
- (6) J. Eisert, M. Cramer, and M. Plenio, Rev. Mod. Phys. 82, 277 (2010).
- (7) M. Stoneham, Physics 2, 34 (2009).
- (8) J. Hubbard, Proc. R. Soc London A 276, 238 (1963).
- (9) H. Tasaki, J. Phys.: Condens. Matter 10 4353 (1998).
- (10) T. Paiva and R. R. dos Santos, Phys. Rev. B 58, 9607 (1998).
- (11) M. Machida, S. Yamada, Y. Ohashi, and H. Matsumoto, Physica C 445-448, 90 (2006).
- (12) O. Derzhko, J. Richter, A. Honecker, M. Maksymenko, and R. Moessner, Phys. Rev. B 81, 014421 (2010).
- (13) S.-J. Gu, S.-S. Deng, Y.-Q. Li, H.-Q. Lin, Phys. Rev. Lett. 93, 086402 (2004). S.-J. Gu, Chin. Phys. Lett. 26, 026401 (2008).
- (14) D. Larsson and H. Johannesson, Phys. Rev. Lett. 95, 196406 (2005).
- (15) V. V. França and K. Capelle, Phys. Rev. A 74, 042325 (2006).
- (16) V. V. França and K. Capelle, Phys. Rev. Lett. 100, 070403 (2008).
- (17) L. D. Carr, M. L. Wall, D. G. Schirmer, R. C. Brown, J. E. Williams, and C. W. Clark, Phys. Rev. A 81, 013613 (2010).
- (18) J. P. Coe, V. V. França and I. D’Amico, Phys. Rev. A 81, 052321 (2010).
- (19) J. P. Coe, V. V. França and I. D’Amico, arXiv:1007.5441 (2010), to appear in Europhys. Lett.
- (20) F. Buscemi, P. Bordone, and A. Bertoni, Phys. Rev. A 75, 032301 (2007).
- (21) J. P. Coe, A. Sudbery, and I. D’Amico, Phys. Rev. B 77, 205122 (2008).
- (22) S. Abdullah, J. P. Coe, I. D’Amico, Phys. Rev. B 80, 235302 (2009).
- (23) P. Zanardi, Phys. Rev. A 65, 042101 (2002).
- (24) A. Osterloh, L. Amico, G. Falci, and R. Fazio, Nature 416, 608 (2002).
- (25) S. S. Deng, S. J. Gu, and H. Q. Lin, Phys. Rev. B 74, 045103 (2006).
- (26) V. V. França and K. Capelle, Phys. Rev. A 77, 062324 (2008).
- (27) S. Ghosh, T.F. Rosenbaum, G. Aeppli. S.N. Coppersmith, Nature 425, 48 (2003).
- (28) V. Vedral, Nature 425, 28 (2003).
- (29) A. M. C. Souza and F. A. G. Almeida, Phys. Rev. A 79, 052337 (2009).
- (30) E. Altman, E. Demler, and M. D. Lukin, Phys. Rev. A 70, 013603 (2004).
- (31) F. Mintert, A. M. Rey, I. I. Satija, and C. W. Clark, Europhys. Lett. 86, 17003 (2009).
- (32) P. Hohenberg and W. Kohn, Phys. Rev. 136, B864 (1964).
- (33) E. H. Lieb and F. Y. Wu, Phys. Rev. Lett. 20, 1445 (1968).
- (34) This is in contrast to the common belief that the strong interaction among particles would destroy local entanglement easily, one of the reasons for which the investigation of single-site, or more generally, local entanglement in many-body systems was somewhat left behind.
- (35) M. Imada, A. Fujimori, and Y. Tokura, Rev. Mod. Phys. 70, 1039 (1998).
- (36) R. Jördens, N. Strohmater, K. Günter, H. Moritz, and T. Esslinger, Nature 455, 204 (2008).
- (37) W. S. Bakr, J. I. Gillen, A. Peng, S. Fölling, and M. Greiner, Nature 462, 74 (2009).
- (38) J. F. Sherson, C. Weitenberg, M. Endres, M. Cheneau, I. Bloch, S. Kuhr, Nature 467, 68 (2010).
- (39) M. Greiner, O. Mandel, T. Esslinger, T. W. Hänsch, and I. Bloch, Nature 415, 39 (2002).
- (40) M. Bartenstein, A. Altmeyer, S. Riedl, S. Jochim, C. Chin, J. H. Denschlag, and R. Grimm, Phys. Rev. Lett. 92, 203201 (2004).
- (41) R. Jördens et. al, Phys. Rev. Lett. 104, 180401 (2010).
- (42) This procedure ensures the continuity of our approximated functional.
- (43) For all calculations open boundary conditions were used and the densities were obtained via DFT with LDA for the exchange-correlation energy.
- (44) V. L. Campo, Jr. and K. Capelle, Phys. Rev. A 72, 061602(R) (2005).
- (45) M. Rontani, F. Rossi, F. Manghi, and E. Molinari, Phys. Rev. B 59, 10165 (1999).
- (46) M. Rigol, A. Muramatsu, G. G. Batrouni, and R. T. Scalettar, Phys. Rev. Lett. 91, 130403 (2003).
- (47) O. Morsch and M. Oberthaler, Rev. Mod. Phys. 78, 179 (2006).