Entropy–energy inequalities for qudit states
Abstract
We establish a procedure to find the extremal density matrices for any finite Hamiltonian of a qudit system. These extremal density matrices provide an approximate description of the energy spectra of the Hamiltonian. In the case of restricting the extremal density matrices by pure states, we show that the energy spectra of the Hamiltonian is recovered for and . We conjecture that by means of this approach the energy spectra can be recovered for the Hamiltonian of an arbitrary finite qudit system. For a given qudit system Hamiltonian, we find new inequalities connecting the mean value of the Hamiltonian and the entropy of an arbitrary state. We demonstrate that these inequalities take place for both the considered extremal density matrices and generic ones.
1 Introduction
Recently [1, 2, 3, 4] an approach was established to study the ground state properties of algebraic Hamiltonians. This approach follows closely the algorithm established in [5, 6]. In particular, the approach was applied to describe the ground state of even–even nuclei within the interacting boson model [1]. In quantum optics, the procedure was used to determine the phase diagrams of the transitions between the normal regime to the superradiant behavior of the ground states of two and threelevel systems interacting with a onemode radiation field [2, 3, 4]. This approach evaluates the mean value of the Hamiltonian with respect to variational test coherent states associated to the corresponding algebraic structures of the Hamiltonian. There exists a tomographic approach, which also uses mean values of density operators in an ensemble of bases to get information on the state of the system [7, 8, 9, 10, 11]. For continuous variables, the tomographic approach has been introduced in [7, 8] in the form of optical tomography. The symplectic tomography is established in [9], and a recent review of these tomograms is given in [10]. The discrete spin tomography has been introduced in [12, 13], while the kernel for product of spin tomograms is presented in [14, 15]. The squeezed tomography is discussed in [11], which is a fair probability distribution of a discrete random variable.
One of the aims of this work is to extend the approach mentioned above to have information of the complete energy spectrum by considering the mean values of the Hamiltonian with respect to extremal density matrices [16, 17]. This is achieved by writing the mean value of the Hamiltonian as a function of the variables of a general finitedimensional density matrix [18, 19, 20, 21, 22, 23] together with the parameters of the Hamiltonian. To guarantee the positivity of the density matrix, we need to include parameters related to the purity of the density matrix [19, 22].
Another goal of this work is to obtain new inequalities connecting entropy and mean value of energy for this qudit system. We show that there exists a bound for the sum of energy and entropy determined by the partition function taken for a particular value of its argument. The method to obtain these inequalities is based on known property of positivity of the relative entropy involving two density matrices of the system states [24]. Analogous mathematical inequalities have been discussed in [25, 26]. The results obtained are valid for generic quantum states (qudits).
The main contribution of our work is to demonstrate the new approach related to the determination of the extremal points of mean values of the Hamiltonian by considering a general parametrization of the density matrices for qudit systems and to test the new entropy–energy inequality. This inequality contains the bound determined by the partition function [27]. The formulated results can be generalized to study the relations between the entropy and an arbitrary hermitian operator describing an observable.
2 Unitary parametrization of the Hamiltonian operator
The Hamiltonian operator can be expanded in terms of the set of operators that form a basis of and the identity operator as follows [19]:
(1) 
with the definitions and . The generators of satisfy the relations
(2) 
They are completely characterized by means of the commutation and anticommutation relations given in terms of the symmetric and antisymmetric structure constants of the special unitary group in dimensions [19].
In a similar form, the density matrix can be expanded, i.e.,
(3) 
because Tr, and in this case one defines
(4) 
Our purpose is to find the extreme values for the variables of the density matrix by taking the expectation value of the Hamiltonian operator. To guarantee the positivity of the density matrix, it is necessary to introduce parameters. Therefore, the extremes are obtained by means of the definition of a new function depending on variables with , Lagrange multipliers with , parameters of the Hamiltonian with , and real constants with characterizing the purity of the density matrix
(5) 
where are nonholonomic constrictions from the characteristic polynomial of , which can be obtained by means of the recursive relation [22]
(6) 
where , , and . The parameters are constants. To find the extrema, we derive the function with respect to obtaining algebraic equations regarding the independent variables of the density matrix. Then by substituting expressions (1) and (3) into (5), one arrives at
(7) 
plus differential equations regarding Lagrange multipliers
(8) 
with , , and we have used the properties of the generators of the unitary group in dimensions. These sets of algebraic equations determine the extremal values of the density matrix, i.e., and for which the expressions (7) and (8) are satisfied.
3 Extremal density matrices for and
3.1 Case
One has three generators with , which can be realized in terms of the Pauli matrices. Therefore, the density matrix can be written in the form
(9) 
and similarly an arbitrary Hamiltonian matrix is given by
(10) 
Substituting the last expressions into Eqs. (5), we obtain
(11) 
yielding, by means of expressions (7) and (8), the system of equations
(12) 
with and . Solving this system of equations, one obtains the results
(13) 
with and we defined the parameters and . Therefore, we have two solutions and substituting them into the expression for the density matrix, we obtain
(14) 
Therefore, the extremal density matrices depend on the parameter whose value is bounded, and, if the density matrix represents a pure or a mixed state, it is determined. For , one has that and the extremal density matrix takes the form
(15) 
corresponding to a mixed state with maximum entropy with the expectation value of the Hamiltonian given by .
For , one gets and the expectation values of the Hamiltonian are given by
(16) 
which corresponds exactly to the eigenvalues of the arbitrary matrix Hamiltonian in (10). The corresponding pure states are given by the density matrices in Eq. (14) replacing the value of the parameter . They are orthogonal projectors as it can be proved by multiplying the corresponding density matrices and taking the trace operation. Therefore, we have reconstructed the Hamiltonian matrix by finding the extremal values of the expectation value of the Hamiltonian.
In Fig. 1, we plot the general behavior of the expectation value of a given Hamiltonian operator as a function of the parameter . Notice that if , the mean values of the energy with respect to mixed states are close to the maximum and minimum values of the energy spectra. Additionally, we can observe that for each value of there are two solutions for the expectation values of the Hamiltonian, except when , where one gets only one expectation value associated to the mixed state with maximum entropy.
This procedure gives the eigenvalues and eigenvectors of the Hamiltonian and besides one gets information on the mean values associated to mixed states, as function of .
3.2 Case
We consider a particular Hamiltonian describing a twomode Bose–Einstein condensate [28, 29],
(17) 
where denotes the th component of the angular momentum operator. The parameter corresponds to the difference in the chemical potentials between the wells, represents the atom–atom interaction, and is related to the atom tunneling parameter.
For the qutrit case, one substitutes the threedimensional representation of the angular momentum operators, and the generators , with can be realized in terms of the GellMann matrices. Thus, an arbitrary density matrix is denoted by
(18) 
and one has a similar expression for the Hamiltonian. Comparing with the angular momentum representation, one obtains that the parameters different from zero are
(19) 
Substituting the previous expressions for the density matrix and the Hamiltonian into Eq. (5), we obtain the function in the form
(20) 
where we define . The extrema of the previous function with respect to , are obtained numerically by establishing different values of the set of parameters . Notice that the and are not independent parameters, and they are keeping track of the purity of the extremal density matrices. If we define and , then following [30], it is straightforward to get the compatible region of the parameters and . For the set , one gets the extremal density matrix, , for the mixed state with maximum entropy and .
To illustrate the procedure, we consider three different choices of the mentioned set of parameters. The obtained expectation values of the Hamiltonian are plotted in Fig. 2 as functions of , , . The pure states are associated to , and we find that there are extremal solutions for the parameters of the density matrix which are indicated by black continuous lines. Comparing these different mean values of the Hamiltonian with the corresponding exact diagonalizations, one finds a complete agreement with the energy spectra.
For and different from zero, additionally one gets information on the mean values of the energy. We obtain in this case that the number of mean values of mixed states is related to the symmetric group of dimensions; thus, for , there are colored continuous lines indicating the expectation values of the Hamiltonian of the possible mixed states.
4 Entropy and energy relations
It is known [24] that for two given density matrices and , the positivity condition for their relative entropy is given by
(21) 
The matrices and satisfy the properties , , , and . On the other hand, two arbitrary matrices which have these properties satisfy inequality (21).
We are going to use this inequality considering the Hamiltonian matrix of a qudit system, by defining the matrix
(22) 
which has the properties Tr and .
Let us write the positivity conditions of relative entropy of two matrices and
(23) 
Using the definition of the von Neumann entropy associated with the density matrix of the state, which reads , one can rewrite inequality (23) in the form
(24) 
where the mean energy is given by .
For an arbitrary Hamiltonian , one introduces the partition function
(25) 
where is interpreted as a temperature, and the matrix can be interpreted as the density matrix of the system in thermal equilibrium state. It is known that the partition function determines all thermodynamic properties of the system (cf. [27]). Then inequality (24) can be rewritten in the form
(26) 
Thus, we got for an arbitrary quantum system a bound for the sum of the mean energy value and the von Neumann entropy, and this bound is determined by the partition function evaluated in a particular value of temperature, .
For example, for a twolevel atom in the state determined by the density matrix given in Eq. (9) with the condition , and the Hamiltonian matrix given in Eq. (10), inequality (26) reads
(27) 
Notice that the last expression is invariant under unitary transformations and thus one can write , where the energy levels and are the solutions of the secular equation .
For any state, one has that the inequality (4) can be written as follows
(28) 
where we define the function in terms of and . By substituting the extremal parameters of the density matrix given in expression (13) the function takes the form
(29) 
where were defined after Eq. (13). This function is displayed in Fig. 3 as a function of the and . Notice that always it is larger than zero. The corresponding contour levels are also positive.
Similar inequalities can be obtained for any Hermitian operator, and the obtained result for the Pauli matrix with , can be written in the form
(30) 
with . For , the result is displayed in Fig. 4.
One can obtain another important inequality between the energy and the partition function by considering
(31) 
which can be written for any state in the form
(32) 
Also inequality (26) is accompanied by the inequality , which means that the partition function provides the bound for the difference of entropy and energy.
5 Summary and Conclusions
We suggested to study the properties of a qudit system considering the mean values of its Hamiltonian with respect to a generic density matrix, which includes the special cases of pure states. The number of parameters which determine the mean values of the Hamiltonian in this case is equal to the number of parameters determining the density matrix. Due to this, the information contained in the mean values of the Hamiltonian is sufficient to reconstruct the Hamiltonian, including both its spectrum and its eigenvectors. We demonstrated how the Hamiltonian spectrum is recovered finding all the extreme density operators which minimize the mean values of the Hamiltonian. The suggested approach for finding all the properties of the Hamiltonian has common features with quantum tomography. In the quantum tomography approach, one determines the density operator by making measures in all the reference frames of the information contained in its mean value. Therefore, the mean values are calculated for an ensemble of basis determined by sufficient number of parameters associated to unitary matrices transforming the density operator (or equivalently, changing the basis where the means are calculated). Thus, we found that the description of the Hamiltonian by the mean values approach and the description of density operators by tomographic approach have common features. We demonstrated that for and 3; the Hamiltonian spectrum is completely recovered using the mean value approach. We surmise on the basis of the number of parameters of the test variational density matrix that the complete recovering of the spectrum takes place for arbitrary dimensions. The previous conjecture is based in the following calculation: For pure states a test density matrix with 4 real parameters can be used and then one can get the exact energy spectrum of the system while with the extremal density matrices associated to the CS one has only approximate results for the energy levels. The problem of the convenience of our approach in comparison with standard diagonalisation of the Hamiltonian has to be clarified. We point out that our approach can be applied to arbitrary hermitian matrices, i.e., to arbitrary observables. The known properties of density matrices and von Neumann entropy were used to formulate and verify a new inequality for a qudit system between the energy and the entropy. We mapped the Hamiltonian matrix onto a densitylike nonnegative matrix, and applied the property of nonnegativity of relative entropy connecting two density matrices. The obtained result for the sum (difference) of entropies and mean value of the Hamiltonian is that this sum (difference) is bound. The bound is determined by the partition function evaluated at a specific value of its argument. We verified the new inequality on an example of a qubit system. The obtained result is valid for all the states, including entangled states of multiqudit systems. One can get analogous inequalities for other hermitian operators corresponding to physical observables. We will study this aspect of the inequalities in future publications.
6 Acknowledgements
This work was partially supported by CONACyTMéxico and PAPIITUNAM (IN110114).
References
References
 [1] LópezMoreno E and Castaños O, 1996 Phys. Rev. C 54 2374.
 [2] Castaños O, LópezPeña R, Hirsch J G, and LópezMoreno E, 2005 Phys. Rev. B 72 012406.
 [3] Castaños O, NahmadAchar E, LópezPeña R, and Hirsch J G, 2011 Phys. Rev. A 83 051601.

