Quantum Mechanics, Volume 3. Claude Cohen-Tannoudji

Quantum Mechanics, Volume 3

this family the operator that yields the lowest value for Φ. This operator will then be the optimal operator within this family. Furthermore, this operator yields an upper value for the grand potential, with an error of second order with respect to the error made on ρ.

2. Approximation for the equilibrium density operator

We now use this variational principle with a family of density operators that leads to manageable calculations.

2-a. Trial density operators

The Hartree-Fock method is based on the assumption that a good approximation is to consider that each particle is independent of the others, but moving in the mean potential they create. We therefore compute an approximate value of the density operator by replacing the Hamiltonian Ĥ by a sum of independent particles’ Hamiltonians :

(25)

We now introduce the basis of the creation and annihilation operators, associated with the eigenvectors of the one-particle operator :

(26)

The symmetric one-particle operator can then be written, according to relation (B-14) of Chapter XV:

(27)

where the real constants are the eigenvalues of the operator .

We choose as trial operators acting in the Fock space the set of operators that can be written in the form corresponding to an equilibrium in the grand canonical ensemble – see relation (42) of Appendix VI. We then set:

(28)

where is any symmetric one-particle operator, the constant β the inverse of the temperature defined in (4), μ a real constant playing the role of a chemical potential, and the trace of :

(29)

Consequently, the relevant variables in our problem are the states , which form an arbitrary orthonormal basis in the individual state space, and the energies . These variables determine the ak as well as , and we have to find which of their values minimizes the function:

(30)

Taking (27) and (28) into account, we can write:

(31)

The following computations are simplified since the Fock space can be considered to be the tensor product of independent spaces associated with the individual states ; consequently, the trial density operator (28) can be written as a tensor product of operators each acting on a single mode k:

(32)

2-b. Partition function, distributions

Equality (32) has the same form as relation (5) of Complement B_XV, with a simple change: the replacement of the free particle energies ek = ħ²k²/2m by the energies , which are as yet unknown. As this change does not impact the mathematical structure of the density operator, we can directly use the results of Complement B_XV.

α. Variational partition function

The function only depends on the variational energies , since the trace of (32) may be computed in the basis {}, which yields:

(33)

We simply get an expression similar to relation (7) of Complement B_XV, obtained for an ideal gas. Since for fermions nk can only take the values 0 and 1, we get:

(34)

whereas for bosons nk varies from 0 to infinity, so that:

(35)

In both

Скачать книгу