{"id":"39bfe696-6fba-419e-9c64-ed92eb1af9f7","arxiv_id":"2502.00582","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"For consensus-based optimization on a bounded domain, the particle system's empirical distribution is within O(1/N) of its mean-field limit uniformly in time, and it concentrates about a limit point as N and t grow.","lead":"An optimization algorithm that uses a swarm of particles without gradients is shown to approximate its infinite-particle limit at rate 1/N for all time, not just short runs. This separates the choice of particle count from the running time in consensus-based optimization.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Proposition 2.4's exponential contraction is an unproven adaptation of CCTT18 to the cut-off dynamics, and every uniform-in-time estimate in the paper depends on it.","rationale":"The paper is a serious attempt to extend Delarue-Tse's uniform-in-time methodology to consensus-based optimization. The master-equation decomposition and the control of the second-order derivatives are internally coherent, and the proof structure is plausible. The most fragile point, however, is the imported contraction estimate. Proposition 2.4 is not proved in the paper; Remark 2.5 merely says it is an adaptation of CCTT18's Theorem 4.1. The cutoff phi changes the diffusion coefficient: it vanishes on the boundary of the search domain, making the noise degenerate exactly at the boundary. Whether the exponential convergence to a unique Dirac persists under this modification, with constants uniform in the initial measure, is not established. Every core lemma (3.1, 3.3, 3.4, and the lemmas in Section 2) explicitly uses the decay (Erg) to control remainder terms and to pass from delta_{tilde x_mu} to the evolving law mu_t; without it, the master-equation integrals in the proof of Theorem 2.6 would not be uniformly bounded in time. This agrees with the reader's weakest_assumption, though I sharpen it to the cut-off-specific gap. The abstract overclaim about convergence to the global minimizer is real but secondary: the paper itself later notes that the limit point is alpha-dependent, so the fix is editorial. The requested remedy—a firmer proof of Lemma 2.15 and a corrected abstract—should include an explicit verification of Proposition 2.4 for the truncated dynamics. Thus the conditional-acceptance verdict is appropriate.","tokens_in":37731,"tokens_out":14118,"duration_ms":136257,"concrete_test":"Provide a self-contained proof of Proposition 2.4 for the cut-off mean-field SDE (2.2): for every mu0 in P(B(c0,2rcut)) and lambda satisfying the stated lower bound, prove that there exist kappa, C > 0, independent of mu0, such that W2(mu_t, delta_{tilde x_mu}) <= C e^{-kappa t}, explicitly handling the boundary-vanishing cutoff phi in the diffusion estimate. A full derivation from CCTT18's Theorem 4.1, with a detailed verification that the presence of phi does not break the contraction argument, would also settle the concern.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Theorem 2.6's uniform-in-time O(N^{-1}) bound rests on Proposition 2.4, which postulates exponential decay of the mean-field cut-off process (2.2) to a distribution-dependent Dirac delta_{tilde x_mu}. The paper does not prove this: Remark 2.5 calls it an adaptation of Theorem 4.1 of CCTT18. But CCTT18's theorem is for the standard CBO dynamics without the compactly supported cutoff phi in the diffusion coefficient. In (2.2), phi vanishes on the boundary of B(c0,2rcut), changing the degeneracy of the noise exactly where particles are confined. It is not automatic that the same contraction holds with constants C and rate kappa uniform in mu0, nor that the invariant Dirac is unique for the cut-off dynamics. This is not a harmless technicality: Lemmas 3.1, 3.3, 3.4 and 2.12-2.15 all use the exponential decay (Erg) to make remainder terms integrable and to replace delta_{tilde x_mu} by mu_t; Theorem 2.6 then uses those bounds to make the master-equation integrals finite uniformly in t. If Proposition 2.4 fails, or if its rate depends on initial data or time, the argument yields only finite-horizon propagation of chaos, not the stated sup_t >= 0 bound. Since the contraction is imported rather than derived for the modified dynamics, the central uniform-in-time claim currently rests on an unverified assumption.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the cut-off consensus-based optimization particle system (2.1) on a bounded domain and establishes a uniform-in-time weak propagation of chaos result: for test functionals Φ satisfying (2.4)-(2.5), it claims sup_{t≥0} |E[Φ(ν_t^N)] - Φ(ν̄_t)| ≤ C_main/N for all N ≥ 2, with C_main independent of the initial distribution. The proof follows the Delarue-Tse methodology: decompose the weak error via the master equation, express second-order derivatives of U through linearized Fokker-Planck equations, establish exponential decay of their Sobolev dual norms using the ergodicity of the mean-field CBO flow, and conclude by Grönwall arguments. Corollaries give convergence of the centered empirical measure to δ_0 in Fourier-Wasserstein and Wasserstein distances with rates in N and t separately.","tokens_in":38002,"tokens_out":4854,"duration_ms":52615,"significance":"If the technical gaps are filled, this is a valuable contribution: it provides one of the first uniform-in-time propagation of chaos results for CBO on a bounded domain, with an explicit O(N^{-1}) rate and independent choices of N and t. The paper is substantial and well-structured, adapting the DT25 machinery with detailed GBM estimates in Appendix B, and it contains no fitted parameters. The explicit nature of the estimates and the careful master-equation decomposition are strengths. However, the central uniform-in-time claim rests on an imported contraction result for the cut-off dynamics and on a partially sketched proof of Lemma 2.15; both are load-bearing and need to be addressed before the result can be considered established.","major_comments":[{"comment":"Proposition 2.4 asserts for the cut-off mean-field SDE (2.2) exponential convergence of the mean and consensus to a point x̃_μ with rate κ = 2(λ - dσ²e^{9αc_E r_cut² - αE}) and a constant C uniform in μ0, but no proof is given; Remark 2.5 calls it an adaptation of Theorem 4.1 of CCTT18. That theorem is stated for the standard CBO dynamics without the compactly supported cut-off φ in the diffusion coefficient, and the cut-off changes the noise degeneracy exactly on the boundary of the confinement. Since the exponential decay (4.1) is used throughout Lemmas 3.1, 3.3, and 3.4 and ultimately in Theorem 2.6 to obtain integrability uniformly in t, this is a load-bearing point. The authors should either prove Proposition 2.4 for (2.2) with the stated uniformity, or identify a published result that applies verbatim to the cut-off dynamics; a short citation is not sufficient.","section":"Section 2.2, Proposition 2.4"},{"comment":"The proof of Lemma 2.15 is incomplete. Section 3.3 provides only a sketch, and Appendix A contains a 'Substitute proof of Lemma 2.15' that says 'By replacing all m(1)(t; μ, δz) with d(1)_j(t; μ, z) in the above proof, we get exactly (3.7) for Lemma 2.15.' This does not verify that the remainder expansion, the cancellation steps, and the associated Sobolev-norm bounds of the Lemma 2.14 proof carry over to d(2), particularly the treatment of the initial condition q(1)_{j,∞}·∇∂_{x_j}δ_{z_1} and the separate bounds (3.7) and (3.8). Lemma 2.15 is used in the proof of Theorem 2.6 to obtain the exponential decay of ∂_{z1_j}∂_{z2_j}δ²U/δm², which is essential for the uniform-in-time bound in (2.9). The authors should provide a complete, self-contained proof.","section":"Section 3.3 and Appendix A, Lemma 2.15"},{"comment":"In the treatment of the case q0 = ∂²_{x_j}δ_z, the displayed formula for γ^rem_j(T0; z) writes the first integral with upper limit T instead of T0, while the subsequent Cauchy-sequence argument treats T1,T2 → ∞. If the integral is genuinely truncated at T, the convergence argument does not apply; if it is a typo and should read T0, the correction should be made. In addition, the bound for the second-derivative terms is asserted to follow from Lemma B.1, but the display records only the expectation identities and not the uniform-in-x estimates needed for the dual-norm conclusion; this step should be expanded explicitly.","section":"Section 4, Lemma 3.1, Step 3"}],"minor_comments":[{"comment":"The proof relies on Lemma 4.11 of DT25 for the identities (2.11), but it does not explicitly check that the functionals satisfying (2.4)-(2.5) lie in the domain of the required measure derivatives of the flow map P_t. The authors should add a short justification or cite the precise result ensuring the needed regularity of U.","section":"Section 2.4, Proof of Theorem 2.6"},{"comment":"The terminology 'centered Wasserstein distance' is used for W2(ν_t^N ∘ τ_{⟨Id,ν_t^N⟩}, δ0)^2, which is the variance functional; the name may confuse readers, since it is not a metric between the centered law and δ0 in the usual Wasserstein sense. A brief clarification would help.","section":"Corollary 2.8"},{"comment":"The statement of Lemma 2.12 includes 'for every j ∈ [d]' although the quantity m(1) involves no index j; this is a copy-paste artifact and should be removed.","section":"Lemma 2.12"},{"comment":"The manuscript contains numerous typographical errors and OCR artifacts, e.g., 'Prop agation' in the title/abstract, 'W e' and 't he' in the abstract, and 'Lemmata' vs 'lemmas'. A careful proofreading pass is needed.","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":"The paper is a serious and substantial application of the Delarue-Tse uniform-in-time propagation of chaos methodology to consensus-based optimization. The main reason for major revision is the unproven reliance on Proposition 2.4 for the cut-off dynamics, which is central to all uniform-in-time estimates, and the incomplete proof of Lemma 2.15. Both issues are fixable in scope but require real mathematical work, not merely exposition. I recommend major revision rather than rejection because the overall architecture is coherent and the result would be significant if these gaps are closed."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This paper is worth reading: it gives the first uniform-in-time weak propagation of chaos bound for consensus-based optimization on a bounded domain, with rate O(1/N). The adaptation of Delarue–Tse to a non-toroidal, non-fully-diffusive setting is nontrivial, and the master-equation/linearized-Fokker–Planck machinery is handled carefully. The proof that the second derivative of the value function decays exponentially is clever, and the O(1/N) rate is genuinely derived, not fitted. The authors are also honest in the body about the finite-α bias, the cut-off, and the need for exponential contraction.\n\nThe main soft spot is Proposition 2.4, on which everything rests. It postulates exponential decay of the mean-field cut-off dynamics toward a point mass, with constants uniform in the initial law, citing CCTT18. But CCTT18's theorem is for the uncut CBO dynamics; adding the compactly supported φ changes the diffusion coefficient and the degeneracy structure exactly where particles are confined. The remark 'this is an adaptation' is not enough, because every uniform-in-time estimate, from Lemma 3.1 to the final Theorem 2.6, uses this contraction. A referee should demand a proof or at least a careful statement of how the CCTT18 estimates transfer to the cut-off system.\n\nSecond, Lemma 2.15 is proved by 'substitute proof': replace m with d in the proof of Lemma 2.14. That is too quick for a key estimate, especially since Lemma 2.14 is already a lengthy arrangement of cancellations. The extra argument for the linear-in-t bound (3.8) is helpful, but the exponential decay (3.7) needs its own write-up.\n\nThird, the abstract says the empirical distribution converges to the Dirac at the global minimizer. The theorem actually gives concentration around some limit point that may depend on the initial data and carries an α-dependent bias. The centered corollaries remove the mean, but they do not put the particles at x*.\n\nThese are fixable issues. The mathematical core is plausible and the methodology is relevant. This is a paper for anyone working on mean-field limits of optimization algorithms or on uniform-in-time PoC beyond fully diffusive settings. It deserves a serious referee, though I would recommend major revision: prove or precisely state Proposition 2.4, expand Lemma 2.15, and correct the abstract.","headline":"A serious and largely convincing uniform-in-time PoC result for cut-off CBO, but the central contraction estimate is imported without proof and Lemma 2.15 is sketched by analogy; the abstract also overclaims for the global minimizer.","tokens_in":38548,"tokens_out":4820,"would_cite":false,"duration_ms":53760,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35Q89","37N40","93D50","82C31","90C26"],"pacs":[],"model":"deepseek-v4-flash","headline":"Consensus-based optimization particles track their mean-field limit with error at most C/N, uniformly for all times.","keywords":["consensus-based optimization","propagation of chaos","uniform-in-time","mean-field limit","Fokker-Planck equation","master equation","Sobolev spaces","global optimization"],"falsifier":"Run the cut-off CBO particle system with a drift $\\lambda$ below the threshold in Proposition 2.4 on a non-convex objective and compute $\\sup_{t\\ge 0} |\\mathbb{E}[\\Phi(\\nu_t^N)] - \\Phi(\\bar\\nu_t)|$ for the smooth translation-invariant centered Fourier–Wasserstein functional; if the supremum grows with $t$ for fixed $N$, the main theorem is false. Equivalently, check whether the mean-field consensus point $M(\\mu_t)$ converges to a single point exponentially fast with a rate independent of the initial measure, since Lemma 3.3 uses exactly this to make the remainder $R^{1,i}_t$ decay.","tokens_in":37512,"feed_emoji":"🎯","tokens_out":12975,"duration_ms":151640,"temperature":0.7,"pith_summary":"Consensus-based optimization (CBO) is a gradient-free algorithm in which $N$ particles drift toward a weighted consensus point to minimize a possibly non-convex objective. This paper proves that, on a bounded search domain cut off by a smooth function, the empirical distribution of the CBO particle system stays within weak error $O(N^{-1})$ of its mean-field limit uniformly over the infinite time horizon, with a constant independent of both time and the initial distribution. The number of particles can therefore be chosen without regard to how long the algorithm runs, eliminating the finite-horizon trade-off in earlier estimates. As a consequence, the empirical measure converges, simultaneously in $N$ and $t$, to the Dirac measure at the global minimizer in centered Wasserstein-type metrics.","feed_headline":"Particle swarm's weak error stays O(1/N) forever","feed_subtitle":"Uniform-in-time bound frees the number of particles from the running time in consensus-based optimization.","key_machinery":"The machinery is the linearized Fokker–Planck equation (L-FPE) for fluctuations of the empirical measure around the mean-field flow, studied through its backward adjoint equation. The proof obtains estimates of the form $\\|q_t - q_\\infty \\cdot \\nabla \\delta_{\\tilde x_\\mu}\\|_{(n,\\infty)'} \\le C e^{-\\kappa_0 t}$ for the L-FPE solutions, an ergodicity property that converts finite-horizon estimates into uniform-in-time ones. Because the CBO generator is not fully diffusive (its second-order coefficient vanishes at the consensus point), the backward equation does not decay directly; instead the paper uses the Feynman–Kac formula to show exponential decay of derivatives of the adjoint solution along geometric Brownian motions, at rate $e^{-\\lambda t}$, and combines this with the exponential contraction of the mean-field CBO flow to a Dirac measure (Proposition 2.4).","core_discovery":"The central discovery is a uniform-in-time weak propagation of chaos for the cut-off consensus-based optimization system. For any smooth, translation-invariant functional $\\Phi$ of the empirical measure, the paper establishes $\\sup_{t\\ge 0} |\\mathbb{E}[\\Phi(\\nu_t^N)] - \\Phi(\\bar\\nu_t)| \\le C_{\\mathrm{main}}/N$ for every $N \\ge 2$, with $C_{\\mathrm{main}}$ independent of the initial law. The proof identifies the second-order derivative of the solution operator $U(t,\\mu)=\\Phi(\\mu_t)$ as the key object, decomposes it through the master equation into solutions of linearized Fokker–Planck equations, and shows that those solutions are uniformly bounded and decay exponentially in Sobolev dual norms to a term of the form $q_\\infty \\cdot \\nabla \\delta_{\\tilde x_\\mu}$. This exponential decay is what prevents errors from accumulating over long time intervals.","pith_inferences":["A natural testable extension is to remove the cutoff function $\\varphi$: if the objective grows fast enough at infinity to control the exponential moments appearing in the estimates, the same $O(N^{-1})$ uniform bound might hold on the whole space.","The translation-invariance condition (2.5) suggests that functionals sensitive to the mean, such as uncentered Wasserstein distances, would not satisfy the theorem; one could test whether their weak error grows with $t$.","The exponential-decay estimates on L-FPE solutions could be reused to prove uniform-in-time propagation of chaos for related consensus algorithms (minimax CBO, constrained CBO), as long as their mean-field dynamics contract exponentially to a Dirac measure.","If the contraction rate $\\kappa$ in Proposition 2.4 can be bounded below explicitly in terms of $\\lambda, \\sigma, \\alpha$ and the objective's parameters, the constant $C_{\\mathrm{main}}$ becomes computable and the result turns into a practical resource-allocation rule for CBO users."],"forward_implications":["For any prescribed error tolerance, the number of particles $N$ can be fixed independently of the running time $t$, eliminating the finite-horizon trade-off in which $N$ had to grow exponentially with $t$.","The empirical measure of the CBO system converges jointly in $N$ and $t$: $\\mathbb{E}\\|\\nu_t^N \\circ \\tau_{\\langle \\mathrm{Id}, \\nu_t^N\\rangle} - \\delta_0\\|^2_{-s,2} \\le C_{FW}(N^{-1} + e^{-2\\kappa t})$ and $\\mathbb{E} W_2^2(\\nu_t^N \\circ \\tau_{\\langle \\mathrm{Id}, \\nu_t^N\\rangle}, \\delta_0) \\le C_W(N^{-1} + e^{-\\kappa t})$.","The uniform-in-time estimate holds for any initial distribution supported in the search ball, so the same $N$ works for any restart or initialization.","The proof transfers a recently developed uniform-in-time weak propagation methodology from torus settings to compactly supported, degenerate-diffusion CBO dynamics."],"supporting_citations":[{"why":"Supplies the master-equation and linearized Fokker–Planck methodology (Propositions 3.3 and 3.4, Lemma 4.11) that the paper adapts to the CBO setting.","marker":"[DT25]"},{"why":"Provides the exponential contraction of the mean-field CBO flow to a Dirac measure used as Proposition 2.4 and in Corollary 2.8.","marker":"[CCTT18]"},{"why":"Gives Theorem 2.14, the second-order expansion on the space of measures that bounds the static part of the error in (2.7).","marker":"[CST22]"},{"why":"Furnishes the higher-order regularity of nonlinear Fokker–Planck equations with respect to the measure argument, underlying the L-FPE analysis.","marker":"[Tse21]"},{"why":"Supplies the global-convergence and exponential behaviour of the CBO mean-field limit supporting Proposition 2.4 and the corollaries.","marker":"[FKR24]"}],"fun_headline_variants":["Uniform-in-time chaos: weak error stays O(1/N)","CBO weak error stays O(1/N) for all time","Consensus optimization: uniform-in-time error bound O(1/N)","Weak chaos propagation: 1/N error forever, no time drift","Particle swarm error remains O(1/N) uniformly in time"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole argument rests on the exponential contraction of the mean-field CBO flow toward a single Dirac measure, imported under a large-drift condition on $\\lambda$; if that contraction rate were zero, negative, or time-dependent, the uniform-in-time $O(N^{-1})$ bound would collapse to a finite-horizon estimate.","fun_headline_variants_meta":{"raw":{"variants":["Uniform-in-time chaos: weak error stays O(1/N)","CBO weak error stays O(1/N) for all time","Consensus optimization: uniform-in-time error bound O(1/N)","Weak chaos propagation: 1/N error forever, no time drift","Particle swarm error remains O(1/N) uniformly in time"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000831,"raw_usage":{"total_tokens":3593,"prompt_tokens":877,"completion_tokens":2716,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":493,"completion_tokens_details":{"reasoning_tokens":2636}},"tokens_in":493,"tokens_out":2716,"duration_ms":19871,"temperature":1.0,"reasoning_tokens":2636,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-09T18:25:55.776641+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the cut-off CBO particle system with a drift $\\lambda$ below the threshold in Proposition 2.4 on a non-convex objective and compute $\\sup_{t\\ge 0} |\\mathbb{E}[\\Phi(\\nu_t^N)] - \\Phi(\\bar\\nu_t)|$ for the smooth translation-invariant centered Fourier–Wasserstein functional; if the supremum grows with $t$ for fixed $N$, the main theorem is false. Equivalently, check whether the mean-field consensus point $M(\\mu_t)$ converges to a single point exponentially fast with a rate independent of the initial measure, since Lemma 3.3 uses exactly this to make the remainder $R^{1,i}_t$ decay.","supporting_citations":[],"review_version":1}