🔬 mesoscale physics

Accurate and efficient simulation of photoemission spectroscopy via Kohn-Sham scattering states

This paper introduces an efficient, first-principles framework that computes photoelectron states as Kohn-Sham scattering solutions to enable accurate, transparent, and widely compatible simulations of angle-resolved photoemission spectroscopy (ARPES), as demonstrated by excellent agreement with experimental data for graphene and WSe $_2$ .

Original authors: Gian Parusa, Sotirios Fragkos, Samuel Beaulieu, Michael Schüler

Published 2026-02-02

📖 4 min read☕ Coffee break read

CC BY 4.0

Original authors: Gian Parusa, Sotirios Fragkos, Samuel Beaulieu, Michael Schüler

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

Imagine you are trying to take a perfect photograph of a bustling city at night. You want to see not just where the buildings are, but exactly how the light reflects off their windows, how shadows fall, and how the city's layout changes the way the light looks from your specific camera angle.

In the world of physics, Angle-Resolved Photoemission Spectroscopy (ARPES) is like that camera. Scientists shoot high-energy light (photons) at a material, knocking electrons out of it. By measuring the speed and direction of these flying electrons, they can map out the material's internal "city plan"—its electronic structure.

However, there is a catch. The photo you get isn't just a direct picture of the city; it's a complex mix of the city's layout and how the light bounces around, hits walls, and interferes with itself before reaching your camera. For a long time, simulating this process on a computer was like trying to solve a massive, tangled knot of string. Existing methods were either too rigid (only working with specific types of materials) or too slow to be practical for modern, complex experiments.

The New "Camera Lens" Approach
This paper introduces a new, efficient way to simulate these photos. The authors, Gian Parusa and his team, developed a method that treats the escaping electrons like waves hitting a shore.

Instead of using complicated, specialized software that only works for certain materials, they built a tool that works with the standard "blueprints" (computer codes) most scientists already use. Their method solves a specific math problem (the Kohn-Sham equation) with special rules at the edges (boundary conditions) that tell the computer: "Imagine these electrons are running away from the material into empty space."

Why is this better?
Think of it like this:

Old Methods: Like trying to build a house by hand-crafting every single brick from scratch. It works, but it's slow and you can't easily change the design later.
This New Method: Like using a high-quality, pre-fabricated wall system that fits into any standard house plan. It's fast, flexible, and lets you see exactly how the light hits the walls before you even build the house.

The "Ghost" in the Machine: Pseudopotentials
One of the biggest hurdles in these simulations is dealing with the heavy atomic cores (the nucleus and inner electrons). To save computer power, scientists often use "pseudopotentials"—which are like simplified masks that represent the heavy atoms without calculating every tiny detail.

The team tested if these "masks" were accurate enough to predict how high-speed electrons scatter. They found that:

Simple masks work well for many materials, provided the mask is high-quality.
However, for heavy atoms (like Tungsten in WSe2), the mask needs to include "deep secrets" (semicore states). If you leave these out, the simulation gets the "shadows" wrong, leading to a distorted photo. It's like wearing a mask that covers your eyes but forgets to cover your ears; you can see, but you miss crucial sound cues that change how you react to the world.

The Proof: Graphene and WSe2
To prove their method works, they simulated two materials:

Graphene (a single layer of carbon): They predicted how the light patterns (called circular dichroism) would look. Their simulation matched real-world experiments perfectly, even predicting subtle "nodal lines" (places where the signal disappears) that other methods missed.
WSe2 (a bulk crystal): They showed that including those "deep secrets" (semicore states) in their masks was essential to get the correct patterns. Without them, the simulation looked like a blurry, wrong version of the real experiment.

The Bottom Line
This paper doesn't just offer a faster way to do math; it offers a clearer window into how light and matter interact. By calculating the exact "flight path" of the escaping electrons, scientists can now:

Understand why certain patterns appear in their experiments.
Distinguish between the material's true nature and the "optical illusions" caused by the measurement process.
Study complex materials and even materials in motion (like those being pumped with lasers) using standard, widely available computer tools.

In short, they have given researchers a sharper, more flexible lens to see the invisible electronic world inside solids.

Technical Summary: Accurate and Efficient Simulation of Photoemission Spectroscopy via Kohn-Sham Scattering States

Problem Statement
Angle-resolved photoemission spectroscopy (ARPES) is a primary tool for mapping the electronic structure of solids, yet a significant gap exists between intrinsic material properties (initial states) and experimentally observed spectra. The observed intensity depends not only on the initial Bloch states but also critically on the final photoelectron states, which are influenced by experimental geometry, light polarization, intra- and inter-atomic scattering, and multiple scattering within the crystal potential. While the "one-step" theory of photoemission provides a rigorous framework, existing implementations—particularly the mature Korringa-Kohn–Rostoker (KKR) multiple-scattering method—are often tied to specialized electronic-structure formalisms. This limits their integration with widely used plane-wave or real-space density functional theory (DFT) codes, hindering the application of advanced exchange-correlation functionals and modern workflows. Furthermore, the accuracy of pseudopotentials (PP), essential for efficient DFT calculations, in describing high-energy photoelectron states remains an open question.

