Landauer-Datta-Lundstrom Generalized Transport Model for Nanoelectronics

The Landauer-Datta-Lundstrom electron transport model is briefly summarized. If a band structure is given, the number of conduction modes can be evaluated and if a model for a mean-free-path for backscattering can be established, then the nearequilibrium thermoelectric transport coefficients can be calculated using the final expressions listed below for 1D, 2D, and 3D resistors in ballistic, quasiballistic, and diffusive linear response regimes when there are differences in both voltage and temperature across the device. The final expressions of thermoelectric transport coefficients through the Fermi-Dirac integrals are collected for 1D, 2D, and 3D semiconductors with parabolic band structure and for 2D graphene linear dispersion in ballistic and diffusive regimes with the power law scattering.


Introduction
The objectives of this short review is to give a condensed summary of Landauer-Datta-Lundstrom (LDL) electron transport model [1][2][3][4][5] which works at the nanoscale as well as at the macroscale for 1D, 2D, and 3D resistors in ballistic, quasiballistic, and diffusive linear response regimes when there are differences in both voltage and temperature across the device.
Appendices list final expressions of thermoelectric transport coefficients through the Fermi-Dirac integrals for 1D, 2D, and 3D semiconductors with parabolic band structure and for 2D graphene linear dispersion in ballistic and diffusive regimes with the power law scattering.

