Symmetry Aspects of the Band Structure and Motion Equations Applied in Calculating the Cyclotron Frequency of Electrons in Metals

Cyclotron frequency of a crystal electron is, in general, not an easily accessible parameter. Nevertheless, its calculation can be simplified when the symmetry properties of the band structure and those of the motion equations in the magnetic field are simultaneously taken into account. In effect, a combined symmetry of the electron Hamiltonian and that of the Lorentz equation provide us with a non-linear oscillator problem of high symmetry. In the next step, the kinetic energy of the oscillator can be separated from the whole of electron energy and applied in a new kind of calculation of the cyclotron frequency which is much more simple than before. In consequence, a detailed approach to the electron circulation, also in more complex band structures, becomes a relatively easy task. For different crystal lattices of cubic symmetry taken as examples the cyclotron frequency of the present and a former method are compared numerically giving the same results.


Introduction. Outline of the Symmetry Properties Concerning the Crystal Electrons Moving in an External Magnetic Field
In theoretical physics one of the principal aims is to calculate the observed data with a possibly high accuracy.Simultaneously, another aim is that this approach to experiment should be of a not too complicated nature.In the solid state theory these two goals are strictly coupled together by the use of the properties of symmetry.
Conventionally, the symmetry in solids is connected with the crystal geometry given by a regular spatial arrangement of the atoms.This is a well-known simplification entering the potential part of the electron Hamiltonian.However, when an external field is acting on a crystal, much of the original symmetry due to the crystal potential is usually lost.In principle, this situation can lead to difficulties, see e.g., [1].But if the potential perturbation, for example that of the magnetic field, is chosen in some special direction with respect to the crystal lattice, a loss of symmetry in the geometry of the crystal potential is-to a large extent-compensated by a new kind of symmetry, namely that connected with the dynamics of electrons moving in the magnetic field.
The purpose of the present paper is to demonstrate that this kind of symmetry, which can be next applied in calculating the cyclotron frequency in solids, is due to both the properties of the Lorentz equation governing the electron motion in the presence of the magnetic field and the wave-packet character of electrons contained in the band of states of a solid.The final effect of such a combination is a pair of equations much similar to a pair of the Hamilton equations possessed by a one-dimensional anharmonic oscillator.The variables of the oscillator, which can be associated with the components of the electron wave vector, have-on one hand-the symmetry properties of the dynamical variables of position and momentum of the oscillator.On the other hand, the remainder of the symmetry details are dictated by the crystal Hamiltonian.
In fact, the Hamilton equations coming from the band structure combined with the action of the magnetic field have been derived some time ago by Wannier [2] and Suhl [3].But these equations were, in our opinion, not sufficiently exploited in describing the motion of the crystal electrons due to the presence of the magnetic field.An attempt to calculate the cyclotron resonance for such electrons has been done in [3], but it was limited to the s-electrons contained in only one, namely simple cubic, crystal lattice.The restriction is perhaps because the formalism applied in [3] to the simple cubic case could not be extended to other lattices.But the experimental interest is focused on crystals rather different than the simple cubic one, moreover the band structure is based on a different kind of the atomic orbits than s, see [4,5].This provides us with a practical aim of the present paper to give a method directed to bridge a definite computational gap.
Cyclotron frequency Ω is a fundamental electron parameter entering practically almost any examination of a conducting solid at its Fermi surface in case when an external magnetic field is applied.This frequency is usually easy to assess when the surface is near to a spherical, or ellipsoidal, shape.For, in this very special circumstance, the influence of the crystalline background exerted on the electron states is rather small and the frequency formula is regulated mainly by a suitable choice of the effective electron mass, see e.g., [6].An extremely simple example is that of free electrons (see below).
In general, when the motion is represented in the space of the wave vector k, a conventional attempt is to calculate the planar area S enclosed by a trajectory of an electron having a constant energy in that space, and next differentiate that area with respect to the energy E possessed by the circulating electron.This gives the reciprocal value equal to the cyclotron frequency [6,7]: For free electrons we have S = πk 2 and E = 2 2m k 2 , thus the frequency (1) becomes which is a well-known result.A characteristic point is that ( 2) is a constant independent of the electron energy and k.But a circulation of electrons in the crystal lattices does not provide us with so simple results.Consequently, in order to approach that problem in a possibly simple manner, a different formalism than that represented in (1) is applied.This takes into account not only the symmetry of the crystal Hamiltonian but also the properties of symmetry characteristic for the motion equations in the magnetic field.

