Abstract

Kubo formula was derived for the electron gas conductivity tensor on the nanotube surface in longitudinal magnetic field considering spatial and time dispersion. Components of the degenerate and nondegenerate electron gas conductivity tensor were calculated. The study has showed that under high electron density, the conductivity undergoes oscillations of de Haas-van Alphen and Aharonov-Bohm types with the density of electrons and magnetic field changes.

1. Introduction

Electromagnetic waves in the systems with cylindrical geometry have been studied for a long time [13]. The authors of the work [4] had solved the Maxwell’s equation taking into account the retardation for the conducting cylinder dipped in a dielectric medium. They considered time dispersion of the dielectric permittivity of conductor and dielectric. Also they had solved the dispersion equation for the wave spectrum and obtained new branches in the spectrum of surface polaritons. The paper [5] deals with polaritons in the magnetic wire with uniaxial anisotropy, dipped into a dielectric medium. Wire's dielectric and magnetic permittivity tensor frequency dependence also was studied. Polaritons spectra were found. Constant, parallel to the magnetic wire axis, external magnetic field effect was studied in paper [6]. In this paper, how wave characteristics allow one to obtain data about the material structure was shown.

A necessary of semiconductor nanotubes waveguide characteristics researching is caused by growing interest [710] to these systems. Plasma and magnetoplasma waves propagation on the nanotube surface were studied in [1118]. They focused on cylindrical geometry systems calculations. These calculations were made by means of matter dielectric permittivity simple models. The problem is that authors of [1118] have not used dielectric permittivity and conductivity tensors exact expressions. Also they have not taken into account the time and spatial dispersion. However, these exact formulas may significantly influence the tube waveguide characteristics expressions.

In this paper, we present conductivity tensor components calculation and their wave vector and frequency dependencies of the following system: electron gas on the nanotube surface, affected by a parallel to the tube axis external magnetic field. In Section 2, Kubo formula was obtained for the conductivity tensor. In Sections 3 and 4 we considered degenerate electron gas, and in Section 5, nondegenerate gas.

2. Dynamic Conductivity Tensor

Confined by the nanotube’s cylindrical surface, electron gas imbedded into a parallel to the cylinder axis magnetic field has axial symmetry. Therefore, it is suitable to characterize the electron state by the conserved quantities: angular momentum projection 𝑚=0,±1,, impulse projection 𝑘 of electron on tube axis 𝑧, and spin quantum number 𝜎=±1. Electron stationary state |𝑚𝑘𝜎 wave function and its energy spectrum are equal [1114] toΨ𝑚𝑘𝜎1(𝜑,𝑧,𝛼)=𝑆𝑒𝑖(𝑚𝜑+𝑘𝑧)𝜒𝜎𝜀(𝛼),(1)𝑚𝑘𝜎=𝜀0Φ𝑚+Φ02+𝑘22𝑚+𝜎𝜇𝐵𝐵,(2) where 𝑚 and 𝜇𝐵 are the effective mass and the electron spin magnetic moment, 𝜑 and 𝑧 are cylindrical coordinates, 𝜒𝜎(𝛼) is a spin wave function, 𝜀0=1/2𝑚𝑎2 is a rotational quantum, 𝑎 is the tube radius, Φ=𝜋𝑎2𝐵 is the magnetic induction 𝐵 flux through the tube section, Φ0=2𝜋𝑐/𝑒 is flux quantum, 𝑆=2𝜋𝑎𝐿 is the nanotube lateral surface area, and the length of nanotube is 𝐿. Hereinafter, the quantum constant is considered to equal unity. The spectrum (2) is a set of one-dimensional subzones, whose boundaries 𝜀𝑚0𝜎 are not equidistant.

