{"id":"9275b2df-590a-4e9e-b7c3-7f0396947182","arxiv_id":"2607.07983","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.5,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"State-of-the-art MCDF calculations plus a Lorentzian-Gaussian overlap model reproduce Cu Kα satellite intensity evolution near threshold, attributing below-edge intensity to Cu(I)/Cu(II) oxide resonances.","lead":"Ab initio multiconfiguration Dirac-Fock calculations reproduce the energy-dependent satellite intensity of Cu Kα lines near the K-edge, matching experiment once oxide resonances are added. The work shows that standard atomic codes plus a simple intensity-modulation integral can track the adiabatic-to-sudden transition without empirical shake parameters.","discovery_kind":"extension","skeptic_critique":{"model":"grok-4.5","headline":"Near-threshold agreement and oxide-peak assignment rest on four process-dependent width multipliers plus oxide intensity normalizations fitted to the same sparse Galambosi data set.","rationale":"The Reader correctly isolates the four width multipliers and oxide normalizations as the weakest link. The high-energy sudden-regime ratio is robust and essentially parameter-free, but the central narrative—that standard ab-initio methods already capture the adiabatic-to-sudden evolution and that the below-threshold intensity is cleanly identified as oxide resonances—depends on those fitted scalings. Freezing the multipliers to unity is a direct, reproducible test that either confirms the physical content of the modulation integral or shows that uncontrolled flexibility is required for the claimed agreement. No stronger internal inconsistency is present, so the verdict remains CONDITIONAL.","tokens_in":20902,"tokens_out":555,"duration_ms":5934,"concrete_test":"Re-run the full spectral simulation of Figs. 9–10 with every width multiplier fixed to 1.0 (pure theoretical partial widths) and only the three oxide intensity normalizations free; if the reduced-χ^{2} for the 8975–9050 eV window rises by more than a factor of two relative to the published fit, or if peaks B and C can no longer be placed within 2 eV of the experimental features, the claim that the modulation is essentially parameter-free fails.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The paper’s strongest claim is that parameter-free MCDF rates plus a physically motivated modulation (overlap of effective Lorentzian level profiles with a Gaussian beam, §II G) already reproduce the satellite-intensity evolution, with oxide resonances only filling residual below-threshold intensity. In practice the high-energy µ∞ is indeed close to the ab-initio total-shake probability (0.260 vs 0.262), but the shape that rises from the edge and the three labeled A/B/C peaks are obtained only after nested-sampling optimization of four independent width multipliers (Md≈5, Mo≈7, Mu=1, Me≈0.54) and three oxide intensity normalizations (Table IV, §IV A). These multipliers rescale the very partial widths that enter the modulation integral (Eqs. 11–13), so they are not mere experimental broadenings; they alter the predicted energy dependence of each channel. Because the experimental points are few and carry large uncertainties, the fit can absorb residual solid-state or cross-section effects into the multipliers, making the “good agreement within a good margin of error” and the oxide-origin claim less ab-initio than the abstract asserts.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.5","summary":"The manuscript computes Cu K-shell diagram and satellite transitions with MCDFGME (energies, radiative/non-radiative rates, and sudden-approximation shake probabilities, including a single-shake binomial correction), then builds synthetic Kα1,2 spectra. Near-threshold intensity evolution is modeled by an overlap integral of a Heaviside-modified Lorentzian level profile with a Gaussian beam profile (Eqs. 11–13). High-energy satellite-to-total ratios are compared with Galambosi et al. and with the Thomas and Roy models; residual below-edge intensity is attributed to resonant 1s→3d (Cu(II)) and 1s→4p (Cu(I)) excitations whose strengths and energy shifts are adjusted by nested-sampling fits that also introduce process-dependent width multipliers (Table IV, §IV A).","tokens_in":21317,"tokens_out":1356,"duration_ms":12554,"significance":"A fully ab-initio MCDF treatment of both shake-off and shake-up channels for Cu Kα, together with an explicit single-shake formula and a transparent modulation integral, is a useful contribution to high-resolution X-ray metrology. The high-energy asymptotic ratio µ∞ ≈ 0.262 matches the calculated total shake probability (0.260) and lies inside the experimental uncertainty, confirming that SA shake probabilities remain a good proxy once all decay channels are open. The oxide-resonance assignment is physically motivated by XANES literature and offers a concrete alternative interpretation of the Galambosi below-edge structure. These strengths are real even if the near-threshold shape is not parameter-free.","major_comments":[{"comment":"Abstract and §VI claim that “standard state-of-the-art ab initio methods” already achieve good agreement “within a good margin of error” and that below-threshold intensity “was found to originate” from oxide resonances. In practice the rise shape and the three labeled A/B/C peaks are obtained only after nested-sampling optimization of four independent width multipliers (Md ≈ 5, Mo ≈ 7, Mu = 1, Me ≈ 0.54) and three oxide intensity normalizations (Table IV, §IV A). These multipliers rescale the partial widths that enter the modulation integral (Eqs. 11–13), so they alter the predicted energy dependence of each channel rather than acting as mere experimental broadenings. With sparse, high-uncertainty experimental points the fit can absorb residual solid-state or cross-section effects, making the central claim less ab-initio than stated. The abstract and conclusions should be rewritten to di","section":null},{"comment":"§III and Table I: approximately 2.9 % of the total shake-up probability is missing because only low-n channels were computed; the missing weight is redistributed proportionally among the existing channels before simulation. Because the redistributed intensity inherits the thresholds of the low-n excitations, the energy dependence of the first rise above the K edge is altered by construction. The manuscript should either (i) quantify the change in the simulated ratio when the missing weight is omitted or placed at a continuum threshold, or (ii) demonstrate that the redistribution does not shift the fitted energy thresholds outside their quoted uncertainties.","section":null},{"comment":"§V B and Table V: the oxide species are modeled as free Cu(I)/Cu(II) ions with a single 4p or 3d spectator. The fitted intensity multipliers for peaks B and C reach factors of order 10 (Table IV), which the text attributes either to normalization of the overlap integral or to resonant cross-section enhancement. Without an independent estimate of oxide thickness or of the resonant excitation cross section relative to K-shell ionization, the assignment remains under-constrained. At minimum the paper should report the absolute intensity scale of the oxide contribution relative to bulk Cu and discuss whether a surface oxide of realistic thickness (~100 nm) can produce the observed peak heights.","section":null}],"minor_comments":[{"comment":"Eq. (1) (Thomas model) and the subsequent fit discussion: the numerical prefactor 15.32 is left unexplained; a brief derivation or reference would help readers reproduce the fit.","section":null},{"comment":"Fig. 5 versus Fig. 9: the same experimental points appear with and without oxide contributions; a single multi-panel figure with a clear legend would reduce confusion.","section":null},{"comment":"Table I: the two 2J columns are not defined in the caption; a short note that they refer to the total angular momentum of the Cu ground state would improve readability.","section":null},{"comment":"Several typographical inconsistencies appear (e.g., “ab initio” vs “abinitio”, missing spaces after commas in author affiliations). A careful copy-edit is needed.","section":null},{"comment":"The Roy-model fit (Table II) reports binding energies that differ from theory by many σ for the experimental data set; a short remark on whether this indicates missing physics or simply the limited number of free parameters would be useful.","section":null}],"recommendation":"major_revision","confidential_remarks":"The atomic-structure core is solid and publishable; the over-selling of “parameter-free” near-threshold agreement is the main obstacle. Once the abstract/conclusions are tempered and the sensitivity of the multipliers is shown, the paper becomes a useful contribution. Scope is appropriate for a specialized atomic-physics or X-ray spectroscopy journal; novelty relative to Galambosi et al. (2003) is incremental but real because of the full ab-initio shake-up treatment and the oxide reinterpretation."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"The useful part of this paper is the complete MCDF calculation of Cu Kα diagram + shake-off + shake-up channels, including the single-electron shake formula, and the resulting synthetic spectra. At high energy the satellite ratio lands at ~0.262, essentially the ab-initio total shake probability (~0.260). That asymptotic result is parameter-free and matches both the Thomas model and the Galambosi data within the large experimental error bars. They also give a transparent modulation integral (overlap of an effective Lorentzian level profile with a Gaussian beam) that recovers the adiabatic-to-sudden rise without inventing new physics.\n\nWhat is new is the continuous energy-dependent ratio that separates the individual shake channels and the explicit assignment of the three below-edge peaks (A/B/C) to 1s\to3d (Cu(II)) and 1s\to4p (Cu(I)) oxide resonances. That assignment is physically plausible from XANES literature and improves the low-energy description once the oxides are added.\n\nThe soft spot is real but limited to the near-threshold region. Four process-dependent width multipliers (Md~5, Mo~7, Mu=1, Me~0.54) plus oxide intensity normalizations are optimized by nested sampling against the same sparse, high-uncertainty Galambosi points. Those multipliers rescale the partial widths that enter the modulation integral, so they are not pure experimental broadenings; they absorb residual solid-state and cross-section effects. The high-energy claim stays ab-initio; the detailed shape of the rise and the oxide-peak intensities do not. The oxide species themselves are treated as simple atomic ions, which is a further approximation.\n\nMath and rates look clean, citations are appropriate, and the calculation is fully described even without released code. This is for people who do high-resolution K-edge fluorescence or HERFD-XANES on 3d metals and need a realistic satellite library. It is incremental, not transformative, but it is careful work that a referee should see. I would accept it for peer review and would cite the high-energy shake numbers and the oxide assignment if I were working on similar spectra.","headline":"Solid ab-initio MCDF satellite library for Cu Kα with a clean high-energy limit; near-threshold shape and oxide peaks still need several fitted multipliers.","tokens_in":21848,"tokens_out":554,"would_cite":true,"duration_ms":5921,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.5","headline":"Ab initio calculations of copper K-shell spectra recover the measured rise of satellite intensity from the adiabatic to sudden regime and show that near-threshold intensity comes from oxide resonances, not shake processes.","keywords":["satellite intensity","K-shell transitions","shake probability","sudden approximation","multiconfiguration Dirac-Fock","copper oxides","X-ray fluorescence","near-threshold evolution"],"falsifier":"A new high-resolution, oxide-free copper-foil fluorescence scan that still shows residual intensity between 8975 eV and 8998 eV would falsify the claim that the below-threshold signal is solely oxide resonance; conversely, an oxide-free surface whose ratios match the pure-ionization simulation would confirm it.","tokens_in":21847,"feed_emoji":"⚛️","tokens_out":670,"duration_ms":6477,"temperature":0.7,"pith_summary":"When an X-ray photon knocks a copper atom's deepest electron out of its shell, weaker satellite lines appear beside the main K-alpha emission. Their strength grows as the photon energy climbs above the ionization edge, marking the change from adiabatic to sudden shake-off and shake-up. This paper shows that ordinary multiconfiguration Dirac-Fock calculations of every transition rate and every shake probability, followed by a simple overlap of Lorentzian level widths with a Gaussian beam profile, already reproduce that measured rise within experimental error. Below the metallic copper edge the same calculation finds no intensity; the residual signal matches resonant 1s-to-3d and 1s-to-4p excitations of surface Cu(I) and Cu(II) oxides. The result means that high-resolution fluorescence spectra near threshold can be simulated from first principles without free intensity parameters, and that apparent satellites can be oxide fingerprints rather than multi-electron shake.","feed_headline":"Copper K satellites rise ab initio; near-edge signal is oxide","feed_subtitle":"First-principles rates recover the adiabatic-to-sudden climb; below-threshold peaks match Cu(I)/Cu(II) resonances","key_machinery":"An energy-dependent intensity modulation obtained by numerical overlap of each ionized level's Lorentzian (or Heaviside-modified) natural profile with a Gaussian beam-energy distribution; this factor multiplies the ab initio shake-weighted rates so that channels open continuously from threshold to the sudden-approximation limit.","core_discovery":"Standard state-of-the-art ab initio methods achieve good agreement with experiment and enable simulation of the intensity evolution near ionization thresholds within a good margin of error; the below-threshold satellite intensity originates from resonant 1s to 3d and 1s to 4p excitations in Cu(I) and Cu(II) oxide phases, which were included in the simulations.","pith_inferences":[],"forward_implications":[],"fun_headline_variants":["Ab initio recovers Cu K satellite climb near threshold","Oxide resonances drive below-threshold Cu K satellites","First-principles match Cu Kα intensity evolution at edge","Cu K-shell satellites simulated from ab initio rates","Resonant 1s excitations explain near-edge Cu satellite rise"],"cache_read_input_tokens":16512,"weakest_assumption_plain":"That four process-dependent width multipliers and three oxide intensity normalizations, obtained by fitting the calculated profiles to the published experimental ratios, correctly absorb solid-state and instrumental broadening without introducing uncontrolled bias.","fun_headline_variants_meta":{"raw":{"variants":["Ab initio recovers Cu K satellite climb near threshold","Oxide resonances drive below-threshold Cu K satellites","First-principles match Cu Kα intensity evolution at edge","Cu K-shell satellites simulated from ab initio rates","Resonant 1s excitations explain near-edge Cu satellite rise"]},"model":"grok-4.5","effort":"low","cost_usd":0.005058,"raw_usage":{"total_tokens":1379,"prompt_tokens":708,"num_sources_used":0,"completion_tokens":82,"cost_in_usd_ticks":50580000,"prompt_tokens_details":{"text_tokens":708,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":589,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":708,"tokens_out":82,"duration_ms":5162,"temperature":1.0,"reasoning_tokens":589,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-10T14:14:50.297306+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"A new high-resolution, oxide-free copper-foil fluorescence scan that still shows residual intensity between 8975 eV and 8998 eV would falsify the claim that the below-threshold signal is solely oxide resonance; conversely, an oxide-free surface whose ratios match the pure-ionization simulation would confirm it.","supporting_citations":[],"review_version":1}