{"id":"61bb4dcb-4797-40a0-9154-2e25fd16739d","arxiv_id":"2601.02300","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"The electron's effective magnetic moment in a very small Dirac atom shrinks to ~eR/2, softening the hyperfine attraction to a harmless 1/R form and preventing collapse.","lead":"This paper shows that in a Dirac treatment of hydrogen and positronium, the effective magnetic moment of an electron shrinks to zero with the assumed atomic size, so the normally divergent hyperfine attraction cannot collapse the atom. The result resolves a long-standing puzzle and offers a framework for treating diquarks as relativistic Coulombic systems.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Small-R magnetic-moment softening depends on minimax trial states in a regime the paper admits is not a rigorous upper bound; a fixed-R Dirac calculation is needed to confirm μ_eff≈eR/2.","rationale":"The reader's weakest_assumption identifies the minimax method in the R < 1/m regime as the central uncertainty, and I agree that this is the load-bearing point. The paper's entire resolution of the hyperfine puzzle rests on the small-R behavior of μ_eff, but it is derived from a variational method that the authors concede is not a rigorous upper bound in that regime. The concrete fixed-R calculation I propose would directly test whether the predicted μ_eff ≈ eR/2 is a property of Dirac dynamics or an artifact of the restricted trial family. The proton compositeness concern raised by the reader is real but less decisive: even if Eq. (47) is not derived, the kinematic relativistic suppression of a magnetic moment in a state of size R < 1/M_p is expected to hold for any charged spin-1/2 fermion, and form factors would likely suppress the proton moment further rather than enhance it. Therefore the primary uncertainty is the minimax trial-state reliability, not the proton structure. The paper is otherwise well supported: for R > 1/m it recovers known Dirac hydrogen and positronium results, and the Breit corrections match standard expressions. The conditional verdict is appropriate, and my concern does not move it.","tokens_in":13430,"tokens_out":12887,"duration_ms":144131,"concrete_test":"Perform a fixed-size Dirac variational calculation with a more flexible trial family: upper component f(r) = N1 r^{γ1−1} e^{-r/R} and lower component g(r) = N2 r^{γ2−1} e^{-r/R}, optimizing N1, N2, γ1, γ2, and X for each R. Compute the effective magnetic moment from <e γ0 γ·A> for a uniform field and for a point-dipole field, and compare with eR/2 for R < 1/m. Alternatively, solve the Dirac equation numerically with a hard wall at r = R (or a confining potential that fixes the size) and extract the magnetic moment from the Zeeman splitting in a weak uniform field as R is varied. If the optimized μ_eff(R) still approaches eR/2, the concern is resolved; if it approaches e/2m or shows a different scaling, the paper's softening mechanism fails.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The entire collapse-avoidance argument hangs on Eq. (39): for R < 1/m, the effective electron magnetic moment is claimed to be μ_eff ≈ eR/2, softening the Fermi 1/R^3 hyperfine attraction to 1/R. This result is obtained by extremizing X in the one-parameter trial family (10)-(13) using the minimax procedure. But the paper itself states in Sec. II that the minimax method 'does not in general necessarily provide an exact upper bound' to the energy, and this is precisely the regime R < 1/m where X → 1 and the trial state contains an O(1) lower-component amplitude, i.e., roughly equal admixture of positive- and negative-energy free Dirac components. The energy functional (21) and the matrix element (33)-(39) are then evaluated at a saddle point of the zero-field Hamiltonian, not at a variational minimum. There is no guarantee that a different trial family—e.g., allowing independent radial profiles for upper and lower components, or using a no-pair projected basis—would yield the same μ_eff(R) scaling. Since the conclusion that the hyperfine interaction remains bounded by the kinetic energy depends quantitatively on this scaling in the window 1/M_p < R < 1/m (and for positronium, 0 < R < 1/m), the central claim is not secure without an independent check. The proton effective moment (Eq. 47) is similarly asserted with 'one readily finds' and inherits this concern; however, proton compositeness would likely only strengthen the suppression, so the minimax issue is the more load-bearing gap.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper addresses the 'hyperfine puzzle': why the attractive −1/R³ hyperfine interaction does not cause collapse of hydrogen or positronium. Using the Gordon solution of the Dirac–Coulomb problem and a minimax variational trial function with parameters (R, X, γ), the authors derive the energy of hydrogen (Eq. 21) and positronium (Eq. 30), recovering the known ground-state energies. For the magnetic response, they compute the expectation value of γ⁰γ·eA in the variational state and obtain an effective electron magnetic moment μ_eff/e = R/(2√((mRγ)²+1)) (Eq. 39), which approaches eR/2 as R→0 instead of the free-field e/2m. They argue that this softens the hyperfine interaction from 1/R³ to 1/R, so that the kinetic energy always dominates at small R, preventing collapse. A similar suppression is claimed for the proton (Eq. 47), and the framework is extended to two fermions of arbitrary masses in Appendix C.","tokens_in":13778,"tokens_out":2888,"duration_ms":30778,"significance":"If the central claim is correct, it resolves a long-standing apparent instability in the Dirac–Coulomb description of hydrogen and positronium and provides a new variational tool for relativistic two-body Coulomb systems, with possible applications to diquarks. The paper contains several genuine strengths: it builds on the exact Gordon solution; the variational calculation reproduces the exact hydrogen ground-state energy (Eq. 23) and the known positronium hyperfine result (Eq. 45); and the derivation is analytic and transparent, with no fitted constants. The proposed R-dependent magnetic moment is a concrete, falsifiable prediction. However, the small-R softening result rests on the minimax procedure in a regime the paper itself concedes is not a rigorous upper bound, and the proton effective moment is asserted without derivation. These load-bearing points need independent checks before the collapse-avoidance claim can be considered secure.","major_comments":[{"comment":"The central result μ_eff/e = R/(2√((mRγ)²+1)) is obtained by substituting the energy-extremization condition (19) into the magnetic moment matrix element (33). But in the regime R < 1/m, where the softening matters, the trial state has X→1 and contains an O(1) admixture of lower components, and the paper states in Sec. II that the minimax method 'does not in general necessarily provide an exact upper bound.' The energy functional (21) is evaluated at a saddle point in X, not a variational minimum. A different trial family—for example, allowing independent radial profiles for upper and lower components, or using a no-pair projected basis—could alter the R-scaling of μ_eff. Since the collapse-avoidance claim depends quantitatively on this scaling in the window 1/M_p < R < 1/m (and for positronium, 0 < R < 1/m), an independent fixed-R Dirac calculation or a different variational ansatz is n","section":"Sec. II and Eq. (39)"},{"comment":"The proton effective magnetic moment, μ_eff^p = (g_p/2) e X_p R/(1+X_p²), is introduced with 'one readily finds' and no derivation. This is a load-bearing element for hydrogen stability: the claim that the hyperfine energy is bounded by the electron kinetic energy for all R down to the proton scale requires the proton moment to be suppressed as ~eR for R < 1/M_p. The proton, however, is a composite object with an anomalous magnetic moment, and treating it as a point Dirac fermion down to R ≈ 0.05 fm is not justified. The derivation should be supplied, and the validity of the point-Dirac treatment in this regime should be discussed, including possible form-factor or compositeness effects.","section":"Sec. IV, Eq. (47)"},{"comment":"The extension to two fermions of arbitrary masses rests on the factorized trial wavefunction (C2) with a common spatial factor r^{γ−1}e^{−r/R}. The paper itself notes in Sec. C2 that the presence of a common exponent γ 'has been argued against in a detailed analysis in Ref. [14].' Yet the proton moment and the equal-mass positronium results are tied to this factorization. The heuristic nature of the factorization should be acknowledged in the main text, and the relation to Ref. [14] should be clarified. If Ref. [14] shows that the factorized form is not reliable, that directly undermines the quantitative small-R predictions for hydrogen and positronium.","section":"Appendix C, Eq. (C13)"}],"minor_comments":[{"comment":"Several typos and corrupted symbols appear: 'minumum' (Abstract), 'proceduce' (Sec. II), 'groiund' (Sec. VI), and '⣨1/r⟩' in Eq. (14) should be ⟨1/r⟩. The reference in the bibliography appears corrupted ('Ko/suppress lakowska').","section":"Throughout"},{"comment":"The Breit factors B_o (Eq. 36) and B_i (Eq. 42) are introduced and then set to unity in the subsequent equations; the order of the approximation should be stated more precisely. Also, Eq. (41) for the hyperfine interaction in hydrogen contains an explicit 1/R³ factor before the X-dependent suppression, which is clear but would benefit from a sentence emphasizing that the R³ divergence is cancelled by the R in Eq. (47).","section":"Sec. IV"},{"comment":"The comparison with the Landau–Lifshitz bound for an r^{-2} potential is used to argue that the hyperfine potential has no bound states; however, the condition μ_f < μ_e/8 is quoted without derivation. A brief derivation or reference to the exact condition would improve self-containedness.","section":"Sec. V"}],"recommendation":"major_revision","confidential_remarks":"The paper is thought-provoking and the electron part is elegant, but the central collapse-avoidance claim currently hinges on an unverified small-R variational regime and an underived proton result. These are not merely presentation issues; they are load-bearing. I would encourage the authors to add a fixed-R numerical or analytic check of Eq. (39) with a more general trial family, and to derive Eq. (47) explicitly. The speculative diquark remarks are not yet developed enough to weigh heavily, but they do not affect the core assessment."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague, you should know this paper is worth reading but is not yet the last word. The central claim—that for R < 1/m the effective electron magnetic moment becomes eR/2, softening the Fermi 1/R^3 hyperfine attraction to 1/R—is a genuinely new result. The derivation via the minimax variational method from the Gordon solution is elegant, and the paper recovers the exact hydrogen ground state and the known positronium hyperfine shift (Eq. 45), which gives me confidence the machinery is sound.\n\nThe best part is that the effective moment formula (39) is not a fitted curve; it follows from extremizing the same trial wavefunction used for the energy. That self-consistency is a real strength. The extension to unequal masses in Appendix C, while heuristic, gives a simple and plausible energy formula that reduces correctly to hydrogen and positronium limits.\n\nNow the soft spots. The reader and the stress-test note are right to worry about the small-R regime. The minimax method is not a rigorous upper bound there—the paper admits this in Sec. II—and the trial family with a single exponent γ may not faithfully represent the Dirac dynamics. If a different trial family changes the μ_eff(R) scaling, the collapse-avoidance argument weakens. This is not a fatal objection because the electron part is anchored to the exact solution for all R, but an independent fixed-R Dirac calculation would settle it.\n\nThe proton effective moment, Eq. (47), is asserted with “one readily finds” and no derivation, and treating a proton as a point Dirac fermion down to 0.05 fm is physically dubious. The authors note proton compositeness likely strengthens the suppression, which is plausible, but a referee should push for the derivation. The two-fermion common-exponent assumption in Appendix C is explicitly disputed by Ref. [14]; the paper acknowledges this but proceeds anyway. That is acceptable for a variational starting point, but it limits the confidence in the diquark framework.\n\nFootnote 1’s pseudopotential explanation also deserves reconciliation with the new Dirac-moment picture. The two mechanisms might be complementary, but the paper leaves the connection loose.\n\nOverall: this is a serious paper with a clear new result and an honest set of caveats. It deserves peer review, not a desk reject. I would send it back with requests for the proton derivation and an explicit discussion of the trial-family dependence. For my own work, I would not yet cite the small-R scaling as established, but the positronium result and the variational method are useful.\n\nRecommendation: accept for review, expect major revision.","headline":"A plausible variational resolution of the hyperfine collapse puzzle, clean for the electron and positronium, but the proton part and the small-R regime need more work.","tokens_in":14282,"tokens_out":2581,"would_cite":false,"duration_ms":26048,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["31.30.jf","32.10.Fn"],"model":"deepseek-v4-flash","headline":"This paper claims that a bound fermion's effective magnetic moment shrinks to eR/2 once an atom's assumed radius R falls below the Compton wavelength, softening the hyperfine attraction from 1/R^3 to at most 1/R and resolving why hydrogen a","keywords":["hyperfine puzzle","Dirac equation","minimax variational method","effective magnetic moment","Compton wavelength","positronium","two-fermion bound states","diquarks"],"falsifier":"Solve the two-body Dirac or Salpeter equation for a hydrogenic or positronium-like state constrained to R = 0.1 ħ/mc and compute the magnetic response: if the effective moment remains e/2m rather than falling to roughly eR/2, the predicted softening is an artifact of the trial family. A scan of alternative trial wavefunctions (e.g., r^ν e^{-r/R} with ν ≠ γ-1) would reveal whether the small-R scaling is robust.","tokens_in":13269,"feed_emoji":"⚛️","tokens_out":5414,"duration_ms":57248,"temperature":0.7,"pith_summary":"This paper claims that the apparent 1/R^3 hyperfine attraction that would naively make hydrogen and positronium collapse is an artifact of treating the fermions' magnetic moments as fixed. In a full Dirac treatment, the effective magnetic moment of an electron (or proton) in a trial state of radius R is not e/2m but eR/(2√((mRγ)²+1)); once R drops below the Compton wavelength, the moment itself shrinks to about eR/2. As a result the hyperfine term scales at worst as -1/R, never outrunning the +1/R kinetic energy, so the atom is stable. The same mechanism is shown to apply to positronium and to any two-fermion Coulombic bound state, and the authors give simple variational energy formulas for these systems.","feed_headline":"Magnetic moments shrink at small R, saving hydrogen","feed_subtitle":"A Dirac variational calculation shows the effective moment falls to eR/2 below the Compton wavelength, so hyperfine attraction is bounded by","key_machinery":"The engine is the minimax variational principle applied to the Dirac equation: for the hydrogenic trial state (φ,χ) ∝ (1, iX σ·r̂) r^{γ-1} e^{-r/R}, one first extremizes the energy with respect to the lower-component mixing parameter X—which yields a maximum, not a minimum—and then minimizes with respect to the radius. This reproduces the exact hydrogenic ground state at the physical radius. When a magnetic field is added, the same X controls the magnetic response, and the identity X/(1+X²) = 1/(2√((mRγ)²+1)) converts the expectation value of the electron's magnetic coupling into the R-dependent effective moment of Eq. (39).","core_discovery":"The central result is that the effective magnetic moment of a fermion in a Coulombic bound state of assumed radius R is μ_eff/e = R/(2√((mRγ)²+1)) (with γ=√(1-α²) for hydrogen), so that for R ≪ 1/m the moment tends to eR/2 instead of its free-field value e/2m. When the radius is at its energy minimum the moment returns to e/2m, but for smaller R the 1/R^3 divergence of the hyperfine energy is replaced by at most a 1/R interaction. Applied to the proton (g_p=5.58), the same suppression gives an effective proton moment of g_p eR/4, so the hyperfine energy stays below the electron kinetic energy at all R. For positronium, the hyperfine energy becomes approximately 2α/(3((mRγ)²+1)R), likewise bo","pith_inferences":["If the mechanism generalizes to any Coulomb-like confining interaction, it predicts that tightly bound color-electric diquarks do not collapse in spin-singlet states even for large effective α; the paper's stability criterion, 2 > α(1 - 2⟨σ₁·σ₂⟩/(3(2γ-1))), could be tested in lattice QCD for (qq) diquarks.","One could test the claim experimentally by probing the hyperfine splitting of hydrogen under strong external compression or in a cavity that restricts the effective radius: the prediction is a strong suppression of the hyperfine coefficient once the confinement scale drops below the electron Compton wavelength.","Because the minimax method lacks a rigorous upper bound for R<1/m, the effective-moment scaling is a property of the chosen variational family rather than a proven operator identity; a direct Bethe-Salpeter or lattice calculation of the magnetic response of a Coulombic pair at R<1/m would settle whether the eR/2 softening is exact.","The same R-dependent moment suppression could influence interpretations of the proton-radius puzzle: in very compact states the Dirac proton moment would be suppressed, making the hyperfine contribution to the Lamb shift smaller than the naive Fermi estimate."],"forward_implications":["Hydrogen, muonium, and positronium are stable against hyperfine collapse at all scales: the hyperfine energy is bounded by kinetic energy even as R→0.","The missing delta-function bound states of the contact hyperfine pseudopotential are explained: at sub-Compton radii the effective moment shrinks, removing the singular attraction.","The Dirac variational method gives analytic positronium ground-state energies in agreement with existing numerical variational Coulomb two-body calculations.","The same formalism extends to arbitrary-mass two-fermion Coulombic atoms, yielding E = √(m_1²-1/R²)+√(m_2²-1/R²).","For diquarks modeled as relativistic Coulombic systems with color-electric and color-magnetic interactions, the analogous boundedness suggests that spin-singlet diquark configurations do not collapse at small size."],"fun_headline_variants":["Magnetic moments shrink at small R, preventing collapse","Effective magnetic moment falls to eR/2, stabilizing atoms","Dirac solution saves hydrogen from hyperfine collapse","Suppressed magnetic moments explain hydrogen's stability","The hyperfine puzzle: why small atoms don't collapse"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The load-bearing premise is that the minimax variational trial states with R smaller than the Compton wavelength faithfully represent the Dirac dynamics of the bound fermions; the paper states that in this regime the method does not in general provide an exact upper bound, so the predicted softening rests on variational states rather than a proven bound.","fun_headline_variants_meta":{"raw":{"variants":["Magnetic moments shrink at small R, preventing collapse","Effective magnetic moment falls to eR/2, stabilizing atoms","Dirac solution saves hydrogen from hyperfine collapse","Suppressed magnetic moments explain hydrogen's stability","The hyperfine puzzle: why small atoms don't collapse"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001762,"raw_usage":{"total_tokens":6834,"prompt_tokens":828,"completion_tokens":6006,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":572,"completion_tokens_details":{"reasoning_tokens":5930}},"tokens_in":572,"tokens_out":6006,"duration_ms":37417,"temperature":1.0,"reasoning_tokens":5930,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T12:32:28.543506+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Solve the two-body Dirac or Salpeter equation for a hydrogenic or positronium-like state constrained to R = 0.1 ħ/mc and compute the magnetic response: if the effective moment remains e/2m rather than falling to roughly eR/2, the predicted softening is an artifact of the trial family. A scan of alternative trial wavefunctions (e.g., r^ν e^{-r/R} with ν ≠ γ-1) would reveal whether the small-R scaling is robust.","supporting_citations":[],"review_version":1}