One can expand the current electron density operator in field 𝐵 into series of cylindrical harmonics (1):𝐉(𝜑,𝑧)=𝑚𝑘𝐉(𝑚,𝑘)Ψ𝑚𝑘(𝜑,𝑧),(3) where 𝐉(𝑚,𝑘)=𝑒𝑚1𝑘1𝑚2𝑘2𝜎𝑚1𝑘1||||𝑚𝐕(𝑚,𝑘)2𝑘2𝑎+𝑚1𝑘1𝜎𝑎𝑚2𝑘2𝜎.(4) Here 𝑎𝑚𝑘𝜎 and 𝑎+𝑚𝑘𝜎 are annihilation and creation operators of electrons in state |𝑚𝑘𝜎,and 𝑚1𝑘1|𝐕(𝑚,𝑘)|𝑚2𝑘2 are operator 𝐕 matrix elements 1𝐕(𝑚,𝑘)=2𝑚𝑒𝑖𝑐𝐀Ψ𝑚𝑘(𝜑,𝑧)+Ψ𝑚𝑘𝑒(𝜑,𝑧)×𝑖𝑐𝐀(5) in the functional basis (1), 𝐁=rot𝐀. The operator’s (4) cylindrical components are equal to𝐽𝜑𝑒(𝑚,𝑘)=𝑚𝑎𝑆𝑚𝑘𝜎𝑚+ΦΦ0+𝑚2𝑎+𝑚𝑘𝜎𝑎(𝑚+𝑚)(𝑘+𝑘)𝜎,𝐽𝑧𝑒(𝑚,𝑘)=𝑚𝑆𝑚𝑘𝜎𝑘+𝑘2𝑎+𝑚𝑘𝜎𝑎(𝑚+𝑚)(𝑘+𝑘)𝜎.(6)

Within the framework of electron gas linear response theory [19, 20] to the weak electric field 𝐄, changes with wave vector 𝑞 and frequency 𝜔, it is possible to obtain the current density𝑗𝜇(𝑚,𝑞,𝜔)=𝜈𝜎𝜇𝜈(𝑚,𝑞,𝜔)𝐸𝜈(𝑚,𝑞,𝜔),(7) where𝜎𝜇𝜈𝑒(𝑚,𝑞,𝜔)=𝑖2𝑛𝑚𝜔𝛿𝜇𝜈+1𝜔0𝑑𝑡𝑒𝑖𝜔𝑡×𝐽𝜇(𝑚,𝑞,𝑡),𝐽𝜈,(𝑚,𝑞,0)(8) Kubo formula for conductivity tensor. Here 𝐽𝜇(𝑚,𝑞,𝑡) is a component of Heisenberg operator (6), 𝑛 is a surface electron density, [𝑎,𝑏]=𝑎𝑏𝑏𝑎, and angle brackets mark Gibbs average.

