{"id":"6d96f5ef-1ad1-4695-99c4-6bde197cba75","arxiv_id":"2501.04396","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Linear systems of moment differential equations with analytic coefficients have unique local analytic solutions for strongly regular moment sequences, with the solution radius governed by a factorial-growth dichotomy.","lead":"This paper proves local existence and analyticity for linear systems of moment differential equations with analytic coefficients, where the moment derivative generalizes ordinary, q-difference, and fractional derivatives. It also characterizes when such systems can be converted into higher-order moment equations and why they cannot be translated away from the origin.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Lemma 2's key comparison c_p ≤ c̃_p is false for α<1: c_p uses Γ(1+αp)/Γ(1+α(p+1)) while c̃_p uses 1/(p+1), and for the basic case m_p=Γ(1+p/2), A(z)=1/(1-z), the intended c_1 exceeds c̃_1. This breaks the printed proof of Theorem 1 in the Assumption (B) case.","rationale":"The reader's weakest-assumption analysis identifies exactly the same load-bearing defect: the comparison between c_p and c̃_p in Lemma 2 breaks for α<1 because of mismatched recursion factors. My reading of the full proof confirms this is not a cosmetic typo: the formal solution of (11) cannot be given by (8), the Caputo reformulation is therefore inapplicable to ĉ, and the final chain ‖y_p‖≤...≤c̃_p has no valid middle step. The concrete numerical check with the simplest fractional moment sequence and A(z)=1/(1-z) shows c_1>c̃_1, making the failure immediate. I also note a related but secondary defect: the recursion defining c_p does not specify c_1, since it is stated for p≥1 and produces c_{p+1}; this reinforces that the printed argument is incomplete. The main theorem is not thereby disproved—there is a natural repair by using the true Caputo majorant recursion (11), and the intended bound likely follows from the integral equation (13). Because the flaw is in the proof of a key lemma rather than a demonstrated counterexample to the theorem, the appropriate disposition remains the reader's CONDITIONAL verdict. The later issues (e.g., Lemma 5 missing an invertibility hypothesis) are less load-bearing for the central existence claim, so they do not alter this assessment.","tokens_in":19344,"tokens_out":22160,"duration_ms":189858,"concrete_test":"Recompute the Lemma 2 comparison in the concrete case m_p=Γ(1+p/2), A(z)=1/(1-z), with K=2, c=C=1, c0=1: verify that c_1=1/Γ(3/2)≈1.128 > c̃_1=1, so the printed inequality c_p ≤ c̃_p fails at p=1. Then analyze the repaired Caputo majorant with coefficients d_{p+1}=Γ(1+αp)/Γ(1+α(p+1)) c C̃ ∑_{k=0}^p K^{p-k} d_k: (i) prove by induction that d_p ≥ c_p for every p≥0; and (ii) show the series ∑ d_p z^p has a positive radius of convergence, for example by comparing d_p with the explicit solution of the integral equation (13) and its known convergence on [0,r1). If both (i) and (ii) hold, Lemma 2 is repairable and Theorem 1 stands; if d_p does not dominate c_p or the repaired series has zero radius, the Assumption (B) existence proof collapses.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central existence theorem (Theorem 1) rests on Lemma 2 for strongly regular sequences whose moment growth is a fractional power of factorial, i.e. Assumption (B). In the proof of Lemma 2, the sequence c_p is defined with the factor Γ(1+αp)/Γ(1+α(p+1)), while c̃_p is defined by recursion (8), which uses the classical factor 1/(p+1). The paper then states that the series ĉ(z)=∑ c̃_p z^p solves the moment-differential problem (11); this is false, since the coefficient recursion for (11) requires the Γ-ratio, not (p+1)^{-1}. Consequently the passage through the Caputo problem (12)-(13) does not apply to ĉ, and the final estimate ‖y_p‖ ≤ c_p ≤ c̃_p is unsupported. The inequality c_p ≤ c̃_p is in fact false: for m_p=Γ(1+p/2) (α=1/2, C=1), A(z)=1/(1-z), take c=1, K=2>1, c0=1. The intended definition of c_1 gives c_1=Γ(1)/Γ(3/2)=1.128..., while (8) with C̃=1 gives c̃_1=1. Thus c_1>c̃_1. The recursion as printed also leaves c_1 undefined, since c_{p+1} is only specified for p≥1. Because the radius bound for y depends on c̃_p, the proof of Lemma 2 does not establish the existence of a positive radius in the α<1 regime. This is a genuine gap in the written mathematics, though not a demonstrated falsity of Theorem 1: a repaired majorant should be the formal solution of (11), whose coefficients satisfy d_{p+1}=Γ(1+αp)/Γ(1+α(p+1)) c C̃ ∑_{k=0}^p K^{p-k} d_k, and one would need to show d_p ≥ c_p and that ∑ d_p z^p has positive radius.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper studies first-order linear systems of moment differential equations with variable analytic coefficients, ∂_m y = A(z)y + b(z), on a disc centered at the origin. The main result, Theorem 1, asserts that for a strongly regular sequence M admitting a pair of kernel functions for M-summability, and for its associated moment sequence m, the homogeneous Cauchy problem has a unique solution with components holomorphic in a neighborhood of the origin. The proof splits into Assumption (A), where the moment ratios grow at least linearly in p and the full disc D(0,r) is preserved, and Assumption (B), where the ratios behave like Γ(1+αp)/Γ(1+α(p−1)) with 0<α<1 and only a possibly smaller disc is obtained. The paper also relates moment differentiation to Caputo fractional derivatives, shows that translation equivariance forces m=(C p!), converts systems to higher-order equations under a cyclic-vector condition, and extends constant-coefficient results from [15].","tokens_in":19811,"tokens_out":12677,"duration_ms":119476,"significance":"If correct, the paper gives a clean structural statement: local analytic solvability for variable-coefficient moment systems holds for all strongly regular moment sequences, and the radius is controlled by whether the moment growth is at least factorial or only fractional-factorial. The comparison with Caputo fractional equations is explicit and the two regimes are illustrated by Example 1. The proofs are based on explicit majorant sequences rather than fitted parameters, and the paper draws on published external theory (strongly regular sequences, generalized summability, [15]) rather than on assumptions invented for the paper. The main obstacle is the majorant comparison in Lemma 2: the written proof does not justify the estimate that produces the convergence radius in the α<1 case. The theorem may still be true, but the proof needs repair.","major_comments":[{"comment":"The comparison c_p ≤ c̃_p is not established and is in fact false as printed. The sequence c_p is defined with the factor Γ(1+αp)/Γ(1+α(p+1)), while c̃_p is defined by (8) with the classical factor 1/(p+1); consequently ĉ(z)=∑ c̃_p z^p is not a formal solution of the moment-differential problem (11), whose recursion requires the Γ-ratio. For m_p=Γ(1+p/2) (α=1/2, C=1), A(z)=1/(1−z), K=2, c=c0=1, the intended definition gives c_1=Γ(1)/Γ(3/2)=1.128..., whereas (8) with C̃=1 gives c̃_1=1, so c_1>c̃_1. Since the final display of Lemma 2 uses c_p≤c̃_p to bound ‖y_p‖ and obtain a positive radius, the printed proof does not establish the conclusion in the Assumption (B) case. A repaired proof should take the majorant to be the formal solution of (11) and show directly that its coefficients dominate c_p and have positive radius of convergence.","section":"Section 3, Lemma 2 (p. 9)"},{"comment":"The recursion c_{p+1}=... is written 'for all p≥1', so c_1 is never defined, although the induction claim for p=1 uses c_1. The intended inequality Y_1 Γ(1)/Γ(1+α)≤c_1 requires a separate definition of c_1 from the p=0 contribution. This is a small but real gap in the written induction.","section":"Section 3, Lemma 2 (definition of c_p)"}],"minor_comments":[{"comment":"For α=1, the displayed solution y(z)=y0/(r−z) does not satisfy (17) with A(z)=r/(r−z) unless r=1; the correct solution is y(z)=y0(1−z/r)^{−r}. The conclusion that the radius remains r is still correct, but the explicit formula should be fixed.","section":"Example 1, p. 11"},{"comment":"The sentence 'Under Assumption (B), one achieves the recursion formula (2 1)' is misleading because (21) is the classical-derivative majorant problem; the fractional majorant is (22)–(23), and the Γ-ratio recursion should be stated explicitly to avoid conflating the two cases.","section":"Section 3.1, pp. 14-15"},{"comment":"The intermediate Laurent-type functions Φ_{k-1}(z)=z^{-k}∑ C_{jk}∆^{k-1}E(λ_j,z) are not holomorphic at the origin, and the step 'applying classical derivation recurrently' is not written precisely. The subsequent Vandermonde argument on the first r Taylor coefficients gives the independence directly and could replace this part.","section":"Theorem 4, proof of linear independence, p. 20"},{"comment":"The abstract contains the typo 'stablished'; also the notation M_α=(α!^α) in the paragraph after Definition 1 should presumably be M_α=(p!^α). These should be corrected in the final version.","section":"Abstract and Section 2.1"}],"recommendation":"major_revision","confidential_remarks":"The paper's main gap is confined to the proof of Lemma 2 and appears repairable by using the formal majorant solution of (11) directly. I saw no circularity or hidden fitted parameters. The heavy self-citation is consistent with the paper being a continuation of [15]."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Candid take: the paper extends your constant-coefficient work to analytic variable coefficients and gives two genuinely new pieces — Theorem 2 (translation invariance forces the factorial sequence) and Theorem 3 (transforming a system into a higher-order moment equation under a cyclic-vector hypothesis). The central existence claim, Theorem 1, is plausible and fits the expected shape, but the proof of Lemma 2 has a real gap. The majorant sequence c̃_p is defined by the ordinary-derivative recursion (8), then claimed to solve the moment problem (11), whose coefficient recursion requires the factor Γ(1+αp)/Γ(1+α(p+1)). Those are different for α<1, and the comparison c_p ≤ c̃_p fails — for m_p=Γ(1+p/2), A=1/(1-z), the intended c_1 exceeds c̃_1. The stress-test note is on target. This is fixable: define the majorant as the formal solution of (11) and then apply the Caputo integral-equation argument to that sequence. As printed, the proof does not establish the radius bound.\n\nThere is also a missing invertibility hypothesis in Lemma 5: A(z)=(∂_m Y)Y^{-1} requires det Y(0)≠0, and without it the lemma is false (Y(z)=zI is a counterexample). The recursion for c_p also leaves c_1 undefined since it is stated for p≥1. These are smaller, easy repairs.\n\nWhat the paper does well: the setup is careful, the relation to Caputo fractional systems and q-difference moments is explicit, Theorem 2 is a neat observation, and the worked examples help. The constant-coefficient adaptation in Section 4 is a reasonable echo of [15].\n\nWho is this for: specialists in moment summability and fractional differential equations. For that audience it deserves referee time — the main theorem is likely true and the gaps are specific, not fatal. I would send it out, with a request that the author fix the Lemma 2 comparison before acceptance. I would not cite it in my own work until the correction appears.","headline":"Plausible and new variable-coefficient existence result, but Lemma 2's majorant comparison is wrong as printed; the paper deserves revision and then a serious referee.","tokens_in":20331,"tokens_out":6649,"would_cite":false,"duration_ms":57567,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["30D15","34M03","34A08"],"pacs":[],"model":"deepseek-v4-flash","headline":"Every linear system of moment differential equations with analytic coefficients has a unique local analytic solution, with the radius controlled by the growth of the moment sequence.","keywords":["moment derivative","moment differential equations","strongly regular sequences","generalized summability","Caputo fractional derivative","local analytic solution","linear systems"],"falsifier":"For the scalar example (17) with $0<\\alpha<1$, compute the formal coefficients and the two bounding sequences $c_p$ and $\\tilde c_p$ used in Lemma 2: if $c_p \\le \\tilde c_p$ fails while the solution still has the claimed radius, the theorem survives but the printed proof does not; if the radius differs from the claimed value $r_0$, the theorem's bound is wrong.","tokens_in":1673,"feed_emoji":"🧮","tokens_out":6045,"duration_ms":118449,"temperature":0.7,"pith_summary":"This paper studies first-order linear systems of moment differential equations $\\partial_m y = A(z)y + b(z)$ whose coefficients are analytic near the origin. Its main claim is that, whenever the moment sequence is strongly regular and arises from a pair of kernel functions for $M$-summability, the Cauchy problem has a unique solution whose components are analytic on a disc around the origin. The proof separates two growth regimes: sequences that grow at least like a factorial power preserve the radius of convergence of the data, while sequences growing like a fractional power of the factorial (the Caputo-type case) give a smaller radius determined by the data and the parameter $\\alpha$. This unifies existence results for ordinary, $q$-difference, and fractional Caputo differentiation, which are all cases of moment differentiation.","feed_headline":"Moment differential systems have unique local analytic solutions","feed_subtitle":"Radius of solution disc is set by how fast the moment sequence grows.","key_machinery":"The central object is the moment derivative $\\partial_m$, defined by $\\partial_m\\bigl(\\sum_{p\\ge 0} a_p z^p/m_p\\bigr) = \\sum_{p\\ge 0} a_{p+1} z^p/m_p$. The argument writes the unknown as a formal series, derives the coefficient recursion, and bounds the coefficients by a scalar majorant series. For factorial-type sequences the majorant solves an ordinary differential equation; for fractional-type sequences it solves the Caputo fractional integral equation on the real axis, using the relation between moment differentiation and Caputo differentiation for $m = (\\Gamma(1+\\alpha p))_{p\\ge0}$. The strongly regular condition, together with kernel functions for $M$-summability, guarantees that the moment sequence is comparable to a factorial power and thus falls into one of the two growth regimes.","core_discovery":"The central claim is Theorem 1: if $M$ is a strongly regular sequence admitting a pair of kernel functions for $M$-summability, and $m$ is the associated moment sequence, then the Cauchy problem $\\partial_m y = A(z)y$, $y(0) = y_0$, with $A$ analytic near $0$, has a unique solution whose components are analytic on some disc $D(0,R)$. The proof splits into two regimes: sequences whose successive quotients grow at least linearly (Assumption (A)) give $R$ equal to the radius of $A$, while sequences of fractional-factorial growth (Assumption (B), $0 < \\alpha < 1$) give a possibly smaller radius computed from $\\alpha$ and the data. The same conclusion extends to nonhomogeneous systems. The paper also shows that, unless the moment sequence is essentially factorial, the moment derivative does not commute with translation, so the local result at $0$ cannot be transported to other base points by the classical shift.","pith_inferences":["The printed comparison of bounding sequences in Lemma 2 has a gap; the intended route through the Caputo integral equation would repair it, but an extra argument is required.","The absence of a Leibniz rule for moment derivatives (noted in the paper) makes a closed-form variation-of-constants formula unlikely; the nonhomogeneous result is obtained by majorants, not by such a formula.","Theorem 2 suggests a quantitative measure of 'momentness': the defect in the translation identity for a non-factorial sequence could be compared with the radius loss in Assumption (B), giving a testable relation between algebraic commutativity and analytic radius.","The cyclic-vector reduction in Theorem 3 is likely sharp: if the matrix $A(0)$ has no cyclic vector, no companion-form reduction to one scalar moment equation should exist, by the standard linear-algebra obstruction."],"forward_implications":["Every strongly regular moment sequence falls into Assumption (A) or (B), so the local existence theorem covers all such sequences and not only the explicit Gevrey-type examples.","Under Assumption (A) the solution's disc of convergence is at least the disc of analyticity of $A$ and $b$; no radius loss occurs.","Under Assumption (B) the radius shrinks to an explicit value determined by $\\alpha$ and the data, matching the picture for Caputo fractional systems.","For constant-coefficient homogeneous equations the solution space has an explicit basis of entire functions built from the kernel function $E$, with growth in generalized order and type governed by the roots of the characteristic polynomial.","Translation-commutation fails except for the factorial moment sequence, so the origin-centered theory cannot be moved to arbitrary points by the classical shift trick."],"supporting_citations":[{"why":"Supplies the definition and properties of strongly regular sequences, including equivalence between the sequence $M$ and its moment sequence $m$.","marker":"[25]"},{"why":"Guarantees the existence of proximate orders and kernel function pairs for strongly regular sequences.","marker":"[12]"},{"why":"Preceding constant-matrix study whose solution space and $\\Delta^h E$ functions are extended and adapted here.","marker":"[15]"},{"why":"Provides the lower factorial-type bound that lets Theorem 1 split into Assumption (A) or (B).","marker":"[24]"},{"why":"Gives the identity identifying moment derivation with Caputo fractional differentiation for $m=(\\Gamma(1+\\alpha p))$.","marker":"[21]"},{"why":"Supplies the fractional calculus tools that convert the Caputo majorant problem into an integral equation.","marker":"[13]"}],"fun_headline_variants":["Moment differential systems have unique local analytic solutions","Unique local analytic solutions proved for moment differential systems","Moment differential systems: unique analytic solutions exist locally","Local analytic solvability proven for moment differential systems","Moment derivative systems: unique analytic solution near zero"],"cache_read_input_tokens":22272,"weakest_assumption_plain":"The printed proof of the fractional case relies on comparing the coefficient bound to a classical-derivative recursion, a comparison that is not justified as written; the existence claim may still hold via the Caputo majorant equation, but that repair is not the printed argument.","fun_headline_variants_meta":{"raw":{"variants":["Moment differential systems have unique local analytic solutions","Unique local analytic solutions proved for moment differential systems","Moment differential systems: unique analytic solutions exist locally","Local analytic solvability proven for moment differential systems","Moment derivative systems: unique analytic solution near zero"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000203,"raw_usage":{"total_tokens":1295,"prompt_tokens":762,"completion_tokens":533,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":378,"completion_tokens_details":{"reasoning_tokens":459}},"tokens_in":378,"tokens_out":533,"duration_ms":5493,"temperature":1.0,"reasoning_tokens":459,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T21:38:11.920965+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"For the scalar example (17) with $0<\\alpha<1$, compute the formal coefficients and the two bounding sequences $c_p$ and $\\tilde c_p$ used in Lemma 2: if $c_p \\le \\tilde c_p$ fails while the solution still has the claimed radius, the theorem survives but the printed proof does not; if the radius differs from the claimed value $r_0$, the theorem's bound is wrong.","supporting_citations":[{"cited_title":"Sanz, Flat functions in Carleman ultraholomorphic classes via pro ximate orders, J","cited_arxiv_id":null,"evidence_quote":"Supplies the definition and properties of strongly regular sequences, including equivalence between the sequence $M$ and its moment sequence $m$."},{"cited_title":"Jim´ enez-Garrido, J","cited_arxiv_id":null,"evidence_quote":"Guarantees the existence of proximate orders and kernel function pairs for strongly regular sequences."},{"cited_title":"Lastra, Entire solutions of linear systems of moment differential eq uations and related asymptotic growth at infinity, Differ","cited_arxiv_id":null,"evidence_quote":"Preceding constant-matrix study whose solution space and $\\Delta^h E$ functions are extended and adapted here."},{"cited_title":"Petzsche, On E","cited_arxiv_id":null,"evidence_quote":"Provides the lower factorial-type bound that lets Theorem 1 split into Assumption (A) or (B)."},{"cited_title":"Michalik, Multisummability of formal solutions of inhomogeneous lin ear partial differen- tial equations with constant coefficients , J","cited_arxiv_id":null,"evidence_quote":"Gives the identity identifying moment derivation with Caputo fractional differentiation for $m=(\\Gamma(1+\\alpha p))$."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the fractional calculus tools that convert the Caputo majorant problem into an integral equation."}],"review_version":1}