The Formalism
A planar electron motion performed by an electron in a constant magnetic field B is governed by the Lorentz equation.Its essence is that if the action of an external electric field is negligible [8] we have This equation describes the motion of the electron wave packet labeled by the wave vector k and having in a real space the velocity v.It is valid for not too strong B, so the band structure is not essentially perturbed by its degeneracy into Landau levels.In fact, (3) is equivalent to a pair of the Hamilton equations of classical mechanics.This is so because for the field B directed, say, along the axis z, Equation (3) can be presented as an equation pair: The components k x , k y of the wave vector behave in (4), (5) like position and momentum variables, say because of the well-known property [7] where E is the electron energy.It should be noted that the one-electron Hamiltonians for crystals in the magnetic field based on the components of vector k seem to be developed first in [9].For the sake of simplicity the constant parameters entering (3)-( 7) can be abbreviated to with the substitutions ( 6) and ( 8), Equations ( 4) and ( 5) become: Equations ( 9) and ( 10) simplify essentially the problem, because the energy expressions E(k x , k y ) for an electron circulating in a (k x , k y )-plane (component k z behaves like a constant parameter) are rather easily accessible quantities.For example, in the tight-binding approach to the s-electrons in cubic crystal lattices, we have the energies [10] for the simple cubic (sc) lattice the nearest neighbors of the central site in (0, 0, 0) are here in positions (±a latt , 0, 0), (0, ±a latt , 0), (0, 0, ±a latt ) for the body-centered cubic (bcc) lattice the nearest neighbors are in positions (±a latt , ±a latt , ±a latt ), for the face centered cubic (fcc) lattice [the nearest neighbors of the central site are in positions (±a latt , ±a latt , 0), (0, ±a latt , ±a latt ), (±a latt , 0, ±a latt )].
The β are multiples of the hopping integrals β latt between the nearest atomic neighbors in the crystal lattices, for example in the sc lattice β = 2β sc , the a latt is the edge length of the elementary cube of a crystal cell.The constant terms of 1 and 3 are introduced in order to make (11), (12) and (13) at small k x , k y and k z proportional to the free-electron expression: In this case the constant term 2 /2m entering the free-electron energy in (1) is replaced by The parameter β has dimension of energy.A substitution of (14) into (1) replaces the frequency (2) by In Sections 3-8 we discuss and apply the Hamiltonian expressions of ( 11)-(13) assuming

Crystal Hamiltonians as the Oscillator Problems
When the motion is limited, for the sake of definiteness, to a plane the energies (11), (12) and (13) become, respectively, equivalent to the following Hamiltonian expressions on the condition that the substitutions presented in ( 6) and ( 17) are taken into account.
A characteristic point of ( 19)-( 21) is that beyond the symmetry properties dictated by the crystal potential [see ( 11)-( 13)], a supplementary symmetry in the dynamic variables x and p x , namely that due to ( 9) and (10), occurs in the formalism.In effect of that symmetry, the Hamiltonians (19), ( 20) and (21) provide us with an anharmonic oscillator problem [11].The oscillators of this kind are one-dimensional, because x and p x refer to the motion in a single dimension, say, that extended along the axis x.In consequence, the electron circulation along the edge of a planar cross-section area of the Fermi surface has been replaced by an oscillation between two opposite (symmetrical) points of the axis x: the axis is lying in a cross-section plane of the Fermi surface labeled by (18).The distance a 0 in (22) is called the amplitude of the oscillator motion.A property of the oscillator is that when x = k x = ±a 0 we have p x = k y = 0, and vice versa: for x = k x = 0 we have p x = k y = ±a 0 .
In course of its oscillatory motion the electron has a constant energy equal to that possessed by an electron circulating in the (x, y)-plane.Consequently, this energy can be considered, for the sake of simplicity, at some special points of the oscillator trajectory, for example those quoted in (22).The values of x, p x entering that points, when substituted into (19), (20) and (21), give respectively the planar Hamiltonians: for each of the considered lattice cases.Our task becomes to calculate the oscillation frequencies for the electron motions represented by the Hamiltonians given in (23), ( 24) and (25).

