Pith. sign in

REVIEW 3 major objections 5 minor 29 references

Photoemission from Semi-Infinite Crystals: An \emph{Ab Initio} Scattering-State Approach

T0 review · 3 major / 5 minor · reviewed 2026-08-01 · deepseek-v4-flash

Pith's one-line read A parameter-free Wannier + Green-function embedding constructs continuum-normalized time-reversed LEED final states for semi-infinite crystal-vacuum interfaces, and the resulting Ag(111) photoemission observables land on the experimental sc

desk verdict A real advance: parameter-free semi-infinite scattering states with absolute-scale Ag(111) agreement; the occupied-state truncation needs convergence data before the quantitative claim fully lands. read the letter →

arxiv 2607.19551 v1 pith:EUIE6JC7 submitted 2026-07-21 cond-mat.mtrl-sci physics.acc-ph

classification cond-mat.mtrl-sciphysics.acc-ph
keywords photoemissionsemi-infinitecrystaltime-reversedLEEDWannierfunctionsGreen-functionembeddingquantumefficiencymeantransverseenergyAg(111)
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

Photoemission is an escape problem, but standard first-principles calculations enclose the electron in a periodic box. This paper removes that artificial boundary by building scattering states for a semi-infinite crystal-vacuum interface directly from Wannier Hamiltonians and Green-function embedding. For each transverse momentum, the problem becomes a one-dimensional open-boundary scattering problem whose solution is a continuum-normalized time-reversed LEED final state with microscopic quasiparticle attenuation. Applied to Ag(111), the method predicts absolute quantum efficiency, the vectorial photoelectric effect, and mean transverse energy on the experimental scale, without rescaling or empirical damping. This matters because it offers a route to fully ab initio one-step photoemission and other open-boundary observables for realistic materials.

What carries the argument

The load-bearing mechanism is the hybrid Wannier representation combined with Green-function embedding. Fourier transforming maximally localized Wannier functions over lattice vectors parallel to the surface gives a short-ranged effective one-dimensional chain along the surface normal for each transverse momentum. The semi-infinite crystal and vacuum half-spaces are integrated out via retarded surface Green's functions, leaving a finite interface problem; solving it with a source term that fixes the incoming Bloch mode yields the complete scattering wave function — the continuum-normalized time-reversed LEED final state. Continuum normalization and microscopic attenuation are built in by the

What would settle it

Repeat the Ag(111) calculation with two, four, and eight additional bulk principal layers in the occupied-state construction and recompute the absolute quantum efficiency and mean transverse energy at a fixed excess photon energy; if the results shift by more than the experimental scatter, the open-boundary final states are not the sole determinant of the reported observables.

Watch

Extended reading notes

Core claim

The central claim is that a parameter-free hybrid-Wannier plus Green-function embedding construction yields the full scattering wave function of the semi-infinite crystal-vacuum system, not merely transmission probabilities. The final states are continuum-normalized time-reversed LEED states; adding a first-principles complex self-energy to the Wannier Hamiltonian produces inelastic attenuation inside the crystal while leaving the vacuum channel propagating. Using these states in one-step photoemission matrix elements, with the Fresnel-transmitted optical field and its surface divergence included, the method reproduces the absolute scale of Ag(111) quantum efficiency, the vectorial photoelec

Load-bearing premise

The occupied initial states are obtained by diagonalizing a finite Hamiltonian with an artificial termination, and the paper gives no convergence test against the number of appended bulk layers; if that termination shifts surface-state energies or occupations relative to the Fermi level, the absolute emission yield and the surface/bulk balance would change.

Editorial extensions