Generalized Model for Current
The generalized model for current can be written in two equivalent forms: where "broadening" () relates to transit time for electrons to cross the resistor channel: density of states () with the spin degeneracy factor   = 2 is included; () is the integer number of modes of conductivity at energy ; the transmission where () is the mean-free-path for backscattering and  is the length of the conductor; Fermi function is indexed with the resistor contact numbers 1 and 2;   is the Fermi energy which as well as temperature  may be different at both contacts.
Journal of Nanoscience Equation (3) can be derived with relatively few assumptions and it is valid not only in the ballistic and diffusion limits, but in between as well: Diffusive:  ≫ ;  =   ≪ 1, Ballistic:  ≪ ;  → 1, Quasi-ballistic:  ≈ ;  < 1. ( The LDL transport model can be used to describe all three regimes.
It is now clearly established that the resistance of a ballistic conductor can be written in the form where ℎ/ 2 is fundamental Klitzing constant and number of modes () represents the number of effective parallel channels available for conduction.This result is now fairly well known, but the common belief is that it applies only to short resistors and belongs to a course on special topics like mesoscopic physics or nanoelectronics.What is not well known is that the resistance for both long and short conductors can be written in the form ) .
Ballistic and diffusive conductors are not two different worlds, but rather a continuum as the length  is increasing.Ballistic limit is obvious for  ≪ , while for  ≫  it reduces into standard Ohm's law: Indeed we could rewrite () above as with a new expression for specific resistivity: which provides a different view of resistivity in terms of the number of modes per unit area and the mean-free-path.Number of modes is proportional to the width  of the resistor in 2D and to the cross-sectional area  in 3D; ⟨V +  ()⟩ is the average velocity in the + direction from contact 1 to contact 2.
where  is the area of the 2D resistor, Ω is the volume of the 3D resistor, and ( −   ) is the Heaviside step function.
Then number of modes is where  V is the valley degeneracy.
Figure 1 shows qualitative behavior of the density of states and number of modes for resistors with parabolic band structure.
For linear dispersion in graphene, where +sign corresponds to conductivity band with   > 0 (-type graphene) and −sign corresponds to valence band with   < 0 (-type graphene), Density of states in graphene is and number of modes is Two equivalent expressions for specific conductivity deserve attention, one as a product of () and the diffusion coefficient (): where with () being the mean free time after which an electron gets scattered, and the other as a product of () and (): where the three items in parenthesis correspond to 1D, 2D, and 3D resistors.Although (19b) is not well known, the equivalent version in (19a) is a standard result that is derived in textbooks.Both (19a) and (19b) are far more generally applicable compared with traditional Drude model.For example, these equations give sensible answers even for materials like graphene whose nonparabolic bands make the meaning of electron mass somewhat unclear, causing considerable confusion when using Drude model.In general we must really use (19a) and (19b) and not Drude model to shape our thinking about conductivity.
These conceptual equations are generally applicable even to amorphous materials and molecular resistors.Irrespective of the specific () relation at any energy, the density of states (), velocity ](), and momentum () are related to the total number of states () with energy less than  by the fundamental relation where  is the number of dimensions.Being combined with (19a), it gives one more fundamental equation for conductivity: where electron mass is defined as For parabolic () relations, the mass is independent of energy, but in general it could be energy-dependent as, for example, in graphene the effective mass where  0 () is the equilibrium Fermi function, and an applied bias is small enough.Using Taylor expansion under constant temperature condition and property of the Fermi function one finds The derivative of the Fermi function multiplied by  to make it dimensionless is known as thermal broadening function and shown in Figure 2. If one integrates   over all energy range, the total area so that we can approximately visualize   as a rectangular pulse centered around  =   with a peak value of 1/4 and a width of ∼4.
The derivative (− 0 /) is known as the Fermi conduction window function.Whether a conductor is good or bad is determined by the availability of the conductor energy states in an energy window ∼±2 around the electrochemical potential  0 , which can vary widely from one material to another.Current is driven by the difference  1 −  2 in the "agenda" of the two contacts which for low bias is proportional to the derivative of the equilibrium Fermi function (26).With this near-equilibrium assumption for current (1b), we have with conductivity known as the Landauer expression which is valid in 1D, 2D, and 3D resistors, if we use the appropriate expressions for ().
For ballistic limit () = 1.For diffusive transport () is given by (3).For a conductor much longer than a meanfree-path the current density equation for diffusive transport is where the electrochemical potential   is also known as the quasi-Fermi level.For a 2D conductor the surface specific conductivity is or in a different form where differential specific conductivity Similar expressions can be written for 1D and 3D resistors.Another way to write the conductance is the product of the quantum of conductance, times the average transmission, times the number of modes in the Fermi windows: Yet another way to write the conductance is in terms of the differential conductance   () as It is important from practical point of view that both expressions-for the electric current (1a) and (1b) and thermocurrent (36)-are suitable for analysis of conductivity of any materials from metals to semiconductors up to modern nanocomposites.
When there are differences in both voltage and temperature across the resistor, then we must the Fermi difference ( 1 −  2 ) expands to Taylor series in both voltage and temperature and get where Deriving a general near-equilibrium current equation is now straightforward.The total current is the sum of the contributions from each energy mode: where the differential current is Using (37) we obtain where is the differential conductance and is the Soret coefficient for electrothermal diffusion in differential form.Note that    () is negative for modes with energy above   (-resistors) and positive for modes with energy below   (-resistors).Now we integrate (39a) over all energy modes and find with three transport coefficients, namely: conductivity given by (35a) and (35b); the Soret electrothermal diffusion coefficient and the electronic heat conductance under the short circuit conditions (Δ = 0) where current  is defined to be positive when it flows in conductor from contact 2 to contact 1 with electrons flowing in opposite direction.The heat current   is positive when it flows in the + direction out of contact 2. Equations (40a), (40b), (40c), and (40d) for long diffusive resistors can be written in the common form used to describe bulk transport as with three specific transport coefficients These equations have the same form for 1D and 2D resistors, but the units of the various terms differ.
The inverted form of (40a), (40b), (40c), and (40d) is often preferred in practice, namely: where In this form of the equations, the contributions from each energy mode are not added; for example,  ̸ = ∫ ().Similarly, the inverted form of the bulk transport equations (41a), (41b), (41c), (41d), and (41e) becomes with transport coefficients In summary, when a band structure is given, number of modes can be evaluated from (11a), (11b), and (11c) and if a model for the mean-free-path for backscattering () can be chosen, then the near-equilibrium transport coefficients can be evaluated using the expressions listed above.

