🔬 mesoscale physics

Dielectric function in WSe2

This paper presents a Hartree-Fock numerical method to compute the static dielectric function of a fully spin-polarized Wigner crystal in monolayer WSe2, offering a theoretical tool for investigating screening and interaction effects in low-density two-dimensional systems.

Original authors: Tiberius O. Cheche, Yia-Chung Chang

Published 2026-07-07

📖 4 min read☕ Coffee break read

CC BY 4.0

Original authors: Tiberius O. Cheche, Yia-Chung Chang

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

The Big Picture: Electrons Turning into a Crystal

Imagine a crowded dance floor where everyone is moving randomly. This is how electrons usually behave in a metal or a semiconductor—they are a chaotic "gas" of particles zipping around.

However, the paper explores a special situation where the electrons are very far apart (low density). In this scenario, the electrons stop dancing randomly and decide to stand still in a perfect, organized pattern, like soldiers in a formation. This is called a Wigner crystal.

The authors of this paper created a computer program to figure out exactly what this "electron crystal" looks like inside a specific material called monolayer WSe2 (a super-thin sheet of a semiconductor). They wanted to know two things:

What are the energy levels of these electrons in this crystal? (The "Band Structure")
How does this crystal react if you try to push or pull on it? (The "Dielectric Function")

The Method: A Self-Checking Simulation

To solve this, the authors used a method called Hartree-Fock. Think of this as a game of "hot potato" played by a computer to find the perfect arrangement of electrons.

The Guess: The computer starts with a guess of where the electrons are (like a blurry photo of a crystal).
The Calculation: It calculates how the electrons push and pull on each other based on that guess.
The Update: It updates the positions of the electrons to see if they can lower their energy.
The Loop: It repeats this over and over until the picture stops changing and settles into a stable, perfect crystal shape.

In this specific study, they assumed all the electrons were spinning in the same direction (fully spin-polarized), which simplifies the math and mimics a specific experimental setup.

The Results: What Happens as the Crowd Thins Out?

The researchers ran their simulation for different "crowd densities," represented by a number called $r_s$ .

Low $r_s$ (Crowded): The electrons are closer together. The "crystal" is weak and fuzzy. The electrons are still moving around a bit, and the pattern isn't very sharp.
High $r_s$ (Sparse): The electrons are very far apart. The "crystal" becomes very strong and sharp. The electrons lock tightly into their spots, like ice crystals forming in very cold water.

Key Findings:

The Pattern: As the electrons get further apart, the charge density (the "weight" of the electrons) becomes highly concentrated in specific spots, creating a clear triangular lattice pattern.
The Energy: The energy levels of the electrons change in a way that confirms they are trapped in this crystal structure.
The Shielding (Dielectric Function): This is a fancy way of asking, "If I poke the crystal, how well does it shield the poke?"
- In a dense crowd, the electrons are good at shielding each other; they rearrange quickly to cancel out the push.
- In the sparse Wigner crystal (high $r_s$ ), the electrons are stuck in their spots. They can't move much to shield the poke. The paper found that as the crystal gets more "frozen," its ability to shield external forces drops significantly. It becomes less like a fluid that flows to protect itself and more like a rigid wall.

Why This Matters (According to the Paper)

The paper doesn't claim to have built a new device or cured a disease. Instead, it provides a theoretical tool.

Think of the paper as a detailed blueprint or a weather forecast for a specific type of electronic storm. The authors built a mathematical model that predicts exactly how these electron crystals behave in WSe2.

They suggest that future experiments could use this blueprint to check if they have actually created a Wigner crystal. For example, if scientists shine specific colors of light on a WSe2 sheet and see the exact energy patterns the computer predicted, it would be proof that they successfully created this electron crystal.

Summary

The Problem: How do electrons behave when they are forced to stand still in a perfect grid?
The Solution: A computer simulation that calculates the energy and shielding properties of this grid in a specific material (WSe2).
The Discovery: As the electrons get further apart, they form a sharper, more rigid crystal, but they become worse at shielding each other from outside forces.
The Goal: To give experimental scientists a "cheat sheet" so they know what to look for when they try to create these crystals in a lab.

Technical Summary: Dielectric function in WSe2 monolayer Wigner crystal

Problem Statement
The paper addresses the theoretical characterization of the Wigner-crystal state in two-dimensional (2D) transition-metal dichalcogenide (TMD) monolayers, specifically focusing on monolayer WSe2. While the existence of Wigner crystallization in low-density electron gases is well-established theoretically and recently observed experimentally in TMD heterostructures, there is a need for a microscopic framework to compute the electronic band structure and static dielectric response of these crystalline states. The challenge lies in accurately modeling the interplay between kinetic energy, Coulomb interactions (specifically the Keldysh interaction relevant to 2D materials), and exchange effects in a broken-symmetry state where continuous translational symmetry is lost.

