Information diagrams in the study of entanglement in symmetric multi-quDit systems and applications to quantum phase transitions in Lipkin-Meshkov-Glick D-level atom models

In this paper we pursue the use of information measures (in particular, information diagrams) for the study of entanglement in symmetric multi-quDit systems. We use generalizations to U(D) of spin U(2) coherent states and their adaptation to parity (multicomponent Schr\"odinger cats) and we analyse one- and two-quDit reduced density matrices. We use these correlation measures to characterize quantum phase transitions occurring in Lipkin-Meshkov-Glick models of D=3-level identical atoms and we propose the rank of the corresponding reduced density matrix as a discrete order parameter.


I. INTRODUCTION
Information diagrams were introduced to discuss the relation between two different information measures, like von Neumann entropy and error probability [1], or von Neumann and linear entropies [2]. In the particular case of linear (L) and von Neumann (S) entropies, pairs (L(ρ), S(ρ)) are usually plotted for any valid probability distribution ρ.
Here ρ can also represent the density matrix of a quantum system (or rather a vector with its eigenvalues), and this is our main interest in this paper. Special attention is paid to the boundaries of the resulting information diagram region, where the associated probability distributions (or density matrices) will be denoted as "extremal". In Ref. [3], a comparison is made between both entropies in the case of two qubits (see also [4] for the case of the ion-laser interaction). In [5], a detailed study of information diagrams is carried out for arbitrary pairs of entropies. There it is proved that, for certain conditions (satisfied by linear, von Neumann and Rényi entropies), the extremal density matrices are always the same. Counterexamples are given but, in general, the deviation will be very small and we can safely assume that these extremal density matrices have universal character.
In this paper we shall use information diagrams to obtain global qualitative information of particle entanglement in symmetric multi-quDit systems described by generalized "Schrödinger cat" (multicomponent DCAT) states (first introduced in [6] as two-component, even and odd, states for an oscillator). These DCAT states turn out to be a Z D−1 2 parity adaptation of U(D)-spin coherent (quasi-classical) states and they have the structure of a quantum superposition of weakly-overlapping (macroscopically distinguishable) coherent wave packets with interesting quantum properties. For that purpose we make use of one-and two-quDit reduced density matrices (RDM), obtained by extracting one or two particles/atoms from a composite system of N identical quDits described by a cat state, and tracing out the remaining system. It is well known (see [3] and references therein) that the entropy of these RDMs provides information about the entanglement of the system. We shall plot the information diagrams associated to these RDMs and extract qualitative information about one-and two-quDit entanglement, and also about the rank of the corresponding RDM, which also provides information on the entanglement of the original system [7].
We shall apply these results to the characterization of quantum phase transitions (QPT) occurring in Lipkin-Meshkov-Glick models of 3-level identical atoms, complementing the results of [8]. In particular, we have seen that the rank of the one-and two-quDit RDMs can be considered as a discrete order parameter precursor detecting the existence of QPTs.
The paper is organized as follows. Section II reviews the notion of information diagram, describing its main properties, particularly with respect to the rank. Section III reviews the concept of U(D)-spin coherent states and their Z D−1 2 parity adapted version, the DCAT. In Section IV we compute one-and two-quDit RDMs for the 2CAT and the 3CAT, their Linear and von Neumann entropies, plotting them and constructing the associated information diagrams. In Section V we use information diagrams to provide qualitative information about QPTs in Lipkin-Meshkov-Glick (LMG) models. Section VI is devoted to conclusions.

II. INFORMATION DIAGRAMS
To determine the boundaries in information diagrams [2] we need to show that, for two different measures of entropy (or information) E 1 and E 2 , there are maximum and minimum possible values of E 1 (resp. E 2 ) for a given value of E 2 (resp. E 1 ) [5]. That is, the region ∆ given by the image of the map ρ → (E 1 (ρ), E 2 (ρ)) is a bounded set in the plane, where ρ denotes all possible probability distributions (or density matrices) for a given dimension d.
Since usual measures of entropy for density matrices are based on the trace, they are invariant under changes of basis. Hence, the only relevant information of a density matrix is contained in its eigenvalues, thus in this paper we shall identify density matrices ρ with their eigenvalues (λ 1 , λ 2 , . . . , λ d ), the order being irrelevant. Therefore, for our purposes, we can identify probability distributions and density matrices using a vector notation in terms of eigenvalues, referring to both of then as density matrices for short. Notwithstanding, we shall continue to treat density matrices as matrices in some situations.
In [5] it was proved that, under rather general assumptions on the convexity/concavity of the entropy measures, the maximum and minimum values are always attained in two standard forms of density matrices respectively, where k = 1, . . . , d − 1. Let us write the previous equations as (convex) sums of density matrices, that in turn can be seen as lower dimensional density matrices. For that purpose denote by ρ k the maximally mixed density matrix (or equal probabilities distribution) in dimension k, ρ k = ( 1 k , (k) . . ., 1 k ) = 1 k I k , where I k is the identity matrix in dimension k. Then we have: where 0 k is the null matrix (or vector) in dimension k and k = 1, . . . , d − 1. The relation between and λ is for eqns. (1,3) and λ = 1− k for eqns. (2,4). In most cases, the pair of entropies (L, S) is considered, where L and S denote linear and von Neumann entropies, respectively. We shall consider here normalized linear and von Neumann entropies, i.e: in such a way that both entropies range from 0 (pure states) to 1 (maximally mixed states). The values of both entropies for each family of curves are: and In Figure 1a the curves ρ → (L(ρ), S(ρ)) are shown for ρ equal to ρ max ( ) and ρ   min ( ) for k = 1, . . . , d − 1, which are traced from right to left. Note that in the case k = 1 the associated curve is the same as in (a), but traced backwards. Also, for k = d − 1 the associated curve coincides with that of ρ max . (c) Plot of the asymptotic curves (8,9) for density matrices near a pure state (bottom-left, red and pink, respectively) and the asymptotic curve (10) near the maximally mixed state (upper-right, green) Note that the density matrices (3) can be seen (for small ) as the maximally mixed density matrix ρ d perturbed by a rank-1 density matrix, while those of (4) can be seen as maximally mixed density matrix of dimension k, ρ k , perturbed by a (orthogonal) rank-1 density matrix, for k = 1, . . . , d − 1.
It should be stressed that the range of the parameter in the curves ρ  min . Each subregion ∆ k contains density matrices of rank greater than k. Density matrices of rank 1 (pure states) lie on the origin, while density matrices of rank 2 lie on the curve ρ (1) min =ρ (1) min . See Figure 2 for a plot of a sample of 20000 density matrices of dimension d = 5 randomly generated following a χ 2 distribution for the eigenvalues where the color of the corresponding point in the information diagram is associated to its rank (warmer colors correspond to higher rank).
From the expression of the extremal density matrices (1,2), or their alternative expressions (3,4), and the expression of the inner curvesρ (k) min ( ), it is clear that, for a given value of the linear entropy and a fixed rank k + 1, the extreme values of the von Neumann entropy are reached for k identical eigenvalues. If the remaining eigenvalue is larger than the rest (i.e. we are inρ min : respectively. Near the maximally mixed density matrix (upper right of the information diagram), both ρ max and ρ (d−1) min collapse into the same curve, with asymptotic expression: See Figure 1c for a plot of these asymptotic curves in an information diagram with d = 5.
Once we have explained what the information diagrams are, and their main features, we shall use them in the study of one-and two-quDit entanglement of generalized Schrödinger cat states, which arise as a parity adaptation of U(D)-spin (symmetric multi-quDit) coherent states.

III. U(D)-SPIN COHERENT STATES AND THEIR ADAPTATION TO PARITY IN SYMMETRIC MULTI-QUDIT SYSTEMS
In this section we introduce the main ingredients and notation required to define parity adapted U(D)-spin coherent states in symmetric multi-quDit systems. These kind of states where introduced long ago in [6] as nonclassical (even and odd) states of light. We shall particularize to D = 2 and D = 3 for practical cases. See [8] for a more detailed study of the general case. We consider a system of N identical (indistinguishable) quDits, namely, D-level identical atoms. Denoting by a † i (resp. a i ) the creation (resp. annihilation) operator of an atom in the i-th level (namely, i = 1, 2 for ground and excited -or spin up and down-in the case D = 2, or i = 1, 2, 3 for a 3-level atom in the case of D = 3), the collective U(D)-spin operators can be expressed (in the fully symmetric representation) as bilinear products of creation and annihilation operators as (Schwinger representation) which generate the unitary symmetry U(D). The operator S ii represents the number of quDits in the level i, whereas S ij , i = j are raising and lowering (tunneling) operators. The fully symmetric representation space of U(D) is embedded into Fock space, with Bose-Einstein-Fock basis (| 0 denotes the Fock vacuum) when fixing n 1 + · · · + n D = N (the linear Casimir C 1 = S 11 + · · · + S DD ) to the total number N of quDits. Collective U(D)-spin operator (11) matrix elements are given by The expansion of a general symmetric N -particle state ψ in the Fock basis will be written as where is a shorthand for the restricted sum. Among all symmetric multi-quDit states, we shall pay special attention to U(D)-spin coherent states (DSCSs for short), which adopt the multinomial form 1 and are labeled by complex points z = (z 2 , . . . , z D ) ∈ C D−1 . These DSCSs can be seen as Bose-Einstein condensates (BECs) of D modes, generalizing the spin U(2) (binomial) coherent states of two modes introduced by [9] and [10] long ago. If we order levels i = 1, . . . , D from lower to higher energies, the state |z = 0 would be the ground state, whereas general |z could be seen as coherent excitations. Coherent states are sometimes called "quasi-classical" states and we shall see in Section V that |z turns out to be a good variational state that reproduces the energy and wave function of the ground state of multilevel LMG atom models in the thermodynamic (classical) limit N → ∞.
Expanding the multinomial (15), we identify the coefficients c n of the expansion (14) of the DSCS |z in the Fock basis as where we have written |z| ≡ (z · z) 1/2 = (1 + D i=2 |z i | 2 ) 1/2 for the "length" of z. Note that DSCS are not orthogonal (in general) since is not zero, in general. However, contrary to the standard (canonical, harmonic oscillator) CSs, they can be orthogonal when z · z = 0.
In [8] we have shown that DSCSs are separable and exhibit no quDit entanglement (although they do exhibit interlevel entanglement). In fact they can be written as a tensor product of 1-quDit coherent states: where we added the superscript (N ) to the N -particle coherent state (15), and |z i denotes the one-particle coherent state for the i-th quDit. Note that this state is explicitly symmetric under the interchange of quDits and therefore there is no need to symmetrize it.
The situation changes when we deal with parity adapted DSCSs, sometimes called "Schrödinger cat states", since they are a quantum superposition of weakly-overlapping (macroscopically distinguishable) quasi-classical coherent wave packets. These kind of cat states arise in several physical situations and display interesting nonclassical properties. The case of even parity cat states is particularly important since they turn out to be good variational states [10], reproducing the energy of the ground state of quantum critical models in the thermodynamic limit N → ∞. In [8], the even parity multi-quDit cat state DCAT have been constructed for general D, and here we shall reproduce the construction to fix notation.
The parity operators are defined as Note that Π −1 i = Π i and Π 1 . . . Π D = (−1) N , a constraint that says that the parity group for symmetric quDits is not instead. Therefore, we can discard in our discussion one of the parity operators, and we select Π 1 (since we will use level 1 as reference level in Sec. V).
Parity operators are conserved when the Hamiltonian scatters pairs of particles conserving the parity of the population n j in each level j = 1, . . . , D, like in the D-level LMG model considered in Sec. V. Using the multinomial expansion (15), it is easy to see that the effect of parity operators on symmetric DSCSs |z is then The projector onto the even parity subspace is given by: where the binary string b = (b 2 , . . . , b D ) ∈ {0, 1} D−1 labels the elements of the parity group Z D−1
Let us define the even parity generalized Schrödinger cat state where |z b ≡ |(−1) b2 z 2 , . . . , (−1) b D z D and we are using b as a shorthand for b∈{0,1} D−1 . The DCAT is just the projection of a DSCS onto the even parity subspace. The normalization factor is given by We shall also use the alternative notation As an illustration, let us provide the particular expressions of |DCAT for D = 2 and D = 3. Denoting by |z = |z 2 = |α the coherent state (15) for D = 2, the corresponding even parity 2CAT state is given by with normalization factor with L ± = 1 ± |α| 2 . Note that the overlap α| − α = (L − /L + ) N N →∞ −→ 0 for α = 0, which means that |α and | − α are macroscopically distinguishable wave packets for any α = 0 (they are orthogonal for |α| = 1).

IV. ENTROPIC MEASURES ON REDUCED DENSITY MATRICES TO QUANTIFY ENTANGLEMENT
One of the most important applications of entropy measures is to quantify the entanglement of the state of a system. For that purpose we define several types of bipartition of the whole system, computing the corresponding RDMs and entanglement measures for symmetric multi-quDit states ψ in terms of linear L and von Neumann S entropies. We shall focus on one-and two-quDit entanglement, computing the one-and two-particle RDMs (ρ 1 and ρ 2 ) for a single and a pair of particles extracted at random from a symmetric N -quDit state ψ. The procedure is straightforwardly extended to ρ M for an arbitrary number M ≤ N/2 of quDits. However, as we shall see, it is not necessary to go beyond two particles since the two-particle RDMs provides enough information for small values of D. Actually, in the particular case of D = 2, the one-particle RDM contains all necessary information about the entanglement of the system.
In [8] we gave the general expression of the one-quDit RDM of any normalized symmetric N -quDit state ψ like (14), expressed in terms of expectation values of U(D)-spin operators S ij as where E ij represent D 2 , D ×D-matrices with entries (E ij ) lk = δ il δ jk (1 in row i, column j, and 0 elsewhere). Likewise, the two-particle RDM of a symmetric state ψ of N > 2 quDits is written as [8] The matrices E ij are the generalization to arbitrary D of standard Pauli matrices for qubits (D = 2), namely E 12 = σ + , E 21 = σ − , E 11 − E 22 = σ 3 and E 11 + E 22 = σ 0 (the 2 × 2 identity matrix). Actually, the one-and two-qubit RDMs for D = 2 were already considered time ago by Wang and Mølmer in [15]. Here we shall consider both cases, D = 2 (qubits) and D = 3 (qutrits), in order to discuss the similitudes and differences.
A. One-quDit reduced density matrices For the case of a DSCS |z , the linear and von Neumann entropies of ρ 1 (z) are zero, i.e. there is no entanglement between quDits in a DSCS. This is because a DSCS is eventually obtained by rotating each quDit individually.
The situation changes when we deal with parity adapted DSCSs or "Schrödinger cat states" like (25)- (26). Indeed, the one-quDit RDM ρ 1 (DCAT) does not correspond now to a pure state since it has the expression (we provide its eigenvalues) for an N -qubit system and for an N -qutrit system. Note that, for α = 0 in the case of ρ N 1 (2CAT), and α = 0 or β = 0 in the case ρ N 1 (3CAT), the corresponding one-quDit RDM has rank greater that 1. That is, unlike |z , the Schrödinger cat |DCAT is not separable in the tensor product Hilbert space [C D ] ⊗N . In addition, ρ N 1 (3CAT) has rank 2 if α = 0 or β = 0 and has rank 3 if both are different from zero. See below for a more detailed discussion on this point.
In Figures 3a and 3b we represent contour plots of linear and von Neumann entanglement entropies in the limit N → ∞ of the one-qutrit RDM, ρ ∞ 1 (3CAT), of the 3CAT in Eq. (26), as a function of the phase-space CP 2 coordinates (α, β) [actually, they just depend on the moduli]. Both entropies are again normalized to 1. They attain their maximum value of 1 at the phase-space point (α, β) = (1, 1) corresponding to a maximally mixed RDM. This behavior of the entropies, and therefore of entanglement (together with squeezing, see [8]) parallels that of the standard harmonic oscillator cat states where the maximun entanglement and squeezing takes place for relatively small values of the coherent state parameter [16]. The diference here in the D = 3 case is that for large values of the parameter there can still be entanglement (and squeezing), dependent on the angle of the directional limit (see eq. (37)). These Figures also show (in magenta color) the values of the entropies along the stationary curve (α(λ), β(λ)) in Eq. (49), that we already mentioned before the Eq. (34). For high interactions λ → ∞ we have (α(λ), β(λ)) → (1, 1), which means that highly coupled quDits are maximally entangled in a cat-like ground state (we shall come back again to this discussion later in Section V). In Figures 3c and 3d the asymptotic behavior for large |α| and β| is shown, where contours of linear and von Neumann entropies coincide with the (isentropic) lines θ =constant, according to the asymptotic behavior of ρ ∞ 1 (3CAT) in (37). In Figure 4a we plot the information diagram for the family of 1-qutrit RDMs for a 3CAT (31) in the limit N → ∞, for all values of |α| and |β|. It can be seen that they fill completely the region ∆. Also, the stationary curve (49) is shown, starting at the origin (zero entropy and therefore no entanglement), moving on the curve ρ (1) min and through the region ∆ 2 , to finish at the maximally mixed state, indicating that this state is maximally entangled.
In Figure 5a-5b, we represent contour plots of normalized linear and von Neumann pairwise entanglement entropies in the thermodynamic limit N → ∞ for the two-qutrit RDM, ρ ∞ 2 (3CAT), of a U(3) Schrödinger cat (26), as a function of the phase-space CP 2 complex coordinates (α, β) [they just depend on the moduli]. As for the one-qutrit case, they attain their maximum value at the phase-space point (α, β) = (1, 1) ("high coupling limit"); however, unlike the one-qutrit case, pairwise entanglement entropies do not attain the maximum value of 1 at this point, but L 2 = 5/6 0.833 and S 2 0.623 for large N . As already commented, variational (parity adapted spin coherent) approximations to the ground state of the LMG 3-level atom model [discussed later in Section V] recover this maximum entanglement point (α, β) = (1, 1) at high interactions λ → ∞, as can be seen in the already discussed stationary curve (49). In Figures 5c and 5d, the asymptotic behavior for large |α| and β| is shown, where contours of linear and von Neumann entropies coincide with the lines θ =constant, according to the asymptotic behavior of ρ ∞ 2 (3CAT) in (44).  (26) for N → ∞, as a function of the phase-space coordinates α, β (they just depend on moduli). The meaning of the magenta curve is the same as in the Figure 3. The asymptotic behaviour of (c) L ∞ 2 and (d) S ∞ 2 for large values of |α| and |β| displays isentropic curves θ =constant, according to the expression of ρ ∞ 2 (3CAT) in Eq. (44).
We also plot in Figure 4b the information diagram for the 2-qutrit RDM of the 3CAT in the thermodynamic limit N → ∞. It is clear that the ∆ region is not completely filled; only the subregions ∆ 2 and ∆ 3 are partially filled, the reason being that ρ 2 (3CAT) has rank 1, 2 or 4. The stationary curve (α 0 (λ), β 0 (λ)) in Eq. (49) is also shown, with a behaviour similar to the case of the 1-qutrit RDM, with the difference that it ends near the maximally mixed RDM of rank 4, more precisely, at the point 1 3 , 2 9 , 2 9 , 2 9 , 0, 0, 0, 0, 0 . It is important to notice that the stationary curve is most of the time at the inferior boundary of the set of 2-qutrit RDMs. This means that, from all 2-qutrit RDMs of the 3CAT with a given linear entropy, it has the minimum allowed value of von Neumann entropy. We conjecture that this is due to the variational character of the ground state and the universality of the extremal states lying at the boundaries of the region ∆.

V. INFORMATION DIAGRAMS AND QUANTUM PHASE TRANSITIONS IN LIPKIN-MESHKOV-GLICK MODELS OF 3-LEVEL IDENTICAL ATOMS
Now we apply the previous results to the study of QPTs of D-level Lipkin-Meshkov-Glick atom models. The standard case of D = 2 level atoms has already been studied in the literature (see e.g. [17]). We shall restrict ourselves to D = 3 level atoms for practical calculations, although the procedure can be easily extended to general D. In particular, we propose the following LMG-type Hamiltonian written in terms of collective U(3)-spin operators S ij . Hamiltonians of this kind have already been proposed in the literature [18][19][20][21][22] [see also [23] for the role of mixed symmetry sectors in QPTs of multi-quDit LMG systems]. We place levels symmetrically about i = 2, with intensive energy splitting per particle /N . For simplicity, we consider equal interactions, with coupling constant λ, for atoms in different levels, and vanishing interactions for atoms in the same level (i.e., we discard interactions of the form S ij S ji ). Therefore, H is invariant under parity transformations Π j in (19), since the interaction term scatters pairs of particles conserving the parity of the population n j in each level j = 1, . . . , D. Energy levels have good parity, the ground state being an even state. We divide the two-body interaction in (46) by the number of atom pairs N (N − 1) to make H an intensive quantity, since we are interested in the thermodynamic limit N → ∞. We shall see that parity symmetry is spontaneously broken in this limit. As already pointed long ago by Gilmore and coworkers [10,24], coherent states constitute in general a powerful tool for rigorously studying the ground state and critical properties of some physical systems in the thermodynamic limit. The energy surface associated to a Hamiltonian density H is defined in general as the coherent state expectation value of the Hamiltonian density in the thermodynamic limit. In our case, the energy surface acquires the following form where we have used the parametrization z = (α, β), as in eq. (26), for U(3)-spin coherent states |z . Note that this energy surface is invariant under α → −α and β → −β, which is a consequence of the inherent parity symmetry of the Hamiltonian (46) and the transformation (20) of |z under parity. The minimum energy is attained at the stationary (real) phase-space values α ± 0 = ±α 0 and β ± 0 = ±β 0 with In Figures 3 and 5 we plotted (in magenta color) the stationary-point curve (α 0 (λ), β 0 (λ)) on top of one-and twoqutrit entanglement entropies, noting that (α 0 (λ), β 0 (λ)) → (1, 1) for λ → ∞ (high interactions). Inserting (49) into (47) gives the ground state energy density at the thermodynamic limit Here we clearly distinguish three different phases: I, II and III, and two second-order QPTs at λ (0) I↔II = /2 and λ (0) II↔III = 3 /2, respectively, where ∂ 2 E0( ,λ) ∂λ 2 are discontinuous. In the stationary (magenta) curve (α 0 (λ), β 0 (λ)) shown in Figures 3a, 3b, 4, 5a, 5b, and 6, the phase I corresponds to the origin (α 0 , β 0 ) = (0, 0) (square point), phase II corresponds to the horizontal part β 0 = 0 up to the star point, and phase III corresponds to β 0 = 0.
Note that the ground state is fourfold degenerated in the thermodynamic limit since the four U(3)-spin coherent states |z ±± 0 = | ± α 0 , ±β 0 have the same energy density E 0 . These four coherent states are related by parity transformations and, therefore, parity symmetry is spontaneously broken in the thermodynamic limit. In order to have good variational states for finite N , to compare with numerical calculations, we have two possibilities: 1) either we use the 3CAT (26) as an ansatz for the ground state, minimizing 3CAT|H |3CAT , or 2) we restore the parity symmetry of the coherent state |α 0 , β 0 for finite N by projecting on the even parity sector. Although the first possibility offers a more accurate variational approximation to the ground state, it entails a more tedious numerical minimization than the one already obtained in (48) for N → ∞. Therefore, we shall use the second possibility which, despite being less accurate, it is straightforward and good enough for our purposes. That is, we shall use the 3CAT (26), evaluated at α = α 0 and β = β 0 , as a variational approximation |3CAT 0 to the numerical (exact) ground state |ψ 0 for finite N .
Let us apply the tools developed in previous sections to this model and draw the main conclusions. Firstly, in Figure  6a and 6b, we have added to the information diagrams for 1 and 2 qutrits RDMs already shown in Figure 4, the curves (as a function of λ) of the numerically computed ground states of the 3-level LMG model for different values on N (in green colors), together with the already shown analytical variational curve (in magenta) (α 0 (λ), β 0 (λ)) for N → ∞. We can conclude that they do not lie in the inferior part of the region ∆, as the variational one, but as N grows the numerical curves approach the analytical one. Secondly, suggested by the results about the rank of 1 and 2 quDits RDMs of Section IV, we plot in Figure 7 the rank of the RDMs as a function of λ for both variational (N → ∞) and numerical (N = 50) solutions for the ground state of Hamiltonian (46). The QPT critical points λ (0) I↔II = /2 and λ (0) II↔III = 3 /2, are clearly marked in the case of the variational curve, with a jump from rank 1 to rank 2 at λ = /2 and another jump from rank 2 to rank 4 (3 in the case of 1 qutrit RDMs) at λ = 3 /2. In the case of the numerical curve, where a small threshold has been applied to the eigenvalues to suppress spurious oscillations, the first jump continues to be at λ 0.5, whereas the second jump takes place at slightly larges values of λ = 1.5 (in = 1 units). This behaviour is the same as with other precursors of QPTs, like susceptibility of fidelity in the 3-level LMG model [23].
From this, it is clear that the rank of the RDMs is a good precursor of a QPT, with the advantage of being a discrete parameter.

