Skip to content

Drude Theory

Drude theory is the single-relaxation-time model of charge-current response. It assumes that a carrier’s mean drift momentum accelerates under an electric field and relaxes exponentially toward equilibrium:

dpddt=qE−pdτ.\frac{d\mathbf p_d}{dt} = q\mathbf E - \frac{\mathbf p_d}{\tau}.

For an isotropic parabolic band, pd=m∗vd\mathbf p_d=m^*\mathbf v_d. The resulting dc conductivity is

σ0=nq2τm∗.\sigma_0 = \frac{nq^2\tau}{m^*}.

This formula looks classical, but its successful modern use is usually Drude–Sommerfeld: carrier density comes from Fermi statistics and band filling, m∗m^* from a band or quasiparticle dispersion, and τ\tau from momentum-relaxing scattering near the Fermi surface. The model does not predict those inputs. It organizes them into a low-frequency response.

This page owns the scalar and tensor Drude equations, relaxation-time interpretation, dc and ac conductivity, Hall response, optical spectral weight, fitting workflow, and failure tests. Transport Coefficients Preview owns generic transport limits and hydrodynamic classification, Kubo Formula owns the exact linear-response construction, and Effective Mass owns curvature, conductivity, cyclotron, and optical mass distinctions.

Required background. Effective Mass supplies the response-specific mass distinctions used here without rederivation.

Helpful background. Fermi Surface supplies low-energy phase space, Metals, Insulators, and Semiconductors supplies carrier and filling classification, and Transport Coefficients Preview supplies the bulk, ballistic, hydrodynamic, and limit-order taxonomy.

Let qq be the signed carrier charge. Thus q=−eq=-e for an electron and q=+eq=+e for a hole. The microscopic carrier velocity can be large even in equilibrium, but opposite directions cancel. vd\mathbf v_d denotes only the small nonequilibrium drift average.

For an isotropic band and zero magnetic field,

m∗dvddt=qE(t)−m∗vdτ.m^* \frac{d\mathbf v_d}{dt} = q\mathbf E(t) - \frac{m^*\mathbf v_d}{\tau}.

The relaxation term asserts:

  • linear response in the drift velocity;
  • local, memoryless damping;
  • one momentum-relaxation time τ\tau;
  • a fixed carrier density and effective mass;
  • spatially uniform response;
  • no distinction among directions, bands, or points on the Fermi surface.

Those are model assumptions, not definitions of a metal.

After transients from the remote past have decayed,

