Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: ground states of quantum systems with many constituents are typically entangled . in a similar manner as one can identify characteristic length scales of correlation functions , quantum correlations are expected to exhibit some general scaling behavior , beyond the details of a fine - grained description . such characteristic features provide a physical picture that goes beyond the specifics of the underlying microscopic model . a central question of this type is the following : if one distinguishes a certain collection of subsystems , representing some spatial region , of a quantum many - body system in a pure ground state , the state of this part will typically have a positive entropy , reflecting the entanglement between this region and the rest of the system @xcite . this degree of entanglement is certainly expected to depend on the size and also on the shape of the region . yet , how does the degree of entanglement specifically depend on the size of the distinguished region ? in particular , does it scale as the volume of the interior  which is meant to be the number of degrees of freedom of the interior ? or , potentially as the area of the boundary , i.e. , the number of contact points between the interior and the exterior ? this paper provides a detailed answer to the scaling behavior of the entanglement of regions with their exterior in a general setting of harmonic bosonic lattice systems and provides a comprehensive treatment of upper bounds on these quantities . we find that in arbitrary spatial dimensions the degree of entanglement in terms of the von  neumann entropy scales asymptotically as the area of the boundary of the distinguished region . this paper significantly extends the findings of ref.@xcite on harmonic bosonic lattice systems : there , the area - dependence of the geometric entropy has been proven for cubic regions in non - critical harmonic lattice systems of arbitrary dimension with nearest - neighbor interactions , corresponding to discrete versions of klein - gordon fields . in this paper we extend our analysis to a general class of finite - ranged harmonic interactions and also take regions of arbitrary shape into account . for thermal gibbs states , the entropy of a reduction is longer a meaningful measure of entanglement . instead , an area - dependence for an appropriate mixed - state entanglement measure , the distillable entanglement , is established . also , an analogous statement holds for classical correlations in classical systems . the area - dependence is even found in certain cases where one can prove the divergence of the two - point correlation length . this demonstrates that this previously conjectured dependence between area and entanglement is valid under surprisingly general conditions . the presented analysis will make use of methods from the quantitative theory of entanglement in the context of quantum information science @xcite . it has been become clear recently that on questions about scaling of entropies and degrees of entanglement  albeit often posed some time ago  new light can be shed with such methods @xcite . in this language , quantum correlations are sharply grasped in terms of rates that can be achieved in local physical transformations . to assess quantum correlations using novel powerful tools from quantum information and to relate them to information - theoretical quantities constitutes an exciting perspective .    in the context of quantum field theory , such questions of scaling of entropies and entanglement have a long tradition under the keyword of geometric entropy . in particular , work on the geometric entropy of free klein - gordon fields was driven in part by the intriguing suggested connection @xcite to the bekenstein - hawking black hole entropy @xcite . in seminal works by bombelli et al . @xcite and srednicki @xcite the relation between the entropy and the boundary area of the region has been suggested and supplemented with numerical arguments . this connection has been made more specific using a number of different methods . in particular , for half spaces in general and intervals in the one - dimensional case , the problem has been assessed employing methods from conformal field theory , notably @xcite , based on earlier work by cardy and peschel @xcite , and by cardy and calabrese @xcite .    in one - dimensional non - critical chains , one observes a saturation of the entanglement of a distinguished block , as was proven analytically for harmonic chains @xcite and was later observed for non - critical spin chains numerically @xcite and analytically @xcite . in turn , in critical systems systems , one often  but not always  finds a logarithmically diverging entropy @xcite . in case of a model the continuum limit of which leads to a conformal field theory , the factor of the logarithmically diverging term is related to the central charge of the conformal field theory . the findings of the present paper and of the above - mentioned results motivate further questions concerning the general area - dependence of the degree of entanglement , for example in fermionic systems @xcite . in particular , the connection between the geometric entropy fulfilling no area law and the correlated quantum many - body system being critical is not fully understood yet . this is particularly true for the interesting case of more than one - dimensional quantum systems . this paper is structured as follows . we start , in section [ major ] , with a presentation of the major results of the paper . in section  [ prelims ] we define our notation and recall some basics on harmonic lattice systems . section [ bounds ] provides a general framework of upper and lower bounds for entanglement measures , expressed in terms of the spectrum of the hamiltonian and two - point correlation functions . this analysis is performed both for the case of the ground state as well as for mixed gibbs states . of particular interest are hamiltonians with finite - ranged interactions in section [ sectionfiniteranged ] . for such hamiltonians , we first study the behavior of the two - point correlation functions , then the entanglement bounds . in section  [ nomore ] , discussing the gibbs state case , we determine temperatures above which there is no entanglement left . a class of examples of hamiltonians which exhibit a divergent two - point correlation length in their ground state , but an area - dependence of the entanglement is presented in section  [ divergent ] and expressed in analytical terms . the specific case of hamiltonians the interaction part of which can be expressed as a square of a banded matrix is discussed in section [ thasquared ] . in this case , very explicit expressions for entanglement measures can be found . we then consider , in section [ classical ] , the case of classical correlations in classical harmonic systems with arbitrary interaction structure . interestingly , this case is related to the quantum case for squared interactions in the sense of section  [ thasquared ] . finally , we summarize what has been achieved in the present paper , and present a number of open questions in this context . throughout the whole paper we will consider harmonic subsystems on a @xmath0-dimensional cubic lattice @xmath1^{\times d}$ ] . the system thus has @xmath2 canonical degrees of freedom . the central question of this work will be : how does the degree of entanglement of a distinguished region @xmath3 with the rest @xmath4 scale with the size and shape of @xmath5 ? ^{\times 2}$ ] with a distinguished ( shaded ) region @xmath5 , consisting of @xmath6 degrees of freedom  . oscillators representing the exterior @xmath4 are shown as  , whereas pairs belonging to the surface @xmath7 of the region are marked by lines ( ) . ] we define the volume @xmath6 and surface @xmath7 of the distinguished region @xmath5 as @xmath8 where @xmath9 defines the @xmath10-norm distance for vectors @xmath11 that specify the position of oscillators on the @xmath0-dimensional lattice . more specifically , @xmath7 is the number of contact points , that is , the number of pairs of sites of @xmath5 and @xmath12 that are immediately adjacent . note that the distinguished region @xmath5 may have arbitrary shape and does not have to be contiguous .    for the pure ground state @xmath13 , we will study the entropy of entanglement of @xmath3 given by @xmath14 where @xmath15 $ ] is the von  neumann entropy of a state @xmath16 and @xmath17 $ ] denotes the reduced state associated with the degrees of freedom of the interior @xmath5 . for the pure ground state , this entropy of entanglement is identical to both the distillable entanglement and the entanglement cost . for pure states , it is indeed the unique asymptotic measure of entanglement . for gibbs states , in turn , we have to study a mixed - state entanglement measure , as the entropy of a subsystem no longer meaningfully quantifies the degree of entanglement . we will bound the rate at which maximally entangled pairs can be distilled , i.e. , the distillable entanglement @xmath18 . clearly , @xmath19 for zero temperature . we will subsequently suppress the index @xmath20 for notational clarity . we derive the following properties of @xmath21 and @xmath22 :    _ _    * for @xmath0-dimensional harmonic lattice systems , we derive general upper and lower bounds to @xmath23 for pure ground states and to @xmath24 for gibbs states with respect to a temperature @xmath25 . these bounds are expressed entirely in terms of the potential matrix and stated in eqs . ( [ upperboundfinitet ] ) , ( [ upperboundzerot ] ) , and ( [ lowerbound ] ) . a necessary condition for the following results to hold is that the spectral condition number @xmath26 of the coupling matrix @xmath27 satisfies @xmath28 for some @xmath29 independent of @xmath5 and @xmath12 .    _ _    * for nearest - neighbor interactions and the ground state , the entropy of entanglement @xmath21 scales as the surface area @xmath7 of @xmath5 . more specifically , there exist numbers @xmath30 independent of @xmath12 and @xmath5 such that @xmath31 this is stated in eqs . ( [ finiteupperbound ] ) and ( [ nearestlowerbound ] ) . note that the specific case of cubic regions @xmath32^{\times d}$ ] has already been proven in the shorter paper ref.@xcite . * for general finite - ranged harmonic interactions  meaning arbitrary interactions which are strictly zero after a finite distance  the entropy of entanglement @xmath21 of the ground state scales at most linearly with the surface area of @xmath5 : there exists a @xmath29 independent from @xmath12 and @xmath5 such that @xmath33 this is expressed in eq . ( [ finiteupperbound ] ) . * for gibbs states with respect to a temperature @xmath25 and general finite - ranged interactions the distillable entanglement @xmath34 scales at most linearly with the surface area of @xmath5 , so @xmath35 with @xmath36 being independent of @xmath12 and @xmath5 , see again eq . ( [ finiteupperbound ] ) , which applies for any temperature . we also determine temperatures above which the entanglement is strictly zero , see eq . ( [ nomoretemp ] ) . * we construct a class of hamiltonians the ground state of which exhibit an infinite two - point correlation length , for which the entropy of entanglement @xmath21 scales provably at most linearly in the boundary area . this is stated in eq . ( [ divergenceexample ] ) . * for interactions that are specified by potential matrices @xmath27 ( see eq . ( [ potmat ] ) ) that can be written as @xmath37 with @xmath38 corresponding to a finite - ranged interaction , we determine the entropy of entanglement . this is made specific in particular for @xmath39 , see eq . ( [ squaredonedim ] ) . * for classical systems of harmonic oscillators with nearest - neighbor interactions prepared in a thermal state , the mutual information @xmath40 measuring classical correlations scales linearly with the surface area of the region . here , @xmath41 , @xmath42 , and @xmath43 are the discrete classical entropies in phase space of the interior , the exterior , and the entire system with respect to an arbitrary coarse graining as defined in eq . ( [ boltzmann ] ) , for a classical gibbs state at temperature @xmath25 . the mathematical formulation of the problem is intimately related to case  ( vi ) . +    the remainder of the paper is now concerned with the detailed discussion and rigorous proofs of these statements . we consider quantum systems on @xmath0-dimensional lattices , where each site is associated with a physical system . the starting point is the following hamiltonian which is quadratic in the canonical coordinates , @xmath44 the matrix @xmath27 is the potential matrix . @xmath45 denotes the potential energy contained in the degree of freedom labeled with @xmath46 . for @xmath47 the element @xmath48 describes the coupling between the two oscillators at @xmath46 and @xmath49 . we assume that @xmath50 is a real symmetric positive matrix . in this paper , we will consider finite - ranged and nearest - neighbor interactions . furthermore , interactions for which the correlation length diverges are considered . we will subsequently study properties of both the ground state @xmath51 as well as of gibbs states @xmath52}\ ] ] corresponding to some nonzero temperature @xmath53 . note that we have set the boltzmann constant @xmath54 . the questions we ask are essentially those of the geometric entropy ( @xmath55 ) , and those on the distillable entanglement ( @xmath25 ) with respect to a distinguished region @xmath5 of the lattice @xmath20 . the specific case of nearest - neighbor interactions , considered in refs.@xcite corresponds to the one of a discrete version of free klein - gordon fields in flat space - time . here , we consider regions of arbitrary shape and allow for more general interactions . the system we consider embodies @xmath2 canonical degrees of freedom , associated with a phase space @xmath56 , where the @xmath57-matrix @xmath58 , @xmath59,\ ] ] specifies the symplectic scalar product , reflecting the canonical commutation relations between the @xmath60 canonical coordinates @xmath61 of position and momentum . instead of considering states we will refer to their moments , considering that both the ground state as well as the gibbs states are quasi - free ( gaussian ) states . the second moments can be collected in the @xmath57 real symmetric covariance matrix @xmath62 ( see , e.g. , ref . @xcite ) , defined as @xmath63)=      2 \re \langle o_j o_k\rangle_{\varrho}.\ ] ] here , it is assumed that the first moments vanish , which does not restrict generality as they can be made to vanish by local unitary displacements that do not affect the entanglement properties . following ref . @xcite , we find from symplectic diagonalization that for the ground state at zero temperature @xmath64 i.e. , there is no mutual correlation between position and momentum , and the two - point vacuum correlation functions are given by the entries of @xmath65 respectively @xmath65 , @xmath66_{{\boldsymbol{i}},{\boldsymbol{j } } } ,      \label{twopoint1}\\      h_{{\boldsymbol{i}},{\boldsymbol{j } } }      & { \ensuremath{\;\raisebox{.06ex}{:}\!\!=}}\langle\text{gs}|p_{{\boldsymbol{i}}}p_{{\boldsymbol{j}}}|\text{gs}\rangle=[v^{1/2}]_{{\boldsymbol{i}},{\boldsymbol{j}}}. \label{twopoint2}\end{aligned}\ ] ] here , @xmath67 denotes the state vector of the ground state . for the thermal gibbs state at finite temperature the covariance matrix takes the form @xmath68 where we define @xmath69 note that the additional term @xmath70 is the same for the position as well as for the momentum canonical coordinates . as in the zero temperature case position and momentum are not mutually correlated and @xmath71 ( @xmath72 ) are the two - point correlation functions of position ( momentum ) , now with respect to the gibbs state . the entropy of entanglement can be expressed in terms of the symplectic eigenvalues of the covariance matrix corresponding to a reduction . let the @xmath73-matrix @xmath74 denote the covariance matrix associated with the interior @xmath5 ; this is the principle submatrix of @xmath75 associated with the degrees of freedom of the interior . the symplectic spectrum of @xmath74 is then defined as the spectrum @xmath76 of the matrix @xmath77 for the situation at hand we find that @xmath78 where we denote by @xmath79 the principle submatrix of @xmath80 associated with the interior and we index the entries of the submatrix as @xmath81_{{\boldsymbol{i}},{\boldsymbol{j}}}$ ] with vectors @xmath82 . given the direct - sum structure of @xmath83 and counting doubly degenerate eigenvalues only once , we find that the entropy of entanglement can be evaluated as @xmath84,\ ] ] where @xmath85 , and the @xmath86 are the square roots of the regular eigenvalues @xmath87 of the @xmath88 matrix @xmath89 , @xmath90 this reflects the fact that the entropy of a harmonic system is the sum of the entropies of uncoupled degrees of freedom after symplectic diagonalization . + using vectors @xmath82 to label the entries of @xmath89 , we find due to @xmath91 , @xmath92_{{\boldsymbol{i}},{\boldsymbol{j}}}&=\left[\left.v^{-1/2}\right|_i\left.v^{1/2}\right|_i\right]_{{\boldsymbol{i}},{\boldsymbol{j}}}\\ & = \delta_{{\boldsymbol{i}},{\boldsymbol{j}}}-\sum_{{\boldsymbol{k}}\in o}\left[v^{-1/2}\right]_{{\boldsymbol{i}},{\boldsymbol{k } } }      \left[v^{1/2}\right]_{{\boldsymbol{k}},{\boldsymbol{j}}}. \end{split}\ ] ]    [ lb]@xmath93 [ lb]@xmath94 [ lb]@xmath95 [ lb]@xmath96   in one dimension , @xmath39 , for a finite - ranged coupling matrix @xmath27 , @xmath97 $ ] and @xmath98\cup[65, ... ,70]$ ] , yielding @xmath99 . bars show @xmath100 , color - encoded surface depicts @xmath101 . all units are arbitrary . note that the entries decay exponentially away from the boundary of @xmath5.,title="fig : " ]    this form of @xmath102 hints at an area theorem in the following way ( see fig . [ r ] , where we depict the entries of @xmath103 ) : if the entries of @xmath80 decay fast enough away from the main diagonal , i.e. , if the correlation functions decay sufficiently fast , the main contribution to the entropy comes from oscillators inside a layer around the surface of @xmath5 as for all the others the product @xmath81_{{\boldsymbol{i}},{\boldsymbol{k } } } [ v^{1/2}]_{{\boldsymbol{k}},{\boldsymbol{j}}}$ ] will be very small ( here vectors @xmath104 are inside @xmath5 , @xmath105 outside the region @xmath5 ) . thus , the matrix @xmath103 has an effective rank proportional to @xmath7 , and the number of symplectic eigenvalues contributing to the sum in eq . ( [ entropy ] ) is approximately proportional to @xmath7 . much of the remainder of the paper aims at putting this intuition on rigorous grounds .    for mixed states , such as thermal gibbs states , the von  neumann entropy no longer represents a meaningful measure of the present quantum correlations , and has to be replaced by concepts such as the distillable entanglement . the distillable entanglement is the rate at which one can asymptotically distill maximally entangled pairs , using only local quantum operations assisted with classical communication . for pure states , it coincides with the entropy of entanglement  the entropy of a reduction  giving an operational interpretation to this quantity @xcite .    from entanglement theory we know that an upper bound for the distillable entanglement is provided by the logarithmic negativity @xcite . note that this is the case even in this infinite - dimensional context for gaussian states can be immediately verified on the level of second moments and a single - mode description . moreover , the entropy of entanglement indeed still has the interpretation of a distillable entanglement , albeit the fact that distillation protocols leave the gaussian setting . this is true as long as one includes an appropriate constraint to the mean energy in the distillation protocol @xcite . the logarithmic negativity @xcite is defined as @xmath106 where @xmath107 denotes the trace norm , and @xmath108 is the partial transpose of @xmath109 with respect to the split @xmath5 and @xmath4 . again following ref . @xcite , we find after a number of steps @xmath110,\ ] ] where @xmath111 labels the @xmath2 eigenvalues of @xmath112 , @xmath113 and matrices @xmath114 are defined as @xmath115 which become @xmath116 for zero temperature following eq . ( [ wt ] ) . the diagonal matrix @xmath117 , defined as @xmath118 is the matrix that implements time reversal in the subsystem corresponding to the inner part @xmath5 , reflecting partial transposition @xmath119 on the level of states . in this section we will derive upper and lower bounds for the entropy of entanglement and the distillable entanglement of the distinguished region @xmath5 with respect to the rest of the lattice . these bounds only depend on the geometry of the problem , i.e. , the region @xmath5 and on properties of @xmath27 , namely its minimal and maximal eigenvalue , which we define as @xmath120 and @xmath121 , respectively , its condition number @xmath122 and the entries of @xmath123 respectively the two - point correlation functions @xmath124 and @xmath125 as in eqs . ( [ twopoint1],[twopoint2 ] ) . in later sections these bounds will be made specific for a wide range of interaction matrices @xmath27 . an upper bound for the entropy of entanglement and the distillable entanglement is provided by the logarithmic negativity as in eq . ( [ logneg ] ) . utilizing this fact , we derive upper bounds using @xmath126-norms @xcite in this section . for brevity and clarity , we will present the case for @xmath55 and finite temperature in a single argument . a direct calculation shows that the matrix @xmath112 is given by @xmath127 introducing the matrix @xmath128 with entries @xmath129 therefore , we can bound the eigenvalues of @xmath112 according to @xmath130 where we denote by @xmath131 the smallest eigenvalue of @xmath70 which is given by @xmath132 hence , we can write @xmath133          \nonumber          \\          & \le&\frac{1}{\ln(2 ) }          \sum_{{\boldsymbol{i}}\in c }          \max\left\{0,\lambda_{\text{min}}(w(t))^{-2}-1 +          2\left|          \lambda_{{\boldsymbol{i}}}(x\omega^+)\right|\right\},\end{aligned}\ ] ] i.e. , for @xmath134 , there is no longer any bi - partite entanglement in the system . we will later see that there is a temperature @xmath135 above which this happens . but for now , we use the fact that @xmath136 to bound the logarithmic negativity , and relate it to the @xmath126 norm @xmath137 @xcite of @xmath138 . we have @xmath139 and therefore @xmath140 the @xmath126 norm is defined as the sum of the absolute values of all matrix entries . inserting the definition of @xmath128 , we find @xmath141 for finite temperature , and @xmath142_{{\boldsymbol{i}},{\boldsymbol{j}}}\right|         \label{upperboundzerot}\ ] ] for zero temperature . these constitute upper bounds on the entropy of entanglement and the distillable entanglement . both only depend in the distinguished region @xmath5 , the maximum eigenvalue of @xmath27 , and the entries of @xmath65 , i.e. , the two - point correlation function with entries @xmath143 . this will be the starting point to derive explicit upper bounds for special types of interaction matrices @xmath27 .      to achieve lower bounds for the zero temperature case , we consider the entropy of entanglement directly . starting from the general expression eq . ( [ entropy ] ) , we use the fact that the symplectic eigenvalues are never smaller than @xmath10 , i.e. , that the eigenvalues of @xmath89 are contained in the interval @xmath144 $ ] , with @xmath145 being the maximal eigenvalue of @xmath89 . using the pinching inequality @xcite , we find @xmath146 . thus we can bound the entropy of entanglement as follows , @xmath147      = \sum_{i=1}^{v(i ) } \frac{\log [ \lambda_i]}{2 }      \ge      \frac{\log[\alpha_i]}{2(\alpha_i-1)}\sum_{i=1}^{v(i)}\left(\lambda_i-1\right)\\      & \ge &      \frac{\log[\sqrt{\kappa}]}{2(\sqrt{\kappa}-1)}\operatorname{tr}\left[a_i-{{\mathbbm{1}}}\right].\end{aligned}\ ] ] using vectors @xmath82 to label the entries of the @xmath88 matrix @xmath89 , we finally arrive at the lower bound @xmath148}{2(\sqrt{\kappa}-1 ) }      \sum_{\substack {        { \boldsymbol{i}}\in i \\        { \boldsymbol{j}}\in o } }      \left[v^{-1/2}\right]_{{\boldsymbol{i}},{\boldsymbol{j } } }      \left[v^{1/2}\right]_{{\boldsymbol{j}},{\boldsymbol{i}}}.\;\;\;\ ;      \label{lowerbound}\end{aligned}\ ] ] this lower bound depends only of the geometry of @xmath5 , the spectral condition number @xmath149 of @xmath27 and the entries of @xmath123 , i.e. , the two - point correlation functions . ( [ lowerbound ] ) is difficult to evaluate in general but for special cases of the interaction matrix @xmath27 it can nevertheless be made specific . in section  [ nearestneighbor ] , for example , we give an explicit expression for the important case of nearest - neighbor interactions . we expect eq . ( [ lowerbound ] ) to be a convenient starting point to derive such lower bounds also for more general cases of @xmath27 . a lower bound for the finite temperature case is difficult to obtain . generally speaking , a lower bound to the distillable entanglement ( two - way distillable entanglement in our case ) is given by the hashing inequality @xcite , @xmath150 where @xmath151 $ ] and @xmath152 $ ] . yet , naively applied , this inequality will vanish , as generally @xmath153 and @xmath154 : this can be made intuitively clear from the following argument . the above analysis demonstrates that both the interior @xmath5 and the exterior @xmath12 can be approximately disentangled with local unitaries up to a layer of the thickness of the two - point correlation length . that is , degrees of freedom associated with @xmath155 for which @xmath156 is sufficiently large for every @xmath157 , can to a very good approximation be decoupled and unitarily transformed into a thermal gibbs state . each such degree of freedom will therefore contribute a constant number to @xmath158 , such that @xmath153 for sufficiently large @xmath159 . similarly , one can argue to arrive at @xmath154 . in order to establish a lower bound to the distillable entanglement , however , we may start with any protocol involving only local quantum operations and classical communication , and apply the hashing inequality on the resulting quantum states . this first step can include in particular local filterings . a bound linear in the boundary area is expected to become feasible if one first applies an appropriate unitary both in @xmath5 and @xmath12 , and then performs a local filtering involving degrees of freedom associated with @xmath160 and @xmath161 for which @xmath162 . this option will be explored elsewhere . in this section we will make the upper bound on the entropy of entanglement and the distillable entanglement explicit for symmetric finite - ranged interaction matrices @xmath27 , i.e. , matrices for which @xmath163 where @xmath156 denotes again the @xmath10-norm distance . denoting as before the maximum and minimum eigenvalue of @xmath27 as @xmath164 and @xmath165 respectively , we require that the spectral condition number @xmath166 be strictly less than infinity independent of @xmath20 , i.e. , independent of @xmath159 . note that we do not require any further assumptions on the matrix @xmath27 . we will make use of a result of ref . @xcite concerning the exponential decay of entries of matrix functions . after generalizing this result to matrices @xmath27 with the properties specified above it enables us to bound the entries of @xmath167 as follows ( see appendix  [ decay ] ) @xmath168 and @xmath169 where @xmath170 for zero temperature we then find @xmath171 this shows that off - diagonal terms of @xmath65 decay exponentially , cf . [ sqrtinvv ] . [ cl]@xmath172 [ cb]@xmath172 [ cc]@xmath172 [ cc]@xmath173 [ cc]@xmath174 [ cc]@xmath175    [ c]@xmath176    _ { { \boldsymbol{i}},{\boldsymbol{j}}}$ ] for a finite - ranged coupling matrix @xmath27 in two dimensions , @xmath177 , @xmath178 . bars show @xmath81_{{\boldsymbol{i}},{\boldsymbol{j}}}$ ] and the color - encoded surface shows @xmath179_{{\boldsymbol{i}},{\boldsymbol{j}}})$ ] . the inset depicts the same for @xmath180 . all units are arbitrary . note the exponential decay away from the main diagonal @xmath181 . ]    substituting eqs . ( [ entrybounds ] ) into the general result in ( [ upperboundfinitet ] ) , we find @xmath182 where @xmath183 note that @xmath184 coincides with the definition of the surface area of @xmath5 . we find @xmath185 ( see appendix  [ abzaehlerei ] ) , i.e. , we have @xmath186 thus we finally arrive at the desired upper bound linear in the surface area of @xmath5 for both the entropy of entanglement and the distillable entanglement , @xmath187 which becomes @xmath188 for zero temperature , i.e. , these upper bounds depend solely on the maximal and minimal eigenvalue of the interaction matrix @xmath27 and the surface area of the region @xmath5 . this result demonstrates that indeed , an area - bound of the degree of entanglement holds in generality for bosonic harmonic lattice systems . this shows that the previously expressed intuition can indeed be made rigorous in form of an analytical argument . notably , the derived results hold also for systems in which the coupling coefficients are not identical , but independent realizations of random variables . if the coefficients @xmath189 of the real symmetric matrix @xmath27 are taken from a distribution with a carrier @xmath190 $ ] , such that the interaction is ( i ) finite - ranged , and ( ii ) the carrier is chosen such that @xmath191 then the same result holds true . this follows immediately from the considerations in appendix  [ decay ] , where it is not assumed that the hamiltonian exhibits translational symmetry . ( [ finiteupperbound ] ) is thus valid also for disordered harmonic lattice systems . in section  [ sectionupperbound ] we found that for @xmath192 there is no entanglement between the regions @xmath5 and @xmath12 . for finite - range interactions @xmath27 we found in section [ sectionfiniteranged ] @xmath193 i.e. , we have @xmath194 on the right - hand side only @xmath195 and @xmath196 depend on the temperature , both are decreasing in @xmath53 and go to zero as @xmath53 goes to infinity , i.e. , @xmath197 gives an implicit equation for the temperature @xmath135 above which there is no bi - partite entanglement left in the system . [ nearestneighbor ] in this section we consider nearest - neighbor interactions and periodic boundary conditions in @xmath0 spatial dimensions , i.e. , block - circulant matrices @xmath27 . for @xmath198 this is a special case of the matrices considered in section  [ sectionfiniteranged ] as boundary conditions become irrelevant in this limit , i.e. , the upper bound coincides with the one derived for finite - ranged interactions . for a tighter bound see appendix  [ nnentries ] . we will now make use of the circulant structure of @xmath27 to show that for matrices of this kind it is possible to also derive a lower bound on the entropy of entanglement that it is proportional to the surface area of @xmath5 . we write @xmath199 for the circulant matrix @xmath200 whose first block column is specified by the tupel of matrices @xmath201 . we can then recursively define @xmath202 via @xmath203 where @xmath204 . we choose an energy scale in which @xmath205 and we demand @xmath29 @xcite . this circulant structure leads to the following properties of @xmath27 : the eigenvalues of @xmath27 are given by @xmath206 in particular the maximum eigenvalue is given by @xmath207 , where equality holds for @xmath159 even . the minimum eigenvalue reads @xmath208 , i.e. , positivity of @xmath27 demands @xmath209 . note that the assumption that @xmath209 independent of @xmath5 and @xmath12 is essential for the argument . if we allow for an @xmath159-dependence of @xmath210 as it arises , e.g. , in the field - limit where @xmath211 , then in one spatial dimension one will encounter an area law up to a logarithmically divergent correction . this behavior  resembling the behavior of critical spin chains and quadratic fermionic models  will be studied in more detail elsewhere . furthermore the circulant nature of @xmath27 yields the following explicit expressions for the entries of @xmath123 @xmath212_{{\boldsymbol{i}},{\boldsymbol{j}}}= \frac{1}{n^d}\sum_{{\boldsymbol{k}}\in c}{\mathrm{e}}^{2\pi{\mathrm{i}}{\boldsymbol{k}}({\boldsymbol{i}}-{\boldsymbol{j}})/n}\left(\lambda_{{\boldsymbol{k}}}(v)\right)^{\pm 1/2}.\ ] ] this reasoning leads to the following properties for @xmath47 that are crucial for the present proof ( see appendix  [ nnentries ] ) , @xmath213_{{\boldsymbol{i}},{\boldsymbol{j}}}[v^{1/2}]_{{\boldsymbol{j}},{\boldsymbol{i } } } = \left|[v^{-1/2}]_{{\boldsymbol{i}},{\boldsymbol{j}}}[v^{1/2}]_{{\boldsymbol{j}},{\boldsymbol{i}}}\right|\\ \ge \left(\frac{c}{2}\right)^{2d({\boldsymbol{i}},{\boldsymbol{j } } ) } \big(\frac{d({\boldsymbol{i}},{\boldsymbol{j}})!}{(2d({\boldsymbol{i}},{\boldsymbol{j}})-1)\prod_{\delta=1}^d|i_\delta - j_\delta|!}\big)^2.\end{gathered}\ ] ] substituting these results into the general lower bound ( [ lowerbound ] ) and keeping only terms with @xmath162 immediately yields @xmath214c^2}{8(\sqrt{\kappa}-1)}s(i ) , \ ] ] where @xmath215 in this case . this generalizes the result of ref.@xcite to regions @xmath5 of arbitrary shape . we expect that these lower bounds can be generalized to other interactions @xmath27 and numerical results suggest that these bounds hold quite generally . a simple special case is related to a certain kind of interaction : this the one where the interaction matrix @xmath27 can be written as @xmath216 with a real symmetric matrix @xmath200 corresponding to a finite - range interaction . then , the covariance matrix associated with the interior is nothing but @xmath217 in this case , the symplectic spectrum can be determined in a fairly straightforward manner . we mention this case also as it appears to be an appropriate toy model as a starting point for studies aiming at assessing the symplectic spectrum of the reduction itself and therefore the spectrum of the reduced density matrix @xcite . we are looking for the spectrum of the @xmath88 matrix @xmath89 as in eq . ( [ amatrix ] ) which now takes the form ( recall that we label the entries of @xmath89 by vectors @xmath218 ) @xmath219_{{\boldsymbol{i}},{\boldsymbol{j}}}=\delta_{{\boldsymbol{i}},{\boldsymbol{j}}}-\sum_{{\boldsymbol{k}}\in o}[m^{-1}]_{{\boldsymbol{i}},{\boldsymbol{k } } }      m_{{\boldsymbol{k}},{\boldsymbol{j } } }      = \delta_{{\boldsymbol{i}},{\boldsymbol{j}}}-r_{{\boldsymbol{i}},{\boldsymbol{j}}}.\ ] ] it is now the central observation that the number of rows of the matrix @xmath103 that are nonzero  and therefore also the number of eigenvalues that are nonzero  is proportional to the surface area of @xmath5 , as @xmath200 is a banded matrix . that is , @xmath220 for @xmath221 , meaning that @xmath222 if @xmath223 for all @xmath224 . as the eigenvalues of @xmath89 are bounded through the pinching inequality by @xmath10 and @xmath225 , an upper bound linear in the surface area of the interior follows immediately from the fact that the number of eigenvalues entering the sum in eq . ( [ entropy ] ) is proportional to @xmath7 . this argument demonstrates that in this simple case , one immediately arrives at bounds that are linear in the number of contact points . the task of finding the eigenvalues of @xmath89 explicitly is now reduced to finding the eigenvalues of the sparse matrix @xmath103 , @xmath226 this case is particularly transparent in the one - dimensional case , @xmath39 , and for @xmath200 being a circulant matrix with first row @xmath227 for @xmath228 . the potential matrix @xmath27 corresponds then to nearest - neighbor interactions , together with next - to - nearest - neighbor interactions . in this one - dimensional setting , we set @xmath32 $ ] . the matrix @xmath103 then takes the simple form ( cf . [ sq_r ] ) @xmath229_{i , k}m_{k , j}\\      & = & -c\left([m^{-1}]_{i , m+1}\delta_{j , m}+[m^{-1}]_{i , n}\delta_{j,1}\right).\end{aligned}\ ] ]    [ cl]@xmath172 [ bc]@xmath172    [ lt]@xmath230 [ lt]@xmath231 [ lt]@xmath232 [ lt]@xmath233 [ lt]@xmath234 [ cl]@xmath235   for the case @xmath39 , @xmath97 $ ] , @xmath236 $ ] , @xmath237 , where @xmath200 is a nearest - neighbor circulant matrix . note that for this particular form of @xmath27 the entries @xmath100 are exactly zero for @xmath238 . this is in contrast to the general form of @xmath103 depicted in fig . the inset shows the entries @xmath239 as bars and @xmath240 in blue ( cf . ( [ rdecay])).,title="fig : " ]    to find its eigenvalues , we calculate @xmath241 , which is straightforward for a matrix of this form . @xmath242 reflecting the fact that @xmath243 eigenvalues are zero , i.e , the number of nonzero eigenvalues is @xmath244)=2 $ ] . we find for the non - vanishing eigenvalues @xmath245 symmetry of @xmath246 yields @xmath247_{1,n}\pm [ m^{-1}]_{m , n}\right).\ ] ] from the circulant structure of @xmath200 we have @xmath248_{i , j}&= & \frac{1}{n}\sum_{k=1}^n\frac{{\mathrm{e}}^{2\pi{\mathrm{i}}k(i - j)/n}}{1 - 2c\cos(2\pi k / n)}\\ & = & \frac{1}{n}\sum_{k=1}^n\frac{\cos(2\pi k(i - j)/n)}{1 - 2c\cos(2\pi k / n)},\\ \ ] ] i.e. , for large @xmath159 we arrive at @xmath248_{1,n}&=&\frac{1}{2c\sqrt{1 - 4c^2}}-\frac{1}{2c},\nonumber\\\ [ m^{-1}]_{m , n}&=&\frac{1}{\sqrt{1 - 4c^2}}\left(\frac{1-\sqrt{1 - 4c^2}}{2c}\right)^m.\label{rdecay}\end{aligned}\ ] ] note that these expressions are asymptotically independent of @xmath249 as shown in the inset of fig . [ sq_r ] . in this limit we finally arrive at @xmath250 this expression specifies the symplectic spectrum of the reduction in a closed form . in this section we will analytically demonstrate that there exist hamiltonians for which the ground state two - point correlation functions diverge whereas the geometric entropy is still bounded by an area - law . in spin systems with a long - range ising interaction such a behavior has been observed in the one - dimensional case and sketched for higher dimensions @xcite . here , we present a class of examples of this type valid in arbitrary dimension @xmath0 , for which one can prove the validity of an upper bound linear in the boundary area . moreover , this set of examples is not restricted to cubic regions . the interaction is here a suitable harmonic long - range interaction . this analysis shows that a divergent two - point correlation function alone is no criterion for a saturating block entropy in the one - dimensional case , and for an area - dependence in higher dimensions . consider the matrix @xmath200 with entries @xmath251 for some @xmath252 . now set @xmath253 . this choice implies that the correlation function @xmath254_{{\boldsymbol{i}},{\boldsymbol{j}}}$ ] decays only algebraically . we will show that despite this fact one still has an upper bound linear in the boundary area of @xmath5 for appropriate values of @xmath255 . firstly , we have to make sure that the maximum eigenvalue of @xmath27 can be bounded from above independent of @xmath159 . from gershgorin s theorem ( see , e.g. , ref . @xcite ) we know that for every eigenvalue @xmath256 there exists an @xmath257 such that @xmath258 i.e. , @xmath259=1,\ ] ] and therefore @xmath260 substituting this and the specific form of @xmath200 into the general expression for the entropy , we obtain @xmath261^\alpha } = \sum_{l=1}^\infty \frac{1}{l^\alpha}n_l\\ & \le & 2^{2d}s(i)\sum_{l=1}^\infty l^{2d-\alpha-1},\end{aligned}\ ] ] which converges for @xmath262 to @xmath263 where @xmath264 is the riemann zeta function . note that these bounds are not necessarily tight in the sense that even for smaller values of @xmath255 , such a behavior can be expected . steps towards tightening these bounds seem particularly feasible in case of cubic regions , where we conjecture that for @xmath265 we arrive at an area - dependence . this analysis shows that for long - range interactions , an area - law in the degree of entanglement can be concomitant with divergent two - point correlation functions . in previous sections we have considered the entanglement between some region and the rest of a lattice of interacting quantum harmonic oscillators with hamiltonians that are at most quadratic in position and momentum . we have demonstrated that both for the ground state and the thermal state of the entire lattice the quantum correlations , i.e. , the amount of entanglement , between the region and the remainder of the lattice is bounded by quantities proportional to their boundary surface area . this suggests similar questions concerning the classical correlations between a region and the rest of the lattice in the corresponding classical systems when the lattice system as a whole is prepared in a thermal state . in this section we will demonstrate that indeed analogous area laws hold . we will note furthermore that there is a quite striking intimate relation between the _ classical system _ with potential matrix @xmath266 and the _ quantum system _ with a potential matrix @xmath267 in this context .      for the following considerations we use the classical equivalent of the quantum mechanical hamiltonian ( [ potmat ] ) , namely again @xmath268 where now @xmath269 and @xmath270 are vectors of classical position and momentum variables , respectively , and @xmath27 denotes again the potential matrix . the state of the classical system is characterized by a phase space density @xmath271 , a classical probability distribution , where @xmath272 denotes all the canonical coordinates in phase space . for nonzero temperatures @xmath25 this phase space distribution is given by the boltzmann distribution @xmath273 where @xmath274 ( as before , @xmath54 ) , and @xmath275    given this density in phase space we will encounter a familiar ambiguity when defining the entropy of the system . using the discrete classical entropy @xcite , we split the phase space into cubic cells each with a volume @xmath276 , where @xmath277 and @xmath278 being an arbitrary constant . from the phase space density we obtain the probability associated with each of these cells which in turn can be used to determine the entropy function of this probability distribution . we will now make use of the multiple indices @xmath279 and @xmath280 , assuming that the cell corresponding to @xmath281 is centered around @xmath282 and @xmath283 . that is , for each degree of the @xmath2 degrees of freedom , the phase space is discretized . the contribution of each cell to the discretized probability density is then given by @xmath284 as usual , the discrete classical entropy is then defined as the corresponding shannon entropy @xmath285 we will denote the discrete classical entropy with respect to the degrees of freedom of the interior with @xmath286 , the entropy of the exterior with @xmath287 .    the value for the entropy will depend on the choice of @xmath288 and in the limit @xmath289 this entropy definition will diverge due to a term proportional to @xmath290 . in classical statistical mechanics this problem is avoided with the help of the third law of thermodynamics . the entropy itself is however not the quantity that we wish to compute but rather the mutual information @xmath291 between the interior @xmath5 and the rest of a lattice , denoted as before by @xmath12 . this classical mutual information meaningfully quantifies the classical correlations between the inner and the outer . in that case we find that the limit @xmath289 exists and that the mutual information @xmath5 can be defined as @xmath292 following these preparations we are now in a position to determine the mutual information between a region and the rest of the lattice explicitly when the lattice as a whole is in a thermal state .      for the evaluation of the mutual information we need to determine the entropy of the total lattice @xmath20 , as well as the entropy determined by the reduced densities describing the two regions @xmath5 and @xmath12 . to this end we carry out the partial summation over all degrees of freedom of region @xmath12 in order to find the reduced phase space density @xmath293 describing region @xmath5 only . employing the schur complement we find that the reduced density @xmath293 is described by the boltzmann distribution corresponding to the same temperature and the hamiltonian @xmath294_{{\boldsymbol{i}},{\boldsymbol{j } } }      x_{{\boldsymbol{j } } }      +      \sum_{{\boldsymbol{i}}\in i }      p_{{\boldsymbol{i}}}^2\right ) \ ] ] an analogous result holds for the reduced phase space density of region @xmath12 .    for a thermal phase space distribution eq . ( [ density ] ) corresponding to a classical hamiltonian function of the form eq . ( [ classicalham ] ) we can compute the entropy straightforwardly , to find @xmath295 @xmath296 , which increases with temperature as expected . for the mutual information we find @xmath297 which is , perhaps surprisingly , independent of temperature . using jacobi s determinant identity @xmath298 this expression can be rewritten as @xmath299=      \frac{1}{2 } \log      \det\left [ { { \mathbbm{1}}}-r\right].\ ] ] it is now advantageous to notice the close connection of this expression , in particular of the matrix @xmath103 @xmath300_{{\boldsymbol{i}},{\boldsymbol{k}}}v_{{\boldsymbol{k}},{\boldsymbol{j}}}\ ] ] with those that arise in the quantum mechanical problem that we have treated previously in section  [ thasquared ] . indeed , the classical problem for a system with potential matrix @xmath301 is related to the quantum mechanical system with the squared potential matrix @xmath302 . this formal similarity arises because in section  [ thasquared ] we have shown that for a lattice of quantum harmonic oscillators with potential matrix @xmath303 in its ground state the symplectic eigenvalues of the covariance matrix describing region i alone are exactly the standard eigenvalues of the matrix @xmath304 . the properties of these eigenvalues have already been discussed in detail in section  [ thasquared ] . this allows us now straightforwardly to establish the area theorem for the mutual information in a classical system employing the result for the corresponding quantum system . this establishes in particular that the classical correlations as measured by the mutual information @xmath305 between the distinguished region @xmath5 and the rest of the lattice @xmath4 satisfy @xmath306 for large @xmath6 and appropriate constants @xmath307 independent of @xmath12 and @xmath5 . ( compare also the assessment of the thermodynamical entropy of parts of classical fluids in ref.@xcite . ) in summary we have seen that the area - dependence of correlations is not restricted to quantum systems , as long as one replaces the notion of entanglement  representing quantum correlations  by the notion of classical correlations in a classical system . in this paper , we have considered the question of the area - dependence of the geometric entropy and the distillable entanglement in general bosonic harmonic lattice systems of arbitrary dimension . the question was the general scaling behavior of these measures of entanglement with the size of a distinguished region of a lattice . such an analysis generalizes assessments of block entropies in the one - dimensional case . using methods from entanglement theory , we established bounds that allow for a conclusion that may be expressed in a nutshell as : _ in surprising generality , we find that the degree of entanglement scales at most linearly in the boundary area of the distinguished region . _ this analysis shows that the intuition that both the interior and the exterior can be approximately disentangled up to a layer of the thickness of the two - point correlation length by appropriate local unitaries carries quite far indeed .    for cubic regions @xmath308^{\times d}$ ] the area law can be formulated as @xmath309 where @xmath310 is the landau theta . such area - laws are expected to have an immediate implication on the accuracy to which ground states can be approximated with matrix - product states and higher dimensional analogues in classical simulations of the ground states of quantum many - body systems @xcite . after all , the failure of dmrg algorithms close to critical points can be related to the logarithmic divergence of the block entropy in the one - dimensional case . the findings of the present paper raise a number of interesting questions . notably , in general quantum many - body systems on a lattice ( fermionic or bosonic ) , what are necessary and sufficient conditions for an area - law in the above sense to hold ? clearly , as we have seen above , the divergence of two - point correlation functions alone is not in one - to - one correspondence with an area law . it would be interesting to consider and possibly decide the conjecture that a one - to - one relationship between a system being critical and not satisfying a law of the form as in eq . ( [ thearea ] ) holds if one ( i ) restricts attention to systems in arbitrary dimension with nearest - neighbor interactions , and ( ii ) grasps criticality in terms of two - point correlation functions with algebraic decay , concomitant with a vanishing energy gap . note that the latter two criteria of criticality have not been rigorously related to each other yet and may indeed not be simultaneously satisfied in lattice systems . the general relationship is still awaiting rigorous clarification .    as steps towards such an understanding of a relationship between criticality and properties of ground state entanglement in more than one - dimensional fermionic and bosonic systems , it seems very interesting to study models different from the ones considered in this paper . for example , complementing our bosonic analysis , area - laws in fermionic critical systems have been addressed @xcite , where logarithmic corrections have been found also in higher - dimensional settings . other settings in the bosonic case , corresponding to field theories beyond this quasi - free setting , are also still not clarified . in particular , the scaling behavior of the geometric entropy in general bosonic theories in higher dimensions is far from clear . then , the case of finite - size effects in harmonic lattice systems where the correlation length is larger than the full system , resembling the critical case , will be presented elsewhere . in less generality , it seems also feasible to identify the prefactors of the leading and next - to - leading order terms in an area - law of the geometric entropy . it is the hope that the present work can contribute to an understanding of genuine quantum correlations in quantum many - body systems and inspire such further considerations . we would like to thank c.  h.  bennett , h.  j.  briegel , j.  i.  cirac , w.  dr , m.  fleischhauer , b.  reznik , t.  rudolph , n.  schuch , r.  g.  unanyan , r.  f.  werner , and m.  m.  wolf for discussions . this work has been supported by the epsrc qip - irc ( gr / s82176/0 ) , the eu thematic network quprodis ( ist-2002 - 38877 ) , the dfg ( spp 1078 and spp 1116 ) , the leverhulme trust ( f/07 058/u ) and the european research councils ( euryi ) . the result concerning the exponential decay of entries of matrix functions of ref.@xcite relies on the fact that the @xmath311-th power @xmath312 of a @xmath313-banded matrix @xmath314 is @xmath315-banded , i.e. , @xmath316_{i , j}=0 $ ] for @xmath317 for a matrix @xmath318 with @xmath319 for @xmath320 . for the purposes of the present paper , we will need a generalization of the result of ref . @xcite to block banded matrices . we refer to @xmath321 as being @xmath313-banded , if @xmath322 for @xmath221 . it can be proven by induction over @xmath311 that the @xmath311-th power of @xmath27 is @xmath315 banded in this sense . this enables us to formulate the general form of ref.@xcite as follows . let @xmath321 be a @xmath313-banded symmetric matrix , i.e. , @xmath322 for @xmath323 . define @xmath164 , @xmath165 , @xmath166 , @xmath324 and @xmath325 as an ellipse with foci in @xmath326 and @xmath10 and half axes @xmath255 , @xmath327 , @xmath328 . + now let @xmath329 be such that @xmath330 is analytic in the interior of the ellipse @xmath325 , @xmath331 , and continuous on @xmath325 . furthermore suppose @xmath332 for @xmath333 . then there exist constants @xmath334 and @xmath335 , @xmath336 , @xmath337 such that @xmath338_{{\boldsymbol{i}},{\boldsymbol{j}}}\right| \le k q^{\sum_{\delta=1}^d|i_\delta - j_\delta|},\ ] ] where @xmath339 to bound the entries of @xmath167 , we apply the above theorem to the function @xmath340 for @xmath341 , we have that @xmath342 is analytic in the interior of the ellipse @xmath325 and continuous on @xmath325 , i.e. , @xmath167 satisfy the assumptions of the above theorem . to make @xmath334 and @xmath335 specific , we choose @xmath343 which yields @xmath344 where @xmath345 and for zero temperature @xmath195 reduces to @xmath346 these findings explicitly relate the spectral properties of the hamiltonian to the two - point correlation functions . in this appendix , we make the evaluation of the entries of @xmath347 and @xmath65 of the important case of nearest - neighbor interactions specific . the power series expansion of the square root is given by @xmath348 which is valid for @xmath349 , i.e. , the positivity constraint @xmath350 allows us to write @xmath351 using the multinomial theorem , we have @xmath352 where the sum runs over all @xmath353 with @xmath354 . now , applying the binomial theorem @xmath355 substituting all the above into eq . ( [ entries ] ) , we find @xmath356_{{\boldsymbol{i}},{\boldsymbol{j}}}= \delta_{{\boldsymbol{i}},{\boldsymbol{j}}}\mp\sum_{l=1}^\infty a_l^\pm c^l \sum l!\prod_{\delta=1}^df(d_\delta , n_\delta),\ ] ] where @xmath357 and @xmath358 to the sum over @xmath359 only terms with @xmath360 for some @xmath361 contribute . we thus arrive at the following expression for the entries of @xmath123 @xmath356_{{\boldsymbol{i}},{\boldsymbol{j}}}=\\      \delta_{{\boldsymbol{i}},{\boldsymbol{j}}}\mp\sum_{l=1}^\infty a_l^\pm c^l\!\!\!\ !      \sum_{\sum_\delta n_\delta = l}\!\!\ !      l!\prod_{\delta=1}^d\!\!\sum_{\substack{r=0\\ d_\delta+n_\delta-2r = zn}}^{n_\delta }      \!\!\!\!\!\!\!\!\!\!\!\!\!\frac{1}{(n_\delta - r)!r!},\ ] ] i.e. , for nearest - neighbor interactions the product @xmath362_{{\boldsymbol{i}},{\boldsymbol{j}}}[v^{+1/2}]_{{\boldsymbol{i}},{\boldsymbol{j}}}$ ] is always positive for @xmath47 . this does not hold in general and makes it difficult to obtain explicit bounds on the entries of @xmath123 for more general interactions @xmath27 . to obtain a lower bound we keep only the term @xmath363 and @xmath364 . the restriction on @xmath359 is then fulfilled for @xmath365 ( @xmath366 ) if @xmath367 ( @xmath368 ) , yielding @xmath369 it is also possible to obtain an upper bound that is tighter than the one derived for finite - ranged interactions : the elements of @xmath123 are symmetric under @xmath370 , where @xmath371 is a unit vector along dimension @xmath372 . thus , we can demand @xmath373 . then we find that @xmath353 has to be larger or equal to @xmath374 otherwise the restriction on @xmath359 can not be fulfilled . this in turn means that @xmath375 has to be larger or equal to @xmath156 . we then obtain an upper bound by summing all terms in the sum over @xmath359 regardless of the given restriction , yielding @xmath376_{{\boldsymbol{i}},{\boldsymbol{j}}}| & \le & \sum_{l\ge      d({\boldsymbol{i}},{\boldsymbol{j } } ) } a_l^\pm c^l \sum_{\sum_\delta n_\delta = l }      l!\prod_{\delta=1}^d \frac{2^{n_\delta}}{n_\delta ! } \\      & = & \sum_{l\ge d({\boldsymbol{i}},{\boldsymbol{j } } ) } a_l^\pm ( 2cd)^l\\      & = & \sum_{l=0}^\infty a_{l+d({\boldsymbol{i}},{\boldsymbol{j}})}^\pm      ( 2cd)^{l+d({\boldsymbol{i}},{\boldsymbol{j}})}\\      & \le & a_{d({\boldsymbol{i}},{\boldsymbol{j}})}^\pm ( 2cd)^{d({\boldsymbol{i}},{\boldsymbol{j } } ) } \sum_{l=0}^\infty      ( 2cd)^l\\      & = & \frac{a_{d({\boldsymbol{i}},{\boldsymbol{j}})}^\pm}{1 - 2cd } ( 2cd)^{d({\boldsymbol{i}},{\boldsymbol{j}})}.\end{aligned}\ ] ] as in appendix  [ abzaehlerei ] . as before , oscillators belonging to the distinguished ( shaded ) region @xmath5 are marked  , the outside ones are shown as  . @xmath377 is shown as the orange shaded area . @xmath378 for a certain oscillator   is shaded green , its surface oscillators @xmath379 are depicted by  . ] we start by identifying the set of oscillators that can contribute to @xmath380 , @xmath381 , @xmath382 oscillators @xmath383 can only contribute if their distance to the boundary @xmath377 is not larger than @xmath384 , @xmath385 thus , we can restrict the sum over @xmath12 to the set @xmath386 @xmath387 i.e. , we can write @xmath388 where @xmath389 is the number of surface oscillators of a ball with radius @xmath375 within the metric @xmath390 , i.e. , @xmath391 for @xmath392 . using the fact that @xmath393 , @xmath394 can now be bounded from above in the following way @xmath395 where @xmath396 is the volume of a ball with radius @xmath375 within the metric @xmath390 , i.e. , @xmath397 . to summarize , we have @xmath398                g.  vidal , j.  i.  latorre , e.  rico , and a.  kitaev , phys . * 90 * , 227902 ( 2003 ) , http://arxiv.org/abs/quant-ph/0211074[quant-ph/0211074 ] ; j.  i.  latorre , e.  rico , and g.  vidal , quant . comp .  * 4 * , 048 ( 2004 ) , http://arxiv.org/abs/quant-ph/0304098[quant-ph/0304098 ] . k.  zyczkowski , p.  horodecki , a.  sanpera , and m.  lewenstein , phys . a * 58 * , 883 ( 1998 ) ; j.  eisert and m.  b.  plenio , j.  mod . opt .  * 46 * , 145 ( 1999 ) ; j.  eisert ( phd thesis , potsdam , february 2001 ) ; g.  vidal and r.  f.  werner , phys . a * 65 * , 032314 ( 2002 ) ; k.  audenaert , m.  b.  plenio , and j.  eisert , phys . * 90 * , 027901 ( 2003 ) ; m.  b.  plenio , http://arxiv.org/abs/quant-ph/0505071[quant-ph/0505071 ] . note that the above choice of @xmath29 does not restrict generality : clearly , there exists a diagonal matrix @xmath399 , the elements of which are @xmath400 , such that @xmath401 and @xmath402 . hence , there exists a local symplectic transformation @xmath403 relating the hamiltonian with @xmath404 to the one with @xmath29 .
Abstract:

Answer: we demonstrate that the entropy of entanglement and the distillable entanglement of regions with respect to the rest of a general harmonic lattice system in the ground or a thermal state scale at most as the boundary area of the region . this area law is rigorously proven to hold true in non - critical harmonic lattice system of arbitrary spatial dimension , for general finite - ranged harmonic interactions , regions of arbitrary shape and states of nonzero temperature . for nearest - neighbor interactions corresponding to the klein - gordon case  upper and lower bounds to the degree of entanglement can be stated explicitly for arbitrarily shaped regions , generalizing the findings of [ phys . rev . lett .  * 94 * , 060503 ( 2005 ) ] . these higher dimensional analogues of the analysis of block entropies in the one - dimensional case show that under general conditions , one can expect an area law for the entanglement in non - critical harmonic many - body systems . the proofs make use of methods from entanglement theory , as well as of results on matrix functions of block banded matrices . disordered systems are also considered . we moreover construct a class of examples for which the two - point correlation length diverges , yet still an area law can be proven to hold . we finally consider the scaling of classical correlations in a classical harmonic system and relate it to a quantum lattice system with a modified interaction . we briefly comment on a general relationship between criticality and area laws for the entropy of entanglement .