🔬 materials science

Unified description of exciton, phonon and plasmon dispersions in 2D materials from an optical conductivity approximation

This paper introduces the optical conductivity approximation (OCA), a unified framework that accurately describes the linear low-momentum dispersions of excitons, phonons, and plasmons in 2D materials by leveraging q=0 optical conductivity to reveal that these behaviors stem primarily from the 2D macroscopic Coulomb interaction rather than band structure details.

Original authors: Alberto Guandalini, Andrea Ferretti, Daniele Varsano, Paolo Barone, Francesco Mauri

Published 2026-09-25

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Alberto Guandalini, Andrea Ferretti, Daniele Varsano, Paolo Barone, Francesco Mauri

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

In the microscopic world of materials science, scientists study how energy moves through solids, much like watching ripples spread across a pond. When a material is squeezed into a sheet just one atom thick, the rules change. In these two-dimensional layers, the forces that hold atoms together and allow electrons to move behave differently than they do in bulk materials. This shift creates unique ripples of energy called excitations. Some of these ripples are vibrations of the atoms themselves, known as phonons. Others are pairs of electrons and the holes they leave behind, called excitons, or collective waves of electrons known as plasmons. Understanding how the speed and energy of these ripples change as they travel across the material is crucial for designing faster electronics and more efficient solar cells. However, measuring these changes is difficult because the forces involved are long-range and complex, requiring immense computational power to track every tiny step of the journey.

A team of researchers from Italy has developed a new way to map these energy ripples without needing to calculate every single step. They focused on materials like hexagonal boron nitride and graphene, which are famous for their strength and electrical properties. The scientists discovered that the way these energy ripples speed up or slow down is not primarily driven by the specific arrangement of the atoms or the complex paths electrons take through the material's internal structure. Instead, the shape of the energy curve is dictated almost entirely by the unique way electric forces behave in a flat, two-dimensional world. In three-dimensional space, electric forces weaken in a specific way as you move away from a source, but in a flat sheet, they weaken more slowly. This difference creates a distinct pattern in how the energy ripples travel.

To prove this, the researchers introduced a method they call the optical conductivity approximation. Imagine trying to predict how a wave moves across a lake. Usually, you would need to measure the water depth, the wind speed, and the shape of the shore at every single point. This new method suggests that if you know how the water responds to a gentle push at a single spot, you can accurately predict how the wave will behave across the entire lake, provided you account for the shape of the lake itself. In this case, the "lake" is the two-dimensional material, and the "push" is the interaction between light and the material's electrons. By using a simplified measurement taken at a single point where the momentum is zero, the team was able to reconstruct the full energy map of the ripples.

When they applied this method to hexagonal boron nitride, a material often used as an insulator in electronics, the results were striking. The team accurately reproduced the energy and intensity of the longitudinal optical phonons, which are specific vibrations of the atoms, as well as the bright excitons, which are the electron-hole pairs that interact strongly with light. They also successfully mapped the behavior of the plasmons in graphene, a single layer of carbon atoms. In every case, the method showed that the linear increase in energy as the ripples moved was a direct consequence of the two-dimensional electric forces, not a complex variation in the material's internal band structure. The researchers found that the energy of these ripples rises in a straight line as they gain momentum, a behavior that their new framework captured perfectly using only the data from the zero-momentum point.

The study also clarified why previous attempts to model these materials were so computationally expensive. To get accurate results before, scientists had to calculate the response of the material at a dense grid of momentum points, a process that required massive computing resources. The new approach shows that this heavy lifting is unnecessary for understanding the low-momentum behavior. By focusing on the optical conductivity, which changes very little as the momentum shifts, the researchers could bypass the need for dense calculations. This not only saves time and energy but also provides a clearer physical picture of what is happening. It confirms that the unique behavior of these two-dimensional materials is a fundamental result of their geometry and the nature of electric forces, rather than a quirk of their specific atomic makeup.

Furthermore, the researchers demonstrated that this method works in reverse. Just as they could predict the energy ripples from the optical data, they showed that one could take experimental data from electron energy-loss spectroscopy and extract the full optical properties of the material. This offers a powerful tool for experimentalists who want to understand the optical behavior of a material without needing large, perfect crystals or high-energy light sources that are difficult to produce. The method was tested on both single-layer and double-layer sheets of hexagonal boron nitride, showing that it holds up even as the material gets slightly thicker. The findings suggest that this unified framework can be applied to a wide range of semiconducting and semi-metallic systems, providing a simpler, more efficient way to understand the fundamental physics of the two-dimensional materials that are shaping the future of technology.

Technical Summary: Unified Description of Exciton, Phonon, and Plasmon Dispersions in 2D Materials

Problem Statement
In two-dimensional (2D) materials, the modified long-range behavior of the Coulomb interaction ( $v_{2D}(q) \sim 2\pi/q$ ) fundamentally alters the dispersion of longitudinal elementary excitations, including excitons, phonons, and plasmons. Unlike in bulk 3D crystals where longitudinal-transverse (LO/TO) splitting occurs at the $\Gamma$ point, 2D materials exhibit a collapse of this splitting at $q=0$ . Consequently, longitudinal modes display a non-analytic linear dispersion ( $\propto q$ ), while transverse modes retain a parabolic dispersion ( $\propto q^2$ ).