We have substituted (6) in (8) and obtained the relation between the conductivity tensor and Fourier time component of two-electron Greens function [19, 20]:𝐾𝑡𝑎(1,2;3,4)=𝑖𝜃(𝑡)+1(𝑡)𝑎2(𝑡),𝑎+3(0)𝑎4,(0)(9) where 1=(𝑚1,𝑘1,𝜎1),, 𝜃(𝑡) is a Heaviside’s function. In case of free electron gas, function (9) is expressed by one-electron Green function. In the result, the conductivity tensor components (8) are equal to𝜎𝜑𝜑𝑒(𝑚,𝑞,𝜔)=𝑖2𝑛𝑚𝜔+𝑖𝑒2𝑚2𝑎2𝜔𝑆𝑚𝑘𝜎𝑓𝜀𝑚𝑘𝜎𝑚+Φ/Φ0+𝑚/22𝜀𝑚𝑘𝜎𝜀(𝑚+𝑚)(𝑘+𝑞)𝜎𝑚+𝜔+𝑖0+Φ/Φ0𝑚/22𝜀(𝑚𝑚)(𝑘𝑞)𝜎𝜀𝑚𝑘𝜎,𝜎+𝜔+𝑖0𝜑𝑧(𝑚,𝑞,𝜔)=𝜎𝜑𝑧=(𝑚,𝑞,𝜔)𝑖𝑒2𝑚2𝑎𝜔𝑆𝑚𝑘𝜎𝑓𝜀𝑚𝑘𝜎𝑚+Φ/Φ0+𝑚/2(𝑘+𝑞/2)𝜀𝑚𝑘𝜎𝜀(𝑚+𝑚)(𝑘+𝑞)𝜎𝑚+𝜔+𝑖0+Φ/Φ0𝑚/2(𝑘𝑞/2)𝜀(𝑚𝑚)(𝑘𝑞)𝜎𝜀𝑚𝑘𝜎,𝜎+𝜔+𝑖0𝑧𝑧𝑒(𝑚,𝑞,𝜔)=𝑖2𝑛𝑚𝜔+𝑖𝑒2𝑚2𝜔𝑆𝑚𝑘𝜎𝑓𝜀𝑚𝑘𝜎(𝑘+𝑞/2)2𝜀𝑚𝑘𝜎𝜀(𝑚+𝑚)(𝑘+𝑞)𝜎+𝜔+𝑖0(𝑘𝑞/2)2𝜀(𝑚𝑚)(𝑘𝑞)𝜎𝜀𝑚𝑘𝜎,+𝜔+𝑖0(10) where 𝑓(𝜀) is Fermi function. These components are included into the dispersion equation for electromagnetic waves, propagating along the tube.

In case of 𝑎, one can obtain from (10) the well-known Drude-Lorentz expression, which describes two-dimensional electron gas with the magnetic field absence: 𝜎𝜑𝜑=𝜎𝑧𝑧=𝑖𝑒2𝑛/𝑚𝜔, 𝜎𝜑𝑧=0. The real part of conductivity tensor in this case is equal to zero. It is caused by not taking into account the electrons collisions. That effect can be considered by means of 𝜔𝜔+𝑖𝜈𝑐 replacement in (10). Here 𝜈𝑐 is the collisions frequency.

3. Degenerate Gas

At zero temperature, integrals over 𝑘, included into (10), can be calculated exactly. Thus, real and imaginary parts of conductivity (10) are equal to Re𝜎𝜑𝜑𝑒(𝑚,𝑞,𝜔)=24𝜋𝑚||𝑞||𝑎3𝜔×𝑚𝜎𝑚+ΦΦ0+𝑚22×𝜃𝐴++𝐴𝜃+𝑚+ΦΦ0𝑚22×𝜃𝐴++𝐴𝜃+,Re𝜎𝜑𝑧𝑒(𝑚,𝑞,𝜔)=24𝜋(𝑞𝑎)2𝜔×𝑚𝜎𝑚+ΦΦ0+𝑚2𝜔Ω+×𝜃𝐴++𝐴𝜃+Φ𝑚+Φ0𝑚2𝜔Ω×𝜃𝐴++𝐴𝜃+,Re𝜎𝑧𝑧(𝑚𝑚,𝑞,𝜔)=𝑒2||𝑞||4𝜋3×𝑎𝜔𝑚𝜎𝜔Ω+2𝜃𝐴++𝐴𝜃+𝜔Ω2𝜃𝐴++𝐴𝜃+,(11)Im𝜎𝜑𝜑𝑒(𝑚,𝑞,𝜔)=2𝑛𝑚𝜔𝑒24𝜋2𝑚𝑞𝑎3𝜔×𝑚𝜎Φ𝑚+Φ0+𝑚22||||𝐴ln++𝐴+||||Φ𝑚+Φ0𝑚22||||𝐴ln++𝐴+||||,Im𝜎𝜑𝑧𝑒(𝑚,𝑞,𝜔)=2𝑚𝑛𝑚𝑒𝑞𝑎𝜔24𝜋2(𝑞𝑎)2𝜔×𝑚𝜎Φ𝑚+Φ0+𝑚2𝜔Ω+||||𝐴×ln++𝐴+||||Φ𝑚+Φ0𝑚2×𝜔Ω||||𝐴ln++𝐴+||||,Im𝜎𝑧𝑧(𝑚,𝑞,𝜔)=2𝑒2𝜀0𝑚2𝑛𝑞2𝜔𝑚𝑒24𝜋2𝑞3×𝑎𝜔𝑚𝜎𝜔Ω+2||||𝐴ln++𝐴+||||𝜔Ω2||||𝐴ln++𝐴+||||,(12) where𝐴++=𝑞𝜐𝜎𝑚𝜔+Ω+,𝐴+=𝑞𝜐𝜎𝑚𝜔+Ω+,𝐴++=𝑞𝜐𝜎𝑚𝜔++Ω,𝐴+=𝑞𝜐𝜎𝑚𝜔++Ω,𝜐𝜎𝑚=2𝑚𝜀𝐹𝜀𝑚𝜎1/2(13) is the electron motion along the tube axis velocity maximum in subzone (𝑚,𝜎), 𝜀𝑚𝜎=𝜀𝑚0𝜎 is the subzone boundary, 𝜀𝐹 is the Fermi energy, 𝜔±=𝜔±𝑞2/2𝑚, andΩ±(𝑚,𝑚)=𝜀0±Φ𝑚+Φ0±𝑚2Φ𝑚+Φ02(14) is the frequency of the “vertical” electron transitions between the spectrum (2) subzones. While obtaining (12) out of (10), the identity below was used 𝐵2𝐶𝐵𝐶=𝐵+𝐶+2.𝐵𝐶(15)

