Advances in Condensed Matter Physics
Volume 2015 (2015), Article ID 397630, 4 pages
http://dx.doi.org/10.1155/2015/397630

Research Article

Entanglement Area Law in Disordered Free Fermion Anderson Model in One, Two, and Three Dimensions

Mohammad Pouranvari, Yuhui Zhang, and Kun Yang

National High Magnetic Field Laboratory and Department of Physics, Florida State University, Tallahassee, FL 32306, USA

Received 12 October 2014; Revised 11 February 2015; Accepted 12 February 2015

Academic Editor: Jan A. Jung

Copyright © 2015 Mohammad Pouranvari et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

We calculate numerically the entanglement entropy of free fermion ground states in one-, two-, and three-dimensional Anderson models and find that it obeys the area law as long as the linear size of the subsystem is sufficiently larger than the mean free path. This result holds in the metallic phase of the three-dimensional Anderson model, where the mean free path is finite although the localization length is infinite. Relation between the present results and earlier ones on area law violation in special one-dimensional models that support metallic phases is discussed.

1. Introduction

Recent years have witnessed tremendous progress in the study of entanglement in condensed matter/many-body physics. Among these studies, free fermion systems play a very special role [1]. Simple as they may seem, fermions are intrinsically nonlocal, due to the anticommutation relation fermion operators satisfy, no matter how far apart they are. Such nonlocality shows up as enhanced entanglement in the ground state; for example, for many years Fermi sea states were the only known ground states whose block entanglement entropy (EE) violates the area law satisfied by most ground states above 1D [2–4]. It is only recently shown that a similar violation occurs in interacting fermion systems in the Fermi liquid phase [5] and bosonic models with excitation spectra that vanish on (extended) Bose surfaces [6]. The existence of sharp Fermi or Bose surfaces is crucial for the area law violation in translationally invariant systems.

Comparatively much less effort has been devoted to studies of fermions in the presence of disorder potential. In a recent work [7] we studied two very special (one-dimensional) 1D models that exhibit free fermion metal-insulator transition (MIT) and found area law violation in the metallic phase, despite the presence of disorder, and thus absence of sharp Fermi surface (actually points in 1D). It was conjectured [7] that as long as the system is metallic, namely, states are delocalized at the Fermi energy, there will be area law violation. In the present work we test this conjecture by performing detailed numerical studies of the Anderson model [8] in one, two, and three dimensions. We find that the area law is actually respected in all cases, including the metallic phase in 3D. We do observe an enhancement (beyond area law) as systems sizes increase while below the mean free path; such enhancement disappears once the system size becomes sufficiently bigger than the mean free path. The origin of the difference between the Anderson model studied here and the special models studied earlier [7] will be discussed.

The remainder of the paper is organized as follows. In Section 2 we introduce our model and numerical method for calculating EE. Results of our calculations are presented in Section 3. Section 4 offers a summary and discussions on our results.

2. Model and Basic Considerations

Anderson model in dimension is a model with constant nearest neighbor hopping term and random on-site energy :

where summation is over all sites in dimensional hyper cubic lattice (with lattice constant set to be 1) and is a vector connecting a site to its nearest neighbor. ’s are uniformly distributed between and [9]. The Fermi energy is set to be 0 (so the lattice is half-filled) in all cases, while in 2D we also study to avoid the van Hove singularity at the band center. We consider cubic-shaped finite-size systems with linear size and open-boundary conditions. We then divide them into two equal subsystems and with size and calculate the disorder-averaged entanglement entropy as detailed below.

For a system in a pure state , the density matrix is . Reduced density matrix of each subsystem ( or ) is obtained by tracing over degrees of freedom of the other subsystem: . Block EE between the two subsystems is . For a single Slater-determinant ground state,

are characterized by free fermion entanglement Hamiltonians

where is determined by the normalization condition . We calculate EE using the method of [10] by diagonalizing correlation matrix of subsystem

and find its eigenvalues ’s. Then EE takes the form

where is number of sites in subsystem .

In one and two dimensions, all states are localized with any finite disorder. However there is an important difference between them: in 1D the localization length is of the same order as mean free path , while in 2D we have for weak disorder. In 3D there is a metal-insulator transition (MIT) at a critical value of disorder strength [11], where diverges. The focus of our numerical calculation is the interplay of the three different length scales, mean free path (calculated perturbatively in the Appendix), localization length , and (sub)system size , and their effects on entanglement.

3. Results

3.1. Anderson Model in One and Two Dimensions

In these two cases, all states are localized as long as .

