{"id":"a66e6dd9-e02f-4303-9acb-31f7c7623be1","arxiv_id":"2504.14758","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A relativistic one-step photoemission model reproduces fine Kikuchi patterns for Ge(100) and Si(100), including the first observed and simulated Si 1s circular dichroism with up to 31% asymmetry.","lead":"This paper models the fine Kikuchi diffraction patterns created by photoelectrons leaving silicon and germanium crystals, and compares them with experimental images. It reports the first observed and simulated circular dichroism in the angular pattern of silicon 1s core electrons, reaching up to 31% asymmetry.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Agreement may be post-hoc: per-energy G_hkl and V0i are tuned with no convergence criterion or quantitative similarity metric, so the claimed validation of the one-step model and Si 1s CDAD is not yet independent.","rationale":"The reader's weakest assumption identifies exactly the load-bearing issue: the agreement rests on hand-chosen numerical parameters and visual comparison. I agree with that assessment and did not find a more serious internal contradiction. The one-step model itself is physically well motivated, and the qualitative pattern geometry in several figures is genuinely suggestive; independent support appears in prior work by the same group. However, the present paper's strongest claims—first Si 1s CDAD reproduction up to 31% and a robust simulation tool across 106-4174 eV—are not supported by quantitative validation. The per-energy G_hkl variation is particularly concerning because Fig. 2 shows that changing it visibly alters the fine structure, and the final comparison uses 193 G_hkl even though the convergence sequence stops at 145. Similarly, V0i is varied until the blur resembles the experiment, so the inelastic scattering treatment is being used as a fit parameter rather than a tested prediction. These features do not falsify the model, but they do prevent the paper's evidence from independently confirming it. The appropriate outcome remains CONDITIONAL, matching the reader's verdict, with the condition being a quantitative, parameter-stability or hold-out based validation.","tokens_in":17619,"tokens_out":5103,"duration_ms":48628,"concrete_test":"Re-run the Si 1s CDAD comparison for G_hkl in {145, 193, 249} and V0i in {1, 2, 3, 5} eV at both 1440 eV and 4174 eV, then compute a normalized cross-correlation or structural-similarity metric between simulated and experimental I_TOT and A_CDAD images. Perform a simple cross-validation: fit V0i and G_hkl on the 1440 eV patterns only, then measure the metric on the 4174 eV patterns (and vice versa). If the held-out metric is comparable to the in-sample metric and stable across neighboring parameter values, the agreement is predictive; if it peaks only at the chosen parameters or degrades sharply at neighboring settings, the reported validation is tuned to the displayed data.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that the relativistic one-step model reliably reproduces Kikuchi diffraction features and Si 1s CDAD. For this to hold, the experimental agreement must be predictive rather than post-hoc. The manuscript's evidence is entirely visual—'matches almost perfectly' and 'near-perfect agreement'—with no similarity metric, error bar, or hold-out comparison. Two tunable ingredients are set to make the displayed agreement: (i) V0i is explicitly increased from 1 to 5 eV until the simulated blur 'closely resembles' the experimental diffractogram (Fig. 1); (ii) G_hkl is set per energy with no convergence test, and Fig. 2 shows that increasing 45 to 89 to 145 materially adds umklapp channels and changes the fine structure, yet the final Si 1s comparison is run at 193, a value not shown in the convergence sequence. lmax=4 is carried over from earlier work without a new check. Because every visual match is made on the same data used for assessment, a favorable pattern can be found even if the model has systematic errors. The paper itself notes that 'systematic trial tests must be carefully taken into account' and that no single G_hkl works across energies, which concedes the absence of a fixed criterion. This is not an internal contradiction, but it is the load-bearing weak spot in the validation: without a quantitative, transferable criterion, the claimed 'robust tool' and the 31% CDAD reproduction are not independently established.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript presents a computational study of photoelectron Kikuchi diffraction in Ge(100) and Si(100) using the fully relativistic one-step photoemission model implemented in SPRKKR. The authors compute core-level diffractograms for Ge 2p3/2 and 3d and Si 1s and 2p3/2 over kinetic energies from 106 eV to 4174 eV, compare them with momentum-microscopy measurements, and analyze how the imaginary inner potential V0i broadens the patterns and how the G_hkl truncation controls fine structure. The main claims are (i) first observation and successful simulation of circular dichroism in the angular distribution (CDAD) of Si 1s, with asymmetries up to 31%, and (ii) that the one-step model reproduces both bulk and more surface-sensitive Kikuchi diffraction features across a wide energy range, including fine details that earlier cluster-based methods could not capture.","tokens_in":17897,"tokens_out":7171,"duration_ms":63472,"significance":"If the claims are correct, the paper would be a valuable methodological advance: it extends the validated application of the relativistic one-step model to hard-X-ray, emitter-site-specific Kikuchi diffraction, and it reports a new CDAD case for a deep core level. The use of a well-established formalism, the comparison with independent experimental data rather than synthetic inputs, the systematic discussion of inelastic broadening, and the broad energy range are genuine strengths. The main weakness is that the validation evidence is entirely qualitative: agreement is assessed by visual inspection, the numerical parameters controlling agreement are not fully reported, and no convergence criterion is given for the basis-set truncations. As a result, the current manuscript establishes plausibility rather than independent validation. No machine-checked proofs or public code repository are provided; the code availability statement is limited to 'upon request', so the reproducibility of the simulations rests on the completeness of the parameter reporting.","major_comments":[{"comment":"The central claim that the one-step model reproduces the experimental Kikuchi patterns is supported only by visual feature matching. The text uses qualitative statements such as 'quantitatively reasonable' (Fig. 3), 'near-perfect agreement' (Section III B), and 'match almost perfectly' (Section III C) without reporting any numerical similarity statistic. Because the paper's goal is to validate the robustness of the model, the authors should add a quantitative comparison, for example a 2D cross-correlation or a normalized R-factor between measured and computed patterns, together with line-profile comparisons for selected features. Without such a metric, the displayed agreement cannot be distinguished from a favorable visual reading of patterns that share the same overall symmetry, especially since the parameters V0i and G_hkl are chosen with reference to the displayed data.","section":"Section III B, Fig. 3; also Figs. 4–5"},{"comment":"The reciprocal-lattice truncation is not shown to be converged for the final comparisons. Fig. 2 demonstrates that increasing G_hkl from 45 to 145 materially changes the fine structure, but the Si 1s comparison in Fig. 3 is computed with 193 G_hkl, which is not included in the convergence sequence. Fig. 5 similarly uses 45, 69, 97, and 57 G_hkl without convergence evidence. The Methods paragraph acknowledges that 'systematic trial tests must be carefully taken into account' and that no single truncation works for all energies. The authors should either provide a convergence test per reported energy showing that the chosen G_hkl lies in a plateau of a pattern-similarity metric, or present a sensitivity analysis demonstrating that the claimed features do not depend on the truncation.","section":"Section III B, Figs. 2–3; Methods II A"},{"comment":"The imaginary inner potential V0i is the main broadening parameter and is reported only for Fig. 1. The text of Section III A shows that increasing V0i from 1 eV to 5 eV substantially changes the appearance of the diffractograms, yet the simulations shown in Figs. 3–5 do not state the V0i values used. Since the manuscript claims agreement in line broadening and contrast, the revision must report V0i(E) for every calculation, or provide an explicit energy-dependent V0i(E) curve, so that the comparison is reproducible and the role of this parameter in the agreement is transparent.","section":"Section III A and Figs. 3–5"},{"comment":"The central numerical claim of 'CDAD asymmetries up to 31%' is not substantiated in the text or figure captions. No value or location of the maximum A_CDAD is given, and it is not stated whether the 31% figure is experimental or calculated. The authors should report the maximum measured and computed A_CDAD values, their angular coordinates, and the associated uncertainties, and should clarify how the 31% is obtained from the data shown in Fig. 3(g,h).","section":"Abstract and Section III B"}],"minor_comments":[{"comment":"The phrase 'cut from a waver' should read 'cut from a wafer'.","section":"Section II B"},{"comment":"The text states that the polar angle for the azimuthal scan is 13.41 degrees, while the Fig. 1 caption gives 2.68 degrees for the same scans; please reconcile this inconsistency.","section":"Section III A and Fig. 1 caption"},{"comment":"The sentence 'IRCP and IRCP stand for the intensity of RCP and LCP light in turn' should read 'IRCP and ILCP'.","section":"Section III A"},{"comment":"The phrase 'Computational results are performed with 193 ⃗Ghkl' is awkward; a clearer wording would be 'Calculations were performed with 193 ⃗Ghkl vectors.'","section":"Section III B"},{"comment":"The caption entry for EF inal = 1036 eV repeats '(e,f)'; it should refer to panels (g,h).","section":"Fig. 5 caption"}],"recommendation":"major_revision","confidential_remarks":"The paper is a plausible and potentially significant contribution, and I do not see an internal inconsistency that would justify rejection. My recommendation of major revision is driven by the absence of quantitative validation and incomplete parameter reporting, which are fixable within the scope of the manuscript. The revision should be judged on whether the authors add a quantitative similarity metric, convergence tests for G_hkl, and complete V0i reporting, and on whether they substantiate the 31% CDAD claim."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The real news here is the first observation and simulation of Si 1s circular dichroism in photoelectron diffraction, with asymmetries up to 31%, and a systematic energy sweep from 106 eV to 4174 eV across several core levels. That is a concrete step beyond the earlier Si 2s/2p work and the previous model paper, and the one-step layered multiple-scattering approach does avoid cluster-size convergence problems. Credit where due: the broad Kikuchi geometry, band positions, and many fine features are reproduced across multiple energies and emitters, which is not trivial.\n\nBut the validation is weaker than the prose suggests. The paper leans on visual matching—'matches almost perfectly', 'near-perfect agreement'—with no quantitative similarity metric and no error bars. The two main tunable knobs, V0i and the reciprocal-lattice cutoff G_hkl, are adjusted per energy with no convergence criterion. The paper itself concedes that 'systematic trial tests must be carefully taken into account' and that no single G_hkl works across energies. That makes the agreement partly post-hoc: you can often find a blur level and a truncation that makes a simulated pattern look like the data. The Si 1s comparison is even run at G_hkl=193, a value not shown in the convergence sequence of Fig. 2. Thermal effects and Debye-Waller are neglected, which the authors acknowledge, and the calculation uses lmax=4 carried over from earlier work without a fresh check.\n\nNone of this breaks the central claim. The comparison is broad enough and the features specific enough that I do not think the agreement is manufactured. But the paper would be much stronger with at least one quantitative similarity metric (e.g., a correlation coefficient or R-factor) and a stated convergence test for G_hkl. The authors also keep code and data 'upon reasonable request', which makes independent checking harder.\n\nWho is this for? Experimentalists using hard X-ray momentum microscopy who want a reference for what one-step theory can do, and anyone working on Kikuchi diffraction simulation. It deserves a serious referee—the first Si 1s CDAD result is significant enough that a careful review, likely requesting quantitative validation and convergence details, is warranted. I would not desk-reject it; I would send it out with the message that the validation needs hardening.","headline":"First Si 1s CDAD observation plus broad energy-dependent Kikuchi comparison is worth a careful referee; the validation is visual and partly tuned, so the robustness claims should be read with caution.","tokens_in":18528,"tokens_out":1165,"would_cite":true,"duration_ms":12522,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper claims to have observed and, with a fully relativistic one-step photoemission model, reproduced circular dichroism in Si(100) 1s photoelectron Kikuchi diffraction, with asymmetries up to 31%, and to match experimental bulk and…","keywords":["photoelectron diffraction","Kikuchi diffraction","one-step photoemission model","circular dichroism","hard X-ray photoemission","multiple scattering","inelastic scattering","core-level spectroscopy"],"falsifier":"At a fixed kinetic energy (e.g., Si 1s at 4174 eV), recompute the total-intensity and dichroism maps while systematically increasing the number of reciprocal-lattice vectors beyond 193 and the angular-momentum cutoff beyond lmax=4. If the maps change substantively, or if the damping parameter V0i has to be re-adjusted to preserve the match, the claimed agreement is an artifact of the chosen truncations. A second check: measure the Si(100) 1s CDAD at intermediate photon energies and compare the predicted sign and positions of the asymmetry lobes.","tokens_in":17402,"feed_emoji":"🔬","tokens_out":9359,"duration_ms":79833,"temperature":0.7,"pith_summary":"This paper tries to establish that a fully relativistic one-step photoemission model, built on multiple-scattering theory, can reproduce fine Kikuchi diffraction—the band-and-line patterns photoelectrons form by scattering inside a crystal on their way to the detector—from Si(100) and Ge(100) core levels in hard X-ray photoemission. It reports the first observation and simulation of circular dichroism in the angular distribution of Si(100) 1s photoelectrons, with asymmetries up to 31%. It also shows the same calculation tracks how the Kikuchi network changes with kinetic energy from 106 eV to 4174 eV, matching experimental bulk and near-surface features. If the claim holds, the model becomes a practical tool for separating diffraction artifacts from angle-resolved photoemission data and for emitter-site-specific structural studies, avoiding the cluster-size convergence problems of earlier methods.","feed_headline":"Si 1s photoelectron Kikuchi diffraction shows 31% circular dichroism","feed_subtitle":"A fully relativistic one-step model reproduces the faint patterns across 106–4174 eV, making hard-X-ray PES cleaner to interpret.","key_machinery":"The load-bearing machinery is the one-step photoemission model with a time-reversed LEED final state, evaluated through a fully relativistic multiple-scattering Green's-function formalism. The final state is expanded in partial waves up to lmax=4 and coupled across atomic layers by a truncated set of reciprocal lattice vectors G_hkl, with the truncation set between 45 and 193 depending on energy. Inelastic scattering enters through a complex inner potential V0i(E), whose imaginary part attenuates the photoelectron wave field and thus controls how much the Kikuchi bands blur. This machinery replaces cluster-based multiple-scattering calculations, avoiding cluster-size convergence problems and allowing site-specific core-level emission to be simulated at hard X-ray energies.","core_discovery":"The central discovery, stated on the paper's own terms, is that the one-step model reproduces experimentally observed Kikuchi patterns in core-level hard-X-ray photoelectron diffraction across a wide energy range, and that the same calculation captures circular dichroism in the angular distribution (CDAD) of Si 1s emission. For Si(100) 1s, the normalized difference between right- and left-circularly polarized intensities reaches 31%, and computed total-intensity and CDAD maps agree with measured diffractograms. The model also reproduces distinct bulk Si 2p patterns at 3180–3374 eV kinetic energy and more surface-sensitive Ge 3d patterns at 106–1036 eV, where the inelastic mean free path ranges from 0.5 to 2.2 nm. The imaginary part of the inner potential, V0i, is the parameter controlling broadening: increasing it from 1 eV to 5 eV blurs the computed pattern until it resembles the measured one. The paper concludes that the model accounts for the energy-driven transition from surface-sensitive to bulk Kikuchi networks and supports hard-X-ray PES investigations.","pith_inferences":["A systematic convergence study at one fixed energy, varying the number of reciprocal-lattice vectors and the angular-momentum cutoff while measuring pattern change, would show whether the hand-chosen parameters are sufficient or whether missing physics such as thermal motion is being absorbed into the damping parameter.","Because CDAD appears in a non-magnetic, non-chiral material, the same model should predict dichroic asymmetries in heavier 5d or 4f compounds where spin-orbit coupling is stronger, with emitter-site specificity isolating inequivalent lattice sites.","The strong energy dependence of the Kikuchi network implies that multi-energy diffractograms carry depth information; inverting the energy-driven pattern evolution could yield a depth-resolved structural probe by analogy with photoelectron holography.","If the damping parameter is made time-dependent, the model could simulate ultrafast pump-probe photoemission, where photoexcitation transiently changes the inelastic mean free path and alters the diffraction pattern."],"forward_implications":["Core-level Kikuchi patterns can be simulated emitter-site-specifically at hard X-ray energies without cluster-size convergence issues, so the method can address buried interfaces and site-selective structure.","The energy-dependent transition between surface-sensitive and bulk networks gives a practical route to separate bulk from surface contributions in angle-resolved photoemission, cleaning up hard X-ray ARPES data.","Deep core levels such as Si 1s can be used for CDAD-based structural analysis, with predicted asymmetries up to 31% comparable to or larger than those of shallower levels.","Matching computed blur to measured patterns via V0i offers an empirical determination of inelastic attenuation in photoelectron diffraction.","The approach is positioned to extend to site-specific effects in alloys and to time-resolved photoemission after excitation, where inelastic electron-hole generation changes the diffraction features."],"supporting_citations":[{"why":"Supplies the layered multiple-scattering approach to hard X-ray photoelectron diffraction that this paper extends and applies to Si and Ge.","marker":"[12]"},{"why":"One of the few prior successful simulations of hard X-ray photoelectron diffraction in a momentum microscope, used as benchmark for reproducing fine Kikuchi structure.","marker":"[13]"},{"why":"Gives the earlier experimental and theoretical study of circular dichroism in hard X-ray photoelectron diffraction that the present CDAD analysis builds on.","marker":"[15]"},{"why":"Provides the many-beam dynamical Kikuchi-band theory baseline for simulating high-energy photoelectron diffraction.","marker":"[18]"},{"why":"Supplies the experimental bulk Si 2p Kikuchi diffractograms that the paper reproduces in its bulk-sensitive section.","marker":"[19]"},{"why":"Establishes how localized inelastic scattering shapes Kikuchi bands, the basis for the paper's V0i broadening analysis.","marker":"[36]"},{"why":"Supplies model calculations and the reciprocity/fine-structure interpretation used for high-energy photoelectron diffraction.","marker":"[37]"},{"why":"Provides the earlier photoelectron-diffraction theory of circular dichroism in Si(001), which this work extends to the 1s core level.","marker":"[42]"},{"why":"Sets out the fully relativistic one-step description of angle-resolved photoemission that the simulations use.","marker":"[49]"},{"why":"Provides the time-reversed LEED final-state theory used to represent the photoelectron final state in the one-step model.","marker":"[51]"}],"fun_headline_variants":["Si 1s Kikuchi CDAD up to 31% reproduced by one-step model","Hard-X-ray PED: theory tracks Kikuchi patterns from 106 to 4174 eV","Bulk-to-surface Kikuchi transition simulated for Ge(100) and Si(100)","Energy-driven Kikuchi fine structure explained by relativistic photoemission","First CDAD in Si 1s photoelectron diffraction: 31% asymmetry"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the agreement between calculation and experiment is real and not produced by adjusting the numerical truncations and damping parameter to the displayed data for each energy.","fun_headline_variants_meta":{"raw":{"variants":["Si 1s Kikuchi CDAD up to 31% reproduced by one-step model","Hard-X-ray PED: theory tracks Kikuchi patterns from 106 to 4174 eV","Bulk-to-surface Kikuchi transition simulated for Ge(100) and Si(100)","Energy-driven Kikuchi fine structure explained by relativistic photoemission","First CDAD in Si 1s photoelectron diffraction: 31% asymmetry"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000366,"raw_usage":{"total_tokens":2041,"prompt_tokens":1091,"completion_tokens":950,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":707,"completion_tokens_details":{"reasoning_tokens":838}},"tokens_in":707,"tokens_out":950,"duration_ms":8285,"temperature":1.0,"reasoning_tokens":838,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T11:41:38.170793+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"At a fixed kinetic energy (e.g., Si 1s at 4174 eV), recompute the total-intensity and dichroism maps while systematically increasing the number of reciprocal-lattice vectors beyond 193 and the angular-momentum cutoff beyond lmax=4. If the maps change substantively, or if the damping parameter V0i has to be re-adjusted to preserve the match, the claimed agreement is an artifact of the chosen truncations. A second check: measure the Si(100) 1s CDAD at intermediate photon energies and compare the predicted sign and positions of the asymmetry lobes.","supporting_citations":[{"cited_title":"Simulation of high energy pho- toelectron diffraction using many-beam dynamical kikuchi-band theory","cited_arxiv_id":null,"evidence_quote":"Provides the many-beam dynamical Kikuchi-band theory baseline for simulating high-energy photoelectron diffraction."},{"cited_title":"Emitter-site specificity of hard x-ray photo- electron kikuchi-diffraction","cited_arxiv_id":null,"evidence_quote":"Supplies the experimental bulk Si 2p Kikuchi diffractograms that the paper reproduces in its bulk-sensitive section."},{"cited_title":"Influence of localized inelastic scattering on kikuchi bands in photoelectron diffraction patterns","cited_arxiv_id":null,"evidence_quote":"Establishes how localized inelastic scattering shapes Kikuchi bands, the basis for the paper's V0i broadening analysis."},{"cited_title":"High-energy photoelectron diffraction: model calculations and future possibilities","cited_arxiv_id":null,"evidence_quote":"Supplies model calculations and the reciprocity/fine-structure interpretation used for high-energy photoelectron diffraction."},{"cited_title":"Circular dichroism in the angular distribution of core photoelectrons from si (001): A photoelectron- diffraction analysis","cited_arxiv_id":null,"evidence_quote":"Provides the earlier photoelectron-diffraction theory of circular dichroism in Si(001), which this work extends to the 1s core level."},{"cited_title":"Correlation, temperature and disorder: Recent developments in the one-step description of angle-resolved photoemission","cited_arxiv_id":null,"evidence_quote":"Sets out the fully relativistic one-step description of angle-resolved photoemission that the simulations use."},{"cited_title":"Theory of photoemission","cited_arxiv_id":null,"evidence_quote":"Provides the time-reversed LEED final-state theory used to represent the photoelectron final state in the one-step model."}],"review_version":1}