Methodology
The authors introduce a first-principles framework that computes final photoelectron states by directly solving the Kohn-Sham (KS) equation with scattering boundary conditions.

Formalism: The approach assumes a two-dimensional periodic system and employs an in-plane Laue representation for wavefunctions. The KS equation is transformed from an eigenvalue problem (for bound states) into a linear equation for scattering states by imposing asymptotic boundary conditions: outgoing waves at $z \to +\infty$ and incoming/standing waves at $z \to -\infty$ .
Equivalence to Lippmann-Schwinger: The authors demonstrate that solving the KS equation with these boundary conditions is formally equivalent to the Lippmann-Schwinger (LS) formulation. However, the direct KS approach allows for the use of highly effective preconditioning strategies derived from iterative diagonalization, resulting in significantly faster convergence compared to solving the LS integral equation directly.
Implementation: The method is implemented within standard plane-wave and real-space DFT codes. It utilizes discrete $z$ -grids with higher-order discretization (up to $\mathcal{O}(h^5)$ ) and iterative linear solvers.
Pseudopotential Analysis: A key component of the study is the systematic evaluation of pseudopotentials. The authors compare all-electron (AE) calculations with generalized norm-conserving pseudopotentials (PP), specifically investigating the necessity of including semicore states and non-local potentials to accurately reproduce high-energy scattering properties.

Key Contributions

Efficient Framework: The paper presents a method that is fully compatible with standard DFT workflows (plane-wave and real-space), enabling the use of advanced functionals and large-scale simulations that were previously difficult with specialized one-step codes.
Explicit Wavefunctions: Unlike methods that yield only intensities, this approach provides explicit photoelectron wavefunctions. This allows for a transparent analysis of matrix-element effects, multiple scattering, and the disentangling of intrinsic initial-state properties from extrinsic final-state effects.
Pseudopotential Validation: The study establishes that high-quality pseudopotentials, particularly those with optimized non-local parts and the inclusion of semicore states, can accurately describe photoelectron states in the VUV to XUV regime, making the method practical for routine use.

Results
The framework was validated through simulations of circular-dichroism ARPES (CD-ARPES) for two systems:

Monolayer Graphene: Theoretical results for energy-integrated circular dichroism in the angular distribution (CDAD) were compared against experimental data (Ref. [61]). The method successfully reproduced nodal lines, overall CDAD textures, and the sign reversal of CDAD at 75 eV. Comparisons between all-electron and pseudopotential methods showed excellent agreement up to 100 eV, provided the pseudopotential included non-local contributions. The non-local part was found crucial for reproducing scattering properties, though its effect varied with photon energy and angular momentum channels (e.g., negligible for $d$ -like partial waves in carbon at 52 eV).
Bulk 2H-WSe $_2$ : New experimental data acquired at the CELIA facility were compared with theoretical simulations. The inclusion of semicore states ( $5s, 5p$ for W and $3d$ for Se) in the pseudopotential was shown to be essential. Simulations with semicore states achieved quantitative agreement with experimental valley-integrated CDAD across all $K$ -points, whereas simulations without semicore states deviated significantly, particularly at $K$ -points 1 and 2. This confirmed that accurate scattering phases around atomic sites are necessary to capture the interferometric nature of photoemission.

Significance and Claims
The authors claim that this work establishes a robust and accessible route toward quantitative ARPES modeling. By bridging the gap between widely used ab initio platforms and photoemission theory, the method opens the door to systematic studies of orbital textures, many-body effects, and nonequilibrium phenomena.

Versatility: The explicit calculation of matrix elements allows for post-processing of arbitrary experimental geometries and light polarizations.
Extensibility: The method can be combined with Green's function techniques (e.g., $GW$) to treat correlation effects beyond standard DFT. It is also positioned as a foundation for simulating time-resolved ARPES (trARPES) in pump-probe scenarios, offering a potential alternative to KKR-based trARPES implementations which are currently limited to standard DFT.
Limitations: The authors modestly note that treating surfaces with very large unit cells may be computationally challenging, and truly semi-infinite systems might still require Green's function formulations like KKR. Additionally, at very high photon energies, core state contributions may necessitate the reconstruction of all-electron wavefunctions from pseudowavefunctions.

In summary, the paper provides a computationally efficient, formally equivalent alternative to the Lippmann-Schwinger formalism that integrates seamlessly with modern DFT codes, validated by high-accuracy agreement with experimental CD-ARPES data.

Technical Summary: Accurate and Efficient Simulation of Photoemission Spectroscopy via Kohn-Sham Scattering States

More like this