Symmetry of Oscillators with Respect to the Dynamic Variables x and p x and Separation of Energy into the Kinetic and Potential Parts
A total electron energy of a one-dimensional oscillator is where the first term on the left-hand side is a kinetic energy, the second one is a potential energy, and C latt (a 0 ) is a total energy of the oscillator which is a parameter constant in time.A differentiation of (26) with respect to time gives dx dt from which we obtain the well-known relation The left-hand side of (28) can be easily calculated on the basis of the Hamilton equations presented in ( 9) and (10).These give the acceleration of x: Symmetrically, we have the acceleration of p x : A characteristic point obtained from (30) and ( 31) is that The Expressions ( 29)-( 30) can be calculated in each lattice case.In the next step, on the basis of the energy Equations ( 23)-( 25), the terms dependent on p x entering (29) can be replaced by those dependent on x, and the constant terms enter also the energy expressions.In effect, the acceleration term on the left-hand side of (29) becomes a function of x alone and the constant terms.According to (28), its integral over the variable x leads to V (x) with the accuracy to a constant potential term V C .
The term V C can be obtained from the condition that the points are the turning points of the oscillator.In this case the electron velocity dx dt as well as the kinetic energy are equal to zero.Therefore, This completes the calculation of the potential energy.Examples are presented below separately for each cubic lattice.
But the same reasoning can be applied also to the acceleration expression for the momentum variable p x in (30).In this case the terms dependent on x are replaced by those dependent on p x combined with the constant terms entering the energy Equations ( 23)-( 25).This gives the acceleration terms dependent only on p x and the constant terms.We find that both functional dependencies obtained for (29) and (30) are the same.
In effect, we can write for (29) and (30) respectively and This kind of symmetry, represented by the same functional dependence of both accelerations ( 35) and (36) on x and p x , allows us to derive another equation of motion of a one-dimensional nonlinear oscillator than that given in (26).The new equation concerns the behavior of the momentum variable p x and has its postulated form similar to (26): The potential V depends solely on the momentum variable p x and the constant terms.The amplitude a 0 labels now a maximum value of |p x | which can be admitted by the oscillator.The time differentiation of (37) gives Examples given below (see Sections 5-7) demonstrate that a functional dependence of V (p x ) is the same as V (x).The constant terms entering V (p x ) and V (x) should be the same because a 0 is equal to the same amplitude for the oscillations of x and p x .
The kinetic energy of the oscillator in the momentum space is given by the formula: Examples of (40) for special lattice cases are given in Section 5.

Example 1: The Simple Cubic Lattice
the acceleration term for electrons oscillating in that lattice is : because of ( 10) and (19).From Equation (23) we have This gives The constant term V C is obtained from the condition of vanishing of the kinetic energy at the turning points (33) [see (34)]: The validity of calculations can be checked by substituting the kinetic energy derived from the first of the Hamilton Equations ( 9) and V (x) from ( 45) and ( 46) into (26).We obtain This is a constant result expected on the basis of the energy Equation ( 23).It shows that relations hold for any x of the oscillator.The same energy expression in the momentum space is: where V (p x ) depends on p x in the same way as V (x) on x and the constant V C remains unchanged.
In Figure 1 parts (a), (b), (c) we demonstrate the phase space (x, p x ) calculated for three forms of the electron Hamiltonian for the sc lattice.Figure 1(a) is representing a constant energy C sc = 1 − cos a 0 calculated for two values of cos a 0 in the case of the Hamiltonian (26) which is fully symmetrical in x and p x .Figure 1(b) plots the same energy of C sc in the case of the Hamiltonian formula given in the first step of (47), Figure 1(c) is a plot of C sc obtained when the kinetic energy is presented according to (48a).