Bipolar Conduction.
Let us consider a 3D semiconductor with parabolic dispersion.For the conduction band and for the valence band The conductivity is provided with two contributions: for the conduction band and for the valence band The Seebeck coefficient for electrons in the conduction band follows from (41a), (41b), (41c), (41d), and (41e): Similarly, for electrons in the valence band we have but the sign of   will be positive.What is going on when both the conduction and valence bands contribute to conduction?This can happen for narrow bandgap conductors or at high temperatures.In such a case, we have to simply integrate over all the modes and will find moreover, we do not have to be worried about integrating to the top of the conduction band or from the bottom of the valence band because the Fermi function ensures that the integrand falls exponentially to zero away from the band edge.What is important is that in both cases we integrate the same expression with the appropriate  3 () and () over the relevant energy difference  2 −  1 .Electrons carry current in both bands.Our general expression is the same for the conduction and valence bands.There is no need to change signs for the valence band or to replace  0 () with 1 −  0 ().
To calculate the Seebeck coefficient when both bands contribute let us be reminded that in the first direct form of the transport coefficients (41a), (41b), (41c), (41d), and (41e) the contributions from each mode are added in parallel so the total specific Soret coefficient then, the Seebeck coefficient for bipolar conduction Since the Seebeck coefficients for the conduction and valence bands have opposite signs, the total Seebeck coefficient just drops for high temperatures and the performance of a thermoelectrical device falls down.In summary, given a band structure dispersion, the number of modes can be evaluated and if a model for a mean-free-path for backscattering can be established, then the near-equilibrium transport coefficients can be calculated using final expressions listed above.

Heat Transfer by Phonons
Electrons transfer both charge and heat.Electrons carry most of the heat in metals.In semiconductors electrons carry only a part of the heat but most of the heat is carried by phonons.
The phonon heat flux is proportional to the temperature gradient with coefficient   known as the specific lattice thermal conductivity.Such an exceptional thermal conductor like diamond has   ≈ 2 ⋅ 10 3 W/m ⋅ K while such a poor thermal conductor like glass has   ≈ 1 W/m ⋅ K.Note that electrical conductivities of solids vary over more than 20 orders of magnitude, but thermal conductivities of solids vary over a range of only 3-4 orders of magnitude.We will see that the same methodology used to describe electron transport can be also used for phonon transport.We will also discuss the differences between electron and phonon transport.For a thorough introduction to phonons use classical books [6][7][8][9].
To describe the phonon current we need an expression like for the electron current (1b) written now as For electrons the states in the contacts were filled according to the equilibrium Fermi functions, but phonons obey Bose statistics; thus, the phonon states in the contacts are filled according to the equilibrium Bose-Einstein distribution Let temperature for the left and the right contacts be  1 and  2 .As for the electrons, both contacts are assumed ideal.Thus the phonons that enter a contact are not able to reflect back, and transmission coefficient  ph () describes the phonon transmission across the entire channel.
It is easy now to rewrite (51) to the phonon heat current.Electron energy  we replace by the phonon energy ℏ.In the electron current we have charge  moving in the channel; in case of the phonon current the quantum of energy ℏ is moving instead; thus, we replace  in (51) with ℏ and move it inside the integral.The coefficient 2 in (51) reflects the spin degeneracy of an electron.In case of the phonons we remove this coefficient, and instead the number of the phonon polarization states that contribute to the heat flow lets us include to the number of the phonon modes  ph (ℏ).Finally, the heat current due to phonons is In the linear response regime by analogy with (26), where the derivative according to (52) is with Now (53) for small differences in temperature becomes where the thermal conductance Equation ( 57) is simply Fourier's law stating that heat flows down to a temperature gradient.It is also useful to note that the thermal conductance (58) displays certain similarities with the electrical conductance The derivative known as the Fermi window function that picks out those conduction modes which only contribute to the electric current.The electron windows function is normalized: In case of phonons the term in square brackets of (58) acts as a window function to specify which modes carry the heat current.After normalization, thus finally with known as the quantum of thermal conductance experimentally observed first in 2000 [10].
Comparing (59) and ( 63) one can see that the electrical and thermal conductances are similar in structure: both are proportional to corresponding quantum of conductance times an integral over the transmission times the number of modes times a window function.
The thermal broadening functions for electrons and phonons have similar shapes and each has a width of a few .In case of electrons this function is given by (27) or with  ≡ (−  )/ and is shown on Figure 2.This function for phonons is given by (62) or with  ≡ ℏ/.Both functions are normalized to a unity and shown together on Figure 3.
Along with the number of modes determined by the dispersion relation, these two window functions play a key role in determining the electrical and thermal conductances.