Figure 1 shows 1D EE as a function of system size for different values of . As size increases, EE grows logarithmically for as expected. For , EE grows with in a manner similar to the disorder free case up to some point and then saturates, indicating area law is obeyed for sufficiently large system sizes. We find substantial deviation (from case) starts when the system size reaches the mean free path , and saturation occurs around . We note in 1D that we have the localization length ; it is thus not immediately clear at this point which of the two controls is the crossover.

Figure 1: Entanglement entropy of one-dimensional Anderson model in log-linear scale for different value of disorder strength . . Mean free path corresponding to selected ’s is indicated as a vertical line. Horizontal axis is the linear size of the system, , varying from up to . The number of samples at each point of is . The strange break in the small ’s including case is a consequence of small size effect in calculating average EE. We see that for bigger sizes the behavior is more smooth.

In Figures 2 and 3, 2D EE divided by boundary length (to account for the area law contribution) as a function of for different values of disorder strength is plotted for and , respectively. We see again that, for , there is a logarithmic growth as expected [2, 3], while for such growth stops beyond certain length scale, indicating area law behavior. This is particularly clear in Figure 2, for ; for smaller the deviation from the behavior and tendency toward saturation are also very obvious. In this case is not well defined perturbatively due to van Hove singularity, but the localization length is much bigger than the system sizes studied here (numerical calculations of [12] show that is in order of for and in order of for and exponentially bigger for smaller ’s), indicating plays no role in the size dependence of EE.

Figure 2: Entanglement entropy divided by linear system size of two-dimensional Anderson model in log-linear scale for different value of disorder strength . . Horizontal axis is the linear size of the system, . The total system has sites. The number of samples at each point of is .

Figure 3: Entanglement entropy divided by linear system size of two-dimensional Anderson model in log-linear scale for different value of disorder strength . . Mean free path corresponding to selected ’s is indicated as a vertical line. Horizontal axis is the linear size of the system, . The total system has sites. The number of samples at each point of is .

For (Figure 3), we again find deviation from the behavior, and tendency toward saturation starts when system size reaches the mean free path . We note in this case that we have for small ; for example, at we expect !. We thus again find that while controls the extensiveness of the fermion wave function at the Fermi level, it does not control the size dependence of EE.

3.2. Anderson Model in Three Dimensions

In three dimensions, there is a critical value of disorder strength, , where a metal-insulator transition occurs. States at the Fermi level are delocalized for . In Figure 4, EE divided by boundary area (to account for the area law contribution) as a function of for different values of disorder strength is plotted. Similar to the 1D and 2D cases, we find a logarithmic growth for , while for we find deviation from such growth, and tendency toward saturation starts when the system size reaches mean free path . It is particularly worth noting that there is no obvious change of the behavior of EE near , which is highlighted in the figure; area law behavior is clearly seen on both sides of . We thus conclude that entanglement area law is respected in both the metallic and insulating phases of the 3D Anderson model, as long as there is finite disorder strength.

Figure 4: Entanglement entropy divided by square of linear system size of three-dimensional Anderson model in log-linear scale for different value of disorder strength . . Mean free path corresponding to selected ’s is indicated as a vertical line. Data points corresponding to are indicated by blue square. Horizontal axis is the linear size of the system, . The total system has sites. The number of samples is for small sizes and for large sizes.

4. Discussion and Summary

In this work we find, through explicit numerical calculations, that entanglement entropy of free fermion systems obeys the area law in the ground states of one-, two-, and three-dimensional Anderson models, for subsystems whose linear size is sufficiently bigger than the mean free path. This result holds in the metallic phase of the three-dimensional case, where states at the Fermi energy are extended. Localization length, either finite or infinite, appears to play no special role in the size dependence of entanglement entropy.

This behavior is in sharp contrast to those of two special 1D models studied earlier [7], where we found violation of area law in the metallic phase. We now discuss the source of difference in these models. The first one is the random dimer model [13], which has a very special feature that there is no back scattering at a special resonance energy, as a result of which the system is metallic when the Fermi energy matches the resonance energy. In this case the mean free path is infinity (and the fermion motion is ballistic instead of diffusive as in ordinary diffusive metals), resulting in the area law violation, consistent with findings of the present work.

The second special model studied in [7] involves power-law long-range hopping [14], and the system is metallic when the decay exponent is small enough. In this case the fermion can hop over arbitrarily large distances; thus mean free path is not well defined. More importantly, in the metallic phase fermion motion is super diffusive [14], different from that in the 3D Anderson model.