Example 2: The Body-Centered Cubic Lattice
The velocity of the oscillator representing the circulating electron is and The corresponding acceleration term is: In this special (bcc) case ( 51) is a function of x alone.The potential energy obtained from the integral of the negative acceleration term is: The constant V C can be derived from the energy balance at the turning points at which the oscillator velocity is zero because of p x = 0; see (49).This gives the condition: A control of the calculations requires to add the kinetic energy to the potential one calculated in (52)-(54).In (55) we eliminate the variable p x with the aid of the relation cos x cos p x = cos a 0 (56) obtained from (24).A sum of (55) and (52) supplemented by the result in (54) for V C gives which is a constant energy term entering Equation (24).In the momentum space a counterpart of (55) is See (55) and a similar change of the variable x into p x concerns the expression of V (x) in (52).The constant V C remains unchanged.In Figure 2 we plot the constant energy C bcc = 1 − cos a 0 in the phase space (x, p x ) for the case of a fully symmetrical Hamiltonian formula (24).But an identical plot in (x, p x ) is obtained for C bcc when E kin = E kin (x, p x ) calculated from the first step in (55) and V (x) from the second step in (57) are taken into account.By the property of symmetry the calculation of C bcc done with the aid of the kinetic energy (55a) and V (p x ) taken instead of V (x) gives a result identical to that of the former two approaches, viz., (24), or (55) and (57).

Example 3: The Face-Centered Cubic Lattice
Here and Figure 2. The (x, p x ) plot of a constant energy C bcc of the body centered cubic lattice; two values applied for cos a 0 are: 1/2 and 1/4.The plot obtained on the basis of (24) degenerates with the energy plot having E kin taken from the first step in (55) and V (x) taken from the second step in (57), as well as with the plot based on E kin taken from (55a) and V (p x ) considered instead of V (x).A further calculation requires to eliminate the dependence on p x on the basis of the energy formula in (25).We obtain which can be substituted into (60).This gives: The expression in (34) taken with a minus sign should be next integrated over the variable x.The result is where The constant V C is obtained from the requirement of the energy conservation at the turning points see (25): The kinetic energy in the fcc lattice is from (58): This energy added to V (x) gives a constant value equal to C fcc (a 0 ) see ( 25): Symmetrical expressions to those in (66) and (67) can be obtained in the momentum space by a replacement-in the first step of (66) and in (67)-of the variable p x by x, and the variable x by p x .In Figure 3 parts (a), (b), (c) we plot the constant energy C fcc = 1 − cos a 0 in the phase space (x, p x ) calculated: (i) from (25) (see Figure 3(a)), (ii) from E kin obtained in the first step of (66) and V (x) is taken from (67), (iii) from the formula (37) specified for the case of the fcc lattice.The degeneracy of the plots observed for the bcc case does not exist for the fcc lattice.3(b) is calculated on the basis of E kin obtained in the first step of (66) and V (x) taken from (67); Figure 3(c) is a plot of the energy formula (37) specified for the case of the fcc lattice.

