ISRN Mathematical Physics

Volume 2014 (2014), Article ID 264956, 19 pages

http://dx.doi.org/10.1155/2014/264956

## Entanglement Dynamics of Second Quantized Quantum Fields

NanoScience Technology Center, Department of Physics, University of Central Florida, Orlando, FL 32826, USA

Received 26 November 2013; Accepted 16 December 2013; Published 28 January 2014

Academic Editors: B. Bagchi, G. Berginc, and M. Sebawe Abdalla

Copyright © 2014 Mikhail Erementchouk and Michael N. Leuenberger. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

We study the entanglement dynamics in the system of coupled boson fields. We demonstrate that there are different natural notions of locality in this context leading to inequivalent notions of entanglement. We concentrate on the particle picture, when entanglement of one particle is determined by one-particle density matrix. We study, in detail, the effect of interaction preserving populations of individual one-particle states. We show that if the system is initially in a disentangled state with the definite total number of particles and the dimension of the one-particle Hilbert space is more than two, then only potentials of the special form admit complete entanglement, which is shown to be reached at NOON states. If the system is initially in Glauber’s coherent state, complete entanglement is not reached despite the presence of two entangling channels in this case. We conclude with studying the time evolution of entanglement of photons in a cavity with multiple quantum dots in the limit of large number of photons. We show that in a relatively short time scale the completely entangled states belong to the class of graph states and are formed due to the interaction with dots in resonance with the cavity modes.

#### 1. Introduction

Entanglement [1] is a quintessentially quantum feature. It signifies that different parts of a compound system may form a new entity, a complex; for example, when neither of two particles can be characterized by a definite state so that instead of two particles one has to consider a pair and so on. If one would attempt to access a particular part of the complex by performing a local measurement, this would unavoidably modify the state of other parts even if the direct interaction between the parts is absent or negligible. Such departure from the classical properties makes entanglement the central object in various contexts, from the perspective of application in quantum informatics [2] to understanding the physics of quantum phase transitions [3]. As a problem of special interest, therefore, the problem of preparation of a system in an entangled state stands out. Among different appearances of this problem, entangled states of quantized electromagnetic field, as perhaps the most accessible and the most flexible object, presents a significant importance on its own. Nowadays, the most developed and widely used method of generating entangled photons is the parametric down conversion [4–6], which is based on the two-photon radiative decay of material excited states. This method, however, suffers from well recognized intrinsic limitations such as very low yield and rescaling the wavelength of the emitted photons [2, 7, 8]. Therefore, there is the constant search for alternative sources of entangled light [9–19], which motivates a thorough consideration of entanglement of quantum many-body states.

The characteristic feature of the process of entangling photons in the course of interaction with matter is the nonconserving number of photons, or involved particles in general. In the few particle limit or in the case when the typical time scales are well separated, these processes can be reduced to a more or less standard quantum mechanical situation. This makes the process of entangling fit into the well developed description of entanglement. Indeed, in this case one can specify time periods when the system is either in the state of excited matter and no (relevant) photon, or in the state when matter is in the ground state and there are emitted photons. The general problem of solid based sources of entangled light, however, requires addressing the more general situation, when the photon states coexist with material excitations and the processes of reabsorption and reemission may play an important role. In this case one has to incorporate the nonconserving number of particles fully into consideration. The quantum field description (more precisely, the formalism of second quantization) provides the most natural framework for dealing with this kind of a situation. In this approach particles appear not as predefined entities as, say, qubits within the standard quantum mechanical treatment, but rather as excitations of the respective quantized fields.

As we will demonstrate below, entanglement in the context of quantized fields has its subtleties. The complete description of entanglement in this context is yet to be developed. For our purposes, however, a basic approach is sufficient. First, we will restrict our attention only to the case when the whole system can be characterized by a pure state. Second, entanglement will be considered from the perspective of the problem of sources of entangled light, which advances consideration of particle properties, as will be elaborated below. From this point of view the existence of entanglement implies that a set of particles is in a nonseparable state and, as a result, one particle cannot be described by a state vector; the particle of necessity is in a mixed state. Pure and mixed states, in turn, can be distinguished using the fact that extremal values of observables corresponding to operators with a nondegenerate spectrum are reached at pure states. Let be an operator acting on a finite dimensional Hilbert space and let its spectral decomposition be , where are the eigenvalues, enumerated in the ascending order, , and are the projectors on the respective eigenspaces. Furthermore, let the smallest eigenvalue be nondegenerate, that is, , then the minimal value of (understood as a function of state described by the density matrix ) is equal to and is reached at corresponding to the pure state, . Respectively, the value of the observable on mixed states will always satisfy inequality . Indeed, on the basis, where is diagonal, one has with the equality reached only when all but are zero. Thus, in a sense any operator with nondegenerate extremal eigenvalues allows distinguishing between pure and mixed states. Using this observation one can implement the ideology of witnessing entanglement, which was originally developed in the context of the full compound system [1, 20] and later was extended to one particle description [21].

This approach can be directly applied to the quantum field description. In this case, the role of witnesses is played by respective one-particle operators and following the standard line of arguments (see, e.g., [22]) one can characterize entanglement of one particle with the rest of the system using one particle correlation matrix (OPCM) , where and enumerate one-particle states and and create and destroy a particle in the respective states. Throughout the paper we incorporate the time dependence into the Heisenber, representation of the operators , where is the Hamiltonian describing the whole system. Thus, the time dependence of the OPCM is given by where the average is taken with respect to the initial state .

In order to illustrate the difficulty of producing entangled states let us consider a simple example of a boson field driven by an external source where are -numbers characterizing the external source.

Before applying rigorous methods let us note that this system may look confusing if approached with the help of often employed arguments based on interference of different paths connecting the initial and final states. Indeed, considering that there are different ways to fill some particular state, say, with only two particles, one might expect that these two particles will become entangled, which, of course, is an incorrect conclusion.