VI. CONCLUSIONS
In this paper we have used an information-theoretic tool like the information diagrams to extract qualitative information about the quDit entanglement (and rank) of parity adapted U(D)-spin coherent states (DCATs) using one-and two-quDit reduced density matrices, and we have applied it to the study of atom entanglement in the ground state (both variational, in the N → ∞, and numerical, with finite N ) of the 3-level atom LMG model.
We have shown how the allowed region ∆ of information diagrams is completely filled in the case of one-qutrit RDMs, while only the lower part of it is partially filled in the case of two-qutrits RDMs. This indicates that the maximum pairwise (2-qutrit) entanglement attained in a 3CAT state is smaller that the maximum one corresponding to a maximally mixed RDM or order 3 2 . We have already seen that this maximally entangled 3CAT is attained for the values (α, β) = (1, 1) (or α = 1 for D = 2), and these are precisely the values obtained for the variational analytical approximation to the ground state of a 3-level LMG model in the high coupling regime.
In addition, we have shown that the variational curve (α 0 (λ), β 0 (λ)) practically all the time lies in the inferior part of the information diagram subregion filled by all 3CAT states. We conjecture that this is due to the variational character of these states (minimum of the energy surface (47)) and the universality character of the extremal states lying at the boundary of the region ∆.
Information diagrams also provide qualitative information about the rank of the RDMs. This has motivated us to study with detail their rank for different values of the parameters α and β of 3CAT states (see Section IV), indicating that the one-and two-quDit RDMs have in general lower ranks than the maximal rank allowed by the corresponding dimension. Focusing on the variational analytic curve (α 0 (λ), β 0 (λ)), and in the numerical solution for the ground state for finite N , Figure 7 shows that the rank of one-and two-qutrit RDMs has jumps precisely at the points where QPTs occurs (or near these values in the numerical finite N case). Therefore the rank can be used as a discrete precursor of a QPT in the LMG model, but this conclusion can be probably extended to other critical models.
All these results motivate us to further study the application of information diagrams and rank of RDMs to other parity adapted U(D)-spin coherent states, but with different parity character. Here we have restricted ourselves to the even case, but remember that there are 2 D−1 different parity adapted U(D)-spin coherent states, the even one just being a particular case. For example, odd parity cat states (for D = 2) are known to be well suited as variational states to approximate excited states in, for example, the Dicke model of superradiance [25].
Since the rank of a RDM is equal to the Schmidt number, by the Schmidt decomposition theorem (see, for instance [26]), it would be interesting to study with detail the Schmidt decomposition of parity adapted U(D)-spin coherent states (not only of the even one, but for all 2 D−1 parity invariant states) when we extract 1, 2, or in general M quDits, and find the basis realizing the Schmidt decomposition in the larger factor.