A scheme for entanglement extraction from a solid

Some thermodynamical properties of solids, such as heat capacity and magnetic susceptibility, have recently been shown to be linked to the amount of entanglement in a solid. However this entanglement may appear a mere mathematical artifact of the typical symmetrization procedure of many-body wave function in solid state physics. Here we show that this entanglement is physical demonstrating the principles of its extraction from a typical solid state system by scattering two particles off the system. Moreover we show how to simulate this process using present-day optical lattices technology. This demonstrates not only that entanglement exists in solids but also that it can be used for quantum information processing or for test of Bell's inequalities.

Abstract. Some thermodynamical properties of solids, such as heat capacity and magnetic susceptibility, have recently been shown to be linked to the amount of entanglement in a solid. However this entanglement may appear a mere mathematical artifact of the typical symmetrization procedure of many-body wave function in solid state physics. Here we show that this entanglement is physical demonstrating the principles of its extraction from a typical solid state system by scattering two particles off the system. Moreover we show how to simulate this process using present-day optical lattices technology. This demonstrates not only that entanglement exists in solids but also that it can be used for quantum information processing or for test of Bell's inequalities. Entanglement [1] is the property of quantum systems to show correlations that cannot be explained by classical mechanics [2]. In the last two decades entanglement has been recognized as a resource in quantum information processing, like quantum cryptography [3] and quantum computation, [4,5]. Furthermore in recent years the role that entanglement plays in the properties of condensed matter systems has received an increasing attention. For instance in a typical magnetic solid the electron orbitals of different atoms overlap giving rise to an exchange interaction between the electron spins which can be described by a Heisenberg Hamiltonian, whose thermal states are entangled [6].
It is a common belief that entanglement cannot exist on a macroscopic scale since decoherence effects would destroy all quantum correlations. However it has been predicted that macroscopic entanglement not only is related to critical phenomena [7] but also that it can exist in solids in the thermodynamic limit [8] even at high temperatures [9]. Indeed there are many works inferring the existence of macroscopic entanglement at various temperatures up to room temperature [10][11][12]. The evidence for this entanglement comes from experiments measuring different thermodynamic properties such as magnetic susceptibility and heat capacity. The main aim of this article is to show that such entanglement is not an artificial mathematical property but that it can be extracted and therefore used for quantum information processing in the same way as heat can be extracted and used for work in thermodynamics. In particular we will show how such entanglement can be transferred to a pair of particles and subsequently used, in principle, for quantum communication or to test the violation Bell's inequalities. Furthermore we propose a scheme to demonstrate entanglement extraction with presentday technology using optical lattices. In the next two paragraphs we outline the basic idea behind our proposal which is then elaborated in the remaining part of the article.
In the ideal scenario the entanglement between a pair of spins inside a solid can clearly be transferred to another pair of probe particles using only local swap operations. To this end one sends simultaneously a pair of probe particles toward the entangled spin chain in such a way that each probe interacts with a different spin (cf. Fig.1). We emphasize that the two probes do not interact with each other nor do they experience an interaction mediated by the solid (cf. [13]). The entanglement between the probes has been extracted from the spin chain and cannot exist without entanglement in the chain. This is a genuine non local process between the two probes like in the case of entanglement swapping [14]. We mention that in the continuum limit our procedure can be used to extract entanglement from vacuum (some steps in this direction have been taken by [15]). In practice, however, the spin chain is in a mixed entangled state and the scattering interaction between probes and spins in the chain can only partially swap their state. Is entanglement extraction still possible under such realistic conditions? In this article we answer affirmatively to this question. We will consider only common physical interactions like the Heisenberg or the XY ones. We will show that in general, after one scattering event, entanglement swapping is only partial. In this case it is interesting to investigate whether by repeating this collision procedure with the same probes, one can achieve a better entanglement extraction.
The most natural way to extract entanglement from entangled electron spins in solids would be to scatter pairs of neutrons off the solid. We will present a proposal for simulation of this process of entanglement extraction with optical lattices [16] that we believe can be implemented at present. Optical lattices are a very useful tool for simulating many quantum effects in solids and have been employed in demonstrating quantum phase transition [16], and quantum gates [17]. For our purposes we need to simulate both the interaction in the solid system as well as the interaction with the external probes. Hamiltonians of entangled spin chains or ladders can be realized using cold neutral atoms trapped in potential wells generated by counter-propagating lasers [18][19][20]. Bosonic atoms (black in Fig. 2) are loaded onto the lattice in such a way that there is only one atom per lattice site and that the hopping between neighboring sites is inhibited. Long-lived atomic states are used to simulate the electronic spin in a typical solid. The interaction between atoms can be varied by adjusting the laser parameters and by means of electric and magnetic fields. In this way one can produce Figure 1. Schematic drawing of the entanglement extraction process. Two probe particles initially prepared in a product state are sent toward a solid whose spins are entangled. Each probe interacts with only one spin in the solid. During the collision, due to the exchange interaction between probe and spin, their state is swapped. The two probes come out from the solid in an entangled state. a variety of common spin Hamiltonians. As probe particles one can use different atoms Figure 2. Simulation of the entanglement extraction using optical lattices. Cold neutral atoms (black) are loaded onto an optical lattice generated by counterpropagating lasers in such a way that each potential well is occupied by only one atom. These atoms represent the spins degrees of freedom in a typical solid. To the first potential a second one is superimposed (wireframe) onto which different atoms (white), called markers, are loaded. It is possible to move the markers across the solid by moving adiabatically the second potential. This can be efficiently done by varying laser parameters of the second lattice. A collision between a marker and a spin can be produced by moving the marker into the same potential well of the spin.
(white in Fig. 2), called marker qubits [21]. The main idea is to construct first the optical lattice, loading the atoms in a register that constitutes the simulated solid and then to superimpose a second optical lattice loading auxiliary atoms onto it as shown in Fig. 2. The auxiliary marker atoms need not to be of the same atomic species as the atoms in the solid. The marker atoms can be moved through the solid by varying adiabatically laser parameters of the second optical lattice. These operations do not affect the register atoms which remain confined in their potential wells. In this way collisions between probes and spins in the solid can be simulated by moving a marker atom close to a register atom letting their orbitals overlap as shown in Fig. 2. This in principle can be performed with high precision and efficiency as in the case of cold controlled collisions [17]. In this way one can realize the swap operations needed in our protocol. Such an operation can be repeated several, times moving the markers forward to collide with other register atoms.
Let us now introduce the detailed description of the interaction between the probes and the spins in the solid. We stress again that all aspects of the following analysis can be implemented not only with optical lattices described above but also with any other quantum systems capable of coherent manipulation. Two probes are sent to interact for a time τ each with a single spin of an entangled pair in a solid as shown in Fig. 1. We suppose that the interaction Hamiltonian is given by an anisotropic XXZ model between probe i = L, R and spin j = 1, 2: where σ x,y,z are the Pauli matrices, J is the coupling which we assume to be constant during the collision. In the case of a time dependent J(t) our analysis still applies by making the substitution: The total Hamiltonian is therefore 2R . This Hamiltonian comes from the same exchange mechanism that gives rise to the interaction between spins in the solid. This Hamiltonian can be easily implemented with optical lattices and marker qubits in a probabilistic fashion. Here we concentrate on two important limits: λ = 1, which is the Heisenberg (exchange) interaction; λ = 0, which is the XY interaction. We assume that the spin chain is initially in the ground state, however our discussion can also be extended to other outof-equilibrium states. We also assume that during the only dynamics is due to the interaction with the probes. The reduced density matrix of the entangled pair can be written in terms of spin-spin correlation functions by using translational invariance and conservation of angular momentum. The full derivation is presented in [22] and the result is: where g zz = σ 1 z σ 2 z /4 and g xx = σ 1 x σ 2 x /4 are spin-spin correlation functions and can be numerically estimated for the Heisenberg model. The concurrence of ρ, assuming |g xx | ≥ g zz , is: Note that the results given below apply to a large class of mixed entangled states including, for example, Werner states. Let us first consider what happens after a single collision between the probes in a pure product state e.g. |00 and the spins when λ = 1, 0. To illustrate the idea let us suppose that the spins are initially in a maximally entangled state |ψ − = 2 −1/2 (|01 − |10 ). After the collision the global state of the four particle reads (apart from a global phase factor): This state is a superposition of the initial state and of the swapped state. By tracing out the two spins one obtains the reduced density matrix ρ LR of the two probes. We use the concurrence [23] to measure the entanglement extracted. For a density matrix of two qubits ρ let us defineρ . = σ y ⊗ σ y ρ * σ y ⊗ σ y and R = ρρ. The concurrence is defined as C = max{0, λ 1 − λ 2 − λ 3 − λ 4 } where λ i are the square roots of the eigenvalues of R labeled in decreasing order. The concurrence for the two probes is simply C = sin 2 2Jτ and is equal to 1 when Jτ = π/4. Thus no matter how small τ is, it is still possible to extract some entanglement from the singlet state.
The extracted entanglement (i.e. the concurrence between the probes), maximized over the set of initial product states of the probes, is plotted in Fig. 3, for λ = 1, as a function of g xx = g zz , including therefore the case of mixed states of spins in the solid. Note that ρ is entangled for −1/4 ≤ g zz < −1/12 . As it is shown in the plot the entanglement between probes decays when g zz decreases. This is confirms the intuitive expectation that the more the spins in the solid are entangled, the more such entanglement can be transferred to the probes.
We want to stress that for a translational invariant system the correlation functions can never reach the maximum value g xx = g zz = −1/4. Pure maximal entangled states can not exist in a homogeneous spin chain. Nevertheless in Fig. 3 we show the concurrence for all possible values of g zz .
The state of the spins and that of the probes are swapped when Jτ = π/4 and the entanglement extracted is maximum. In an optical lattice implementation this can be achieved in a time comparable to the trapping frequency of the atoms (10µs for 87 Rb, [21]). In the case of solid J is of the order of 1K giving approximately τ ∼ 10 −11 s [11]. We also note that entanglement oscillates with τ , a feature which reflects the fidelity of the swapping operation. The probes are always entangled for g zz −0.16 no matter how small τ is. However, as shown in Fig. 3, for −0.16 < g zz < −1/12, τ must be chosen appropriately in order to extract some entanglement. Notice that obviously no entanglement can be extracted for g zz ≥ −1/12 since the two spins are not entangled. The important point is that as long as some entanglement is shared between two spins in the solid (i.e. for g zz < −1/12) our procedure is capable of transferring part of it to the probes. Similar results applies also to the XY model λ = 0.
Once τ is fixed for a specific experimental setup, the amount of entanglement extracted is less than the entanglement in the pair of spins as long as Jτ < π/4. Is it possible to extract more entanglement by repeating the collision process several times? . Extracted entanglement for an Heisenberg chain, measured by the concurrence C between the two probes, maximized over the set of initial product states of the probes, after one collision as a function of g zz and τ . Notice how C increases with the correlation g zz and oscillates with the interaction time τ .
This can be done in optical lattices by moving the marker qubits so they interact with different identically prepared chains (cf. Fig. 2). Although this approach has some conceptual analogies with the so called homogenization process [24], the two schemes are not equivalent since the transformation induced by the interaction Hamiltonian (1) is not a global partial SWAP transformation. Indeed with a local interaction generated by the spin Hamiltonian (λ = 1) the full transformation of the four spins is: swaps two qubits. It is easy to demonstrate that, for λ = 0 (XY model), the state |ψ + LR |ψ − 12 is an eigenstate of H T . It is thus a fixed point of the evolution transformation U 12LR and the state of the probes approaches the state |ψ + as the number of collisions goes to infinity, independently of the input state and of how small τ is. This is a counter-intuitive result: the state of the probes do not converge to the state of the bath but to an orthogonal state. The reason is that XY interaction does not generate a partial SWAP operation but something similar to a partial iSWAP operation: that swaps the state of the two qubits but with a i factor. The composition of two of these transformation produces a minus sign that transform |ψ + in |ψ − . In Fig. 4 we show how the concurrence grows increasing the number of collision with the spins. For λ = 1 (Heisenberg model) there is not an eigenstate of the total Hamiltonian which is the tensor product of a maximally entangled state of the probes times a maximally entangled state of the two spins.
When the state of the two spins in the solid is mixed, the transformation has still a unique fixed point [25,26]. This follows form the Banach-Cacciopoli theorem for contractive maps. The transformation we are considering is the tensor product of two partial swap transformation which are contractive maps and thus it is itself a contractive map. The entanglement of the fixed point depends on g zz , g xx and τ . For example, for a Heisenberg interaction, when g zz = −0.2 and Jτ = 0.6, the concurrence of the fixed point is C ≃ 0.28. However, when the two spins in the chain are in a mixed state, the entanglement extracted by the probes after repeated collisions is always less than that after only one collision for λ = 1 (Heisenberg interaction) and "nearly" always, i.e. except for mixed states which are very close to a pure state, for λ = 0 (XY interaction). Therefore repeated collisions do not increase the extracted entanglement.
Under different hypothesis it is possible by repeated collisions to homogenize the state of the probes to the state of the spins in the solid. By modifying the interaction between the four spins 1, 2, L, R we can build a Hamiltonian that generates a global SWAP transformation: P SW 1L,2R = e iJt (cos 2Jt 1l 1L ⊗ 1l 2R − i sin 2Jt SW AP 1L ⊗ SW AP 2R ) (9) This can be found by inspection and reads in the computational basis of 12LR: Notice that the elements of the first two rows can be generated by pairwise interactions as in (1). The elements of the third row instead are four spin processes. With this interaction it is proven by the homogeneization process [24] that the probes state will approach the state of the spins no matter the initial state. Until now we concentrated on extracting bipartite entanglement, however our scheme can also be adapted to extract genuine multipartite entanglement. Let us consider n spins in a chain, each scattering a spin probe. Let the initial state of the chain be a W state: √ n (|100 · · · 0 + |010 · · · 0 + · · · + |000 · · · 1 ) (11) and the probes be in state |0 = |00 · · · 0 . W states can be ground states of some Hubbard related Hamiltonians [27], critical spin chains [7,28] and fermionic lattice models of high T c superconductivity [9]. The state of the system at time τ is the analogue for n spins of Eq. (4) with |W n instead of |ψ − . Thus the reduced density matrix for the n probes reads: Notice that when Jτ = π/4 the state of the chain is fully swapped onto that of the probes as in the case with n = 2. For any value of Jτ = kπ/2, where k is an integer, it can be shown that the above state contains genuine multipartite entanglement [28] (i.e. it cannot be written as a mixture of biseparable state). So far we have assumed that each probe interacts with a single spin in the chain. This assumption is, in some cases, unrealistic. For instance it is difficult, with present day technology, to address a single spin in a solid. For example, in neutron scattering experiments, the spin of the incoming neutron interacts with the total angular momentum of a bunch of spins in the solid [29]. It is thus more realistic to analyze a model in which two probes interact with several spins at once. Let us consider a chain of L spins and let us assume that each neutron, being a wavepacket of a certain width, interacts with a different subset of N spins. We will assume that the spin of the probe is coupled to the total angular momentum of the N spins, which is equivalent to assume that each probe is equally coupled to each of the N spins in the chain. This model is equivalent to the ones considered in others contexts [30], and for λ = 0, 1 its eigenstates are known. Let us consider a W L state as the initial state of the chain. Though this is not the most general entangled state of a spin chain, this model can be solved analytically. The concurrence for the two probes for λ = 0 is C = 2N L sin 2 2J √ N τ and reaches the maximum value 2N L for Jτ = π/(4 √ N) (the entanglement for the Heisenberg model λ = 1 is greater than zero but is always less than that of the XY model). Notice that 2N/L is just the fraction of spins interacting with the probes and thus it is the maximum bipartite entanglement that can be extracted. We note that these results have been obtained for a W state of the chain. This state is symmetric under permutations of any number of spins which means that the entanglement is equally distributed among all spins. However for Heisenberg and XY chain models long range correlations are very small and the entanglement extracted should decrease with N and with the distance between the two beams of neutrons. From the above analysis it emerges that a necessary condition to extract entanglement within this realistic model is that segments of the spin chains should be entangled. This analysis can be extended the case of non-equilibrium states and it would be interesting to see if the entanglement extracted could be larger than in the case of the W state.
We have described a scheme to extract entanglement from spins in a solid by scattering probe particles. In this way macroscopic thermal entanglement of solids can be converted into a useful resource for quantum information processing. To illustrate the scheme we have proposed an optical lattice implementation where all aspects of our procedure can be realized with present-day technology. In this context optical lattices offer only a simulation of entanglement extraction. This is why we have discussed how neutron scattering can be used to achieve entanglement extraction from a real solid as well as the physical limitations of this process. We hope that our ideas in the long run will lead to an entangling procedure for neutrons in the same way as parametric down conversion is for creating entangled photons.