In virtue of the discussion above, the absence of entanglement would be manifested by rank one of the OPCM [23, 24], while the rank of the OPCM can be easily investigated. The solutions of the operator equations of motion have the form where are -number functions. Substituting this representation into (1) we find where . Generally the structure of the OPCM (throughout the paper hats denote matrices in one-particle Hilbert space) driven by the external source depends on time nontrivially. The last three terms in (4) may lead to variation of entanglement depending on the structure of the initial state. If, however, the system is initially in the vacuum state, that is, , then , where is the matrix with the elements . In turn, is at most of rank which implies the absence of entanglement.

This example shows that the problem of the dynamics of entanglement should be treated with certain care. First, the naive arguments based on the picture of interference of different paths may be misleading. Second, no matter how complex the internal dynamics of the system described by the spectrum, , and independently on particular time dependence of the external excitation , the states, reached out of vacuum under the action of this excitation, are disentangled. In particular, the initial state, which may lead to a nontrivial time dependence of entanglement, as has been mentioned above, must be created by means other than the external classical excitation.

The rest of the paper is organized as follows. In Section 2, we provide a more detailed description of entanglement in the context of quantized fields. In Section 3 we consider the entanglement dynamics for the system of linearly coupled boson fields. In Section 4 we study entanglement dynamics for the case of self-interacting boson field. Finally, in Section 5 we apply the obtained results to the analysis of the entanglement dynamics of photons in a cavity with multiple quantum dots in the limit of a large number of photons.

#### 2. Entanglement within Field and Particle Pictures

Some results presented in Sections 3 and 4 may seem to contradict some results readily available in the literature. This reflects a certain ambiguity of the notion of entanglement in the context of quantized fields. Therefore, in order to avoid a possible misunderstanding, it is worthwhile to analyze the problem of entanglement in detail.

##### 2.1. Field and Particle Pictures

The ambiguity stems from the fact that entanglement is understood as a relation between a part of the compound system and the whole system, while there are two distinctive notions of the part when quantized fields are considered: fields and particles. Which of these two different entities, fields and particles, appear more naturally is dictated by the physical content of the specific problem. For example, if a system of harmonic or anharmonic oscillators is considered [25], individual oscillators (i.e., fields) stand out as part of the big system, while, say, the dynamics of excitons in semiconductors [10, 12, 14, 16, 17, 19] promotes consideration of particles.

An arbitrary (pure) state of a system of (boson) fields can be specified in terms of degrees of excitation (population numbers) of each field where denotes the population of the -th field and are the respective amplitudes. Equivalently, this state can be presented in terms of the particle creation operators where is vacuum, is the total number of particles, and creates a particle in the -th one-particle state. Representations (5) and (6) are equivalent if the amplitudes satisfy conditions that follow from the commutation of the creation operators and the relation between the individual fields in (5) and one-particle states in (6), .

While representations (5) and (6) are equivalent they imply different notions of “locality.” This follows simply from the following observation. An operator acting on a particular *field* in (5) does not preserve the total number of particles unless it is proportional to the identity operator. At the same time an operator, which changes the state of one *particle* in (6), obviously affects the population numbers of more than one field. In turn, entanglement strongly relies on the notion of locality and, therefore, is sensitive to the choice of eligible local transformations. Therefore, one should distinguish representations (5) and (6). For this reason we will call representation (5) *the field picture* and (6) will be referred to as *the particle picture*. Most of the papers dealing with entanglement in the field context are effectively restricted to one of the pictures implied by the physical situation. This, however, leads to possible ambiguities because results are often formulated in some general terms (entanglement, part, compound system, and so on), which are identical for both pictures. Entanglement, however, does depend on the picture, as is suggested by the nonequivalence of the notions of locality and will be demonstrated below.

Within the field picture, if the general form of operators acting on a particular field is allowed, entanglement is naturally related to separability of amplitudes . In systems with superselection rules [26], the space of allowed operators is “smaller” and the operators, whose mean value is extremal at the particular state, may not be accessible. From the perspective of witness ideology, this means that such state should be considered as entangled. For such systems, therefore, one needs a criteria of entanglement different from the mere separability of the respective amplitudes. This problem, however, goes beyond the scope of the present paper.

In order to quantify entanglement of the -th field with the rest of the system, that is, with the other fields, it is convenient to introduce the reduced density matrix where the summation excludes . If the rank of the reduced density matrix is higher than one, we have entanglement, which can be quantified, for example, by the von Neumann entropy. Thus, entanglement in the field picture fits the canonical quantum mechanical description based on separability of the amplitudes and agrees with the general witness ideology.

Within the particle picture, however, one has to rely upon the witness approach more due to indistinguishability of individual particles, which makes the separability property of amplitudes an inadequate criterium. The quantity of main interest becomes the one-particle correlation matrix which is related to the amplitudes in a rather complex way. The one-particle density matrix (OPDM) is defined as a properly normalized correlation matrix . Then, entanglement can be quantified, for example, by von Neumann entropy (we will return to the question of the entanglement measure below) . One can easily check that when the system immediately admits the standard description (e.g. when all particles are distinguishable, i.e., all in (6) are different, and the field and particle pictures are identical) the approach based on the OPDM yields results consistent with this description. Of course, OPDM may only answer questions regarding entanglement of a single particle with the rest of the system. If one is interested in more subtle details, such as, for instance, entanglement in pairs, one has to look at the density matrices of higher orders. We, however, limit ourselves to studying the basic properties of entanglement and for this purpose it suffices to consider OPDM, which can be shown to yield the upper bound for entanglement in the system.

The necessity to distinguish entanglement within different pictures is illustrated by states with amplitudes (in the field picture) being separable, , and for each field only amplitudes of the same parity are not zero, for example, for all and integral . These states are disentangled within the field picture; they, however, may be completely entangled within the particle picture. Indeed, the off-diagonal elements of the OPDM vanish for such states and if, additionally, the amplitudes are chosen in such manner that the average number of particles in each state (the degree of excitation of each field) is the same, , then the von Neumann entropy takes the maximal value, implying maximal entanglement within the particle picture.