If this is right

  • One-step photoemission calculations can proceed without a periodic supercell, an empirical optical potential, or a fitted inelastic mean free path.
  • Momentum-resolved observables and matrix elements that depend on the full quantum state become directly accessible, not just total transmission probabilities.
  • The same scattering-state construction applies to any lattice-matched or commensurate interface for which a localized Wannier representation exists, covering transport and interface spectroscopy.
  • The inelastic mean free path emerges from the first-principles self-energy, so it can be extracted as a function of energy, momentum, and temperature rather than supplied as an input.
  • The symmetrized optical perturbation, including the surface divergence term, gives a first-principles account of polarization-dependent emission such as the vectorial photoelectric effect.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • A natural test is to compare the energy-dependent inelastic mean free path extracted from these scattering states against independent attenuation-length measurements; the paper treats it as an output rather than an input.
  • The finite-slab treatment of occupied states is the least protected part of the calculation; systematic convergence checks against the number of appended bulk layers would sharpen every absolute number reported.
  • The interface optical field is modeled with a delta-function sheet charge; probing s- versus p-polarized emission over a broader frequency range would test that approximation independently of the electronic-structure machinery.
  • Because explicit wave functions are produced, the same embedding could be used for higher-order processes (for example, phonon-assisted or two-step emission) without changing the boundary-condition construction.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 5 minor

Summary. The paper presents an ab initio one-step photoemission formalism for semi-infinite crystal–vacuum interfaces. It combines hybrid Wannier functions, Green-function embedding of semi-infinite leads, and first-principles G0W0/Fan–Migdal self-energies to construct time-reversed LEED final states with microscopically computed attenuation. For Ag(111), the method is used to predict absolute quantum efficiency, the vectorial photoelectric effect, and mean transverse energy; the results are compared directly with experiment on an absolute scale without rescaling. The central claim is that this procedure is parameter-free and removes the need for a periodic box, an empirical optical potential, or an externally imposed inelastic mean free path.

Significance. If the construction is sound, this is a significant advance for first-principles photoemission: it supplies explicit open-boundary scattering wave functions rather than transmission probabilities, and it replaces phenomenological damping with computed self-energies. The absolute-scale agreement for Ag(111) QE and MTE is a nontrivial, falsifiable test that goes beyond the usual band-structure comparisons. The method's applicability to other interfaces is also plausible. However, the validation depends on two under-documented elements — the finite-slab construction of the occupied initial states and the inputs entering the Fresnel optical field — so the significance is currently conditional.

major comments (3)
  1. [End Matter, Occupied-state construction] The occupied initial states, which dominate the near-threshold emission and the surface/bulk balance, are obtained by diagonalizing a finite Hamiltonian with 'many bulk principal layers appended,' but no layer count, no energy-shift data, and no observable convergence tests are reported. The claim that the states 'were converged with respect to the artificial termination' is therefore unsupported. This is load-bearing: the paper itself attributes the remaining QE shape discrepancies to 'the relative weight of surface-state and bulk-state emission' (Results, QE paragraph). Please provide convergence of the Shockley surface-state energy and occupation, and of the computed QE and MTE, as a function of the number of appended bulk layers. Also specify how the finite set of occupied bulk states samples the near-surface continuum and whether the thresholds in Fig. 3 shift with termination.
  2. [Methods, Optical Transition Matrix Elements and Eq. (19)] The abstract and conclusion call the approach 'parameter-free,' but the optical field entering Eq. (19) is determined by the complex dielectric function of Ag, described as the sum of an independent-particle interband response and a metallic intraband Drude response 'following Ref. [19].' It is not stated whether the Drude plasma frequency, damping rate, or any other dielectric input is computed from first principles or taken from experiment/fitted. Since the absolute QE and the vectorial effect depend sensitively on the Fresnel-transmitted field and on the delta-function surface term, any empirical or fitted dielectric parameter would directly affect the headline claim. Please list all parameters entering ε and state explicitly whether any is adjusted to improve agreement with the photoemission data.
  3. [End Matter, Hybrid Wannier construction] The vacuum sector is an essential part of the open-boundary construction, but the regular close-packed vacuum Wannier lattice is said to be 'developed in Ref. [14],' which appears to be an unpublished companion paper. The brief description in the End Matter (Gaussian trial orbitals projected onto the empty-state manifold) is not sufficient for a reader to reproduce or independently assess the method. Please either expand the description so that the vacuum Wannier construction is self-contained, or provide a published or otherwise available reference with full computational details.
