{"id":"47c029e1-bfce-450b-9925-248455895a9d","arxiv_id":"2411.09669","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"A layer-KKR k-space implementation of the one-step photoemission model simulates hard X-ray photoelectron diffraction and circular dichroism for Si(100) and Ge(100) at 6 keV.","lead":"The authors present a k-space, layer-KKR implementation of photoelectron diffraction in the SPRKKR package and test it at 6 keV on silicon and germanium surfaces. Simulated total intensities, CDAD differences and asymmetries reproduce the main experimental Kikuchi features, though only qualitatively.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The quantitative accuracy claim is unsupported by the paper's own convergence tests: lmax=4 is admitted (Section 4.2) to be insufficient for quantitative intensities, and Fig. 4 shows it is not converged in the MsSpec cross-check.","rationale":"The reader's weakest_assumption correctly identifies the single-site final-state approximation as load-bearing. My stress-test sharpens the same point: the manuscript itself provides the decisive evidence against its central accuracy claim. Section 4.2 explicitly concedes that lmax > 30 is required for quantitative intensities, and Figure 4b shows that the lmax = 4 value used for all headline patterns is not fully converged even for the cross-section test. Since the claimed advantage over real-space cluster codes is precisely the absence of angular-momentum convergence problems, this self-admitted limitation directly undercuts the abstract's '20-8000 eV without angular momentum or cluster size convergence issues' statement. The experiment-theory agreement is also only qualitative: the authors apply Gaussian and Perona-Malik smoothing before comparison, and no error bar or correlation metric is provided. The single-site approximation could be physically reasonable at 6 keV because forward scattering dominates, but that is a physics argument that must be demonstrated, not assumed; the paper provides no full multiple-scattering benchmark for the CDAD patterns themselves. I do not think the method should be rejected: it is a plausible and useful k-space approach, and the visual matches, particularly for Si 2p3/2 and Ge 3p3/2, are encouraging. But the central claim of quantitative accuracy across a broad energy range is not supported by the provided evidence. The appropriate disposition remains CONDITIONAL, as the reader concluded: the manuscript should be published only after a quantitative validation benchmark, public artifacts, and a scoped energy-range statement that acknowledges the lmax limitation. I therefore keep the reader's verdict unchanged, while noting that my concern is slightly more focused on the internal contradiction between Section 4.2 and the abstract than on a purely external failure of the single-site approximation.","tokens_in":24460,"tokens_out":3474,"duration_ms":35789,"concrete_test":"Run the same Si 2p3/2 case at hnu = 6 keV with the MsSpec full multiple-scattering cluster code at lmax = 24 (or 30), using identical phase shifts and geometry, and compare against the LKKR lmax = 4 result for A_CDAD and I_TOT with a quantitative metric such as normalized cross-correlation or mean absolute difference over the measured (kx,ky) range. If the discrepancy exceeds the experimental pattern variability (e.g., RCP/LCP noise after binning), the single-site, lmax = 4 approximation is not quantitatively valid at 6 keV and the abstract's 'accurate...without angular momentum convergence issues' claim fails. For the 20-8000 eV range, repeat at a low-energy point (e.g., 100 eV) where backscattering and multiple scattering are known to matter; if the single-site layer approximation fails there, the energy-range claim must be rescoped.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's strongest claim—efficient and accurate hard-XPD/CDAD simulations without angular-momentum convergence issues—is contradicted by its own convergence data. Section 4.2 states 'we are aware that we would need to approach lmax > 30 ... to make a quantitative comparison between experimental and theoretical intensities,' yet all headline patterns use lmax = 4. The MsSpec cross-check (Fig. 4b) explicitly shows lmax = 4 is 'a bit far from convergence' relative to lmax = 24, with shape differences near phi = 30 and 60 degrees. Therefore the validation rests on (i) an unconverged angular-momentum truncation, (ii) a single-site per-layer final-state approximation whose omitted intra-layer multiple scattering is never benchmarked against a full multiple-scattering calculation, and (iii) qualitative visual matching after Gaussian and Perona-Malik smoothing, with no quantitative agreement metric. These are not independent faults but one load-bearing gap: there is no demonstrated connection between the computed intensities and the true multiple-scattering final state at 6 keV. The 20-8000 eV claim is likewise extrapolated from a 4690-6690 eV window. If the method is intended as a practical qualitative pattern-prediction tool, the abstract and Section 3 should be scoped accordingly.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a k-space implementation of the one-step model of photoemission based on the layer-KKR method in the SPRKKR package, aimed at hard X-ray photoelectron diffraction (XPD) and circular dichroism in angular distributions (CDAD). The method is validated against momentum-microscopy measurements for Si 2p3/2, Ge 2p3/2, and Ge 3p3/2 at 6 keV, with the authors claiming qualitative reproduction of Kikuchi bands, mirror-plane symmetries, and several fine features, while attributing remaining discrepancies to potential choices, inelastic scattering, lattice vibrations, and crystal imperfections. The central claim is that the method covers a wide energy range (20–8000 eV) without the angular-momentum and cluster-size convergence issues of real-space methods.","tokens_in":24581,"tokens_out":5654,"duration_ms":52833,"significance":"If validated, the layer-KKR/one-step approach would offer a practical reciprocal-space alternative to real-space cluster codes, for which lmax grows to roughly 100 at 10 keV, and would extend the SPRKKR package to hard-X-ray CDAD simulations. A clear strength is that the method is compared against independent experimental data without fitting the compared datasets: lmax and the Ghkl expansion are justified by convergence tests, and the imaginary potentials V0i are standard phenomenological inputs. The CDAD antisymmetry and mirror-plane features are reproduced in the simulated patterns. However, the quantitative accuracy claim is currently unsupported by the convergence data, and the validation is largely qualitative. The work at this stage is best viewed as a demonstration of qualitative pattern prediction rather than as a validated quantitative method.","major_comments":[{"comment":"The paper explicitly states near the end of Section 4.2 that 'we would need to approach lmax > 30' for a quantitative comparison of intensities, and Fig. 4b shows that lmax = 4 is 'a bit far from convergence' with shape differences near φ = 30° and 60°. Nevertheless, all headline comparisons in Figs. 6–9 use lmax = 4, while the abstract and Section 3 claim 'efficient and accurate calculations' without angular-momentum convergence issues. This is an internal contradiction at the core of the paper: the validation does not establish that the computed intensities represent the fully multiple-scattered final state. The authors should either provide a benchmark showing that lmax = 4 converges the features used for the experimental comparison, or soften the convergence and accuracy claims to 'qualitative pattern prediction.'","section":"4.2"},{"comment":"The statement 'This study focuses exclusively on single-site scattering when calculating the final state' means that intra-layer multiple scattering is neglected, justified by forward-scattering dominance at multi-keV energies. This assumption is never benchmarked against a full multiple-scattering calculation. The MsSpec cross-check in Fig. 4 is only an lmax convergence test and does not test the single-site-per-layer approximation. Because the correctness of the Kikuchi bands and CDAD features depends directly on this assumption, a quantitative benchmark is needed—for example, a comparison of the layer-KKR intensities with and without intra-layer multiple scattering, or a cluster calculation using the same potential and geometry. Without such a test, the agreement with experiment could be fortuitous or the result of the truncation.","section":"4.2"},{"comment":"The evaluation of experiment–theory agreement in Figs. 6–9 is entirely qualitative: visual inspection after Gaussian (σ = 2) and Perona–Malik (λ = 1) smoothing, with no quantitative agreement metric. The number of smoothing parameters is not justified systematically, and Section 3 states 'The observed outcomes were reproduced in the hard X-ray regime,' which is stronger than what the evidence supports. The authors should add a quantitative comparison measure (e.g., a reliability factor, Pearson correlation, or feature-based metric) for at least the total-intensity and ACDAD patterns, or explicitly limit the claim to reproduction of qualitative symmetry and band positions.","section":"2.3"},{"comment":"The abstract claims the method addresses kinetic energies of 20–8000 eV without convergence problems, but the experimental validation covers only final-state energies in the window 4.69–6.69 keV (Fig. 5 and Section 2.3). No test outside this narrow window is shown, so the wide-range claim is an extrapolation. To support the 20–8000 eV statement, the authors should either add a validation or convergence test at another energy (for example, near 20 eV or at 8 keV) or state explicitly that this range refers to formal applicability rather than to validated performance.","section":"Abstract / 2.2"},{"comment":"Section 2.1 states that with lmax = 4, 'it is possible to achieve satisfactory agreement between simulations and experiments,' only after noting in the same section and in Fig. 4b that lmax = 4 is 'a bit far from convergence' relative to lmax = 16 and 24. These statements are not reconciled. The choice of lmax = 4 is justified in Section 4.2 by 'computational time and memory limitations,' not by convergence. The manuscript should clearly distinguish 'converged' from 'practically usable for qualitative features' and should state the known truncation error separately for each headline figure.","section":"2.1"}],"minor_comments":[{"comment":"The definition 'αk = σκ ⊗ σκ, k = (x, y, z)' appears to be a typo; the Dirac alpha matrices should be defined with specific Pauli matrices, e.g., αx = σ1 ⊗ σ1, and β should be defined accordingly. Please correct the notation.","section":"4.2, Eq. (1)"},{"comment":"The phrase 'result in I non direct transitions (Eq. (4))' is unclear; 'non-direct' or 'incoherent' direct transitions should be explicitly defined.","section":"4.2"},{"comment":"The statement 'The agreement between observed and computed intensity (a-b), CDAD difference (c,d), and ACDAD (e,f) looks quite reasonable' is followed later by 'the consensus between the experimented and simulated ACDAD is far from perfect.' Please ensure that the overall assessment is consistent across the different quantities.","section":"2.3"},{"comment":"The phrase 'convoluted-calculated' is awkward; 'convolved' is standard and should be used in the captions and text.","section":"Figure captions 6, 8, 9"},{"comment":"The sentence 'the position of the photoemission peaks as well as their positions do not change' is redundant; please revise.","section":"2.1"},{"comment":"The number of Ghkl vectors is not consistent across figures (161/113 in Section 2.2, 193/177/241 in Section 2.3); please state how the final values were selected.","section":"2.2 / 2.3"}],"recommendation":"major_revision","confidential_remarks":"The manuscript would benefit from an explicit statement that the method is currently validated for qualitative pattern prediction, not quantitative intensity analysis. The claim 'without angular momentum convergence issues' should be replaced or qualified by 'reduced angular momentum requirements compared to real-space cluster methods,' which is what the data actually support. The lmax and single-site assumptions are the key weaknesses; addressing them in a revision with at least one targeted benchmark would substantially strengthen the paper."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Worth a look if you work on hard X-ray photoelectron diffraction or momentum microscopy. The genuinely new content is a k-space layer-KKR implementation of the one-step model inside SPRKKR, with CDAD calculations for Si 2p and Ge 3p at 6 keV that reproduce main Kikuchi bands, mirror planes, and a fair amount of fine structure. The literature review is thorough, the convergence tests on G_hkl and lmax are documented, and the authors are honest about many limitations: LDA/ASA potentials, phenomenological V0i, lattice vibrations, and the need for lmax > 30 for quantitative intensities.\n\nThe soft spots are real but mostly about scope, not a fatal flaw. The abstract's claim of '20–8000 eV without angular momentum or cluster size convergence problems' is not supported by the 4690–6690 eV tests shown. Figure 4's MsSpec cross-check actually shows lmax = 4 is 'a bit far from convergence' against lmax = 24, so the headline patterns are computed with an unconverged truncation. That is acceptable for qualitative pattern prediction, but the paper does not consistently frame it that way. Second, the final state uses single-site scattering within each layer, with inter-layer multiple scattering in a plane-wave basis. That is a reasonable high-energy approximation, but it is never benchmarked against a full multiple-scattering calculation, so the reader cannot tell how much the missing intra-layer scattering bends the intensities. Third, the experiment–theory comparison is visual, aided by Gaussian and Perona–Malik smoothing, with no numerical agreement metric. The matches shown are suggestive, not demonstrated quantitatively. Fourth, code and data are only available on request.\n\nNone of this kills the paper. The k-space route is a practical alternative to real-space clusters, and the authors are candid about the gaps. What is missing is a scoped statement of what the current implementation can claim (qualitative pattern reproduction at multi-keV, with lmax and single-site caveats), and ideally one quantitative validation case. I would send this to referees; the method is useful and the community would benefit from a corrected version. I would cite it as a methodological reference, though not for quantitative accuracy.","headline":"A useful k-space layer-KKR method for hard-XPD that reproduces Kikuchi patterns qualitatively, but the paper's abstract overclaims convergence and quantitative accuracy relative to its own lmax=4 tests.","tokens_in":25346,"tokens_out":2100,"would_cite":true,"duration_ms":21351,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["79.60.-i"],"model":"deepseek-v4-flash","headline":"A k-space layer-KKR implementation of the one-step photoemission model reproduces measured hard X-ray photoelectron diffraction and circular dichroism for Si(100) and Ge(100) core levels at 6 keV, while avoiding the large angular-momentum…","keywords":["hard X-ray photoelectron diffraction","layer Korringa-Kohn-Rostoker method","one-step photoemission model","circular dichroism in angular distributions","Kikuchi diffraction","momentum microscopy","core-level photoemission","multiple scattering theory"],"falsifier":"Recompute the Si 2p$_{3/2}$ 6 keV CDAD patterns with intra-layer multiple scattering included, or with an angular-momentum cutoff near 30, and compare the Kikuchi-band intensities and the $A_{\\rm CDAD}$ maps with the $l_{\\max}=4$ results; if the intensity ratios shift by more than the experimental noise, the single-site-per-layer approximation is falsified.","tokens_in":24105,"feed_emoji":"🔬","tokens_out":10054,"duration_ms":87606,"temperature":0.7,"pith_summary":"This paper argues that hard X-ray photoelectron diffraction (XPD) and the circular dichroism in its angular distribution (CDAD) can be computed efficiently in reciprocal space by combining the layer-KKR multiple-scattering scheme with the one-step photoemission model. The central claim is that at photoelectron kinetic energies near 6 keV the crystal can be treated as a stack of atomic layers with single-site scattering inside each layer and plane-wave propagation between layers, so Kikuchi-band diffraction is captured with a small angular-momentum cutoff rather than the very large phase-shift sets required by real-space cluster codes. The authors report that simulated total intensities, CDAD differences, and CDAD asymmetries for Si 2p$_{3/2}$ and Ge 3p$_{3/2}$ emitted from (100) surfaces under 6 keV circularly polarized light reproduce the measured momentum-microscope patterns, with Ge 2p$_{3/2}$ in qualitative agreement. If the claim holds, the same implementation can describe photoelectron diffraction from ultraviolet to hard X-ray energies and help separate diffraction structure from valence-band ARPES signals.","feed_headline":"Layer model reproduces hard X-ray photoelectron diffraction at 6 keV","feed_subtitle":"K-space layer-KKR model reproduces Si and Ge core-level CDAD at 6 keV, easing hard-X-ray simulations.","key_machinery":"The load-bearing object is the layer-KKR expansion of the photoelectron final state. The final state is a time-reversed LEED state: an incoming plane wave from the detector direction is scattered by the semi-infinite crystal, with layer-by-layer scattering matched by expanding partial waves inside a layer into plane waves labeled by the reciprocal lattice vectors $\\vec{G}_{hkl}$ between layers. This mixed partial-wave/plane-wave basis decouples the two convergence parameters: $l_{\\max}$ only needs to describe single-site scattering within one layer, while the $\\vec{G}_{hkl}$ expansion carries the long-range interlayer diffraction that produces the Kikuchi structure. The transition is evaluated with the relativistic one-step photoemission formula from a core-level initial state, with small phenomenological imaginary self-energies for initial and final states. This machinery lets a calculation at $l_{\\max}=4$ reproduce patterns for which cluster approaches at 10 keV would need $l_{\\max}$ near 100.","core_discovery":"The paper's central discovery is that the one-step model of photoemission, reformulated on the layer-KKR Green's-function scheme, describes hard X-ray photoelectron diffraction without the convergence bottlenecks of real-space cluster methods. The crystal is approximated by an infinite stack of atomic layers; within each layer the photoelectron final state is expanded in partial waves up to a small cutoff ($l_{\\max}=4$), while the coupling between layers is expanded in the two-dimensional reciprocal lattice vectors $\\vec{G}_{hkl}$, which act as the Umklapp channels that generate Kikuchi lines and bands. With at least 137 $\\vec{G}_{hkl}$ vectors the computed patterns converge in feature positions, and the paper demonstrates agreement with measured total intensity $I_{\\rm RCP}+I_{\\rm LCP}$, the CDAD difference $I_{\\rm RCP}-I_{\\rm LCP}$, and the normalized asymmetry $A_{\\rm CDAD}$ for Si 2p$_{3/2}$ and Ge 3p$_{3/2}$ at $h\\nu=6$ keV, with Ge 2p$_{3/2}$ reproduced qualitatively.","pith_inferences":["The paper validates the single-site-per-layer approximation only at 6 keV; a natural extension is to test the same code at 1–2 keV, where backscattering is stronger, against a full multiple-scattering cluster calculation, since the approximation is most likely to fail there.","If the approximation holds quantitatively, the residual disagreement in the Ge 2p$_{3/2}$ CDAD points to the potential and self-energy choices rather than the scattering geometry; repeating the comparison with an improved potential would isolate that source of error.","The predicted fine Kikuchi lines are sharper than the effective momentum resolution of most measured panels, so a higher-resolution re-measurement of Si 2p$_{3/2}$ could either confirm the predicted line crossings or reveal that the $\\vec{G}_{hkl}$ truncation overproduces fine structure."],"forward_implications":["Hard X-ray PED and CDAD simulations become practical inside a KKR Green's-function code, since the self-consistent potential and the diffraction calculation share the same electronic-structure framework.","The same implementation spans 20–8000 eV, so diffraction patterns can be followed continuously from UV to hard X-ray energies without switching to a cluster model.","Through the alloy-analogy model, finite-temperature and disorder effects can be included in XPD and CDAD simulations, extending the comparisons beyond 0 K perfect crystals.","Because core-level photoelectrons carry element-specific binding energies, the calculated Kikuchi patterns can distinguish chemically different emitter sites in momentum-microscope images.","The computed up-down antisymmetry of $A_{\\rm CDAD}$ and its Kikuchi-grid geometry give a template for interpreting circular-dichroism textures in hard X-ray photoemission experiments."],"supporting_citations":[{"why":"Supplies the relativistic KKR Green's-function framework whose layer scheme the paper extends to photoelectron diffraction.","marker":"[19, 20]"},{"why":"Real-space cluster multiple-scattering package used as an independent cross-check for the angular-momentum convergence tests.","marker":"[18]"},{"why":"Earlier measurement and calculation of W 3d5/2 that motivates the method and provides the CDAD symmetry argument.","marker":"[33]"},{"why":"Dynamical Kikuchi-band simulation approach for high-energy PED that the layer model is contrasted with as a benchmark.","marker":"[22]"},{"why":"Review of the one-step model of photoemission that supplies the time-reversed LEED final-state formalism used here.","marker":"[57]"},{"why":"Alloy-analogy model used to include finite-temperature and disorder effects in the one-step photoemission current.","marker":"[111]"}],"fun_headline_variants":["Layer-KKR model reproduces hard X-ray PED at 6 keV","K-space layer-KKR enables broad-energy hard X-ray PED","Broad-energy XPD without cluster convergence: layer-KKR","One-step layer-KKR reproduces Si and Ge hard X-ray CDAD","Layer-KKR method: broad-energy XPD without convergence issues"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that at multi-keV energies each photoelectron scatters only once inside an atomic layer, with all repeated scattering occurring between layers, so intra-layer multiple scattering can be neglected.","fun_headline_variants_meta":{"raw":{"variants":["Layer-KKR model reproduces hard X-ray PED at 6 keV","K-space layer-KKR enables broad-energy hard X-ray PED","Broad-energy XPD without cluster convergence: layer-KKR","One-step layer-KKR reproduces Si and Ge hard X-ray CDAD","Layer-KKR method: broad-energy XPD without convergence issues"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001666,"raw_usage":{"total_tokens":6674,"prompt_tokens":1069,"completion_tokens":5605,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":685,"completion_tokens_details":{"reasoning_tokens":5513}},"tokens_in":685,"tokens_out":5605,"duration_ms":36622,"temperature":1.0,"reasoning_tokens":5513,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T20:24:29.944534+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the Si 2p$_{3/2}$ 6 keV CDAD patterns with intra-layer multiple scattering included, or with an angular-momentum cutoff near 30, and compare the Kikuchi-band intensities and the $A_{\\rm CDAD}$ maps with the $l_{\\max}=4$ results; if the intensity ratios shift by more than the experimental noise, the single-site-per-layer approximation is falsified.","supporting_citations":[{"cited_title":"Exploring the xps limit in soft and hard x-ray angle-resolved pho- toemission using a temperature-dependent one-step theory","cited_arxiv_id":null,"evidence_quote":"Alloy-analogy model used to include finite-temperature and disorder effects in the one-step photoemission current."}],"review_version":1}