{"id":"596956d4-1780-4091-bdb8-8d28f7974890","arxiv_id":"2607.12579","paper_version":2,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Under Coulomb MMD energy, the Wasserstein gradient flow converges exponentially to uniformly positive targets under PL-type coercivity, exhibits rigidity of critical points, and admits no uniform whole-space convergence rate.","lead":"This paper proves global existence, instantaneous regularization, and long-time convergence results for Wasserstein gradient flows of the squared Coulomb Maximum Mean Discrepancy, including new Polyak–Łojasiewicz inequalities and explicit rate obstructions. It also classifies critical points and proves planar convergence, giving a rigorous basis for kernel-based generative modeling and sampling dynamics.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"No significant objection identified: the acknowledged non-uniqueness for measure initial data limits interpretation but does not invalidate the stated convergence theorems.","rationale":"The paper is careful with its claims: long-time theorems are stated for constructed solutions, and the non-uniqueness is disclosed before §2.1. I traced the key estimates: the ultracontractive L^p bound (Prop. 2.1) is sound; the compactness argument for existence uses standard Aubin–Lions and near–far field bounds; the defective-PL proof of Thm 1.7 has correct variation-of-constants and the defect decay Lemma 3.4 is valid; the variational barrier proof of Thm 1.3, though intricate, has a consistent Euler equation and Hessian-trace contradiction; the radial PL proof's shellwise inequality is correct; the planar convergence proof's tightness/LaSalle argument is standard and the application of the rigidity theorem is legitimate because (ρ̄−μ)^+ ≤ ρ̄. The only substantive caveat is the selection issue, which the reader also identified. I therefore do not see a reason to move the verdict.","tokens_in":44043,"tokens_out":41568,"duration_ms":378574,"concrete_test":"To test the selection issue: on T^2 with µ≡1 and ρ0=δ0, compute the smooth-approximation limit of (1.1) with two different mollifier sequences (e.g., Gaussian vs. box, width ε→0) and compare ρ_t at a fixed positive time in the H^{-1} metric; also compare against a JKO/minimizing-movement limit if feasible. If the limiting trajectories differ at positive times, the constructed-solution caveat is substantive; if they agree numerically, it supports (but does not prove) a canonical semigroup.","verdict_should_be":"UNCHANGED","load_bearing_attack":"I found no internal inconsistency or missing step that would break the central claims. The main soft spot is exactly the one the paper concedes: Theorem 1.2 constructs solutions from arbitrary Borel initial measures by smooth approximation, and uniqueness is only proved in the L∞ class; the paper states 'We do not know whether uniqueness holds' (Section 2, before §2.1). Consequently Theorems 1.7 and 1.22 describe the behavior of the particular constructed (selected) solutions, not an intrinsic semigroup for (1.1). This is a genuine scope limitation: if different regularizations or JKO-type schemes select different limits, the equation has no canonical long-time behavior from singular data. It is not, however, a correctness flaw in the stated results, because those results are explicitly quantified over 'any solution constructed in Theorem 1.2,' and the convergence proofs use only properties (ultracontractivity, energy inequality, positive-time boundedness) shared by every such constructed solution. The variational-barrier proof of Theorem 1.3 and the use of the external level-set theorem [APR20, Thm 1.1] in Theorem 1.21 are intricate, but I could not locate a concrete gap; the first-variation and Euler-equation steps are internally consistent, modulo a harmless typo in the display of (3.34), which the subsequent argument uses correctly.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the long-time behavior of the Wasserstein gradient flow of the squared Coulomb Maximum Mean Discrepancy, i.e. the nonlocal transport equation ∂tρ − div(ρ∇g∗(ρ−µ)) = 0 on the flat torus and on R^d. It constructs global weak solutions from arbitrary Borel initial measures when the target µ is bounded, with instantaneous L∞ regularization and explicit Lp estimates, proves uniqueness in the bounded initial-data class by an Osgood/modulated-energy argument, and shows that Hölder seminorms can grow exponentially. On T^d it proves a global metric PL inequality under a target density-ratio condition, a defective PL inequality yielding exponential decay of the squared MMD for uniformly positive targets, and two counterexamples showing limits of PL coercivity. On R^d it proves a radial PL inequality under source-support inclusion, and whole-space travel-time obstructions showing that no uniform convergence rate or global PL inequality can hold on the unrestricted class. Finally, it proves that every Lagrangian critical point with (ρ−µ)+ absolutely continuous equals the target, and combines this with a logarithmic-capacity tightness estimate to obtain narrow, weak-* L∞, and H^{-α} convergence in dimension two whenever the Coulomb energy is finite at some positive time.","tokens_in":44315,"tokens_out":14907,"duration_ms":146651,"significance":"If the results hold, this is a substantial contribution to the analysis of MMD gradient flows and Wasserstein gradient flows of singular interaction energies. The paper's strengths are notable: the Cauchy theory handles arbitrary Borel initial measures; the torus PL inequality is proved by an intricate variational barrier argument rather than by fitting constants; the defective PL proof gives explicit exponential rates without a lower bound on the evolving density; the radial Euclidean PL result is proved by a clean quantile-shell decomposition; and the planar convergence theorem uses a genuine logarithmic-capacity tightness mechanism. The negative results—PL failure for a target vanishing at one point, non-uniformity over targets with a fixed lower bound, and the travel-time obstruction on R^d—are constructive and clarify the sharpness of the hypotheses. The acknowledged non-uniqueness for measure initial data limits the interpretation of Theorems 1.7 and 1.22 as statements about the particular constructed solutions rather than an intrinsic semigroup, but the paper states this caveat explicitly and the theorems are correctly quantified over constructed solutions. I found no intern","major_comments":[],"minor_comments":[{"comment":"The display reads “2w − εΦ′(ρ̄) = κρ̄ dx-a.e.”, but the derivation from (3.33) gives the density (2w − εΦ′(ρ̄) − κ)ρ̄ = 0, and since ρ̄ > 0 a.e. the conclusion should be “2w − εΦ′(ρ̄) = κ dx-a.e.” The subsequent Euler-equation argument uses the correct form, so this is a harmless typo, but it should be corrected.","section":"Eq. (3.34)"},{"comment":"Definition 1.1 requires ρ_t ∈ L∞ for a.e. t > 0, while Theorem 1.2 states the bound for every t > 0. This is a minor mismatch between the weak-solution definition and the regularity actually proved; the authors may wish to harmonize the wording or explicitly state that the constructed solutions have the stronger pointwise-in-time regularity.","section":"Definition 1.1 and Theorem 1.2"},{"comment":"The non-uniqueness caveat for arbitrary Borel initial data is stated in Section 2, but it is also relevant to the long-time theorems in the introduction. Adding a sentence near Theorems 1.7 and 1.22 reminding the reader that these results concern solutions constructed by smooth approximation would improve the exposition, even though the theorems are already quantified correctly.","section":"Section 2, before §2.1"}],"recommendation":"accept","confidential_remarks":"The paper is well within the scope of the journal and the proofs are detailed and convincing. The only substantive caveat—the lack of uniqueness for measure initial data—is explicitly acknowledged and does not affect the validity of the stated theorems. I recommend acceptance; the minor typo and wording issues can be fixed in proof."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read this if you work on Wasserstein gradient flows or MMD dynamics. It resolves several genuinely open questions: global weak solutions from arbitrary Borel initial measures with instantaneous L∞ regularization, a global metric PL inequality near uniform torus targets, exponential decay for uniformly positive targets via a defective PL inequality, a radial Euclidean PL result under support inclusion, whole-space travel-time obstructions, a clean rigidity theorem for Lagrangian critical points when (ρ−μ)+ is absolutely continuous, and planar convergence in narrow, weak-*, and H^{-α} topologies. That is a lot, and the paper's proofs look real. The variational-barrier proof of Theorem 1.3 and the use of the external level-set theorem [APR20] are intricate, but I found no gap; the stress-test note's one harmless typo in the display near (3.34) does not affect the argument. The main caveat is the one the authors state in Section 2: uniqueness for global weak solutions starting from arbitrary Borel measures is not known, so Theorems 1.7 and 1.22 describe the behavior of the particular solution selected by smooth approximation, not an intrinsic semigroup for the equation. That is a real scope limitation—if different regularizations select different limits, the long-time theorems do not characterize the equation itself—but it is disclosed clearly in the text and in Section 7. It does not invalidate the stated results, which are explicitly quantified over constructed solutions. I would not treat non-uniqueness as a fatal flaw; it is a boundary of the theory, and the paper says so. I also note that the PL inequalities are proven, not fitted, and the self-citations are external published tools, not circular dependencies. The paper is honest, carefully scoped, and technically dense; it deserves a full peer review rather than a desk rejection. I would bring it to our reading group and would cite it as the current reference for these questions. Send it to a serious referee, with the expectation of a demanding but fair report.","headline":"A serious paper that largely delivers on its advertised results; the only significant caveat—non-uniqueness for measure initial data—is openly acknowledged rather than hidden.","tokens_in":44795,"tokens_out":1588,"would_cite":true,"duration_ms":19936,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35A01","35B40","49Q22","35D30"],"pacs":[],"model":"deepseek-v4-flash","headline":"Coulomb MMD flow relaxes to target on torus and in plane, with exponential rates under coercivity","keywords":["Wasserstein gradient flow","Maximum Mean Discrepancy","Coulomb potential","Polyak–Łojasiewicz inequality","long-time behavior","nonlocal transport equation","Lagrangian critical point","planar tightness"],"falsifier":"Take a singular initial measure and two sequences of smooth mollifications that converge to it; if the corresponding global weak solutions differ at positive times or have different long-time limits, the paper's convergence theorems for arbitrary initial data would be regularization-dependent rather than intrinsic. Alternatively, construct a Lagrangian critical point with (ρ−μ)+ absolutely continuous whose Coulomb field vanishes on the positive part but ρ≠μ; this would directly contradict the rigidity theorem.","tokens_in":43912,"feed_emoji":"📉","tokens_out":6982,"duration_ms":63529,"temperature":0.7,"pith_summary":"The paper studies the Wasserstein gradient flow of the squared Maximum Mean Discrepancy (MMD) with a Coulomb kernel, a nonlocal transport equation that pushes a probability density toward a target density. The main message: from arbitrary Borel initial measures the flow exists globally, instantly becomes bounded when the target is bounded, and, under coercivity conditions, relaxes to the target at an exponential rate. On the torus, a defective Polyak–Łojasiewicz inequality yields exponential decay for uniformly positive targets even when the evolving density has vacuum regions. On the whole space the situation is more delicate: radial sources supported inside the target converge exponentially, but no uniform decay rate or global PL inequality holds when the source may start far away. In dimension two, finite Coulomb energy at a positive time forces tightness, and every constructed solution converges to the target in several weak topologies.","feed_headline":"Coulomb MMD flow relaxes to target on torus and in plane","feed_subtitle":"Flow exists from any starting measure, regularizes instantly, decays exponentially for bounded positive targets.","key_machinery":"The Polyak–Łojasiewicz (PL) inequality D(ρ|μ) ≥ c MMD²(ρ,μ), comparing the dissipation D(ρ|μ)=∫|∇g*(ρ−μ)|² dρ to the squared MMD energy, is the main coercivity mechanism; where the classical inequality fails (vacuum regions, distant sources), the paper uses a defective PL inequality with an explicit exponentially decaying defect, and a metric-slope formulation to handle singular measures. The rigidity of Lagrangian critical points — ∇h=0 on (ρ−μ)+ forces ρ=μ when (ρ−μ)+ is absolutely continuous — supplies the final identification in the planar LaSalle argument, together with a logarithmic-capacity estimate that turns finite Coulomb energy into uniform tightness in two dimensions.","core_discovery":"The central discovery is that the squared-MMD Coulomb discrepancy functional, though nonconvex, has a Polyak–Łojasiewicz dissipation-to-energy structure in exactly the regimes where transport is not obstructed. The paper proves a global metric PL inequality on the torus for all finite-energy sources when the target is nearly uniform, and a defective PL inequality for uniformly positive targets whose vacuum defect decays exponentially; these yield exponential decay of the squared MMD for every constructed solution after any waiting time. It also proves rigidity of Lagrangian critical points — if (ρ−μ)+ is absolutely continuous and the Coulomb field vanishes on it, then ρ=μ — and, in the plane","pith_inferences":["The defective PL mechanism suggests that for degenerate targets (vanishing on a set), the relaxation rate may be governed by local Łojasiewicz exponents tied to the order of vanishing; the paper leaves this as an open problem, and a natural next step is to compute such rates for targets vanishing like |x|^k.","Because uniqueness for singular initial data is unresolved, the planar convergence theorem is a statement about the particular approximation-selected solution; if selection is later shown to matter, the theorem would need reinterpretation as a property of a canonical semigroup rather than the equation itself.","The travel-time obstruction at spatial infinity implies that any data-independent convergence rate for the whole-space problem must encode location or tail information; in sampling and generative-modeling applications, this quantifies the cost of initializing far from the target.","The rigidity theorem classifies Lagrangian critical points with absolutely continuous positive part, but singular positive parts (e.g., atoms) can yield non-target critical points for suitable field representatives, so the critical-point landscape outside the absolutely continuous class remains open."],"forward_implications":["On the torus, every solution constructed from an arbitrary probability measure and targeting a bounded, uniformly positive density attains exponential decay of the squared MMD after any positive time, with any rate below 2λ/3 where λ is the target's essential infimum.","Global PL coercivity holds on the torus for all finite-energy sources when the target is close to uniform, and in dimension one for every bounded uniformly positive target; this gives quantitative transport-information bounds independent of the source's lower bound.","In the plane, finite Coulomb energy at one positive time is enough to guarantee convergence of the constructed solution to the target in narrow, L∞ weak-*, and H^{-α} topologies; no moment or support assumptions are needed.","On the whole space, any uniform exponential decay modulus or global PL inequality fails: a localized source separated from a compactly supported target by distance D retains a positive fraction of its initial squared MMD for times proportional to D.","Radial sources with connected target support and source-support inclusion converge exponentially on R^d, d≥2, with explicit constant max{1/λ, d²A/(2λ²)}."],"fun_headline_variants":["Coulomb MMD flow hits global PL on torus, decays exponentially","Exponential relaxation for Coulomb MMD flow on torus and plane","PL inequality unlocks exponential decay in Coulomb MMD flow","Global PL inequality drives Coulomb MMD flow to target","Coulomb MMD flow: from any start, decay to target on torus"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The long-time convergence results for arbitrary Borel initial measures apply to a solution obtained by a specific smoothing procedure, and the paper does not prove that this limiting solution is unique; if a different regularization produced a different trajectory, the convergence statements would not characterize the equation itself.","fun_headline_variants_meta":{"raw":{"variants":["Coulomb MMD flow hits global PL on torus, decays exponentially","Exponential relaxation for Coulomb MMD flow on torus and plane","PL inequality unlocks exponential decay in Coulomb MMD flow","Global PL inequality drives Coulomb MMD flow to target","Coulomb MMD flow: from any start, decay to target on torus"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000185,"raw_usage":{"total_tokens":1216,"prompt_tokens":861,"completion_tokens":355,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":605,"completion_tokens_details":{"reasoning_tokens":266}},"tokens_in":605,"tokens_out":355,"duration_ms":4020,"temperature":1.0,"reasoning_tokens":266,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-02T06:25:06.834944+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take a singular initial measure and two sequences of smooth mollifications that converge to it; if the corresponding global weak solutions differ at positive times or have different long-time limits, the paper's convergence theorems for arbitrary initial data would be regularization-dependent rather than intrinsic. Alternatively, construct a Lagrangian critical point with (ρ−μ)+ absolutely continuous whose Coulomb field vanishes on the positive part but ρ≠μ; this would directly contradict the rigidity theorem.","supporting_citations":[],"review_version":2}