Equations of Motion as a Basis in Calculating the Cyclotron Frequency
This kind of approach replaces that quoted in Equation ( 1).
On the first step, the motion of crystal electrons close to a free-electron behavior attained at small k x , k y , k x is examined; see (14).In this case the energy expressions (23)-(25) of a planar motion (k x = x, k y = p x ) are reduced to the formula This is a typical Hamiltonian of the harmonic oscillator.In accordance with the oscillator property at x = ±a 0 is attained p x = 0, and vice versa, p x = ±a 0 is attained at x = 0.
The electron velocity is On the basis of (69) a half of the time period of the oscillator can be readily calculated from the integral In effect In case the potential V (x) instead of the velocity (69) is considered, we have as a starting point the relation 1 2 The V (x) in (72) is obtained readily from the acceleration formula since the integral of the potential derivative over the variable x gives: At x = ±a 0 we have dx/dt = 0, so the potential energy should be equal there to a total energy (68).This gives the equation so V C = 0.In effect, the expression under the integral with respect to x leading to T becomes identical to that applied in (70).Since the general formula for the circular frequency is we obtain in view of ( 70) This is precisely a result for Ω obtained in (2), on the condition that the values of parameters ( 8) are substituted into (2).
The situation for the crystal electrons becomes much similar to that for the free ones.Because the kinetic and potential energies add together to a constant value for any x, it is inessential whether we apply the kinetic, i.e., the velocity-dependent, energy expression, or the potential energy in the calculations.In the velocity case, since can be applied for any lattice.For example, for the sc lattice we have see ( 41) where the transformation of the dependence on p x into that on x is obtained on the basis of (43).But the potential V (x) can be also used.In this case, because of the relation the integral formula becomes: An equivalence of the both approaches can be easily demonstrated.For the bcc case we have [see (52)] where the second and the last step are obtained because of the formulae (54) and (55).
A similar calculation can be done for the fcc lattice, since from (25), ( 63) and (66) we have:

Cyclotron Oscillation Periods and Frequencies
The effect of the crystal potential is represented by the dependence of Ω on the amplitude a 0 absent in (2), ( 16) and (77).These Ω are reduced to that obtained in (77) only if a 0 is tending to zero, a situation obtained when the crystal electrons approach the behavior of the free-electron particles.
With a substitution of the electron oscillation in the sc lattice can be represented by the formula With another substitution k = sin a 0 2 (88) the integral in (87) becomes: Assuming that 2z = w the last integral becomes With the aid of the formula (76) this leads to the expression for the cyclotron frequency where The expansion of (89) in powers of a 0 is presented in Table 1.This expansion is equal to that calculated earlier with the aid of more conventional, but also more complicated, methods [11].
The bcc problem is equivalent to an elliptic integral similar to that met in the problem of the mathematical pendulum: A substitution of sin x = sin a 0 sin z (92) provides us with the formula Because of (76) the cyclotron frequency in the bcc lattice becomes its power expansion in terms of a 0 is also presented in Table 1.
sc lattice The period of the electron oscillation in the fcc lattice is given by the integral: provides us with the following expression which evidently is also proportional to an elliptic integral.The cyclotron frequency in the fcc lattice is see (76): A power expansion of it is also given Table 1.The accordance of the expansions in Table 1 with those calculated formally with aid of other methods [11] becomes evident.A physical cyclotron frequency is obtained when ω latt of Table 1 is multiplied by a factor presented in (16); see the end of Section 1.
An extension of the frequency calculations for the sc and bcc lattices to the case of k z = 0 is straightforward.If the energies E sc in (11) and C sc in (23) are equal, i.e., they concern the same Fermi surface [see also (17)], we obtain for k x = x, k y = p x the relation where a 0 the amplitude of the oscillation on the plane k z = 0. Equation (41) holds also in this case, therefore the formula (89) for the cyclotron frequency is valid with a 0 replaced by a 0 .In effect, the a 0 entering the expansion for the sc in Table 1 should be replaced by A similar calculation can be done for the bcc lattice.Because of ( 12), ( 17) and (24) we have where is the formula for the amplitude a 0 at k z = 0.In this case (49) becomes In effect, an infinitesimal oscillation time becomes inversely proportional to cos k z .The formula for the oscillation frequency is similar to that in (94), namely where from (101) The calculation of the oscillation frequency at k z = 0 for the fcc lattice is based on the energy formulae (13), ( 17) and (25): This gives from which where with we obtain Another substitution gives By applying the formula [12] π/2 a following elliptic integral defines the time period T : from which the calculation of the cyclotron frequency ω fcc is straightforward.Following (76), any time period of the oscillation can be next transformed into the frequency Ω.These Ω differ from the result of unity for free electrons [see (77)] by the power expansions presented in terms of a dimensionless amplitude parameter a 0 , or a 0 and k z .Physical values of Ω are obtained by multiplying the dimensionless Ω by a factor equal to (16); see the end of Section 2.