The example of opposite situations is presented by single-particle entanglement [21, 27–30], when, say, , while all other amplitudes in (5) are zero. These states are obviously entangled within the field picture while disentangled in the framework of the particle picture.

##### 2.2. Disentangled States in Particle Picture

For a more detailed comparison of different pictures we show that the states disentangled within both field and particle pictures constitute the special class and this class is not particularly rich. Besides the states with only excitations of one field, when trivially there is nothing to entangle (we would like to remind the reader that we consider the case without superselection rules), are canonical coherent states [31–33].

First of all we describe all disentangled states within the particle picture. The main result here is almost obvious: *particles are disentangled within the particle picture if and only if they all are in the same state*. In order to prove this statement (nontrivial in the part that there are no other disentangled states) it is convenient to introduce special notations.

We introduce a vector-operator with components . The commutation relation in terms of the vector-operators can be formally written as . These are vectors in the following sense. The choice of a different set of basis one-particle states in space corresponds to a linear transformation of , which translates into choosing a different set of operators linearly related to the old set or . The commutation relation for new operators can be shown to have the form . Thus, in order to satisfy the boson commutation relation, the vector-operators and must be related through a unitary transformation. These transformations, in particular, preserve the operator of the total number of particles. This can be illustrated by presenting the operator as . In this section denotes the product , which transforms as a tensor.

The reason why we have introduced vector-operators is that OPCM also transforms as a tensor that is, is mapped into . Thus, unitary transformations do not change the spectrum of OPCM and, due to invariance of , they leave entanglement intact.

The important consequence of this geometrical picture is that it immediately provides the description of all disentangled states. OPCM is a Hermitian matrix and, therefore, can be diagonalized by a unitary transformation, implying the existence of the preferred set of operators. In turn, OPCM of *any* disentangled state in the diagonal form has only single nonzero element. Thus, choosing the appropriate transformation we can have only while all other elements of OPCM are zero. The only states yielding such OPCM are of the form
where . Using relation (9) we can expand . As follows from unitarity, the single column of a unitary matrix is a unit vector. Denoting this vector by and defining the “scalar" product , we can parametrize all disentangled states
by a set of amplitudes and a vector on the unit sphere in . There is a one-to-one correspondence between such sphere and the space of one-particle states; therefore, we will sometimes refer to vector as a state. Disentangled states with a definite number of particles (with for all but specific ) will play an important role in the following consideration. Extending the terminology used for the case we will call them *spin coherent states*.

Once we have established the general form of disentangled states in the particle picture we may proceed and find which of those are disentangled in the field picture. In order to find amplitudes we first notice that only the term with the same total number of particles in (12) contributes to and from the polynomial expansion of this term we need only one term with matching population numbers for each field. Thus, we find
with and the important convention that even if . It follows from (13) that there are only two possibilities to have separability of . The first one is when all but one are zero. This corresponds to the trivial case when only the single type of fields is excited. The second possibility is when with some complex numbers and . Substituting these amplitudes into (12) and enforcing the normalization condition we find that all states disentangled in both field and particle pictures can be presented as
These are *Glauber’s coherent states*. It is interesting to emphasize a relation with the example considered in the introduction. Hamiltonian (2) is diagonalized by introducing operators , which are obtained by employing Glauber’s shift operator with . This, in particular, proves that states reached from vacuum under classical excitation are disentangled within both pictures.

##### 2.3. NOON States

The general “inverse” problem of a relation between entanglement and the structure of respective states requires consideration, which goes far beyond the objectives of the present paper. Therefore we limit ourselves to an explicit description of a few sets of completely entangled states.

In Section 4 states of NOON type appear. While these states were introduced for qubits [34], they can be defined in a more general setup as follows. Let vectors , form a basis in . Then NOON states of particles with dimensional one-particle Hilbert space are defined as Thus, taking into account that different rearrangements of produce the same state, the manifold of NOON states is isomorphic to , where is the symmetric group (group of all permutations of elements).

It should be noted that NOON states do not exhaust all completely entangled states. In order to see this it is useful to present NOON states as they appear in a more general context. The natural representation of on the basis has the form for any . Taking any , NOON states can be defined in terms of the orbits where we have taken into account that the orbit on visits each element times as runs over . Representation (16) demonstrates that NOON states are a particular class of entangled coherent states [35, 36]. Namely, NOON states are equally weighted superpositions of orthogonal spin coherent states.

Using the same approach as for (16) another set of completely entangled states can be constructed for . For any pair such that states are completely entangled. This makes an interesting connection between completely entangled states (constituting a straightforward generalization of Dicke states [37, 38]) and irreducible representations of symmetric group .

##### 2.4. One-Particle Entanglement from Lie-Algebraic Perspective

The natural framework for dealing with entanglement in particle picture is provided by the language of irreducible representations of Lie algebras. For example, in the case of a two-dimensional one-particle Hilbert space (e.g., photons characterized by two “+” and “−” polarizations), it is convenient to employ Schwinger’s model of angular momentum [39] and to introduce which satisfy the commutation relation of algebra with being a completely antisymmetric tensor. These operators provide a representation of OPDM. Consider where runs over and are Pauli matrices and is the identity matrix. The sum over describes a deviation of from the identity matrix and hence a deviation of the state from a completely entangled one. This implies that entanglement can be expressed in terms of .

Before formulating a general approach let us consider a two-particle case, which provides a clear connection with the standard quantum-mechanical consideration. On the basis of the population numbers, any two-particle state can be presented as where denotes the state with particles in the “+”-mode and particles in the “−”-mode. Alternatively, the state can be presented as [23] where is a symmetric matrix. Comparing (20) and (21) one finds , and . From (21) one obtains .

