{"id":"34f10969-05c9-44f7-972a-3544d37eea6c","arxiv_id":"1908.07158","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A new multivariable confluent hypergeometric function yields explicit fundamental solutions with r^{2-m} singularity for the Helmholtz equation with n singular Bessel coefficients.","lead":"This paper introduces a new family of multivariable confluent hypergeometric functions and uses them to write explicit fundamental solutions for a generalized Helmholtz equation with Bessel-type singular coefficients. A generalist reader may care because it offers a unified formula for singular solutions of a class of elliptic PDEs and clarifies the order of the singularity.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Formula (5.5) is internally inconsistent with system (5.4): the first parameter of H should be α, not \\tildeα_k, so the claimed fundamental solutions are not shown to solve (1.1).","rationale":"The reader's weakest assumption was about completeness of the separated ansatz and of the hypergeometric solution family. My reading confirms that the derivation from (5.4) to (5.5) is the soft spot, but the issue is more acute than a missing completeness proof: the parameter α in system (5.4) is fixed by the PDE reduction, while formula (5.5) uses \\tildeα_k as the first argument of the confluent hypergeometric function. This is not a notational variant because the first argument enters the Pochhammer symbols and the PDE system (4.2). The prefactor ξ_i^{1-c_i} and the definition ξ_i = -4x_i x0_i/r^2 only account for the radial exponent r^{-2\\tildeα_k}; they do not change the first argument of H. Thus, as written, the q_k in (5.5) do not satisfy (5.4) and hence are not established as solutions of (1.1). The singularity computation in Section 6 uses the same erroneous parameter, so it does not rescue the claim. A single symbolic/numeric substitution for n=1 would settle the matter; the algebraic comparison already shows the residual is nonvanishing generically. Because the central application of the paper depends on these explicit solutions, the verdict should be REJECT rather than CONDITIONAL unless the formulas are corrected and re-verified.","tokens_in":12951,"tokens_out":16439,"duration_ms":155401,"concrete_test":"For n=1, m=3, α_1=1/4, compare (4.2) with (5.4). A direct calculation shows H(a,α_1;2α_1;ξ,η_1) satisfies system (5.4) iff (a-α)(ξ ω_ξ + α_1 ω)=0 and (a-α)ω_{η_1}=0 for all ξ,η_1. Here α=3/4 and \\tildeα_0=1/4, so evaluating at (ξ,η_1)=(1/2,-1) gives a nonzero residual. Equivalently, substitute q_0 from (5.5) into (1.1) with λ=1 at x=(1,1,1), x0=(2,3/2,1), using a truncated series for H; the residual will not vanish. This settles whether (5.5) is a genuine fundamental solution.","verdict_should_be":"REJECT","load_bearing_attack":"The central claim is that all fundamental solutions are explicitly given by (5.5). But (5.5) does not follow from the stated derivation. System (5.4) is exactly system (4.2) with a=α, b_i=α_i, c_i=2α_i. Therefore the solution family (4.3) applied to (5.4) yields ω_k = C_k ∏_{i=1}^k ξ_i^{1-2α_i} H^{(n,1)}_A[α, 1-α_1,...,1-α_k, α_{k+1},...,α_n; 2-2α_1,...,2-2α_k, 2α_{k+1},...,2α_n; ξ,η_1]. Using ξ_i = -4x_i x0_i/r^2, this becomes const · r^{-2\\tildeα_k} ∏_{i=1}^k (x_i x0_i)^{1-2α_i} H^{(n,1)}_A[α, ...]. The radial exponent r^{-2\\tildeα_k} is correct, but the first argument of H must remain α. Formula (5.5) instead puts \\tildeα_k there. The prefactor transformation cannot change the first argument of H. Unless \\tildeα_k = α, which fails generically (e.g. \\tildeα_0 = α - 2∑α_i), the functions in (5.5) are not the solutions of (5.4). Consequently, q_k(x,x0) as written is not shown to solve equation (1.1). The same mismatch appears in (6.1), so the singularity analysis in Section 6 applies to a different function than the one obtained from the reduction. This is a correctness gap in the explicit construction, more basic than the completeness question raised in the reader's verdict.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper introduces a new confluent hypergeometric function H^{(n,p)}_A of several variables, derives a system of partial differential equations satisfied by it, and establishes decomposition formulas. It then uses a separation ansatz to construct fundamental solutions of the generalized Helmholtz equation (1.1) with n singular coefficients, claiming that all such solutions are given explicitly by formula (5.5) and that they have a singularity of order r^{2-m} as r tends to 0 for m>2. A similar statement is made for equation (7.1) with several parameters λ_k, with solutions in formula (7.4).","tokens_in":13406,"tokens_out":15064,"duration_ms":135765,"significance":"If the construction were correct, the paper would provide a unified explicit family of fundamental solutions for a class of multidimensional elliptic equations with several singular coefficients, extending earlier work that is mostly limited to n≤3. The new hypergeometric function and its decomposition formulas are potentially useful, and the paper does give checkable series definitions and a derived PDE system. However, the central explicit formula (5.5) contains a parameter mismatch that breaks the connection to the derived system, and the singularity proof relies on an identity quoted from the author's previous paper without proof. The significance of the contribution is therefore contingent on the corrections described below.","major_comments":[{"comment":"The system (5.4) is exactly the system (4.2) with a=α, b_i=α_i and c_i=2α_i. Applying the solution family (4.3) to (5.4) gives ω_k = C_k ∏_{i=1}^k ξ_i^{1-2α_i} H^{(n,1)}_A[α, 1-α_1,...,1-α_k, α_{k+1},...,α_n; 2-2α_1,...,2-2α_k, 2α_{k+1},...,2α_n; ξ,η_1]. The first argument of H must therefore be α, not \\tildeα_k. Since \\tildeα_k = α + k - 2Σ_{i=1}^k α_i, equality holds only when k=0. Consequently, for k≥1 the functions in (5.5) do not satisfy system (5.4) and are not shown to solve equation (1.1). The same mismatch appears in (7.4). This is a load-bearing error because the claim that all fundamental solutions are given by (5.5) is the main result of the paper; the first parameter should be α, with the constants γ_k in (5.6) adjusted accordingly.","section":"§5, Eq. (5.5); §7, Eq. (7.4)"},{"comment":"The singularity analysis is carried out only for q0, and the decisive limit (6.4) uses the identity [21, p.94, (33)] without proof. Since this identity is quoted from the author's own previous paper and is nontrivial, the claimed singularity order r^{2-m} is not established within the manuscript. For qi with i=1,...,n, the text says 'similarly' without specifying the relevant parameters or verifying that the same summation identity applies after the correction in (5.5). Because the singularity order is a central advertised result, this gap must be addressed, either by proving the quoted identity or by supplying the missing calculations for all k.","section":"§6, Eq. (6.4)"},{"comment":"The derivation from the equation in (5.3) to the system (5.4) is compressed into a single sentence. The coefficients in (5.3) contain ratios x0_i/x_i, and it is not transparent how these cancel to yield the x-independent system (5.4). Since (5.4) is the basis for all explicit solutions, the authors should provide the full computation or a detailed outline showing that the separation ansatz leads exactly to this system.","section":"§5, reduction from (5.3) to (5.4)"},{"comment":"The paper asserts that all solutions of system (4.2) are expressed by (4.3), but this system is a coupled PDE system in several variables, and the list (4.3) contains only n+1 symmetry-grouped solutions. No proof of completeness is given, nor is it shown that no other linearly independent solution can contribute to a fundamental solution of (1.1). The abstract and Section 5 use the phrase 'all fundamental solutions'; this is stronger than what is demonstrated. The authors should either prove the needed completeness or explicitly qualify the statement to 'all fundamental solutions obtained by this separation ansatz.'","section":"§4, formula (4.3), and Abstract"}],"minor_comments":[{"comment":"The phrase 'substituting these solutions in (5)' is unclear because no equation (5) exists; it should refer to the ansatz preceding (5.1) or to equation (1.1).","section":"§5, sentence before Eq. (5.5)"},{"comment":"The text says 'substitute them into equation (1.1)', but it should say equation (7.1).","section":"§7, after Eq. (7.2)"},{"comment":"The summation in (5.3) uses index m while the coefficients are labelled A_k, B_mk, etc.; the notation should be made consistent.","section":"§5, Eq. (5.3)"},{"comment":"The equality H^{(n,p)}_A(...; η_1,...,η_p) = H^{(n,1)}_A(...; η_1+...+η_p) is stated without justification; it follows from the series expansion but should be stated with the appropriate convergence conditions.","section":"§2, after (2.2)"},{"comment":"The manuscript contains many typographical and language issues, such as 'сonsidered', 'F undamental', and inconsistent spacing in formulas; a careful proofreading is needed.","section":"General"},{"comment":"After replacing the first argument of H in (5.5) with α, the constants γ_k in (5.6) must be recomputed; as written they contain Γ(\\tildeα_k) and are inconsistent with a hypergeometric function whose first parameter is α except in the case k=0.","section":"§6, Eq. (5.6) after correction of (5.5)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript relies heavily on the author's own previous work: the solution family (2.5) is credited to [12], and the summation identity used in (6.4) is quoted from [21] without proof. The editor may wish to ask the author to include a proof of this identity or to state it as a lemma with a full derivation. The parameter mismatch in (5.5) is not a mere typo, because it also affects the constants γ_k and the singularity argument; the paper needs a careful reworking of Section 5 and the corresponding parts of Sections 6 and 7 before it can be considered for publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear Colleague,\n\nThe paper has one genuinely new piece: the confluent hypergeometric function H^{n,p}_A, a confluent limit of Erdélyi's function, with a PDE system and a Burchnall–Chaundy decomposition formula. That part is coherent and extends earlier work by Hasanov and Srivastava. The application to the generalized Helmholtz equation with n singular coefficients is the right generalization of the m=n case in [12], and the reduction from the PDE to the hypergeometric system (5.4) is standard and plausible.\n\nThe problem is that the final explicit formula does not follow. System (5.4) is exactly system (4.2) with a=α, b_i=α_i, c_i=2α_i. Formula (4.3) then gives ω_k = ∏ ξ_i^{1-2α_i} H^{n,1}_A[α, ...]. With ξ_i = -4x_i x0_i/r^2, the radial prefactor becomes r^{-2\\tilde α_k}, but the first argument of H remains α. Equation (5.5) puts \\tilde α_k there. For k=0 that is fine because \\tilde α_0=α; for k≥1 the two parameters differ generically, so the functions in (5.5) are not solutions of (5.4). The same issue carries into (7.4). The stress-test note is right about (5.5); its remark about (6.1) is not, since q0 is the k=0 case.\n\nOther soft spots are secondary: the derivative computation in Section 5 is compressed into one sentence, the singularity analysis for q_i is dismissed with 'similarly', and the key limit (6.4) uses identity [21] from the author's own earlier paper without proof. The completeness claim for 'all' fundamental solutions is also asserted rather than proved.\n\nIf the \\tilde α_k in the H brackets is a typo for α, the construction probably works, but as written it is internally inconsistent. This deserves a serious referee because the special-function material is salvageable and potentially useful, but I would not accept the current version, and I would not cite it until the parameters are fixed and the identity [21] is either proved or replaced.","headline":"The new special function and decomposition are worth knowing, but the explicit fundamental solution formula (5.5) does not follow from the hypergeometric system, so the paper needs correction before it can be used.","tokens_in":13903,"tokens_out":7338,"would_cite":false,"duration_ms":71128,"reading_group":"maybe","serious_thinker":"no","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35A08","33C65","33C70","35J05"],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper introduces a new class of many-variable confluent hypergeometric functions and uses it to write all fundamental solutions of the generalized Helmholtz equation with n singular coefficients explicitly, with singularity of order…","keywords":["confluent hypergeometric functions","Lauricella functions","fundamental solutions","generalized Helmholtz equation","singular coefficients","decomposition formula","symbolic operators","singularity order"],"falsifier":"For $m=3$, $n=1$, $\\alpha_1=1/4$ and $\\lambda=0$, compute $q_0$ from (5.5) and test distributionally whether $(\\Delta+(1/(2x_1))\\partial_{x_1})q_0=\\delta(x-x_0)$; a mismatch would show (5.5) is not a fundamental solution. Independently, evaluating the multivariable identity on [21, p.94] at $n=2$ would confirm or refute the predicted $r^{2-m}$ order.","tokens_in":12727,"feed_emoji":"🧮","tokens_out":12176,"duration_ms":106111,"temperature":0.7,"pith_summary":"The paper introduces a new family of confluent hypergeometric functions of many variables, $H_A^{(n,p)}$, defined as a confluent limit of a generalized hypergeometric $H$-function of several variables, and proves that this family supplies explicit fundamental solutions for the generalized Helmholtz equation with $n$ singular coefficients. The central claim is that every fundamental solution of equation (1.1) is captured by formula (5.5), a separated expression $u=P(r)\\omega(\\xi,\\eta_1)$ in which the angular part is one of the $2^n$ solutions of a hypergeometric system written with $H_A^{(n,1)}$. Using decomposition formulas that expand $H_A^{(n,1)}$ into one-variable hypergeometric functions, the paper shows each solution is singular of order $r^{2-m}$ as $r\\to 0$ for $m>2$, and logarithmically for $m=2$. The same construction is extended to equation (7.1), where several Helmholtz parameters $\\lambda_k$ appear, yielding fundamental solutions (7.4) through $H_A^{(n,p)}$. A reader would care because explicit fundamental solutions are the standard route to integral representations and boundary-value solvability for singular elliptic equations.","feed_headline":"One hypergeometric function yields every Helmholtz fundamental solution","feed_subtitle":"Explicit fundamental solutions for all n singular coefficients, blowing up like r^{2-m} as the radius tends to 0.","key_machinery":"The carrying object is the many-variable confluent hypergeometric function $H_A^{(n,p)}(a,b_1,\\ldots,b_n;c_1,\\ldots,c_n;\\xi,\\eta)$, defined by the multiple series (2.2) as a confluent limit of a generalized hypergeometric $H$-function of several variables; when $p=1$ it reduces to (4.1). The mechanism that connects this function to the differential equation is the separated ansatz $u(x,x_0)=P(r)\\omega(\\xi,\\eta_1)$ with $P(r)=(r^2)^{-\\alpha}$, $\\alpha=\\sum_{i=1}^n\\alpha_i-1+m/2$, and with the ratio variables $\\xi_j=(r^2-r_j^2)/r^2$ and the confluent variable $\\eta_1=-\\lambda^2 r^2/4$. Substitution turns equation (1.1) into the hypergeometric system (5.4), whose solutions are taken from the Kummer-type family (4.3); combining these with $P(r)$ yields (5.5). The singularity analysis relies on a symbolic-operator technique and the decomposition formula (4.5), which express $H_A^{(n,1)}$ as a sum of products of Lauricella-type functions and one-variable confluent functions; the limit $r\\to 0$ is then evaluated with the identity quoted from [21, p.94].","core_discovery":"The paper claims that all fundamental solutions of the generalized Helmholtz equation with several singular coefficients can be written explicitly as $q_k(x,x_0)=\\gamma_k\\prod_{i=1}^{k}(x_i x_{0i})^{1-2\\alpha_i} r^{-2\\widetilde{\\alpha}_k} H_A^{(n,1)}[\\widetilde{\\alpha}_k,1-\\alpha_1,\\ldots,1-\\alpha_k,\\alpha_{k+1},\\ldots,\\alpha_n; 2-2\\alpha_1,\\ldots,2-2\\alpha_k,2\\alpha_{k+1},\\ldots,2\\alpha_n;\\xi,\\eta_1]$ for $k=0,\\ldots,n$, where $\\gamma_k$ and $\\widetilde{\\alpha}_k$ are defined by (5.6) and (5.7). The variables are $\\xi_j=(r^2-r_j^2)/r^2$ and $\\eta_1=-\\lambda^2 r^2/4$, with $r_j$ the distance from $x$ to the reflection of $x_0$ in the $j$-th coordinate plane. The proof substitutes the ansatz $u=P(r)\\omega(\\xi,\\eta_1)$, reduces equation (1.1) to the hypergeometric system (5.4), and identifies its solutions with the Kummer-type family (4.3). The same argument for equation (7.1) gives (7.4) through $H_A^{(n,p)}$. As $r\\to 0$, each $q_k$ has a singularity of order $r^{2-m}$ for $m>2$, and the $m=2$ case is logarithmic.","pith_inferences":["A direct numerical or symbolic check of the identity quoted from [21, p.94] for small $n$ would independently confirm the singularity-order step, since the paper does not prove that identity itself.","Because $H_A^{(n,p)}$ reduces to $H_A^{(n,1)}$ by summing the $\\eta_j$ variables, the multi-parameter equation (7.1) is essentially a superposition of single-frequency fundamental solutions, which suggests Duhamel-type arguments for variable or nonlinear $\\lambda$.","The confluent variable is tied to $\\lambda^2 r^2/4$, so allowing complex $\\lambda$ should produce oscillatory or damped kernels directly, potentially useful for dissipative Helmholtz or Klein-Gordon-type equations.","The completeness of the separated ansatz $u=P(r)\\omega(\\xi,\\eta_1)$ is the natural next step to verify; if a fundamental solution requires a different radial factor, formula (5.5) would give only a particular subfamily."],"forward_implications":["Formula (5.5) gives explicit integral kernels for potential theory on the singular half-space $\\mathbb{R}^n_{+m}$, so boundary-value problems for equation (1.1) can be converted into integral equations with known kernels.","The uniform singularity order $r^{2-m}$ for $m>2$ matches the Laplacian's fundamental solution, so classical layer-potential estimates should transfer to this singular elliptic setting.","For $m=2$ the same family yields logarithmic singularity, resolving the borderline case of the construction.","Equation (7.1) with several parameters $\\lambda_1,\\ldots,\\lambda_p$ has fundamental solutions (7.4) through $H_A^{(n,p)}$ with the same singularity type, so the method covers multi-frequency Helmholtz-type problems.","The decomposition formulas (3.13) and (4.5) reduce the new confluent function to one-variable hypergeometric terms, making the kernels computable in practice."],"supporting_citations":[{"why":"Supplies the one-variable hypergeometric function, the summation formula at argument 1, and the multivariable definitions on which the new confluent function is built.","marker":"[1]"},{"why":"Provides the bi-axially symmetric Helmholtz fundamental solutions and the logarithmic-singularity statement for $m=2$ that the present work generalizes.","marker":"[6]"},{"why":"Gives the $m=n$ case and the $2^n$ Kummer-type solution family (4.3) used to solve the hypergeometric system (5.4).","marker":"[12]"},{"why":"Supplies the limiting device $(1/\\varepsilon)^q\\varepsilon^q=1$ used to define $H_A^{(n,p)}$ and a decomposition formula for the multivariable Lauricella function.","marker":"[15]"},{"why":"Introduces the multivariable symbolic operators used to derive the decomposition formula (4.5) for $H_A^{(n,1)}$.","marker":"[18]"},{"why":"Provides the multivariable identity used to evaluate the $r\\to 0$ limit in (6.4), the step that fixes the singularity order $r^{2-m}$.","marker":"[21]"}],"fun_headline_variants":["All Helmholtz fundamental solutions from a single hypergeometric function","Unified explicit solutions for generalized Helmholtz singular coefficients","New confluent hypergeometric class cracks singular Helmholtz equation","Explicit singular Helmholtz solutions from a new hypergeometric family"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The claim that (5.5) gives all fundamental solutions rests on the separated ansatz $u=(r^2)^{-\\alpha}\\omega(\\xi,\\eta_1)$ with $\\alpha=\\sum_i\\alpha_i-1+m/2$ and on the assertion that the Kummer-type family (4.3) exhausts the solutions of the resulting hypergeometric system; the singularity order additionally rests on an identity quoted from [21, p.94].","fun_headline_variants_meta":{"raw":{"variants":["All Helmholtz fundamental solutions from a single hypergeometric function","Unified explicit solutions for generalized Helmholtz singular coefficients","New confluent hypergeometric class cracks singular Helmholtz equation","Explicit singular Helmholtz solutions from a new hypergeometric family"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001314,"raw_usage":{"total_tokens":5370,"prompt_tokens":981,"completion_tokens":4389,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":597,"completion_tokens_details":{"reasoning_tokens":4318}},"tokens_in":597,"tokens_out":4389,"duration_ms":30583,"temperature":1.0,"reasoning_tokens":4318,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T12:23:56.904513+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"For $m=3$, $n=1$, $\\alpha_1=1/4$ and $\\lambda=0$, compute $q_0$ from (5.5) and test distributionally whether $(\\Delta+(1/(2x_1))\\partial_{x_1})q_0=\\delta(x-x_0)$; a mismatch would show (5.5) is not a fundamental solution. Independently, evaluating the multivariable identity on [21, p.94] at $n=2$ would confirm or refute the predicted $r^{2-m}$ order.","supporting_citations":[{"cited_title":"Erdelyi, W","cited_arxiv_id":null,"evidence_quote":"Supplies the one-variable hypergeometric function, the summation formula at argument 1, and the multivariable definitions on which the new confluent function is built."},{"cited_title":"Hasanov, Fundamental solutions bi-axially symmetric Helmholtz equation","cited_arxiv_id":null,"evidence_quote":"Provides the bi-axially symmetric Helmholtz fundamental solutions and the logarithmic-singularity statement for $m=2$ that the present work generalizes."},{"cited_title":"Urinov, T.G","cited_arxiv_id":null,"evidence_quote":"Gives the $m=n$ case and the $2^n$ Kummer-type solution family (4.3) used to solve the hypergeometric system (5.4)."},{"cited_title":"Appell, J.Kampe de Feriet, Fonctions Hypergeometriques et Hyperspheriques; Polynom es d’Hermite, Gauthier - Villars","cited_arxiv_id":null,"evidence_quote":"Supplies the limiting device $(1/\\varepsilon)^q\\varepsilon^q=1$ used to define $H_A^{(n,p)}$ and a decomposition formula for the multivariable Lauricella function."},{"cited_title":"Hasanov, H.M","cited_arxiv_id":null,"evidence_quote":"Introduces the multivariable symbolic operators used to derive the decomposition formula (4.5) for $H_A^{(n,1)}$."},{"cited_title":"Ergashev, The Dirichlet problem for elliptic equat ion with several singular coefficients, e-Journal of Analysis and Applied Mathematics, 2018 (1), p","cited_arxiv_id":null,"evidence_quote":"Provides the multivariable identity used to evaluate the $r\\to 0$ limit in (6.4), the step that fixes the singularity order $r^{2-m}$."}],"review_version":1}