Thermal Conductivity of the Bulk Conductors.
The thermal conductivity of a large diffusive resistor is a key material property that controls performance of any electronic devices.By analogy with the transmission coefficient (3) for electron transport, the phonon transmission It is also obvious that for large 3D conductors the number of phonon modes is proportional to the cross-sectional area of the sample:  ph (ℏ) ∝ . (68) with specific lattice thermal conductivity or substituting (67) into (63), one for the lattice thermal conductivity finally obtains It is useful now to define the average number of phonon modes per cross-sectional area of the conductor that participate in the heat transport: Then where the average mean-free-path is defined now as Thus, the couple of the phonon transport equations ( 69) and (73) corresponds to similar electron transport equations: The thermal conductivity (73) and the electrical conductivity (76) have the same structure.It is always a product of the corresponding quantum of conductance times the number of modes that participate in transport, times the average mean-free-path.These three quantities for phonons will be discussed later.

Specific Heat versus Thermal Conductivity.
The connection between the lattice specific thermal conductivity and the lattice specific heat at constant volume is well known [6][7][8][9].We will show now that corresponding proportionality coefficient is a product of an appropriately-defined meanfree-path ⟨⟨Λ ph ⟩⟩ and an average phonon velocity ⟨V ph ⟩, namely: The total phonon energy per unit volume where  ph (ℏ) is the phonon density of states.By definition, where ( 55) and (62) were used.Next, multiply and divide (71) by ( 79) and obtain the proportionality we are looking for: To obtain final expression (77) and correct interpretation of the proportionality coefficient we need to return to (67).This expression can be easily derived for 1D conductor with several simplifying assumptions.Nevertheless it works very well in practice for a conductor of any dimension.Derivation of (67) is based on the interpretation of the mean-free-path () or (ℏ) as that its inverse value is the probability per unit length that a positive flux is converted to a negative flux.This is why  is often called a mean-freepath for backscattering.Let us relate it to the scattering time .The distinction between mean-free-path and mean-freepath for backscattering is easiest to see for 1D conductor.Let an electron undergo a scattering event.For isotropic scattering the electron can forward scatter or backscatter.Only backscattering is relevant for the mean-free-path for scattering, so the time between backscattering events is 2.Thus the mean-free-path for backscattering is twice the mean-free-path for scattering: It was shown that the proper definition of the mean-freepath for backscattering for a conductor of any dimension [11] is where averaging is performed over angles.For isotropic bands, The scattering time is often approximately written as the power law scattering: where exponent  describes the specific scattering mechanism: for acoustic phonon scattering in 3D conductor with parabolic dispersion  = −1/2 and for ionized impurity scattering  = +3/2 [12].Analogous power law is often used for mean-free-path: For parabolic zone structure V() ∝  1/2 ; thus  =  + 1/2 with  = 0 for acoustic phonon scattering, and  = 2 for ionized impurity scattering.
Coming back to our initial task to derive (77) from (80) for 3D conductor according to (81c) we have where, according to (81a), and finally It was stated above in (11c) that the density of states and number of modes for electrons in 3D are Let us rewrite this formula for phonons.Note that the spin degeneracy for electrons   = 2 is included with the density of states: and for spherical bands in 3D conductor, Collecting (87) up to (89) all together in case of phonons, we have Substituting ( 86) and ( 90) into (80), we obtain Multiplying and dividing (91) by we finally get (77) with proportionality coefficient between   and   as the product of an average mean-free-path as and an average velocity as with the appropriate averaging.Equation ( 77) is often used to estimate the average meanfree-path from the measured   and   , if we know the average velocity, which is frequently assumed to be the longitudinal sound velocity.The derivation above has identified the precise definitions of the ⟨⟨Λ ph ⟩⟩ and ⟨V ph ⟩.If a phonon dispersion is chosen one can always compute the average velocity according to (94), and it is typically very different from the longitudinal sound velocity.Thus, estimates of the average mean-free-path can be quite wrong if one assumes the longitudinal sound velocity [13].