The OPCM is a matrix and, therefore, its eigenvalues are completely determined by and with being the total number of particles. Entanglement, in turn, is determined by the normalized eigenvalues , which are found as where is the concurrence [40]. In order to see the relation with the standard definition of the concurrence in the two-particle case we introduce the “spin flip" transformation , and the spin flip state , then . Thus, the two-particle case completely fits into the canonical quantum-mechanical description and, moreover, concurrence is tightly connected to average angular momentum; for an arbitrary one has .

Now we can turn to a general case. Let with be generators of , which are presented in by traceless Hermitian matrices normalized according to [41]. It can be checked that operators constitute a symmetric representation of on the boson Fock space and this representation is irreducible within each sector with a definite total number of particles.

Any Hermitian matrix can be uniquely presented as a linear combination of and the identity matrix . Taking into account the chosen normalization of we find a representation of OPDM generalizing (19) [42] where and is a Bloch vector with components .

As follows from (24) vector completely determines the spectrum of OPDM and, therefore, any measure of entanglement in the particle picture is a function of or, more precisely, of its invariants with respect to adjoint representation of [42]. It is worth noting that since operators correspond to a set of single-particle observables, this implies that entanglement is an observable.

Entanglement is a global characteristic of the spectrum of OPDM. Qualitatively, it is related to the density of the distribution of eigenvalues of ; if they, including zeros, are far away from each other, entanglement is weak; if they are close, entanglement is strong. Taking into account that the simplest characteristic of the spectrum is the averaged squared distance between the eigenvalues It changes from , when all (complete entanglement), to , when one eigenvalue is and the remaining are (absence of entanglement). By expanding (25) we find This establishes the relation between and the linear entropy

On the other hand, from (24) we find . Thus, changes from for completely entangled states to for disentangled ones. The condition has been considered as a criterion for complete entanglement from different perspectives: maximization of quantum fluctuations [43] and geometric invariant theory [44]. We would like, therefore, to emphasize that within the present approach this condition is a consequence of the physical requirement that an entangled particle cannot be characterized by any definite one-particle state, which allows for a straightforward generalization for the case of indefinite number of particles.

Other measures of entanglement may involve more subtle characteristics of the distribution of the spectrum of OPDM and do not have to coincide with the linear entropy [45]. It should be noted, however, that the choice of essentially different measures is limited. Clearly, any measure should be invariant with respect to basis transformations and, therefore, should be a function of invariants of OPDM only. For matrices there are invariants determined by the exact positions of eigenvalues of OPDM constrained by the requirement . Thus, any measure can be expressed in terms of with . In particular, for there is only one invariant and, therefore, any measure of entanglement is a function of a single parameter , which in this case has the meaning of the magnitude of average angular momentum. For example, in the case, entanglement without any loss of generality can be quantified by the linear entropy Moreover, if one is interested mostly in limiting situations of complete entanglement and disentanglement, a description employing only is sufficient also for the case of general . Therefore, in the following consideration, we will use the linear entropy for characterization of entanglement. It is especially convenient because , which effectively allows one to calculate without dealing with the problem of specific representations of [46, 47].

#### 3. Entanglement Transfer between Coupled Fields

We begin our analysis of dynamics of entanglement by considering a simple but important case of linearly coupled fields described by the Hamiltonian where and enumerate the modes of the fields; , , and are the spectra of the fields and the coupling constants between them, respectively, and the operators and are assumed to obey the boson commutation relations.

The important feature of the time evolution of entanglement in this system is that the total entanglement remains constant and is solely determined by the initial state. For a more precise formulation, instead of the operators and , let us introduce the combined operators . That is instead of two fields with two sets of modes and we consider the single field, whose modes are the direct sum . Thus, if and if . For field we define the OPCM as When both and belong to, say, the respective matrix elements give the OPCM for field and so on. Furthermore, let have the spectral representation where and and are the eigenvalues and the unit eigenvectors and denotes the tensor product, which is defined as . Then the total entanglement can be characterized by , where .

The important feature of systems with linear coupling between fields is that the entanglement evolution is restricted purely to its redistribution between the fields, while the total entanglement remains constant, . This result while specific for the particle picture (thus contrasting the results obtained, e.g., in [25, 48]) is in a more general context agreeable with the description of entanglement in terms of irreducible representations of Lie groups. Therefore, we prove it for a system described by the Hamiltonian , with Hermitian and *pairs* of field operators obeying the commutation relation . The total entanglement is determined by the spectrum of , which satisfies . The spectrum of matrices, whose time dependence is governed by such equations with the commutator in the r.h.s., does not change with time. Indeed, the solution has the form . Thus, due to unitarity of , the spectral representation of is given by (31), with constant , and only the eigenvectors are functions of time. This implies that entanglement is an integral of motion.

For example, if initially there was only a single nonzero eigenvalue (i.e., entanglement was zero) it remains the only one later on, implying no production of entanglement. The same result holds for the OPCM corresponding to either fields or . Indeed, the OPCM of field is obtained from the OPCM of the combined field applying the respective projection operators, , so that Using for its spectral representation one can see that such projection cannot increase the rank of the OPCM.

While inability of linear coupling to entangle initially disentangled fields can be expected [49], it should be noted that the evolution of entanglement of initially entangled states is nontrivial. As an example, let us consider the situation of small total entanglement, more specifically, when there are only two nonzero terms in (31) corresponding to with . According to (32) the OPCM for field is given by where and . Thus, the value of entanglement depends on the magnitude of the projections , which is determined by the internal dynamics of the coupled fields. In particular, if one may have , resulting in the significant entanglement of particles , either with each other or with particles . For more concrete information one has to take into account that, generally speaking, (33) may not be the spectral representation of the matrix because the vectors and are not necessarily orthogonal. Let . Then the eigenvalues of are found to be Thus, if at some particular instant one has , then, depending on the details of the dynamics of the eigenvectors, one may have . Strong entanglement is produced when the “weak” component of the OPCM is transferred more effectively than the major component , whose only small part is moved into field during the evolution. At the same time, as can be seen from (34) after a slight change of notations, such spike of entanglement between particles is accompanied with disentanglement of particles .

