Plasmonic properties and correlation energies from a compact multipole representation of the dielectric response in 2D metals

This paper generalizes the Multipole-Padé approximant framework to create a compact, symmetry-conserving, and anisotropic representation of the inverse dielectric function for 2D metals, enabling efficient and accurate calculation of plasmonic properties and correlation energies across the full Brillouin zone while bridging *ab initio* calculations with analytical models.

The Gist

Imagine a 2D metal (like a single layer of atoms) as a giant, bustling dance floor. When you tap on the floor, the dancers (electrons) don't just move individually; they ripple and wave together in a coordinated pattern. In physics, these collective waves are called plasmons, and the way the material responds to these waves is described by something called the dielectric function.

For a long time, scientists have had two ways to study this dance floor:

1. The "Brute Force" Method: They use supercomputers to calculate the movement of every single dancer at every single spot on the floor. This gives a massive amount of data—like a video recording with billions of frames. It's accurate, but it's huge, hard to read, and impossible to use for making new predictions quickly.
2. The "Simple Model" Method: They try to describe the whole dance with a simple rule, like "everyone moves in a circle." This is easy to use, but often too simple to capture the complex, real-life choreography of different materials.

What this paper does:
The authors, Dario A. Leon, Claudia Cardoso, and Kristian Berland, have created a new "smart summary" tool that sits perfectly between these two extremes. They call it a Multipole-Padé Approximant (MPA).

Think of their tool as a musical synthesizer.

• Instead of recording the entire orchestra (the brute force data), they figure out that the complex sound of the orchestra can be recreated perfectly by just a few specific notes played on a few specific instruments.
• In their case, they found that the complex "dance" of electrons in 2D metals can be described accurately by just a handful of collective modes (their "notes").

How it works (The Analogy):
Imagine you are trying to describe the shape of a bumpy hill (the electron response) to someone who has never seen it.

• The old way: You hand them a map with 1,000,000 dots showing the exact height at every single point. It's accurate, but they can't hold the map, and they can't easily guess what the hill looks like between the dots.
• The new way (This paper): You give them a smooth, flexible wire frame. You only need to bend this wire at a few key points (the "poles" or "modes") to make it match the hill perfectly. Once they have the wire frame, they can instantly see the shape of the hill from any angle, even in places where they didn't put a dot.

What they found:

1. It works for many different "dance floors": They tested this on seven different types of 2D metals, ranging from simple ones (like Sodium) to complex ones with multiple types of dancers (like Magnesium Boride).
2. Few notes are enough: Even for the complex materials, they only needed between 1 and 6 "notes" (modes) to recreate the entire dance floor's behavior perfectly.
3. It fills in the gaps: Because their model is a smooth mathematical formula (not just a list of dots), it can predict what happens in the "gaps" between the data points. This is crucial for calculating the correlation energy (a measure of how much energy the dancers save by moving together). Their method calculates this energy much faster and more accurately than the old "brute force" method, especially when looking at very small movements.

Why it matters:
This paper doesn't just give a pretty picture; it builds a bridge. It connects the heavy, slow, supercomputer calculations (the "brute force" data) with fast, easy-to-use mathematical models. Now, scientists can take the massive data from supercomputers, compress it into this "wire frame" summary, and use it to quickly predict how new materials will behave without needing to run the supercomputer again.

In short: They found a way to turn a million-page instruction manual on how electrons dance into a simple, 5-step recipe that works just as well.

Technical Summary

Technical Summary: Compact Multipole Representation of Dielectric Response in 2D Metals

Problem Statement
Two-dimensional (2D) metals exhibit complex, momentum-dependent screening and plasmonic excitations arising from reduced dimensionality and anisotropic electronic structures. While first-principles calculations based on the Random Phase Approximation (RPA) can access the dielectric response over wide frequency ranges, they typically generate large numerical datasets. These datasets are often difficult to interpret and impractical to use as a starting point for further modeling, particularly for properties requiring integration over the full Brillouin zone (BZ). Existing approaches using multipole-Padé approximants (MPA) have successfully captured plasmon dispersions for bulk metals and along one-dimensional momentum paths, but they lack the capability to evaluate properties requiring full 2D BZ integration.

