pith. machine review for the scientific record. sign in

arxiv: 2605.00028 · v1 · submitted 2026-04-22 · ⚛️ physics.chem-ph

Recognition: 8 theorem links

· Lean Theorem

A Noble-Gas-Centered Coordinate for Within-Period Atomic Property Trends

Elshad Allahyarov, Jonathan Washburn, Megan Simons

Authors on Pith 3 claimed

Pith reviewed 2026-05-05 22:54 UTC · model claude-opus-4-7

classification ⚛️ physics.chem-ph PACS 31.10.+z31.15.-p32.10.Hq
keywords periodic tableionization energyelectron affinityelectronegativitychemical hardnessnoble-gas coordinategolden ratiowithin-period trends
4
0 comments X

video Unifying the Periodic Table - The Noble-Gas-Centered Coordinate

The pith

A single golden-ratio cosh function on a noble-gas-centered coordinate reproduces within-period trends in ionization energy, electron affinity, electronegativity, and chemical hardness across periods 2–6.

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

The paper proposes that within-period chemistry can be read off one dimensionless landscape: a hyperbolic-cosine function of the noble-gas-centered coordinate ρ = d/Lp, with the golden ratio φ entering as the rescaling factor. From this one function the authors derive proxies for first ionization energy (an outward one-electron step), electron affinity and Pearson hardness (the same inward gap to the period boundary, with two independent scales), and Mulliken electronegativity (the half-sum). They report that this single shape reproduces the full noble-gas-to-alkali ordering of IE1 across periods 2–6, with all 8 of 34 upward deviations sitting on the textbook anomaly sites {p³, d⁵, f⁷, s², d¹⁰}. Two ratio identities — heavy noble-gas IE1 ratios near φ^(1/4) and halogen-to-alkali IE1 ratios near φ² — match NIST data to roughly 1% and 5%. Because EA and η share the same kernel, their ratio collapses to a per-period constant; on period-4 data the nine-atom mean EA/η = 0.180 matches the predicted 0.182 to within 1%, although per-atom scatter is large. The paper presents this as a compact analytical baseline against which shell-specific and relativistic corrections can be measured, not as an ab initio derivation.

Core claim

The authors claim that one dimensionless function on a single periodic-table coordinate — distance to the nearest noble gas, normalized by row length — organizes four central within-period observables at once. The function is J(ρ) = cosh(ρ ln φ) − 1 with φ the golden ratio. Its outward one-electron step serves as the ionization-energy proxy; its inward gap to ρ=1 serves simultaneously as the electron-affinity proxy and the chemical-hardness proxy; Mulliken electronegativity is half their sum. With one fitted energy scale per period, this single shape reproduces noble-gas-to-alkali ordering across periods 2–6, localizes upward bumps exactly at the textbook half-filled and closed-subshell site

What carries the argument

The dimensionless landscape J_chem(ρ) = cosh(ρ ln φ) − 1, evaluated on ρ = (G_p − Z)/L_p ∈ [0,1). Two derived quantities do all the work: an outward step ΔJ⁺ between consecutive ρ values (ionization energy proxy) and an inward gap ΔJ⁻ = J_chem(1) − J_chem(ρ) (shared kernel for both electron affinity and Pearson hardness). One per-period energy scale E_p converts the dimensionless shape to eV; the dimensionless ordering and the EA/η ratio are independent of E_p.

If this is right

  • Cross-period comparison of IE1, EA, electronegativity, and hardness can be done on a single dimensionless axis, with the per-period scale E_p as the only adjustable energy.
  • The shared kernel forces EA/η to be a period-only constant on regular-Aufbau interior atoms, giving a parameter-free benchmark (period-4 mean predicted 0.182, observed 0.180).
  • Deviations from the smooth landscape become diagnostic: anomalies localize at {p³, d⁵, f⁷, s², d¹⁰}, and relativistic shifts (Og, Cn, Ts, Fr) appear as quantified offsets against a defined nonrelativistic reference.
  • Future derivation of E_p from first principles would promote the empirical φ^(1/4) heavy-noble-gas and φ² halogen/alkali ratios from regularities to parameter-free predictions.
  • The single-coordinate kernel currently fails to recover the within-class halogen ordering F > Cl > Br > I, identifying a class-dependent mixing of ΔJ⁺ and ΔJ⁻ as the next refinement.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Because Eq. (14) reduces the φ^(1/4) noble-gas identity to a statement about the ratio E_p/L_p² between consecutive heavy periods, the golden-ratio content is effectively a smoothness condition on the fitted period scales rather than a property of the cosh form itself.
  • The shared EA/η kernel makes a sharper testable claim than the paper foregrounds: the period-mean EA/η should be a fixed function of L_p alone, so periods 5 and 6 ought to reproduce a similar 1%-level agreement if the picture is right.
  • The post-hoc exclusion list (period 1, period 7, halogens in EA, anomaly sites, first-row screening pairs) is doing real work; a principled within-model criterion that predicts which atoms are 'regular-Aufbau interior' would convert the framework from a fit baseline into a selection rule.
  • The φ^(1/4) prediction sets a clean nonrelativistic zeroth-order target for IE1(Og); a future high-precision Og measurement would either confirm the ~7.6% spin-orbit shift assignment or refute the heavy-noble-gas identity.

