{"id":"d3a57e64-1537-4a93-9e5d-d62fd9566dcb","arxiv_id":"2608.12909","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A triangular-lattice finite-difference model is constructed that is maximally superintegrable, exactly solvable by Krawtchouk and dual Hahn polynomials, and whose continuum limit reproduces the continuous Smorodinsky-Winternitz II system with Hermite-Laguerre eigenfunctions.","lead":"The authors construct a finite lattice version of the Smorodinsky-Winternitz II oscillator with exact energy levels and eigenfunctions, then show it tends to the continuous system in a scaling limit. The construction connects Hahn-type symmetry algebras on the lattice to the Laguerre-Heun algebra of the continuous model.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Krawtchouk–Hermite limit (B.14) is inconsistent with (B.13): at n=1 the LHS is -x√(2(1-p)/p), not the printed RHS, and with (3.23) the Krawtchouk factor vanishes under p=N^{-θ}; (4.6) does not follow.","rationale":"The strongest claim includes both maximal superintegrability and recovery of the continuous SW II system. The finite-difference part is explicit and plausible, but the continuum eigenfunction limit is load-bearing and rests on (B.14), which is demonstrably inconsistent at n=1. This is not a question of outside consensus; it is an internal mismatch between displayed equations. The reader flagged the continuum limit as a condition; I sharpen this to a concrete normalization error in the Krawtchouk factor. The rest of the paper, including the cubic algebra and Tratnik solution, may be salvageable, so I keep the verdict conditional rather than reject outright.","tokens_in":18065,"tokens_out":45622,"duration_ms":378604,"concrete_test":"Specialize to n2=0, n1=1. First, evaluate the N→∞ limit in (B.14) using \\hat K_1 = 1 - x/(Np) and compare with the printed RHS; second, compute the left side of (3.20) for P_{1,0} from (3.23): if it is not 1, the normalization is wrong. A corrected derivation should replace the exponent n1 in (3.23) by n1/2 and re-run (4.6); the Hermite factor will emerge only after that change.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"Equations (B.13), (3.23) and (B.14) do not support the eigenfunction limit (4.6), already in the n2=0 sector. For \\hat K_n defined in (B.13), the n=1 case of (B.14) reads LHS = sqrt(N)(1 - (pN + sqrt(2p(1-p)N)x)/(Np)) = -x sqrt(2(1-p)/p), while the printed RHS is -2x sqrt(2) p/(1-p); these are not equal for fixed p. A correct Krawtchouk–Hermite limit carries a factor ((1-p)/p)^{n/2} or its equivalent, not (p/(1-p))^n. Inserting the printed (B.14) into (3.23), whose K-factor already multiplies sqrt(...) by (p/(1-p))^{n1}, gives a limiting factor (p/(1-p))^{2n1}; under p=N^{-θ} this tends to zero for n1>0. Even with a corrected (B.14), the printed exponent in (3.23) leaves a residual (p/(1-p))^{n1/2} or similar, again vanishing. Equivalently, the P_{1,0} eigenfunction from (3.23) has squared norm p/(1-p) with respect to the stated weight, contradicting the orthonormality (3.20). The finite-difference Hamiltonian and ladder structure may be correct, but the claimed convergence to Hermite–Laguerre wavefunctions is not established by the manuscript as written.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript constructs a finite-difference analogue of the two-dimensional Smorodinsky–Winternitz II superintegrable system on a triangular region of the square lattice. The Hamiltonian is written as the sum of two commuting number operators with a local weight function; the paper claims a pair of additional integrals of motion, a cubic symmetry algebra admitting a Hahn-algebra presentation, and an exact spectral solution in terms of bivariate Krawtchouk and dual Hahn polynomials of Tratnik type. A continuum limit is then claimed to recover the continuous Smorodinsky–Winternitz II Hamiltonian with Hermite–Laguerre eigenfunctions, while the Hahn presentation degenerates and the limiting symmetry algebra is identified with the Laguerre–Heun algebra. The exact-discrete part of the construction is explicit and largely self-contained, but the continuum-limit claims rest on a specific asymptotic relation for Krawtchouk polynomials that is incorrect as stated.","tokens_in":18361,"tokens_out":17149,"duration_ms":167191,"significance":"If the continuum-limit claims can be made correct, the paper would be a valuable contribution: it gives a fully explicit finite model with exact ladder operators and closed-form eigenfunctions, and it proposes a discrete origin for the Laguerre–Heun algebra of the continuous Smorodinsky–Winternitz II system. The exact discrete solution, the identification with Tratnik polynomials, and the algebraic presentation of the symmetry are attractive and likely of independent interest. However, the paper's advertised central claim—that the continuum limit recovers the continuous system with Hermite–Laguerre eigenfunctions—is not established by the present derivations, because the Krawtchouk–Hermite limit relation used in the proof is wrong and the normalization of the eigenfunctions is internally inconsistent. The paper therefore needs substantive revision before the main claims can be accepted.","major_comments":[{"comment":"Equation (B.14) is incorrect. For n=1, \\hat K_1(x;p,N)=1-x/(Np), so with z=pN+\\sqrt{2p(1-p)N}\\,x the left-hand side of (B.14) equals \\sqrt{N}\\hat K_1(z;p,N)=-x\\sqrt{2(1-p)/p}. The printed right-hand side is (-1)^1 H_1(x)\\sqrt{2}\\,(p/(1-p)) = -2\\sqrt{2}\\,x\\,p/(1-p). These two expressions are not equal for fixed p. The correct fixed-p Hermite limit carries a factor ((1-p)/p)^{n/2}/\\sqrt{2^n n!}, not (p/(1-p))^n. Since (B.14) is the only ingredient supplied for the first limit in (4.6), the claimed recovery of the Hermite factor in the eigenfunctions (4.15) is not established.","section":"Appendix B.4, Eq. (B.14)"},{"comment":"The normalization in (3.23) is inconsistent with the orthonormality stated in (3.20). For n_1=1, n_2=0, equations (3.22)-(3.24) give P_{1,0}(x_1,x_2)=-\\sqrt{N}\\,(p/(1-p))\\,\\hat K_1(x_1;p,N)/\\sqrt{\\Gamma(\\alpha+1)}. Using the binomial variance for \\hat K_1 and \\sum_{x_2} w_2=\\Gamma(\\alpha+1), the weighted squared norm in (3.20) equals p/(1-p), not 1. Moreover, even if (B.14) is corrected, the factor (p/(1-p))^{n_1} in (3.23) multiplies the Hermite limit by an additional factor (p/(1-p))^{n_1/2}, which tends to zero for n_1>0 under p=N^{-\\theta}; hence (4.6) cannot follow from the printed definitions. The prefactor in (3.23) appears to be missing a compensating factor ((1-p)/p)^{n_1/2}, and the normalization of the Krawtchouk factor must be reconsidered.","section":"Section 3.3, Eq. (3.23)"},{"comment":"The limits in (4.5) are asserted after applying the Taylor expansion (4.4) to the number operators (2.9) and the ladder operators (3.3)-(3.4), but no expansion or cancellation is displayed. The coefficient functions in (2.9) and Appendix A contain several powers of N in numerators and denominators, and the expansion (4.4) produces contributions of different orders in N. The claimed convergence to the differential operators in (4.5) therefore requires an explicit leading-order computation showing that all singular contributions cancel. Since (4.5) is the basis for the continuum Hamiltonian in (4.14), this verification should be included.","section":"Section 4.1, Eq. (4.5)"}],"minor_comments":[{"comment":"Reference [4] is listed as \"arXiv:XXXX.XXXXX, 2026\"; a placeholder arXiv identifier is not acceptable in a submitted manuscript and should be replaced before publication.","section":"References, [4]"},{"comment":"The Jacobian determinant in (4.8) is negative with the stated orientation of the change of variables; the absolute value should be taken before multiplying by the weight function.","section":"Section 4.1, Eq. (4.8)"},{"comment":"The text says \"what follows is self-contained,\" but the definition of superintegrability used throughout is deferred to the companion paper [4]; either include a precise definition or soften the self-containment claim.","section":"Section 1"},{"comment":"The phrase \"respectively as follow\" should read \"respectively as follows\" in the sentences before (B.8) and (B.12).","section":"Appendix B.2 and B.3"}],"recommendation":"major_revision","confidential_remarks":"The continuum-limit claim is central to the abstract and conclusions, and the identified errors in (B.14) and (3.23) are load-bearing. They appear correctable within the manuscript's framework, so I recommend major revision rather than rejection. The editor may also wish to verify that the companion papers [4] and [5] are available in final form, since the manuscript relies on them for framing and for the identification of the limiting Laguerre–Heun algebra."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's the short version. The paper builds a finite-difference model of the Smorodinsky–Winternitz II system on a triangular lattice with two commuting number operators, a ladder structure, and eigenfunctions identified as Tratnik-type Krawtchouk/dual Hahn polynomials. That part is substantial and, as far as I can tell from the manuscript, the operators are explicit enough to verify. The Hahn-algebra presentation and the singular degeneration to the Laguerre–Heun algebra are also interesting and novel. So the paper is not a mirage.\n\nThe soft spot is not small, though. The stress-test note is right. Equation (B.14) is inconsistent with the definition (B.13) even at n=1: with x = pN + sqrt(2p(1-p)N) x, the LHS of (B.14) tends to -x sqrt(2(1-p)/p), while the printed RHS is -2 sqrt(2) x p/(1-p). The correct Krawtchouk–Hermite limit carries a factor ((1-p)/p)^{n/2} (or equivalent), and (B.14) as printed cannot be fixed by a mere sign. Worse, the normalization in (3.23) is off: for (n1,n2)=(1,0), the squared norm of P with respect to the weight (2.12) is p/(1-p), not 1, contradicting the stated orthonormality (3.20). This is not a cosmetic detail. The claimed limit (4.6) uses (B.14) and (3.23); under p = N^{-\\theta} the prefactor leaves a vanishing factor (p/(1-p))^{n1} (or (p/(1-p))^{2n1} if you use their printed B.14) for every n1>0, so the Hermite–Laguerre wavefunctions are not recovered. In other words, the continuum limit of the eigenfunctions—one of the two headline results—does not follow from the manuscript as written. The Hamiltonian and ladder limits in (4.5) may well be fine; the problem is in the matching to the continuous wavefunctions.\n\nI don't see fitted-parameter circularity. The model is constructed exactly, and the exact solution is imported from the independent Tratnik/Leonard-pair literature. That part is honest. The symmetry algebra and Casimir are asserted as 'direct computation' with no intermediate steps; that is a referee annoyance rather than a demonstrated error, but it compounds the difficulty of checking the paper.\n\nThe companion reference [4] is a placeholder with no arXiv identifier, which is a practical obstacle for a paper that leans on it.\n\nRecommendation: this deserves a serious referee, not a desk reject, because the construction is real and the error looks localized to a normalization/limit formula. But as submitted the central claim about Hermite–Laguerre eigenfunctions is not established, and the orthonormality contradiction is load-bearing. I'd send it out, expecting major revision. I wouldn't cite it in its current form.","headline":"The finite SW II model is a real construction, but the continuum limit as written is broken by a Krawtchouk–Hermite normalization error, so the central eigenfunction claim does not hold.","tokens_in":18942,"tokens_out":14006,"would_cite":false,"duration_ms":110532,"reading_group":"maybe","serious_thinker":"no","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["81R12","33C45","39A70"],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper constructs a finite-difference model on a triangular lattice that is maximally superintegrable with Krawtchouk-dual-Hahn eigenfunctions and whose continuum limit is the Smorodinsky-Winternitz II system.","keywords":["superintegrable systems","Smorodinsky-Winternitz II","finite-difference Hamiltonian","Krawtchouk polynomials","dual Hahn polynomials","Hahn algebra","Laguerre-Heun algebra","continuum limit"],"falsifier":"Compute the limits directly: substitute the explicit ladder coefficients of Appendix A into the difference operators, apply the scaling (4.1) with $p=N^{-\\theta}$, and keep the Taylor expansions through all orders that contribute to $O(1/p)$ and $O(1)$; if any non-cancelling term remains beyond the operators listed in (4.5), or if the gauge-transformed operators (4.13) differ from the stated Hermite and Laguerre forms, the continuum-limit claim is refuted.","tokens_in":17821,"feed_emoji":"🔺","tokens_out":11099,"duration_ms":101139,"temperature":0.7,"pith_summary":"The paper constructs a finite quantum system defined on a triangular region of the two-dimensional square lattice and claims it is an exact discrete analogue of the Smorodinsky-Winternitz II superintegrable system. The Hamiltonian is the sum of two commuting finite-difference number operators, so the model has three independent integrals of motion and is maximally superintegrable. Its eigenfunctions are given in closed form as products of Krawtchouk and dual Hahn polynomials of Tratnik type. Under a rescaling of the lattice variables with $p=N^{-\\theta}$, the Hamiltonian, ladder operators, and eigenfunctions converge to the continuous Smorodinsky-Winternitz II Hamiltonian and its Hermite-Laguerre eigenfunctions. The paper further shows that the discrete symmetry algebra's Hahn presentation becomes singular in this limit, while the limiting symmetry operators close on the Laguerre-Heun algebra.","feed_headline":"A triangular lattice hosts a discrete Smorodinsky-Winternitz II system","feed_subtitle":"Three commuting integrals, exact Krawtchouk-dual-Hahn eigenfunctions, and a continuum limit to Hermite-Laguerre.","key_machinery":"The load-bearing machinery is the pair of commuting finite-difference number operators $N_1,N_2$ of (2.9), whose commutativity and self-adjointness with respect to the factorized weight (2.12) force a finite spectrum on $\\mathcal{R}(N)$. Around them are built the ladder operators $(a_1,a_2,a_1^\\dagger,a_2^\\dagger)$ with structure functions (3.5), which satisfy the dynamical algebra (3.6) and give the raising/lowering actions (3.18) on energy eigenstates; the symmetry operators $C_1=N_2+\\alpha+1/2$ and $C_2=\\{a_1,a_2^\\dagger\\}+\\{a_1^\\dagger,a_2\\}$ are formed from these and close on the cubic algebra (3.9), which is brought to the Hahn algebra (3.15) by the renormalization (3.13). The continuum limit is carried by the Taylor expansion (4.4) of finite differences at shifted lattice points under the scaling $x_1=pN+\\sqrt{2p(1-p)N}\\,x$, $x_2=N-x_1-y/\\sqrt{p}$, $p=N^{-\\theta}$, followed by the gauge factor $g=\\pi^{-1/4}e^{-x^2/2}y^{\\alpha+1/2}e^{-y^2/4}$, which converts the limiting difference operators into Hermite and Laguerre differential operators.","core_discovery":"The central discovery is that the operator $H=N_1+N_2+\\alpha/2+1$ with $N_1,N_2$ the finite-difference number operators of (2.9) acting on the triangular lattice $\\mathcal{R}(N)$ and self-adjoint with respect to the weight (2.12) is a maximally superintegrable two-dimensional system whose symmetry algebra admits a Hahn-algebra presentation. The additional integrals are $C_1=N_2+\\alpha+1/2$ and $C_2=\\{a_1,a_2^\\dagger\\}+\\{a_1^\\dagger,a_2\\}$, built from the ladder operators; together with $H$, they close on the cubic algebra (3.9), which the redefinition (3.13) brings to Hahn form. The energy eigenfunctions are exactly the bivariate Krawtchouk and dual Hahn polynomials of Tratnik type. Under the scaling (4.1) with $p=N^{-\\theta}$ and the gauge transformation (4.12), the finite-difference operators converge to the differential operators (4.5) of the continuous Smorodinsky-Winternitz II system, and the polynomial eigenfunctions converge to Hermite-Laguerre products, so the finite model is a genuine realization rather than a finite-difference approximation. The Hahn presentation does not survive this limit; the limiting symmetry algebra is the Laguerre-Heun algebra.","pith_inferences":["Because the finite Hilbert space has dimension $(N+1)(N+2)/2$ and the ladder operators act by explicit square-root factors, the model offers a ready-made exactly solvable finite register; using it as a quantum-information platform is an extension not pursued in the paper.","The pattern seen here, regular operator limits coexisting with a singular algebraic presentation, could be a general feature of finite superintegrable models; a natural test is whether other bivariate orthogonal-polynomial systems attached to the Askey scheme exhibit the same presentation degeneration.","One could attempt to invert the construction: start from the Laguerre-Heun algebra of the continuous system and look for other finite-difference realizations whose limit is the same algebra, which would sharpen the sense in which this lattice model is canonical."],"forward_implications":["The model is maximally superintegrable: the three independent integrals $H,C_1,C_2$ on a two-dimensional configuration space imply maximal degeneracy of the finite spectrum.","The spectral problem is exactly solvable: eigenfunctions are explicit products of Krawtchouk and dual Hahn polynomials with orthogonality governed by the factorized weight (2.12).","The continuum limit recovers the continuous Smorodinsky-Winternitz II Hamiltonian, its Hermite-Laguerre eigenfunctions, and its separation of variables, so the finite model is a genuine realization rather than a finite-difference approximation of the equations of motion.","The Hahn presentation of the discrete symmetry algebra is not stable under the continuum limit: it becomes singular, while the limiting symmetry operators close on the Laguerre-Heun algebra.","The construction supplies a finite-dimensional origin for the Laguerre-Heun algebra, because the regular limit of the discrete symmetry algebra is exactly that continuous algebra."],"supporting_citations":[{"why":"Provides the factorized $A_2$-Leonard pair whose difference equations and orthogonal polynomials are exactly the spectral problem solved here.","marker":"[6]"},{"why":"Supplies the Tratnik-type bivariate Krawtchouk and dual Hahn polynomials that form the model's eigenfunctions.","marker":"[28]"},{"why":"Identifies the continuous Smorodinsky-Winternitz II potential that the continuum limit is designed to recover.","marker":"[11]"},{"why":"Supplies the review definition of the continuous system, the superintegrability framework, and the separation-of-variables setup used in Section 4.","marker":"[24]"},{"why":"Identifies the Laguerre-Heun algebra of the continuous system, which the discrete symmetry algebra converges to in the limit.","marker":"[5]"},{"why":"Gives the Hahn algebra presentation in which the discrete symmetry algebra (3.15) is cast.","marker":"[10]"},{"why":"Is the source of the hypergeometric, Hermite, Laguerre, Krawtchouk, and dual Hahn formulas used in the continuum limit and Appendix B.","marker":"[21]"},{"why":"States the exact-solvability framework and the Hermite-Laguerre eigenfunctions of the continuous system that are matched in the limit.","marker":"[27]"}],"fun_headline_variants":["Maximally superintegrable discrete SWII on triangular lattice","Exact eigenfunctions and Hahn algebra for a discrete SWII system","Triangular lattice discrete SWII: Hahn symmetry and Hermite-Laguerre limit","Discrete SWII with continuum limit to Hermite-Laguerre","Krawtchouk-dual Hahn eigenfunctions for discrete SWII"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the Taylor expansion (4.4) of the finite differences at shifted lattice points, together with the scaling $p=N^{-\\theta}$, produces exact cancellation of all singular terms of order $1/p$ and order $1$, so that $N_1,N_2$ and the ladder operators converge to the differential limits (4.5); if those cancellations fail, the claimed recovery of the continuous Smorodinsky-Winternitz II Hamiltonian and its Hermite-Laguerre eigenfunctions fails.","fun_headline_variants_meta":{"raw":{"variants":["Maximally superintegrable discrete SWII on triangular lattice","Exact eigenfunctions and Hahn algebra for a discrete SWII system","Triangular lattice discrete SWII: Hahn symmetry and Hermite-Laguerre limit","Discrete SWII with continuum limit to Hermite-Laguerre","Krawtchouk-dual Hahn eigenfunctions for discrete SWII"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001467,"raw_usage":{"total_tokens":5936,"prompt_tokens":1018,"completion_tokens":4918,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":634,"completion_tokens_details":{"reasoning_tokens":4831}},"tokens_in":634,"tokens_out":4918,"duration_ms":31657,"temperature":1.0,"reasoning_tokens":4831,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T20:49:22.802483+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the limits directly: substitute the explicit ladder coefficients of Appendix A into the difference operators, apply the scaling (4.1) with $p=N^{-\\theta}$, and keep the Taylor expansions through all orders that contribute to $O(1/p)$ and $O(1)$; if any non-cancelling term remains beyond the operators listed in (4.5), or if the gauge-transformed operators (4.13) differ from the stated Hermite and Laguerre forms, the continuum-limit claim is refuted.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the factorized $A_2$-Leonard pair whose difference equations and orthogonal polynomials are exactly the spectral problem solved here."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the Tratnik-type bivariate Krawtchouk and dual Hahn polynomials that form the model's eigenfunctions."},{"cited_title":"Friˇ s, V","cited_arxiv_id":null,"evidence_quote":"Identifies the continuous Smorodinsky-Winternitz II potential that the continuum limit is designed to recover."},{"cited_title":"42, 423001","cited_arxiv_id":null,"evidence_quote":"Supplies the review definition of the continuous system, the superintegrability framework, and the separation-of-variables setup used in Section 4."},{"cited_title":"The 2D Smorodinsky--Winternitz II system and the Laguerre--Heun algebra","cited_arxiv_id":"2606.00903","evidence_quote":"Identifies the Laguerre-Heun algebra of the continuous system, which the discrete symmetry algebra converges to in the limit."},{"cited_title":"14, 1531–1535","cited_arxiv_id":null,"evidence_quote":"Gives the Hahn algebra presentation in which the discrete symmetry algebra (3.15) is cast."},{"cited_title":"Lesky, and Ren´ e F","cited_arxiv_id":null,"evidence_quote":"Is the source of the hypergeometric, Hermite, Laguerre, Krawtchouk, and dual Hahn formulas used in the continuum limit and Appendix B."},{"cited_title":"Turbiner, and Pavel Winternitz,Exact solvability of superintegrable systems, Journal of Mathematical Physics42(2001), no","cited_arxiv_id":null,"evidence_quote":"States the exact-solvability framework and the Hermite-Laguerre eigenfunctions of the continuous system that are matched in the limit."}],"review_version":1}