Debye Model.
For 3D conductors there are three polarization states for lattice vibrations: one for atoms displaced in the direction of propagation (longitudinal/) and two for atoms displaced orthogonally to the direction of propagation (transverse/).The low energy modes are called acoustic modes/: one  mode analogous to sound waves propagating in air and two  modes.Near  → 0 dispersion of acoustic modes is linear: and is known as Debye approximation.The Debye velocity V  is an average velocity of the  and  acoustic modes.In case of  mode V  is simply the sound velocity V  ∝  −1/2 with  being an effective mass of vibrating atom.Typically, V  ≈ 5 × 10 3 m/s, about 20 times slower that the velocity of a typical electron.The bandwidth of the electronic dispersion is typically BW ≫ , so only states near the bottom of the conduction band where the effective mass model is reasonably accurate are occupied.For phonons the situation is much different; the bandwidth BW ≈ , so states across the entire Brillouin zone are occupied.The widely used Debye approximation (95) fits the acoustic branches as long as  is not too far from the center of the Brillouin zone.
With the Debye approximation (95) it is easy to find the density of the phonon states: where the factor of three is for the three polarizations.Then one can obtain the number of phonon modes per crosssectional area from (90): Since all the states in the Brillouin zone tend to be occupied at moderate temperatures, we are to be sure that we account for the correct number of states.For a crystal there are 3/Ω states per unit volume.To find the total number of states we have to integrate the density of states: with the upper limit as ℏ times the so-called Debye frequency to produce the correct number of states, namely: The Debye frequency defines a cutoff frequency above which no states are accounted for.This restriction can also be expressed via a cutoff wave vector   or as a Debye temperature: For  ≪   , only states with  → 0 for which the Debye approximation is accurate are occupied.Now we can calculate the lattice thermal conductivity by integrating (71) to the Debye cutoff energy and estimate  ph (ℏ) according to (71).The integral can be taken numerically or analytically if appropriate expression for the mean-free-path is used.This is how the lattice thermal conductivities were first calculated [14,15].The theory and computational procedures for the thermoelectric transport coefficients were developed further in [11,13,16].

Phonon Scattering.
Phonons can scatter from defects, impurity atoms, isotopes, surfaces and boundaries, and electrons and from other phonons.Phonon-phonon scattering occurs because the potential energy of the bonds in the crystal is not exactly harmonic.All higher order terms are treated as a scattering potential.Two types of phonon scattering are considered.In the normal process two phonons interact and create a third phonon with energy and momentum being conserved: The total momentum of the phonon ensemble is conserved; thus this type of scattering has little effect on the heat flux.
In a second type of scattering, umklapp/-scattering, the two initial phonons have larger momentum; thus the resulting phonon would have a momentum outside the Brillouin zone due to unharmonic phonon-phonon as well as electron-phonon interactions.The -scatterings are the basic processes in the heat transport especially at high temperatures.Scattering on defects/ and on boundaries/ is also important.Scattering rates are additive; thus the total phonon scattering rate is or alternatively, in terms of the mean-free-path (84), Expressions for each of the scattering rates are developed [17].For scattering from point defects, known as the Rayleigh scattering which is like the scattering of light from the dust.For boundaries and surfaces, where  is the shortest dimension of the sample.A commonly used expression for -scattering is With this background we are now able to understand the temperature dependence of the lattice thermal conductivity.

Lattice Thermal Conductivity versus
Temperature.The temperature dependence of the lattice thermal conductivity   is illustrated for bulk Si on Figure 4.According to (73)   is proportional to the number of the phonon modes that are occupied ⟨ ph /⟩ and to the average value of the phonon mean-free-path ⟨⟨ ph ⟩⟩.The curve   () can be explained by understanding how ⟨ ph ⟩ and ⟨⟨ ph ⟩⟩ vary with temperature.) Figure 4: The experimental [18] and calculated [13] thermal conductivity of bulk Si as a function of temperature.
It can be shown using ( 72) that, at low temperatures, so the initial rise in thermal conductivity is due to the fact that the number of populated modes rises quickly with temperature.At low temperatures boundary scattering is important.As the temperature increases more short-wavelength phonons are produced.These phonons scatter from point defects, so defect scattering becomes more and more important.As the temperature approaches the   , all of the phonon modes are populated and further increases in temperature do not change ⟨ ph ⟩.Instead, the higher temperatures increase the phonon scattering by -processes, and the thermal conductivity drops with increasing temperatures.

