Quantum corrections are obtained for the single-particle density matrix in a semiclassical ensemble where the distribution is unrestricted. A form for the density matrix containing explicitly a function of the wave-mechanical Hamiltonian operator is devised; a formalism is then developed to decompose this operator function into effectively classical and nonclassical parts. The classical part corresponds to the semiclassical density matrix and the quantum corrections are obtained from the nonclassical terms.The quantum-corrected density matrix for a spherically symmetrical and a one-dimensional system have been evaluated. The densities for a linear harmonic oscillator and a Coulomb potential for a Fermi-Dirac distribution have also been examined in some detail.
Quantum corrections are obtained for the single-particle density matrix in a semiclassical ensemble where the distribution is unrestricted. A form for the density matrix containing explicitly a function of the wave-mechanical Hamiltonian operator is devised; a formalism is then developed to decompose this operator function into effectively classical and nonclassical parts. The classical part corresponds to the semiclassical density matrix and the quantum corrections are obtained from the nonclassical terms.The quantum-corrected density matrix for a spherically symmetrical and a one-dimensional system have been evaluated. The densities for a linear harmonic oscillator and a Coulomb potential for a Fermi-Dirac distribution have also been examined in some detail.
The semiclassical $\ensuremath{\hbar}$ expansion of the one-particle density matrix for a two-dimensional Fermi gas is calculated within the Wigner transform method of B. Grammaticos and A. Voros [Ann. Phys. (N.Y.) 123, 359 (1979)], originally developed in the context of nuclear physics. The method of Grammaticos and Voros has the virtue of preserving both the Hermiticity and idempotency of the density matrix to all orders in the $\ensuremath{\hbar}$ expansion. As a topical application, we use our semiclassical expansion to go beyond the local-density approximation for the construction of the total dipole-dipole interaction energy functional of a two-dimensional, spin-polarized dipolar Fermi gas. We find a finite, second-order gradient correction to the Hartree-Fock energy, which takes the form $\ensuremath{\varepsilon}{(\mathbf{\ensuremath{\nabla}}\ensuremath{\rho})}^{2}/\sqrt{\ensuremath{\rho}}$, with $\ensuremath{\varepsilon}$ being small ($|\ensuremath{\varepsilon}|\ensuremath{\ll}1$) and negative. We test the quality of the corrected energy by comparing it with the exact results available for harmonic confinement. Even for small particle numbers, the gradient correction to the dipole-dipole energy provides a significant improvement over the local-density approximation.
For systems which obey quantum mechanics, the formulation of the problem of treating large numbers of particles is, of course, somewhat different than it is for classical systems. The microscopic description of the system is provided by a wave function which (in the absence of spin) is a function of the classical coordinates { q i }. The mathematical model is provided by a Hamiltonian operator H which is often obtained from the corresponding classical Hamiltonian by the replacement p i → ( ħ /i)(∂/∂ q i ). In other cases the form of the Hamiltonian operator is simply postulated. The microscopic dynamics are provided by the Schrödinger equation i ħ (∂Ψ/∂ t ) = H Ψ which requires as initial condition the knowledge of the wave function Ψ({ q i }, t ) at some initial time t 0 . (Boundary conditions on Ψ({ q i }, t ) must be specified as part of the description of the model as well.) The results of experiments in quantum mechanics are characterized by operators, usually obtained, like the Hamiltonian, from their classical forms and termed observables. Operators associated with observables must be Hermitian. In general, the various operators corresponding to observables do not commute with one another. It is possible to find sets of commuting operators whose mutual eigenstates span the Hilbert space in which the wave function is confined by the Schrödinger equation and the boundary conditions. A set of such (time independent) eigenstates, termed ψ ν ( q ), is a basis for the Hilbert space.
We state in the following a structure theory of the density matrices from the standpoint of secind quantization. The difficulties involved in the quantum-mechanical many body problem have so far proved to be so tremendous that as yet some have been obliged to stop before exact treatments and others have retired to phenomenological model theories. There are two main tools to approach the problem, the theory of second quantization and the theory of density matrix. The former, in spite of its generality and compactness, has never revealed us its adaptability to the practical problems of quantum statistical mechanics. In fact, about the exchange effect, namely the statistical correlation of particles, the latter informs us the detailed features of it by application of the group theory, while the former lacks, at present, such an exposition. For a long time a more general and simpler formulation than the current group theoretic method has been desired in the theoretical researches of spectroscopy. A distinguished contribution has been made by P. Jordan (1934) and H. Ostertag (1937) from the standpoint of second quantization. H. Ostertag made use of the theory of hypercomplex ring and introdeced the original concept of the elementary algebra. (Part II) Since the latter part of his formulation, however, seems to be complicated, it will be better, for clearness' sake, to avoid this part. Accordingly the concepts of the elementary algebra and the symmetric algebra will be introduced as a remedy for it. Then we shall complete our theory to include the whole theory of symmetric group. These concepts will serve to clarify and generalize the previous formulations and to unify the several group theoretic standpoints of the spectroscopy. In §1 and §2 the representation of the quantized density of Bose or Fermi particles will reveal itself as the density matrix of Bose or Fermi particles or as the density matrix of many electrons with a certain term of spin multiplets. Indeed, the theory of the statistical operator, or the density matrix of J. von Neumann and P.A.M. Dirac, the equilibrium properties of which have been fully examined by one of the authors (1940), is the natural formalism of dealing with quantum gases. A contribution has been made by W. Kofink (1938) to introduce the intrinsic degrees of freedom of particles. This attempt has resulted in a fair seccess. His formulation consists in making use of the initial stages of the theory of spectroscopy, which will be reexamined in the second section and Part II. In §3, as an example at zero temperature, the Hartree-Fock equation for a prescribed orbital and spin configuration will be derived, and in §4 some properties of the reduced density matrix will be illustrated.
We state in the following a structure theory of the density matrices from the standpoint of secind quantization. The difficulties involved in the quantum-mechanical many body problem have so far proved to be so tremendous that as yet some have been obliged to stop before exact treatments and others have retired to phenomenological model theories. There are two main tools to approach the problem, the theory of second quantization and the theory of density matrix. The former, in spite of its generality and compactness, has never revealed us its adaptability to the practical problems of quantum statistical mechanics. In fact, about the exchange effect, namely the statistical correlation of particles, the latter informs us the detailed features of it by application of the group theory, while the former lacks, at present, such an exposition. For a long time a more general and simpler formulation than the current group theoretic method has been desired in the theoretical researches of spectroscopy. A distinguished contribution has been made by P. Jordan (1934) and H. Ostertag (1937) from the standpoint of second quantization. H. Ostertag made use of the theory of hypercomplex ring and introdeced the original concept of the elementary algebra. (Part II) Since the latter part of his formulation, however, seems to be complicated, it will be better, for clearness' sake, to avoid this part. Accordingly the concepts of the elementary algebra and the symmetric algebra will be introduced as a remedy for it. Then we shall complete our theory to include the whole theory of symmetric group. These concepts will serve to clarify and generalize the previous formulations and to unify the several group theoretic standpoints of the spectroscopy. In §1 and §2 the representation of the quantized density of Bose or Fermi particles will reveal itself as the density matrix of Bose or Fermi particles or as the density matrix of many electrons with a certain term of spin multiplets. Indeed, the theory of the statistical operator, or the density matrix of J. von Neumann and P.A.M. Dirac, the equilibrium properties of which have been fully examined by one of the authors (1940), is the natural formalism of dealing with quantum gases. A contribution has been made by W. Kofink (1938) to introduce the intrinsic degrees of freedom of particles. This attempt has resulted in a fair seccess. His formulation consists in making use of the initial stages of the theory of spectroscopy, which will be reexamined in the second section and Part II. In §3, as an example at zero temperature, the Hartree-Fock equation for a prescribed orbital and spin configuration will be derived, and in §4 some properties of the reduced density matrix will be illustrated.
Abstract In this chapter, we continue our parallel exploration of statistical and computational physics, but now move to the field of quantum mechanics, where the density matrix generalizes the classical Boltzmann distribution. The density matrix constructs the quantum statistical probabilities from their two origins: the quantum mechanical wave functions and the Boltzmann probabilities of the energy levels. The density matrix can thus be defined in terms of the complete solution of the quantum problem (wave functions, energies), but all this information is available only in simple examples, such as the free particle or the harmonic oscillator. A simple general expression exists for the density matrix only at high temperature.
The density matrix for a many-electron system has been examined. A formalism has been arrived at which facilitates its evaluation in terms of a basis corresponding to a discrete spectrum of eigenvalues. This allows a representation to be employed associated with a suitably chosen approximation to the Hamiltonian. Thereby, reasonably accurate estimates of the properties of many-electron systems may be anticipated for low orders of approximation.In developing the formalism, attention was focused upon a means for approximating the quantity exp(zH), where H is the many-particle Hamiltonian and z is a complex number. It was found possible to represent this quantity exactly as an infinite product of exponential factors, each of which depends upon a portion of the Hamiltonian alone. The relation of this result to statistical mechanical applications is indicated. An approximation procedure is described in terms of which various finite products of exponential factors approach the desired quantity in the limit of indefinitely increasing numbers of factors; this limit is approached in a manner that permits any given approximation to contain all the terms of the previous approximation.The first order of approximation was employed to approximate the many-particle density matrix. The resulting theory was applied to the helium-like systems: He, Li+, Be++. The energies of the ground 1S state of helium and the three lowest states of Li+ and Be++ were calculated. The agreement with experimental term values is reasonably satisfactory, the discrepancy between calculations and experiments not exceeding one percent. The triplet-singlet splitting for Li+ was in the order observed; the magnitude calculated was about one-half that observed.The relation between the present theory and the Thomas-Fermi theory is discussed. It is pointed out that the former is the analog of the latter when the density matrix is expressed in other than eigenfunctions of momentum. The quasi-classical character of the approximation is discussed.