{"id":"2553f8ff-e5bf-472c-a9d5-3bc0d4160502","arxiv_id":"1909.00492","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Every nonnegative solution of the static Schrödinger-Hartree-Maxwell equations with higher-order or fractional Laplacians is either zero, or in the critical case an explicit rescaled bubble.","lead":"Researchers classified all nonnegative solutions to a family of equations mixing fractional Laplacians with Hartree-type convolutions. In the critical case the only nonzero solutions are explicit rescaled bumps; in the subcritical and high-order critical cases only the zero solution exists.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The proof of Theorem 1.3 does not cover the m=0 case that the theorem statement includes, so the PDE classification in Corollary 1.5 is not established for the m=0 range as written.","rationale":"The integral-equation classification in Theorem 1.4 is largely self-contained and has independent support, including the explicit bubble computation via formula (4.47). The soft spot is the PDE-IE equivalence, exactly as the reader identified: the m=0 case is included in the statement of Theorem 1.3 and Corollary 1.5, but the proof in Section 3 starts from iterates u_i defined only for m>=1 and never writes the direct m=0 Green-Poisson/Liouville argument. This is addressable and likely fixable, but until it is written, Corollary 1.5 is not fully proved in the stated range. The CONDITIONAL verdict is therefore appropriate. I did not find a more severe structural flaw in the Theorem 1.4 moving-spheres proof: the Hardy-Littlewood-Sobolev exponents in (4.17) are consistent once the kernel convention of Lemma 4.1 is applied, and the subcritical-case wording 'WLOG tau>0 and mu>0' is a minor strict-inequality patch rather than a genuine counterexample.","tokens_in":32267,"tokens_out":35100,"duration_ms":439417,"concrete_test":"Add to Section 3 a self-contained proof for m=0: for a nonnegative classical solution u of (-Delta)^{alpha/2}u = (|x|^{-sigma} * u^p) u^q, set f = (|x|^{-sigma} * u^p) u^q and define v_R(x) = int_{B_R} G_R^alpha(x,y) f(y) dy; use Lemma 3.1 to show u >= v_R, let R -> infinity to get u >= I_alpha(f), then Lemma 3.2 gives u = I_alpha(f) + C; use integrability int u^p/|x|^sigma < infinity to force C=0 and recover (1.5). If this proof reproduces (3.45) for m=0, the concern is resolved; if the additive constant cannot be eliminated, the m=0 classification fails.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Theorem 1.3 is the bridge that turns the integral-equation classification (Theorem 1.4) into the PDE classification (Corollary 1.5). Its statement includes m=0 with alpha in (0,2] and p>0, but the proof in Section 3 defines u_i := (-Delta)^{i-1+alpha/2}u for i=1,...,m and proves the representation (3.1) for u_m; for m=0 these objects do not exist. The later iteration (3.19)-(3.29) also runs over m-1 steps. The final Green-Poisson/Liouville step (3.30)-(3.45) could in principle be applied directly with f_1(u) as the right-hand side, but that direct m=0 proof is never written. Consequently, the PDE-IE equivalence for m=0 is an unproved claim in the manuscript, not a purely notational variant. This matters because Corollary 1.5 includes m=0 and thereby claims classification for a parameter range not covered by the written proof. The fractional higher-order m>=1 case additionally relies on the same-group preprint [5] for part of the super poly-harmonic property (Remark 1.2), but the decisive gap is the explicit m=0 omission in the equivalence proof.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies nonnegative classical solutions of the higher-order/fractional static Schrödinger-Hartree-Maxwell equation (1.1). It establishes super poly-harmonic properties (Theorem 1.1), proves an equivalence between the PDE and the integral equation (1.5) (Theorem 1.3), classifies all nonnegative continuous solutions of the integral equation (Theorem 1.4), and transfers this classification to PDEs (Corollary 1.5). It also proves a Liouville theorem in critical and super-critical order cases (Theorem 1.8) and derives the best constant for a Hardy-Littlewood-Sobolev inequality (Corollary 1.7). The central classification result states that, for 0<s<n/2 in the critical exponent case, every nontrivial solution has the explicit form μ^{(n-2s)/2} Q(μ(x-x0)); in the subcritical cases only the zero solution exists.","tokens_in":32400,"tokens_out":19188,"duration_ms":171286,"significance":"The integral-equation classification is the core achievement and is largely convincing: the moving-sphere argument is detailed, the small-radius starting estimates (4.27)-(4.30) are standard, and the final profile is verified by direct substitution using the explicit identity (4.47). If the PDE-to-IE bridge is completed for all stated parameter ranges, this substantially generalizes earlier results by Liu, Cao-Dai, Dai-Fang-Qin, Dai-Liu, and others, and gives a complete classification over the full range of n, s, σ, p, q. The explicit best-constant formula in Corollary 1.7 is a useful byproduct. The main weakness is that Theorem 1.3, the bridge used for the PDE classification, omits the m=0 case in its written proof; this is a fixable but load-bearing gap.","major_comments":[{"comment":"The case m=0 is included in the statement of Theorem 1.3 and in Corollary 1.5, but the proof begins with the definition u_i := (-Δ)^{i-1+α/2}u for i=1,...,m and then proves the representation (3.19) for u_m, followed by the iteration (3.27)-(3.29) over k=1,...,m-1. For m=0 none of these objects is defined, so the PDE-to-IE representation is not proved in the range m=0. The direct Green-Poisson step that would handle m=0, namely applying the argument of (3.30)-(3.45) with f_1(u) as the right-hand side, is not written. Because Corollary 1.5 explicitly covers m=0, this is a load-bearing gap and not a purely notational issue.","section":"Section 3 (Theorem 1.3)"},{"comment":"In the subcritical cases the text says 'Without loss of generality, suppose that τ>0 and μ>0.' This is not a WLOG reduction: if p is subcritical while q is critical then μ>0 and τ=0, and if q is subcritical while p is critical then τ>0 and μ=0. The displayed strict inequalities in (4.45) use both (λx0/|z-x0|)^μ - 1 > 0 and (λx0/|y-x0|)^τ - 1 > 0; when one exponent is critical, one of these factors vanishes. The mixed cases should be proved explicitly, for example by keeping only the surviving positive factor, since the subcritical classification is the conclusion being established.","section":"Section 4 (subcritical case, Eq. (4.45))"}],"minor_comments":[{"comment":"The abstract states n≥1, but the PDE classification Corollary 1.5 assumes n≥2, Theorem 1.1 assumes n≥2, and Theorem 1.3 assumes n≥2; the abstract should be aligned with the statements.","section":"Abstract and Corollary 1.5"},{"comment":"The proof of Theorem 1.3 establishes the PDE-to-IE direction only; the asserted converse IE-to-PDE direction is not proved, though it is standard via Riesz potential properties. A sentence or short argument should be added for completeness.","section":"Section 3"},{"comment":"Proposition 4.3 states the identity u_{x0,λx0}(x)=u(x) only for x in B_{λx0}(x0)\\setminus{x0}, but Eq. (4.43) uses it for all x∈R^n\\setminus{x0}. The extension by Kelvin reflection should be justified either in the proposition or in the proof of Lemma 4.2.","section":"Section 4 (Proposition 4.3 and Eq. (4.43))"},{"comment":"For 0<α<2 the proof of Theorem 1.1 refers to the same-group preprint [5] for the key integral estimates leading to (2.18), and the induction for the remaining layers is only summarized by 'through a similar argument'. Since this property is used in the PDE-to-IE equivalence, the proof should either be fully self-contained or the precise result from [5] should be quoted with all hypotheses in force.","section":"Section 2 (Remark 1.2)"},{"comment":"There are several minor typographical issues, such as inconsistent use of |u|^p versus u^p in places where u≥0, and notational overload of the constant C in different estimates; these do not affect the mathematics but should be cleaned up.","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":"The m=0 gap in Theorem 1.3 is real and directly affects the PDE classification in Corollary 1.5, so I cannot recommend acceptance in the present form. The missing argument is short and standard, so a careful revision should be able to close it without changing the paper's scope. I would also ask the editor to monitor the status of the same-group preprint [5], on which part of the fractional higher-order super poly-harmonic property relies."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This paper is a serious extension of prior classification results for static Schrödinger-Hartree-Maxwell type equations to the full parameter range of (m, α, σ, p, q). The main payoff is Theorem 1.4, which classifies nonnegative solutions of the integral equation (1.5) via moving spheres, and Corollary 1.5, which transfers that to the PDE (1.1) through the equivalence in Theorem 1.3. The proof of Theorem 1.1 for integer higher-order cases (α = 2) is new and clean, and the explicit bubble verification plus the derived HLS best constant give real independent support. The moving-sphere argument in Theorem 1.4 is standard in structure and largely sound as far as I can tell.\n\nThe soft spots are real but not fatal. The stress-test is correct: the proof of Theorem 1.3 does not cover m = 0. The iteration defines u_i for i = 1,...,m and runs over m−1 steps; for m=0 there are no such functions, and the direct Green–Poisson step to (1.5) is never written. Since Corollary 1.5 explicitly includes m=0, the PDE classification for that parameter range is not established as written. This is a missing paragraph rather than a deep flaw, but it must be fixed. Also, for 0<α<2 and m≥1, the super poly-harmonic property (Theorem 1.1) and the Liouville ingredient in Theorem 1.8 are imported from the same-group preprints [5,6]. That is a self-containment problem, not circularity, but a referee will want those dependencies either proved here or clearly stated as external results. One minor imprecision: in the subcritical case the authors say \"without loss of generality, suppose τ>0 and μ>0\" when only one needs to be positive; the argument still works, so this is harmless but should be tightened.\n\nWho gets value: researchers working on classification for nonlocal elliptic equations, sharp HLS inequalities, or ground states for Hartree-type equations. I would send this to peer review rather than desk reject it. The main results are likely correct and important, but the m=0 case and the preprint dependencies need to be resolved before the paper is fully acceptable.","headline":"A strong full-range classification paper with a real but fixable gap: the PDE-to-integral-equivalence proof skips the m=0 case, and fractional-order ingredients lean on a same-group preprint.","tokens_in":33068,"tokens_out":3149,"would_cite":true,"duration_ms":26852,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35B53","35J30","35J91","35B06"],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper proves that in the critical case all nonnegative solutions of the fractional Schrödinger–Hartree–Maxwell equation are scaled translates of a single explicit profile, and that subcritical nonnegative solutions vanish identically.","keywords":["higher-order fractional Laplacians","Schrödinger-Hartree-Maxwell equations","classification of nonnegative solutions","super poly-harmonic properties","method of moving spheres","integral equations","Hardy-Littlewood-Sobolev inequality","Liouville theorem"],"falsifier":"Search for a nonnegative classical solution of (1.1) in the subcritical range that is not identically zero; even a single numerical example with admissible parameters (for instance $n=4$, $s=1$, $\\sigma=2$, $p=2$, $q=1$) would contradict the classification. Alternatively, compute $( -\\Delta)^{i+\\alpha/2}Q$ for the explicit $Q$ and check for a sign change, which would invalidate the super poly-harmonic premise behind the PDE-to-integral-equation transfer.","tokens_in":31947,"feed_emoji":"⚛️","tokens_out":8678,"duration_ms":73710,"temperature":0.7,"pith_summary":"This paper establishes a complete classification of nonnegative solutions to a family of higher-order and fractional Schrödinger–Hartree–Maxwell equations in which a fractional Laplacian is coupled to a Hartree-type nonlinearity through a Riesz potential. The authors first show that every nonnegative classical solution has all intermediate fractional Laplacians nonnegative, the super poly-harmonic property, and then prove that the PDE is equivalent to a single integral equation. Applying the method of moving spheres in integral form, they show that in the critical exponent range every nonzero solution is a scaled and translated copy of one explicit profile, while in the subcritical range the only nonnegative solution is zero. This classification transfers back to the PDE, and a parallel argument gives a Liouville theorem in critical and super-critical order cases. The result pins down the extremal functions and best constants for the associated Hardy–Littlewood–Sobolev inequality.","feed_headline":"One explicit profile classifies every critical Hartree solution","feed_subtitle":"Subcritical cases have only the zero solution; critical cases reduce to one explicit profile and give the sharp constant.","key_machinery":"The load-bearing object is the Kelvin transform $u_{x_0,\\lambda}(x)=(\\lambda/|x-x_0|)^{n-2s}u(x_\\lambda)$ centered at $x_0$, with $x_\\lambda=x_0+\\lambda^2(x-x_0)/|x-x_0|^2$, applied to the integral equation. The method of moving spheres starts from small $\\lambda$ where $u_{x_0,\\lambda}\\ge u$ in the ball $B_\\lambda(x_0)$, increases $\\lambda$ up to a critical scale $\\lambda_{x_0}$, and uses a contradiction argument to prove that at a finite critical scale the Kelvin transform coincides with $u$, forcing the explicit conformal profile; if no finite critical scale exists, the solution must be constant, which the integrability condition rules out. A calculus lemma classifies functions invariant under all such Kelvin transforms as the one-parameter family of $Q$. The earlier sections supply the two ingredients that make this legitimate for the PDE: the super poly-harmonic sign conditions and the equivalence theorem that converts the higher-order fractional PDE into the integral equation by iterated Riesz potentials.","core_discovery":"The central discovery is that the integral equation (1.5), and hence the PDE (1.1) through the equivalence theorem, admits only the explicit one-parameter family in the critical case. For $n\\ge 1$, $0<s:=m+\\alpha/2<n/2$, $0<\\sigma<n$, every nonnegative continuous solution with $p=(2n-\\sigma)/(n-2s)$ and $q=(n+2s-\\sigma)/(n-2s)$ is either identically zero or $$u(x)=\\$mu^{{\\frac{n-2s}}${2}}Q(\\mu(x-x_0)),\\qquad Q(x)=\\left(\\frac{1}{R_{2s,n}I(\\$\\sigma$/2)I((n-2s)/2)}\\right)^{\\frac{n-2s}{2(n+2s-\\$\\sigma$)}}\\left(\\frac{1}{1+|x|^2}\\right)^{\\frac{n-2s}{2}},$$ with $\\mu>0$ and $x_0\\in\\mathbb{R}^n$. In the subcritical cases $0<p<(2n-\\sigma)/(n-2s)$ or $0<q<(n+2s-\\sigma)/(n-2s)$, the only nonnegative solution is $u\\equiv 0$. The proof routes through three steps: super poly-harmonic inequalities, the equivalence between the PDE and the integral equation, and the moving-spheres classification of the integral equation; the explicit $Q$ then yields the best constant of the corresponding Hardy–Littlewood–Sobolev inequality.","pith_inferences":["If the equivalence and classification hold, the same explicit $Q$ should control sharp constants in weighted Hardy–Littlewood–Sobolev inequalities across the full admissible parameter range, not only the special cases computed before.","The moving-spheres proof provides a template for classifying solutions of systems of coupled integral equations with different exponents $p,q$, a natural next step for multi-component Hartree systems.","A direct check suggested by the proof is whether the explicit $Q$ satisfies the super poly-harmonic inequalities $( -\\Delta)^{i+\\alpha/2}Q\\ge 0$ for every intermediate order $i$; computing these for representative parameters would test the consistency of the equivalence theorem with the classification."],"forward_implications":["In the critical case, every nonzero nonnegative solution of the PDE is a scaling-translation of the explicit profile $Q$; in the subcritical case the zero solution is unique.","The explicit $Q$ is the unique extremal function of the associated Hardy–Littlewood–Sobolev inequality, and the best constant $S_{\\sigma,s,n}$ can be computed in closed form.","For $s\\ge n/2$, no nonzero nonnegative classical solution exists for any admissible $p,q$ and $\\sigma<n$ (Liouville theorem).","The PDE–integral-equation equivalence means classification results proved for integral equations automatically transfer to the original fractional Laplacian problem, including the nonlocal-nonlocal interaction.","Earlier classifications for special parameter choices are subsumed by one statement covering the full range of $n$, $s$, $\\sigma$, $p$ and $q$."],"supporting_citations":[{"why":"Supplies the super poly-harmonic property for fractional higher-order cases, used in Theorem 1.1.","marker":"[5]"},{"why":"Supplies the Liouville and representation arguments used to prove the PDE–integral-equation equivalence and Theorem 1.8.","marker":"[6]"},{"why":"Supplies the moving-spheres calculus lemma classifying functions invariant under Kelvin transforms, used in Theorem 1.4.","marker":"[47]"},{"why":"Provides the Riesz potential constants and the composition formula used to collapse iterated convolutions into (1.5).","marker":"[54]"},{"why":"Provides identity (4.47), which fixes the constant $C$ in the explicit profile $Q$.","marker":"[27]"},{"why":"Prior classification of a second-order Hartree case that the paper extends.","marker":"[43]"},{"why":"Prior classification of a bi-harmonic Hartree case that the paper extends.","marker":"[4]"},{"why":"Identifies the Hardy–Littlewood–Sobolev inequality and its extremals, which the explicit profile saturates.","marker":"[39]"}],"fun_headline_variants":["Critical Hartree solutions: one explicit family only","Subcritical zero, critical one explicit profile","All critical solutions reduce to a single explicit form","Hartree equation: only zero or explicit profile","Explicit profile classifies all critical Hartree solutions"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The classification for PDEs rests on the super poly-harmonic property—every nonnegative classical solution has $( -\\Delta)^{i+\\alpha/2}u\\ge 0$ for all intermediate orders; if that property fails in any admissible parameter range, the integral-equation classification no longer transfers to the PDE.","fun_headline_variants_meta":{"raw":{"variants":["Critical Hartree solutions: one explicit family only","Subcritical zero, critical one explicit profile","All critical solutions reduce to a single explicit form","Hartree equation: only zero or explicit profile","Explicit profile classifies all critical Hartree solutions"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000254,"raw_usage":{"total_tokens":1711,"prompt_tokens":1232,"completion_tokens":479,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":848,"completion_tokens_details":{"reasoning_tokens":408}},"tokens_in":848,"tokens_out":479,"duration_ms":64528,"temperature":1.0,"reasoning_tokens":408,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T05:47:51.394732+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Search for a nonnegative classical solution of (1.1) in the subcritical range that is not identically zero; even a single numerical example with admissible parameters (for instance $n=4$, $s=1$, $\\sigma=2$, $p=2$, $q=1$) would contradict the classification. Alternatively, compute $( -\\Delta)^{i+\\alpha/2}Q$ for the explicit $Q$ and check for a sign change, which would invalidate the super poly-harmonic premise behind the PDE-to-integral-equation transfer.","supporting_citations":[{"cited_title":"Super poly-harmonic properties, Liouville theorems and classification of nonnegative solutions to equations involving higher-order fractional Laplacians","cited_arxiv_id":"1905.04300","evidence_quote":"Supplies the super poly-harmonic property for fractional higher-order cases, used in Theorem 1.1."},{"cited_title":"Li and L","cited_arxiv_id":null,"evidence_quote":"Supplies the moving-spheres calculus lemma classifying functions invariant under Kelvin transforms, used in Theorem 1.4."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the Riesz potential constants and the composition formula used to collapse iterated convolutions into (1.5)."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides identity (4.47), which fixes the constant $C$ in the explicit profile $Q$."},{"cited_title":"Liu, Regularity, symmetry, and uniqueness of some integral type quasilinear equations, Nonlinear Anal., 71 (2009), 1796-1806","cited_arxiv_id":null,"evidence_quote":"Prior classification of a second-order Hartree case that the paper extends."},{"cited_title":"Cao and W","cited_arxiv_id":null,"evidence_quote":"Prior classification of a bi-harmonic Hartree case that the paper extends."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Identifies the Hardy–Littlewood–Sobolev inequality and its extremals, which the explicit profile saturates."}],"review_version":1}