Abstract

Here we use the extended phase space formulation of quantum statistical mechanics proposed in an earlier work to define an extended lagrangian for Wigner's functions (WFs). The extended action defined by this lagrangian is a function of ordinary phase space variables. The reality condition of WFs is employed to quantize the extended action. The energy quantization is obtained as a direct consequence of the quantized action. The technique is applied to find the energy states of harmonic oscillator, particle in the box, and hydrogen atom as the illustrative examples.

1. Introduction

Among different methods proposed for the quantum statistical mechanics, one may employ the phase space distribution functions to calculate the average quantities. Wigner [1] pioneered in this direction. Based on the quantum state functions he proposed distribution functions with interesting properties and wide range of applicability. Since then many authors have examined Wigner’s distribution functions, and an amazingly large amount of literature has been amassed. Wigner himself used them to make quantum corrections for thermodynamic equilibrium. Moyal [2] and Husimi [3] developed Wigner’s technique extensively. They also showed that averaging by WFs is equivalent, in quantum mechanics, by an averaging in which noncommutating functions of operators ̂𝑝 and Μ‚π‘ž are ordered according to Weyl’s ordering rule [4]. Further improvements were made by Margenau and Hill [5], Glauber [6], Mehta [7], O’Connell and Wigner [8], Hillery et al. [9], and Agarwal [10]. Another way of studying the quantum statistical mechanics is proposed by Sobouti and Nasiri [11, 12] by generalizing the principle of least action in configuration space to one in phase space. In this method the phase space coordinates are treated as independent variables on the virtual system point trajectories in phase space. This further permits a definition of extended phase space (EPS), extended lagrangians and hamiltonians. A subsequent canonical quantization procedure or using the path integral technique [13] through these extended entities gives a formulation for quantum statistical mechanics. To demonstrate the simplicity and advantages of the EPS formulation, the technique is applied to solve a few problems. Treatment of Bloch’s theory, partition functions for simple harmonic, and linear potentials are worked out by Nasiri et al. [12]. Obtaining a unique expression for conductivity of a dissipative system of charged particles in different gauges is done using EPS technique by Khademi and Nasiri [14]. Considerable assistance of the EPS technique for studying the dissipative quantum systems is shown by Nasiri and Safari [15] and Razavy [16]. Here we use the later method to obtain an extended hamiltonian governing the evolution of WFs by an appropriate canonical transformation in EPS. An arbitrary extended lagrangian is then obtained by Legender’s transformation. The extended action obtained by this lagrangian is quantized upon the reality condition of WFs that directly yields the quantization of energy. The layout of the paper is as follows. In Section 2, a brief review of the EPS formalism is introduced. In Section 3 the extended hamiltonian and lagrangian are introduced for WFs, and the quantization of extended action and energy are illustrated by a few examples. Section 4 is devoted to conclusions.

2. Background

In this paper we use a direct approach to quantum statistical mechanics that are proposed by extending the conventional phase space and by applying the canonical quantization procedure to the extended quantities in this space (the EPS formalism).

2.1. EPS Formalism

In this formalism an extended lagrangian is introduced as follows [11]:

β„’(𝑝,π‘ž,̇𝑝,Μ‡π‘ž)=βˆ’Μ‡π‘π‘žβˆ’Μ‡π‘žπ‘+𝐿𝑝(𝑝,̇𝑝)+πΏπ‘ž(π‘ž,Μ‡π‘ž),(1) where𝐿𝑝(𝑝,̇𝑝)and πΏπ‘ž(π‘ž,Μ‡π‘ž)are lagrangians in p and q space, and 𝑝={𝑝𝑖},π‘ž={π‘žπ‘–} are the collection of generalized momenta and coordinates specifying a system, respectively. Assuming the phase space coordinates p and q to be independent variables on the virtual trajectories allows one to define momenta πœ‹π‘ and πœ‹π‘ž, conjugate to p and q, as follows:

πœ‹π‘=πœ•β„’=πœ•Μ‡π‘πœ•πΏπ‘πœ•Μ‡π‘βˆ’π‘ž,(2a)πœ‹π‘ž=πœ•β„’=πœ•Μ‡π‘žπœ•πΏπ‘žπœ•Μ‡π‘žβˆ’π‘.(2b) Then, using these extended momenta, one may define an extended hamiltonian as follows:

ℋ𝑝,π‘ž,πœ‹π‘,πœ‹π‘žξ€Έ=πœ‹π‘Μ‡π‘+πœ‹π‘žξ€·Μ‡π‘žβˆ’β„’(𝑝,π‘ž,̇𝑝,Μ‡π‘ž)=𝐻𝑝+πœ‹π‘žξ€Έξ€·,π‘žβˆ’π»π‘,π‘ž+πœ‹π‘ξ€Έ.(3) Using the quantization rule [11], one gets the following dynamical equation:

π‘–β„πœ•πœ’=ξ‚»π»ξ‚΅πœ•πœ•π‘‘=β„‹πœ’π‘βˆ’π‘–β„ξ‚Άξ‚΅πœ•πœ•π‘ž,π‘žβˆ’π»π‘,π‘žβˆ’π‘–β„πœ’=ξ“πœ•π‘ξ‚Άξ‚Όπ‘›=1(𝑖ℏ)π‘›ξ‚»πœ•π‘›!π‘›π»πœ•π‘π‘›πœ•π‘›πœ•π‘žπ‘›βˆ’πœ•π‘›π»πœ•π‘žπ‘›πœ•π‘›πœ•π‘π‘›ξ‚Όπœ’,(4) where πœ’(𝑝,π‘ž,𝑑)=πœ“(π‘ž,𝑑)πœ‘βˆ—(𝑝,𝑑)π‘’βˆ’(π‘–π‘π‘ž/ℏ) is distribution function in EPS where βˆšπœ“(π‘ž,𝑑)=(1/∫2πœ‹β„)π‘‘π‘πœ‘(𝑝,𝑑)π‘’π‘–π‘π‘ž/ℏ, πœ‹π‘=βˆ’π‘–β„(πœ•/πœ•π‘), and πœ‹π‘ž=βˆ’π‘–β„(πœ•/πœ•π‘ž). The averaging rule for an observable𝑂(𝑝,π‘ž), a c-number operator in this formalism specified by antistandard ordering rule, is

ξ€œβŸ¨π‘‚(𝑝,π‘ž)⟩=𝑂(𝑝,π‘ž)πœ’βˆ—(𝑝,π‘ž)π‘‘π‘π‘‘π‘ž.(5) For more details, one may consult Hillery et al. [9] and Sobouti and Nasiri [11].

2.2. The Extended Canonical Transformation and WFs

The quantization procedure of Section 2 assumes p and q as c-numbers, similar to Hamilton’s formulation of classical mechanics. Therefore, the rules of linear canonical transformations in classical mechanics are directly applicable to EPS quantities. These linear canonical transformations correspond to the unitary transformations in quantum mechanics [17].

A simple infinitesimal extended canonical transformation on 𝑝,π‘ž,πœ‹π‘,and πœ‹π‘ž gives WFs and the corresponding evolution equation. Let

π‘βŸΆπ‘+π›Ώπœ†πœ‹π‘ž,πœ‹π‘βŸΆπœ‹π‘,(6a)π‘žβŸΆπ‘ž+π›Ώπœ†πœ‹π‘,πœ‹π‘žβŸΆπœ‹π‘ž,(6b) where π›Ώπœ†β‰ͺ1, the generator of this infinitesimal transformation is

𝐺=πœ‹π‘πœ‹π‘ž.(7) The corresponding unitary transformation for finiteπœ† is

π‘ˆ=π‘’π‘–πœ†πΊ/ℏ=π‘’π‘–πœ†β„πœ•2/πœ•π‘πœ•π‘ž,π‘ˆπ‘ˆ+=1.(8) The above transformation, forπœ†=βˆ’1/2, transforms the distribution πœ’(𝑝,π‘ž,𝑑), to those of the Wigner [1]. That is,

=ξ€œπœ“ξ‚€1π‘Š(𝑝,π‘ž,𝑑)=π‘ˆπœ’(𝑝,π‘ž,𝑑)π‘ž+2ξ‚πœ“β„πœβˆ—ξ‚€1π‘žβˆ’2ξ‚π‘’β„πœπ‘–π‘πœπ‘‘πœ,(9) and the evolution equation (4) changes into the Wigner equation:

πœ•π‘Šπœ•π‘‘=βˆ’π‘πœ•π‘Š+ξ“πœ•π‘žπ‘›=01ℏ(2𝑛+1)!2𝑖2π‘›πœ•2𝑛+1π‘‰πœ•π‘ž2𝑛+1πœ•2𝑛+1π‘Šπœ•π‘2𝑛+1.(10) A similar technique could be used to obtain other distribution functions and their corresponding ordering rules [18–20].

3. The Extended Hamiltonian and Extended Lagrangian for WFs

By transformation (8) the extended hamiltonian of (3) becomes

ℋ𝑀𝑝,π‘ž,πœ‹π‘,πœ‹π‘žξ€Έξ‚€1=𝐻𝑝+2πœ‹π‘ž1,π‘žβˆ’2πœ‹π‘ξ‚ξ‚€1βˆ’π»π‘βˆ’2πœ‹π‘ž1,π‘ž+2πœ‹π‘ξ‚.(11) The subscript 𝑀denotes the quantities for WFs. For 𝐻=𝑝2/2π‘š+𝑉(π‘ž) (11) gives

ℋ𝑀=1ξ‚€12π‘šπ‘+2πœ‹π‘žξ‚2ξ‚€1+π‘‰π‘žβˆ’2πœ‹π‘ξ‚βˆ’1ξ‚€12π‘šπ‘βˆ’2πœ‹π‘žξ‚2ξ‚€1βˆ’π‘‰π‘ž+2πœ‹π‘ξ‚.(12) The extended lagrangian, ℒ𝑀, for WFs may be defined by the following Legendre transformation:

ℋ𝑀𝑝,π‘ž,πœ•β„’π‘€,πœ•Μ‡π‘πœ•β„’π‘€ξ‚Άπœ•Μ‡π‘ž=Μ‡π‘πœ•β„’π‘€πœ•Μ‡π‘+Μ‡π‘žπœ•β„’π‘€πœ•Μ‡π‘žβˆ’β„’π‘€(𝑝,π‘ž,̇𝑝,Μ‡π‘ž).(13) Equation (13) serves as a differential equation to determine ℒ𝑀. Note that the extended lagrangians are not unique and they just serve to obtain the correct equations of motion through the Euler-Lagrange equations for configuration and momentum spaces.

The solutions of the SchrΓΆdinger equation can be written as [21]

πœ“(π‘ž,𝑑)=π‘…π‘ž(π‘ž,𝑑)π‘’π‘–π‘†π‘ž(π‘ž)/ℏ,(14) where π‘†π‘ž(π‘ž) is the ordinary action in q space defined by

π‘†π‘ž(ξ€œπΏπ‘ž)=π‘ž(π‘ž,Μ‡π‘ž,𝑑)𝑑𝑑.(15) By the same manner one may do in p space:

πœ‘(𝑝,𝑑)=𝑅𝑝(𝑝,𝑑)π‘’βˆ’(𝑖𝑆𝑝(𝑝)/ℏ),(16) where

𝑆𝑝(ξ€œπΏπ‘)=𝑝(𝑝,̇𝑝,𝑑)𝑑𝑑(17) is the action in p space. Substituting for πœ“ and πœ‘, from (14) and (16) in the expression given for πœ’, one gets

πœ’=πœ“πœ‘βˆ—π‘’βˆ’(π‘–π‘π‘ž/ℏ)=π‘…π‘žπ‘…π‘π‘’π‘–(π‘†π‘ž+π‘†π‘βˆ’π‘π‘ž)/ℏ.(18) Now, let the extended action be defined as

π‘†πœ’(ξ€œπ‘,π‘ž)=β„’(𝑝,π‘ž,̇𝑝,Μ‡π‘ž,𝑑)𝑑𝑑.(19) Using (1), (15), and (17) one gets

π‘†πœ’(𝑝,π‘ž)=π‘†π‘ž+π‘†π‘βˆ’π‘π‘ž,(20) which is the same as appeared in the phase of (18). Hence,

πœ’=π‘…πœ’π‘’π‘–π‘†πœ’(𝑝,π‘ž)/ℏ,(21) whereπ‘…πœ’(𝑝,π‘ž,𝑑)=𝑅𝑝(𝑝,𝑑)π‘…π‘ž(π‘ž,𝑑). Therefore, in general, any phase space distribution function,𝜌(𝑝,π‘ž,𝑑) can be written in the form of an amplitude multiplied by a phase factor [13] as follows:

𝜌(𝑝,π‘ž,𝑑)=π‘…πœŒ(𝑝,π‘ž,𝑑)π‘’π‘–π‘†πœŒ(𝑝,π‘ž)/ℏ,(22) whereπ‘…πœŒ(𝑝,π‘ž,𝑑)is the amplitude, and π‘†πœŒ(𝑝,π‘ž)is the corresponding extended action of 𝜌(𝑝,π‘ž,𝑑). Both of π‘…πœŒ and π‘†πœŒ are real functions. In the same manner for WFs, as a possible distribution functions in the phase space, we formally have

π‘Š(𝑝,π‘ž,𝑑)=𝑅𝑀(𝑝,π‘ž,𝑑)𝑒𝑖𝑆𝑀(𝑝,π‘ž)/ℏ,(23) where

𝑆𝑀(ξ€œβ„’π‘,π‘ž)=𝑀(𝑝,π‘ž,̇𝑝,Μ‡π‘ž,𝑑)𝑑𝑑,(24) and ℒ𝑀is given by (13). However, it is well known that the WFs are real functions of the phase space variables. Thus, the phase 𝑆𝑀(𝑝,π‘ž)in (23) must satisfy the following constraint:

𝑆sin𝑀(𝑝,π‘ž)ℏ=0,(25) or

𝑆𝑀(𝑝,π‘ž)=π‘›πœ‹β„,𝑛=0,1,2,….(26) Equation (26) shows the quantization of extended action. As will be shown later, due to the uncertainty principle, the values of 𝑛must be correctly chosen for a given problem. Note that the extended action in EPS is a function of phase space variables and is basically different from ordinary action π‘†π‘ž(π‘ž)defined in configuration space. As an interesting result, one has cos{𝑆𝑀(𝑝,π‘ž)/ℏ}=cos(π‘›πœ‹)=(βˆ’1)𝑛. Thus, the well known negativity of WFs, 𝑀(𝑝,π‘ž,𝑑)=(βˆ’1)𝑛𝑅𝑀(𝑝,π‘ž,𝑑), seems to have a quantum origin [22].

3.1. Energy Quantization

In this section we show that the quantization of the extended action of (26) is directly related to the quantization of energy. To see this we consider a few examples.

3.1.1. Harmonic Oscillator

For harmonic potential, 𝐻=𝑝2/2π‘š+(1/2)π‘šπœ”2π‘ž2 (12) gives

ℋ𝑀𝑝,π‘ž,πœ‹π‘,πœ‹π‘žξ€Έ=π‘π‘šπœ‹π‘žβˆ’π‘šπœ”2π‘žπœ‹π‘.(27) By (13) one gets

π‘π‘šπœ•β„’π‘€πœ•Μ‡π‘žβˆ’π‘šπœ”2π‘žπœ•β„’π‘€πœ•Μ‡π‘=Μ‡π‘πœ•β„’π‘€πœ•Μ‡π‘+Μ‡π‘žπœ•β„’π‘€πœ•Μ‡π‘žβˆ’β„’π‘€(𝑝,π‘ž,̇𝑝,Μ‡π‘ž).(28) The extended lagrangian as the solution of differential equation (28) may be obtained as

β„’π‘€πœ‹(𝑝,π‘ž,̇𝑝,Μ‡π‘ž)=πœ”ξ‚€π‘π‘šΜ‡π‘+π‘šπœ”2ξ‚π‘žΜ‡π‘ž.(29) It can be shown that ℒ𝑀 satisfies the Euler-Lagrange equations in both configuration and momentum space. Using (24) and (29) one obtains

π‘†π‘€πœ‹(𝑝,π‘ž)=πœ”ξ‚΅π‘2+12π‘š2π‘šπœ”2π‘ž2ξ‚Ά=πœ‹πœ”ξ‚€1Hβˆ’2,β„πœ”(30) where H is the quantum mechanical Hamiltonian assuming the uncertainty principle, implicitly [23]. Equations (26) and (30) give

𝐸𝑛1(𝑝,π‘ž)=𝑛+2ξ‚β„πœ”,𝑛=0,1,2,3,…,(31) as the quantization of energy of a harmonic oscillator.

3.1.2. Particle in the Box

