InTechOpen uses cookies to offer you the best online experience. By continuing to use our site, you agree to our Privacy Policy.

Engineering » Chemical Engineering » "Thermodynamics - Interaction Studies - Solids, Liquids and Gases", book edited by Juan Carlos Moreno-Pirajan, ISBN 978-953-307-563-1, Published: November 2, 2011 under CC BY 3.0 license. © The Author(s).

Chapter 17

Statistical Mechanics That Takes into Account Angular Momentum Conservation Law - Theory and Application

By Illia Dubrovskyi
DOI: 10.5772/19645

Article top

Statistical Mechanics That Takes into Account Angular Momentum Conservation Law - Theory and Application

Illia Dubrovskyi1

1. Introduction

The fundamental problem of statistical mechanics is obtaining an ensemble average of physical quantities that are described by phase functions (classical physics) or operators (quantum physics). In classical statistical mechanics the ensemble density of distribution is defined in the phase space of the system. In quantum statistical mechanics the space of functions that describe microscopic states of the system play a role similar to the classical phase space. The probability density of the system detection in the phase space must be normalized. It depends on external parameters that determine the macroscopic state of the system.

An in-depth study of the statistical mechanics foundations was presented in the works of A.Y. Khinchin (Khinchin, 1949, 1960). For classical statistical mechanics an invariant set was introduced. It would be mapped into itself by transforming with the Hamilton equations. The phase point of the isolated system remains during the process of the motion at the invariant set at all times. If the system is in the stationary equilibrium state, this invariant set has a finite measure. The Ergodic hypothesis asserts that in this case the probability dP (R) to detect this system at any point R of the phase space is:


where ΩΣ - the measure (phase volume) of the invariant set ∑; φΣ(R)- the characteristic function of the invariant set, which is equal to one if the point R belongs to this set, and is equal to zero in all other points of the phase space; dΓ=i=1Ndpidri- the phase space volume element. The number of distinguishable states in a phase space volume element dΓ is[(2π)3NN!]1. The system that will be under consideration is a collection of Nstructureless particles. The averaged value of a phase function F(R) isF¯=ΓF(R)dP(R). Here the integral goes over all phase spaceΓ. This is microcanonical distribution. A characteristic function often would be presented asδ[f(R)z], where f(R) is a phase function and z is it’s fixed value.

A hypersurface in a hyperspace is a set with zero measure. Therefore the invariant set is determined as a thin layer that nearly envelops the hypersurface in the phase space. The determining equations of this hypersurface are the equalities that fix the values of controllable motion integrals. A controllable motion integral is a phase function, the value of which does not vary with the motion of the system and can be measured. An isolated system universally has the Hamiltonian that does not depend on the time explicitly, and is the controllable motion integral. A fixed value of the Hamiltonian is the energy of the system. The kinetic energy of majority of systems is a positive definite quadric form of all momenta. It determines a closed hypersurface in the subspace of momenta of the phase space. If motions of all particles are finite, the hypersurface of the fixed energy is closed and the layer that envelops it has the finite measure. Then this hypersurface can determine the invariant set of the system. A finiteness of motions of particles as a rule is provided by enclosing the system in an envelope that reflects particles without changing their energy, if the system is considered as isolated. It is common in statistical mechanics to consider the layer enveloping the energy hypersurface as the invariant set. But A.Y. Khinchin (Khinchin, 1949) shows that other controllable integrals of the system, if they exist, must be taken into account. In the general case an isolated system can have another two vector controllable integrals. That is the total momentum of the system, and the total angular momentum relative to the system’s mass centre. The total momentum is a sum of all momenta of particles. If the volume of the system is bounded by an external field or an envelope, the total momentum does not conserve. In the absence of external fields the total momentum conservation cannot make particle motions finite. Therefore the total momentum cannot be a controllable motion integral that determines the invariant set.

The angular momentum is another case. A vector of angular momentum relative to the mass center always is conserved in an isolated system. If this vector is nonzero, a condition should exist that provides a limitation of a gas expansion area. For example, nebulas do not collapse because they rotate, and do not scatter because of the gravitation. In the system of charged particles in a uniform magnetic field the conservation of the angular momentum provides a limitation of a gas expansion area (confinement of plasma). If a gas system is enclosed into envelope, and total system has nonzero angular momentum, the vector of the angular momentum should be conserved. However an envelope can have the non-ideal form and surface. That is the cause of the failure to consider the angular momentum of the gas as a controllable motion integral (Fowler, & Guggenheim, 1939). But if the cylindrical envelope rotates and the gas rotates with the same angular velocity deviations of the angular momentum of the gas from the fixed value as the result of reflections of particles from the envelope should be small and symmetric with respect to a sign. These fluctuations are akin to energy fluctuations for a system that is in equilibrium with a thermostat. Therefore the angular momentum conservation in specific cases can determine the invariant set and the thermodynamical natures of the system together with the energy conservation. Taking into account all controllable motion integrals is the necessary condition of the validity of the Ergodic hypothesis (Khinchin, 1949).

There is a contradiction in physics at the present time. Firstly, it has been proven that in the equilibrium state a system spin can exist only if the system is rigid and can rotate as a whole (Landau, & Lifshitz, E.M., 1980a). Therefore a gas, which supposed not be able to rotate as a whole, cannot have any angular momentum and spin. Based on this reasoning R.P. Feynman proves that an electron gas cannot have diamagnetism (the Bohr – van Leeuwen theorem) (Feynman, Leighton, & Sands, 1964). On the other hand, it is well known that density of a gas in a rotating centrifuge is non-uniform. This effect is used for the separation of isotopes (Cohen, 1951). The experiment by R. Tolman, described in the book (Pohl, 1960), is a proof of the existence of the electron gas angular momentum. In this experiment a coil was rotated and then sharply stopped. An electrical potential was observed that generated a moment of force, which decreased to zero the angular momentum of electron gas.

The contradiction described above requires creation of statistical mechanics for non-rigid systems taking into account the nonzero angular momentum conservation. This statistical mechanics differs from common one in many respects. If the angular momentum relative to the axis that passes through the mass centre conserves, the system is spatially inhomogeneous. This means that passage to the thermodynamical limit makes no sense, a spatial part of the system is not a subsystem that similar to the total system, specific quantities such as densities or susceptibilities have no physical meaning.

