REVIEW 4 major objections 5 minor 23 references
Effects of gold cluster intercalation in graphene: stationary waves and modified QPI features
T0 review · 4 major / 5 minor · reviewed 2026-07-31 · grok-4.5
Pith's one-line read Gold atoms under graphene hollow sites act as six-carbon ring scatterers that create M-centered QPI ellipses and nearly stationary standing waves.
desk verdict Solid, focused T-matrix explanation of the Au-cluster M-centered QPI and stationary waves; ring geometry is the real contribution, with U and multi-cluster coherence left soft. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The six-site ring T-matrix: the hollow-site impurity is projected onto the six neighbouring carbon pz orbitals, the full multiple-scattering T-matrix is solved in that six-dimensional space, and the result is embedded back into graphene Bloch states, generating the hollow-site form factor that selects M-centered double ellipses.
What would settle it
An FT-STS map of the same Au-cluster phase that shows only extended line-like QPI features instead of compact double ellipses around M, or an independent structural measurement that places the intercalated gold under carbon atoms rather than hollow sites, would directly falsify the central mechanism.
Extended reading notes
Core claim
Gold intercalation below graphene hollow sites induces a local electrostatic potential on the six surrounding carbon atoms; within a single-impurity T-matrix treatment this ring-like scatterer produces elliptical QPI structures centered near the M points that naturally generate the nearly stationary real-space standing waves observed in the Au-cluster phase, while the LDOS contrast on and off a small cluster is strongly energy-dependent.
Load-bearing premise
The claim rests on gold atoms sitting under hollow sites and acting, at the energies of interest, as a roughly constant electrostatic shift on a six-carbon ring, with inter-cluster multiple scattering neglected when computing the local density of states.
Editorial extensions
If this is right
- Ring-like hollow-site scatterers, not ordinary onsite defects, are required to produce the observed M-centered elliptical QPI.
- Standing-wave periods stay nearly constant over a bias window because the relevant wave-vectors remain pinned near M by the upper van-Hove warping.
- LDOS contrast between cluster and background must reverse sign multiple times between the lower and upper van-Hove singularities.
- A minimal double-lobe model of the M-centered QPI already reproduces the coherence length of the experimental standing waves.
Reading between the lines
- The same hollow-site ring geometry should appear in other noble-metal intercalants that occupy hollow sites, offering a route to engineer stationary interference patterns by choice of adsorbate.
- If coherent inter-cluster scattering is what enhances the experimental M-weight relative to small-q background, deliberately ordered cluster lattices could further amplify or suppress selected QPI channels.
- Energy-dependent LDOS contrast on individual clusters could be used as a local spectroscopic fingerprint of hollow-site versus top-site intercalation.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript argues that the unusual M-centered elliptical QPI features and nearly stationary real-space standing waves in the Au-cluster (“ostrich-leather”) phase of graphene/SiC arise from Au atoms sitting under graphene hollow sites and acting as six-carbon ring scatterers. Using a sixth-nearest-neighbour tight-binding graphene model and a single-impurity T-matrix for an effective constant electrostatic shift U on the ring (Eqs. 13–19), the authors obtain compact double-lobe structures near the M points (Fig. 2), in contrast to the extended line-like QPI of single-site impurities. An analytic double-lobe Fourier model then produces coherent, weakly energy-dependent standing waves. A real-space LDOS calculation for a three-ring cluster further yields strongly energy-dependent on/off-cluster contrast with multiple sign changes, qualitatively matching STM observations.
Significance. If correct, the work supplies a concrete microscopic link between intercalant registry (hollow-site ring geometry) and a distinctive class of FT-STS signatures that differ from standard intra-/intervalley graphene QPI. The ring-versus-single-site comparison in Sec. III.B and Fig. 2 is a clear, falsifiable diagnostic, and the double-lobe construction transparently connects momentum-space lobe geometry to real-space coherence length and near-stationarity. That is a useful addition to the intercalated-graphene and QPI literature. Strengths include an explicit T-matrix formulation for a finite-range ring impurity, a controlled 6NN band structure with stated hoppings, and direct side-by-side theory–experiment figures for QPI, filtered standing waves, and LDOS contrast.
major comments (4)
- [Sec. III.A, Eq. (14); Ref. 16] The load-bearing structural premise—that Au occupies hollow sites and, at the relevant energies, reduces to V_eff ≃ U I_6 on a six-carbon ring—is imported from the companion microscopic study Ref. 16, which is still listed as “in preparation” (Abstract, Sec. I, Sec. III.A). For a standalone journal article this is insufficient: either the essential structural/DFT evidence and the justification for dropping the hybridization term in Eq. (14) must be included (or summarized with data) here, or publication should await Ref. 16 so that the premise is citable and checkable.
- [Sec. III.A–B, Eqs. (14), (17); Fig. 2] U = 3 eV is introduced only as “of the same order of magnitude as that of Ref. 16” (Sec. III.A), with no sensitivity analysis. Because T(ω) = V_eff [I_6 − g_0 V_eff]^{-1} is nonlinear in U, the compactness, lobe separation 2Δq, and weak energy drift of the M-centered features (Fig. 2, top row)—and therefore the stationarity argument built from the analytic double-lobe model—can change if U differs appreciably or if the off-resonant Au term in Eq. (14) is not negligible. A short scan over U (and a check with a weakly energy-dependent V_eff) is needed to show that the double-ellipse topology is robust rather than tuned.
- [Sec. IV (paragraph on small-q vs M-centered weight); Fig. 3 vs Fig. 2] Sec. IV states that the same single-impurity / incoherent three-ring LDOS is dominated by small-q scattering, so M-centered weight appears only after Fourier filtering, whereas experiment shows strong unfiltered M-centered intensity. The text invokes possible coherent inter-cluster multiple scattering to reconcile this, but that physics lies outside the T-matrix used for the central QPI claim and is not computed. Either a multi-impurity coherent calculation (or a controlled argument why ring form factors alone suffice once disorder-averaged) should be provided, or the claim that the single-ring T-matrix “explains” the experimental FT-STS weight distribution should be narrowed accordingly.
- [Sec. II; Sec. III.B; Figs. 1–3, 5–6] Theory–experiment energy alignment is acknowledged as only qualitative (Dirac point and VHS positions uncertain; experimental maps at ~0.8–1 eV vs theory near the upper VHS at ~1.74 eV). Given that the upper-VHS warping of the 6NN model is essential to the extended spectral weight near M (Fig. 1), the manuscript should either (i) show QPI maps over a broader energy window including a shifted μ consistent with experiment, or (ii) state more sharply which observables (lobe topology, stationarity, LDOS sign changes) are predicted to be robust under μ/VHS shifts and which are not.
minor comments (5)
- [Sec. II; Figs. 2, 5] The chemical potential is set so the Dirac point is at E_F (μ = −0.377 eV), while experiment may have the Dirac point down to −0.5 eV; a single sentence in the figure captions of Figs. 2 and 5 stating the assumed μ would help readers compare panels.
- [Sec. II, after Eq. (11)] Broadening η = 0.2 eV is large; the claim that M-centered ellipses persist for η ∼ 0.1–0.3 is useful but only stated in text—adding one supplemental or inset panel at a second η would make that check visible.
- [Sec. III.C, Eq. (21)] In Eq. (21) the double-lobe model parameters (σ1, σ2, Δq) are fitted per energy; briefly tabulating them (as already partly done in the text) would aid reproducibility of the bottom row of Fig. 2.
- [Fig. 6] Fig. 6 axis labels mix “Tension (V)” and energy; standardize units and clarify whether the experimental curve is raw bias or an estimated energy scale.
- [Sec. II; Sec. III.C; Sec. IV] Typos / wording: “thenth” → “the nth” (Sec. II); “approximatively” → “approximately” (several places); “resonatorlike” hyphenation; ensure consistent notation for M vs M points.
Circularity Check
Hollow-site ring premise is load-bearing self-citation to unfinished Ref. 16; QPI/standing-wave math is independent forward calculation, not tautological.
-
self citation load bearing
[Introduction; Sec. III.A (Hollow-site Au impurity model); Ref. 16]
"In a recent work16 we developed an effective microscopic description of the cluster phase which allowed us to conclude that the Au atoms sit underneath the graphene hollow sites and couple to the six neighbouring carbon atoms either through a local electrostatic potential or through Au–C hybridization. ... Since the QPI features considered here lie far from the relevant Au hybridization levels (∼ -2eV – -3eV)16, we use the minimal approximation Veff(ω)≃UI6; the fitted value of U may then also absorb the slowly varying off-resonant contribution of the Au orbitals. In what follows we will take U"
The load-bearing structural premise—Au under hollow sites acting as a six-carbon ring scatterer, and hybridization reduced to a constant electrostatic U—is justified only by citation to the authors’ own unfinished companion (Ref. 16, ‘in preparation’), not by independent derivation or external measurement in this paper. The entire QPI/standing-wave claim is conditioned on that imported geometry and mechanism.
-
fitted input called prediction
[Sec. II (η); Sec. III.A–B (U=3 eV); Fig. 2 comparison to experiment]
"with a broadening η = 0.2 eV. We use this relatively large broadening to mimic, at the level of the effective model, the disorder and finite lifetime associated with the Au-cluster phase. ... In what follows we will take U =3eV, a value which is of the same order of magnitude as that of Ref. 16. ... These two-lobe features reproduce the characteristic structure of the experimental FT-STS maps shown in Fig. 3."
U and η are not predicted from first principles in this work; they are chosen by order-of-magnitude match to Ref. 16 and to mimic cluster-phase disorder so that the T-matrix maps resemble experiment. Because T(ω)=V_eff[I−g0 V_eff]−1 is nonlinear in U, lobe compactness and separation (and thus the stationarity narrative) are partly tuned by these inputs rather than fixed a priori—mild fitted-input circularity, not full tautology.
full rationale
The paper’s central computational chain—6NN graphene Green’s function, six-site ring T-matrix, δρ(q,ω), and Fourier connection to real-space waves—is a genuine forward calculation once the impurity geometry and U are fixed. It is not self-definitional: elliptical M-centered lobes are not built into the definition of V_eff, and the single-site control (Fig. 2, last column) shows the ring form factor is doing real work. Circularity is limited to (i) importing the hollow-site placement and the reduction V_eff≃U I_6 from the authors’ own companion manuscript Ref. 16 (still “in preparation”), and (ii) choosing U=3 eV and η∼0.2 eV by order-of-magnitude / phenomenological match rather than predicting them here. The analytic double-lobe model is fitted to the authors’ own T-matrix QPI maps to illustrate standing-wave coherence; that is pedagogical, not a fit-to-experiment renamed as prediction. Overall this is moderate self-citation dependence on an unfinished companion for the structural premise, with independent content in the QPI and LDOS calculations—score 3, not a by-construction derivation.
Assumptions & free parameters
free parameters (5)
- U (ring electrostatic shift) =
3 eV
- η (spectral broadening) =
0.2 eV
- μ (chemical potential) =
−0.377 eV
- Double-lobe Gaussian parameters (σ1, σ2, Δq) =
e.g. 2Δq ~ 0.64→0.42 Å⁻¹, σ1 ~ 0.33→0.41 Å⁻¹, σ2 ~ 0.15→0.11 Å⁻¹
- 6NN hopping set (t1…t6) =
t1=−2.937, t2=0.249, t3=−0.260, t4=0.025, t5=0.050, t6=−0.024 eV
assumptions (6)
- domain assumption Single-impurity T-matrix on a projected six-site ring embedded in Bloch graphene states correctly captures the FT-STS signal of the cluster phase.
- domain assumption Away from Au hybridization resonances (~−2 to −3 eV), Veff(ω) ≈ U I6 on the six carbons is sufficient.
- ad hoc to paper Au atoms occupy graphene hollow sites and couple primarily to the six neighboring pz orbitals.
- domain assumption Sixth-nearest-neighbor tight-binding with given hoppings adequately represents the warped upper van-Hove contours.
- ad hoc to paper Linear superposition of three single-ring LDOS responses approximates a compact three-Au cluster (inter-ring multiple scattering negligible).
- standard math Standard retarded Green’s function and LDOS/QPI trace formulas in the sublattice basis.
invented entities (1)
-
Hollow-site six-carbon ring scatterer (effective Au impurity model)
Cite this review
Pith. "Pith review of Effects of gold cluster intercalation in graphene: stationary waves and modified QPI features." pith.science (2026). https://pith.science/paper/QJMMZAA7
@misc{pith2026260728297,
author = {Pith},
title = {Pith review of: Effects of gold cluster intercalation in graphene: stationary waves and modified QPI features},
year = {2026},
howpublished = {\url{https://pith.science/paper/QJMMZAA7}},
note = {Machine review of arXiv:2607.28297}
}
read the original abstract
Gold intercalation beneath epitaxial graphene on SiC produces a cluster phase with unusual standing waves and quasiparticle-interference (QPI) features concentrated near the graphene M points. We show that this can be explained by Au intercalation below graphene hollow sites, which induces a local scattering potential on the six surrounding carbon atoms. Within a T-matrix treatment, this ring-like scatterer produces elliptical QPI structures centered near M, in agreement with the experimental FT-STS measurements. We further show that these QPI features naturally generate the nearly stationary standing-wave patterns observed in real space. Finally, we compute the local-density-of-states contrast on and off a small cluster and show that its sign and magnitude are strongly energy dependent, consistent with the experimental observations.
Figures
Figures from the paper (2 more)
Reference graph
Works this paper leans on
-
[1]
Cranney, F
M. Cranney, F. Vonau, P. B. Pillai, E. Denys, D. Aubel, and L. Simon, ``Superlattice of resonators on monolayer graphene created by intercalated gold nanoclusters,'' Europhys. Lett. 91, 66004 (2010)
2010
-
[2]
Riedl, C
C. Riedl, C. Coletti, T. Iwasaki, A. A. Zakharov, and U. Starke, ``Quasi-free-standing epitaxial graphene on SiC obtained by hydrogen intercalation,'' Phys. Rev. Lett. 103, 246804 (2009)
2009
-
[3]
Gierz, T
I. Gierz, T. Suzuki, R. T. Weitz, D. S. Lee, B. Krauss, C. Riedl, U. Starke, H. H\"ochst, J. H. Smet, C. R. Ast, and K. Kern, ``Electronic decoupling of an epitaxial graphene monolayer by gold intercalation,'' Phys. Rev. B 81, 235408 (2010)
2010
-
[4]
Forti, K
S. Forti, K. V. Emtsev, C. Coletti, A. A. Zakharov, C. Riedl, and U. Starke, ``Large-area homogeneous quasifree standing epitaxial graphene on SiC(0001): Electronic and structural characterization,'' Phys. Rev. B 84, 125449 (2011)
2011
-
[5]
Briggs, Z
N. Briggs, Z. M. Gebeyehu, A. Vera, T. Zhao, K. Wang, A. De La Fuente Duran, B. Bersch, T. Bowen, K. L. Knappenberger, and J. A. Robinson, ``Epitaxial graphene/silicon carbide intercalation: a minireview on graphene modulation and unique 2D materials,'' Nanoscale 11, 15440--15447 (2019)
2019
-
[6]
Premlal, M
B. Premlal, M. Cranney, F. Vonau, D. Aubel, L. Cattin, F. Scheidt, C. Deville-Cavellin, and L. Simon, ``Surface intercalation of gold underneath a graphene monolayer on SiC(0001) studied by scanning tunneling microscopy and spectroscopy,'' Appl. Phys. Lett. 94, 263115 (2009)
2009
-
[7]
M. N. Nair, M. Cranney, F. Vonau, D. Aubel, and L. Simon, ``High van Hove singularity extension and Fermi velocity increase in epitaxial graphene functionalized by intercalated gold clusters,'' Phys. Rev. B 85, 245421 (2012)
2012
-
[8]
M. N. Nair, M. Cranney, T. Jiang, S. Hajjar-Garreau, D. Aubel, F. Vonau, A. Florentin, E. Denys, M.-L. Bocquet, and L. Simon, ``Noble-metal intercalation process leading to a protected adatom in a graphene hollow site,'' Phys. Rev. B 94, 075427 (2016)
2016
Show all 23 references
-
[9]
Rosenzweig, H
P. Rosenzweig, H. Karakachian, S. Link, K. K\"uster, and U. Starke, ``Tuning the doping level of graphene in the vicinity of the van Hove singularity via ytterbium intercalation,'' Phys. Rev. B 100, 035445 (2019)
2019
-
[10]
Rosenzweig, H
P. Rosenzweig, H. Karakachian, D. Marchenko, K. K\"uster, and U. Starke, ``Overdoping graphene beyond the van Hove singularity,'' Phys. Rev. Lett. 125, 176403 (2020)
2020
-
[11]
ohr, K. K\
S. Link, S. Forti, A. St\"ohr, K. K\"uster, M. R\"osner, D. Hirschmeier, C. Chen, J. Avila, M. C. Asensio, A. A. Zakharov, T. O. Wehling, A. I. Lichtenstein, M. I. Katsnelson, and U. Starke, ``Introducing strong correlation effects into graphene by gadolinium intercalation,'' ...
2019
-
[12]
Zaarour, V
A. Zaarour, V. Malesys, J. Teyssandier, M. Cranney, E. Denys, J. L. Bubendorff, A. Florentin, L. Josien, F. Vonau, D. Aubel, A. Ouerghi, C. Bena, and L. Simon, ``Flat band and Lifshitz transition in long-range-ordered supergraphene obtained by erbium intercalation,'' Phys. Rev...
2023
-
[13]
G. M. Rutter, J. N. Crain, N. P. Guisinger, T. Li, P. N. First, and J. A. Stroscio, ``Scattering and interference in epitaxial graphene,'' Science 317, 219--222 (2007)
2007
-
[14]
Brihuega, P
I. Brihuega, P. Mallet, C. Bena, S. Bose, C. Michaelis, L. Vitali, F. Varchon, L. Magaud, K. Kern, and J.-Y. Veuillen, ``Quasiparticle chirality in epitaxial graphene probed at the nanometer scale,'' Phys. Rev. Lett. 101, 206802 (2008)
2008
-
[15]
Bena, ``Effect of a single localized impurity on the local density of states in monolayer and bilayer graphene,'' Phys
C. Bena, ``Effect of a single localized impurity on the local density of states in monolayer and bilayer graphene,'' Phys. Rev. Lett. 100, 076601 (2008)
2008
-
[16]
Kumari, A
P. Kumari, A. Zobelli, I. de Melo Froldi, A. Crepieux, L. Simon, and C. Bena, ``An iterative method bridging DFT, disorder averaging, and experiment in intercalated materials: application to Au-intercalated graphene,'' in preparation
-
[17]
Reich, J
S. Reich, J. Maultzsch, C. Thomsen, and P. Ordej\'on, ``Tight-binding description of graphene,'' Phys. Rev. B 66, 035412 (2002)
2002
-
[18]
G. H. Wannier, ``The structure of electronic excitation levels in insulating crystals,'' Phys. Rev. 52, 191--197 (1937)
1937
-
[19]
Kaasbjerg, ``Atomistic T -matrix theory of disordered two-dimensional materials: Bound states, spectral properties, quasiparticle scattering, and transport,'' Phys
K. Kaasbjerg, ``Atomistic T -matrix theory of disordered two-dimensional materials: Bound states, spectral properties, quasiparticle scattering, and transport,'' Phys. Rev. B 101, 045433 (2020)
2020
-
[20]
P. Kot, V. Fodje M\'ela, J. F. Annett, and S. Rossi, ``Band dispersion of graphene with structural defects,'' Phys. Rev. B 101, 235116 (2020)
2020
-
[21]
A. A. Mostofi, J. R. Yates, Y.-S. Lee, I. Souza, D. Vanderbilt, and N. Marzari, ``wannier90: A tool for obtaining maximally-localised Wannier functions,'' Comput. Phys. Commun. 178, 685--699 (2008)
2008
-
[22]
Giannozzi, S
P. Giannozzi, S. Baroni, N. Bonini, et al., ``QUANTUM ESPRESSO: a modular and open-source software project for quantum simulations of materials,'' J. Phys.: Condens. Matter 21, 395502 (2009)
2009
-
[23]
Giannozzi, O
P. Giannozzi, O. Andreussi, T. Brumme, et al., ``Advanced capabilities for materials modelling with Quantum ESPRESSO,'' J. Phys.: Condens. Matter 29, 465901 (2017)
2017
Reviewed July 31, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.