Combining earlier results with the present ones, we conclude that while violation of entanglement area law does not necessarily require a sharp Fermi surface (which is destroyed by disorder), it does not occur in diffusive metals. It can occur, however, in metallic phases of disordered free fermion models where the fermion motion is super diffusive.

While in this paper we studied noninteracting fermions, a recent paper studied the role of localization length in the case of interacting fermions [15]. Also another paper [16] proved the area law in the localized regime of the Anderson model, although our numerical calculations demonstrate area law in the metallic phase as well.

Note Added. While the present paper is being written up, a related preprint [17] appeared on the arxiv reaching a very similar conclusion. The numerics of that work are limited to quasi-1D strips.

Appendix

Calculation of Mean Free Path as a Function of Random Potential Strength

We explain the calculation of mean free path here. We assume that ’s are uncorrelated random numbers with uniform distribution, that is, a white noise uniform distribution. First we need to calculate the ensemble averaged transition rate using Fermi’s golden rule:

in which stands for ensemble average and is the dimension. We calculate the mean free path at the Fermi level, which is chosen to be for 1D and 3D and for 2D. Also we know

where for a uniform distribution we have

Thus, finally we have

Then, , where is the average velocity of electron on Fermi surface. Calculated mean free path in 1, 2, and 3 dimensions at specific Fermi energy, , is

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

This research is supported by DOE Grant no. DE-SC0002140.

References

I. Peschel and V. Eisler, “Reduced density matrices and entanglement entropy in free lattice models,” Journal of Physics A. Mathematical and Theoretical, vol. 42, no. 50, Article ID 504003, 30 pages, 2009.

D. Gioev and I. Klich, “Entanglement entropy of fermions in any dimension and the widom conjecture,” Physical Review Letters, vol. 96, Article ID 100503, 2006.

M. M. Wolf, “Violation of the entropic area law for fermions,” Physical Review Letters, vol. 96, Article ID 010404, 2006.

B. Swingle, “Entanglement entropy and the Fermi surface,” Physical Review Letters, vol. 105, no. 5, Article ID 050502, 2010.

W. Ding, A. Seidel, and K. Yang, “Entanglement entropy of fermi liquids via multidimensional bosonization,” Physical Review X, vol. 2, no. 1, Article ID 011012, 18 pages, 2012.

H.-H. Lai, K. Yang, and N. E. Bonesteel, “Violation of the entanglement area law in bosonic systems with Bose surfaces: possible application to Bose metals,” Physical Review Letters, vol. 111, no. 21, Article ID 210402, 2013.

M. Pouranvari and K. Yang, “Maximally entangled mode, metal-insulator transition, and violation of entanglement area law in noninteracting fermion ground states,” Physical Review B, vol. 89, Article ID 115104, 2014.

P. W. Anderson, “Absence of diffusion in certain random lattices,” Physical Review, vol. 109, no. 5, pp. 1492–1505, 1958.

N. Laorencie, “Scaling of entanglement entropy in the random singlet phase,” Physical Review B, vol. 72, no. 14, Article ID 140408(R), 4 pages, 2005.

I. Peschel, “Calculation of reduced density matrices from correlation functions,” Journal of Physics A: Mathematical and General, vol. 36, no. 14, pp. L205–L208, 2003.

A. MacKinnon and B. Kramer, “One-Parameter scaling of localization length and conductance in disordered systems,” Physical Review Letters, vol. 47, article 1546, 1981.

A. MacKinnon and B. Kramer, “The scaling theory of electrons in disordered solids: additional numerical results,” Zeitschrift für Physik B: Condensed Matter, vol. 53, no. 1, pp. 1–13, 1983.

D. H. Dunlap, H.-L. Wu, and P. W. Phillips, “Absence of localization in a random-dimer model,” Physical Review Letters, vol. 65, no. 1, pp. 88–91, 1990.

A. D. Mirlin, Y. V. Fyodorov, F.-M. Dittes, J. Quezada, and T. H. Seligman, “Transition from localized to extended eigenstates in the ensemble of power-law random banded matrices,” Physical Review E—Statistical Physics, Plasmas, Fluids, and Related Interdisciplinary Topics, vol. 54, no. 4, pp. 3221–3234, 1996.

R. Berkovits, “Entanglement entropy in a one-dimensional disordered interacting system: the role of localization,” Physical Review Letters, vol. 108, no. 17, Article ID 176803, 5 pages, 2012.

L. Pastur and V. Slavin, “Area law scaling for the entropy of disordered quasifree fermions,” Physical Review Letters, vol. 113, Article ID 150404, 2014.

A. C. Potter, “Boundary-law scaling of entanglement entropy in diffusive metals,” http://arxiv.org/abs/1408.1094