The main condition for this geometric effect is the smallness of the projection of the respective eigenvector of the total OPCM. As a result, the characteristic feature of the OPCM of strongly entangled states in this case is . That is, the states with developed from the states with low total entanglement, are characterized by low excitation, while the strongly excited field, say the field in the considered example, for which one has , remains only weakly entangled.

The fact that entanglement of the specific particles is determined by the projections of the total OPCM may lead not only to increased entanglement of the particles but also to disentanglement. Indeed, if for the OPCM given by (31) with at some instant the vectors and belong to different subspaces (say, and ) then entanglement of particles of both types, and , are zero. One can understand this effect introducing the isospin quantum number, so that one value of isospin corresponds to the particle and another one stands for the particle . Thus, non-zero total entanglement coexisting with zero entanglement of specific particles can be interpreted as stored in isospin.

These effects are of a general nature and do not depend on the number of particles and may play an important role in the dynamics of entanglement in specific systems, for example, the time dependence of entanglement of initially entangled light in a leaky cavity.

#### 4. Entanglement Produced by Interaction

As has been demonstrated above, simple linear coupling of quantum fields is not sufficient for the fields to become entangled. We show in this section that the situation is different when there is an interaction in the system. In order to illustrate some important features of the dynamics of entanglement, we consider evolution the of initially disentangled states in systems whose dynamics preserves populations of one-particle states, so that the Hamiltonian can be written as , where is the operator of the number of particles in -th one-particle state, are energies of these states, and are the matrix elements of the interaction potential.

Such systems demonstrate nontrivial entanglement dynamics while admitting an exact solution. First, we consider a simple example of a system with 2d one-particle Hilbert space, which we will refer to as an model, then we will generalize this consideration to an arbitrary dimension, model.

##### 4.1. -Model

In the normal form the Hamiltonian preserving population of individual one-particle states has the form Here the first two terms describe the internal dynamics of the modes and the last term represents the interaction between the modes. Following the discussion in Section 2 it is convenient to rewrite in terms of generators of defined in (18). Consider where and depend on the total number of particles , and The dynamics governed by Hamiltonian (36) was a subject of numerous studies [50–52]. It is constructive, however, to consider its properties from the perspective of the formalism described in Section 2 and establish a connection with entanglement in the particle picture. In particular, we show that the time evolution of entanglement essentially depends on the structure of the initial disentangled state; if initially the system is in spin coherent state it may evolve into a completely entangled state; however, if the system is initially in Glauber’s coherent state (i.e., disentangled in both field and particle pictures) it cannot be completely entangled. Moreover, the structure of maximally entangled states is qualitatively different for these initial states.

From it follows that the dynamics of entanglement is subject to the general constraint . Thus, the variation of entanglement is determined by the change with time of the “transversal” component , which is unaffected by rotations around the -axis. This imposes an upper limit on the value of entanglement produced by the interaction, . The symmetry with respect to rotations around the -axis implies that dynamics of entanglement is identical for states related through such rotations. An additional simplification comes from the fact that preserves the total number of particles. We use this circumstance restricting the main analysis to the case when states are characterized by a definite number of particles, that is, when the initial state is spin coherent states . Expansion in terms of creation operators yields where in the last equality we have used spherical coordinates for spin coherent states and .

In order to find the transversal component , we first notice that solutions of the operator equations of motion can be presented as The effect of the term on the initial state while calculating reduces to a rotation of the vector at the right end. Thus, we obtain and or in the spherical coordinates . It follows from (40) that for or the total angular momentum oscillates with the period . Instants with integer correspond to minimal and, hence, maximal entanglement reached from the particular initial state.

In order to better understand the effect of the number of particles on the time evolution of entanglement we consider the case , that is, the case of initially disentangled states with symmetrically populated modes (). Up to rotation around -axis these are the only states (among initially disentangled ones) that yield maximum entanglement in the course of time evolution. Such states are the object of special interest in our consideration and, therefore, it is worth noting that they constitute a manifold (defined by ) in the Fock space of all states, including superposition of states with different total number of particles.

From (40) one finds that for these states The time dependence of entanglement is particularly simple in this case. When the number of particles increases, entanglement considerably changes (practically from to and back) within vicinities of with . When the number of particles increases these regions narrow down. In the limit one can approximate Thus, away from the points where , it can be regarded as a train of Gaussian bumps. As a result, when one can consider the system as spending most of the time in states with high entanglement, .

The Gaussian decay of is similar to the Gaussian decay of coherence of central system [53, 54]: two spins coupled to the bath. In particular, the same dependence of the decay rate on the number of particles in the environment [notice in (42)] should be emphasized. There are, however, two important differences between this situation and our case. First, the Gaussian decay for the case of central system appears in the limit of slow dynamics of the environment. In the opposite limit, the decay follows the Lorentzian law and in the intermediate case both types of decays present at different time scales [55]. For the two-mode boson field it is meaningless to separate particular particles and environment due to indistinguishability; however, it is worth noting that does not depend on the single-particle energies. The second important difference is that the concurrence of the two-mode boson field exhibits oscillations while the loss of coherence of the central system is irreversible.

Despite the complex structure of the manifold of all completely entangled states, the states reached in the course of evolution of spin coherent states are simple and belong to the class of NOON states. As a result, we can describe these states completely despite the circumstance that the general problem of the state structure refers to the full hierarchy of many-particle density matrices and therefore is not suited for the one-particle approach developed in Section 2 (see, e.g., [56]).