4. Longitudinal Conductivity

This section concentrates on the analysis of degenerate electron gas longitudinal conductivity, (11) (12), assuming 𝑚=0. The numerical calculations were carried out also for 𝑚=1 case. In this case, electrons make only intra-subzone transitions without spin change. Longitudinal conductivity is included into the propagating along the nanotube axis intra-subzone magnetoplasma waves dispersion equation [13], because it is related to the polarization operator 𝑃 by the relation belowIm𝜎𝑧𝑧𝑒(0,𝑞,𝜔)=2𝜔𝑞2Re𝑃(0,𝑞,𝜔).(16) In case of 𝑚=0, we obtain the following equations from (11) and (12): Re𝜎𝑧𝑧𝑒(𝑞,𝜔)=2𝑚𝜔||𝑞||4𝜋3𝑎×𝑚𝜎𝜃𝑞𝜐𝜎𝑚𝜔𝜃𝑞𝜐𝜎𝑚𝜔𝜃𝑞𝜐𝜎𝑚𝜔++𝜃𝑞𝜐𝜎𝑚𝜔+,Im𝜎𝑧𝑧𝑒(𝑞,𝜔)=2𝑚𝜔4𝜋2𝑞3𝑎𝑚𝜎||||ln𝑞𝜐𝜎𝑚𝜔𝑞𝜐𝜎𝑚𝜔||||||||ln𝑞𝜐𝜎𝑚𝜔+𝑞𝜐𝜎𝑚𝜔+||||.(17) Suppose that the flux ratio in (2) should be equal to Φ/Φ0=𝑀+𝜂, where 𝑀=0,1, is an integer part of Φ/Φ0and 𝜂 is a fractional part (0𝜂<1). One can use the achievable in experiments values of carbon and semiconductor nanotube radii and obtain that the value of 𝑀 in the accessible magnetic fields is low, unless the tube radius is too large. Then the energy of electron (2) is minimal in the subzone (𝑚=𝑀, 𝜎=1):𝜀min=𝜀0𝜂2𝜇𝐵𝐵.(18) In this case, the electron rotation energy is nearly compensated by the rotation in the magnetic field. Usually spin level splitting (2) is small, that is,2𝜇𝐵𝐵<𝜀0(1+2𝜂),(19) and the electron density corresponds to inequality𝑛<𝑚𝜇𝐵𝐵𝜋2𝑎.(20) Then electrons partially fill only the lower subzone (𝑀,1). Therefore, (17) leads to the following:𝜋Re𝛿𝜎=2𝑘0𝑞𝑥2,𝑘Im𝛿𝜎=0𝑥2|||||2𝑞ln1+𝑞/2𝑘02𝑥21𝑞/2𝑘02𝑥2|||||,(21) where𝛿𝜎=𝜎𝑧𝑧𝑚(𝑞,𝜔)𝜔𝑒2𝑛𝜔,𝑥=𝑞𝜐0,𝜐0=𝑘0𝑚=2𝑚𝜀𝐹+𝜇𝐵𝐵𝜀0𝜂21/2(22) is the electrons velocity maximum in subzone (𝑀,1). Figures 1 and 2 show the 𝑥=𝜔/𝑞𝜐0 dependencies of Re𝛿𝜎 and Im𝛿𝜎 under 𝜂=0.01,𝑞=𝑘0,and 𝑘0𝑎=1 in quantum limit, when Fermi energy is located in subzone (𝑀,1).

