- About this Journal
- Abstracting and Indexing
- Aims and Scope
- Article Processing Charges
- Articles in Press
- Author Guidelines
- Bibliographic Information
- Citations to this Journal
- Contact Information
- Editorial Board
- Editorial Workflow
- Free eTOC Alerts
- Publication Ethics
- Reviewers Acknowledgment
- Submit a Manuscript
- Subscription Information
- Table of Contents
Advances in Mathematical Physics
Volume 2012 (2012), Article ID 375182, 15 pages
Creation of Two-Particle Entanglement in Open Macroscopic Quantum Systems
1Department of Mathematics and Statistics, Memorial University of Newfoundland, St. John’s, NL, Canada A1C 5S7
2Theoretical Division, Los Alamos National Laboratory, MS B213, Los Alamos, NM 87545, USA
3Dipartimento di Matematica e Fisica, and Interdisciplinary Laboratories for Advanced Materials Physics, Università Cattolica, Via Musei 41, 25121 Brescia, Italy
4Istituto Nazionale di Fisica Nucleare, Sezione di Pavia, Via Bassi 6, 27100 Pavia, Italy
5Department of Applied Physics, Polytechnic Institute of NYU, 6 MetroTech Center, Brooklyn, NY 11201, USA
Received 6 March 2012; Accepted 17 May 2012
Academic Editor: P. G. Kevrekidis
Copyright © 2012 M. Merkli et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
We consider an open quantum system of not directly interacting spins (qubits) in contact with both local and collective thermal environments. The qubit-environment interactions are energy conserving. We trace out the variables of the thermal environments and qubits to obtain the time-dependent reduced density matrix for two arbitrary qubits. We numerically simulate the reduced dynamics and the creation of entanglement (concurrence) as a function of the parameters of the thermal environments and the number of qubits, . Our results demonstrate that the two-qubit entanglement generally decreases as increases. We show analytically that, in the limit , no entanglement can be created. This indicates that collective thermal environments cannot create two-qubit entanglement when many qubits are located within a region of the size of the environment coherence length. We discuss possible relevance of our consideration to recent quantum information devices and biosystems.
In open many-body systems, such as solid-state and biological ones, quantum behavior reveals itself in many ways. Often the quantitative parameter used to measure “quantumness” (possibly even of macroscopic order) is entanglement. The presence of entanglement implies that the wave function (or the reduced density matrix) cannot be represented as a product of the corresponding objects for the individual qubits. It is important to note that to produce and to measure entanglement in such systems, one does not necessarily need to know much detail about the system, possibly not even its Hamiltonian . The questions then are how useful is entanglement as a measure of quatumness, what can it add to our knowledge of system properties and behavior, and how can it be utilized? Indeed, just knowing that the system is entangled (knowledge of complicated quantum behavior) is not sufficient to imply that its quantum properties are useful for specific applications. Fortunately, however, in some situations entanglement provides very useful properties, including additional exponential resources for quantum computation  and a possible enhancement of photosynthesis in bio-systems  (See also references therein.)
Entanglement could be produced by direct interaction between the qubits. This interaction should be of a “conditional” nature, mixing an initial product state (disentangled) in such a way that the final state becomes correlated in a quantum way. Many aspects of entanglement creation are widely discussed in the literature (See, e.g., [4–11] and references therein.)
More recently, interest appeared in the possibility to create entanglement in the absence of direct interactions between qubits (or when the latter are very small). Entanglement can then still be created merely by the indirect interaction of noninteracting qubits through a collective thermal bath. In  this situation was considered for a model of two noninteracting spins (qubits) interacting only with a collective thermal bosonic environment. It was demonstrated numerically in  that for some initially unentangled two-qubit states and under some conditions on the thermal bath, measurable entanglement between the two qubits is created for intermediate times. The model of  is energy conserving, ignoring relaxation processes for the qubits and including only the effects of decoherence. In  these results were extended to a more general model having (i) both local and collective thermal environments (at the same temperature) and (ii) energy-conserving and energy exchange interactions between qubits and their environments. The conditions for entanglement creation were discussed and analyzed numerically in . It was concluded that, in spite of the competition between the local thermal environments (which destroy entanglement) and the collective thermal environment (tending to create entanglement), the creation of measurable entanglement can be realized for some finite times. In both papers, [4, 11], only two qubits are analyzed hence no direct connection to many-body systems was made. As was recently shown in , the presence of a large number of indirectly interacting qubits interacting only with their common collective thermal environment could significantly modify the effective single-qubit characteristics including their relaxation and decoherence rates.
In the present paper, we consider a model of not directly interacting spins (qubits) placed in a constant effective magnetic field (oriented in the -direction). The qubits interact with both local and collective thermal environments (all at the same temperature). The collective interaction introduces an indirect qubit interaction. In the total density matrix of all qubits and environments, we trace over the variables of the environments and qubits. This gives us the time-dependent reduced density matrix for two arbitrary qubits. In the -representation it is represented by a time-dependent matrix. It is important to notice that the matrix elements, , , depend not only on the parameters of the thermal environments but also on the total number of qubits, . We study numerically the concurrence of the reduced two qubit density matrix and its dependence on the parameters of the system and on . To realize and study this situation in an experiment, one must have access to the two selected qubits (such as their particular frequencies), in order to manipulate them and prepare the initial state. Our main result is that the amplitude of concurrence, , generally decreases as increases. This means that one should not expect that the collective thermal environment can create by itself measurable entanglement even of two qubits, in the presence of many other qubits within the range of the collective environment coherence length.
1.1. Outline of Main Results
The initial state of the entire system is disentangled, a product state in which each of the spins is in a state , ; all local reservoir states are thermal and so is that of the collective reservoir, at a fixed common temperature (A generalization to a nonequilibrium situation where each local and the collective reservoir have different individual temperatures is immediate.).
1.1.1. Analytic Results
(i) Explicit Dynamics
As the spins interact with the reservoirs via energy-conserving couplings only, the reduced two-spin dynamics can be calculated explicitly; see Proposition 2.1. Consequences of the energy conservation are that populations, that is, the diagonal density matrix elements, are time independent and that the off-diagonal elements evolve independently. As an example, we discuss here the dynamics of the (1, 2) matrix element: The other matrix elements have similar behavior. Each factor on the r.h.s. has an interpretation.(i) is the initial condition of the matrix element in question. None of the other initial matrix elements are involved (energy-conserving coupling).(ii) is the uncoupled dynamics (no interaction with environments);(iii) is a dephasing factor with a time-dependent phase becoming linear for large (for the considered infrared behavior of the coupling constants in three dimensions, see (2.10)); it represents a “Lamb shift” contribution to the real part of the effective energy; this term is generated by the collective reservoir, but it is independent of the presence of the traced-out spins (the term would be the same if only two spins were coupled to the reservoirs).(iv) is a decaying factor with time-dependent decay rates, , becoming linear for large (see (2.11)). Both the local and collective reservoirs contribute; however, the term is independent of the traced-out spins (again, it would be the same if only two spins were coupled to the reservoirs);(v) is a product of oscillating terms encoding the effect of all the traced-out spins (see (2.9)). It is important to notice that only depends on the diagonal density matrix elements of the initial states of the traced-out spins. (This is so since the dynamics is energy-conserving, and tracing out any of the spins involves only the diagonal of the initial (time zero) density matrix. (See also Remark 2.2 after Proposition 2.1.)). Consequently, the two-qubit state does not depend on the initial off-diagonal density matrix elements of the traced-out “background” spins. Typically, we expect those spins to be initially in (close to) equilibrium, corresponding to vanishing off-diagonals.
Some general properties of can be explained easily for the case in which all spins are initially in the high-temperature equilibrium state . Then and its magnitude oscillate between zero and one. Due to the large power , the peaks of the function , centered at the discrete times satisfying , are of very narrow width for large . Consequently, in the limit , with held fixed, is zero for all except for , where . But the density matrix becomes very simple if , because many entries vanish (c.f. Proposition 2.1) and the corresponding concurrence is zero. It follows that, in the large limit, concurrence is zero for all times (except possibly for some isolated instances, ).
(ii) -Dependent Scaling of the Interaction
The above analysis suggests that one cannot generate two-spin entanglement for large at fixed interaction strength . However, the width of the peaked function which is of order becomes appreciable if . Hence we consider an -dependent scaling of the coupling, replacing by , for some . According to the above discussion, the borderline case is .
Starting from the explicit expressions (Proposition 2.1) and using the scaling , we calculate the limit of , for fixed. The analytic expressions we obtain for and show that the limiting dynamics does not create entanglement, for any time . While we are able to obtain explicit expressions for concurrence in the regime of , we are not so for finite. (The reduced density matrix is given explicitly for all and all , but calculating from it explicitly the concurrence is more difficult.) However, since no entanglement is generated in the limiting case, , but we know entanglement is created for (see, e.g., [4, 11]), we expect that entanglement creation decays with increasing . We study this decay numerically.
1.1.2. Numerical Results
We introduce , the highest frequency at which spin-reservoir interactions occur and call it the cut-off frequency. In the simulations, we take of the order of the thermal frequency . In the infra-red regime, our coupling is proportional to (see after (2.11)).
For , concurrence creation is maximal if both spins start out in their high-temperature state ; see Figure 1. Consequently, in the subsequent simulations, we take initial states of the two not traced-out qubits very close to this state, and we take the diagonals of the initial states of the traced-out quibts to be constant (remember that the off-diagonals of these qubits do not influence the dynamics at all). In Figure 4 we modify the initial state of the two not traced-out qubits and check that maximal concurrence is indeed obtained when both qubits are in the above state, even for large .
For general , entanglement evolves according to a rescaled time , see Figure 2. This figure shows that a reduction of diminishes the created concurrence in a moderate way. For instance, decreasing by a factor only decreases concurrence by less than .
In Figure 3 we show that the maximum of created concurrence decays with increasing . For intermediate values of (with the current parameters N~10–150) the decrease is exponential; for smaller and larger values of , it is superexponential.
In the same Figure 3(c), we study the dependence of the maximal time, , (before recurrence) at which the concurrence is not zero. We have found that this time decays exponentially in the number of spins, , for sufficiently large .
Results on the rescaled model are shown in Figure 5. We find a decrease of maximal concurrence with increasing for all . The critical value, (see analytic results above), divides the concurrence decay into two regimes. In the range, , the maximal concurrence decreases exponentially in , for intermediate values of (between 10 and 180), with a universal decay rate (i.e., not depending on ). For the decay is superexponential and varies with . We conclude that no scaling can compensate the decay of created concurrence for large .
2. Model and Reduced Density Matrix
The full Hamiltonian of the noninteracting spins coupled by energy conserving interactions to local and collective bosonic heat reservoirs is given by
Below we use dimensionless variables and parameters. To do so, we introduce a characteristic frequency, , typically of the order of spin transition frequency. The total Hamiltonian, energies of spin states, and temperature are measured in units . The frequencies of spins, , bosonic excitations, (where is the speed of light), the wave vectors of bosonic excitations are normalized by , and all constants of interactions are measured in units of . The dimensionless time is defined as .
In (2.1) and (2.2), is the frequency of spin , and denotes the of spin . is the Hamiltonian of the bosonic collective reservoir, and is that same Hamiltonian for the th individual reservoir. For a square-integrable form factor , , is given by The real numbers, and , are coupling constants, measuring the strengths of the energy-conserving collective coupling and the energy-conserving local coupling, respectively.
Since the spins interact with the reservoirs only through energy-conserving channels, this model is exactly solvable. For simplicity of exposition, we take We also take, for simplicity, all local form factors equal () and all collective ones also ().
Fix any pair of spins, and (re-)label their frequencies by and ; see (2.1). We write the reduced density matrix, , of the two fixed spins as a matrix in the ordered energy basis (Another equivalent notation is: .) where . For instance, .
The initial state of the spins is a product state of the form , where with and . The upper bound on the off-diagonal guarantees that the eigenvalues of are nonnegative.
We introduce the quantities:
The integrals in (2.10) and (2.11) are made to converge introducing a suitable cut-off wavenumber, , or cut-off frequency, . (Here we use dimensionless units.) For instance, for numerical simulations, we choose as form factor the function , where if and otherwise.
We also define to be the same as , but with replaced by . With this notation, we have the following result.
Proposition 2.1 (explicit dynamics of the reduced density matrix). The evolution of the density matrix is given by and the populations are constant, , for and .
Remark 2.2. (1) The effect of spins is contained entirely in the factors and . They only depend on the initial populations , (see (2.9)), but not on the off-diagonals, . This is explained by the fact that when tracing over a single spin, , we perform the operation , where are operators commuting with (energy-conserving interactions only!). Clearly the latter trace only involves the diagonal of .
(2) The oscillatory phases, , in (2.12) represent the free, uncoupled dynamics of the spins. Consider the modified two-spin density matrix: (“interaction picture” dynamics of ). Because and are related by conjugation of a unitary operator of the product form, , the concurrences of and are the same. In other words, when examining concurrence of , we may use formulas (2.12) with .
Recall that the concurrence of the reduced density matrix, , is unchanged when we pass to the interaction picture (see the remark explaining (2.13)). In the basis (2.7), the evolution of (2.13) is given by (2.12) with and where the initial condition is ():
2.1.1. Variation of Dependence
For homogeneous initial conditions, for , we have Unless , oscillates in between its minimum value (when ) and its maximum value (when ) (We have ) The width of the oscillations becomes very narrow with increasing . In the limit of large , is zero for all times, except for the discrete set of satisfying , in which case .
This implies that, for large , all off-diagonal density matrix elements of vanish with the exception of (and , of course) for almost all values of . This suppression of off-diagonals comes from the large number of particles and is mediated through the collective energy-conserving interaction. (For we have .)
In order to try to have a nontrivial dynamics for large , one may scale the collective conserving coupling constant as Then (2.15) becomes An expansion in large yields Thus as ,
Remark 2.3. (1) By replacing, in these limits, by , we obtain the corresponding limits for .
(2) For , vanishes as .
(3) The rapid oscillating phases disappear if (any ) or (any ).
2.1.2. Asymptotic Concurrence ()
The reduced two-spin density matrix (in the interaction picture) at time is from which we obtain the concurrence Remember that . This shows that the asymptotic dynamics cannot create entanglement at any time.
Call the r.h.s. of (2.19) (a quantity still depending on unless ). By replacing by in (2.19), we obtain the limit of , which we call . For we have the relation . The reduced two-spin density matrix (in the interaction picture) at time is, for , where The density matrix (2.22) is of the product form: This shows that the asymptotic dynamics is factorizable and cannot create entanglement at any time.
3. Numerical Results
Let us first consider the case of two spins only, . In (2.12) we put for simplicity and regard and as two independent parameters. Taking both spins initially in the same state given by , , see (2.8), we examine the maximal concurrence, as a function of and , for arbitrary fixed values of and .
We find that, for fixed , , the maximal concurrence is given at . Having such values fixed and plotting the concurrence as a function of , the maximal concurrence is realized when ; see Figure 1, where a plot of the concurrence as a function of is shown. Maximal generation of concurrence is thus obtained starting from pure-state initial conditions for each spin.
Let us now consider the case of spins. For concreteness we choose for all spins, (the traced-out ones and the two not traced-out ones). For the two not traced-out spins we take off-diagonals . (Then is close to which favors larger entanglement creation.) Recall that the dynamics is independent of the off-diagonals of the traced-out spins (i.e., we do not have to specify the of the traced-out spins).
As mentioned after (2.11), we choose the form factor, , with the cut-off frequency equal to the thermal frequency, , at room temperature, K.
In Figure 2 we investigate the effects of an increase in the coupling parameter, . The first effect is a time shift for the concurrence evolution, described by a scaling, ; see Figure 2(a), where is the cut-off frequency. The second effect is a reduction of the maximal concurrence in a smooth way; see Figure 2(b). As one can see, the effective decrease in amplitude for not too strong coupling strength, , is relatively small. For instance, changing for by one order of magnitude from 0.04 to 0.4 changes the amplitude by only 27%. The percentage change is almost the same for larger values; see Figure 2(b).
In Figure 3 we show that the creation of concurrence decreases with the number of spins. In (a) we plot the concurrence as a function of (rescaled) time for various values of . As one can see, the same time rescaling is also valid for . Moreover, the maximum concurrence created, , reported in (b), decreases exponentially in in the range and faster than exponentially outside this range. For larger , the concurrence decays superexponentially in . For exceeding , the concurrence becomes too small to be significant (of the order ).
It is also interesting to note that the graph of concurrence shows collapses and revivals and that the revival times for are always less than the revival time for . It is also interesting to consider how the collapse time, , defined as the first time at which concurrence drops abruptly to zero, depends on the number of spins, . This study has been reported in Figure 3(c) and shows that the rescaled collapse time decays exponentially with the number of spins: , where . (See dashed line in Figure 3(c).)
One can also vary the initial conditions for the spins by choosing independent and , while all other spins have the same value , (their off-diagonals do not influence the dynamics at all).
In order to simplify the problem, we also set and and consider the maximal concurrence as a function of two independent parameters, and , only.
An example of the 3D plot obtained is reported in Figure 4: the maximal concurrence is realized for , independently of (in the picture , but similar plots are obtained for other values of ).
The numerical analysis of the rescaled model with replaced by shows that the concurrence is always a decreasing function of and that the maximum of the created concurrence is a universal function of the number of spins , independent of for .
Results are shown in Figure 5, where the dashed line is the best exponential fit , with , the best fitting value, for the cases . The same figure shows that when , the decay is superexponential and no universality occurs.
This suggests that no power law scaling with of the coupling strength can compensate the rapid decay of concurrence with the number of spins.
Qualitatively similar results, not reported here, can be obtained by changing the ratio between the thermal and cut-off frequency in the range .
We have analyzed the two-qubit (two effective spin) entanglement in an -qubit open system, where individual qubits interact through a collective thermal environment. We have demonstrated that concurrence (a measure of two-qubit entanglement) quickly decays with increasing number of surrounding qubits. It follows from our consideration that, for creation of entanglement by implementing our approach, one has to use a small number of qubits collectively interacting with the thermal environment, preferably only two qubits. In this paper, we consider a quantum noise generated by a thermal bosonic environment, interacting in a purely dephasing way with qubits, but no Markovian approximation has been made. In contrast, it has been shown [16–18] that an -qubit register can be driven to entangled target states (within a wide class of states), under an appropriately engineered Markovian dissipative dynamics in Lindblad form and that the convergence rate to the target state is independent of . This entirely different noise has thus a completely different effect on entanglement.
The results of this paper are important for a better understanding and a characterization of the collective thermal environment and its ability to create and destroy entanglement in many-particle open quantum systems. Any quantum information processor could be used for the study of entanglement induced by a thermal environment, as, for example, two (or more) ions in an ion trap quantum computer  or two (or more) superconducting qubits  or an NMR quantum information processor  operating at room temperature. In the latter case, the experiment is performed with an ensemble of many identical molecules, each having a similar number of weakly interacting nuclear spins. One can study in this system, for example, how the collective vibrational modes of a molecule (thermal environment) influence the two-spin “effective entanglement” in the ensemble of molecules, depending on the number of spins in a molecule. Important and relevant systems are represented by photosynthetic complexes in which quantum effects for exciton dynamics and primary charge separation have been recently discovered [20–23]. Here, the effective spins are exciton or electron sites (pigments), and the thermal environment represents the vibrational protein modes.
This work Supported by NSERC under Discovery Grant 205247 (M. Merkli). This work was carried out under the auspices of the NNSA of the U. S. DOE at LANL under Contract No. DEAC52-06NA25396. Supported by the LDRD program at LANL (G. P. Berman).
- M. Cramer, M. B. Plenio, and H. Wunderlich, “Measuring entanglement in condensed matter systems,” Physical Review Letters, vol. 106, no. 2, Article ID 020401, 2011.
- M. A. Nielsen and I. L. Chuang, Quantum Computation and Quantum Information, Cambridge University Press, Cambridge, UK, 2000.
- M. Sarovar, A. Ishizaki, G. R. Fleming, and K. B. Whaley, “Quantum entanglement in photosynthetic light-harvesting complexes,” Nature Physics, vol. 6, no. 6, pp. 462–467, 2010.
- D. Braun, “Creation of entanglement by interaction with a common heat bath,” Physical Review Letters, vol. 89, no. 27, pp. 2779011–2779014, 2002.
- T. Yu and J. H. Eberly, “Finite-time disentanglement via spontaneous emission,” Physical Review Letters, vol. 93, no. 14, Article ID 140404, pp. 1–4, 2004.
- T. Yu and J. H. Eberly, “Sudden death of entanglement,” Science, vol. 323, no. 5914, pp. 598–601, 2009.
- T. Yu and J. H. Eberly, “Sudden death of entanglement: classical noise effects,” Optics Communications, vol. 264, no. 2, pp. 393–397, 2006.
- B. Bellomo, R. Lo Franco, and G. Compagno, “Non-Markovian effects on the dynamics of entanglement,” Physical Review Letters, vol. 99, no. 16, Article ID 160502, 2007.
- J. H. Huang and S. Y. Zhu, “Sudden death time of two-qubit entanglement in a noisy environment,” Optics Communications, vol. 281, no. 8, pp. 2156–2159, 2008.
- J. P. Paz and A. J. Roncaglia, “Dynamics of the entanglement between two oscillators in the same environment,” Physical Review Letters, vol. 100, no. 22, Article ID 220401, 2008.
- M. Merkli, G. P. Berman, F. Borgonovi, and K. Gebresellasie, “Evolution of entanglement of two qubits interacting through local and collective environments,” Quantum Information & Computation, vol. 11, no. 5-6, pp. 390–419, 2011.
- M. Merkli, G. P. Berman, and A. Redondo, “Application of resonance perturbation theory to dynamics of magnetization in spin systems interacting with local and collective bosonic reservoirs,” Journal of Physics A, vol. 44, no. 30, Article ID 305306, 2011.
- M. Merkli, I. M. Sigal, and G. P. Berman, “Resonance theory of decoherence and thermalization,” Annals of Physics, vol. 323, no. 2, pp. 373–412, 2008.
- M. Merkli, G. P. Berman, and I. M. Sigal, “Dynamics of collective decoherence and thermalization,” Annals of Physics, vol. 323, no. 12, pp. 3091–3112, 2008.
- M. Merkli, I. M. Sigal, and G. P. Berman, “Decoherence and thermalization,” Physical Review Letters, vol. 98, no. 13, Article ID 130401, 2007.
- F. Verstraete, M. M. Wolf, and J. Ignacio Cirac, “Quantum computation and quantum-state engineering driven by dissipation,” Nature Physics, vol. 5, no. 9, pp. 633–636, 2009.
- B. Kraus, H. P. Büchler, S. Diehl, A. Kantian, A. Micheli, and P. Zoller, “Preparation of entangled states by quantum Markov processes,” Physical Review A, vol. 78, no. 4, Article ID 042307, 2008.
- S. Diehl, A. Micheli, A. Kantian, B. Kraus, H. P. Büchler, and P. Zoller, “Quantum states and phases in driven open quantum systems with cold atoms,” Nature Physics, vol. 4, no. 11, pp. 878–883, 2008.
- J. M. Martinis, “Superconducting phase qubits,” Quantum Information Processing, vol. 8, no. 2-3, pp. 81–103, 2009.
- H. Lee, Y. C. Cheng, and G. R. Fleming, “Coherence dynamics in photosynthesis: protein protection of excitonic coherence,” Science, vol. 316, no. 5830, pp. 1462–1465, 2007.
- G. Panitchayangkoon, D. Hayes, K. A. Fransted et al., “Long-lived quantum coherence in photosynthetic complexes at physiological temperature,” Proceedings of the National Academy of Sciences of the United States of America, vol. 107, no. 29, pp. 12766–12770, 2010.
- E. Collini, C. Y. Wong, K. E. Wilk, P. M. G. Curmi, P. Brumer, and G. D. Scholes, “Coherently wired light-harvesting in photosynthetic marine algae at ambient temperature,” Nature, vol. 463, no. 7281, pp. 644–647, 2010.
- G. L. Celardo, F. Borgonovi, M. Merkli, V. I. Tsifrinovich, and G. P. Berman, “Superradiance transition in photosynthetic light-harvesting complexes,” Submitted, http://arxiv.org/abs/1111.5443.