The microcanonical distribution is seldom used directly when the computations and the justifications of thermodynamics are done. The more usable Gibbs distribution can be deduced from microcanonical one (Krutkov, 1933; Zubarev, 1974). The Gibbs assembly describes a system that is in equilibrium with environment. These systems do not have motion integrals because they are non-isolated. All elements of the Gibbs assembly must have equal values of parameters that are determined by the equilibrium conditions. In usual thermodynamics this parameters are the temperature and the chemical potential. The physical interpretation of these parameters is getting by statistical mechanics. A rotating system can be in equilibrium only with rotating environment. The equilibrium condition in this case is apparent. That is equality of the both angular velocities of the system and of the environment. The Gibbs assembly density of distribution and thermodynamical functions in the case of a rotating classical system will be obtained in the second section of this work. It was done (Landau, & Lifshitz E.M., 1980a) but an object, to which this distribution is applied, is incomprehensible, because an angular velocity of an equilibrium gas has not been determined.

In quantum statistical mechanics the invariant set is the linear manifold of the microscopic states of the system in which the commutative operators that correspond to the controllable motion integrals have fixed eigenvalues. The phase volume of system in this case is the dimension of the manifold, if this dimension is limited. It directly determines the number of distinguishable microstates of the system that are accessible and equiprobable. The role of the angular momentum conservation in quantum statistical mechanics is similar to one in classical statistical mechanics. The method of computing this phase volume will be also proposed in the second section of this work. The Gibbs assembly density of distribution and thermodynamical functions in the case of a rotating quantum system also will be obtained.

In the third section of this work statistical mechanics of an electron gas in a magnetic field is considered. This question was investigated by many during the last century. Many hundreds experimental and theoretical works were summarized in the treatises (Lifshits, I.M. et al., 1973; Shoenberg, 1984). However, together with successful theoretical explanations of many experimental effects some paradoxes and discrepancies with observed facts remain unaccounted.

“Finally, it is shown that the presence of free electrons, contrary to the generally adopted opinion, will not give rise to any magnetic properties of the metals”. This sentence ends a short report on the presentation “Electron Theory of Metals” by N. Bohr, given at the meeting of the Philosophical Society at Cambridge. It was well-known that a charged particle in a uniform magnetic field moves in a circular orbit with fixed centre in such a way that the time average value of the magnetic moment, generated by this motion, is directed opposite to the magnetic field and equal to the derivative of the kinetic energy with respect to the magnetic field. N. Bohr computed the magnetic moment of an electron gas by statistical mechanics with the density of distribution that is determined only by a Hamiltonian. Zero result of this theory (Bohr – van Leeuwen theorem) is the first paradox. Many attempts of derivation and explanation of this were summarized in the treatise (van Vleck, 1965). The most widespread explanation was that the magnetization generated by the electrons moving far from the bound is cancelled by the near-boundary electrons that reflect from the bound. But this explanation is not correct because, when formulae are derived in statistical mechanics, any peculiarities of the near-boundary states shall not be taken into account. Another paradox of the common theory went unnoticed. It is well known that a uniform magnetic field restricts an expanse of a charged particles gas in the plane perpendicular to the field. But from common statistical mechanics it follows that the gas uniformly fills all of the bounded area. The diamagnetism of some metals also was left non-explained.

L.D. Landau (Landau, 1930) explained the diamagnetism of metals as a quantum effect. He solved the quantum problem of an electron in a uniform magnetic field. The cross-section of the envelope perpendicular to the magnetic field is a rectangle with the sides 2Lx and2Ly. The solutions are determined by three motion integrals. The first is energy that takes the valuesεnp=ωс(n+1/2)+p2/2m, where ωс is the cyclotron frequencyωс=eH/m, (e)is the charge and m is the mass of an electron, His the magnetic induction, nis a positive integer or zero. The second is the zcomponent of the momentump. The third motion integral is the Cartesian coordinate of the centre of the classical orbit. It takes the valuesyj=(π/eHLx)j, wherej=0,±1,±2,...±eHLxLy/π. The thermodynamical potential with this energy spectrum, when the spin degeneracy is taken into account, is:


Here kB is the Boltzmann constant, Tis the temperature, S=4LxLy, μis the chemical potential. If in this formula the summation over nis changed to the integration, the result Ω0does not depend onH, and the magnetic momentM=(Ω0/H)=0. That agrees to the classical and paradoxical Bohr – van Leeuwen theorem. L.D. Landau uses the Euler – Maclaurin summation formula in the first order and obtains the amendment that depends on the magnetic field. In the limit T0 the thermodynamical potential has appearance (Abrikosov, 1972):


wherepF=2mζ=(3π2)1/3(N/V)1/3, and ζ is the Fermi energy, V=4LxLyLzis the volume. This result cannot be correct because the magnetic moment does not depend on the Plank constant and thus it cannot be a quantum effect. This problem is simpler for a two-dimensional gas. In this case the common formula of the thermodynamical potential has the form:

. WhenT0, this sum can be computed without to change the summation to the integration.

It is suggested that the Fermi level is filled. In this case the magnetic moment does not depend not only on the Plank constant, but also on the number of the electrons. Therefore the fundamental formula for the thermodynamical potential is incorrect.

In the third section of this work the diamagnetism of an electron gas is investigated with taking into account the conservation of zero value of the total angular momentum in classical and quantum statistical mechanics. The paradoxes described above are eliminated; however many other theories should be reconsidered.

2. Statistical mechanics of rotating gas

For the computation of average values of macroscopic quantities it is necessary to derive a formula of the phase volume as a function of macroscopic parameters. This function is called “structural function” by Khinchin (Khinchin, 1949) and “number of accessible states (or complexions)” by Fowler (Fowler, & Guggenheim, 1939). It determines the normalizing factor in the probability density of the microcanonical distribution (1). In usual theory this function is essential to the derivation of formulae that connect statistical physics with thermodynamics.

The system that will be considered is a collection of Nstructureless particles. If forces of interaction between particles manifest themselves only at distances considerably smaller than the average distance between particles, the interaction energy of particles is essential only in a small fraction of the phase volume. Therefore, the interaction of particles can be neglected or be taken into account as a perturbation in calculating the phase volume and the average values (Uhlenbeck & Ford, 1963). Otherwise, if particles interact by a long-range force, this interaction needs to be considered using a mean field method. This is a model of an ideal gas under an external field. Meanwhile, this external field can be also a periodical crystal field. In the commonly considered cases the Hamiltonian and other phase functions of the system can be presented as the sum of identical terms, each of which depends on the coordinates and momenta of a single particle. Such phase function is said to be a summatory function.

