Pith. sign in

REVIEW 3 major objections 4 minor 17 references

Liouvillian solutions for second order linear differential equations with polynomial coefficients

T0 review · 3 major / 4 minor · reviewed 2026-08-14 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read For each fixed degree, the equations admitting Liouvillian solutions form a countable union of disjoint algebraic varieties of codimension at most n, and quasi-solvable Schrödinger potentials have at most d+1 Liouvillian energy levels.

desk verdict The spectral-set description and the d+1 eigenvalue bound are sound; the flagged gap in Theorem 3.3 is a terse proof, not a flaw. read the letter →

arxiv 1908.07666 v2 pith:4KNFFS6L submitted 2019-08-21 math-ph math.CAmath.MP

Primitivo B. Acosta-Humánez , David Blázquez-Sanz , Henock Venegas-Gómez This is my paper · ORCID
classification math-phmath.CAmath.MP MSC 34M1581Q35
keywords AnharmonicoscillatorsAsymptoticiterationmethodKovacicalgorithmLiouvilliansolutionsparameterspacequasi-solvablemodelsSchrödingerequationspectralvarieties
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

This paper studies the family of all second-order linear differential equations $u''+P(x)u'+Q(x)u=0$ with polynomial coefficients, parameterized by the coefficients of $P$ and $Q$. It characterizes, within that parameter space, the subset of equations that admit a Liouvillian solution—one built by exponentials, integrals, and algebraic operations. The main theorem states that for each fixed polynomial degree this 'spectral set' is a countable union of pairwise disjoint algebraic varieties $L_{2n,d}$, each of codimension at most $n$, and any compact part of the parameter space meets only finitely many of them. A reader should care because this turns a transcendental integrability question into finite algebraic equations and yields a concrete bound on the number of energy values for which a polynomial Schrödinger potential has Liouvillian eigenfunctions.

What carries the argument

The load-bearing object is the sequence of universal differential polynomials $\Delta_d\in\mathbb{Q}\{a,b\}$ generated by the asymptotic iteration method. Applied to the auxiliary equations (7) and (8) of Theorem 2.2, the recurrence (21) with initial values $\ell_0^\pm=\mp 2A$, $r_0^\pm=B\mp A'$ produces obstructions $\delta_d^\pm = r_d^\pm\ell_{d-1}^\pm-\ell_d^\pm r_{d-1}^\pm$, and $\Delta_d(A,B)$ is the specialization giving $\delta_d^+$ while $\Delta_d(-A,B)$ gives $\delta_d^-$. The vanishing of their product, together with the arithmetic condition $b_{n-1}^2=(n+2d)^2$, exactly detects the degree-$d$ polynomial factor $P_d$ in a Liouvillian solution $P_d e^{\pm\int A\,dx}$; the same polynomials work for every $n$, which is what makes the description degree-independent.

What would settle it

Take an explicit point on the variety $V^+_{4,2}$ from Table 2, say $b_1=6$ and $b_0^3+16a_0b_0-16=0$ with a generic $a_0$, and compute the AIM sequences $\ell_0^+$, $\ell_1^+$, $\ell_2^+$ and the obstruction $\delta_2^+$ from the recurrence (21) for the auxiliary equation (7). The theorem asserts $\delta_2^+=0$ exactly when a degree-2 polynomial solution exists; if any point on this variety has $\ell_2^+\ell_1^+\equiv 0$, or has $\delta_2^+=0$ without a polynomial solution, the algebraic criterion of Theorem 3.3 is false.

Watch

Extended reading notes

Core claim

On the paper's own terms, the central discovery is that Liouvillian integrability of (1) is exactly membership in one of the varieties $L_{2n,d}$ defined by the existence of a polynomial-hyperexponential solution of polynomial degree $d$, i.e. a solution $P_d(x)e^{\int A_k(x)\,dx}$ with $P_d$ of degree $d$ and $A_k$ a polynomial of degree $k$. Theorem 1.3 assembles these varieties into a countable disjoint union $L_{2n}=\bigcup_{d=0}^\infty L_{2n,d}$, with the codimension bound and the local finiteness property. The computational engine is Theorem 3.3: after writing the coefficient in the trace-free equation as $M(x)=A(x)^2+B(x)$, the equation $y''=M(x)y$ has a polynomial-hyperexponential solution of degree $d$ if and only if $b_{n-1}^2=(n+2d)^2$ and $\Delta_d(A,B)\,\Delta_d(-A,B)=0$, where $\Delta_d$ is a universal differential polynomial independent of $n$. This algebraic description supports Theorem 4.2: for a quasi-solvable polynomial Schrödinger potential with arithmetic condition $d$, at most $d+1$ values of the energy parameter admit a Liouvillian solution.