First we notice that due to the structure of maximally entangled states is determined solely by a nonlinear term in the Hamiltonian, which we consider in detail and we discuss the effect of the linear terms later. In order to present the state as a NOON state we rewrite where . A rationale behind such reformulation, which allows one to treat cases of even and odd on the equal ground [51, 57], will be provided in the next subsection.

The operator is represented on the basis of population numbers as an integer number . Hence, the action of the last operator on the right in (43) is reduced to a multiplication by a phase factor , a function of an integer variable with the period . Thus, using the discrete Fourier transform, we obtain where .

Expressing in terms of and , one can see that at the initial spin coherent state turns into a superposition of two spin coherent states related to each other through rotation by around the -axis. For example, if the initial state belongs to then the coordinate system can be chosen in such a way so that the maximally entangled state can be written as , or, choosing the -axis as the new quantization axis, as . The particular orientation of the NOON state and the phase are determined by the terms linear in in (36) and (43) and by the number of particles.

Two important circumstances should be emphasized with this regard. First, while the time dependence of entanglement is a periodic function, however, the time dependence of the state, generally speaking, is not periodic. For example, at entanglement vanishes, and therefore the system returns to a spin coherent state. The form of the state, however, is not necessarily the same as of the initial state but is given by that is, it is rotated around the -axis. Second, the angle of rotation depends on the number of particles, unless with an integer . As a result, states with a different number of particles pass through different sequences of spin coherent states.

This suggests that the time evolution of entanglement of states, which cannot be characterized by a definite number of particles, demonstrates features peculiar from the perspective of the dynamics of the spin coherent states. A general form of a disentangled state is given by (12). Using that we find where is the average taken within the -particle sector. Despite an independent contribution of each sector into the average components of the angular momentum, their effect on entanglement is more complex due to . For concreteness, we limit ourselves to considering a case of special interest when the initial state is disentangled within both field and particle pictures. That is is Glauber’s coherent state and, hence, , where is the average number of particles. Moreover, we restrict ourselves to the case when and obtain for the normalized angular momentum where .

First of all, it should be noted that the maximal value of entanglement is limited from above by , where . Thus, Glauber’s coherent states do not pass through completely entangled states. However, the “dephasing" of dynamics within sectors corresponding to a different total number of particles brings interesting features to the time dependence of entanglement.

It follows from (47) that, unless is a rational number, entanglement vanishes only at . In other words, for a generic interaction potential, the loss of coherence of Glauber’s coherent states is irreversible (while, of course, nonmonotonous). The initial loss is described by a Gaussian law with . This expression demonstrates that for the case of superposition of states with different numbers of particles, one has two cooperating channels of entanglement: the internal dynamics within each sector and the dephasing between different sectors.

Another manifestation of existence of an additional channel is the formation of entangled states even if , when spin coherent states remain disentangled. In this case the variation of entanglement is due to , that is, due to dephasing of states with different number of particles. It follows from (47) that the maximal entanglement in this case is reached at with . It can be seen that this corresponds to the “antiferromagnetic” arrangement of spin coherent states with different total number of particles (see Figure 1). The one-particle states in sectors with the numbers of particles differing by 1 are characterized by opposite transversal components of the average angular momentum.

Clearly, the same result holds when but . The instant when NOON states are formed for each component with the definite number of particles (when ) does not correspond to maximal entanglement because of the uncompensated angular momentum of the one-particle component. In fact, is reached at , when again one has spin coherent states within different sectors but states with odd numbers of particles are oppositely oriented, compared to orientation determined by the linear dynamics.

It follows from this consideration that when both channels contribute, in order to reach it is necessary to have commensurate and . Thus, for a generic interaction, potential states with maximal entanglement are inaccessible and for the inequality holds. The overall dependence of on and has a simple structure (see Figure 2). In order to describe its main features we introduce and consider the normalized angular momentum as a function of two variables , where . The function has a lattice of resonances , where and are integers of the same parity. The resonances are of the Gaussian form with the principal axes along vectors and on the -plane and the respective decay parameters and , where and we have neglected corrections of the order .

We would like to emphasize that the reason why Glauber’s coherent states cannot be completely entangled is that the dynamics governed by Hamiltonian (36) is unable to cancel the contribution of the one-particle component when the distribution of amplitudes corresponds to Glauber’s coherent state. If one discards the requirement for the initial state to be disentangled in the field picture, it is not difficult to construct such sets of , which produce states evolving through completely entangled states of various structures.

##### 4.2. -Model

A general case of -dimensional one-particle Hilbert space can be treated using the same ideas as those that have been applied in the previous subsection. As has been discussed in Section 2, the most convenient framework for studying entanglement in the particle picture is provided by the language of Lie algebras of Lie groups describing transformations of the basis in one-particle Hilbert space. From this perspective it would be, perhaps, more natural to approach the problem in the spirit of studies of dynamical symmetry and to present the Hamiltonian in a form emphasizing the role of as a spectrum generating algebra [58], for example, in terms of Cartan generators (the elements of the Cartan subalgebra) of . We, however, find it more convenient to assume that a basis in the Hilbert space is chosen in such a way that the Hamiltonian can be written solely in terms of the operators of the number of particles in individual modes Hamiltonians of this form emerge from a general Hamiltonian with two-particle interaction, when there are linearly independent one-particle integrals of motion. This corresponds to an off-resonant dynamics of many-level atoms in a cavity [59]. A similar situation may take place in highly symmetrical molecular systems [60, 61] and the special interest the model presents in the context of the Bose-Hubbard model [62], where it corresponds to the strong coupling limit and describes fast on-site dynamics.

The structure of Hamiltonian (48) allows one to find the evolution of time with a general ladder operator where is the order of the operator. Indeed, noticing that we find the solution of in the form where , and .