For integration characteristic functions over a phase space the method by Krutkov will be used. The main idea of this method is to make the Laplace transformation of the δfunction with respect the value of the summatory function. Then the product of the N like exponents from the terms of the summatory function would be integrated over variables of the phase space. The inverse transformation would be made by the saddle-point method with using the large parameterN.

Let us write several equalities with a characteristic function. If the system can be divided into two independent subsystems described by non-overlapping groups of phase variables, so thatΣ=Σ1Σ2,R=R1+R2, and the determining functions possess the values independently, then


Here the multiplier (N!)1 is not taken into account because it cannot be introduced logically in classical statistical mechanics. Considering the fact that a density of distribution for a system is equal to the product of densities of distribution for subsystems, a conclusion is drawn in the treatise (Landau, & Lifshitz, E.M., 1980a) that the logarithm of the density of distribution should be an additive motion integral and, hence, it should be a linear combination of the additive controllable motion integrals, such as the energy, the momentum and the angular momentum. However, as it follows from the formula (6) this is incorrect for the microcanonical distribution, since the logarithm of the characteristic function is meaningless. A system in a thermostat does not have any motion integrals. If the invariant set is determined by some conservation laws, its characteristic function is


where φi(R) is the characteristic function that is determined by the conservation law numberi. Let us denote a set, at which the phase function A(R)is equal toa, byΣAa, its characteristic function byφAa(R), and its measure byΩAa. It is supposed that this measure is limited and is not equal to zero. Then


Here Ais the range of values of the functionA(R). The prevalent formula φAa=Δaδ(Aa) (Landau, & Lifshitz, E.M., 1980a; Uhlenbeck, & Ford, 1963) satisfies to the equalities (7) and (8) but does not satisfy to the equality (6). It is more frequently considered the separation on subsystems that conserves the total value of the function. IfA(R)=A1(R1)+A2(R2)=a, thenφAa=AφA1axφA2xdx, and the prevalent formula is correct.

2.1. Classical statistical thermodynamics of rotating gas

The formula for average values, when the conservation of the angular momentum is taken into account, has the form:


In these formulae the axis Z is parallel to the angular momentumL¯. The angular momentum of the gas isL(p,r)=i=1Nl(pi,ri). The effective Hamiltonian of the gas with the fixed angular momentum in the cylindrical coordinates can be obtained from the usual formula by substitution l1=L¯i=2Nliwith reduction of the quadratic form to the standard appearance. It is:


where U(r2+z2)is potential energy that confines the particles in bounded volume, and the second term of the Hamiltonian is the potential energy of the centrifugal force that leads to a collapse of rotating nebula into disk. This Hamiltonian is not a summatory function. Therefore the Krutkov’s method cannot be used for the subsequent computations.

Let us consider equilibrium of a gas with a rotating rigid body. The rigid body can be determined as the body in which the rotatory degree of freedom can not transfer energy and angular momentum to the internal degrees of freedom. This possibility arises when this body is a cylindrical rotating envelope with non-ideal surface filled by a gas. The state of the gas is characterized by two parameters: the temperature Tand the angular velocityω. For introducing the statistical parameter that corresponds to the thermodynamical temperature it is necessary to deduce the canonical Gibbs distribution for a system that is in equilibrium with a thermostat. That can be done, for example, by a method developed by Krutkov (Krutkov, 1933; Zubarev,1974). The conditions of the equilibrium between the rotating envelope and the gas are apparent. Those are the equalities of the temperature and the angular velocity.Let us determine the angular velocity of a gas. An angular velocity of a particle is a stochastic quantity with an average valueω¯. The sum i=1Nωi is the Gaussian random variable with the average valueNω¯=ωg. That is the angular velocity of a gas. The conditions of the equilibrium between the rotating envelope and the gas are the equality of the temperatures and


The total system can be considered as motionless if it will be described in the rotating reference frame, when the right part of the equality (11) is zero. The hollow cylinder is the envelope, the thermostat, and it keeps the gas spin. It should be named “termospinstat”. The potential energy of the centrifugal force Ucf=(mω2r2/2) is added in the Hamiltonian of the gas particle in the rotating reference frame (Landau, & Lifshitz, E.M., 1980a). The average angular momentum of the gas mωi=1Nri2¯=ωNmri2¯=L¯ depends on the angular velocity nonlinearly because the gas moment of inertia I¯(ω)=Nmri2¯ is the function of the angular velocity. This function can be obtained from the Gibbs distribution for a gas in the system of reference that rotates with the angular velocityω:


where hz and R are the dimensions of the envelope, and U0=NmR2ω2/2 is an appending constant that does the energy positive. Going to the thermodynamics (it rather can be entitled by “thermospindynamics”) it is naturally to consider σ=ω2/2as an external parameter and the moment of inertia I=mi=1Nri2as a characteristic of the gas. Then


The formulae of the isotopes separation (Cohen, 1951) can be obtained from the distribution (12). If m(ωR)2<<kBT the formula (12) can be presented as:


where F0 - is the free energy of the ideal gas that does not rotate. Hence it follows that the parametershz, R2, and correspondinglyPz, PSshould be introduced instead of the volume V and the pressureP. Other thermodynamical equations are changed also. The parameter of expansion in the formula (14) can be of the order of unity whenm=1025kg,R=1m,T=100K,ω=100s-1.

2.2. Quantum statistical thermodynamics of rotating gas

The characteristic function of the invariant set that takes into account conservation of the angular momentum in quantum statistical mechanics can be presented as a set of diagonal elements of the operator:


in the space of microstates of the system. Here H^=i=1Nh^i is Hamiltonian of gas; L^=i=1Nl^iis the operator of the total angular momentum of gas;E, L¯are values of these quantities for the considered macroscopic state;ϑ, ηare real numbers which will be defined below. As usually, let us assume that energies of one-particle states, and, hence, both eigenvalues of the operator H^ and gas energyE, are expressed by the dimensionless positive integers. This formula would be generalized by the transition to representation of secondary quantization. In this representation function of one-particle states are eigenfunctions of the one-particle Hamiltonian h^i and angular momentum. These functions are numbered using index Ψwhich consists of a pair of quantum numbers(t,l), where ε(t) is energy and l(Ψ) is the angular momentum at the stateΨ. Let us suppose that only two quantum numbers determine the state. Both the production and annihilation operators are determined by the kind of statistics. The operator (15) should be replenished by multiplier that describes the conservation of the particle number. Then the measure of the invariant set for N-particle system is:


Here the following variables are entered:


Thus, integrals are rearranged into integrals along contours which enclose the origin of coordinates. If L¯0 then the axis Z is parallel to the angular momentum L¯ andL¯>0,η>0. The lower operator or sign should be taken for the Bose statistics, and upper ones should be taken for the Fermi statistics. The expression of the measure of the invariant set in the case of quantum statistics is obtained from the apparent formula of the characteristic function (15). This expression is similar to the initial one in the Darwin – Fowler method (Fowler, & Guggenheim, 1939), where it was substantiated as a mathematical device. These contour integrals can be computed by the saddle-point method, whenL¯0. The saddle-points determine the values ofχ,ϑ,η. It can be shown that taking into account the conservation of the angular momentum does not change statistic mechanics of this model of an ideal gas whenL¯=0.

Let us describe a quantum gas in a termospinstat with the temperature T and the angular velocityω. The potential of the centrifugal force would be introduced in the Hamiltonian only in classical statistical mechanics (Landau, & Lifshitz, E.M., 1980a). Other method should be used for quantum theory. Wave functions in a system of reference that rotates with the angular velocity ω depend on the time, and should be determined from the Schrödinger equation:


Dependence of the wave function on the time should beψ=exp{i[zpz/+l(φωt)]}. Thenψ=exp{i[pz/+l(φωt)]}J|l|(r2m(εωl)pz2/)and J|l|(x)is the Bessel function. The boundary condition ψ(R)=0 determines the energy spectrum:


where j(|l|,ν) is the null of the Bessel function J|l|(x)with number in the order of increasingν. Whenν>1, j(|l|,ν)π(ν+|l|/2)andj(|l|,1)|l|. This spectrum is quasicontinuous because a distance between nearest-neighbor levels is proportional toR2. The lowest level has value (mω2R2/2)whenl=(mR2ω/). Then the reference point of energy should be altered by this value. It conforms to the appearance of the centrifugal force potential in the classical system. Energies of states with l>0 are lower than the ones of states with equal |l| andl<0. Then in the gas the part of particles with l>0should be more than half, and as result a circular current of the probability density should exist. This describes rotating of the system. With increasing argument modulus of extremes of the Bessel functions decrease. If the values of energy εand positive angular momentum are fixed the value of the null number ν in the rotating system should be lower than this value in the motionless system. Therefore, the gas density increases with distance from the axis in rotating system.

Let us compute the thermodynamical potential Ω=±kBTpνlln{1exp[(μεpνl)/kBT]}of ideal rotating gas whenexp(μ/kBT)<<1. Nulls of Bessel functions would be approximated by formulaj(|l|,ν)π(ν+|l|/2), but when ν=1 thenj(|l|,1)|l|. The computation is performed by passing from summation over ν and l to integration over j andl. This approximation for non-rotating gas leads to the result that differs from the common result by the multiplier(4/π). The result of computation for the rotating gas is:


If this result is compared with that of Eq. (14), it can be shown that amendments differ only by coefficients.

3. Statistical mechanics of electron gas in magnetic field

The review of the current status of this theory is in the paper (Vagner et al., 2006). There are some inaccuracies in this problem consideration besides disregard of the angular momentum conservation. To clarify the problem, in the first subsection we consider formulations of the one-particle problem in classical and quantum mechanics and its simplest application to the statistical mechanics. For simplicity, we will restrict ourselves to the case of a two-dimensional gas on a plane perpendicular to the uniform magnetic fieldH=(0,0,H). As will be shown the magnetization of electron gas is nonuniform. We will suppose that the magnetization is small as against the uniform field, and will not regard effect of it. Then magnetic induction His proportional to the external magnetic field strength by the coefficientμ0. Where it is needed, we imply the plane to be of finite “thickness”δz, and, for example, the equations of electrodynamics are written for three-dimensional space.

3.1. Two-dimensional electron ideal gas in uniform magnetic field

This problem traditionally is considered in quasiclassical theory (Lifshitz, I.M. et al., 1973; Shoenberg, 1984). Some corrections will be inserted in this consideration in the section 3.1.1. In this section classical statistic mechanics of ideal gas will be discussed. In the next section the new correction will be obtained from the consistent quantum theory.

3.1.1. Classical statistical mechanics of ideal gas in magnetic field

The Hamiltonian of an electron in a magnetic field has the form:


where m and (e)are the mass and the charge of an electron, pis momentum, and Ais vector potential of the magnetic field:


This Hamiltonian does not have the translation symmetry. This symmetry, seemingly, should be, if the magnetic field is uniform at an unlimited plane. But a uniform magnetic field at unlimited plane is impossible because an electrical current that generates it according to Maxwell equation should envelope a part of this plane. It is asserted (Landau, & Lifshitz, E.M. 1980b; Vagner et al., 2006) that the Hamiltonian (21) with the vector potential (22) would be converted by gauge transformationAA+f(x,y,z). If the functionf=eHxy/2, then Hamiltonian will be in the Landau form :


and will have the translation symmetry in the direction of the axis X in return for axial symmetry. That is strange assertion because the symmetry is the physical property of the system rather than of a method of it description. In fact the transformation to the Hamiltonian (23) in classical mechanics is result of the canonical transformation of the Hamilton variables with the generating functionΦ=pxx+pyyeHxy/2. Then px=px(1/2)eHy, py=py(1/2)eHx,x=x,y=y. Therefore the px,py(in factpx,py) in the Landau Hamiltonian (23) are not the momentum components in the Cartesian coordinates and the absence of the xcoordinate in this Hamiltonian does not lead to the momentum xcomponent conservation. In quantum mechanics the unitary transformation with operator exp[ieHxy/2] is equivalent to this canonical transformation. The boundary is created by the line of intersection of the plane with a solenoid that generates the magnetic field. Electrons, orbits of which transverse this boundary, will be extruded from the area, and the gas will evaporate. The more realistic problem is gas in the area with a reflecting boundary.

An isolated electron has three motion integrals. Those are the angular momentum relative to the centre of area and two coordinates of the centre electron orbit:


Two motion integrals that have the physical importance would be created from it: energy Eand squared centre electron orbit distance from the centre of areaR2:


Here ρ is the radius of the electron orbit. The motion integral R2 is proportional to the Hamiltonian with opposite direction of the magnetic field. The values of R2 and εshould be discrete by the rules of the quasiclassical quantization:


Here λ is magnetic length and ωc is cyclotron frequency. Every energy level En is degenerated because the integer k that determines the position of the orbit can take any values from zero tok1, where k1 is determined by the condition that electron orbit does not have common points with the boundary. That condition would be expressed by formula when the boundary is a circle with radiusR. Then Rkmax+ρn<R orλ(2k1+1+2n+1)<R. Those are the “ordinary” states. Whenk1<k<k2, whereλ(2k2+12n+1)<R, the states are nominated “near-boundary state “. Their energies are not described by the formula (26). The instant magnetic moment is determined by the formula:


The first determination is valid for any negative charged particle with and without an external magnetic field. The second equality is valid when the vector potential has the form (22), or when the equality vi=h/piis transformed by canonical transformation. For the ordinary states the averaged over the orbit magnetic moment isμz=εn/H. The trajectories of the near-boundary states are composed from arcs and envelop the all area. Their magnetic moment is positive. If in the area exists a potential field U(r)the orbit centers of the ordinary states also move along equipotential lines and the energy values depend onk. The degeneration of the energy levels goes off. The angular momentum does not represent the electron motion, but rather the position of the electron orbit. From Eqs. (25-26) it follows thatlz=(eH/2)(ρ2R2)=(nk). Then lz>0 if the orbit envelops the centre, and lz<0if the area centre locates out the orbit. If ρn<<Ra large share of the states with energy εn have angular momentumlz<0.

Going to consideration of the ideal gas with electron-elektron collisions, let us suppose that the interaction does by a central force. Then total energy and angular momentum are conserved. It is generally believed that the area is filled by uniform and motionless positive charged background that neutralizes the electrostatic interaction. This assumption is inconsistently. If electron gas is in equilibrium with motionless background, its angular momentum should be equal to zero. But then it should be nonuniform as is evident from the foregoing consideration. It should be regarded more comprehensively. Let us go to classical statistic mechanics for a gas of charged particles in magnetic field. The characteristic function of the total system (gas and background) is:


Here the indexes gand b denote the quantities those relating to the gas and to the background,H is a Hamiltonian, Eis the value of the total energy. To provide of Gibbs distribution for the gas it is needed to integrate the function φ(Rg+Rb)over the phase space of the background. Then δ(Lg(Rg))is factored out from integral. The function from Hg(Rg)can be factored from integral by the method Krutkov (Krutkov, 1933; Zubarev,1974) asexp(Hg(Rg)/kBT). Then:


The Hamiltonian Hg(Rg)would be represented in the form:


Let us substitute Hamiltonian (30) to the formula (29) and take into account formula i=1Nlzg=Lg and the Eqs. (8). Then it is obtained:


The integration over lz1 leads to change lz1i=2Nlziand after reduction of the quadratic form to the standard appearance the Hamiltonian of ideal gas that is collection of N1 particles in harmonic potential field is obtained. Obviously, that taking into account conservation of the zero value of the angular momentum eliminates the paradoxes that were mentioned at the Introduction. The magnetic momentμz=(/H)(kBTlnZN)0. The gas is confined by the magnetic field. But that confinement provides to the inconsistency of the model that neglects of the electrostatic interaction because the uniform background cannot neutralize it. The model that regards this interaction will be considered below.

3.1.2. Quantum problem of electron in magnetic field at bounded area

The quasiclassical description of an electron in a magnetic field would not give the correct picture of the probability density distribution and the current density. It would not also describe the alternation of the energy spectrum when a perturbation does the classical motion nonperiodical. But that problem in quantum mechanics also is considered insufficiently. As suggested in the paper (Vagner et al., 2006) the density of the probability current of the wave function ψ(r)=exp[iΦ(r)]|ψ| is


The eigenfunctions ψn(r) cannot be chosen so that j=0 because any vector potential Acannot be equal to a gradient of a continuous function. Then it is necessary for stationarity of the states that the current lines to be closed in the area under consideration. The boundary condition best suited to the research of this problem is null of the function on the circumference that bounds the area:

This condition retains the greatest possible symmetry. The current lines in this case are concentric circumferences. The density of the current can be zero only at separate circumferences. Therefore the magnetization is nonuniform. The magnitudes of the eigenfunctions should have the axial symmetry. The localization of the electron cannot coincide with any classical orbit because the uncertainties of values of the orbit centre coordinates (24) should satisfy to Heisenberg uncertainty relation.

The Hamiltonian is as follows:


The operators h^0 and l^zcommutate with each other and each of them with the Hamiltonianh^. The potential energy U(r) is created by the interaction with other electrons and with a neutralizing background. The nonuniform magnetization M(r)should be neglected. If the potential energy U(r)also will be neglected the eigenfunctions of the Hamiltonians h^ and h^0 will have the form:


Here l is eigenvalue of the operator of the angular momentum, Φ(a,c;x)is the degenerate hypergeometric function, Ais normalizing factor that depends on |l| andα. The eigenvalues of energy are expressed by α andl: for the Hamiltonian h^0 that isε0=ωc[α+(|l|+1)/2], and for h^ that isε=ωc[α+(|l|+1)/2]+ωc(l/2). The permissible values of αare determined by the boundary condition. In the common theory (Landau, & Lifshitz, E.M., 1980b; Vagner et al., 2006) that is normability of the eigenfunctions at an infinity plane. Then it is necessary thatα=nr, where nr is integer or zero.Φ(nr,|l|+1;x)=[nr!|l|!/(|l|+nr)!]Lnr|l|(x), where Lnr|l|(x)is the Laguerre polynomial. The energy spectrum of the operator h^0 isε0(n0)=(ωc/2)(2nr+|l|+1)=(ωcn0/2), where n0is integer. The energy spectrum of the operator h^ isε(n)=(ωc/2)(2nr+l+|l|+1)=, =ωc(n+1/2)where n is integer. There are two kinds of a degeneracy of energy levels of those Hamiltonians. The degeneracy of the first kind arise as result of the formulae:n0=2nr+|l|+1, and n=nr+l whenl>0. Every level is degenerated with multiplicity that equals to its number. The perturbation U(r) would be described as a power series without the linear term. Then the term that is proportional to r2 in the spectrum of the Hamiltonian h^0 should change the distance between levels without elimination of the degeneration. In the spectrum of the Hamiltonian h^ the level would be split. The other terms of potential field series eliminate the degeneration in the first order of the perturbation theory. The second kind of the degeneracy is inherent only to the spectrum of the Hamiltonian h^ in absence of any perturbation. Every level is degenerated with infinity multiplicity because when l<0 the energy value does not depend onl. At a bounded area the multiplicity would be limited but very large and would depend on the magnetic field. That will be shown below. This degeneration eliminates by any potential field. The modulo l would have many integer value, and whenl|U|lωc, if that’s possible at somel, the spectrum would be quasicontinuous. That is essential because the explanations of many experimental effects rely on the discreteness of the Landau spectrum. In so doing ones do not study the model stability relative to the electrostatic interaction that is unavoidable perturbation. The other inconsistence of the common model is the large-scale negative total angular momentum of the ground state when all levels have the identical and large-scale multiplicity of degeneracy.