Load-bearing premise

The load-bearing premise is that, for every admissible coefficient pair $(A,B)$, the asymptotic-iteration entries $\ell_d^\pm$ and $\ell_{d-1}^\pm$ are nonzero polynomials; the paper asserts this 'comes easily' from a leading-degree comparison but gives no detailed argument, so if those entries can vanish identically for some admissible pair, the equivalence between $\Delta_d=0$ and existence of a degree-$d$ polynomial solution could fail.

Editorial extensions

If this is right

  • For fixed $n$, the spectral set $L_{2n}$ is a singular analytic submanifold of the parameter space $P_{2n}$, a countable disjoint union of algebraic components of codimension at most $n$.
  • Any bounded region of parameter space involves only finitely many possible solution degrees $d$, so a search for Liouvillian solutions can be truncated after checking finitely many $\Delta_d$ conditions.
  • For a quasi-solvable polynomial Schrödinger potential with integer arithmetic condition $d$, the number of energy values with a Liouvillian eigenfunction is at most $d+1$; explicit spectral systems for $x^6-(4J+1)x^2$ are computed through $d=13$.
  • The spectral varieties are nonempty for infinitely many degrees: $L_{2n,d}\neq\varnothing$ whenever $d\equiv 0$ or $1\pmod{n+1}$.
  • The varieties $L_{2n,d}$ and $L_{2n,k}$ are disjoint for $d\neq k$, so a given equation cannot admit polynomial-hyperexponential solutions of two different polynomial degrees.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The $d+1$ bound is likely sharp in the computed examples (e.g. $d=3$ gives two $\lambda$-values), so one could test numerically on other quasi-solvable potentials whether the bound is generically attained; a counterexample would point to an unstated assumption in the asymptotic-iteration criterion.
  • Because $\Delta_d$ is independent of $n$, the method suggests a purely algebraic elimination procedure for quasi-solvability: compute $\Delta_d(A(x),B(x)+\lambda)$ once, eliminate $x$, and obtain polynomial equations in the potential coefficients and $\lambda$.
  • The D'Alembert reduction preserves the polynomial degree of polynomial-hyperexponential solutions, so the same variety description should extend to equations with rational coefficients after clearing denominators; the paper does not develop this extension.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 4 minor

Summary. The paper studies the general second-order linear differential equation with polynomial coefficients, y'' + P(x)y' + Q(x)y = 0, under a non-degeneracy condition. It defines the spectral set L2n of equations admitting a Liouvillian solution and spectral subvarieties L2n,d according to polynomial-hyperexponential solutions of polynomial degree d. The main result, Theorem 1.3, states that L2n is the countable union of pairwise disjoint algebraic varieties L2n,d of codimension at most n in P2n, and that any compact subset of P2n meets only finitely many of them. The proof uses D'Alembert reduction to monic trace-free form, a dichotomy theorem from the authors' earlier work (Theorem 2.2) characterizing Liouvillian solutions via polynomial solutions of two auxiliary equations, and the Asymptotic Iteration Method to express solvability as the vanishing of universal differential polynomials Δ_d. A further result, Theorem 4.2, bounds by d+1 the number of energy parameters for which an algebraically quasi-solvable polynomial Schrödinger potential admits a Liouvillian solution, where d is an arithmetic condition read from the potential's coefficients.

Significance. If correct, the results provide a genuinely algebraic, degree-independent description of Liouvillian integrability for a large family of second-order equations, together with explicit equations for spectral varieties in low degree (Tables 2-3) and a quantitative bound on quasi-exactly-solvable eigenvalues of polynomial Schrödinger potentials (Theorem 4.2). The paper relies on a published dichotomy theorem [1] and on an external AIM criterion [15], so there is no circular dependence on the new results. The concrete computations, including the Turbiner potential example in Example 4.3, give explicit and falsifiable spectral predictions. However, two proof steps that are load-bearing for Theorem 1.3 are insufficiently detailed: the non-vanishing of the AIM sequences in Theorem 3.3 and the codimension argument in Proposition 3.6. The contribution is potentially valuable, but the central claims are not yet fully demonstrated.