One can calculate what included into (17) sum by 𝑚 using a Poisson formula, assuming 𝜀𝐹𝜀0, that is, a large quantity of subzones are filled. Within long wavelength assumption 𝑞𝜐𝐹𝜔 (𝜐𝐹 is a Fermi velocity), we obtain the following formula from (17): Im𝜎𝑧𝑧=𝑒2𝑛𝑚𝜔21+𝜋2𝜀0𝜀𝐹3/4𝑙=11𝑙3/2×sin2𝜋𝑙𝜀𝐹𝜀0𝜋4Φcos2𝜋𝑙Φ0.(23) Here we neglected the spin level splitting (2). Monotonous part of conductivity (23) under the given electron density does not depend on the tube radius and coincides with the two-dimensional electron gas conductivity [21]. It is within the 𝑎, when the tube is cut along its generatrix and unrolled onto a plane. Conductivity (23) experiences the Aharonov-Bohm and de Haas-van Alphen oscillation types. The reason of the de Haas-van Alphen oscillation types existance is that the electron states density has root peculiarities on the subzone bounaries. Change of the tube radius or electron density causes these peculiarities to pass throw Fermi energy and de Haas-van Alphen oscillation types to appear. If we study conductivity dependence (23) on 𝜀𝐹1/2=(𝑛/𝜈)1/2 (𝜈 is the two-dimensional electron gas states density), then the oscillation period is 1/2𝑚𝑎. If we study the radius dependence, the period of conductivity oscillations is equal to 1/2𝑚𝜀𝐹. Aharonov-Bohm oscillations are caused by the magnetic flux through the tube cross-section change. Their period is equal to the flux quantum. Figure 3 shows the dependence of the (23) first conductivity harmonic on (𝜀𝐹/𝜀0)1/2 under Φ/Φ0=0.01.

5. Nondegenerate Electron Gas