Extensions of the Method to More Complicated Hamiltonians
The Hamiltonians discussed above were symmetrical with respect to an interchange of the variables x and p x .The amplitude a 0 can be referred to a crystal state (k x , k y ) lying in a plane k z = 0 by the formula where latt = sc, bcc, fcc.
More complicated Hamiltonians can be obtained, for example, when more atomic layers than that of the nearest atomic neighbors are considered.Let the second atomic layer interacting with a central atom be assumed as also important for the band energy in the sc lattice.The atomic distribution in that layer is similar to a layer of the nearest atomic neighbors in the fcc lattice.In this case the oscillation amplitude is referred to the electron energy by the formula where β I is characteristic for the nearest neighbors in the sc lattice and β II for the nearest neighbors in the fcc lattice.Since the interaction between atoms decreases strongly with an increase of the interatomic distance, we may assume that The oscillator Hamiltonian for two layers taken together into account becomes This gives the velocity which vanishes at p x = 0, and x = a 0 is a new amplitude characteristic for this case.The oscillation problem and calculation of the cyclotron frequency become similar to those considered for the Hamiltonians discussed above.
But not only the spherical s-like atomic orbits define the band structure of numerous crystals.In the case of orbits different than s, an asymmetry of the Hamiltonian describing the electron motion in the (k x , k y )-plane should be usually taken into account.For example, the wave functions based on p x -like atomic orbitals (this symbol should not be confused with the electron momentum) can be considered, see e.g., [10].In the sc lattice we obtain then different interaction between the nearest atomic neighbors lying in two sites (±a latt , 0, 0) than with the neighbors lying in two sites (0, ±a latt , 0); in the nearest neighbors approximation the interaction with the remaining two sites (0, 0, ±a latt ) becomes inessential for the planar motion in the (x, y)-plane.In this case the Hamiltonian for three dimensions becomes Here β x refers to the interaction between the nearest atomic neighbors in the x-direction and the ratio refers the β x interaction to the interaction of the central atom with the neighbors placed in the y and z directions.Assuming the oscillator motion along axis x in the plane k z = 0 is considered, the Hamiltonian for that motion which replaces (19) becomes : since β in (17) can be replaced by β x .The last step in ( 125) is obtained because for p x = 0, and in this case x = a 0x .The kinetic energy of the motion calculated with the aid of ( 125) and ( 126) is: The potential energy obtained from the acceleration formula The integration constant V C can be obtained from the condition that at x = a 0x the total energy in ( 125) is equal to the potential energy.This gives A check of calculations is a sum which is equal to a constant energy value entering (125).This makes possible to calculate the oscillation period of time by an integration of the electron velocity obtained either from the formula (127), or from the energy balance represented by (129) and (131).The integration should be performed between two turning points, −a 0x and a 0x , of the oscillator.But an oscillator along the axis y can be equally considered.This changes the amplitude of the oscillation because of the asymmetry of the Hamiltonian in (123).For the sake of convenience this Hamiltonian can be written (k x = p y , k y = y, k z = 0): We have dy dt = ∂C sc(y) px orbits ∂p y = γ sin p y = 0 (133) attained at the amplitude a 0y different than a 0x because the right-hand side of (132) is different than that in (125).The kinetic energy of the motion is: and the acceleration in direction y: The corresponding potential energy is an integral of (135): with the integration constant obtained from the condition that at y = a 0y the potential energy (136) should be equal to the total energy in (132).A check of calculations is a sum which gives for any y the result of a constant energy entering (132).The cyclotron frequency can be calculated equally from the oscillator motions performed in x and y directions.The integral limits in the second case are (−a 0y , a 0y ).