major comments (3)
  1. [Section 3.2, proof of Theorem 3.3] The proof must justify that ℓ±_{d-1}ℓ±_d ≠ 0 for the auxiliary equations (7)-(8) before applying Theorem 3.2(ii), which requires exactly that condition to conclude that δ_p = 0 implies a polynomial solution. The sentence 'This comes easily from the fact that ℓ±_0 = ±2A(x) is of bigger degree than r0 = B ± A′' addresses only the initial step of a recurrence. A complete proof should include an induction showing, for instance, that deg ℓ±_j = (j+1)n and deg r±_j ≤ (j+1)n - 1 for every j, so that ℓ±_j is never the zero polynomial. Without this, the equivalence in Theorem 3.3 is not established, and Theorem 1.3(b)-(d) and Theorem 4.2 inherit the gap.
  2. [Section 3.3, Proposition 3.6 and Remark 3.7] The rank assertion rk(M±_{d-1,n}(A,B)) = d is stated without proof, and the passage from that rank to the conclusion codim(L′_{2n,d}, M2n) ≤ n is only sketched. The text does not identify precisely which n equations are locally independent, nor why the hyperplane condition b_{n-1} = ±(2d+n) together with the determinant conditions gives codimension at most n rather than n+1. Since Theorem 1.3(b) depends on this codimension bound, a rigorous local argument, including a transversality statement around a smooth point of L′_{2n,d}, is required.
  3. [Section 3.2, Example 3.4] The displayed identity Δ_d(x,b_0) = Δ_d(-x,b_0) = 2^{d+1} ∏_{k=0}^d (d-k) is incorrect as written: the factor k=d makes the product zero for every d, and the equality of Δ_d(x,b_0) with Δ_d(-x,b_0) is false in general (for d=0, Δ_0(x,b_0)=b_0-1 and Δ_0(-x,b_0)=b_0+1). Since this example is meant to illustrate the general method, the formula should be corrected; the final displayed form of L′_{2,d} is, however, correct.
minor comments (4)
  1. [Section 4, Theorem 4.2 proof] The phrase 'which is a polynomial in x of λ' is imprecise: Δ_d(A(x),B(x)+λ) is a polynomial in x and λ, and for a fixed λ the Liouvillian condition is that this polynomial vanish identically in x, i.e., all its x-coefficients vanish. The bound on the number of λ-values follows from the degree in λ of the gcd of these coefficient polynomials, which is at most d+1 by Lemma 4.1; this should be stated explicitly.
  2. [Section 3.1, Theorem 3.1 proof] There is a typo in the sentence 'v,β are solutions of u′ = αu, v′ = ℓ0 v y β′ = u2v'; it should read 'u′ = αu, v′ = ℓ0 v, and β′ = u²v'.
  3. [Section 3.2, Lemma 4.1] The bounds 'small or equal to (d+1)/2' and 'small or equal to (d+2)/2' should use integer floors or ceilings to avoid the impression that non-integer degrees are allowed.
  4. [Throughout] The manuscript contains many typographical errors and infelicities, for example 'inpendently' in the Introduction, 'Corolary 2.6', 'Liovullian' in the Section 2 heading, 'posibilities' in Proposition 2.4, and inconsistent punctuation in the display of Tables 2-3. A thorough language and copy-editing pass is needed.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the spectral variety description follows from independent Kovacic-based and AIM-based criteria; the only flagged weakness is a proof gap in the nonvanishing of AIM sequences, not a circular reduction.

full rationale

