{"id":"549545b9-622d-4a88-aad4-40f06ef6822f","arxiv_id":"1908.05719","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"For the subcritical fractional NLS with harmonic potential, ground states exist, are radial and decreasing, and are orbitally stable if the Cauchy problem is well posed.","lead":"The paper proves that a nonlinear fractional Schrödinger equation with a harmonic trapping potential has ground-state standing waves, and that these waves are orbitally stable when the initial-value problem is well posed. It also runs numerical simulations showing how the fractional power, the nonlinearity, and the trap shape the wave and its stability.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Theorem 2.1 rests on Lemma 3.1, whose proof invokes a false compact embedding and omits the s >= N/2 case; the subsequence extraction is unsupported as written.","rationale":"The reader's weakest-assumption analysis identifies exactly the load-bearing gap: Theorem 2.1's proof of attainment of the infimum in (8) requires strong convergence of a minimizing sequence in L^{2σ+2}, and that convergence is supplied only by Lemma 3.1. The lemma's proof contains a false statement about compactness of Hs on R^N, and it also fails to address the case s >= N/2, which is not excluded by the theorem's hypotheses. I therefore agree with the reader's conditional verdict: the main conclusions are plausible and likely correct, but the existence proof as written is not sound. I would not move to REJECT because the gap is standard and repairable: the weighted space with |x|^2 should restore compactness, and interpolation can extend p = 2 compactness to the needed range. I also considered the conditional well-posedness assumption in Theorem 4.1; while it is a genuine limitation, it is explicitly stated as a hypothesis, so it is less load-bearing than the invalid Lemma 3.1. The absence of code and data in the numerical section is a reproducibility concern but does not bear on the central mathematical claim. Thus the appropriate final verdict remains CONDITIONAL, matching the reader's verdict.","tokens_in":15316,"tokens_out":8390,"duration_ms":85739,"concrete_test":"Perform an independent proof check of Lemma 3.1. Write out the p = 2 compactness argument fully using the tail estimate (14) and a diagonal extraction over balls; then, for s < N/2, use Riesz-Thorin interpolation with the continuous embedding Hs(R^N) -> L^{2N/(N-2s)} to cover every 2 < p < 2N/(N-2s), and for s >= N/2 adapt the argument with finite q > p. Verify that the needed exponent p = 2σ+2 is covered for every 0 < σ < 2s/N allowed by Theorem 2.1. If the compactness proof cannot be completed, Theorem 2.1 lacks a valid proof.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The existence theorem 2.1 depends on Lemma 3.1, and that lemma's proof is invalid. The proof states that 'Hs(R^N) can be compactly embedded into Lp(R^N) for 2 < p < 2N/(N-2s)'; this is false on the unbounded domain R^N, where translations prevent compactness. For p = 2 the separate diagonal/tail argument in (14) can be made correct, and for p > 2 compactness would follow from compactness in L^2 combined with interpolation and the continuous embedding Hs -> L^{2N/(N-2s)}, but the paper supplies no such argument. There is also a missing case distinction: for s >= N/2, the expression 2N/(N-2s) is negative or undefined, yet Theorem 2.1 states 0 < s < 1 with no restriction, so Lemma 3.1's stated range is empty in that case. As written, the step in the proof of Theorem 2.1 where a minimizing sequence is extracted as strongly convergent in L^{2σ+2} is unsupported. The lemma is likely repairable, so I would not reject the theorem, but the proof needs correction before the existence claim is accepted.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the fractional nonlinear Schrödinger equation with a harmonic potential, Eq. (1). It introduces the weighted Hilbert space Σ_s(R^N), Eq. (10), and studies the constrained minimization of the energy functional J on the mass sphere S_c, Eq. (8). Theorem 2.1 claims that for 0<σ<2s/N this problem admits a nonnegative, radial, radially decreasing minimizer; Theorem 2.2 claims orbital stability of the associated ground states under a uniqueness/well-posedness assumption on the Cauchy problem. The second half of the paper develops split-step Fourier spectral and normalized-gradient-flow numerical methods, reports numerical ground states over a wider parameter range, and studies dynamics and stability numerically.","tokens_in":15556,"tokens_out":10289,"duration_ms":95797,"significance":"If the existence and stability theorems are established, the paper would provide a natural variational treatment of fractional NLS with a harmonic trap, extending known results for the Gross-Pitaevskii equation (s=1) and for fractional NLS without potential. The numerical section is broad and gives useful heuristics on the role of s, σ, and the trapping potential. The central existence proof, however, has a gap in Lemma 3.1: the compactness argument rests on a false Sobolev embedding statement, and no alternative proof is supplied. The result is likely repairable, and the stability theorem is explicitly conditional on an unproved well-posedness assumption, so the manuscript's claims are defensible as conditional statements but not yet fully proved as written.","major_comments":[{"comment":"The proof of compactness for p>2 relies on the assertion that H^s(R^N) is compactly embedded into L^p(R^N) for 2<p<2N/(N-2s). This is false on the unbounded domain R^N, where translation invariance prevents compactness. Moreover, the stated range 2≤p<2N/(N-2s) is empty when N≤2s, while Theorem 2.1 allows such values (e.g., N=1, s=0.8, σ=1). The lemma is used in the proof of Theorem 2.1 to extract a subsequence converging strongly in L^{2σ+2}, and again in Theorem 4.1. As written, these subsequence extractions are unsupported. The compactness of the weighted space Σ_s(R^N) is likely true, but a correct proof must be supplied, and the statement must be modified to handle the case N≤2s.","section":"Section 3, Lemma 3.1"},{"comment":"The paper claims existence of ground states for the range 2s/N ≤ σ < 2s/(N-2s) via the constrained problem (62)-(63), and states 'Similar to Lemma 3.2 and Theorem 2.1, there exists a local minimizer'. This is not a proof; no rigorous argument is given that the auxiliary problem has a minimizer or that it yields a solution to (65). If the paper intends to claim existence in this range, a complete proof is needed. As it stands, those results are numerical evidence only and should be labeled as such.","section":"Section 6, Eqs. (62)-(63)"}],"minor_comments":[{"comment":"The exponent θ in the fractional Gagliardo-Nirenberg inequality appears to be misprinted as θ = Ns/(2s(σ+1)), which simplifies to N/(2(σ+1)). The subsequent definitions p = 1/(θ(1+σ)) and q = 1/(1-θ(1+σ)) are consistent only if θ = Nσ/(2s(σ+1)). Please correct this typo.","section":"Section 3, Lemma 3.2, Eq. (15)"},{"comment":"The orbital-stability theorem is conditional on an assumed uniqueness and well-posedness result for the Cauchy problem (1). The paper is explicit about this, but the abstract and introduction should state clearly that stability is proved only under this assumption, and it would be helpful to cite or prove local well-posedness in Σ_s(R^N).","section":"Theorems 2.2 and 4.1"},{"comment":"The Lagrange multiplier argument is presented formally. Since Theorem 2.1 only asserts existence of a minimizer of (8), this is not a flaw for the theorem as stated, but the paper repeatedly calls these minimizers 'ground state solutions'; a precise statement connecting minimizers to solutions of (6) would strengthen the paper.","section":"Section 2, Eqs. (11)-(13)"},{"comment":"There are numerous typographical and notational errors: in Lemma 3.2 the integral set is written as Ω but should be R^N; 'Figure 7.2' in Section 7.2 should presumably be 'Figure 10'; the caption of Figure 10 says 'δ = 1' but should be 'σ = 1'; and 'Abosolute value' should be 'Absolute value'. The manuscript would benefit from a careful proofreading pass.","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":"The core variational proof is likely repairable, and the numerical work is substantial, but the manuscript as written contains a genuine gap in the compactness argument of Lemma 3.1 and an unproved existence claim in Section 6. I recommend major revision rather than rejection. The paper cites several works by the second author; these citations appear to be relevant and standard, and I did not see a circularity problem."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nShort version: this paper fills a natural gap—fractional NLS with a harmonic trap—and the main results are probably true, but the proof of Theorem 2.1 has a real gap in its compactness lemma, and the critical-range existence is only asserted. I'd accept it for peer review, but revision is required.\n\nWhat's genuinely new is the combination of [9] (s=1 with trap) and [15] (fractional without trap). The weighted space Σ_s is the right setting, and the existence/stability structure follows the standard variational pattern. The numerical work is the most original part: it explores ground states for different s, non-radial potentials, the critical case, and dynamics, and the figures show the expected behavior.\n\nThe stress-test concern about Lemma 3.1 lands. The proof says H^s(R^N) compacts into L^p(R^N) for 2<p<2N/(N−2s); that is false on unbounded domains because translations break compactness. The lemma's statement (compactness of the weighted Σ_s) is true, but the supplied proof doesn't establish it. The p=2 argument can be fixed, the p>2 case is missing, and for s≥N/2 the exponent 2N/(N−2s) is meaningless. As written, the strong convergence of a minimizing subsequence in L^{2σ+2} is unsupported. This is repairable via standard compactness/interpolation arguments for weighted spaces, but it's a genuine gap, not a typo.\n\nSecond, Section 6 introduces the constrained problem (62)-(63) and claims existence of a local minimizer 'similar to Lemma 3.2 and Theorem 2.1,' without proof. That's a separate unproven assertion, especially since the functional K and the projection step differ from the subcritical case.\n\nThird, the numerics have no code, data, or error analysis. The stability measure D(s,t) in Figure 9(b) isn't obviously connected to the orbital stability norm, so the conclusion that smaller s gives worse stability is informal.\n\nOn citations: the paper leans on several co-author papers for the variational framework and symmetrization inequalities. That's legitimate; there's no circularity. The stability theorem is conditional on uniqueness and conservation, which is standard in this literature but not proved here.\n\nThe core message: the main ideas are sound and the theorems are likely correct, but the write-up is not. A serious referee should ask for a corrected Lemma 3.1 and either a full proof or a clearly labeled conjecture for the critical range. I wouldn't cite it in its current form; after careful repair I might. It's a good reading-group example of a false Sobolev embedding claim slipping into a proof.\n\nI'd send it to review, with the expectation of major revision.","headline":"A plausible extension of ground-state theory to fractional NLS with a harmonic trap, but the proof of the main existence theorem rests on a faulty compactness lemma.","tokens_in":16073,"tokens_out":6265,"would_cite":false,"duration_ms":56464,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35Q55","35J60","47J30"],"pacs":[],"model":"deepseek-v4-flash","headline":"For 0<σ<2s/N, the fractional Schrödinger equation with a harmonic trap has nonnegative, radial, decreasing ground-state standing waves, orbitally stable when the Cauchy problem is unique.","keywords":["fractional Schrödinger equation","fractional Laplacian","harmonic potential","standing waves","ground state solutions","orbital stability","constrained minimization","normalized gradient flow"],"falsifier":"Translate a fixed compactly supported bump: $u_n(x)=\\varphi(x-ne_1)$ is bounded in $H^s(\\mathbb{R}^N)$ but has no convergent subsequence in $L^p(\\mathbb{R}^N)$, so the proof's appeal to compactness of $H^s$ on the whole space is demonstrably false; a correct proof must use the $\\|xu\\|_{L^2}$ weight, for example by showing that $\\|u_n\\|_{\\Sigma_s}\\to\\infty$ for such translations.","tokens_in":15086,"feed_emoji":"⚛️","tokens_out":9643,"duration_ms":84746,"temperature":0.7,"pith_summary":"The paper establishes that the fractional nonlinear Schrödinger equation with a harmonic trapping potential, $i\\psi_t=(-\\Delta)^s\\psi+|x|^2\\psi-|\\psi|^{2\\sigma}\\psi$, has standing-wave solutions when the nonlinearity is subcritical, $0<\\sigma<2s/N$. The proof works by minimizing the associated energy on the sphere of fixed $L^2$ mass within the weighted space $\\Sigma_s(\\mathbb{R}^N)$, where the potential restores compactness lost on the whole space. The resulting minimizer is nonnegative, radial, and radially decreasing, and it gives a ground-state solution through a Lagrange multiplier $\\lambda$. Under an additional assumption that the initial-value problem has unique solutions conserving mass and energy, these standing waves are orbitally stable. The paper also supplies numerical evidence for ground states, stability, and dynamics, including the critical and supercritical ranges via a second constrained problem.","feed_headline":"Fractional Schrödinger ground states exist in a harmonic trap","feed_subtitle":"Subcritical nonlinearities yield radial, decreasing minimizers of the constrained energy; stability follows conditionally.","key_machinery":"The central object is the weighted Sobolev space $\\Sigma_s(\\mathbb{R}^N)=\\{u\\in H^s(\\mathbb{R}^N): \\|u\\|_{L^2}+\\|\\nabla^s u\\|_{L^2}+\\|xu\\|_{L^2}<\\infty\\}$, whose norm includes the harmonic potential. Its compact embedding into $L^p$ for $2\\le p<2N/(N-2s)$ is what lets a minimizing sequence for $I_c$ converge strongly, overcoming the lack of compactness on $\\mathbb{R}^N$. The fractional Gagliardo\\textendash Nirenberg inequality supplies the $L^{2\\sigma+2}$ bound that keeps the energy finite below the threshold $\\sigma<2s/N$. Schwarz symmetrization then converts any minimizer into a radial, radially decreasing one without raising the energy. For stability, the machinery is a standard compactness-based variational argument: the set of minimizers at fixed mass is shown to be stable under the flow whenever the Cauchy problem is well posed with conserved mass and energy.","core_discovery":"The central claim is Theorem 2.1: for $0<\\sigma<2s/N$, the constrained minimization problem $I_c=\\inf\\{J(u): u\\in\\Sigma_s(\\mathbb{R}^N),\\ \\|u\\|_{L^2}=c\\}$ admits a minimizer, and that minimizer is a nonnegative, radial, radially decreasing ground state solution of $(-\\Delta)^s u+|x|^2u-|u|^{2\\sigma}u=\\lambda u$. The mechanism is that the harmonic potential makes the embedding $\\Sigma_s(\\mathbb{R}^N)\\hookrightarrow L^p(\\mathbb{R}^N)$ compact for $2\\le p<2N/(N-2s)$, so minimizing sequences converge; the fractional Gagliardo\\textendash Nirenberg inequality bounds the energy below precisely when $\\sigma<2s/N$. Theorem 2.2 then asserts orbital stability of the set of minimizers, conditional on uniqueness and on conservation of mass and energy for the Cauchy problem. Numerical sections complement the existence result by solving the constrained problems with normalized gradient flow and split-step Fourier methods, and by probing the critical range $2s/N\\le\\sigma<2s/(N-2s)$ through an alternative minimization with fixed $L^{2\\sigma+2}$ norm.","pith_inferences":["An editorial repair of Lemma 3.1: a correct compactness proof should split $\\mathbb{R}^N$ into a ball and its complement, use the tail bound $\\|u\\|_{L^2(|x|>R)}\\le R^{-1}\\|xu\\|_{L^2}$, then take a diagonal subsequence; the false step in the paper is a proof gap, not evidence against the theorem.","The convergence $\\lambda_c\\to\\lambda_0$ as $c\\to0$ seen numerically is consistent with $\\lambda_0$ being the lowest eigenvalue of the linear operator $(-\\Delta)^s+|x|^2$; proving this would connect the nonlinear ground states to the linear spectrum.","The non-radial numerical ground states for asymmetric potentials suggest the radial symmetry theorem depends essentially on the potential being radial; extending existence to general trapping potentials would require new rearrangement tools."],"forward_implications":["For every mass $c>0$ and exponent $0<\\sigma<2s/N$, the constrained problem (8) has a minimizer, so the elliptic equation (6) has a nonnegative, radial, radially decreasing ground-state solution for some Lagrange multiplier $\\lambda$.","Because the minimizer set is compact in $\\Sigma_s(\\mathbb{R}^N)$, the standing waves $e^{-i\\lambda t}u(x)$ form an orbitally stable family whenever the Cauchy problem is well posed with conserved mass and energy.","The threshold $\\sigma=2s/N$ is sharp for the original minimization: above it the energy is unbounded below, and the paper's alternative fixed-$L^{2\\sigma+2}$ formulation gives standing waves numerically up to the $L^2$-supercritical threshold $\\sigma<2s/(N-2s)$.","Numerical experiments indicate that as $s$ decreases toward $\\sigma N/2$, ground states become more peaked, the energy diverges to $-\\infty$ at fixed mass, and orbital stability degrades."],"supporting_citations":[{"why":"Supplies the fractional Gagliardo\\textendash Nirenberg inequality used to bound the $L^{2\\sigma+2}$ term and prove the energy is bounded below when $\\sigma<2s/N$.","marker":"[12]"},{"why":"Provides the variational approach to standing waves, the symmetrization inequality for the potential term, and the $\\lambda_c$ convergence comparison.","marker":"[9]"},{"why":"Supplies the Schwarz symmetrization facts that preserve $L^p$ norms and lower the fractional gradient, converting minimizers into radial decreasing ones.","marker":"[14]"},{"why":"Establishes the $C^1$ regularity and variational stability framework used in Lemma 4.1 for the energy functionals.","marker":"[3]"},{"why":"Gives the compactness-based stability argument for standing waves used in Theorem 4.1.","marker":"[13]"},{"why":"Defines the notion of orbital stability of the minimizing set used in Definition 4.1 and Theorem 4.1.","marker":"[4]"},{"why":"Supplies the normalized gradient flow method used numerically to compute ground states.","marker":"[2]"},{"why":"Supplies the split-step, mass-conservative Fourier spectral method used to simulate the time dynamics.","marker":"[7]"},{"why":"Provides the numerical baseline for fractional NLS without potential; the paper compares dynamics, convergence, and blow-up phenomena against it.","marker":"[15]"}],"fun_headline_variants":["Subcritical nonlinearity yields ground states in fractional trap","Harmonic trap ensures compactness for fractional ground states","Radial decreasing ground states for fractional Schrodinger","Conditional stability of standing waves in fractional trap"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument's load-bearing premise is Lemma 3.1's claim that the weighted space $\\Sigma_s(\\mathbb{R}^N)$ sits compactly inside the Lebesgue spaces $L^p(\\mathbb{R}^N)$; as written, the proof of that lemma appeals to a compactness of $H^s$ on the whole space that is false, so the extraction of a convergent minimizing sequence is not justified until that point is repaired.","fun_headline_variants_meta":{"raw":{"variants":["Subcritical nonlinearity yields ground states in fractional trap","Harmonic trap ensures compactness for fractional ground states","Radial decreasing ground states for fractional Schrodinger","Conditional stability of standing waves in fractional trap"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000867,"raw_usage":{"total_tokens":3720,"prompt_tokens":874,"completion_tokens":2846,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":490,"completion_tokens_details":{"reasoning_tokens":2785}},"tokens_in":490,"tokens_out":2846,"duration_ms":21591,"temperature":1.0,"reasoning_tokens":2785,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:05:50.049246+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Translate a fixed compactly supported bump: $u_n(x)=\\varphi(x-ne_1)$ is bounded in $H^s(\\mathbb{R}^N)$ but has no convergent subsequence in $L^p(\\mathbb{R}^N)$, so the proof's appeal to compactness of $H^s$ on the whole space is demonstrably false; a correct proof must use the $\\|xu\\|_{L^2}$ weight, for example by showing that $\\|u_n\\|_{\\Sigma_s}\\to\\infty$ for such translations.","supporting_citations":[{"cited_title":"Hajaiej, L","cited_arxiv_id":null,"evidence_quote":"Supplies the fractional Gagliardo\\textendash Nirenberg inequality used to bound the $L^{2\\sigma+2}$ term and prove the energy is bounded below when $\\sigma<2s/N$."},{"cited_title":"Hadj Selem, H","cited_arxiv_id":null,"evidence_quote":"Provides the variational approach to standing waves, the symmetrization inequality for the potential term, and the $\\lambda_c$ convergence comparison."},{"cited_title":"Hajaiej and C","cited_arxiv_id":null,"evidence_quote":"Supplies the Schwarz symmetrization facts that preserve $L^p$ norms and lower the fractional gradient, converting minimizers into radial decreasing ones."},{"cited_title":"Carles and H","cited_arxiv_id":null,"evidence_quote":"Establishes the $C^1$ regularity and variational stability framework used in Lemma 4.1 for the energy functionals."},{"cited_title":"Hajaiej and C","cited_arxiv_id":null,"evidence_quote":"Gives the compactness-based stability argument for standing waves used in Theorem 4.1."},{"cited_title":"Cazenave and P.-L","cited_arxiv_id":null,"evidence_quote":"Defines the notion of orbital stability of the minimizing set used in Definition 4.1 and Theorem 4.1."},{"cited_title":"Bao and Q","cited_arxiv_id":null,"evidence_quote":"Supplies the normalized gradient flow method used numerically to compute ground states."},{"cited_title":"Duo and Y","cited_arxiv_id":null,"evidence_quote":"Supplies the split-step, mass-conservative Fourier spectral method used to simulate the time dynamics."},{"cited_title":"Klein, C","cited_arxiv_id":null,"evidence_quote":"Provides the numerical baseline for fractional NLS without potential; the paper compares dynamics, convergence, and blow-up phenomena against it."}],"review_version":1}