Load-bearing premise

The choice of the golden ratio φ as the argument-rescaling factor and the cosh form of the landscape is a modeling ansatz rather than something derived from chemistry, so the central identities depend on a kernel selection the authors acknowledge is not unique.

What would settle it

Refit the same observables with alternative smooth kernels — ρ², or cosh(cρ) − 1 with a freely fitted c — on the same period 2–6 NIST data set. If those alternatives match or beat the φ-cosh form on within-period ordering, the φ^(1/4) and φ² ratio identities, and the shared-kernel period-4 EA/η ≈ 0.182 prediction, then the golden-ratio content is empirically unsupported and the central claim collapses.

Figures

Figures reproduced from arXiv: 2605.00028 by Elshad Allahyarov, Jonathan Washburn, Megan Simons.

Figure 1
Figure 1. Figure 1: He 2 Ne 10 Ar 18 Kr 36 Xe 54 Rn 86 Og 118 0 4 8 12 16 Atomic number Z δ(Z) view at source ↗
Figure 2
Figure 2. Figure 2: The dimensionless landscape function Jchem(ρ) = cosh(ρ ln φ) − 1 for a period of length Lp = 8 (periods 2 and 3). Black dots mark the noble-gas (ρ = 0), halogen (ρ = 1/8), p 3 (ρ = 3/8), s 2 (ρ = 6/8), and alkali (ρ = 7/8) positions. The mathematical uniqueness of J is established in Refs.[5,6] and has no chemistry input: it is the unique reciprocal cost on positive ratios under standard smoothness conditi… view at source ↗
Figure 3
Figure 3. Figure 3: Absolute percentage deviation of consecutive noble-gas view at source ↗
Figure 4
Figure 4. Figure 4: Absolute percentage deviation of halogen/alkali view at source ↗
Figure 5
Figure 5. Figure 5: Normalized first ionization energy IEnorm 1 (Z) = IE1(Z)/IE1(Gp) versus displacement ρ, one panel per period (2 blue, 3 red, 4 green, 5 orange, 6 violet, 7 brown). Light fill: regular elements; thick black border: anomaly sites {p 3 , d5 , f 7 , s2 , d10}. Dashed gray: rescaled landscape Jchem(ρ) from Eq. (3). Dotted black: linear guide between the period endpoints. Data: NIST[9] for periods 2–6 and for th… view at source ↗
Figure 6
Figure 6. Figure 6: Electron affinity vs. displacement ρ = d/Lp for periods 4 (red squares), 5 (blue triangles), and 6 (green diamonds), with one-parameter fits from Eq. (7) (solid curves). EA values: NIST. [9] Data and exclusions. The same rules are applied to all three periods. Halogens are dropped because the monotone form of Eq. (7) underestimates the steep shell-closure binding (the period 4–6 halogens lie roughly 2 eV a… view at source ↗
Figure 7
Figure 7. Figure 7: Parity test: NIST Mulliken χM vs. the single-coordinate kernel χstruct (Eq. (9)) for the 15 atoms of view at source ↗
Figure 8
Figure 8. Figure 8: Pearson chemical hardness η = (IE1 − EA)/2 for the non-anomalous elements of periods 2–4 vs. valence position v(Z), with one-parameter landscape fits η = C (p) κRS (Eq. (20)); omitted anomaly points are in view at source ↗
read the original abstract

We introduce a single dimensionless landscape function $J_{\rm chem}(\rho) = \cosh(\rho \ln \varphi) - 1$, $\varphi = (1+\sqrt{5})/2$, on the noble-gas-centred coordinate $\rho = d/L_p \in [0,1)$, and show that it organizes four central atomic observables: first ionization energy \IE$_1$, electron affinity EA, Mulliken electronegativity $\chi_M$, and Pearson chemical hardness $\eta$, on one periodic-table axis. The outward step $\Delta J_{\rm chem}^{+}$ delivers IE$_1$; the inward gap $\Delta J_{\rm chem}^{-} = J_{\rm chem}(1) - J_{\rm chem}(\rho)$ delivers EA and $\eta$; $\chi_M$ follows by Mulliken's identity. Three results establish the empirical content. (i) The within-period IE$_1$ envelope reproduces the full noble-gas-to-alkali ordering across periods 2--6: of 34 atoms compiled across periods 2-4, 26 lie on the predicted monotone descent and the 8 upward deviations occur exactly at the textbook anomaly sites $\{p^3, d^5, f^7, s^2, d^{10}\}$. (ii) Two golden-ratio identities, ${\rm IE}_1(G_p)/{\rm IE}_1(G_{p+1}) \approx \varphi^{1/4}$ on three heavy noble-gas pairs and ${\rm IE}_1(\text{halogen})/{\rm IE}_1(\text{alkali}) \approx \varphi^2$ on four within-period pairs, agree with NIST data to MAD $\approx 1\%$ and $\approx 5\%$, respectively. (iii) The shared kernel $\Delta J_{\rm chem}^{-}$ provides single-parameter analytical fits to EA across periods 4--6 (MAE $0.3$--$0.4$~eV), to Pearson hardness $\eta$ across periods 2--4 (MAE $\sim 1$~eV on noble-gas maxima up to $10.8$~eV), and to Mulliken $\chi_M$ across a 15-atom four-class benchmark ($R^2 = 0.73$). \edit{At the period-averaged scale level, the shared-kernel relation ${\rm EA}/\eta \approx C^{(p)}_{\mathrm{EA}}/C^{(p)}_{\eta}$ is supported on period-4 NIST data: the empirical nine-atom mean $\overline{{\rm EA}/\eta} = 0.180$ agrees with the predicted constant $0.182$ to better than $1\%$, although individual-atom scatter ($\sigma \approx 0.13$) is much larger.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

5 major / 9 minor

Summary. The manuscript proposes a one-parameter analytical framework for within-period atomic property trends. A dimensionless landscape J_chem(ρ) = cosh(ρ ln φ) − 1 is defined on the noble-gas-centred coordinate ρ = d/L_p, with φ = (1+√5)/2. The outward step ΔJ⁺ is assigned to IE₁ and the inward gap ΔJ⁻ is reused as a shared kernel for EA and Pearson η; χ_M follows by Mulliken's identity. Three IE₁ "predictions" (heavy noble-gas ratio φ^(1/4), halogen/alkali ratio φ², and a within-period monotone-with-anomaly envelope across periods 2–6) and three single-parameter fits (EA on periods 4–6, η on periods 2–4, and a 15-atom Mulliken χ_M benchmark with R²=0.73) are reported. A "shared-kernel" benchmark EA/η ≈ C_EA^(p)/C_η^(p) is tested on period-4 data, with predicted ratio 0.182 vs empirical mean 0.180.

Significance. If accepted on its stated phenomenological terms, the paper provides a compact analytical baseline for organizing four standard within-period trends on a single coordinate, with the smooth landscape acting as a backdrop against which textbook anomaly sites and relativistic shifts can be quantified. The within-period IE₁ envelope (Prediction 3) — 26/34 atoms on the predicted monotone descent across periods 2–4 with all 8 upward deviations on the textbook anomaly set {p³,d⁵,f⁷,s²,d¹⁰} — is a non-trivial coordinate-level statement that does not require a fitted scale. The acknowledgement that φ is a modeling ansatz, the explicit per-period exclusion lists, and the auxiliary mathematical uniqueness results (with reported Lean 4 formalization) are strengths. The contribution is presentational/organizational rather than predictive in the strict parameter-free sense; the abstract should reflect this more cleanly.

major comments (5)
  1. [§5.7, Eq. (21) and Table 6] The headline claim that the empirical period-4 mean EA/η = 0.180 matches the 'predicted constant' C_EA/C_η = 0.182 is close to a regression identity rather than an independent test. For OLS-through-origin slopes on overlapping atoms, C_EA/C_η = (Σ EA_i κ_i)/(Σ η_i κ_i) is a κ-weighted mean of EA_i/η_i, while Table 6 reports the unweighted mean of the same nine ratios. Agreement at the ~1% level between a weighted and unweighted mean of nine numbers is largely automatic. The genuine per-atom test of the shared-kernel hypothesis (constancy of EA_i/η_i) fails: σ ≈ 0.13, CV ≈ 70%, range 0.022–0.370. This is acknowledged in the body but the abstract leads with the <1% framing. Please reframe the abstract, and present the per-atom scatter as the primary diagnostic. A more informative test would compare a joint shared-kernel fit (one κ scale) to two independent fits via likelihood ratio or AIC.
  2. [§5.1, Eq. (14)] The paper itself shows that the φ^(1/4) heavy-noble-gas identity reduces, after the cosh small-argument expansion, to (E_p/L_p²)/(E_{p+1}/L_{p+1}²) ≈ φ^(1/4), where E_p is the per-period scale fitted from data (Appendix B). The empirical content is therefore a regularity in fitted scale ratios on three retained pairs (Ar/Kr, Kr/Xe, Xe/Rn) after excluding He/Ne and Ne/Ar. This is stated honestly in §5.1, but the abstract continues to call it a 'golden-ratio identity ... agreeing with NIST data to MAD ≈ 1%.' Please harmonize the abstract and conclusion language with the §5.1 caveat: this is an empirical regularity in E_p/L_p², not a parameter-free derivation.
  3. [§3 (kernel choice)] The choice of φ as the argument-rescaling factor and cosh as the kernel is load-bearing for the title and the framing, yet the paper notes (§3) that ρ² or cosh(cρ)−1 with free c would produce qualitatively similar shapes, and provides no quantitative robustness comparison. Given that the EA, η, and χ_M fits each have a free per-period (or global) scale, alternative kernels are likely to fit nearly as well. Please include a quantitative comparison — e.g., refit EA on periods 4–6 with cosh(cρ)−1 letting c vary, and report MAE/RMSE/R² alongside the φ-fixed fit — so readers can judge how much of the empirical content is specifically 'golden-ratio.'
  4. [§5.4, Eq. (18)] The period-5 EA fit reports R² = 0.07 on n = 8, i.e. the kernel explains essentially none of the variance after exclusions; period 6 has R² undefined and σ_EA comparable to MAE. Aggregating these as 'MAE 0.3–0.4 eV across periods 4–6' in the abstract masks that two of three periods are not actually discriminated from a constant. Please separate the period-4 result (where R²=0.29 is at least non-trivial) from the period-5/6 results in the abstract and conclusion, and state explicitly that the period-5 fit is not better than the period mean.
  5. [§5 (Benchmarking philosophy) and Tables 2, 3, 8] The exclusion set is substantial and is justified kernel-by-kernel (period 1; Ne/Ar; F/Li; period 7; halogens in the EA fit; the anomaly classes {p³,d⁵,f⁷,s²,d¹⁰s²}; relativistic/borderline-Aufbau atoms in period 6). Each exclusion is physically motivated, but the rules differ across observables, and most are introduced in the same section that reports the metrics. Please provide a single, pre-stated exclusion table applied uniformly where possible, and report the metrics on (a) the retained subset and (b) the full data without exclusions, so readers can assess sensitivity. The leave-one-out check in §5.7 is a good model; please extend it to Predictions 1–2 (the φ^(1/4) MAD is computed on three pairs only, so it is sensitive to inclusion of Rn/Og and Ne/Ar).
minor comments (9)
  1. [Abstract] 'golden-ratio cosh coordinate' overstates the parameter-free content given the §3 caveat and Eq. (14). Consider 'one-parameter cosh-form coordinate (golden-ratio rescaling)' or similar.
  2. [Eq. (14)] State explicitly that the φ does not appear on the right-hand side; the empirical claim is (E_p/L_p²)/(E_{p+1}/L_{p+1}²) ≈ φ^(1/4). A one-line remark would prevent the reader from over-reading the identity.
  3. [Table 6] The column 'EA/η − 0.182' with mean −0.002 numerically demonstrates the regression-identity point. Please add a column for the κ-weighted mean of EA_i/η_i so readers can see directly that this equals C_EA/C_η.
  4. [Appendix B, Eq. (25)] Exponents 2.1 and 0.72 are fitted on five noble-gas anchors. State the residual statistics and uncertainties of the fit, and clarify that the pilot has 2 free parameters fit to 5 points.
  5. [Figure 5] The period-7 panel shows Cn above the Og reference (IE_norm > 1). Consider drawing the IE_norm = 1 line explicitly and labelling the violation, since this is discussed in the text as the lone counterexample.
  6. [References / dating] Several references are dated 2026 with future arXiv identifiers (e.g. arXiv:2603.13553, Refs. [5–8]). Please confirm dates and make the supporting-paper status (preprint vs. accepted) clear, since the uniqueness claim used in §3 and Appendix A leans on Refs. [5,6].
  7. [§4.4, Eq. (11)] κ_RS is defined as algebraically identical to ΔJ⁻ from Eq. (7). Two distinct symbols for the same object risk confusing readers; consider unifying or stating the equality at the point of definition.
  8. [§5.5, Table 5 / Figure 7] The χ_struct values for halogens collapse to ~0.060 across F/Cl/Br/I, so the within-class ordering F>Cl>Br>I is not recovered. This limitation is acknowledged, but Figure 7 risks reading as four nearly-coincident points; consider annotating which halogens lie on top of each other.
  9. [§2.1 / Eq. (1)] δ(Z) is introduced for visual purposes and then unused; the body uses ρ(Z)=d(Z)/L_p. Consider deferring δ to a footnote or removing Figure 1, which currently advertises a quantity that does not enter the analysis.

Simulated Author's Rebuttal

5 responses · 2 unresolved

We thank the referee for a careful and constructive report that engages substantively with the empirical content of the manuscript. We agree with the referee's overall framing — that the contribution is presentational/organizational rather than parameter-free predictive — and with the specific diagnosis that several abstract phrases overstate this. We will revise accordingly. The five major comments fall into three classes: (1) abstract/wording corrections that we accept (Comments 1, 2, 4); (2) a substantive request for a quantitative kernel-robustness comparison, which we will add (Comment 3); and (3) a request for a unified, pre-stated exclusion table with full-data sensitivity metrics, which we will add as a new appendix table and extend the leave-one-out analysis to Predictions 1–2 (Comment 5). None of the requested revisions undermine the empirical claims that survive the referee's scrutiny — Prediction 3 (the within-period IE₁ envelope) and the period-4 EA fit (R²=0.29) — but they do require us to reframe the abstract and conclusion in line with §5.1, §5.4, and §5.7 of the body, where these caveats already appear.

read point-by-point responses
  1. Referee: [§5.7, Eq. (21)] The 0.180 vs 0.182 agreement is close to a regression identity (κ-weighted vs unweighted mean of the same nine ratios). The genuine per-atom test (constancy of EA_i/η_i) fails: σ≈0.13, CV≈70%. Reframe the abstract; present per-atom scatter as primary diagnostic; add joint-vs-independent fit comparison (LRT/AIC).

    Authors: We accept this point. The referee's algebraic observation is correct: with OLS-through-origin slopes, C_EA/C_η is a κ-weighted mean of the same ratios that Table 6 averages unweighted, so ~1% agreement on nine numbers is largely structural rather than an independent test. The body of §5.7 already states σ≈0.13 and gives the IQR and leave-one-out range, but the abstract leads with the <1% framing, which we agree is misleading. Revisions: (i) The abstract sentence on EA/η will be rewritten to lead with the per-atom scatter (σ≈0.13, CV≈70%, range 0.022–0.370) and present the period-mean agreement as a secondary, structural observation rather than a test. (ii) The conclusion paragraph will be harmonized. (iii) We will add a model-comparison subsection to §5.7 reporting AIC and a likelihood-ratio test comparing a joint shared-kernel fit (one κ scale shared between EA and η on period 4) against two independent fits, so readers can judge whether the shared-kernel constraint is statistically supported relative to the unconstrained alternative. revision: yes

  2. Referee: [§5.1, Eq. (14)] The φ^(1/4) identity reduces (after cosh expansion) to (E_p/L_p²)/(E_{p+1}/L_{p+1}²) ≈ φ^(1/4), with E_p fitted from data. This is honestly stated in §5.1 but the abstract still calls it a 'golden-ratio identity ... MAD ≈1%.' Harmonize abstract with §5.1 caveat.

    Authors: Accepted. The referee correctly identifies that the empirical content of Prediction 1 — once the small-step expansion of cosh is applied — is the regularity (E_p/L_p²)/(E_{p+1}/L_{p+1}²) ≈ φ^(1/4) on the three retained heavy noble-gas pairs, and that E_p is a per-period scale fitted from data (Appendix B). This is stated in §5.1 but, as the referee notes, the abstract and conclusion phrasing is inconsistent with that caveat. Revision: We will rewrite the relevant abstract sentence to read approximately: 'On the three heavy noble-gas pairs Ar/Kr, Kr/Xe, Xe/Rn, the consecutive IE₁ ratio takes the value φ^(1/4) to MAD≈1%; this is an empirical regularity in the fitted per-period scale ratio E_p/L_p², not a parameter-free derivation.' The conclusion will be amended in parallel. We will retain the language 'golden-ratio identity' only in contexts where the per-period-scale caveat is adjacent. revision: yes

  3. Referee: [§3 (kernel choice)] The choice of φ and cosh is load-bearing for the title and framing, but the paper notes ρ² or cosh(cρ)−1 with free c would give similar shapes, with no quantitative comparison. Refit EA on periods 4–6 with cosh(cρ)−1 letting c vary; report MAE/RMSE/R² alongside the φ-fixed fit.

    Authors: Accepted, and we agree the missing comparison is a substantive gap rather than only a presentational one. We will add a new subsection (§5.4.1, 'Kernel robustness') reporting refits of the EA data on periods 4–6 under three kernels: (a) the present cosh(ρ ln φ)−1 with fixed φ; (b) cosh(cρ)−1 with c a free parameter alongside the per-period scale (i.e., a two-parameter fit per period); and (c) a pure ρ² kernel with one per-period scale. We will report MAE, RMSE, and R² for each, together with AIC penalising the extra parameter in (b). We will also report the best-fit c and assess whether c is consistent across periods and whether c=ln φ is preferred or merely admissible. We expect — consistent with the referee's intuition — that (b) and (c) will perform comparably to the φ-fixed kernel on EA, and we will state this honestly in the revised text. The empirical case for φ specifically rests primarily on the ratio-level Predictions 1–2, not on the EA fit quality, and we will sharpen this attribution in the conclusion. revision: yes

  4. Referee: [§5.4, Eq. (18)] Period-5 EA fit has R²=0.07 (n=8); period-6 R² is undefined. Aggregating as 'MAE 0.3–0.4 eV across periods 4–6' in the abstract masks that two of three periods are not discriminated from a constant. Separate period-4 from period-5/6 in abstract and conclusion; state explicitly that period-5 is no better than the period mean.

    Authors: Accepted. The body of §5.4 already notes 'the period-5 fit is the weakest of the three: R²=0.07 ... comparable to fitting the period mean,' but the abstract and conclusion aggregate the three periods under a single MAE band, which we agree is misleading. Revision: The abstract will distinguish period 4 (R²=0.29, MAE 0.28 eV; the only period where the kernel discriminates from a constant) from periods 5–6, where we will state explicitly that the fit reproduces sign and trend but does not improve on the period mean in R². The conclusion will be amended to match. We will also tabulate period-5 and period-6 'fit vs. constant' MAE side-by-side in the §5.4 diagnostics block so readers can read off the comparison directly. revision: yes

  5. Referee: [§5, Tables 2,3,8] Exclusions are physically motivated but kernel-by-kernel and introduced alongside the metrics. Provide a single pre-stated exclusion table applied uniformly where possible; report metrics on (a) retained subset and (b) full data without exclusions. Extend the leave-one-out check to Predictions 1–2.

    Authors: Accepted. We will add a new appendix table ('Exclusion protocol') listing every excluded atom or pair, the kernel/observable affected, and the physical justification (period 1; first-row screening Ne/Ar and F/Li; period 7 throughout; halogens in the EA fit; anomaly classes {p³,d⁵,f⁷,s²,d¹⁰s²}; relativistic borderline-Aufbau atoms in period 6). The table will be referenced from the start of §5, before the metrics, so the protocol is pre-stated rather than co-introduced with the results. For each fit (EA periods 4–6, η periods 2–4, χ_M 15-atom, EA/η period-4) we will report metrics on both (a) the retained subset and (b) the full unexcluded data, so the sensitivity to the exclusion set is visible. We will also extend the leave-one-out analysis (currently only in §5.7) to Predictions 1 and 2: for Prediction 1 we will report the MAD when Rn/Og is added back and when each of Ar/Kr, Kr/Xe, Xe/Rn is dropped in turn (the φ^(1/4) MAD is computed on three pairs and is, as the referee notes, sensitive); for Prediction 2 we will do the same on the four retained periods. revision: yes

standing simulated objections not resolved
  • The referee's central methodological concern in Comment 1 — that the EA/η <1% agreement is essentially a regression identity — is correct as stated, and no rewriting can convert it into an independent per-atom test. We can add the AIC/LRT model comparison the referee requests, but if it disfavors the shared-kernel constraint we will need to weaken the §5.7 claim accordingly; the outcome of that comparison cannot be guaranteed in advance.
  • The referee's Comment 3 anticipates that cosh(cρ)−1 with free c will fit EA as well as the φ-fixed kernel. If this turns out to be the case quantitatively (which we consider likely), the empirical case for φ specifically will rest entirely on the ratio-level Predictions 1–2, not on the fit-quality results. We accept this consequence and will state it plainly, but it does narrow the load-bearing empirical content of the 'golden-ratio' framing relative to the current manuscript.

Circularity Check

5 steps flagged

Two of three headline "predictions" reduce by the paper's own algebra (Eq. 14) or by overlapping-fit construction (§5.7) to relabelings of fitted per-period scales; the kernel/φ ansatz is imported from same-author work.

specific steps
  1. self definitional [Section 5.1, Eq. (14)]
    "Equation (14) contains no factor of φ explicitly; rather, the cosh expansion reduces the IE ratio to the ratio of Ep/L²_p between consecutive heavy periods. The empirical content of Prediction 1 is the observation that this ratio takes the golden-ratio value... Ep enters here as a per-period scale fitted from data (Appendix B); a closed-form derivation of Ep that would promote this empirical regularity to a parameter-free derivation is left to future work."

    The headline 'φ^{1/4} prediction' for heavy noble-gas IE ratios is shown by the paper's own derivation to be algebraically equivalent to a statement about ratios of fitted per-period scales Ep/L²_p. By Appendix B Eq. (24), E^(G)_p = IE1(Gp)/[cosh((lnφ)/Lp)−1], so Ep is anchored to IE1(Gp) itself. The ratio IE1(Gp)/IE1(Gp+1) reproduces a golden-ratio number iff the data already do. φ has no causal role in P1; the paper itself admits this is empirical, not derived.

  2. fitted input called prediction [Section 5.7, Eq. (21) and Table 6]
    "For period 4 the fitted constants give C^(p=4)_EA/C^(p=4)_η ≈ 7.49/41 ≈ 0.182. ... The empirical period-4 mean EA/η = 0.180 agrees with the predicted shared-kernel constant C^(p=4)_EA/C^(p=4)_η = 0.182 to within ∼1%. This is a period-mean benchmark, not a per-atom proportionality law: the per-atom standard deviation σ ≈ 0.13 (CV ≈ 70%) is much larger."

    C_EA = 7.49 eV (Eq. 17) and C_η = 41 eV (Eq. 20) are linear least-squares slopes of EA and η on the same κ_RS coordinate over essentially the same period-4 atoms. Their ratio equals Σ EA_iκ_i / Σ η_iκ_i, a κ²-weighted mean of EA_i/η_i; 0.180 is the unweighted mean of the same nine ratios. Agreement at 1% between weighted and unweighted means of identical samples is a near-identity of the fitting procedure, not a falsifiable test. The genuine per-atom shared-kernel test is refuted by the paper's own σ ≈ 0.13.

  3. ansatz smuggled in via citation [Section 1, Section 3, Appendix A]
    "The function is the cosh form J_chem(ρ) = cosh(ρ ln φ) − 1, with φ = (1+√5)/2, established in Refs. [5,6] as the unique reciprocal cost on positive ratios under standard smoothness conditions ... (ii) The argument-rescaling factor φ = (1+√5)/2 is presently a modeling ansatz, not a derived chemical constant."

    Refs. [5,6,7,8] are all by the same author group (Washburn, Zlatanović, Allahyarov, Simons, Pardo-Guerra). The cosh kernel's 'uniqueness' is imported from these self-citations to license the chemistry application. The Lean 4 formalization (footnote 1) legitimately certifies the math uniqueness of J on positive ratios but cannot certify that x = φ^ρ is the correct chemical substitution. The paper itself concedes φ is an ansatz and that ρ² or cosh(cρ)−1 with free c would yield similar shapes — so the load-bearing chemical premise rests on a same-author choice.

  4. uniqueness imported from authors [Section 1 footnote 1; Appendix A]
    "The mathematical uniqueness of J is established in Refs. [5,6] and has no chemistry input: it is the unique reciprocal cost on positive ratios under standard smoothness conditions (Appendix A). ¹The uniqueness statement and the φ-power identities used below are also formalized in Lean 4 in the Recognition Science library."

    Mild instance, partially mitigated by Lean 4 formalization. 'Uniqueness' is invoked to forbid alternative kernel forms but is only uniqueness of J on R>0 under chosen regularity conditions, not uniqueness of the chemical embedding ρ ↦ φ^ρ. The paper elsewhere admits alternatives (ρ², cosh(cρ)−1) would 'produce qualitatively similar within-period shapes,' undermining the implied forced-choice rhetoric. Score impact small because machine-checked math is genuine evidence for what it covers.

  5. renaming known result [Section 5.3, Prediction 3]
    "the only upward departures from that monotone sequence sit at the textbook half-filled or completed-subshell sites {p³, d⁵, f⁷, s², d¹⁰}. ... The physics behind the five anomaly sites is textbook and we do not re-derive it here."

    Prediction 3's content is within-period IE descent plus anomaly localization at textbook sites. Both are standard inorganic chemistry; the paper acknowledges this. Novel content is the ρ = d/Lp coordinate organization. Not strongly circular — the ordering test against NIST data is genuinely falsifiable — but the 'predicted' anomaly sites are imported from textbooks rather than derived from J_chem, making the headline closer to organized restatement than derivation.

full rationale

The paper's empirical content rests on three pillars: (P1) heavy-noble-gas IE ratio ≈ φ^(1/4); (P2) halogen/alkali IE ratio ≈ φ^2; (P3) within-period IE ordering with anomaly localization; plus single-parameter EA/η/χ fits and the EA/η ≈ 0.182 "shared-kernel" benchmark. Two distinct circularities are visible in the paper's own equations: (a) Self-definitional via Eq. (14). The paper derives IE1(Gp)/IE1(Gp+1) ≈ (Ep/Lp²)/(Ep+1/Lp+1²); it then states "Equation (14) contains no factor of φ explicitly... The empirical content of Prediction 1 is the observation that this ratio takes the golden-ratio value... Ep enters here as a per-period scale fitted from data." Combined with Appendix B Eq. (24), Ep is anchored to IE1(Gp) itself, so the φ^(1/4) "prediction" is a relabeling of fitted quantities. (b) Fitted-input-called-prediction in §5.7. The "shared-kernel constant" 0.182 is C_EA/C_η = 7.49/41 from two LSQ slopes on the same period-4 atoms; comparing to the unweighted mean EA/η = 0.180 of the same nine atoms is approximately a κ²-weighted vs unweighted average — a regression near-identity. The paper itself reports σ ≈ 0.13 (CV ≈ 70%), undercutting the headline 1% match. (c) Ansatz/uniqueness imported via self-citation. The cosh form and golden-ratio rescaling are licensed by Refs. [5–8], all same author group. A footnote claims Lean 4 formalization, which partially mitigates the math-uniqueness claim but cannot certify x = φ^ρ is the correct chemical substitution; the paper concedes φ is "a modeling ansatz." Mitigating factors: the paper is unusually candid about each of these — explicitly acknowledging Eq. (14) drains φ-content from P1, that φ is an ansatz, that period-5 EA fits explain only 7% of variance, and that EA/η scatter is huge. Prediction 3 (ordering + anomaly sites across 34 atoms vs NIST) is independent of the per-period fits and is genuine empirical content, even though the anomaly sites themselves are textbook. Net: the parameter-free headline claims (φ^(1/4), EA/η ≈ 0.182) reduce to fitted-scale relabelings; one prediction has real independent content. Score 5–6.

Axiom & Free-Parameter Ledger

5 free parameters · 3 axioms · 1 invented entities

The central claim rides on a chosen functional form (cosh) at a chosen scale (φ^ρ) plus several fitted per-period (and per-class) energy scales. The φ choice is admitted to be an ansatz; the cosh form is justified by a uniqueness theorem from the authors' prior work. Anomaly atoms and entire periods are excluded by stated rules. The closest invented-entity is the broader 'Recognition Science' reciprocal-cost framework supplying J; in this paper it is treated as background, but it is not standard in chemistry and lacks chemistry-side independent evidence.

free parameters (5)
  • Per-period IE scale E_p = Ne 1.19e4, Ar 8.71e3, Kr 3.92e4, Xe 3.39e4, Rn 9.51e4 eV (Table 7)
    Fitted from noble-gas IE; required to convert dimensionless landscape to eV; closed-form derivation is explicitly deferred.
  • Pilot fit exponents in E_p^pilot = E_0 L_p^2.1 p^(-0.72) = 2.1, 0.72, E_0 ≈ 2.49e2 eV
    Stated as 'phenomenological, fitted to the five noble-gas E_p^(G) values, and not predicted by RS' (Appendix B).
  • Per-period EA scale C^(p=4,5,6) = 7.49, 9.64, 6.15 eV
    One free scale per period absorbs all EA magnitude information; the dimensionless kernel only supplies shape.
  • Per-period hardness scale C^(p=2,3,4) = ≈68, 48, 41 eV
    One free scale per period for η fits.
  • Per-class and global χ_M scale C = Global 117.4 eV; per-class 92–138 eV
    Fitted to 15-atom benchmark; supplies absolute eV scale for χ_struct.
axioms (3)
  • ad hoc to paper Reciprocal cost J(x) = cosh(ln x)−1 is the unique solution to a functional equation under stated regularity.
    Imported from authors' prior work (refs [5,6]).
  • ad hoc to paper Argument-rescaling factor on the chemical landscape is φ.
    Acknowledged ansatz.
  • domain assumption Anomaly atoms can be excluded from smooth-landscape fits.
    Standard chemistry justification, but used to gate benchmarks.
invented entities (1)
  • Recognition Science reciprocal-cost framework as a chemistry input no independent evidence
    purpose: Supplies the cosh form and motivates the φ-rescaling on which the entire periodic-table parameterization rests.
    Imported from the authors' own institutional program (refs [5–8]) with a claimed Lean formalization; outside that program it has no independent chemistry-side falsification handle. Treated as background here, but readers should note it is not standard atomic-physics machinery.

pith-pipeline@v0.9.0 · 45271 in / 7606 out tokens · 255718 ms · 2026-05-05T22:54:39.291975+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.