Boltzmann-Bloch Equation Approach to the Theory of the Optical Inter- and Intraband Response in Noble Metals

This paper introduces a momentum-resolved Boltzmann-Bloch equation framework that incorporates many-body interactions and anisotropic electronic dispersion to provide a microscopic description of the temperature-dependent linear optical response in gold, offering deeper insight into the interplay of intra- and interband processes than traditional phenomenological models.

The Gist

Imagine a bustling city made of gold. This city is the metal, and the people running around are electrons. Some of these people live in the "basement" (the valence band), and some live in the "penthouse" (the conduction band). Usually, the penthouse is only half-full, which is what makes gold a conductor of electricity.

When you shine a light on this gold city, you are essentially sending a wave of energy through the streets, causing the people (electrons) to move, bump into each other, and sometimes jump from the basement to the penthouse.

This paper is about building a super-detailed map to understand exactly how these people move when the light hits them.

The Old Map vs. The New Map

The Old Way (Drude-Lorentz Model):
For a long time, scientists used a simplified map. They treated the electrons like a single, giant crowd of people running through a foggy street. They knew the crowd slowed down because of collisions (friction), but they didn't really know why or how the collisions happened. They just guessed the rules to make the math work. It was like saying, "The traffic slows down because of 'traffic'," without looking at the individual cars.

The New Way (The Boltzmann-Bloch Approach):
The authors of this paper built a new, microscopic map. Instead of treating the crowd as one blob, they tracked individual electrons. They looked at:

1. Intra-band movement: People running around inside the penthouse (conduction band).
2. Inter-band movement: People jumping from the basement to the penthouse.
3. The Bumps: How electrons bump into other electrons or bump into the vibrating walls of the city (phonons/heat).

They call this the Metal Boltzmann-Bloch Equation (MBBE). Think of it as a high-tech GPS that tracks every single electron's speed, direction, and who they bumped into.

The Shape of the City (Anisotropic Dispersion)

Here is the tricky part: The gold city isn't a perfect sphere. It's shaped like a complex, multi-faceted gemstone with specific "highways" and "dead ends."

• The Problem: If you assume the city is a perfect sphere (isotropic), your map is wrong. The electrons don't move the same way in every direction.
• The Solution: The authors used an anisotropic model. Imagine the city is made of 14 different cone-shaped districts (around specific points called X and L). In some districts, the roads are wide and fast; in others, they are narrow and winding.
• Why it matters: When they tried to predict how gold reacts to light using a simple sphere, the colors were wrong. But when they used the complex, cone-shaped map, the predicted colors matched the real gold perfectly. It turns out, the weird shape of the city is the secret sauce to why gold looks the way it does.

The Collisions (Relaxation and Dephasing)

When the light hits the gold, the electrons get excited. But they don't stay excited forever; they crash into things and calm down. The paper breaks down these crashes into two types:

1. Bumping into the Walls (Electron-Phonon): The city walls are vibrating because of heat. As electrons run, they bump into these vibrating walls.
   • The Finding: This is the main reason the electrons slow down. It's like running on a treadmill that's shaking. The hotter the room (higher temperature), the more the walls shake, and the more the electrons get slowed down.
2. Bumping into Each Other (Electron-Electron): The electrons also crash into each other.
   • The Finding: This happens too, but it follows a different rule (it gets worse much faster as you heat things up, like a quadratic curve).

The Temperature Twist

The authors tested their map at different temperatures, from freezing cold to very hot.

• Cold: The walls stop shaking. The electrons can run faster, and the "friction" drops.
• Hot: The walls vibrate wildly. The electrons get hammered, and the gold absorbs more light (it gets darker in the infrared).

Their new map predicted exactly how the gold's color and reflectivity would change as it got hotter, matching real-world experiments they did in their lab.

The Big Takeaway

This paper is a victory for details over guesses.

By creating a microscopic map that accounts for the weird, cone-shaped geometry of the gold atoms and the specific ways electrons bump into each other and the heat, the authors can explain the optical properties of gold without needing to "fudge" the numbers.

In simple terms: They stopped guessing how gold behaves and started actually simulating the traffic jam of electrons inside it. This helps us design better solar cells, faster computer chips, and more sensitive medical sensors, because now we know exactly how light interacts with the "people" inside the metal.

Technical Summary

1. Problem Statement

Plasmonic applications rely heavily on the accurate description of the dielectric function $\epsilon(\omega)$ of noble metals. Historically, the optical response of these metals has been modeled using phenomenological approaches:

• Drude Model: Describes intraband (conduction band) processes but treats electrons as a free gas, ignoring band structure details.
• Lorentz Oscillators: Used to phenomenologically fit interband transitions (valence to conduction band) but lacks a microscopic foundation, hiding the underlying quantum mechanical mechanisms.

The Gap: There is a lack of a unified, microscopic framework that simultaneously describes intra- and interband processes, including many-body relaxation (electron-electron and electron-phonon scattering) and dephasing, while accounting for the complex, anisotropic geometry of the Fermi surface in noble metals. Existing models often fail to explain temperature-dependent spectral features or the specific shape of the interband absorption edge without arbitrary fitting parameters.

2. Methodology

The authors introduce a Momentum-Resolved Metal Boltzmann-Bloch Equation (MBBE) framework. This approach extends the Semiconductor Bloch Equations (SBE), typically used for semiconductors with empty conduction bands, to noble metals with partially filled conduction bands.