The paper's central claim (Theorem 1.3) is derived from two sources that are independent of the present result. Theorem 2.2, quoted from the authors' [1], is a parameter-free specialization of Kovacic's algorithm for y'' = M(x)y with the stated assumption that M is a monic polynomial of even degree; it does not assume that the spectral set L_{2n,d} has the algebraic description being proved, so citing it is real evidence rather than circular. The AIM criteria in Theorems 3.1 and 3.2 are explicitly translations of external results [15], and the universal differential polynomials Δ_d are computed from the AIM recurrence (21), not fitted to the spectral varieties. Theorem 3.3 is the point where the paper uses AIM to replace 'there exists a polynomial solution of degree d' by δ_d = 0; this is a mathematical equivalence supplied by Theorem 3.2, not a renaming of the conclusion. The only genuine weakness is that the proof of Theorem 3.3 asserts without a detailed argument that the AIM sequences ℓ^±_d and ℓ^±_{d−1} are nonzero for admissible A,B: 'We need only to check that ℓ±_{d−1}ℓ±_d ≠ 0 for the auxiliar equations. This comes easily from the fact that ℓ±_0 = ±2A(x) is of bigger degree than r0 = B ± A′.' That is a proof gap: if the nonvanishing failed, the 'if' direction of the equivalence could fail. However, a gap in justification is not circularity, because the required statement is an analytic property of the recurrence and is not imposed as an input equivalent to the target result. Likewise, the bound in Theorem 4.2 is obtained from Lemma 4.1, a degree estimate for Δ_d in b, not from fitting the number of eigenvalues to d. Consequently the derivation chain is self-contained modulo standard external theorems, and no predicted quantity reduces by construction to a fitted parameter or to the paper's own definitions.

Assumptions & free parameters 0 free parameters · 5 assumptions · 0 invented entities

No free parameters are fitted; the paper is pure mathematics. The load-bearing inputs are published theorems from differential Galois theory (Kovacic, Martinet-Ramis), the AIM criterion from Saad-Hall-Ciftci, and the authors' own Theorem 2.2 from [1], which is itself a parameter-free classification.

assumptions (5)
  • standard math Kovacic's algorithm provides a complete classification of Liouvillian solutions for second-order linear differential equations.
    Used implicitly to justify that any Liouvillian solution is polynomial-hyperexponential after a gauge transformation; cited as [11].
  • domain assumption Theorem 2.2 from [1]: for y''=M(x)y with M monic even degree 2n, the differential Galois group is either SL2(C) or C*⋉C; the latter iff ±b_{n-1}-n is a non-negative even integer 2d and an auxiliary equation has a polynomial solution.
    This is the load-bearing premise for the entire spectral set description. It is a prior result by two of the current authors, quoted as [1, Theorem 2.5].
  • standard math AIM criteria (Theorems 3.1 and 3.2 from [15]): the vanishing of the obstruction δ_p is necessary and sufficient (with non-vanishing ℓ_p ℓ_{p-1}) for a polynomial solution.
    Used in Theorem 3.3 to translate the auxiliary equation condition into algebraic equations on the coefficients.
  • standard math Martinet-Ramis theorem on Liouvillian solutions of the Whittaker equation.
    Used in Propositions 2.4 and 2.5 to characterize integrability of the canonical equation y''=(x^{2n}±(2d+n)x^{n-1})y.
  • domain assumption The D'Alembert transformation u=exp(-1/2∫P dx)y preserves the property of having a polynomial-hyperexponential solution and its polynomial degree.
    Used in Section 2.1 to reduce the parameter space from P2n to K2n.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Liouvillian solutions for second order linear differential equations with polynomial coefficients." pith.science (2026). https://pith.science/paper/4KNFFS6L

@misc{pith2026190807666,
  author       = {Pith},
  title        = {Pith review of: Liouvillian solutions for second order linear differential equations with polynomial coefficients},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/4KNFFS6L}},
  note         = {Machine review of arXiv:1908.07666}
}
read the original abstract

In this paper we present an algebraic study concerning the general second order linear differential equation with polynomial coefficients. By means of Kovacic's algorithm and asymptotic iteration method we find a degree independent algebraic description of the spectral set: the subset, in the parameter space, of Liouiville integrable differential equations. For each fixed degree, we prove that the spectral set is a countable union of non accumulating algebraic varieties. This algebraic description of the spectral set allow us to bound the number of eigenvalues for algebraically quasi-solvable potentials in the Schr\"odinger equation.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