minor comments (5)
  1. [Title] The title contains a typo: 'AnAb Initio' should be 'An Ab Initio.'
  2. [Fig. 4] The y-axis label 'Absolute QE 1e 5' is ambiguous; it appears to indicate a multiplicative factor of 10^5, but the figure would be clearer with explicit units (e.g., 'QE (×10^−5)').
  3. [Abstract/Conclusion] The word 'parameter-free' is used prominently, but the manuscript does not provide a complete list of numerical inputs (plane-wave cutoff, k-mesh densities, number of principal layers, Gaussian trial orbital widths, etc.). A concise input table or statement of defaults would make the claim verifiable.
  4. [Results, QE paragraph] The text lists several sources of remaining discrepancy but does not quantify them (e.g., how much the calculated QE rises earlier at low excess energy, or the magnitude of the underestimate at higher excess energy). Quantitative error measures would strengthen the comparison and make future improvements trackable.
  5. [End Matter, Optical matrix elements, Eq. (19)] The surface delta-function term assumes a sharp interface at z=0 with a discontinuous dielectric function. For a real Ag(111) surface the profile has finite width; please comment on the sensitivity of the vectorial effect to this idealization and, if possible, test a broadened profile.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the scattering-state construction and observables are computed from first-principles inputs, not from the target data.

full rationale

I walked the derivation chain from the Wannier Hamiltonian and Green-function embedding through the scattering states, quasiparticle self-energy, optical matrix elements, and final QE/MTE observables. The final states are obtained by solving the embedded interface problem (Eq. 6) with lead self-energies from iterative Green functions and quasiparticle self-energies from G0W0/Fan-Migdal (End Matter). The occupied initial states are obtained by diagonalizing a finite Hamiltonian in the same Wannier representation, which is an approximation with convergence asserted but not quantified, but it is not a circular step: it does not use the experimental QE/MTE or the paper's conclusions as inputs. The optical response uses a dielectric function from Ref. [19] and Fresnel equations, with no parameter fitted to the photoemission observables. Self-citations to Refs. [14] and [21] supply vacuum Wannier functions and DFT machinery as building blocks, while Refs. [11] and [20] provide experimental comparisons; none is an unverified premise that contains the paper's conclusion. The finite-slab termination concern is a correctness/missing-support risk, not an equation-level reduction to inputs, so it does not raise the circularity score.

Assumptions & free parameters 1 free parameters · 5 assumptions · 1 invented entities

No fitted parameter enters the visible equations; the two main non-trivial assumptions are the transferability of bulk G0W0/Fan-Migdal self-energies to the interface and the convergence of the occupied-state finite-Hamiltonian truncation. The classical Fresnel/delta-function field profile is another non-trivial modeling assumption. The vacuum-Wannier basis is the method's enabling construct and carries the weight of the whole approach.

free parameters (1)
  • none identified
    The paper claims parameter-free; the complex dielectric function/refractive index of Ag is taken from a first-principles response (Ref. [19]) and the Fresnel fields are computed from it; the self-energies are computed from G0W0 and Fan-Migdal, not fit. No free parameter such as an inelastic mean free path or an empirical optical potential is introduced in the equations shown. The claim is therefor
assumptions (5)
  • domain assumption The hybrid Wannier representation remains short-ranged and converged with a tractable number of layers, including vacuum Wannier functions placed on a close-packed grid (Eqs. 1-3).
    The entire reduction to a one-dimensional open-boundary problem depends on the localization of the Wannier basis and the vacuum-Wannier grid of Ref. [14]; the paper states that without the regular close-packed grid 'random Gaussian guesses would lead to large numerical instabilities.'
  • domain assumption The G0W0 and Fan-Migdal self-energies, computed in bulk or slab, remain valid as local-in-energy and momentum-dependent attenuation when inserted into the Wannier Hamiltonian of the interface (Eqs. 8, 14-17).
    The diagonal quasiparticle approximation Σαβ ≈ δαβ Σαk(E) and the Wannier interpolation of the self-energy presume transferability of the bulk self-energy to the near-surface region.
  • domain assumption The classical Fresnel field profile with a sharp interface, including the delta-function surface term (Eqs. 18-20), adequately represents the microscopic optical field at a metal-vacuum surface.
    The surface term is a classical electrodynamics idealization; the paper treats the field discontinuity as a delta function at z=0, which is an approximation for the real microscopic screening profile.
  • ad hoc to paper Direct diagonalization of a large finite extended Hamiltonian yields converged occupied initial states, including surface states, for the purpose of absolute photoemission amplitudes.
    End Matter 'Occupied-state construction': states are obtained by 'diagonalizing a large finite Hamiltonian' with appended bulk principal layers, with no convergence tests shown. This is a modeling choice that the central results depend on.
  • domain assumption The plane-wave DFT/GGA electronic structure plus Wannier interpolation is an adequate starting point for the near-threshold photoemission matrix elements.
    GGA underpredicts band gaps and can misplace surface-state energies; the paper does not discuss quasiparticle corrections to the initial states, only to the final-state attenuation.