Difference between Lattice Thermal and Electrical Conductivities.
We have already noted the similarity between phonon transport equations ( 69) and (73) and electron transport equations ( 75) and (76).The average electron and phonon mean-free-paths are of the same order of magnitude.Why then does the electrical conductance vary over many more orders of magnitude while the lattice thermal conductance varies only over a few?The answer lies in the corresponding windows functions (60) and (62).For both electrons and phonons, higher temperatures broaden the window function and increase the population of states.For electrons, however, the position of the Fermi level has a dramatic effect on the magnitude of the window function.
By controlling the position of the Fermi level, the electrical conductivity can be varied over many orders of magnitude.
For the phonons the width of the window function is determined by temperature only.
Another key difference between electrons and phonons relates to how the states are populated.For a thermoelectric device   ≈   , and the electron and phonon window functions are quite similar.However, for electrons the BW of the dispersion is very large, so only a few states near the bottom of the conduction band are populated: the effective mass approximation works well for these states, and it is easy to obtain analytical solutions.For phonons, the BW of the dispersion is small.At moderate temperatures states all across the entire Brillouin zone are occupied: simple analytical approximations do not work, and it is hard to get analytical solutions for the lattice thermal conductivity.

3.7.
Lattice Thermal Conductivity Quantization.By analogy with the electronic conduction quantization, over 30 years ago Pendry [19] stated the existence of the quantum limits to the heat flow.In fact, if  → 0 in (63), then the phonon window  ph (ℏ) is sharply peaked near ℏ = 0: For a bulk conductor  ph (ℏ) → 0 as ℏ → 0, but for nanoresistors like a nanowire or nanoribbon one can have a finite number of phonon modes.For ballistic phonon transport,  ph = 1 and one can expect that Exactly this result was proved experimentally using 4mode resistor at  < 0.8 K [10]; thermoconductivity measurements agree with the predictions for 1D ballistic resistors [20][21][22].
The quantum of thermal conductance represents the maximum possible value of energy transported per phonon mode.Surprisingly, it does not depend on particle statistics: the quantum of thermal conductance is universal for fermions, bosons, and anyons [23][24][25].

Conclusions
In summary, we see that the LDL concept used to describe electron transport can be generalized for phonons.In both cases the Landauer approach generalized and extended by Datta and Lundstom gives correct quantitative description of transport processes for resistors of any nature and any dimension and size in ballistic, quasiballistic, and diffusive linear response regimes when there are differences in both voltage and temperature across the device.We saw that the lattice thermal conductivity can be written in a form that is very similar to the electrical conductivity, but there are two important differences.The first difference between electrons and phonons is the difference in bandwidths of their dispersions.For electrons, the dispersion BW ≫  at room temperature, so only low energy states are occupied.For phonons, BW ≈ , so at room temperature all of the acoustic modes across the entire Brillouin zone are occupied.As a result, the simple Debye approximation to the acoustic phonon dispersion does not work nearly as well as the simple effective mass approximation to the electron dispersion.
The second difference between electrons and phonons is that for electrons the mode populations are controlled by the window function which depends on the position of the Fermi level and the temperature.For phonons, the window function depends only on the temperature.The result is that electrical conductivities vary over many orders of magnitude as the position of the Fermi level varies, while lattice conductivities vary over only a few orders of magnitude.
Finally, we also collect below the thermoelectric coefficients for parabolic band semiconductors and for graphene [5,26].

Figure 1 :
Figure 1: Comparison of the density of states () and number of modes () for 1D, 2D, and 3D resistors with parabolic dispersion.

Figure 2 :
Figure 2: Fermi function and the dimensionless normalized thermal broadening function.

Figure 3 :
Figure 3: Broadening function for phonons compared to that of electrons.