For a particle in a one dimensional box of side L located at the right side of origin, one has

⎧βŽͺ⎨βŽͺβŽ©π‘‰(π‘ž)=∞,π‘ž<0,0,0<π‘ž<𝐿,∞,π‘ž>𝐿.(32) For 𝐻=𝑝2/2π‘š one gets the following from (12):

ℋ𝑀𝑝,π‘ž,πœ‹π‘,πœ‹π‘žξ€Έ=π‘π‘šπœ‹π‘ž.(33) By legendre’s transformation on the above extended hamiltonian, one gets for an arbitrary extended lagrangian;

ℒ𝑀(𝑝,π‘ž,̇𝑝,Μ‡π‘ž)=𝐿̇𝑝.(34) Using (24) and (34) the extended action becomes

𝑆𝑀(𝑝,π‘ž)=𝐿𝑝=𝐿(2π‘šπ»)1/2.(35) Then, by (26) and squaring (35) one gets

𝐸𝑛=𝑛2πœ‹2ℏ22π‘šπΏ2,𝑛=1,2,…,(36) where 𝑛 cannot assume zero due to the uncertainty principle. Equation (36) shows the energy levels of a particle in the box.

3.1.3. Hydrogen Atom

For the electron in the hydrogen-like atom, as a bound system with Z protons, the hamiltonian is 𝑝2π‘Ÿ/2π‘š+𝑙2/2π‘šπ‘Ÿ2βˆ’(𝑧𝑒2/π‘Ÿ) [24]. By (12) one gets

β„‹π‘€ξ€·π‘π‘Ÿ,π‘Ÿ,πœ‹π‘π‘Ÿ,πœ‹π‘Ÿξ€Έ=π‘π‘Ÿπ‘šπœ‹π‘Ÿ+2β„“2π‘Ÿπœ‹π‘π‘Ÿ+2π‘šπ‘π‘’2π‘Ÿ2πœ‹π‘π‘Ÿ+ξ€·π‘šπ‘§π‘’2ξ€Έπœ‹/23π‘π‘Ÿ2π‘š(π‘Ÿβˆ’(1/2)πœ‹π‘π‘Ÿ)2(π‘Ÿ+(1/2)πœ‹π‘π‘Ÿ)2.(37) A suitable extended lagrangian can be obtained by (25) to be

β„’π‘€ξ€·π‘π‘Ÿ,π‘Ÿ,Μ‡π‘π‘Ÿξ€Έ=,Μ‡π‘Ÿπ‘šπ‘2√2πœ‹π‘’2ξ‚΅π‘π‘Ÿπ‘šΜ‡π‘π‘Ÿβˆ’β„“2Μ‡π‘Ÿπ‘šπ‘Ÿ3+𝑒2Μ‡π‘Ÿπ‘Ÿ2×𝑝2π‘Ÿ+β„“2π‘š22π‘šπ‘Ÿ2βˆ’π‘§π‘’2π‘Ÿξ‚Άβˆ’3/2.(38) The extended action becomes

π‘†π‘€ξ€·π‘π‘Ÿξ€Έ=ξ‚™,π‘Ÿπ‘š2πœ‹π‘π‘’2π»βˆ’1/2.(39) Then, by (37) and squaring (39) one gets

||𝐸𝑛||=π‘šπ‘2𝑒42ℏ2𝑛2=𝑍2𝑒22π‘Ž0𝑛2,𝑛=1,2,…,(40) where π‘Ž0 is the Bohr radius, and (40) gives the energy levels of the hydrogen atom. The same as the previous case 𝑛 starts from 1 due to the uncertainty principle.

4. Conclusion

The EPS formulation of quantum statistical mechanics permits the definition of extended hamiltonians and extended lagrangians. A canonical quantization in EPS gives a dynamical equation for the phase space distribution functions. Wigner’s distribution functions and their corresponding evolution equation can be obtained by a simple extended canonical transformation in EPS. The extended hamiltonian through Legendre’s transformation provides a differential equation for the corresponding extended lagrangian. Having the extended lagrangian one may define the extended action as a function of ordinary phase space variables. Applying the reality condition of the WFs yields the quantization of the extended action which in turn gives the energy quantization. The energy levels of the harmonic oscillator, the particle in the box and the hydrogen atom are worked out as a few illustrative examples.