For the analysis of one-particle entanglement, we need to consider a particular case . Its time dependence is given by (51) with Initially, a disentangled state with a definite number of particles is a spin coherent state . Finding follows the same steps as in the case. Operator acting on turns it into another spin coherent state with . Thus, we have for the OPCM This yields where and . The term is an integral of motion and is determined by the population of one-particle states in the initial state. It can be expressed in terms of the normalized average values of Cartan generators , where ; obviously this expression does not depend on a particular choice of the defining representation of the Cartan generators. Thus, plays a similar role as for the -model and limits the maximal value of entanglement, which can be reached for a particular initial state. The states of special interest are those that admit complete entanglement. It can be seen that these are states with for . In more general terms, these are states for which average values of all Cartan generators vanish. This condition defines a manifold in the Fock space. Thus, for a spin coherent state from we have .

The variation of entanglement with time is related to . It can be expressed in terms of the “transversal” component of the Bloch vector with , where the summation runs over the ladder operators of in the Hermitian representation. Using the expressions found for the off-diagonal elements of OPCM we find In particular, we find that for complete entanglement is reached at instants when vanishes for all . It is not difficult to see that not all potentials admit the existence of such instants. Moreover, one can show that in striking contrast to the model such potentials constitute a special and not very rich class.

In order to see this we first notice that can be seen as matrix elements of a matrix , where . The condition of complete entanglement in these terms is simply . Thus, must be a symmetric (due to ) unitary matrix. Next, we notice that the condition remains invariant under the symmetry transformation for an arbitrary set of real . This transformation can be used to turn all elements of the first row (and the first column) of to zeros. Finally, we notice that up to permutations of coordinates of the orthonormal vectors constituting the matrix or, equivalently, up to permutations of the vectors themselves, symmetric unitary matrix with elements of the first row equal to is unique since it is determined by phase factors satisfying inhomogeneous equations. The matrix elements of such matrix are given by Thus, after a proper enumeration of the basis vectors in one-particle Hilbert space, potentials admitting the complete entanglement are of the form with . It follows from this expression that for the model the manifold of such potentials is parametrized by parameters, while the dimension of the manifold of symmetric matrices is . For these dimensions coincide and as a result any interaction potential causes complete entanglement (of states with ) except for , which corresponds to with the factor stemming from the difference in definitions of and and in (35). In the case, the manifold of completely entangling potentials has codimension two and so on. Thus, by increasing the dimensionality of the one-particle Hilbert space completely, entangling potentials are rather exceptional, at least when dynamics preserves populations of individual modes.

The time dependence when increases becomes more complex with a characteristic self-similar pattern emerging (see Figure 3(a)). Its overall form can be found using (59) in (56) for . In the , the time dependence of has a typical form of strong (when ) “resonances” at instants , where the one-particle coherence is briefly restored. In the case, additional instants of partial recovery of coherence emerge. By performing summation over in (56) we obtain When the term under the sum is small unless . Thus, we have a partial recovery of coherence, when with coprime and . It follows then that the number of such points is equal to the number of irreducible fractions of the form and the degree of recovery (the magnitude ) is proportional to the number of pairs satisfying .

Of more general interest is the question of initial formation of entanglement. In the limit , we obtain for small Thus, one has a spectrum of entanglement times (the same holds also for points of partial recovery of coherence). Estimating the entanglement time as the longest one we obtain .

The distribution of instants when complete entanglement is reached is also more complex compared to the case. Within the period , one has multiple instants where . We denote that at the earliest time when the complete entanglement is reached. One can see from (60) that , where with integral such that . From the perspective of determining potentials admitting complete entanglement, this corresponds to permutations of either coordinates of the vectors constituting symmetric unitary matrix or the vectors themselves.

Thus, the number of points , when , depends on prime factorization of . For example, if is prime, for all . For composite ’s, some of these points are excluded and one has the least density of points of complete entanglement when is a factorial. It should be noted, however, that and always correspond to complete entanglement (see Figure 3(b)).

Finally, we turn to the question of the structure of completely entangled states obtained in the course of evolution of initial spin coherent states. The same approach used for the model can be applied here. First, we present the potential as in (59) with , where is the earliest time when a completely entangled state is reached, and study the structure of , where and . Similar to what we had in the previous subsection, the structure of the entangled state is determined by the term , while the term affects only the orientation of the state. Therefore, in order to analyze the structure of the state it suffices to consider . The idea is the same as for the -model. We expand over basis states on the basis of population numbers and use the fact that action of on these states is equivalent to multiplication by an integer number . Thus, can be regarded as a function of an integer variable and with the help of the discrete Fourier transform can be expanded in terms of linear operators in .

The implementation of this Hubbard-Stratonovich type of transformation is slightly different for even and odd , and therefore we consider these cases separately.

For even the period of is and we obtain where Here we have used which follows from Siegel’s reciprocity relation for generalized quadratic Gauss sums [63, 64] and is valid when is even.

When is odd we need to consider instead. Here the second exponential is a periodic function with the period . Thus, using (64), we find where .

Applying (62) and (65) to the case when we can see that the structure of completely entangled states is essentially the same as in the case: they are superpositions of spin coherent states characterized by state vectors , where and . It can be easily seen that . In other words, a completely entangled state is built of spin coherent states corresponding to an orthogonal set of one-particle states. Thus, completely entangled states belong to the class of NOON states.

The structure of completely entangled states reached at with differs only by precise relations between phases of superimposed spin coherent states. In order to see this we consider the more transparent case of even . The expansion of is given by (62) with amplitudes where is defined by or . In (66) only the phase factor in front of the sum is relevant and it describes the difference between completely entangled states reached at different .