Let us study the degeneration eliminating by the boundary condition (33) in the absence of other perturbation. The polynomial Lnr|l|(x)has n simple zeroes, which, ifnr>>1, are expressed by formula:


Here j(|l|,i) is a null of the Bessel function J|l|(x) with number in the order of increasingi. Then the function ψ|l|(nr) has (nr+1) extremes that are decrease modulo. The extreme number (nr+1) has the sign(1)nr, and further the function tends to zero asymptotically. Obviously, this function cannot satisfy to the boundary condition (33). When α=nrγ the degenerate hypergeometric function has the form:


where Lnr|l|(γ;x)is the polynomial, the nulls of which are less than ζ(nr,|l|;i)by quantities that are proportional toγ, and Tnr(x)is the infinite series that at x>>ζ(nr,|l|;nr) is proportional tox(|l|+1+nrr)exp(x). Then the function (35) will have one more null at the large-scale valuex. This null X tends to infinity when γ tends to zero, and it would be shown that


whenγ<<1. Whenγ1,Xζ(nr+1,|l|;nr+1). Then the boundary condition (33) would be satisfied when(R2/2λ2)=X. (In the mathematical handbook (Erdélyi, 1953) it is written that the function Φ(α,|l|+1;x) at α=nrγ has nr+1nulls, but all these nulls are determined by the formula that is like to the formula (36). Then the boundary condition would be satisfied at arbitrarily large values (R2/2λ2) only when energy has also large-scale value). The value of the null ζ(nr+1,|l|;nr+1) increases when the value of |l| increases. The maximal value|l1|, at which the inequality (R2/2λ2)ζ(nr+1,|l1|;nr+1) is fulfilled, determines the number of the eigenstates of the Hamiltonian h^ that have γ<<1 (ordinary states). It can be shown that when |l1|>>nr then|l1|(R2/2λ2)(nr+1)(R/2λ). The quantum number of the energy n is equal to nr because when l<0 the energy depends on l only byγ(n,l). This number of the ordinary states is consistent with the estimate that was obtained in classic mechanics theory in section 3.1.1 and in the work (Landau, 1930). Those ordinary states are quasi degenerated. The multiplicity of this degeneration is proportional to the magnetic field. Then the main term in the thermodynamical potential, as it was computed by Landau (Landau, 1930), should not depend on the magnetic field, and the Bohr – van Leeuwen theorem is vindicated. When |l|>|l1| and l<0 the other nulls that are described by the formula (36) would satisfy the boundary condition. Those are the near-boundary states, and in this case are not restrictions on the values of γ besides0<γ<1. There are n near-boundary states.

In the section 3.2 it will be shown that for statistical mechanics of the electron gas in the magnetic field the Hamiltonian h^0 has fundamental importance. Let as study the spectrum of this Hamiltonian with boundary condition (33). Then


The degenerate levels are transformed in zonule. It follows from formula (39) that the zonule upper edge is determined by minimum value of (n0nr1)!(nr)! that is roughly[(n0/2)!]2. In the vicinity of nr=n0/2 a shift of energy γ most slowly varies withnr. It means that a density of states is the highest in the vicinity of the zonule upper edge. The zonule lower edge is shifted byγmin=γmax[(n0/2)!]2/n0!<<γmax. The zonule width γmax(n0)increases with n0 if(R2/2λ2)>(n0/2). If γmax(n0)<1/2 then number of states on each interval of energy values of width ωc/2 is equal ton0, (the spin will be taken into account subsequently). For γmax(n0)>1/2 zonules overlap, gaps in the spectrum disappear, and the number of states in the interval becomes less thann0, i.e., grows of the density of states is decelerated with energy increase. That is transitive area of energy. For higher values of energy not the greatest but other nulls, which are described by formula (36), satisfy the boundary condition. With the approximate formula for the nulls of the Bessel function j(|l|,ν)π(ν+|l|/2)the second in magnitude null of the function (35) with ν=nr is[π2n02/8(n0+2γ)](π2n0/8). When it will be so that π2n0/8R2/2λ2 the boundary condition should be satisfied only not the greatest but other nulls of the degenerate hypergeometric function. Corresponding values of energy look like:


where ν=nrn/2 and γ is not described by the formula (39). They coincide with eigenvalues energy in a circular potential well with reflecting boundaries (see formula (19)). The spectrum becomes quasicontinuous as distances between the nearest levels are proportional toR2. The density of states does not depend on energy, as well as for two-dimensional gas of free particles. The function (35) would be expanded in series over the Bessel functions (Erdélyi, 1953) and the first term is proportional toJ|l|((r/)2εm). The parameter of this expansion is(n0)1/2. Hence, the wave functions in this approximation also coincide with free electron wave functions. This form of the spectrum would be illustrated by the classical consideration. If the formula (34) would be considered as the classical Hamiltonian, the ultimate energy for which the classical accessible area is determined by the parabolic potential is equal tomωc2R2/8. That within a factor ξ that is close to 1 coincides with the energy of the transitive area and with energy of the transition to quasicontinuous spectrum(2/π2)mωc2R2.

These results can be described as energy spectrum breakdown into two bands. A spectrum lower part is denoted as a magnetic band, and the upper one will be denoted as a conduction band. Bands are not separated by a gap or sharp boundary, but far from transitive area the density of states and wave functions differ substantially. Fine structure of the density of states in the lower part of magnetic band represents the narrow zonule separated by gaps. The total width of the allowed zonule and gap is equal toωc/2. Number of states in an interval with number n0is equal ton0. In the transitive area gaps disappear and in the conduction band the spectrum is quasicontinuous. The maximum of magnitude of a wave function in the magnetic band is localized within a ring of width about λ2 and aboutR=λ2n0 in radius. The density of states in the magnetic band, averaged over the intervalωc/2, is 2n0/ωc=4ε0/(ωc)2 and grows with energy. In the transitive area this grows is decelerated, and in the conduction band the density of states does not depend on energy and is equal to mR2/22(without spin consideration). It is believed that the magnetic band is finished when n0=R2/4λ2=nb or ε0=mωc2R2/8=εb as would be expected from the quasiclassical consideration. Then the density of states will be continuous.

