Thermal state entanglement in harmonic lattices
Abstract
We investigate the entanglement properties of thermal states of the harmonic lattice in one, two and three dimensions. We establish the value of the critical temperature for entanglement between neighbouring sites and give physical reasons. Further sites are shown to be entangled only due to boundary effects. Other forms of entanglement are addressed in the second part of the paper by using the energy as witness of entanglement. We close with a comprehensive diagram showing the different phases of entanglement versus complete separability and propose techniques to swap and tune entanglement experimentally.
I Introduction
Entanglement is a purely quantum phenomenon essential for understanding how information is shared in quantum systems and implementing new computing technologies such as measurementbased computation. The existence of entanglement in a composite system implies strong correlations between its parts; so strong that no classical model can explain them. When classical correlations become strong in a condensed matter system, i.e. when they stretch over large distances, a new ‘phase’ emerges which manifests in qualitatively different behaviour of the material. Phase transitions happen when a parameter, such as the temperature, reaches a critical value. In this paper we will determine the critical temperature below which entanglement is present and at which a quantumcorrelated phase emerges in the harmonic lattice.
The entanglement properties of the bisected onedimensional harmonic chain have been analysed in AEPW () using extensive mathematical formalism to discuss the ground state (see also Botero ()). The authors also provide a remark on the situation for thermal states, however, the treatment there is rather short. In the present paper we illuminate the thermal case in greater detail by first discussing twosite entanglement and then any entanglement for finite size chains as well as in the thermodynamic limit, for varying external potential and in one, two and three dimensions. Our intention is to gain intuition about the physics behind the quantumclassical transition at . We explain how the excitation and mixing of the phononic modes degrades the entanglement unless the lattice is in the intrinsically ‘quantum phase’ below where quantum zeropoint fluctuations dominate. We find that entanglement depends on external parameters as well as size and propose that manipulating these will enable techniques of tuning and swapping entanglement in trapped ions.
The outline of the paper is as follows. In Section II the harmonic lattice model is introduced and we summarise technical methods that build the fundament for the results discussed later. Readers familiar with these techniques may skip to Section III where we derive analytical expressions for the entanglement of the onedimensional harmonic chain in the thermodynamic limit and at zero temperature, and then approximate the critical temperature for nearest neighbour entanglement. In Section IV we present comprehensive numerical results for twosite entanglement in the onedimensional chain. The value of the respective critical temperatures is explained and finite chain sizes are discussed. In Section V we investigate the dependence of the nearest neighbour entanglement in the onedimensional chain when increasing an external trapping potential. Finally, this behaviour is compared with the behaviour of entanglement in the two and threedimensional lattice.
In Section VI a different approach to detect entanglement is taken. Using the energy as a witness for any entanglement (not just between two sites) we derive the minimal energy bound for completely separable thermal states. We then establish the critical temperature for witnessed entanglement for the physical sites and find that it is independent of the dimension of the lattice. In Section VII we show and compare the two temperature thresholds derived in the present paper with the threshold for blockentanglement found in AEPW () and the threshold for full separability derived in Winter07 () in a comprehensive entanglement phase diagram. We draw concluding remarks in Section VIII, propose possible uses of our results in controlling and storing information in trapped ion systems and discuss the notion of an entangled phase.
Ii The model
ii.1 Hamiltonian and phononic spectrum
A harmonic lattice can be realised, for example, as a system of ions that are trapped at sites by a harmonic trapping potential . The ions shall couple to their nearest neighbours with equal strength throughout the lattice. The Hamiltonian for the dimensional lattice is
(1) 
Here is the operator for the deviation from the equilibrium position of each ion. The sites are indexed by with being the lattice constant and the components of , running over . Assuming square geometry of the lattice the length in each spatial direction is and is the effective mass of the ions. is the momentum operator corresponding to , so that for all components and the commutation relations hold . The first sum in Eq.(1) is taken over all sites, and the second sum takes all neighbouring pairs of sites, , that are harmonically coupled. To simplify the calculation we will adopt periodic boundary conditions and additionally require to be an odd number for the sake of definitive calculations.
The Hamiltonian can be diagonalised by applying first a Fourier transformation followed by a Bogoliubov transformation into a new set of position and momentum operators, see for instance Brink (). The frequencies of the vibrational modes of the lattice, the phonons, are given by
(2) 
with for a dimensional square lattice with sites, each having neighbours. It will turn out that the frequency spectrum of the phononic modes is key to the value of the critical temperature for nearest neighbour entanglement and the entanglement properties can be modulated by varying the interaction and the trapping potential. (In the following vector arrows indicating the dimension will be dropped and the dimension of the lattice will be stated explicitly when necessary.)
ii.2 Thermal states and covariance matrix
In this work we are interested in the states of thermal equilibrium, , at inverse temperature . For the harmonic lattice these can be written as a tensor product over all phononic modes,
(3) 
where are the number operators for each phononic mode. Since the Hamiltonian is quadratic the thermal states are Gaussian states (see for example contvarstates (); Andersdipl () and references therein) which are specified uniquely by their first and second moments, i.e. the expectation values of the canonical operators themselves and the expectation values of any combination of two canonical operators. The first moments can be changed by choice of the coordinatesystem, i.e. by applying local, single mode displacements, and their absolute value is of no great physical significance. Only their relative values obey conservation laws. Correlations between modes are determined by higher moments and for Gaussian states the second moments are hence vital for the discussion of entanglement. The second moments of a state with modes 1, 2, …. and canonical operators are collected in the covariance matrix, whose elements are given by
(4) 
is a real, symmetric and positive matrix which reveals significant properties of the state such as the occurrence and amount of entanglement in the state. The advantage of formulating the discussion in terms of the covariance matrix is the substantial reduction of parameters, (number of modes), compared to the for the density matrix.
For the harmonic lattice the covariance matrix for two sites and can be found using the transformation from real space into phononic modes in which the Hamiltonian is diagonal. Then expectation values are easily calculated and one obtains
(5) 
where the entries are given by (for a onedimensional chain, similar expressions hold for higher dimensions)
(6) 
which are independent of because of translational symmetry, and
(7) 
which only depend on the distance between the two sites measured in terms of the lattice constant , . Here the inverse temperature is replaced by a unitfree quantity and the phononic frequencies are rescaled by the interaction strength of nearest neighbours .
ii.3 Twosite entanglement condition
For Gaussian continuous variable states of two modes there exists a necessary and sufficient criterion that decides whether the modes are entangled. (The criterion is necessary and sufficient for Gaussian states of modes, Werner:Wolf ().) It is a variation of the positive partial transpose criterion (PPTcriterion) for discrete systems PPTcrit (); HHH96 (). For separability of two sites and in the lattice the criterion requires the positivity of a matrix inequality,
(8) 
where is the symplectic form. is the covariance matrix with timereversed mode ,
(9) 
i.e. for in Eq. (5) we simply replace by to obtain . Using the symmetry of criterion (8) reduces to requiring the positivity of two temperature dependent functions, ,
(10)  
Here is the separation of the two sites and in the onedimensional chain. If positivity fails the two sites are entangled and the magnitude of their entanglement can be given by the logarithmic negativity , defined in Vidal02 (); JensPhD (). can be calculated directly from the functions as
(11) 
Analytical and numerical results based on criterion (8) for two sites are discussed in Section III and Section IV, respectively. The criterion is also sufficient to detect entanglement between any two blocks of sites as shown in AEPW (). Another approach to reveal other kinds of entanglement, i.e. exclude full separability, is the method of entanglement witnesses, introduced in HHH96 (); Terhal02 (). Entanglement witnesses are observables that can take a particular range of expectation values only for entangled state and therefore ‘witness’ entanglement. They have been used to investigate spins systems Brukner (); Toth (); Dowling () and also Wu (), as well as Bosonic gases Anders06 (). In Section VI we use the same strategy for the harmonic lattice to deduce a temperature below which entanglement must be present. A third approach to investigate entanglement in a harmonic lattice uses the full separability condition for the covariance matrix, , for all sites Werner:Wolf (). The mathematical derivation following this direction is given in Winter07 (). Fig. 4 in the present paper unites the results on nearest neighbour entanglement and witnessed entanglement with these previous works and shows how the different strategies complete each other.
Iii Twosite entanglement (1D)
iii.1 Entanglement in the ground state
Any pair of nearest neighbours in the harmonic chain is coupled and we expect it to be entangled in the ground state. Before going in the discussion of the numerical plots in Section IV, let us here derive the negativity Eq. (11) for zero temperature and in the thermodynamical limit.
Let , i.e. , then Further substituting and and taking the continuum limit to replace the sums by integrals one obtains
(12)  
Simplifying this expression further we let and by using the symmetry of the integrands one obtains which can be easily evaluated for arbitrarily distant neighbours. For nearest neighbours we find
(13)  
(14) 
Finding negative implies that nearest neighbours are entangled and this can be quantified in terms of the negativity, Eq. (11),
(15) 
The behaviour of entanglement between any other pair of two sites in the chain is not obvious. Even though nonnearest neighbours do not interact directly, common nearest neighbours could mediate entanglement. We find that this is not the case when the chain size grows arbitrarily, as both and remain positive for . However, for finite chain sizes, the above calculation is only an approximation and the numerical results show small entanglement between next nearest neighbours for very small chains at low temperatures, see Section IV for discussion. To conclude, at zero temperature and in the thermodynamic limit of large twosite entanglement only exists between nearest neighbours and not between any further separated sites. This result agrees with Cramer05 (), where the authors showed that in gapped harmonic lattice systems the correlation functions decay exponentially in the ground state and thermal states.
iii.2 Critical temperature
When the temperature rises, thermal mixing becomes stronger and correlations between the nearest neighbours decrease. We expect intuitively that any entanglement should vanish at large enough temperatures. The transition into the classical regime is marked by a critical temperature which we here want to approximate. The exact temperature for various parameter ranges can be found numerically and is discussed in Sections IV and V.
Again we take the continuum limit and replace the sums by integrals in Eq.(10). Assuming that the trapping potential vanishes, , we obtain The integrals have no closed expression and we need to further approximate the function to evaluate . The coupling parameter , where is the finite interaction strength between neighbours, is temperature dependent. Assuming high temperatures, , one can expand to first order in : . Thus for nearest neighbours one obtains
(16)  
(17) 
Here we have dropped quadratic and higher powers of for consistency. Indeed can become negative and the point where crosses through zero gives the critical temperature for nearest neighbour entanglement, from follows
(18) 
We remark that this temperature is not as large as the initial approximation assumed and one may wonder about its validity. Higher expansions could be considered to evaluate the critical value more accurately. We will not follow this direction here, instead we will present numerical data confirming the derived critical value in the next section.
For the nextnearest neighbours the factors change slightly and the expressions become and . Both expressions are never negative irrespective of the temperature and no entanglement is expected between nextnearest neighbours in the continuum limit. The numerical data in following section independently confirm this result.
Iv Entanglement for finite temperature and size (1D)
iv.1 Data
Numerical data of nearest neighbour (solid) and nextnearest neighbour (dashed) entanglement in terms of the logarithmic negativity , Eq.(11) are displayed in Fig. 1 over the temperature and for varying size of a onedimensional chain. We observe that in general the entanglement is highest for low temperatures and decreases sharply with increasing until it vanishes abruptly at some critical temperature. For nearest neighbours the critical temperature decreases for increasing and converges to in the thermodynamical limit, which agrees with the analytical result in Eq.(18). Hence, the stronger the interaction between nearest neighbours, the more robust their entanglement with respect to thermal noise. The maximal amount of entanglement at decreases for growing and converges to the threshold value of independent of , in good agreement with the analytical result Eq. (15).
Nextnearest neighbour entanglement (dashed curves) is in general much smaller than the nearest neighbour entanglement. Nonetheless, the numerical data show that a small amount of entanglement exists between nextnearest neighbours when the chain is very small. However, the entanglement as well as the critical temperature are considerably lower than for nearest neighbours and decrease rapidly to zero for increasing .
iv.2 Interpretation
We can understand the behaviour of entanglement with changing temperature in the following way. Without quantum zeropoint fluctuations, the ions in the chain would be completely fixed and entanglement would stretch over all sites at zero temperature. However, due to the fluctuation of the phononic modes of the chain (# of phononic modes is with ) correlations are mixed with anticorrelations and the entanglement in the chain reduces considerably. For nearest neighbours the mixing of vibrations is biased. This is a result of the direct coupling between them leading to an average deviation of the ions from equilibrium position of less than the harmonic oscillator length of the interaction. Indeed for small temperatures entanglement exists between nearest neighbours, see the solid curves in Fig. 1.
Nextnearest neighbours and more distant sites do not have a bias of the phononic modes in the thermodynamic limit and entanglement between them is averaged out even at . However, when the chain size is small, finite size effects occur. So few phononic modes are available that the different swinging motions of the ions are not completely balanced and some correlations survive. This is why we find nextnearest neighbour entanglement for very small lattices, see the dashed curves in Fig. 1. It is also the reason why the entanglement between nearest neighbours is in general slightly higher for smaller lattices, as can be seen for the solid set of curves where the uppermost curve displays the behaviour of a tiny chain with , and the lower curves apply for increasing .
When the temperature is raised, more phonons are excited and the mixing of different vibrations is enhanced. This leads to a rapid decay of the entanglement and when the temperature reaches the magnitude of the phononic modes that drive nearest neighbours against each other the entanglement vanishes completely. This happens when the two degrees of freedom, the position and the momentum of each of the ions, with classical thermal energy per degree of freedom exceeds the energy of the quantum fluctuations of each of the springs that connect nearest neighbours, (for ). We thus expect nearest neighbours to behave classically when is fulfilled. This heuristic argument gives the critical temperature for entanglement, , in agreement with the numerical cutting point in Fig. 1. It also indicates that the energy is a good witness of entanglement, which we will elaborate in Section VI.
V Dependency on frequency spectrum and dimension of space
v.1 Data
The impact of the phononic frequency spectrum on the entanglement properties are studied quantitatively by varying the ratio between trapping potential and nearest neighbour coupling. The numerical results are displayed in Fig. 2. As before, for small the entanglement starts at for and decreases for increasing temperature until it vanishes at . However, when the trapping potential is raised and reaches the magnitude of the interaction the initial entanglement starts off much lower, for . Notable is that this entanglement is more robust and survives increasingly high temperature mixing. For example, when , the critical temperature has increased to .
v.2 Interpretation
When the trapping potential is raised close to and above two effects emerge. Firstly, the onsite potential fixes the ions tightly to their sites and only a small fraction of the available energy can be put into the interaction. Even though they are coupled, the ions cannot follow the motion of their neighbours and entanglement between nearest neighbours starts off at substantially lower values for compared to the case of small , see Fig. 2. The second effect is that higher phononic modes ‘cost’ a higher amount of energy. Thus for growing the bias of phononic modes cannot be compensated even at higher temperatures and nearest neighbours stay entangled up to increasing critical temperatures, as shown in Fig. 2.
For higher dimensions, that is, when the chain becomes a two or threedimensional square lattice, the behaviour of twosite entanglement is qualitatively similar to the onedimensional case, as shown in Fig. 3. Comparing this figure with the graphs for varying in Fig. 2 one can see that raising the dimension leads to similar effects as raising the trapping potential. The higher the dimension the lower is the ground state entanglement, but the higher rises the critical temperature.
The reason behind this behaviour is the growing fixation of the ions to their sites for increasing dimension. Namely, in the thermodynamic limit (and for ) the average displacement of the ions from their equilibrium position, , diverges linearly in 1D and logarithmically in 2D, while converging in 3D. No longrange order can be established in 1D and 2D as the crystal simply ‘melts’ (for a discussion see for example Brink () and references therein). Only in three dimensions do the fluctuations of the ionic positions become smaller than the lattice spacing, longrange order can occur and a solid with a threedimensional crystalline structure can form. To obtain the same effect one could impose a strong trapping potential to hold the ions tightly in place and thereby force crystallisation in low dimensions. Thus, increasing the trapping potential is effectively raising the dimension of the lattice structure and vice versa.
Vi Energy as an entanglement witness in dimensions
Twosite entanglement is a very restricted kind of entanglement. Other kinds of entanglement exist, for example genuine threepartite entanglement. As indicated by the previous numerical results the balance between the available thermal energy compared to the phononic spectrum determines the existence of entanglement implying that the energy may be a useful criterion as an entanglement witness. This idea has been used in a number of works before, for example in Brukner (); Dowling (); Toth (); Anders06 (). Adopting this strategy for the harmonic lattice in dimensions we identify the temperature range where some kind of entanglement must necessarily be present.
vi.1 Derivation of energy bound for separated sites.
The most general, fully separable (into sites) state is of the form
(19) 
where are probabilities with . The energy (c.f. Eq. (1)) can be evaluated for all separable states as
(20) 
The second term can w.l.o.g. be set zero since all first order expectation values of any operator can be shifted locally to zero without affecting the nonlocal entanglement properties. Thus, the energy of the whole lattice in any separable configuration, , becomes the sum of the individual energies for all sites,
(21) 
where is an arbitrary state for site and the effective Hamilton operator for one site alone, where the and are in general dimensional operators. Thus the description of the lattice has become a meanfield model where each site interacts with an effective, identical background produced by the interaction with all other sites. Reversely, a description using only effective onesite Hamiltonians leads to a classical model where entanglement has been averaged out and states are always separable.
can be rewritten as a harmonic oscillator, , where the frequency is identical for all sites . By the presence of quantum zeropoint fluctuations we always have for any state . The possible energy for a separable configuration of the whole lattice is thus bounded from below by
(22) 
Any state having an energy below this bound is necessarily entangled with respect to the sites.
vi.2 Entanglement for different modes.
It may seem puzzling that one could apply the same strategy for other set of modes, for instance, for the phononic modes themselves and find a nonzero energy bound. Indeed, for the normal modes one would find
(23) 
with as given in Eq. (2) for dimensions. This inequality is always fulfilled as the right hand side gives just the ground state energy. The criterion is thus useless here as it can never be violated and is hence unable to detect any entanglement.
By symmetry we expect that the sites in space are ‘most’ entangled, i.e. give the highest energy bound, while the phononic modes remain separable in any thermal state. The situation for ‘intermediate modes’ is between these extremes. To show this we use again the energy witness argument (assuming one dimension for simplicity). For all orthogonal transformations that transform the set of phononic momentum operators, , into some new momentum operators , such that and similarly for the position operators, , so that the commutation relations remain , the Hamiltonian becomes
(24)  
where the interactions are given by the elements of the matrix . With the same argument as before we obtain a minimal energy threshold for all separable states with respect to the new set of modes,
(25) 
In particular, the diagonal frequencies are with the doublystochastic matrix , . With the Birkhoffvon Neumann theorem can be written as a convex combination of permutation matrices of the same order and because the squareroot is concave we only need to consider the two extremal cases which bound from below and above. The first is a single permutation with resulting in a reordering , and the second is the complete mixture of all possible permutations such that and . These extremes are indeed identical to the phononic modes and the sites, respectively. The energy bound for the general set of modes is thus bounded from below and above,
(26) 
being bigger than the zeropoint energy, while remaining smaller than the energy bound for individual sites. Site entanglement thus possesses the highest energy bound.
vi.3 Critical temperature for separated sites.
The energy of the dimensional harmonic lattice in a thermal state can be related to the temperature. The internal energy, , can be explicitly calculated and for high temperatures one finds the Boltzmann’s equipartition law, where each degree of freedom of the system obtains a thermal energy contribution of ,
(27) 
(In this context the relation is identical to the DulongPetit law of constant heat capacity at high temperatures.) Setting the internal energy equal to the energy bound for the sites, Eq. (22), one finds that any fully separable configuration has a temperature above a certain critical value
(28) 
Conversely, any configuration with a lower temperature must necessarily contain entanglement between the sites of the lattice. Remarkably, the critical value is independent of the dimension of the system. Furthermore, the witness temperature for , is close to the nearest neighbour critical temperature in 1D. In higher dimensions the nearest neighbour critical temperature increases and approaches the witnesstemperature, as shown in Fig. 3.
The energy relation Eq. (27) and hence the witness temperature Eq. (28) is initially only valid in the asymptotic limit of high temperatures. For small temperatures the behaviour of the internal energy and its derivative, the heat capacity , has been subject of intense discussion at the beginning of the 20th century (which is covered in many textbooks, see e.g. heatcapacitydebate ()). Introducing a heuristic cutoff frequency , Debye found that in three dimensions the phononic oscillations lead to a heat capacity proportional to , implying for temperatures below . With the same method the internal energy scales with in lower dimensions. Discussing the Debye temperature in detail is beyond the scope of this paper, however, direct conclusions can be drawn just from the properties of the energy function. Because of the monotonicity of the internal energy and its convexity at small temperatures, is below the value of the hightemperature limit, in Eq. (27), in the regime . The witness temperature derived with the high temperature approximation, , is thus a lower bound to the critical temperature that one would obtain when using the low temperature formula. The entangled region detected with the high temperature approximation is thus a subset of the region which one could detect when using the small temperature approximation. This implies that is a valid bound for witnessed entanglement also at small temperatures, being most tight in one dimension and lesser so in higher dimension due to the stronger bend of the law.
Vii Entanglement phase diagram
Fig. 4 is a comprehensive and exact diagram comparing the two derived temperature thresholds for nearest neighbour entanglement and witnessed entanglement with previous results from AEPW () and Winter07 (). As found in Section IV, twosite entanglement resides between nearest neighbours (dark shaded) for temperatures below (for , ) while nonnearest neighbour entanglement vanishes for large . The hatched set are states whose nonseparability is detected by the energy witness. It fully includes all nearest neighbour entanglement and extends further up to the temperature , cf. Eq. (28), where no twosite entanglement exists. A different phase of how entanglement is distributed is thus reached.
For example, entanglement can exist between two neighbouring blocks of sites. Results in AEPW () for (this appears to be the critical value according to Figure 12 in AEPW () after retaining all units) confirm this and their critical temperature for blockentanglement, , is indicated in Fig. 4 by the dark point. When the temperature is raised further the chain passes over into complete separability (light shaded), i.e. no entanglement exists for any possible split. The temperature at which this transition happens is given as (for ) in Winter07 (). For sufficiently small , the witness bound is in good agreement with the exact threshold for full separability, Winter07 (). This indicates that the temperature range where the entanglement survives thermal mixing is identical to the regime where the quantum zeropoint effects are important. The use of the energy as an entanglement witness is thus justified and in general an alternative, easier way to evaluate regions where entanglement is present, c.f. Brukner (); Toth (); Dowling (); Anders06 (). For higher dimension we find that the critical temperature for nearest neighbour entanglement increases towards the dimension independent witness temperature (see Fig. 3) thus leaving less temperature range where nontwosite entanglement is witnessed.
When is large compared to , the frequency spectrum shrinks to essentially one dominating frequency, . The mixing of frequencies occurring at increasing temperatures becomes irrelevant and entanglement persist up to higher temperatures, the higher becomes. The asymptotic behaviour of the witnesstemperature becomes proportional to as . Yet the fullseparability temperature scales with as which implies a growing gap between the witnessed and the true entangled region for growing .
Interestingly in this model any entanglement must be nonPPT Cavalcanti07 (), in other words, there always exists a bipartite split of the chain into two parties such, that partial transposition on one of the parties results in a nonpositive density matrix. Any entanglement can thus potentially be identified by this property and the value of the negative eigenvalue can be used to measure the amount of entanglement.
Viii Conclusion and Outlook
In condensed matter physics the twopoint correlation function correlationfunction () between two sites is an important quantity used to determine the behaviour of the material. If correlations stretch over a wide range of sites, i.e. the correlation function becomes a constant or decreases polynomially for increasing distances, then one speaks of the establishment of an ordered phase. Quantum physics allows for a stronger kind of correlation, that is, entanglement. The entanglement criterion Eq. (10) is reminiscent of the classical order criterion for the twopoint correlation function. Instead of requiring that the fluctuations of the ionic positions, , should be smaller than the lattice spacing, we here have a condition on the product of the relative positions and the relative momentum fluctuations. In other words, instead of requiring order in space the entanglement criterion (10) requires order between two sites in phase space.
We have found that entanglement exists between nearest neighbours in all dimensions, for all lattice sizes, and for all values of , when the system is cooled below some critical temperature, which depends on the former parameters. By changing the temperature the entanglement can be modulated and even cut out. The critical temperature at which the entanglement vanishes is directly proportional to the coupling constant; for the harmonic chain this is the coupling . This relation is in agreement with previous results in spin chains Brukner () and Bose gases Anders06 (); Heaney07 () where also the interaction strength determines the critical temperature.
Entanglement between nextnearest neighbours exists in very small lattices, but is in general much weaker in magnitude and also only observable at much lower temperatures than for the nearest neighbour case. A possible application for trapped ions in lattices with harmonic interaction, would be to tune the value of and such, that only nearest neighbour ions remain entangled. By slowly inserting an additional ion in the lattice one can destroy entanglement to the former nearest neighbours and the new nearest neighbours become entangled. Varying the trapping potential or the interaction strength , one can further change the amount and range of entanglement in the lattice. Such techniques and the possibility to apply additional external potentials that modify the interaction and hence affect the structure of the entanglement may allow to store and manipulate information on the ions. For instance, applying an extra harmonic potential to all ions in a small periodic chain results in a higher frequency of the lowest phononic mode, the centre of mass mode, that could influence the chain to ‘swap’ entanglement from nearest neighbours to opposite ions in the ring by simply raising the temperature.
Even though in the thermodynamic limit only neighbouring pairs are entangled and more distant pairs are not, we can speak of an ‘entangled phase’. This notion is justified because entanglement is much stronger than classical correlations and nearest neighbour entanglement is already sufficient to transport information across the whole lattice. Such correlations enable universal quantum computation in the measurementbased model, as initiated in Raussendorf01 () for spin systems with nearest neighbour interaction and its analogue for the continuous variable case proposed in Menicucci06 (). The entangled phase survives, approximately, as long as the system is in the quantum regime where the zeropoint fluctuations dominate the physics of the lattice. Going beyond that range implies the validity of a meanfield theory where each site can be described as an independent oscillator whose effective frequency includes the contributions from the rest of the system.
Acknowledgements.
J.A. is supported by the Gottlieb Daimler und Karl BenzStiftung. This work was supported in part by the Singapore A*STAR Temasek Grant No. 0121040040 and this research is part of QIP IRC www.qipirc.org (GR/S82176/01). I wish to thank A. Ekert, B.G. Englert, L. C. Kwek, Ch. Miniatura, J. Suzuki, V. Vedral, A. Winter for many insightful discussions on the subject.Footnotes
 email: janet@qipc.org
References
 K. Audenaert, J. Eisert, M. B. Plenio and R. F. Werner, Phys. Rev. A 66, 042327 (2002).
 A. Botero, and B. Reznik, Phys. Rev. A 70, 052329 (2004).
 J. Anders, and A. Winter, Quant. Inf. & Comp. 8 (3&4), 0245 (2008).

van den Brink, J.,
http://www.ilorentz.org/~brink/course/tcm.html.  Quantum Information with continuous variables of atoms and light, Editors: N. J. Cerf, G. Leuchs, and E. S. Polzik, Imperial College Press (2007).
 J. Anders, Diploma Thesis, Universität Potsdam, arXiv:quantph/0610263, (2003).
 R. F. Werner, and M. M. Wolf, Phys. Rev. Lett. 86, 3658 (2001).
 A. Peres, Phys. Rev. Lett. 77, 1413 (1996).
 M. Horodecki, P. Horodecki, and R. Horodecki, Phys. Lett. A 223, 1 (1996).
 G. Vidal, and R. F. Werner, Phys. Rev. A 65, 032314 (2002).
 J. Eisert, PhD Thesis, Universität Potsdam, arXiv:quantph/0610253, (2001).
 B. M. Terhal, J. Th. Comp. Sc. 287(1), 313 (2002).
 Č. Brukner, and V. Vedral, 2004, arXiv:quantph/0406040.
 G. Toth, Phys. Rev. A 71, 010301(R) (2005).
 M. R. Dowling, A. C. Doherty, S. D. Bartlett, Phys. Rev. A 70, 062113 (2004).
 L. A. Wu, S. Bandyopadhyay, M.S. Sarandy, and D.A. Lidar, Phys. Rev. A 72, 032309 (2005).
 J. Anders, D. Kaszlikowski, Ch. Lunkes, T. Oshima, and V. Vedral, New J. Phys. 8, 140 (2006).
 M. Cramer, and J. Eisert, New J. Phys. 8, 71 (2006).
 G. H. Wannier, Statistical Physics, John Wiley & Sons, Inc., New York (1966).
 D. Cavalcanti, A. Ferraro, A. GarciaSaez, and A. Acin, Phys. Rev. Lett. 100, 080502 (2008).
 The twopoint correlation function measures the correlation between the orderparameter at two positions separated by . A solid is characterised by longrange order of the mass density . For the harmonic lattice the twopoint correlation function can be expressed as , which can be calculated explicitly, see for example Brink (); Chaiken ().
 P. M. Chaiken, and T. C. Lubensky, Principles of condensed Matter Physics, Cambridge University Press, Cambridge (UK), (1995).
 L. Heaney, J. Anders, D. Kaszlikowski and V. Vedral, Phys. Rev. A 76, 053605 (2007).
 R. Raussendorf and H.J. Briegel, Phys. Rev. Lett. 86, 5188 (2001).
 N. C. Menicucci, P. van Loock, M. Gu, Ch. Weedbrook, T. C. Ralph, M. A. Nielsen, Phys. Rev. Lett. 97, 110501 (2006).