Accurately describing these dispersions typically requires calculating response functions over dense momentum grids, which is computationally prohibitive, particularly for many-body perturbation theory methods like the Bethe-Salpeter Equation (BSE) and GW approximations. Furthermore, a unified theoretical framework that simultaneously accounts for the distinct behaviors of electronic (excitons, plasmons) and ionic (phonons) longitudinal excitations under the influence of the 2D macroscopic Coulomb field has been lacking.

Methodology: The Optical Conductivity Approximation (OCA)
The authors introduce the Optical Conductivity Approximation (OCA), a unified framework that links the finite-momentum ( $q$ ) longitudinal density-density response function, $\chi(q, \omega)$ , to the optical conductivity at zero momentum, $\sigma(\omega) \equiv \sigma(q=0, \omega)$ .

The derivation relies on the continuity equation and the relationship between the induced current and the total macroscopic electric field. The key theoretical insight is that while the reducible density-density response $\chi(q, \omega)$ exhibits singular $q$ -dependence due to the macroscopic Hartree potential, the conductivity $\sigma(q, \omega)$ is inherently free from this macroscopic Coulomb divergence. Consequently, $\sigma(q, \omega)$ displays a much milder dependence on momentum near the $\Gamma$ point.

The OCA posits that $\sigma(q, \omega) \approx \sigma(0, \omega)$ . By substituting this approximation into the exact relation derived from the continuity equation:
$\chi(q, \omega) = \frac{\sigma(q, \omega)|q|^2}{-i\omega - v_{2D}(q)\sigma(q, \omega)|q|^2}$
the authors obtain the OCA expression:
$\chi_{OCA}(q, \omega) = \frac{\sigma(\omega)|q|^2}{-i\omega - v_{2D}(q)\sigma(\omega)|q|^2}$
where $v_{2D}(q) = 2\pi/q$ . This formulation allows the calculation of the finite- $q$ loss function (proportional to $-\text{Im}[\chi]$ ) using only the $q=0$ optical conductivity. The conductivity is decomposed into electronic ( $\sigma_{el}$ ) and vibrational ( $\sigma_{ph}$ ) contributions, modeled as sums of Lorentzian resonances.

Key Results
The authors applied the OCA to monolayer and bilayer hexagonal boron nitride (hBN) and graphene, comparing the results against ab-initio calculations (BSE for excitons, DFPT for phonons, and GW-BSE for plasmons) performed on finite momentum grids.

hBN Phonons and Excitons: The OCA accurately reproduces the energy and intensity dispersions of the longitudinal-optical (LO) phonon and the low-energy bright excitons (A and B branches) in monolayer hBN. The linear dispersion of these modes and the rapid increase in the intensity ratio between the B and A branches are captured without calculating $\sigma$ at finite $q$ .
Graphene Plasmons: The framework successfully reproduces the linear dispersion of the graphene $\pi$ -plasmon. The results confirm that this dispersion is driven by the $q$ -dependence of the 2D Coulomb interaction acting on the $q=0$ optical conductivity, rather than by intrinsic variations in the band structure dispersion.
Model Analysis: By approximating $\sigma(\omega)$ with a single Lorentzian plus a non-absorptive background, the authors derived analytical expressions for the momentum-dependent energy $\Omega(q)$ and residue $R(q)$ . They demonstrated that the linear dispersion coefficient is directly proportional to the residue of the optical conductivity at $q=0$ .
Limitations: The OCA is exact at $q=0$ and remains accurate for low momentum transfers. However, it cannot describe dark (transverse) excitations, such as dark excitons or transverse optical (TO) phonons, as these do not couple to the macroscopic longitudinal field. Additionally, the method does not capture the onset of the excitation spectrum in graphene, which is determined by the linear band structure of the Dirac cone (a purely band-structure effect), whereas the plasmon dispersion is a collective effect.

Significance and Claims
The paper claims to provide a unified description of 2D longitudinal excitations, demonstrating that their sharp linear dispersions are primarily a consequence of the 2D macroscopic Coulomb interaction form ( $v_{2D} \sim 1/q$ ) rather than the dispersion of the underlying band structure.

The significance of the work is twofold:

Computational Efficiency: The OCA offers an efficient scheme to generate accurate low-momentum Electron Energy-Loss Spectroscopy (EELS) data without the need for expensive calculations on dense momentum grids. It effectively interpolates between ab-initio points using only a single $q=0$ calculation.
Inverse Methodology: The framework serves as an inverse method, allowing researchers to extract the full complex optical conductivity from low-momentum EELS spectra. This enables access to optical properties over broad energy ranges, potentially circumventing limitations of standard optical measurements regarding high-energy reach and sample purity requirements.

The authors note that the scheme is applicable to semiconducting, semimetallic, and metallic systems, providing a consistent physical picture for the dispersion of plasmons, excitons, and phonons in reduced dimensions.

Technical Summary: Unified Description of Exciton, Phonon, and Plasmon Dispersions in 2D Materials

More like this