3.2. Statistical mechanics of electron gas in uniform magnetic field with regard for electrostatic interaction

The quantum-mechanical average value of the magnetic moment in the ordinary eigenstate of the Hamiltonian h^ in absence of the potential energy U(r)is:


It is a negative quantity because the positive term that proportional to (γ/H)is small. In the paper (Landau, 1930) were taken into account only ordinary states. Then the quantum-statistical average value of the gas magnetic moment must be:


Here n(n,l) is the average occupation number of the stateΨ(n,l), D(H)is the multiplicity of degeneracy that in this case does not depend on the state energy and depend on the magnetic field. But in this work the magnetic moment of the gas was computed as


When this result is compared with the formula (42) it is apparent that the thermodynamical potential Ω is determined incorrectly.

Let us obtain the density of distribution for an electron gas that is at equilibrium with thermostat, which is described by classical mechanics. The conservation of the zero value of the angular momentum also will be taken into account. The characteristic function of this system is:


This formula is obtained on a basis of the properties of characteristic functions that was described in the formulae (68), and the quantum characteristic function (see formula (15)). Here E is a total energy, Eis the energy of the electron gas, Hthis the Hamiltonian of the thermostat that is a summatory function of the classical variables. Values of energy are expressed by the dimensionless positive integers, and the angular momentum is measured in. The expression of the Hamiltonian h^ as the sum h^0+(ωc/2)l^z is taken into account. Then the second equality in (44) would be rewritten as:


by using the formula (8). Let us generalize the Krutkov method (Krutkov, 1933; Zubarev, 1974) for this case. We calculate the Laplace transformation with respect the total energyE. To do this let us multiply the function (2π)1IE (formulae (44 - 45) by exp(ϑE)and integrate over E between 0 and. Then the integrations over E and τ can be performed, and the result is:


where L{f(E),ϑ}=Lϑ[f(E)] is Laplace transform of a functionf(E). The integration over the thermostat phase space leads to the formula:


where Γ1 is the phase space of a thermostat particle and hth is its Hamiltonian and Nthis the number of thermostat particles. As result of inverse Laplace transformation we obtain:


This integration can be performed by the saddle-point method because Nthis a large-scale number. The saddle-point ϑ0=(kBT)1 is determined by the first exponent, and the second exponent that depend on variables-operators of the quantum subsystem can be factored from integral at ϑ=(kBT)1as well as in the classical case.The result should be substituted into (44). Then the non-normalized statistical operator is obtained:


Here the first multiplier is the common formula of the statistical operator for the quantum Gibbs distribute (Landau, & Lifshitz, E.M., 1980a; Zubarev, 1974). The second multiplier would be computed by the Darwin – Fowler method (Fowler, & Guggenheim, 1939) like as in the formula (16) and describes the conservation of the particle number. If the grand canonical ensemble is considered, then the statistical operator of particle number


would be obtained from this multiplier by the Krutkov method. The last multiplier in formula (49) cannot be computed by those methods because it does not have any large-scale parameter. This multiplier imposes constraints on ensembles that the total angular momentum equal to zero. If Hamiltonian and an operator that should be averaged have the commutative term that is proportional to the total angular momentum operator, this term should be eliminated when the averaging is performed. That is the reason for the change ε to ε0 in the Hamiltonian Ψε0(Ψ)aΨ+aΨ in the formula (48).

The model that will be considered below is described by the Hamiltonian:


Here U(r) is potential of the self-consisted electric field that describes interaction of an electron with other electrons and neutralizing background that have the form of circle with radiusR.

The electron density in the magnetic field should be distributed in such a way as to shield the external potentialmωc2r2/8. This shielding cannot be perfect. The current that would be generated by a residual potential would be compensated by the diffusion current that generated by an inhomogeneity of the electron density. Because the kinetic energy in the Hamiltonian (51) has the standard form, electrons density distribution and self-consistent potential can be calculated by the density functional method, described in the work (March, 1983). In this approximation, energy of the ground state of electron gas is presented in the form of functional of the gas density, i.e., number of particles per unit of areaρ(r):


Here π2ρ2/2m is density of kinetic energy of two-dimensional degenerated Fermi gas (each state is supposed to be twice degenerated), the second addend describes energy of electrons Coulomb interaction, e0is vacuum inductivity, the third addend is energy of electrostatic interaction with the neutralizing background, the last addend described gas potential energy in the effective harmonic potential depending on magnetic field. The neutralizing background has the charge density eρ0Y(Rr)and generate the potentialeu(r). This potential would be expressed by the elliptical integrals and is approximated by the parabola and a repulsing bound. The step function Y(x)=1at x>0 and Y(x)=0 atx0. The constant eu0would be assigned so that the potential energy in the effective Hamiltonian would be zero atr=0. All integrations are made over the areaC, which is a circle of radiusR. The functional E0[ρ(r)] should be minimized under the following supplementary condition:


The quadratic term of the residual potential is of chief interest. Therefore the minimization would be performed in the quadratic approximation. The test function has the form:


and as result we obtain:


where a0=4πe02/me2 is Bohr radius. The Lagrange multiplierζ, which has meaning of the chemical potential or the Fermi energy is obtained:


where ζ0 is the Fermi energy of two-dimensional electron gas in absence of a magnetic field. The residual potential energy has the form(4a0/3R)(mωc2r2/8). It is expressed conveniently by new frequencyωr=ωc4a0/3R. Then one-particle effective Hamiltonian will have the form that was study in the previous subsection:


The Coulomb interaction and the electron density inhomogeneous that are neglected commonly decrease the frequency ωcωror another way the effective charge eer by factor4a0/3R104. This factor does not depend on the electron density if this density satisfies to the condition of the degeneration. The operator of the total magnetic momentM^=(er/2m)iN(l^i+erHri2/2). When the quantum-statistical average is computed, iNl^i=0, andM=M^=iNh^r/H. The energy magnetic splitting over spin does not alter. Hence it should be considered separately.

The electron gas always interacts with electromagnetic field that in the ordinary circumstance has zero temperature. This leads to the fact that the gas passes to the ground state by the spontaneous photon irradiation. The role of statistical mechanics is that it imposes a constraint on a value of the angular momentum in the ground state. The energy of the ground state with zero angular momentum would be computed by the spectrum of the one-particle states that was described in the previous subsection. The energy levels, the degeneracy multiplicity, and the boundary of the magnetic band are determined only the orbital motion. Then the energy of the magnetic band is:


The number of states in the magnetic band equals to:


The quantity nb is discrete; therefore the formulae (58 – 59) describe a smoothed function. The distance between the neighboring spikes in Emb isωrnb2. If Nb is less than the number of electrons N=πR2ρ0 the energy of the ground state with a glance the formula (56) is:


Here the first term is the energy of the electron gas in the absence of a magnetic field. The second term that makes the main addend in the magnetic moment is the product of the Fermi energy of the electron gas ζ0and the additional term in it, which depends on the magnetic field (formula (56)). The third term in the formula (60) is created by the addendsEmb, the squared additional term of the Fermi energy, and the negative addend in the result of integration(mR2εb2/22). The last term in this formula describes the energy lowering by the spin polarization. It always is smaller than the second term. The magnetic moment is:


whereEp=(eRH)2/4m. The first term is the diamagnetic moment of the orbital motion in the states of the conduction band. It is proportional to the electron density and third power of the radius as distinct from the result of the Landau theory (see to formula (5)). This formula would be obtained by other way. Let us consider the density of electric current (formula (32)) in the ground state with zero angular momentum:


The magnetic moment is:


The first term equals to zero because 2π0R|ψnl(r)|2rdr1 andnln(n,l)l=0. Let us change in the second addendnl|ψnl(r)|2n(n,l)=ρ(r). Then the integration gives both the first and the second terms of the formula (61). If the sum nln(n,l)l will be computed in the common theory, the linear with respect to the magnetic field terms in the formula (63) will cancel.

The orbital diamagnetic susceptibility decreases with increasing of the magnetic field. This function also has a spikes that is caused by the addend2Emb/H2. When Nb will be equal to N=πR2ρ0 the diamagnetic susceptibility will be zero. When the magnetic field is stronger, the Fermi energy remains on the level nF=N of the magnetic band and increases proportionally to magnetic field. The magnetic moment remains constant. The paramagnetic moment also does not increase because the density of states on the Fermi level decreases linearly. The value of the magnetic induction whereby the saturation of the magnetic moment will start is:

Ifρ0=1019m-2, R=0.1m, then Hs=5,15TandMs=2.92102JT-1.

4. Conclusion

The fundamental theory of statistical mechanics requires taking into account the law of the angular momentum conservation. The fulfilment of this requirement does not introduce any essential alterations into statistical thermodynamics, when the angular momentum of the system equal zero and the system Hamiltonian is a positive definite quadric form of all momenta. An equilibrium isolated system would have nonzero angular momentum only if an attraction of particles can resist centrifugal forces, as it is in nebulas. A gas can be in equilibrium with a rotating envelope that is a termospinstat. The condition of this equilibrium is the equality of the average value of sum of particles angular velocities to the angular velocity of the envelope. The Gibbs density of distribution and the thermodynamical functions are generalized for this case. If a system has the angular momentum equal to zero, the conservation of this value is important only when the Hamiltonian or/and an averaged quantity depend on the angular momentum. The problem of an electron gas in a uniform magnetic field is considered with taking into account the conservation of the zero value of the angular momentum. This consideration eliminates the paradoxical statement of the conventional theory that diamagnetic moment of the gas equals zero in classical as well as quantum physics (the Bohr – van Leeuwen theorem). The new formulae for the magnetic moment of the electron gas are obtained. It also leads to the effect of confinement of two-dimensional gas of charged particles by magnetic field. This results in effect of a non-uniform density of a gas, which decreases with distance from a center according both to classical as well as quantum theory. Then the model of noninteracting charged particles does not have areas of application. Many theories should be reconsidered, if they are founded on this model and on the statistical mechanics which does not take into account the angular momentum conservation law.


1 - A. A. Abrikosov, 1972 Introduction to the Theory of Normal Metals, Nauka, Moskow.
2 - K. Cohen, 1951 The Theory of Isotope Separation as Applied to the Large Scale Production of U235, McGraw- Hill, New York.
3 - A. Erdélyi, 1953 Higher Transcendental Functions, based, in part, on notes left by Harry Bateman, 1 Mc Graw-Hill book company, INC, New York Toronto London.
4 - R. P. Feynman, R. B. Leighton, M. Sands, 1964 The Feynman lectures on physics, 2 Addison-Wesley Publishing Company, Inc. Reading, Massachusetts. Palo Alto. London
5 - R. H. Fowler, E. A. Guggenheim, 1939 Statistical Thermodynamics, Cambridge University Press, London.
6 - A. Y. Khinchin, 1949 Mathematical Foundation of Statistical Mechanics, Ed. Dover, NewYork.
7 - A. Y. Khinchin, 1960 Mathematical Foundation of Quantum Statistics, Ed. Dover, NewYork.
8 - Y. A. Krutkov, 1933 Zs. Phys. 81 377 & Supplement by Editor, In translation into Russian of the book by H. A. Lorentz, Statistical Theory in Thermodynamics, (1935), ONTI, Leningrad- Moscow.
9 - L. D. Landau, 1930 Diamagnetism of Metals, In: Collected papers, 1 47 (1969), Nauka, Moscow.
10 - L. D. Landau, E. M. Lifshitz, 1980a Statistical Physics (3rd rev., part I), Pergamon Press, New York.
11 - L. D. Landau, E. M. Lifshitz, 1980b Quantum mechanics. Nonrelativistic theory, Pergamon Press, New York.
12 - I. M. Lifshits, M. Y. Azbel, M. I. Kaganov, 1973 Electron Theory of Metals, Consultant Bureau, 0-30610-873-9 York.
13 - March N. H. , 1983 Origins: theory by Thomas- Fermi, In: Theory of the inhomogeneous electron gas, Ch. 1, S. Lundqvist & N.H. March (Ed.), Plenum Press, New York & London.
14 - R. W. Pohl, 1960 Elektrizitätslehre, Springer- Verlag, Berlin- Gottingen- Heidelberg.
15 - D. Shoenberg, 1984 Magnetic oscillations in metals, Cambridge University Press, Cambridge.
16 - van Vleck J.H. 1965 Theory of Electric and Magnetic Susceptibilities, Oxford University Press, Oxford.
17 - G. Uhlenbeck, G. Ford, 1963 Lectures in Statistical Mechanics, American Mathematical Society, Providence, Rhode Island.
18 - I. D. Vagner, V. M. Gvozdikov, P. Wyder, 2006 Quantum mechanics of electrons in strong magnetic field. HIT Journal of Science and Engineering, 3 1 555 ,Holon Institute of Technology.
19 - D. N. Zubarev, 1974 Nonequilibrium Statistical Thermodynamics, Consultant Bureau, New York.