invented entities (1)
  • hybrid vacuum Wannier functions on a regular close-packed lattice
    purpose: Represent the semi-infinite vacuum continuum as a short-ranged one-dimensional chain (Eqs. 1-3 and End Matter).
    This is a methodological construct, not a physical entity; its validity is established by the agreement of the final results, not by an independent falsifiable handle.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Photoemission from Semi-Infinite Crystals: An \emph{Ab Initio} Scattering-State Approach." pith.science (2026). https://pith.science/paper/EUIE6JC7

@misc{pith2026260719551,
  author       = {Pith},
  title        = {Pith review of: Photoemission from Semi-Infinite Crystals: An \emphAb Initio Scattering-State Approach},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/EUIE6JC7}},
  note         = {Machine review of arXiv:2607.19551}
}
read the original abstract

Photoemission is an escape problem, yet first-principles calculations usually trap the electron in a periodic box. We present a parameter-free \emph{ab initio} framework that removes this artificial boundary by constructing open scattering states for semi-infinite crystal-vacuum interfaces from Wannier Hamiltonians and Green-function embedding. The method gives continuum-normalized time-reversed LEED final states with microscopic quasiparticle attenuation. For Ag(111), it predicts absolute quantum efficiency, vectorial photoemission, and mean transverse energy on the experimental scale.

Figures

Figures reproduced from arXiv: 2607.19551 by the authors.