The obtained results can be generalized to the case of superpositions of states with different total numbers of particles following the same line as for the model. For an initially disentangled state, presented as in (12), we have Thus, one can expect a qualitative similarity between and models. Let the initial state be Glauber’s coherent state characterized by the one-particle state and the average total number of particles . Then, the time dependence of entanglement is given by (54) with First of all it should be noted that ; thus complete entanglement cannot be achieved for any interaction potential. Instead one has to consider the problem of the minimal value that can reach for given potential and initial state. This complex problem is significantly simplified for the case when the interaction potential and the initial state admit complete entanglement within sectors with a definite number of particles, that is, when potential can be presented as in (59) and . In this case we have where . This equation has the same structure as (47) and posseses the same basic features related to the existence of different channels of entanglement in the case of superposition of states with different number of particles. One channel is through the formation of NOON states as has been discussed above and another channel is related to dephasing of states belonging to sectors with a different total number of particles. As well as in the case, initially, at , the effects of these two mechanisms combine leading to faster entanglement compared to the case of spin coherent state.

When there is an interesting dependence of the general profile of on the parity of . In order to see it, we consider the case . The magnitude of in this case is determined by the largest term in the sum over and . Therefore, a special attention should be paid to vicinities of , instants of partial recovery of coherence. The respective contributions to have the form with integral . Thus, when is odd all correspond to partial recovery of coherence. While for even , some of these points in fact correspond to increased entanglement, an effect similar to what we have seen in the case (see Figure 4).

When both entangling channels contribute (i.e., we have and ), the initial time dependence of entanglement is relatively simple: the contribution of terms with vanishes at time scales ~, leading to a reduction of entanglement oscillations. The long time behavior, however, is quite complex and strongly depends on the commensurability of and .

#### 5. Entanglement Dynamics in Two-Mode Jaynes-Cummings Model

In the previous sections we have considered only interaction between boson fields. Here, we extend the analysis for the case when the boson field interacts with the system, whose dynamics is subject to the Pauli exclusion principle. More specifically, we consider the evolution of entanglement of photons in a cavity with quantum dots (QD) modelled by two-level systems admitting transitions characterized by definite helicities.

The excitations of QDs are described by , operators creating electron with spin in the QD at point , and , operators creating holes. A dot in the ground state is excited by creating an electron-hole pair. This is described by acting on the QD’s ground state by . The helicity of the respective transition is given by . The operators dealing with transitions of fixed helicity constitute algebras and, therefore, it is more convenient to introduce the respective generators: , , and so on, where we have assumed that the states in the valence bands are heavy-holes.

Using these notations, we present the Hamiltonian governing the dynamics of the system as a Hamiltonian of two-mode Jaynes-Cummings model [65, 66]. Consider where runs over QDs and are the excitation energies of QDs and . In terms of the electron and hole operators, are written as . The interaction between photons and QDs is given by where are the respective Rabi frequencies. Generally speaking, they are complex numbers depending on the phase of the photon mode. Their argument, however, is of minor importance and, therefore, in what follows we will assume that .

We do not require the parameters of the QDs to be the same. This accounts for the variation of QD’s excitation energies due to the size fluctuations and the variation of the Rabi frequencies due to position of the QDs in the cavity. The dynamics of entanglement in this system is a result of interference of harmonics with incommensurate frequencies [67] and in the long time limit the time dependence of entanglement is very complex. However, when the number of photons is large, then for relatively short time scales the evolution of entanglement demonstrates universal features, which we will analyse in detail. More specifically, we consider the evolution of entanglement for the case when initially we have a disentangled photon state in the cavity and all QDs are in the ground state. This closely corresponds to the situation when, for example, the cavity in the ground state is pumped by an external source.

It is convenient to take into account that the total number of excitations remains constant: , where . The effect of the operator on the dynamics is simple and reduces to rotation of the initial photon state around the -axis. It results in acquiring phase factors by and, therefore, does not affect the entanglement. This circumstance allows us to exclude the single-photon contribution to the energy and to introduce , where is the detuning of the -transition in the QD at from the resonance with the cavity mode.

The nontrivial part of the dynamics is due to the excitation of a QD by absorbing a photon. Let the initial state of the system be characterized by definite numbers of photons with each polarization, , where corresponds to all QDs being their ground states. Then at any instant we can present where and are the QDs state with QDs in excited -states and describe the distribution of the excited states over the QDs. The matrix element of the Hamiltonian between the respective states is proportional, for example, to , which varies in the range from to . In the limit , we can neglect this variation, which is sufficient for the consideration of relatively short time scales, as will be elaborated below. Moreover, if one is interested, for example, only in photon observables (say, the photon OPDM), one can extend this approximation further and neglect the variation of photon states in expansion (72) and set up to terms . Thus, treating consistently the limit one can approximate by replacing photon creation and annihilation operators by , similar to the Bogoliubov approximation in the theory of excitations of the Bose-Einstein condensates. The validity of this approximation for initially disentangled states, however, is not self-evident. Indeed, except the cases when or , disentangled states are not characterized by a definite numbers of photons. Therefore, generally speaking, the condition is not fulfilled for disentangled states. On the other hand, states with a significantly mismatched number of photons with different polarizations are typical for states with , and the time evolution of entanglement of such states is rather simple due to the limitation on the maximum value of entanglement. Similar to previous sections, states of our main interest are those which pass through completely entangled states. These are states with and for simplicity we will consider states with . The contribution of states with mismatched and decreases with fast , where . Thus, if , the contribution of the states, for which the condition does not hold, is negligible.

Adopting these approximations we can write the Hamiltonian in a simple form where with , and so on, and . Hamiltonian (73) has the form of a Hamiltonian of a lattice of spins in a magnetic field depending on the number of photons. It should be noted that “spins" are formed by the generators of the algebras corresponding to each helicity independently.

Different terms in (73) commute and, therefore, the operator of evolution is represented as a product of operators acting on individual QDs. Thus, by expanding the photon state over the basis of population numbers one can evaluate the action of the operator of evolution on the initial state and find the resulting photon OPDM. The final expression turns out to be very cumbersome due to noncommutativity of operators taken for photon states with different . It can be seen, however, that the variation of the directions of leads to corrections to the state