vd(t)=qm∗∫−∞tdt′ e−(t−t′)/τE(t′).\mathbf v_d(t) = \frac{q}{m^*} \int_{-\infty}^{t} dt'\, e^{-(t-t')/\tau} \mathbf E(t').

The response kernel is causal: only earlier fields contribute.

If a constant field E0\mathbf E_0 is switched on at t=0t=0 and the initial drift is zero,

vd(t)=qτm∗E0(1−e−t/τ).\mathbf v_d(t) = \frac{q\tau}{m^*} \mathbf E_0 \left( 1-e^{-t/\tau} \right).

If the field is switched off after steady state is reached, the current decays as e−t/τe^{-t/\tau}. This is the operational meaning of the model’s relaxation time.

τ\tau is a transport momentum-relaxation time. It need not equal:

  • the interval between all microscopic collisions;
  • a single-particle quantum lifetime from spectral broadening;
  • a phase-coherence time;
  • an energy-relaxation time;
  • a recombination time;
  • an escape time from a finite device.

A collision that changes momentum only slightly can strongly broaden a single-particle phase while relaxing current inefficiently. Naming the lifetime is part of reporting it.

In a time-independent field, the steady drift is

vd=qτm∗E.\mathbf v_d = \frac{q\tau}{m^*} \mathbf E.

The current density is

j=nqvd.\mathbf j = nq\mathbf v_d.

Although electron drift is opposite to E\mathbf E, multiplying by the negative electron charge makes the conventional current parallel to E\mathbf E. Therefore

j=σ0E,σ0=nq2τm∗.\mathbf j = \sigma_0\mathbf E, \qquad \sigma_0 = \frac{nq^2\tau}{m^*}.

The resistivity is

ρ0=1σ0=m∗nq2τ.\rho_0 = \frac{1}{\sigma_0} = \frac{m^*}{nq^2\tau}.

Define the positive mobility

μ=∣q∣τm∗.\mu = \frac{|q|\tau}{m^*}.

Then

σ0=n∣q∣μ.\sigma_0 = n|q|\mu.

Conductivity is a bulk coefficient. Conductance also depends on sample geometry and contacts. For a uniform bar of length LL and cross-sectional area AA,

G=σ0AL,G = \sigma_0 \frac{A}{L},

only when current distribution and contact effects permit the elementary geometry reduction.

The original Drude model treated carriers as a classical gas. That picture gives incorrect electronic heat capacity and misses Pauli blocking, band structure, and Fermi-surface geometry. Modern use keeps the momentum-balance equation while replacing its inputs.

nn is obtained from the occupied electronic states, not from a classical Maxwell distribution. In a simple monovalent metal it may be comparable to one carrier per atom; in a doped semiconductor or semimetal it can be many orders of magnitude smaller. In a multiband system there is no single nn unless a reduction is justified.

m∗m^* is a band, optical, or quasiparticle response parameter. It can be anisotropic, density dependent, or renormalized by interactions. Inserting the bare electron mass merely because the charge carrier is an electron is generally unjustified.

In a degenerate metal, the characteristic velocity for transport is the Fermi velocity rather than a classical thermal velocity. The transport mean free path is

ℓtr=vFτtr.\ell_{\mathrm{tr}} = v_{\mathrm F}\tau_{\mathrm{tr}}.

This connects a time-domain Drude parameter with a spatial semiclassical criterion.

Boltzmann Transport owns the distribution function, collision operator, scattering-in terms, thermoelectric moments, and conservation tests. Its simplest weak-field relaxation-time reduction gives

σij=q2∑n∫ddk(2π)d(−∂f∂ε)×vn,i(k)vn,j(k)τn(k).\begin{aligned} \sigma_{ij} = q^2 \sum_n \int \frac{d^dk}{(2\pi)^d} &\left( - \frac{\partial f} {\partial\varepsilon} \right) \\ &\times v_{n,i}(\mathbf k) v_{n,j}(\mathbf k) \tau_n(\mathbf k). \end{aligned}

Only states in the thermal window around the chemical potential contribute through −∂f/∂ε-\partial f/\partial\varepsilon. For one isotropic parabolic band and constant τ\tau, the integral reduces to nq2τ/m∗nq^2\tau/m^*.

For an isotropic degenerate band, the same physics can be written as an Einstein relation,

σ=q2ρ(EF)D,\sigma = q^2 \rho(E_{\mathrm F}) D,

with

D=vF2τd.D = \frac{v_{\mathrm F}^2\tau}{d}.

Here ρ(EF)\rho(E_{\mathrm F}) is the density of states per volume, including the degeneracies used in nn. The Drude and diffusion forms agree only when their masses, state counts, and scattering approximations describe the same carriers.

For a parabolic band with mass tensor M\mathsf M and scalar relaxation time,

dvddt+vdτ=qM−1E.\frac{d\mathbf v_d}{dt} + \frac{\mathbf v_d}{\tau} = q\mathsf M^{-1}\mathbf E.

The dc conductivity tensor is

σ0=nq2τM−1.\boldsymbol\sigma_0 = nq^2\tau \mathsf M^{-1}.

In principal axes,

σii=nq2τmi.\sigma_{ii} = \frac{nq^2\tau}{m_i}.

If τ\tau is also anisotropic or momentum dependent, it cannot in general be represented by multiplying M−1\mathsf M^{-1} by one scalar. A measured conductivity anisotropy can come from mass, Fermi-surface geometry, scattering, carrier density, or all four.

Effective Mass distinguishes the tensor, conductivity, density-of-states, cyclotron, and optical masses used in such reductions.

Use the convention

E(t)=Re⁡[Eωe−iωt].\mathbf E(t) = \operatorname{Re} \left[ \mathbf E_\omega e^{-i\omega t} \right].

The drift equation gives

vω=q/m∗τ−1−iωEω.\mathbf v_\omega = \frac{q/m^*} {\tau^{-1}-i\omega} \mathbf E_\omega.

Therefore

σ(ω)=nq2/m∗τ−1−iω=σ01−iωτ.\sigma(\omega) = \frac{nq^2/m^*} {\tau^{-1}-i\omega} = \frac{\sigma_0} {1-i\omega\tau}.

The pole lies at ω=−i/τ\omega=-i/\tau, in the lower half-plane as required for a retarded response.

The dissipative and reactive parts are

Re⁡σ(ω)=σ01+ω2τ2,\operatorname{Re}\sigma(\omega) = \frac{\sigma_0} {1+\omega^2\tau^2}, Im⁡σ(ω)=σ0ωτ1+ω2τ2.\operatorname{Im}\sigma(\omega) = \frac{\sigma_0\omega\tau} {1+\omega^2\tau^2}.

The cycle-averaged absorbed power density is

P‾=12Re⁡σ(ω)∣Eω∣2.\overline P = \frac{1}{2} \operatorname{Re}\sigma(\omega) |\mathbf E_\omega|^2.

Re⁡σ\operatorname{Re}\sigma is a Lorentzian centered at zero frequency with half width

γ=1τ.\gamma = \frac{1}{\tau}.

For ωτ≪1\omega\tau\ll1, the current nearly follows the field. For ωτ≫1\omega\tau\gg1, collisions are too slow to establish the dc drift during one cycle and the response is mainly reactive:

σ(ω)≈inq2m∗ω.\sigma(\omega) \approx i \frac{nq^2} {m^*\omega}.

Real and imaginary Drude conductivity versus frequency times relaxation time, beside the semicircle traced by the complex conductivity.

For σ/σ0=(1−iωτ)−1\sigma/\sigma_0=(1-i\omega\tau)^{-1}, the dissipative peak has half width 1/τ1/\tau, while the reactive part peaks at ωτ=1\omega\tau=1. At positive frequency the complex conductivity traces the upper semicircle from σ/σ0=1\sigma/\sigma_0=1 toward the origin.

Write

D0=nq2m∗.\mathcal D_0 = \frac{nq^2}{m^*}.

Then

Re⁡σ(ω)=D0γγ2+ω2.\operatorname{Re}\sigma(\omega) = \mathcal D_0 \frac{\gamma} {\gamma^2+\omega^2}.

The positive-frequency intraband area is

∫0∞dω Re⁡σ(ω)=π2nq2m∗.\int_0^\infty d\omega\, \operatorname{Re}\sigma(\omega) = \frac{\pi}{2} \frac{nq^2}{m^*}.

Changing τ\tau changes the width and height but not this Drude area. In the clean limit,

lim⁡γ→0+γγ2+ω2=πδ(ω)\lim_{\gamma\to0^+} \frac{\gamma} {\gamma^2+\omega^2} = \pi\delta(\omega)

on the full frequency axis. Thus

Re⁡σ(ω)⟶πD0δ(ω).\operatorname{Re}\sigma(\omega) \longrightarrow \pi\mathcal D_0 \delta(\omega).

Different communities absorb factors of π\pi into the definition of “Drude weight.” A quoted Drude weight must include its conductivity convention.

In a continuum parabolic model, the full optical ff-sum is governed by the bare mass. In a lattice, low-frequency intraband weight can involve band curvature and interactions, while interband transitions carry the remaining spectral weight. A partial Drude integral is not automatically the full sum rule.

Sum Rules owns exact integrated spectral constraints. The delta-function distinction between an ideal conductor and an insulator is broader than a finite-τ\tau fit and connects to Metals, Insulators, and Semiconductors.

For the e−iωte^{-i\omega t} convention,

ϵ(ω)=ϵ∞+iσ(ω)ϵ0ω.\epsilon(\omega) = \epsilon_\infty + \frac{i\sigma(\omega)} {\epsilon_0\omega}.

Define the Drude plasma frequency

ωp2=nq2ϵ0m∗.\omega_p^2 = \frac{nq^2} {\epsilon_0m^*}.

Then

ϵ(ω)=ϵ∞−ωp2ω(ω+i/τ).\epsilon(\omega) = \epsilon_\infty - \frac{\omega_p^2} {\omega \left( \omega+i/\tau \right)}.

The zero of the longitudinal dielectric function controls a bulk plasma mode after background screening and damping are included. Reflectivity edges, loss-function peaks, and zeros of Re⁡ϵ\operatorname{Re}\epsilon need not coincide when absorption or interband structure is appreciable.

Plasmons Preview owns self-consistent charge oscillations, dimensional scaling, dielectric zeros, and damping.

Include a static magnetic field:

m∗dvddt=q(E+vd×B)−m∗vdτ.m^* \frac{d\mathbf v_d}{dt} = q \left( \mathbf E + \mathbf v_d\times\mathbf B \right) - \frac{m^*\mathbf v_d}{\tau}.

Take B=Bz^\mathbf B=B\hat{\mathbf z} and define the signed cyclotron frequency and field parameter

ωc=qBm∗,βc=ωcτ.\omega_c = \frac{qB}{m^*}, \qquad \beta_c = \omega_c\tau.

In steady state,

σ=σ01+βc2(1βc−βc1).\boldsymbol\sigma = \frac{\sigma_0} {1+\beta_c^2} \begin{pmatrix} 1 & \beta_c \\ -\beta_c & 1 \end{pmatrix}.

The resistivity tensor is

ρ=1σ0(1−βcβc1).\boldsymbol\rho = \frac{1}{\sigma_0} \begin{pmatrix} 1 & -\beta_c \\ \beta_c & 1 \end{pmatrix}.

With the convention

RH=EyjxB=ρyxB,R_H = \frac{E_y} {j_xB} = \frac{\rho_{yx}}{B},

the single-carrier Hall coefficient is

RH=1nq.R_H = \frac{1}{nq}.

It is negative for electrons and positive for holes. The Hall angle has magnitude

∣tan⁡θH∣=∣ωcτ∣.|\tan\theta_H| = |\omega_c\tau|.

The same one-band isotropic model predicts

ρxx(B)=ρ0.\rho_{xx}(B) = \rho_0.

That is, it has zero transverse magnetoresistance despite a field-dependent σxx\sigma_{xx}. Inverting the conductivity tensor matters. Observed magnetoresistance requires additional structure such as multiple carriers, anisotropy, momentum-dependent scattering, open orbits, inhomogeneity, or quantum effects.

For carrier species aa with positive mobility μa\mu_a and sign sa=sgn⁡(qa)s_a=\operatorname{sgn}(q_a),

σxx(B)=∑ana∣qa∣μa1+μa2B2,\sigma_{xx}(B) = \sum_a \frac{ n_a|q_a|\mu_a }{ 1+\mu_a^2B^2 }, σxy(B)=∑asana∣qa∣μa2B1+μa2B2.\sigma_{xy}(B) = \sum_a \frac{ s_a n_a|q_a| \mu_a^2B }{ 1+\mu_a^2B^2 }.

For electrons of density nn and holes of density pp, the weak-field Hall coefficient is

RH=pμh2−nμe2e(pμh+nμe)2.R_H = \frac{ p\mu_h^2 - n\mu_e^2 }{ e \left( p\mu_h+n\mu_e \right)^2 }.

The Hall sign is mobility weighted. A positive RHR_H does not prove that every Fermi-surface sheet is hole-like, and 1/(eRH)1/(eR_H) is not generally the total carrier density. Fermi Surface develops this warning in the context of electron and hole pockets.

The Drude equation does not calculate τ\tau. Candidate momentum-relaxing mechanisms include:

  • impurities, defects, and isotope disorder;
  • phonons;
  • electron–electron Umklapp;
  • boundaries and interfaces;
  • magnetic fluctuations;
  • coupling to other collective modes;
  • spatial inhomogeneity.

For elastic scattering on an isotropic Fermi surface, a representative transport rate is

1τtr=∑k′Wkk′(1−cos⁡θ).\frac{1}{\tau_{\mathrm{tr}}} = \sum_{\mathbf k'} W_{\mathbf k\mathbf k'} \left( 1-\cos\theta \right).

θ\theta is the angle through which the velocity is scattered. The corresponding single-particle or quantum rate lacks the current-relaxation factor:

1τq=∑k′Wkk′.\frac{1}{\tau_q} = \sum_{\mathbf k'} W_{\mathbf k\mathbf k'}.

Forward scattering can make τtr≫τq\tau_{\mathrm{tr}}\gg\tau_q. Comparing mobility and quantum-oscillation damping as though they measured one lifetime can therefore be misleading.

Independent weak scattering channels are often approximated by

1τtr≈∑a1τa.\frac{1}{\tau_{\mathrm{tr}}} \approx \sum_a \frac{1}{\tau_a}.

This is Matthiessen’s rule. It can fail when channels interfere, scattering is strongly energy dependent or anisotropic, the band structure changes with temperature, or one process modifies the distribution sampled by another.

A measured temperature dependence of ρ\rho is not predicted by Drude theory until a microscopic or empirical τ(T)\tau(T) is supplied.

The model survives because it is more general as a momentum-balance structure than its classical origin suggests.

If the uniform current overlaps mainly with one slowly relaxing quantity, its low-frequency response often has one dominant pole:

σ(ω)∼DΓ−iω.\sigma(\omega) \sim \frac{\mathcal D} {\Gamma-i\omega}.

Microscopic quantum mechanics determines the weight D\mathcal D and relaxation rate Γ\Gamma. The line shape can look Drude-like even when carriers are Bloch quasiparticles rather than classical balls between collisions.

Fermi statistics repair the thermodynamics

Section titled “Fermi statistics repair the thermodynamics”

Only states near the Fermi surface respond to a weak field. This removes the classical prediction that every conduction electron contributes a thermal kinetic energy of order kBTk_{\mathrm B}T while preserving the simple conductivity reduction in a parabolic band.

Many collisions rapidly establish local equilibrium while relaxing total momentum only weakly. Uniform conductivity is then governed by the slower momentum-relaxing processes. A single long τ\tau can summarize this bottleneck over a limited frequency range.

If experimental frequencies satisfy ωτ≪1\omega\tau\ll1, many details collapse into the measured ratio nτ/m∗n\tau/m^*. Good agreement with ρ0\rho_0 alone therefore validates that combination, not each microscopic parameter separately.

For a degenerate Fermi gas with sufficiently elastic, weakly energy-dependent scattering,

κeσT⟶L0,L0=π23(kBe)2.\frac{\kappa_e} {\sigma T} \longrightarrow L_0, \qquad L_0 = \frac{\pi^2}{3} \left( \frac{k_{\mathrm B}}{e} \right)^2.

The quantum Sommerfeld calculation supplies the correct Lorenz number. Deviations can arise from inelastic scattering, additional heat carriers, multiband structure, hydrodynamics, or non-Fermi-liquid behavior.

A general causal intraband response can be parameterized as

σ(ω)=ϵ0ωp2M(ω)−iω,\sigma(\omega) = \frac{\epsilon_0\omega_p^2} {M(\omega)-i\omega},

where M(ω)M(\omega) is a memory function. Ordinary Drude theory takes

M(ω)=1τ.M(\omega) = \frac{1}{\tau}.

Given a chosen intraband plasma frequency, optical analyses sometimes define

1τop(ω)=ϵ0ωp2Re⁡[1σ(ω)],\frac{1} {\tau_{\mathrm{op}}(\omega)} = \epsilon_0\omega_p^2 \operatorname{Re} \left[ \frac{1}{\sigma(\omega)} \right], mop(ω)mb=−ϵ0ωp2ωIm⁡[1σ(ω)].\frac{m_{\mathrm{op}}(\omega)} {m_b} = - \frac{\epsilon_0\omega_p^2} {\omega} \operatorname{Im} \left[ \frac{1}{\sigma(\omega)} \right].

For a perfect Drude response these give 1/τ1/\tau and 11. For real data they are derived optical functions, not direct measurements of a microscopic collision clock and mass. Their values depend on the chosen plasma frequency and on how interband, phonon, and collective contributions are subtracted.

The model is unreliable when its one-fluid, local, memoryless assumptions fail.

Different densities, masses, charges, and lifetimes produce several response scales. One Drude peak can fit a limited interval while assigning physically false nn, m∗m^*, and τ\tau.

Current is a Fermi-surface average of velocities and scattering. A scalar mass cannot represent strongly warped surfaces, open orbits, or momentum-dependent τ\tau.

Optical conductivity can contain interband transitions, phonons, excitons, plasmons, density-wave modes, and superconducting condensate weight. Fitting all low-energy structure with one broadened Drude term violates spectral accounting.

Localized carriers do not support the same extended semiclassical drift. Activated hopping, variable-range hopping, weak localization, and Anderson localization require different frequency and temperature dependence.

When

kFℓtr≲1,k_{\mathrm F}\ell_{\mathrm{tr}} \lesssim 1,

the picture of well-defined momentum between separated collisions becomes suspect. This Mott–Ioffe–Regel criterion is a heuristic crossover, not a universal sharp bound. Some correlated systems have broad, non-Drude conductivity without long-lived electron quasiparticles.

Momentum-conserving electron–electron collisions can be much faster than momentum relaxation. The system can form a viscous electron fluid whose geometry-dependent flow is not captured by one local τ\tau.

When device dimensions are shorter than the relevant mean free path, leads, channels, contacts, and transmission probabilities control conductance. A bulk Drude resistivity is no longer the complete protocol.

When Landau quantization, quantum oscillations, magnetic breakdown, or edge states matter, the smooth classical Lorentz-force solution is incomplete.

Berry curvature, skew scattering, side jump, magnetic texture, and quantized Hall response generate transverse currents not reducible to RH=1/(nq)R_H=1/(nq).

A superconductor transfers spectral weight into a zero-frequency condensate response and develops gap-dependent absorption. A narrow normal Drude peak alone cannot represent phase stiffness or the superconducting state.

Large fields can heat or redistribute carriers, while short wavelengths and anomalous skin effects make j(r)\mathbf j(\mathbf r) depend on fields elsewhere. The local linear relation j=σE\mathbf j=\sigma\mathbf E then needs replacement.

From ρ0\rho_0 alone one obtains only

nτm∗=1q2ρ0.\frac{n\tau}{m^*} = \frac{1} {q^2\rho_0}.

Separating nn, m∗m^*, and τ\tau requires independent information. Hall density is independent only in a validated one-carrier regime. Optical weight, quantum oscillations, ARPES, capacitance, chemistry, and band calculations each bring their own assumptions.

Measure both longitudinal and transverse voltages, antisymmetrize the Hall signal in BB, symmetrize the longitudinal signal, and invert the full resistivity tensor when fitting conductivity components. Nonlinearity in ρyx(B)\rho_{yx}(B) is evidence against the weakest one-carrier reduction, but a linear curve alone does not prove it.

Fit complex σ(ω)\sigma(\omega) or an electrodynamically consistent reflectivity/transmission model, not only one visible peak. State:

  • the Fourier convention;
  • frequency range and resolution;
  • background dielectric constant;
  • interband and phonon components;
  • number of Drude terms;
  • fitted plasma frequencies and scattering rates;
  • partial-sum cutoff;
  • Kramers–Kronig or causal-model procedure.

Optics measures Drude weight and width. It does not separately determine nn and m∗m^* without another constraint.

ObservableSimple Drude inferenceRequired warning
ρ0\rho_0m∗/(nq2τ)m^*/(nq^2\tau)only one parameter combination
weak-field RHR_H1/(nq)1/(nq)one isotropic carrier
mobility μ\mu$q
optical width1/τ1/\tautransport and optical lifetime may differ
Drude areaπnq2/(2m∗)\pi nq^2/(2m^*)partial weight and background subtraction
magnetoresistancezerononzero response signals physics beyond scalar one-band Drude

A trustworthy interpretation seeks one set of carriers and conventions that explains several observables, while allowing response-specific masses and lifetimes where theory requires them.

  • Using the bare electron mass without checking the band or optical mass.
  • Calling every collision time a transport relaxation time.
  • Inferring n=1/(eRH)n=1/(eR_H) in a multiband or anomalous-Hall material.
  • Forgetting to invert the conductivity or resistivity tensor in a magnetic field.
  • Fitting only Re⁡σ\operatorname{Re}\sigma without a causal complex response.
  • Treating a numerical broadening as 1/τ1/\tau.
  • Adding scattering rates by Matthiessen’s rule without testing independence.
  • Using a classical thermal velocity instead of vFv_{\mathrm F} for a degenerate metal.
  • Reading a successful dc fit as proof of coherent quasiparticles.
  • Calling a zero-frequency delta function a finite resistivity.
  • Ignoring interband spectral weight when applying an optical sum.
  • Reporting a mobility, plasma frequency, or Hall density without units, geometry, temperature, field, and model.

A uniform field E0\mathbf E_0 is switched on at t=0t=0. If j(0)=0\mathbf j(0)=0, derive j(t)\mathbf j(t).

Solution

The current obeys

djdt+jτ=nq2m∗E0.\frac{d\mathbf j}{dt} + \frac{\mathbf j}{\tau} = \frac{nq^2}{m^*} \mathbf E_0.

The solution is

j(t)=σ0E0(1−e−t/τ),\mathbf j(t) = \sigma_0\mathbf E_0 \left( 1-e^{-t/\tau} \right),

where

σ0=nq2τm∗.\sigma_0 = \frac{nq^2\tau}{m^*}.

At short times the current grows approximately linearly; at long times it reaches the dc value.

Verify the positive-frequency area of the dissipative conductivity.

Solution

Using

Re⁡σ(ω)=σ01+ω2τ2,\operatorname{Re}\sigma(\omega) = \frac{\sigma_0} {1+\omega^2\tau^2},

set x=ωτx=\omega\tau:

∫0∞dω Re⁡σ=σ0τ∫0∞dx1+x2=σ0τπ2=π2nq2m∗.\begin{aligned} \int_0^\infty d\omega\, \operatorname{Re}\sigma &= \frac{\sigma_0}{\tau} \int_0^\infty \frac{dx}{1+x^2} \\ &= \frac{\sigma_0}{\tau} \frac{\pi}{2} \\ &= \frac{\pi}{2} \frac{nq^2}{m^*}. \end{aligned}

The area is independent of τ\tau.

Take

n=5.0×1028 m−3,ρ0=2.0×10−8 Ω m,n = 5.0\times10^{28}\,\mathrm{m^{-3}}, \qquad \rho_0 = 2.0\times10^{-8}\,\Omega\,\mathrm m,

with m∗=mem^*=m_e and vF=1.5×106 m s−1v_{\mathrm F}=1.5\times10^6\,\mathrm{m\,s^{-1}}. Estimate τ\tau, mobility, and mean free path.

Solution

From ρ0=me/(ne2τ)\rho_0=m_e/(ne^2\tau),

τ=mene2ρ0≈3.6×10−14 s.\tau = \frac{m_e} {ne^2\rho_0} \approx 3.6\times10^{-14}\,\mathrm s.

The mobility is

μ=eτme≈6.3×10−3 m2 V−1 s−1,\mu = \frac{e\tau}{m_e} \approx 6.3\times10^{-3}\, \mathrm{m^2\,V^{-1}\,s^{-1}},

or about 63 cm2 V−1 s−163\,\mathrm{cm^2\,V^{-1}\,s^{-1}}. The mean free path is

ℓ=vFτ≈5.4×10−8 m,\ell = v_{\mathrm F}\tau \approx 5.4\times10^{-8}\,\mathrm m,

about 54 nm54\,\mathrm{nm}. These values inherit the assumed one-band mass and density.

For one electron-like carrier, determine the Hall coefficient and longitudinal resistivity in a magnetic field.

Solution

With q=−eq=-e,

RH=1nq=−1ne.R_H = \frac{1}{nq} = -\frac{1}{ne}.

The inverse of the one-band conductivity tensor gives

ρxx(B)=1σ0=m∗ne2τ,\rho_{xx}(B) = \frac{1}{\sigma_0} = \frac{m^*}{ne^2\tau},

independent of BB. Thus the isotropic one-carrier Drude model has a negative Hall coefficient and zero transverse magnetoresistance.

A compensated semimetal has n=pn=p, but μh=2μe\mu_h=2\mu_e. What is the sign of its weak-field Hall coefficient?

Solution

For n=pn=p,

RH=nμh2−nμe2en2(μh+μe)2.R_H = \frac{ n\mu_h^2-n\mu_e^2 }{ en^2 \left( \mu_h+\mu_e \right)^2 }.

Since

μh2−μe2=4μe2−μe2>0,\mu_h^2-\mu_e^2 = 4\mu_e^2-\mu_e^2 > 0,

RHR_H is positive. The material is exactly carrier compensated, yet the more mobile holes control the weak-field Hall sign.

Suppose scattering is concentrated at a small angle θ0≪1\theta_0\ll1. Compare its contributions to 1/τq1/\tau_q and 1/τtr1/\tau_{\mathrm{tr}}.

Solution

The quantum rate counts the event with weight approximately one. The transport rate includes

1−cos⁡θ0≈θ022.1-\cos\theta_0 \approx \frac{\theta_0^2}{2}.

Thus the same forward-scattering channel contributes much less to current relaxation:

τtr−1τq−1∼θ022.\frac{ \tau_{\mathrm{tr}}^{-1} }{ \tau_q^{-1} } \sim \frac{\theta_0^2}{2}.

Consequently τtr\tau_{\mathrm{tr}} can be much longer than τq\tau_q.

An optical conductivity has a narrow zero-frequency peak, a broad mid-infrared band, and a frequency-dependent 1/τop(ω)1/\tau_{\mathrm{op}}(\omega). Is a single Drude fit a complete interpretation?

Solution

No. A defensible analysis should:

  1. fit the full complex response over a declared range;
  2. test whether more than one intraband component is required;
  3. identify possible interband, phonon, collective, or localization contributions to the mid-infrared band;
  4. check spectral-weight transfer as temperature or tuning changes;
  5. state the plasma frequency used to define the extended-Drude functions;
  6. compare dc conductivity with the ω→0\omega\to0 limit;
  7. test Kramers–Kronig consistency and resolution effects;
  8. compare carrier densities and masses with Hall, quantum-oscillation, photoemission, or band information.

The narrow component may be Drude-like, but the broad band and frequency-dependent memory show that one constant τ\tau does not describe the entire spectrum.

  • Terahertz and Infrared Probes develops the field- and power-resolved experimental workflow that turns a sample stack into complex optical conductivity and tests a Drude fit.
  • Transport Measurements shows how terminal records, contact configurations, geometry, reversal, and tensor inversion produce the experimental conductivity to which a Drude model may be fitted.
  • Boltzmann Transport resolves state-dependent velocities, collision integrals, electric and thermal driving, and the limits of a diagonal relaxation time.
  • Disorder in Quantum Matter defines disorder ensembles and distinguishes the transport time controlling the Drude width from quantum and single-particle lifetimes.
  • Weak Localization adds the leading phase-coherent interference correction to this classical diffusive baseline and fixes its magnetoconductance sign convention.
  • Anderson Insulators replaces the extended-state Drude picture by localized orbitals, bath-assisted hopping, and Coulomb-gap transport.
  • Hall Effect develops sign conventions, multiband tensor inversion, anomalous mechanisms, experimental extraction, and the bridge to quantization.
  • Transport Coefficients Preview develops diffusion, Einstein relations, ballistic and hydrodynamic regimes, Green–Kubo formulas, and orders of limits.
  • Kubo Formula derives conductivity from equilibrium current response and keeps electromagnetic contact terms explicit.
  • Density Operators and Current Operators constructs microscopic currents and continuity equations.
  • Lifetime and Spectral Weight distinguishes transport, quantum, coherence, and decay times.
  • Fermi Surface connects carrier pockets, velocities, quantum oscillations, and multiband Hall warnings.
  • Effective Mass explains which response-specific mass belongs in a Drude reduction.
  • Fermi Liquid Theory Preview treats quasiparticle renormalization, backflow, and conservation constraints.
  • Strange Metals audits when a measured linear resistivity can be converted into a Planckian rate and why that conversion is model dependent.
  • Kondo Effect provides a concrete failure of a temperature-independent lifetime through many-body spin-flip scattering and strong-coupling saturation.
  • Quantum Matter Conventions fixes Fourier signs, charge signs, tensor order, and conductivity normalization.
  • Condensed Matter Roadmap places Drude theory between band structure and more general Boltzmann, Kubo, and mesoscopic transport.
  • P. Drude, “Zur Elektronentheorie der Metalle,” Annalen der Physik 306, 566–613 (1900), doi:10.1002/andp.19003060312.
  • A. Sommerfeld, “Zur Elektronentheorie der Metalle auf Grund der Fermischen Statistik. I,” Zeitschrift für Physik 47, 1–32 (1928), doi:10.1007/BF01391052.
  • N. W. Ashcroft and N. D. Mermin, Solid State Physics, Holt, Rinehart and Winston, 1976.
  • J. M. Ziman, Electrons and Phonons: The Theory of Transport Phenomena in Solids, Oxford University Press, 1960; Oxford Classic Texts reissue, 2001, doi:10.1093/acprof:oso/9780198507796.001.0001.
  • M. Dressel and G. Grüner, Electrodynamics of Solids, Cambridge University Press, 2002, doi:10.1017/CBO9780511606168.
  • R. Kubo, “Statistical-Mechanical Theory of Irreversible Processes. I,” Journal of the Physical Society of Japan 12, 570–586 (1957), doi:10.1143/JPSJ.12.570.
  • W. Kohn, “Theory of the Insulating State,” Physical Review 133, A171–A181 (1964), doi:10.1103/PhysRev.133.A171.
  • N. D. Mermin, “Lindhard Dielectric Function in the Relaxation-Time Approximation,” Physical Review B 1, 2362–2363 (1970), doi:10.1103/PhysRevB.1.2362.
  • D. N. Basov, R. D. Averitt, D. van der Marel, M. Dressel, and K. Haule, “Electrodynamics of Correlated Electron Materials,” Reviews of Modern Physics 83, 471–541 (2011), doi:10.1103/RevModPhys.83.471.