Methodology
The authors extend the momentum-dependent multipole-Padé approximant (MPA(q)) framework to the full 2D Brillouin zone. The core of the method involves constructing a symmetry-conserving, anisotropic representation of the inverse dielectric function, $\epsilon^{-1}(q, \omega)$, over a wide energy range.

1. Analytic Formulation: The inverse dielectric function is expressed as a sum over a finite number of poles ($n_p$):
   $$\epsilon^{-1}_{\text{MPA}}(q, \omega) = 1 + \sum_{p}^{n_p} \frac{2R_p(q)\Omega_p(q)}{\omega^2 - [\Omega_p(q)]^2}$$
   where $R_p(q)$ and $\Omega_p(q)$ are complex, smooth functions of momentum $q$.

2. Symmetry and Anisotropy: Unlike previous formulations that used simple polynomials along high-symmetry paths, this work generalizes the model to the full 2D BZ using polar coordinates $(q, \theta)$. The parameters are decomposed into isotropic and anisotropic terms:

   • Isotropic terms: Modeled using third-order polynomials in $q$ to capture smooth momentum dependence and finite plasmon energy at $q \to 0$ (due to interband contributions).
   • Anisotropic terms: Modeled using symmetry-conserving rational forms involving $\cos(n\theta)$, where $n$ is fixed by the lattice's rotational symmetry (e.g., $n=4$ for tetragonal, $n=6$ for hexagonal). This preserves analyticity at $q \to 0$ and captures leading angular dependencies.

3. Parameterization: The coefficients for these functions are obtained by fitting the calculated first-principles inverse dielectric function. The study applies this to seven 2D metals (Na, Mg, K, MgB$_2$, AlB$_2$, NbSe$_2$, Ca$_2$N) using RPA data generated via Quantum Espresso and the Yambo code.

4. Correlation Energy Calculation: The RPA correlation energy ($E_c$) is evaluated by integrating the analytical MPA(q) representation along the imaginary frequency axis. This allows for evaluation on flexible and dense momentum grids, including systematic refinement of the $q \to 0$ limit.

Key Results

• Compact Representation: A small number of dispersive plasmonic modes ($n_p$) suffices to accurately describe the dielectric response across the full Brillouin zone for diverse electronic regimes. Simple metals (Na, Mg, K) required only 1–2 poles, while more complex systems with multiband carriers (MgB$_2$, AlB$_2$, NbSe$_2$, Ca$_2$N) required between 4 and 6 poles.
• Spectral Accuracy: The MPA(q) model successfully reproduces the main spectral features, including plasmonic dispersions and damping, for both the direct and inverse dielectric functions. The model captures the expected lattice symmetries and quantifies the degree of anisotropy.
• Correlation Energy Precision: The deviation of the correlation energy computed via MPA(q) compared to numerical first-principles evaluation on a $24 \times 24 \times 1$ grid was found to be less than 1.3 meV for all systems studied.
• Improved Convergence: The analytical nature of the model allows for the evaluation of $E_c$ on denser momentum grids than computationally feasible with direct numerical integration. Crucially, the parametrization improves the description of the $q \to 0$ limit, significantly accelerating the convergence of the correlation energy toward the dense-grid limit.

Significance and Claims
The paper claims to establish a direct bridge between large-scale ab initio calculations and analytical models of screening. By replacing large numerical datasets with compact analytical forms, the MPA(q) framework enables:

1. Efficient Evaluation: Rapid calculation of quantities involving dynamical screening, such as spectral features and correlation energies, without the computational prohibitions of dense grid sampling.
2. Systematic Refinement: Controlled refinement of momentum sampling and the $q \to 0$ limit, which is essential for metallic and low-dimensional systems where dielectric responses exhibit strong nonlocal behavior.
3. Foundation for Many-Body Methods: The authors posit that this approach facilitates efficient implementations of advanced many-body approaches requiring momentum- and frequency-dependent dielectric functions, such as GW/BSE calculations and coupled quasiparticle frameworks.

The work demonstrates that despite the diversity of electronic regimes in 2D metals, a small set of dispersive plasmonic modes can accurately capture the full momentum and frequency-dependent dielectric response, offering a practical route for evaluating screening observables within RPA and beyond.