Summary
The paper presents a method of calculating the cyclotron frequency of electrons circulating in cubic crystal lattices upon the action of a constant magnetic field.An approximation of the tightly-bound electrons is mainly assumed for the electron Hamiltonians.
In the first step, an analysis is developed on the symmetry properties of the crystal Hamiltonian as well as the similar properties of the Lorentz equation governing the electron motion in an external magnetic field.
For the field directed along the axis z, the planar motion in the space of the wave vector k has the variable k z as a constant parameter, and the variables k x and k y in the motion plane can be considered as canonical position and momentum variables of an anharmonic one-dimensional oscillator.The second step is based on a separation of the total energy of the oscillator into the kinetic and potential parts, which readily provides us with the analytic integral formulae for the time period of the cyclotron oscillation.This way is much more simple than a former one [11] in which the oscillation frequencies, or time periods, are obtained either from the action function, or from a direct solution of the differential equation of an oscillator.Moreover, the motion parameters of a crystal electron oscillating upon the action of the magnetic field, for example the electron velocity dx dt , can be easily examined for any electron position within the oscillation limits (22).This velocity can be obtained for any x on the basis of (26): on the condition that V (x) is a known function of x.In [11] the same velocity is calculated only with the aid of a rather tedious solution of the equation of motion.Moreover, the time t at any instant of the motion can be coupled with position x because of the formula: This dependence of t on x is absent in [11].The V (x) entering ( 139) and ( 140) is calculated for different lattices in (45), ( 52), ( 54) and ( 63), (64).
A similar examination can be done for the dependence of p x on x and t, on the condition that we note a coupling between x and p x due to the energy relations (23)-( 25).The same energy formulae can be applied in calculating the speed of the change of p x , namely cf. here (30).The dependence of the electron acceleration d 2 x dt 2 on x is easily obtainable from V (x) on the basis of (28).The dynamical parameters in more complicated band structures than listed in Section 4 (see Section 10) can be examined in a similar way.
The Hamiltonian plots of crystal electrons in the phase space (x, p x ) are also examined.Beyond the original Hamiltonian H(x, p x ), the Hamiltonians H(x, dx dt ) and H(p x , dpx dt ) are also considered.This kind of plots in the phase space is obtained when dx dt and dpx dt entering H are calculated from the Hamilton equations.
The calculations performed for the s-electrons taken as an example fit the cyclotron frequency data obtained on the basis of more conventional methods; see Table 1.Outlines of an extension of the method to other crystal cases than described by the tight-binding approximation and the s-electron band structure, for example the case of the p x -electron band, are also given.

Figure 1 .
Figure1.The (x, p x ) plot of a constant energy C sc of the simple cubic lattice done for two values of cos a 0 : 1/2 and −1/3.Figure1ais based on the Hamiltonian formula (23), Figure1bis obtained from the first step of (47), Figure1cis based on the energy expression given in (48a) and inferences below of it.x and p x are expressed in radians.

Figure 3 .
Figure 3.The (x, p x ) plot of a constant energy C fcc of the face centered cubic lattice done for two values of cos a 0 : 1/2 and −1/2.Figure 3(a) is based on the Hamiltonian formula (25); Figure3(b) is calculated on the basis of E kin obtained in the first step of (66) and V (x) taken from (67); Figure3(c) is a plot of the energy formula (37) specified for the case of the fcc lattice.

Figure 3 (
Figure 3.The (x, p x ) plot of a constant energy C fcc of the face centered cubic lattice done for two values of cos a 0 : 1/2 and −1/2.Figure 3(a) is based on the Hamiltonian formula (25); Figure3(b) is calculated on the basis of E kin obtained in the first step of (66) and V (x) taken from (67); Figure3(c) is a plot of the energy formula (37) specified for the case of the fcc lattice.