Methodology
The authors develop a self-consistent Hartree-Fock (HF) numerical method to compute the band structure and static dielectric function of a fully spin-polarized 2D Wigner crystal. The key methodological components include:

Hartree-Fock Formulation: The calculation is performed within a restricted HF framework where the ground state is assumed to be fully spin-polarized. The single-particle orbitals are expanded in a Bloch form using a plane-wave basis $|k + G\rangle$ , where $k$ belongs to the first Brillouin zone of the emergent triangular Wigner-crystal lattice and $G$ represents reciprocal lattice vectors.
Interaction Potential: The effective electrostatic interaction between charge carriers is modeled using the Keldysh interaction, which accounts for the dielectric screening of the surrounding medium and the 2D nature of the material.
Self-Consistent Loop: The algorithm initializes with a periodic Gaussian-like charge density. It iteratively updates the Hartree potential (derived from density harmonics) and the Fock exchange matrix (derived from occupied eigenvectors) until convergence. The density harmonics are reconstructed from the updated eigenvectors to close the loop.
Numerical Implementation:
- The Brillouin zone sampling utilizes an irreducible $\Gamma-K-M$ wedge for the target momenta, while the exchange term is evaluated via a full-zone summation reconstructed from symmetry images.
- A circular cutoff in reciprocal space selects the plane-wave basis, preserving the $D_6$ point-group symmetry of the triangular lattice.
- Singularities in the Coulomb/Keldysh potential at $q=0$ are handled by analytic averaging over a small circular cell in reciprocal space.
Dielectric Response Calculation: The static dielectric function $\epsilon(q, 0)$ is computed using a Lindhard-type polarizability formula adapted for the insulating, folded-Brillouin-zone nature of the Wigner crystal. This involves summing over interband transitions between occupied and empty HF bands, utilizing the converged eigenvectors and eigenvalues.

Key Contributions

Development of a Specialized Code: The authors present a numerical code capable of solving the self-consistent HF equations for a 2D Wigner crystal with Keldysh interactions, specifically tailored for TMD monolayers.
Microscopic Dielectric Response: The work provides a method to derive the static dielectric function directly from the HF band structure and eigenfunctions of the crystalline state, rather than relying on classical phonon models or uniform electron gas approximations.
Characterization of Crystalline Order: The study introduces and utilizes specific metrics to quantify the degree of crystallization:
- The first-star order parameter ( $M_1$ ), measuring the amplitude of the fundamental density modulation.
- The density-contrast ratio ( $C_n$ ), quantifying the real-space separation between density maxima and minima.
- The evolution of the band structure at high-symmetry points ( $\Gamma, K, M$ ) as a spectral signature of the triangular lattice symmetry.

Results
The authors applied their method to monolayer WSe2 across a range of density parameters ( $r_s$ ), representing the transition from weakly crystalline to strongly localized states:

Band Structure Evolution: As $r_s$ increases (density decreases), the system transitions from a weakly modulated state to a robust Wigner crystal. At higher $r_s$ , the band structure exhibits clear degeneracy patterns consistent with the hexagonal Brillouin zone symmetry, and the energy bands evolve smoothly along the $\Gamma-K-M$ path.
Density Modulation: The order parameter $M_1$ increases with $r_s$ (e.g., from 0.100 at $r_s=4$ to 0.586 at $r_s=50$ ), indicating stronger charge localization. The density-contrast ratio $C_n$ also increases significantly, confirming the formation of distinct charge density peaks at lattice sites.
Dielectric Function: The static dielectric function $\epsilon_{HF}(q, 0)$ shows a systematic reduction in its peak value as $r_s$ increases. This reflects a decrease in electronic polarizability; as electrons become more localized, they become less efficient at screening external perturbations, causing the screened interaction to approach the bare Keldysh interaction.
Convergence: The self-consistent solutions demonstrate stability, with the final density profiles relaxing from the initial Gaussian guess to a state determined by the balance of kinetic, Hartree, and exchange energies.

Significance and Claims
The paper positions itself as a theoretical tool for investigating screening, band-structure reconstruction, and interaction effects in low-density 2D systems. The authors claim that their method provides a consistent microscopic picture of the evolution from a weakly crystalline solution to a localized Wigner-crystal state.

The significance of the work is framed in the context of recent experimental advances in atomically thin semiconductors. The authors suggest that their calculated spin-conserving, frozen-spin/valley sector could be tested via time-resolved, helicity-resolved optical spectroscopy on gate-tunable WSe2 monolayers, potentially driven into a polarized state by optical pumping. However, they explicitly state that their results serve as "phase-sensitive indicators" of the crystalline character rather than a thermodynamic proof of the liquid-crystal transition, noting that a definitive identification of the transition would require energy comparisons or independent phase-stability analyses. The work connects conventional low-density 2D Wigner crystal theory with recent optical studies of charge order in TMDs, offering a framework adapted to the specific Keldysh interaction and optical response of monolayer TMDs.