[4]
Cordero S, Castaños O, LópezPeña R, and NahmadAchar E,
2013 J. Phys. A: Math. Theor. 46 505302.  [5] Gilmore R 1979 J. Math. Phys. 20 891.
 [6] Gilmore R 1981 Catastrophe Theory for scientists and engineers Wiley (New York).
 [7] Bertrand J and Bertrand D, 1987 Found. Phys. 17 397.
 [8] Vogel K and Risken H, 1989 Phys. Rev. A 40 2847.
 [9] Mancini S, Man’ko V I, and Tombesi P, 1996 Phys. Lett A 213 1.
 [10] Ibort A, Man’ko V I, Marmo G, Simoni A, and Ventriglia F, 2009 Phys. Scr. 79 065013.

[11]
Castaños O, LópezPeña R, Man’ko M A, and Man’ko V I ,
2004 Phys. A: Math. Gen. 37 8529.  [12] Dodonov V V and Man’ko V I, 1997 Phys. Lett. A 229 335.
 [13] Man’ko V I and Man’ko O V, 1997 J. Exp. Theor. Phys. 85 430.
 [14] Castaños O, LópezPeña R, Man’ko M A, and Man’ko V I, 2003, J. Phys. A: Mat. Gen. 36 4677.
 [15] Castaños O, LópezPeña R, Man’ko M A, and Man’ko V I, 2003, J. J. Opt. B: Semiclass. Opt. 5 227.
 [16] Landau L D, 1927 Z. Physik 45 430.
 [17] von Neumann J, 1932 Mathematische Grundlagen der Quantenmechanik, Springer Berlin.
 [18] Tilma T and Sudarshan E C G, 2002 J. Phys A: Math. Gen. 3510467.
 [19] Kimura G, 2003 Phys. Lett. A 314 339.
 [20] Byrd M S and Khaneja N, 2003 Phys. Rev. A 68 062322.
 [21] Jarlskog C, 2006 J. Math. Phys. 47 013507.
 [22] Brüning E, Mäkelä H, Messina A, and Petruccione F, 2012 J. Mod. Opt. 59 1.
 [23] Akhtarshenas S J, 2006 Opt. Spectrosc. 103 411.

[24]
Nielsen M A, and Chuang I L, 2010 Quantum Computation and Quantum Information ,
Cambridge University Press.  [25] Chernega V N, Man’ko O V, and Man’ko V I, 2014 J. Russ. Laser Res. 35 278
 [26] Man’ko M A, and Man’ko V I, 2014 Phys. Scr. T160 014030.

[27]
Kubo R, 1965 Statistical Mechanics: An Advanced Course with Problems and Solutions,
Elsevier Science.  [28] Milburn G J, Corney J, Wright E W, and Walls D F, 1997 Phys. Rev. Lett. A 55 4318.

[29]
PérezCampos C, GonzálezAlonso J R, Castaños O, and LópezPeña R,
2010 Ann. Phys. (N.Y.) 325 325.  [30] Gerdt V P, Khvedlidze A M, and Palii Yu G, 2014 J. Math. Sci. 200 682.