A. Theoretical Framework

• Hamiltonian: The system is described by a two-band Hamiltonian (valence $v$ and conduction $c$) including:
  • Single-particle dispersion.
  • Light-matter interaction (intra- and interband dipole couplings).
  • Many-body interactions: Electron-phonon (including Umklapp processes) and electron-electron (screened Coulomb) scattering.
• Equations of Motion: The authors derive coupled equations for:
  • Electron Occupations ($f^\lambda_k$): Describing intraband dynamics and relaxation.
  • Interband Transitions ($p^{vc}_k$): Describing the coherence between valence and conduction bands.
• Scattering Rates: Using a Born-Markov approximation and correlation expansion, they derive semi-analytical expressions for:
  • Relaxation rates ($\gamma^\lambda_k$): Energy and momentum relaxation of occupations.
  • Dephasing rates ($\gamma^p_k$): Decay of interband coherence.
  • These rates distinguish between Normal (N) and Umklapp (U) scattering processes for both electron-phonon and electron-electron interactions.

B. Anisotropic Band Structure Model

To accurately model the Fermi surface of gold (fcc lattice), the authors employ an anisotropic electronic dispersion model:

• The Brillouin zone is approximated by 14 high-symmetry cones centered at the X and L points.
• The dispersion relation $\epsilon^\lambda_{k}$ is expanded quadratically in radial ($k_\parallel$) and tangential ($k_\perp$) directions relative to these symmetry points.
• This model captures the "necks" at L-points and "domes" at X-points, which are critical for describing the density of states (DOS) and transition energies accurately.

C. Linearization and Calculation

• The MBBE are linearized for weak optical fields to calculate the linear optical response.
• The macroscopic polarization is derived from the microscopic currents and interband coherences.
• The total permittivity $\epsilon(\omega)$ is calculated as the sum of a background term, an intraband (Drude-like) term, and an interband term derived directly from the microscopic band structure and scattering rates.

3. Key Contributions

1. Unified Microscopic Framework: The MBBE provides the first comprehensive microscopic description that treats intra- and interband processes on equal footing, including full many-body scattering dynamics.
2. Anisotropic Dispersion Implementation: The paper demonstrates that the complex geometry of the Fermi surface (anisotropy) is essential for explaining the spectral shape of interband transitions, particularly the broadening of the absorption edge, which isotropic models fail to reproduce.
3. Semi-Analytical Scattering Rates: The derivation of temperature-dependent, momentum-resolved relaxation and dephasing rates for both electron-phonon and electron-electron interactions in a metal context.
4. Parameter Reduction: Unlike Drude-Lorentz models that require fitting multiple oscillator strengths and linewidths, the MBBE approach significantly reduces the number of free parameters, relying primarily on the band structure fit and fundamental material constants.

4. Results

The framework was applied to bulk gold and compared against experimental ellipsometry data across a temperature range of 75 K to 700 K.

• Relaxation and Dephasing Rates:

  • Electron-Phonon: Dominates intraband relaxation. The rate follows a linear temperature dependence ($T$) at high temperatures and a $T^5$ dependence at low temperatures (Bloch-Grüneisen relation).
  • Electron-Electron: Follows a $T^2$ law, consistent with Fermi liquid theory.
  • Dephasing: The interband dephasing rate is dominated by electron-phonon scattering. A notable "dip" near the Fermi edge is observed at low temperatures due to Pauli blocking.

• Optical Spectra (Permittivity):

  • Fit Range Selection: Three different anisotropic band structure fits were tested. The fit considering global band curvature (Range i) provided the best agreement with experimental data and the correct plasma frequency ($\hbar\omega_{pl} \approx 9.7$ eV). Local fits (Ranges ii and iii) overestimated the plasma frequency.
  • Temperature Dependence:
    • The imaginary part of the intraband susceptibility shows strong temperature dependence, vanishing as $T \to 0$ (consistent with reduced scattering).
    • The interband absorption edge shows a spectrally extended broadening. The study proves this broadening is not primarily due to temperature-dependent dephasing (which is weak at low $T$) but is instead a direct consequence of the anisotropic band structure.
  • Comparison with Isotropic Models: An isotropic model failed to reproduce the experimental absorption edge shape at low temperatures, confirming that the anisotropic geometry of the Fermi surface is the dominant factor in the spectral line shape.

5. Significance

• Fundamental Insight: The work bridges the gap between phenomenological Drude-Lorentz models and rigorous many-body theory, revealing that the "broadening" of interband transitions in noble metals is largely a geometric effect of the Fermi surface rather than just a thermal dephasing effect.
• Predictive Power: By reducing the number of fitting parameters, the MBBE approach offers a more predictive tool for designing plasmonic devices, especially those operating under non-equilibrium conditions or at varying temperatures.
• Future Applications: The framework is established as a foundation for studying non-linear optical processes, hot carrier generation, and the dynamics of electrons in plasmonic nanostructures far from equilibrium.

In conclusion, this paper successfully establishes a microscopic, momentum-resolved theory for the optical response of noble metals, demonstrating that accurate modeling of the anisotropic Fermi surface is crucial for explaining experimental optical spectra, particularly the temperature-dependent behavior of interband transitions.