{"id":"09ec2c99-dc0a-4d61-a02d-c0f1deb17b0e","arxiv_id":"2501.00010","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Hydrogen eigenfunctions can be obtained by differentiating harmonic 4D polynomials and substituting t=-ir, giving a coordinate-space version of Fock's theory.","lead":"Hydrogen-atom wavefunctions can be generated from four-dimensional harmonic polynomials by differentiation and a time-to-radius substitution, placing Fock's momentum-space symmetry into ordinary coordinate space. The paper also applies this tensor method to re-derive the Stark effect, the Schwinger resolvent, and ladder operators.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The §VI.D 'discard delta functions' step is asserted, not proved: for the polynomial numerator, the contour integral (64) is not convergent, and no general distributional lemma shows that the residue value (65) is the physical Fourier inverse rather than one of many regularizations.","rationale":"The reader's weakest assumption is exactly the load-bearing formal step: the inverse Fourier transform of Eq. (62) is claimed to reduce to the residue at τ = −ir after delta functions are discarded. I agree that this is the point on which the central algebraic recipe depends. The paper verifies the formula on examples and Appendix D checks a polynomial identity, but neither supplies a general distributional justification of the truncation. Because the polynomial numerator makes the contour integral ill-defined as written, the gap is real; however, it is likely repairable by a standard tempered-distribution argument, so it warrants a conditional verdict rather than rejection. The internal Stark inconsistency (Eq. 43 vs Eq. 49) is secondary and does not affect the main claim, but it strengthens the need for a careful revision. The proposed numerical check would settle whether the truncation is merely under-proved or actually wrong for a nontrivial state.","tokens_in":30901,"tokens_out":19150,"duration_ms":212544,"concrete_test":"Compute for n = 3, l = 0 (k = 2) the left side of Eq. (63) as lim_{ε↓0} ∫ e^{−iτ} P(τ) e^{−ε|τ|} / (τ² + r²) dτ, with P(τ) obtained from the bracket in Eq. (62), and compare at r > 0 with the right side of Eq. (66) using the known R_{30}(r) up to the paper's normalization. If the regularized integral equals Eq. (66), the delta discarding is validated for a nontrivial case; if it differs by terms such as r^{-1} or r^{-3}, the central recipe is incomplete. An independent check is to re-derive Eq. (67) by performing the Fourier integral before differentiating, instead of after.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section VI.D (Eqs. 63–66) is the hinge of the paper. The function in Eq. (62) is a polynomial in τ divided by τ² + r², times e^{−iτ}; its Fourier transform at ω = 0 is a tempered distribution, not an ordinary contour integral. The claim that 'the result is a combination of delta functions, which must be discarded' is not a derivation: closing the contour in the lower half-plane is only allowed if the integrand vanishes on the large semicircle, and it does not when the Gegenbauer numerator has degree k ≥ 2. A correct statement would be that the delta terms are supported at ω = 1, hence vanish in a neighborhood of ω = 0, but this is not written down. Footnote 11 gestures at choosing a polynomial to remove the deltas but does not prove that this choice preserves the physical boundary condition or is equivalent to the original Fock normalization. Appendix D verifies a polynomial identity after the residue substitution, not the truncation itself. Thus the central algebraic map (65)–(66) rests on an unproved distributional regularization. This is a correctness risk, not merely a style issue, because a different regularization can change the result by local terms at r = 0.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a coordinate-space reformulation of Fock's momentum-space treatment of the hydrogen atom. The author extends the Fock eigenfunctions harmonically into a fourth dimension, applies an inverse 4D Fourier transform, and claims that the Schrödinger equation becomes the 4D Laplace equation in a coordinate space with a fictitious time-like coordinate. The central result is an algebraic recipe, Eqs. (65)–(66), that generates hydrogen eigenfunctions from derivatives of harmonic 4D polynomials after the substitution τ = -ir. The paper also derives a differential equation in momentum space, a compact quadratic Stark-effect calculation, the Schwinger resolvent, and vector ladder operators using harmonic-tensor methods.","tokens_in":31183,"tokens_out":12320,"duration_ms":131977,"significance":"If the central map is valid, the paper gives a genuinely coordinate-space realization of Fock's SO(4) symmetry and an efficient algebraic route to hydrogen wavefunctions, with no fitted parameters. The tensor identities in Sec. III and the example checks in Sec. VI and Appendix D are explicit and checkable, and the Stark formula provides a concrete quantitative prediction. The main obstacle is not circularity: the derivation does not fit parameters to the known eigenfunctions. The obstruction is the unproved distributional truncation in Sec. VI.D, which supports the principal claim. The result is therefore currently conditional rather than established.","major_comments":[{"comment":"The step \"the result is a combination of delta functions, which must be discarded\" is the central load-bearing step of the paper, and it is asserted rather than proved. The integrand in Eq. (64) has a polynomial numerator, so the Fourier transform is a tempered distribution; for a Gegenbauer numerator of degree k ≥ 2 the integral over τ is not absolutely convergent, and the large-semicircle contribution in the lower half-plane does not vanish because the endpoint values of e^{-iτ} do not decay. A valid argument must either compute the distributional Fourier transform of P(τ)e^{-iτ}/(τ²+r²) and show that the polynomial contributions are supported at ω = -1 (with the convention of Eq. (63)) and hence do not affect the value at ω = 0, or introduce an explicit regularization and prove that the residue value is the unique physical one. Footnote 11 gestures at removing the delta terms but does not prove that the choice preserves the physical boundary condition or the Fock normalization. Until this is supplied, Eqs. (65)–(66) are not established.","section":"Sec. VI.D, Eqs. (63)–(66)"},{"comment":"The verification in Appendix D applies the residue substitution and then proves a polynomial identity for the resulting Laguerre/Gegenbauer expressions. This confirms the internal consistency of the recipe for individual n, l, k but does not address the truncation step itself. In particular, it does not exclude regularizations that add local terms supported at r = 0, which would change the boundary condition of the claimed physical eigenfunction. A complete proof of Eqs. (65)–(66) must include the missing distributional lemma, not only the post-substitution polynomial identity.","section":"Appendix D, Eqs. (D.1)–(D.10)"}],"minor_comments":[{"comment":"The claimed agreement between Eq. (49) and Eq. (43) for m = n-1 should be displayed explicitly. If the numerator in Eq. (49) is (n+1)(4n+5), the two expressions coincide; if it is (n-1)(4n+5), they do not. The current typography is ambiguous, and a one-line substitution check would remove the ambiguity.","section":"Sec. IV.C, Eq. (49)"},{"comment":"Eq. (63) defines an inverse Fourier transform in ω, but the following paragraph evaluates the integral at ω = 0 without stating that this is the physical zero-frequency component. This convention should be stated explicitly, especially because the distributional support of the discarded terms depends on the Fourier convention.","section":"Sec. VI.D, Eq. (63)"},{"comment":"The manuscript contains numerous transcription-level errors and inconsistent notation, for example \"Decomcoordinate\" in Sec. III.E, \"Acknoledgements\" in the contents, and \"Feinman\" in Ref. [58]. These do not affect the mathematics but make the paper harder to referee and use.","section":"Throughout"},{"comment":"The constants in Eqs. (60)–(66) are described as \"floating\" and the normalization is left unspecified. Since Eqs. (65)–(66) are intended as a calculational recipe, a fixed normalization convention (or an explicit statement that all formulas are up to an n,l-dependent constant) should be stated before the main formula.","section":"Sec. VI.C"},{"comment":"The statement that the substitution τ = -it means a transition to a wave equation with speed 2Ze/n is made without derivation; a one-line derivation from the 4D Laplace equation would help the reader assess this physical interpretation.","section":"Sec. X"}],"recommendation":"major_revision","confidential_remarks":"The central distributional gap is fixable in principle, but the manuscript in its current form does not supply the needed lemma. I recommend major revision rather than rejection because the examples and Appendix D suggest the recipe is correct, and the missing proof can be added without changing the paper's scope. The editor may also wish to ask the author to clarify the status of the momentum-space differential equation in Sec. V, which is presented as a new result but is not used in the main algebraic map."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nQuick take: this paper gives a genuinely clean coordinate-space route from 4D harmonic polynomials to hydrogen wavefunctions, but the hinge of the derivation, the \"discard delta functions\" step in Section VI.D, is not rigorous as written. The paper deserves a serious referee, not a desk reject, and not immediate acceptance.\n\nWhat is actually good: the tensor-calculus apparatus in Sections III and VIII is well built and does real work. The Stark-effect calculation and the Schwinger-resolvent derivation show that the harmonic-polynomial method can compress fairly messy standard computations. Appendix D's Laguerre-Gegenbauer identity, Eq. (D.10), looks genuinely new and is hand-verifiable; the example checks in the body pass. The paper does not fit parameters or normalize its way to the conclusion: it starts from Fock's theory and recovers the standard eigenfunctions. The reliance on the author's earlier tensor papers is not a problem; those identities are standard enough and are re-derived here.\n\nThe soft spot is exactly where the stress-test note puts its finger. Eq. (64) is not a convergent contour integral: the polynomial numerator prevents closing the large semicircle. The delta contributions from the polynomial part are supported away from the frequency at which the inverse transform is being evaluated, so the discarding may well be correct, but that is a distributional lemma the paper never states. Footnote 11 gestures at removing the deltas by adding a polynomial, but it does not prove that the chosen polynomial preserves the physical boundary condition or that the residue value is the true Fourier inverse rather than one regularization among many. Appendix D verifies the algebra after the residue substitution, not the substitution itself. This is a correctness gap, not a style issue, because a different regularization could add local terms at r = 0.\n\nOne secondary point: the reader's claim of an internal inconsistency between Eq. (43) and Eq. (49) does not land cleanly. On the text as typeset, Eq. (49) appears to intend n^4(n+1)(4n+5)/4, which does match Eq. (43) when m = n-1. If the minus sign is literal, the paper is simply wrong there, but given the OCR quality I would not build a rejection on it.\n\nNo new physics is claimed or delivered, and the paper should say more plainly that it is a reformulation. That is not a fatal objection; a genuinely simpler derivation of known results can be worth publishing.\n\nBottom line: send it to peer review, but ask the referee to demand a precise distributional justification for Section VI.D or a downgraded claim that the map is verified for all n,l by Appendix D plus examples. If the author can supply that lemma, I would be comfortable with the paper.","headline":"An elegant coordinate-space repackaging of Fock's SO(4) treatment of hydrogen, but the central inverse-Fourier step in Section VI.D is asserted rather than proved; worth refereeing, not accepting as-is.","tokens_in":31659,"tokens_out":6787,"would_cite":false,"duration_ms":76597,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["02.30.Em","03.65.Db","03.65.Ge"],"model":"deepseek-v4-flash","headline":"The paper claims that in a 4D coordinate space the hydrogen Schrödinger equation becomes the 4D Laplace equation, and that physical eigenfunctions are generated by differentiating harmonic 4D polynomials and setting the extra coordinate…","keywords":["Fock's theory","quantum Coulomb problem","harmonic polynomials","harmonic tensors","SO(4) symmetry","Fourier transform","ladder operators","Stark effect"],"falsifier":"Take a state beyond the paper's worked examples, such as $n=3$, $l=0$, and evaluate the right-hand side of Eq. (66) with and without the delta-function contributions generated by the polynomial numerator. If the two results differ from the standard hydrogen $3s$ wavefunction, the claimed algebraic map is incomplete; if the truncated expression reproduces it, the truncation rule is supported at least for that case.","tokens_in":30698,"feed_emoji":"⚛️","tokens_out":8551,"duration_ms":82493,"temperature":0.7,"pith_summary":"Fock showed that the hydrogen atom's hidden SO(4) symmetry becomes manifest when momentum space is wrapped onto a 3-sphere. This paper tries to establish that the same symmetry can be realized in an ordinary four-dimensional coordinate space, where the Schrödinger equation becomes the 4D Laplace equation and the eigenfunctions are harmonic polynomials. The return to physical space is claimed to be algebraic: differentiate the harmonic polynomial with respect to the extra coordinate, multiply by $e^{-r}$, and then set that coordinate to $t=-ir$. If this recipe is correct, hydrogen wavefunctions and derived quantities such as Stark shifts, the Schwinger resolvent, and ladder operators can be computed without integrals or stereographic projection. The paper also claims a new special-function identity connecting Laguerre polynomials to derivatives of Gegenbauer polynomials.","feed_headline":"Hydrogen eigenstates reduced to 4D harmonic polynomials","feed_subtitle":"A complex-time coordinate turns hydrogen eigenstates into derivatives of 4D harmonic polynomials.","key_machinery":"The load-bearing objects are harmonic 4D polynomials, written invariantly as traceless symmetric tensors whose contraction over any two indices vanishes. The Gegenbauer polynomial $C_k^{l+1}(i t/R)$ carries the radial structure, and the operator $\\partial_t^{n-1}$ converts the polynomial degree into the principal quantum number $n$. The substitution $t=-ir$ after differentiation projects the 4D harmonic polynomial onto physical space; the paper interprets it as passing from the 4D Laplace equation to a wave equation with speed $2Ze/n$. A four-dimensional raising operator $\\hat D$ steps between multipole ranks and underlies the ladder-operator results, while trace-reduction formulas make perturbation-theory contractions algebraic.","core_discovery":"The paper's claim is that Fock's momentum-space construction can be replaced by a coordinate-space construction. Starting from solid 4D spherical functions, that is, harmonic polynomials of degree $n-1$, it performs a 4D Fourier transform, discards the delta-function terms coming from the polynomial numerator, and closes the integration contour to obtain $\\Psi_{nl}(\\mathbf{x})\\propto e^{-r}\\,\\partial_t^{n-1}\\left[Y_l(\\mathbf{x}) C_k^{l+1}(i t/R) R^{n-1}\\right]\\big|_{t=-ir}$. The squared 4D radius $R^2=r^2-t^2$ is set to zero only after differentiation. In this picture the SO(4) symmetry lives in the harmonic tensor itself, and the substitution $t=-ir$ is what hides the symmetry in the familiar coordinate wavefunctions. The same machinery yields a differential equation in momentum space, a compact derivation of the quadratic Stark effect, an electrostatic rederivation of Fock's integral equation, and the Schwinger resolvent as a Gegenbauer-polynomial series.","pith_inferences":["If the delta-function truncation can be made rigorous, the construction would give a fully algebraic derivation of the hydrogen spectrum and a compact route to matrix elements.","The same pattern, inversion plus Fourier transform instead of stereographic projection, might extend to other quantum systems whose hidden symmetry is known in momentum space, such as Dirac hydrogen or dynamical-symmetry problems.","The Appendix D identity could be tested numerically for arbitrary $n,l,k$; if it holds beyond the examples checked, it is an independent generating-function result.","The wave-equation interpretation of $t=-ir$ hints at a time-dependent formulation of Coulomb scattering, though the paper does not develop that direction."],"forward_implications":["Hydrogen eigenfunctions can be constructed by differentiating harmonic 4D polynomials: no Fourier integrals or spherical-function expansions are needed in the final step.","The SO(4) symmetry of the Coulomb problem is realized in a coordinate space whose extra coordinate acts as complex time; after $t=-ir$ the symmetry is concealed in standard solutions.","The quadratic Stark effect for states without a linear effect follows from tensor contractions and Euler's theorem, yielding the dipole-moment formula without parabolic-coordinate separation.","Fock's integral equation can be rederived from electrostatic boundary conditions on spheres in 3D and 4D, and the Schwinger resolvent is obtained as a series in Gegenbauer polynomials.","The polynomial-correspondence identity in Appendix D connects Laguerre polynomials to derivatives of Gegenbauer polynomials, a relation the paper says is absent from standard special-function theory."],"supporting_citations":[{"why":"Fock's original momentum-space treatment that this paper modifies; supplies the stereographic projection and the 3-sphere equation for the Coulomb problem.","marker":"[9-11]"},{"why":"Standard textbook formulation of the hydrogen Schrödinger equation, its eigenfunctions, and the Stark-effect formula that the tensor method reproduces.","marker":"[1]"},{"why":"Source for the 4D spherical-function formalism, the Fock integral equation, and the Schwinger resolvent used as comparisons.","marker":"[3]"},{"why":"Introduces the harmonic-tensor formulas and the raising operator between multipole states that carry the paper's algebraic machinery.","marker":"[26]"},{"why":"Provides the symmetrization notation and multipole-tensor methods used throughout the invariant-tensor calculations.","marker":"[27]"},{"why":"Supplies the perturbation technique of obtaining the next approximation algebraically from the preceding one, used in the Stark-effect section.","marker":"[25]"},{"why":"Schwinger's Coulomb Green's function is the object rederived here as a Gegenbauer-polynomial series.","marker":"[50]"},{"why":"Reference for special-function properties of Gegenbauer and Laguerre polynomials, against which the new polynomial-correspondence identity is claimed to be unknown.","marker":"[32]"}],"fun_headline_variants":["Hydrogen wavefunctions as derivatives of 4D harmonics","Coulomb problem simplified via 4D harmonic tensors","Fock's theory recast in coordinate space via tensors","Hydrogen from 4D Laplace: harmonic polynomials","Coordinate-space twist on Fock's hydrogen theory"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the delta-function pieces produced by the polynomial part of the inverse Fourier transform can be thrown away; the paper verifies the resulting recipe on examples but supplies no general proof that this discarding is legitimate.","fun_headline_variants_meta":{"raw":{"variants":["Hydrogen wavefunctions as derivatives of 4D harmonics","Coulomb problem simplified via 4D harmonic tensors","Fock's theory recast in coordinate space via tensors","Hydrogen from 4D Laplace: harmonic polynomials","Coordinate-space twist on Fock's hydrogen theory"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000697,"raw_usage":{"total_tokens":3139,"prompt_tokens":924,"completion_tokens":2215,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":540,"completion_tokens_details":{"reasoning_tokens":2136}},"tokens_in":540,"tokens_out":2215,"duration_ms":15321,"temperature":1.0,"reasoning_tokens":2136,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T16:33:28.515539+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take a state beyond the paper's worked examples, such as $n=3$, $l=0$, and evaluate the right-hand side of Eq. (66) with and without the delta-function contributions generated by the polynomial numerator. If the two results differ from the standard hydrogen $3s$ wavefunction, the claimed algebraic map is incomplete; if the truncated expression reproduces it, the truncation rule is supported at least for that case.","supporting_citations":[{"cited_title":"D., Lifshitz E","cited_arxiv_id":null,"evidence_quote":"Standard textbook formulation of the hydrogen Schrödinger equation, its eigenfunctions, and the Stark-effect formula that the tensor method reproduces."},{"cited_title":"I., Zel’dovich","cited_arxiv_id":null,"evidence_quote":"Source for the 4D spherical-function formalism, the Fock integral equation, and the Schwinger resolvent used as comparisons."},{"cited_title":"P., Transition operator between multipole states and their tensor structure, Theor","cited_arxiv_id":null,"evidence_quote":"Introduces the harmonic-tensor formulas and the raising operator between multipole states that carry the paper's algebraic machinery."},{"cited_title":"Z., Multipoles and Ellipsoid Fields, Publishing House- MISiS, Moscow (2015), ISBN 978-5-600-01057-4","cited_arxiv_id":null,"evidence_quote":"Provides the symmetrization notation and multipole-tensor methods used throughout the invariant-tensor calculations."},{"cited_title":"P., Muratov R","cited_arxiv_id":null,"evidence_quote":"Supplies the perturbation technique of obtaining the next approximation algebraically from the preceding one, used in the Stark-effect section."},{"cited_title":"U., Vilenkin N","cited_arxiv_id":null,"evidence_quote":"Reference for special-function properties of Gegenbauer and Laguerre polynomials, against which the new polynomial-correspondence identity is claimed to be unknown."}],"review_version":1}