{"id":"7deb290e-08ce-43b8-8b26-717fff2a5513","arxiv_id":"2411.18532","paper_version":1,"verdict":"REJECT","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"high","formal_verification":"none","parameter_count":0,"one_line_summary":"An L2-norm-preserving nonlocal parabolic flow is shown to have global solutions that converge to stationary states, with positive data on a ball converging to the ground state.","lead":"This paper proves that a nonlocal parabolic flow which preserves the L2 norm admits global energy-space solutions on bounded domains and Euclidean space, and that positive solutions on a ball converge to the unique ground state. The global well-posedness result stands in contrast to a closely related norm-preserving flow studied earlier, which can blow up for large nonlinearities.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Lemma 4.1(44) is false: a Dirichlet eigenfunction with eigenvalue -ω is stationary with μ=0, so A≡0 and the strict negativity asserted in (44) never occurs; this invalidates the ω<0 global well-posedness proof in Proposition 4.3 and leaves Theorem 1.2 unsupported for ω<0.","rationale":"Reading in good faith, the paper contains substantial correct work: local well-posedness via contraction and Schauder arguments, global well-posedness for ω≥0, and asymptotic convergence to stationary states and to the ground state for ω>0. However, the abstract and Theorem 1.2 advertise global well-posedness for all ω∈R, and the only proof for ω<0 is Proposition 4.3, which depends on Lemma 4.1(44). The Dirichlet-eigenfunction example is a genuine counterexample to (44), not a missing edge case: at the instant A=0 with ∂_t u=0, the asserted one-sided negativity fails identically for all time. This invalidates the claimed invariance of G and H in Proposition 4.3 and leaves the uniform H^1 bound for data in H unjustified. Since Theorem 1.2's full statement requires exactly this bound, the central claim is unsupported as written. The flaw is localized and plausibly repairable by a first-hitting-time argument that treats A=0 separately and by handling stationary eigenfunctions, which is why the paper's other contributions remain credible; nevertheless, the manuscript in its current form does not prove its advertised main theorem.","tokens_in":17809,"tokens_out":11442,"duration_ms":98410,"concrete_test":"Take Ω to be the unit ball (or any bounded C^2 domain), fix ω<0, and choose u0=φ_k, a Dirichlet eigenfunction satisfying -Δφ_k=λ_k φ_k with λ_k=-ω. Compute μ[u0]=(λ_k||φ_k||^2_2+ω||φ_k||^2_2)/||φ_k||^{2σ+2}_{2σ+2}=0, so u(t)=φ_k solves (1) with ∂_t u=0. Then A(t)=||∇φ_k||^2+ω||φ_k||^2=0 for every t≥0, directly contradicting Lemma 4.1(44). This explicit computation settles the concern without further analysis.","verdict_should_be":"REJECT","load_bearing_attack":"Theorem 1.2 claims global well-posedness for every ω∈R. For ω<0 the proof is Proposition 4.3, which splits H^1_0 into G={A≤0} and H={A>0}, with A(v)=||∇v||^2+ω||u0||^2, and asserts these sets are invariant. The invariance argument relies on Lemma 4.1(44), which says that if A(t0)=0 then A(t0+ε)<0 for some ε>0. This claim is false. Let φ be a nonzero Dirichlet eigenfunction of -Δ satisfying -Δφ=-ω φ, which is possible because ω<0. Take u0=φ. Then ||∇u0||^2 = -ω||u0||^2, so A(0)=0 and μ[u0]=0. The function u(t)≡φ is a solution of (1), so A(t)≡0 for all t≥0, and no ε with A(t0+ε)<0 exists. Thus (44) fails. Consequently, Proposition 4.3's invariance of G and H is not established, and its uniform H^1 bound for data in H—which uses (45) to keep F positive—has no valid proof. Since this is the only mechanism offered for the ω<0 case, Theorem 1.2 lacks a proof for that parameter range. The ω≥0 global result and the ω>0 asymptotic theorems appear unaffected.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the nonlocal parabolic equation (1), in which the nonlocal coefficient mu[u] is chosen so that the L^2 norm of solutions is preserved in time. It claims local well-posedness in H^1_0 on bounded C^2 domains for every subcritical sigma, and on R^d under additional restrictions; it then claims global well-posedness with a uniform H^1 bound for every omega in R (Theorem 1.2). The asymptotic results assert that for omega>0 the omega-limit set consists of stationary states (Theorem 1.3), and that on a ball with positive initial data the solution converges in H^1 to the unique positive ground state (Theorem 1.4). The central mechanism for the omega<0 global bound is Proposition 4.3, which relies on the invariance of the sets G and H established through Lemma 4.1.","tokens_in":18045,"tokens_out":11769,"duration_ms":102556,"significance":"The omega>0 results are interesting: if they are correct, they provide a sharp contrast with the grow-up behavior known for the related norm-preserving model (3), and they give a relatively clean Lyapunov-functional route to convergence to stationary states. The paper also contains useful technical work in the local well-posedness arguments, especially the Schauder fixed-point proof on bounded domains and the space-time estimates on R^d. However, the claimed global well-posedness for every omega in R, which is one of the headline results, is not established: the proof for omega<0 rests on a false statement in Lemma 4.1. The omega>=0 parts and the asymptotic analysis appear to be unaffected by this defect, but the manuscript in its present form does not prove one of its central theorems.","major_comments":[{"comment":"The assertion that A(t0)=0 implies A(t0+epsilon)<0 for some epsilon>0 is false. Let omega<0 and let phi be a nonzero Dirichlet eigenfunction of -Delta on H^1_0(Omega) with eigenvalue -omega. Then mu[phi]=0 and u(t)≡phi solves (1); moreover A(t)=||grad phi||^2+omega||phi||^2=0 for every t>=0. Hence no t0>0 admits an epsilon with A(t0+epsilon)<0. The derivation in the proof only gives A'(t0)=-2||partial_t u(t0)||^2 <=0, which is compatible with A remaining identically zero in the stationary case; the stationary case is not excluded by the hypotheses of the lemma.","section":"Lemma 4.1, eq. (44)"},{"comment":"The invariance of G and H for omega<0, and consequently the uniform H^1 bound for data in H, rests entirely on the false statement (44). The claim that (45) implies F[u0] and F[u(t)] have the same sign for all t is unjustified, because (45) is derived only while the denominator A(s)=||grad u(s)||^2+omega||u0||^2 does not vanish; once A hits zero, the representation (45) breaks down. Since Proposition 4.3 is the only argument given for global well-posedness when omega<0, Theorem 1.2 is not proved for that parameter range. The omega>=0 global result in Corollary 4.2 is not affected.","section":"Proposition 4.3"}],"minor_comments":[{"comment":"The sentence 'sup_{t>=0} ||grad u(t)||_{L^2} <= M < 0' is nonsensical as written; M must be a positive finite constant, so it should read M < infinity.","section":"Lemma 5.1, proof"},{"comment":"The text says 'by (28), (26) and (26)' but the second reference should be (27), since (26) alone gives the self-map property while (27) gives the contraction estimate.","section":"Proposition 3.1, after (27)"},{"comment":"The sentence 'Thus, u satisfies (37)' should refer to equation (1), since (37) is the regularized equation with mu_epsilon and the limit has already been taken.","section":"Proposition 3.8, final sentence"},{"comment":"The set S0 is used without being defined; presumably it denotes the set of stationary states of (7) with the fixed L^2 norm, but the notation should be introduced.","section":"Corollary 5.3"},{"comment":"The notation Q is used inconsistently: the proof repeatedly writes 'u(t_k) -> Q' and 'E[u]=E[Q]' when the intended limit is the ground state Qgs; this should be corrected for clarity.","section":"Theorem 5.4, proof"}],"recommendation":"reject","confidential_remarks":"The defect in Lemma 4.1 is elementary and strikes at one of the main theorems: the omega<0 part of Theorem 1.2. I do not see how to repair the present proof strategy without a substantially different argument for the omega<0 case. If the authors restrict the global well-posedness claim to omega>=0, the remaining results may be publishable after a careful revision, but the manuscript in its current form does not support the advertised range omega in R."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe paper is a mix of solid work and one load-bearing gap. The headline: the abstract advertises global well-posedness for all ω∈R, but the proof for ω<0 rests on Lemma 4.1(44), which is false. A nonzero Dirichlet eigenfunction φ with -Δφ = -ωφ (ω<0) is a stationary solution: μ[φ]=0, so A(t)≡0 for all t, contradicting the claimed strict negativity in (44). That knocks out Proposition 4.3's invariance argument and leaves Theorem 1.2 unsupported for ω<0. The ω≥0 results are unaffected, and the asymptotics for ω>0 stand on their own.\n\nWhat's genuinely new: local well-posedness on bounded domains for the full subcritical range via a Schauder fixed point argument (the earlier Ma–Cheng paper only treated ω=0 on compact manifolds), the treatment of R^d with two separate contraction arguments, and the ground-state convergence on balls for positive data. The Lyapunov structure (Lemma 4.1 up to the false claim, Corollary 4.2) is clean, and the asymptotic analysis in Section 5 is mostly careful. The paper is not a model of exposition—there are typos like 'M < 0' in Lemma 5.1—but the mathematics is mostly well-organized.\n\nThe central question is whether the ω<0 case can be fixed. The stationary eigenfunction sits exactly on the boundary between G and H, so the decomposition into G and H is not invariant in the claimed way. One could try to split trajectories at first hitting time of A=0, but then the H^1 bound for data in H may fail if the solution actually crosses; you'd need a different mechanism. The author should either repair this or restrict Theorems 1.1–1.2 to ω≥0. As stated, the paper overclaims.\n\nWho is this for? People working on norm-preserving parabolic flows and ground-state computation will find the ω≥0 part useful, and the flaw itself is instructive. I'd send it to a referee, but I'd expect a major revision. My recommendation: engage with it, but do not accept the current version.","headline":"Solid ω≥0 results and a clean Lyapunov structure, but the advertised all-ω global well-posedness rests on a false lemma.","tokens_in":18667,"tokens_out":4220,"would_cite":false,"duration_ms":33828,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35K55","35B40","35B45"],"pacs":[],"model":"deepseek-v4-flash","headline":"A nonlocal heat flow that keeps its $L^2$ mass fixed is globally well-posed for every real frequency and subcritical nonlinearity, and on a ball with positive data it converges strongly to the unique ground state.","keywords":["nonlocal parabolic flow","norm-preserving flow","L2-constraint","global well-posedness","Lyapunov functional","ground state","asymptotic convergence","subcritical nonlinearity"],"falsifier":"Take $\\omega<0$ and $u_0=\\varphi$, where $\\varphi$ is a Dirichlet eigenfunction of $-\\Delta$ with eigenvalue $-\\omega$ normalized in $L^2$: the exact solution is $u(t)=\\varphi$, so $\\|\\nabla u(t)\\|^2+\\omega\\|u_0\\|^2\\equiv 0$, contradicting the strict negativity asserted in (44). Initializing instead just inside the set $G$ and numerically integrating (1) would test whether the claimed invariance of $G$, and hence the $\\omega<0$ uniform bound, actually holds.","tokens_in":17476,"feed_emoji":"🔥","tokens_out":16850,"duration_ms":139721,"temperature":0.7,"pith_summary":"This paper studies a heat-type equation modified by a nonlocal term chosen so that every solution keeps its $L^2$ mass exactly constant over time. The central claim is that, for every real frequency parameter $\\omega$ and every subcritical nonlinearity exponent, the flow is global in time: no solution blows up, and the $H^1$ norm (size together with gradient size) stays uniformly bounded, on bounded domains and on the whole space under the stated restrictions. The paper then proves that on bounded domains the asymptotic dynamics are constrained to stationary states, and that on a ball with positive initial data and $\\omega>0$ the entire trajectory converges strongly in $H^1$ to the unique positive ground state. This matters because norm-preserving flows are widely used as numerical tools to compute stationary states, and the closely related constrained gradient flow was previously known to admit growing-up solutions in the intercritical regime.","feed_headline":"A heat flow that pins its mass never blows up","feed_subtitle":"An energy-decreasing quantity keeps solutions bounded and drives them to steady states—on a ball, to the ground state.","key_machinery":"The load-bearing object is the Lyapunov functional $F[u]=(\\|\\nabla u\\|_{L^2}^2+\\omega\\|u\\|_{L^2}^2)/\\|u\\|_{L^{2\\sigma+2}}^{2\\sigma+2}$, together with the nonlocal multiplier $\\mu[u]$ that enforces the fixed-$L^2$ constraint. $F$ decreases along the flow and its explicit decay formula transfers the integrated decay of $\\|\\partial_t u\\|_{L^2}$ into a uniform bound on $\\|\\nabla u\\|_{L^2}$, because Gagliardo-Nirenberg interpolation bounds the denominator $\\|u\\|_{L^{2\\sigma+2}}$ by a power of the gradient. The same functional selects the ground state as its minimizer, which is what turns subsequential convergence into full convergence on a ball.","core_discovery":"The discovery is that the nonlocal multiplier $\\mu[u]=(\\|\\nabla u\\|_{L^2}^2+\\omega\\|u\\|_{L^2}^2)/\\|u\\|_{L^{2\\sigma+2}}^{2\\sigma+2}$, which enforces $\\|u(t)\\|_{L^2}=\\|u_0\\|_{L^2}$ for all $t$, makes the flow gradient-like rather than destabilizing. The quantity $F[u(t)]=(\\|\\nabla u(t)\\|^2+\\omega\\|u_0\\|^2)/\\|u(t)\\|_{L^{2\\sigma+2}}^{2\\sigma+2}$ is non-increasing and satisfies the exact differential identity $\\frac{d}{dt}\\log \\frac{\\sqrt{\\|\\nabla u\\|^2+\\omega\\|u_0\\|^2}}{\\|u\\|_{L^{2\\sigma+2}}}=-\\frac{\\|\\partial_t u\\|_{L^2}^2}{\\|\\nabla u\\|^2+\\omega\\|u_0\\|^2}$. Combined with the Gagliardo-Nirenberg inequality, this identity controls the gradient uniformly for every subcritical $\\sigma$, giving global existence; the same Lyapunov function, together with compactness and uniqueness of the positive stationary solution, gives strong $H^1$ convergence to the ground state on a ball.","pith_inferences":["The paper leaves open whether the $\\omega$-limit set is a singleton outside the ball/positive-data case; a natural next step would be to combine the Lyapunov identity with spectral properties of the linearized operator around the ground state to obtain full convergence and rates.","The same $L^2$-sphere constraint and the same $F$-type functional underlie normalized gradient flows used to compute ground states of Bose-Einstein condensates; the global-bound result suggests these schemes remain stable across the full subcritical range, beyond the intercritical regime where the unmodified flow grows up.","The paper notes that its Aubin-Lions/Schauder route fails on $\\mathbb{R}^d$, leaving a gap for the whole space; a plausible repair is to adapt the density argument with localized compactness, which would remove the extra restrictions on $\\sigma$ currently needed there."],"forward_implications":["Global existence and a uniform $H^1$ bound hold for all $\\omega\\in\\mathbb{R}$ and all subcritical $\\sigma<2/(d-2)_+$ on bounded $C^2$ domains, without any smallness assumption on the initial data.","For $\\omega>0$ on a bounded domain, the $\\omega$-limit set is nonempty, compact, connected, and contained in the set of stationary states with the same $L^2$ norm; the same conclusion holds on $\\mathbb{R}^d$ for radially symmetric data.","When the domain is a ball, $\\omega>0$, and the initial datum is positive, the whole trajectory converges in $H^1$ to the unique positive stationary state $Q_{gs}$, not merely along a subsequence.","The Lyapunov identity implies $\\int_0^\\infty\\|\\partial_t u(t)\\|_{L^2}^2\\,dt<\\infty$, so the time derivative of the solution decays to zero in $L^2$ as $t\\to\\infty$.","In contrast to the earlier constrained gradient flow (3), which exhibits growing-up solutions for $2/d\\le\\sigma<2/(d-2)_+$, the present flow is asserted to remain uniformly bounded for this whole range."],"supporting_citations":[{"why":"Provides the prior constrained gradient flow that admits growing-up solutions, the contrast motivating the global-bound result.","marker":"[1]"},{"why":"Identifies ground states as minimizers of the Lyapunov functional, supplying the target state in the convergence theorem.","marker":"[6]"},{"why":"Supplies the Aubin-Lions compactness lemma used in the Schauder fixed-point argument for bounded domains.","marker":"[8]"},{"why":"Supplies the Gagliardo-Nirenberg interpolation inequality used to turn Lyapunov decay into uniform gradient bounds.","marker":"[12]"},{"why":"Gives symmetry and uniqueness properties of the positive stationary solution used as the ground state.","marker":"[17]"},{"why":"Proves uniqueness of positive solutions in the whole space, extending ground-state identification to $\\mathbb{R}^d$.","marker":"[20]"},{"why":"Studies the same equation with $\\omega=0$ on compact manifolds, the prior setting generalized here.","marker":"[23]"},{"why":"Provides the heat-semigroup smoothing estimates and parabolic well-posedness tools used in the local existence proofs.","marker":"[24]"}],"fun_headline_variants":["Norm-fixed heat flow never blows up, hits ground state","L2-conserving flow converges to ground state on balls","Nonlocal pinning yields global solutions and steady states","Mass-pinned parabolic flow: no blow-up, strong convergence"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is Lemma 4.1's claim that a solution reaching the surface $\\|\\nabla u\\|^2+\\omega\\|u_0\\|^2=0$ must immediately move below it; a stationary solution that is a Dirichlet eigenfunction of $-\\Delta$ with eigenvalue $-\\omega$ stays on that surface for all time, so this premise fails exactly there.","fun_headline_variants_meta":{"raw":{"variants":["Norm-fixed heat flow never blows up, hits ground state","L2-conserving flow converges to ground state on balls","Nonlocal pinning yields global solutions and steady states","Mass-pinned parabolic flow: no blow-up, strong convergence"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000709,"raw_usage":{"total_tokens":3159,"prompt_tokens":877,"completion_tokens":2282,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":493,"completion_tokens_details":{"reasoning_tokens":2213}},"tokens_in":493,"tokens_out":2282,"duration_ms":18163,"temperature":1.0,"reasoning_tokens":2213,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T11:08:21.260658+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take $\\omega<0$ and $u_0=\\varphi$, where $\\varphi$ is a Dirichlet eigenfunction of $-\\Delta$ with eigenvalue $-\\omega$ normalized in $L^2$: the exact solution is $u(t)=\\varphi$, so $\\|\\nabla u(t)\\|^2+\\omega\\|u_0\\|^2\\equiv 0$, contradicting the strict negativity asserted in (44). Initializing instead just inside the set $G$ and numerically integrating (1) would test whether the claimed invariance of $G$, and hence the $\\omega<0$ uniform bound, actually holds.","supporting_citations":[{"cited_title":"Antonelli, P","cited_arxiv_id":null,"evidence_quote":"Provides the prior constrained gradient flow that admits growing-up solutions, the contrast motivating the global-bound result."},{"cited_title":"Berestycki and J","cited_arxiv_id":null,"evidence_quote":"Identifies ground states as minimizers of the Lyapunov functional, supplying the target state in the convergence theorem."},{"cited_title":"Boyer and P","cited_arxiv_id":null,"evidence_quote":"Supplies the Aubin-Lions compactness lemma used in the Schauder fixed-point argument for bounded domains."},{"cited_title":"Cazenave","cited_arxiv_id":null,"evidence_quote":"Supplies the Gagliardo-Nirenberg interpolation inequality used to turn Lyapunov decay into uniform gradient bounds."},{"cited_title":"Gidas, W","cited_arxiv_id":null,"evidence_quote":"Gives symmetry and uniqueness properties of the positive stationary solution used as the ground state."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Proves uniqueness of positive solutions in the whole space, extending ground-state identification to $\\mathbb{R}^d$."},{"cited_title":"Ma and L","cited_arxiv_id":null,"evidence_quote":"Studies the same equation with $\\omega=0$ on compact manifolds, the prior setting generalized here."},{"cited_title":"Quittner and P","cited_arxiv_id":null,"evidence_quote":"Provides the heat-semigroup smoothing estimates and parabolic well-posedness tools used in the local existence proofs."}],"review_version":1}