Figure 1
Figure 1. FIG. 1. Partitioning of the semi-infinite crystal-vacuum sys [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2. Inverse LEED scattering states for a semi-infinite [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3. Near-threshold quantum efficiency of Ag(111). Ab [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figures from the paper (1 more)
Figure 4
Figure 4. Figure 4: FIG. 4. Vectorial photoelectric effect in Ag(111). Abso [PITH_FULL_IMAGE:figures/full_fig_p004_4.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

29 extracted references

  1. [19]

    Claudia Ambrosch-Draxl and Jorge O. Sofo. Linear opti- cal properties of solids within the full-potential linearized augmented planewave method.Computer Physics Com- munications, 175(1):1–14, 2006

  2. [14]

    Vacuum wannier functions for first-principles scattering and photoemission

    Tyler Wu and Tom´ as Arias. Vacuum wannier functions for first-principles scattering and photoemission. 2026

  3. [1]

    Sobota, Yu He, and Zhi-Xun Shen

    Jonathan A. Sobota, Yu He, and Zhi-Xun Shen. Angle- resolved photoemission studies of quantum materials. Rev. Mod. Phys., 93:025006, May 2021

  4. [2]

    C. T. Parzyck, A. Galdi, J. K. Nangoi, W. J. I. DeBenedetti, J. Balajka, B. D. Faeth, H. Paik, C. Hu, T. A. Arias, M. A. Hines, D. G. Schlom, K. M. Shen, and J. M. Maxson. Single-crystal alkali antimonide pho- tocathodes: High efficiency in the ultrathin limit.Phys. Rev. Lett., 128:114801, Mar 2022

  5. [3]

    W. E. Spicer. Photoemissive, photoconductive, and op- tical absorption studies of alkali-antimony compounds. Phys. Rev., 112:114–122, Oct 1958

  6. [4]

    Kevin Nangoi, Siddharth Karkare, Ravishankar Sundararaman, Howard A

    J. Kevin Nangoi, Siddharth Karkare, Ravishankar Sundararaman, Howard A. Padmore, and Tom´ as A. Arias. Importance of bulk excitations and coher- ent electron-photon-phonon scattering in photoemission from pbte(111): Ab initio theory with experimental com- parisons.Phys. Rev. B, 104:115132, Sep 2021

  7. [5]

    Dim- itrov, Jai Kwan Bae, Luca Cultrera, Ivan Bazarov, and Andrei Afanasev

    Oksana Chubenko, Siddharth Karkare, Dimitre A. Dim- itrov, Jai Kwan Bae, Luca Cultrera, Ivan Bazarov, and Andrei Afanasev. Monte carlo modeling of spin-polarized photoemission from p-doped bulk gaas.Journal of Ap- plied Physics, 130(6):063101, 08 2021

  8. [6]

    Ultrabright and ultrafast iii–v semiconductor photocathodes.Phys

    Siddharth Karkare, Laurent Boulet, Luca Cultrera, Bruce Dunham, Xianghong Liu, William Schaff, and Ivan Bazarov. Ultrabright and ultrafast iii–v semiconductor photocathodes.Phys. Rev. Lett., 112:097601, Mar 2014

Show all 29 references
  1. [7]

    J.B. Pendry. Theory of photoemission.Surface Science, 57(2):679–705, 1976

  2. [8]

    J. Braun. The theory of angle-resolved ultraviolet pho- toemission and its applications to ordered materials.Re- ports on Progress in Physics, 59(10):1267–1338, 1996

  3. [9]

    E. E. Krasovskii and W. Schattke. Angle-resolved photoemission from surface states.Phys. Rev. Lett., 93:027601, Jul 2004

  4. [10]

    E. E. Krasovskii. Augmented-plane-wave approach to scattering of bloch electrons by an interface.Phys. Rev. B, 70:245322, Dec 2004

  5. [11]

    Chi- ang, and Howard A

    Siddharth Karkare, Weishi Wan, Jun Feng, Tai C. Chi- ang, and Howard A. Padmore. One-step model of pho- toemission from single-crystal surfaces.Phys. Rev. B, 95:075439, Feb 2017

  6. [12]

    Lippmann- schwinger approach for accurate photoelectron wave 6 functions and angle-resolved photoemission spectra from first principles.Physical Review Letters, 135:056403, 2025

    Ji Hoon Ryoo and Cheol-Hwan Park. Lippmann- schwinger approach for accurate photoelectron wave 6 functions and angle-resolved photoemission spectra from first principles.Physical Review Letters, 135:056403, 2025

  7. [13]

    S. A. Rozhkov, V. V. Bakin, H. E. Scheibler, V. S. Ruset- sky, D. V. Gorshkov, D. A. Kustov, V. A. Golyashov, V. L. Alperovich, and O. E. Tereshchenko. Subthermal mean transverse energies induced by electron refraction on the jump in mass at the surface of multialkali photo- c...

  8. [15]

    Mostofi, Jonathan R

    Nicola Marzari, Arash A. Mostofi, Jonathan R. Yates, Ivo Souza, and David Vanderbilt. Maximally localized wannier functions: Theory and applications.Rev. Mod. Phys., 84:1419–1475, Oct 2012

  9. [16]

    On the surface states associated with a periodic potential.Phys

    William Shockley. On the surface states associated with a periodic potential.Phys. Rev., 56:317–323, Aug 1939

  10. [17]

    M. P. Lopez Sancho, J. M. Lopez Sancho, and J. Ru- bio. Highly convergent schemes for the calculation of bulk and surface green functions.Journal of Physics F: Metal Physics, 15(4):851–858, 1985

  11. [18]

    Chico and L

    L. Chico and L. M. Falicov. Electron scattering at inter- faces: A tight-binding approach.Phys. Rev. B, 52:6640– 6646, Sep 1995

  12. [20]

    Karkare, H

    S. Karkare, H. A. Padmore, G. Adhikari, and W. A. Schroeder. Intrinsic emittance of single crystal cathodes. InProceedings of the 9th International Particle Acceler- ator Conference (IPAC2018), pages 4263–4266, Vancou- ver, BC, Canada, 2018. JACoW Publishing

  13. [21]

    Schwarz, Deniz Gunceler, Yalcin Ozhabes, and T.A

    Ravishankar Sundararaman, Kendra Letchworth- Weaver, Kathleen A. Schwarz, Deniz Gunceler, Yalcin Ozhabes, and T.A. Arias. Jdftx: Software for joint density-functional theory.SoftwareX, 6:278–284, 2017

  14. [22]

    Advanced capabilities for materials modelling with quantum espresso.Journal of Physics: Condensed Mat- ter, 29(46):465901, oct 2017

    P Giannozzi, O Andreussi, T Brumme, O Bunau, M Buongiorno Nardelli, M Calandra, R Car, C Cavaz- zoni, D Ceresoli, M Cococcioni, N Colonna, I Carnimeo, A Dal Corso, S de Gironcoli, P Delugas, R A DiStasio, A Ferretti, A Floris, G Fratesi, G Fugallo, R Gebauer, U Gerstmann, F Gi...

  15. [23]

    Quantum espresso: a modular and open-source software project for quantum simulations of materials.Journal of Physics: Condensed Matter, 21(39):395502, sep 2009

    Paolo Giannozzi, Stefano Baroni, Nicola Bonini, Mat- teo Calandra, Roberto Car, Carlo Cavazzoni, Davide Ceresoli, Guido L Chiarotti, Matteo Cococcioni, Ismaila Dabo, Andrea Dal Corso, Stefano de Gironcoli, Ste- fano Fabris, Guido Fratesi, Ralph Gebauer, Uwe Ger- stmann, Christ...

  16. [24]

    Optimization algo- rithm for the generation of oncv pseudopotentials.Com- puter Physics Communications, 196:36–44, 2015

    Martin Schlipf and Fran¸ cois Gygi. Optimization algo- rithm for the generation of oncv pseudopotentials.Com- puter Physics Communications, 196:36–44, 2015

  17. [25]

    Strubbe, Manish Jain, Marvin L

    Jack Deslippe, Georg Samsonidze, David A. Strubbe, Manish Jain, Marvin L. Cohen, and Steven G. Louie. Berkeleygw: A massively parallel computer package for the calculation of the quasiparticle and optical proper- ties of materials and nanostructures.Computer Physics Communicat...

  18. [26]

    M. S. Hybertsen and S. G. Louie. Electron correlation in semiconductors and insulators: Band gaps and quasi- particle energies.Physical Review B, 34(8):5390–5413, 1986

  19. [27]

    da Jornada, Gabriel Antonius, Tonatiuh Rangel, Steven G

    Mauro Del Ben, Felipe H. da Jornada, Gabriel Antonius, Tonatiuh Rangel, Steven G. Louie, Jack Deslippe, and Andrew Canning. Static subspace approximation for the evaluation ofg 0w0 quasiparticle energies within a sum- over-bands approach.Physical Review B, 99(12):125128, 2019

  20. [28]

    da Jornada, Andrew Canning, Nicole Wichmann, Koushik Raman, R

    Mauro Del Ben, Felipe H. da Jornada, Andrew Canning, Nicole Wichmann, Koushik Raman, R. Sasanka, Chao Yang, Steven G. Louie, and Jack Deslippe. Large-scale gw calculations on pre-exascale hpc systems.Computer Physics Communications, 235:187–195, 2019

  21. [29]

    Electron-phonon interactions from first principles.Rev

    Feliciano Giustino. Electron-phonon interactions from first principles.Rev. Mod. Phys., 89:015003, Feb 2017. 7 End Matter Hybrid Wannier construction.—All density- functional-theory calculations were performed using JDFTx [21] and Quantum ESPRESSO [22, 23], with SG15 optimized...

Pith tools

Reviewed August 1, 2026 · model on record in the stance chip above.