17 extracted references · 17 canonical work pages

  1. [1]

    Acosta-Hum´ anez and D

    P. Acosta-Hum´ anez and D. Bl´ azquez-Sanz. Non-integrability of some hamiltonian sys- tems with rational potential. Discrete and Continuous Dynamical Systems Series B , 10(2–3):265–293, 2008

    work page 2008

  2. [15]

    N. Saad, R. Hall, and H. Ciftci. Criterion for polynomia l solutions to a class of linear differential equations of second order. Journal of Physics A: Mathematical and General , 39(43):13445–13454, 2006

    work page 2006

  3. [2]

    Acosta-Hum´ anez.Galoisian Approach to Supersymmetric Quantum Mechanics

    P.B. Acosta-Hum´ anez.Galoisian Approach to Supersymmetric Quantum Mechanics. PhD thesis, Universitat Polit` ecnica de Catalunya, https://w ww.tdx.cat/handle/10803/22723, 2009

    work page 2009

  4. [3]

    Acosta-Hum´ anez.Galoisian Approach to Supersymmetric Quantum Mechanics

    P.B. Acosta-Hum´ anez.Galoisian Approach to Supersymmetric Quantum Mechanics. Th e integrability analysis of the Schr¨ odinger equation by mea ns of differential Galois theory . VDM Verlag, Dr M¨ uller, Saarbr¨ ucken, Deutschland, 2010. 16

    work page 2010

  5. [4]

    Acosta-Hum´ anez, J.J

    P.B. Acosta-Hum´ anez, J.J. Morales-Ruiz, and J.-A. Weil. Galoisian approach to integra- bility of schr¨ odinger equation.Reports on Mathematical Physics , 67(3):305–374, 2011

    work page 2011

  6. [5]

    Bender and G

    C. Bender and G. Dunne. Quasi-exactly solvable systems a nd orthogonal polynomials. Journal of Mathematical Physics , 37(1):6–11, 1996

    work page 1996

  7. [6]

    Galoisian approach for a sturm-liouville problem on the infinite interval

    David Bl´ azquez-Sanz and Kazuyuki Yagasaki. Galoisian approach for a sturm-liouville problem on the infinite interval. Methods and Applications of Analysis , 19(3):267–288, 2012

    work page 2012

  8. [7]

    Ciftci, R

    H. Ciftci, R. Hall, and N. Saad. Asymptotic iteration met hod for eigenvalue problems. Journal of Physics A: Mathematical and General , 36(47):11807–11816, 2003

    work page 2003

Show all 17 references
  1. [8]

    Ciftci, R

    H. Ciftci, R. Hall, N. Saad, and E. Dogu. Physical applica tions of second-order linear dif- ferential equations that admit polynomial solutions. Journal of Physics A: Mathematical and Theoretical, 43(41):415206–415219, 2010

    2010

  2. [9]

    T. Combot. Integrability of the one dimensional schr¨ od inger equation. Journal of Math- ematical Physics , 59(2):022105, 2018

    2018

  3. [10]

    Duval and M

    A. Duval and M. Loday-Richaud. Kovacic’s algorithm and its application to some families of special functions. Applicable Algebra in Engineering, Communication and Comput ing, 3(3):211–246, 1992

    1992

  4. [11]

    J. Kovacic. An algorithm for solving second order linea r homogeneous differential equa- tions. Journal of Symbolic Computation , 2(1):3–43, 1986

    1986

  5. [12]

    Martinet and J.-P

    J. Martinet and J.-P. Ramis. Theorie de galois differenti elle et resommation. In E. Tournier, editor, Computer Algebra and Differential Equations , pages 117–214. Aca- demic Press, London, 1989

    1989

  6. [13]

    R. Hall N. Saad and H. Ciftci. Sextic harmonic oscillato rs and orthogonal polynomials. Journal of Physics A: Mathematical and General , 39(26):8477–8486, 2006

    2006

  7. [14]

    G. A. Natanzon. Investigation of a one dimensional schr ¨ odinger equation that is generated by a hypergeometric equation. Vestnik Leningrad. Univ , 10:22–28, 1971. in Russian

    1971

  8. [16]

    A. V. Turbiner. Quantum mechanics: problems intermedi ate between exactly solvable and completely unsolvable. Soviet Phys. JETP , 10(2):230–236, 1988

    1988

  9. [17]

    Venegas-G´ omez

    H. Venegas-G´ omez. Enfoque galoisiano de la ecuaci´ on de schr¨ odinger con potenciales polinomiales y polinomios de laurent. Master’s thesis, Uni versidad Nacional de Colombia, sede Medell ´ ın, http://bdigital.unal.edu.co/71580/, 2018. 17

    2018

Pith tools

Reviewed August 14, 2026 · model on record in the stance chip above.