Nondegenerate electron gas emerges if inequality 𝛽|𝜇|1 is true on the semiconductor nanotube surface. Here 𝛽 is a reverse temperature and 𝜇 is a chemical potential. Using Boltzmann’s statistics, the chemical potential is related to the electron density as 1𝑛=𝜋3/2𝑎𝑚𝑒2𝛽𝛽𝜇ch𝛽𝜇𝐵𝐵𝑚=𝑒𝛽𝜀𝑚,(24) where 𝜀𝑚=𝜀0(𝑚+Φ/Φ0)2. The included here summation may be transformed by means of formula [22] 𝑚=𝑒𝑥(𝑚+𝜐)2=𝜋𝑥𝑙=𝑒𝜋2𝑙2/𝑥cos(2𝜋𝑙𝜐),𝑥>0.(25) Real and imaginary tensor (10) parts in the studied case assuming 𝑚=0 are equal toRe𝜎𝜑𝜑=𝑒(𝑞,𝜔)24𝜋𝑚𝛽𝜀03/2||𝑞||𝑎3𝜔𝑒𝛽𝜇ch𝛽𝜇𝐵𝐵×𝑚exp𝛽2𝜔𝑞2𝑚exp𝛽2𝜔+𝑞2×𝑙=𝜋exp2𝑙2𝛽𝜀012𝜋2𝑙2𝛽𝜀0Φcos2𝜋𝑙Φ0,Re𝜎𝜑𝑧=(𝑞,𝜔)𝜋𝑒22𝛽𝜀03/2(𝑞𝑎)2𝑒𝛽𝜇ch𝛽𝜇𝐵𝐵×𝑙=𝜋𝑙exp2𝑙2𝛽𝜀0×𝑚exp𝛽2𝜔𝑞2𝑚exp𝛽2𝜔+𝑞2Φ×sin2𝜋𝑙Φ0,Re𝜎𝑧𝑧=𝑒(𝑞,𝜔)2𝑛𝜔||𝑞||3𝜋𝑚𝛽2×𝑚exp𝛽2𝜔𝑞2𝑚exp𝛽2𝜔+𝑞2,Im𝜎𝜑𝜑=𝑒(𝑞,𝜔)2𝑛𝑚𝜔+𝑒24𝜋𝑚𝛽𝜀03/2𝑞𝑎3𝜔𝑒𝛽𝜇ch𝛽𝜇𝐵𝐵×𝐹𝑚𝛽2𝜔𝑞𝐹𝑚𝛽2𝜔+𝑞×𝑙=𝜋exp2𝑙2𝛽𝜀012𝜋2𝑙2𝛽𝜀0Φcos2𝜋𝑙Φ0,Im𝜎𝜑𝑧=𝑒(𝑞,𝜔)22𝛽𝜀03/2(𝑎𝑞)2𝑒𝛽𝜇ch𝛽𝜇𝐵𝐵×𝑙=𝜋𝑙exp2𝑙2𝛽𝜀0×𝐹𝑚𝛽2𝜔𝑞𝐹𝑚𝛽2𝜔+𝑞Φ×sin2𝜋𝑙Φ0,Im𝜎𝑧𝑧=𝑒(𝑞,𝜔)2𝑛𝜔𝑞3𝑚𝛽2×𝐹𝑚𝛽2𝜔𝑞𝐹𝑚𝛽2𝜔+𝑞,(26) where1𝐹(𝑥)=𝜋𝑝𝑒𝑑𝑦𝑦2.𝑥𝑦(27) Components 𝜎𝜑𝜑 and 𝜎𝜑𝑧 in (26) are expressed through 𝛽 and 𝜇, that is, they refer to an open system of tube surface electrons. They experience only Aharonov-Bohm oscillations with the magnetic field changes. Aharonov-Bohm oscillations are not taking place in 𝜎𝑧𝑧 component while it is described in 𝛽 and 𝑛 terms.

6. Conclusion

Simplified conductivity models are usually used for studying electromagnetic waves propagation in the cylindrical geometry systems, for example, in nanotubes. The metal cylinder conductivity is often believed to be endless, and the dielectric permittivity of the matter, where cylinder is dipped, is considered to be constant or only frequency dependent. Conductivity’s spatial dispersion as a rule is not taken into account. Nonetheless, the electromagnetic field’s nature in the tube and its waveguide characteristics are sensitive to the surface currents. Therefore, the electron gas conductivity tensor components calculation problem, taking into account spatial and time dispersion, is worth consideration. It is given in this work conductivity tensor components may be used in obtaining a dispersion equation for electromagnetic waves in the tube. The example of longitudinal conductivity showed what data of electron gas is possible to obtain by measuring conductivity. Particularly, imaginary part of conductivity experiences de Haas-van Alphen and Aharonov-Bohm oscillations types with the electron density and magnetic field changes. The oscillation periods measurement enables one to identify the electron effective mass and the combination of universal constants included into flux quantum.

Acknowledgment

The authors would like to thank T. Rashba for help with preparation of the paper.