{"id":"7678b40b-a5a3-4c78-8567-bd322d88024d","arxiv_id":"2502.09300","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"For SDEs on R^d, the paper proves existence and uniqueness of the optimal infinitesimal local kernel perturbation maximizing the linear response of an observable, and demonstrates a numerical approximation on a double-well example.","lead":"This paper finds which small, localized change to the transition kernel of a stochastic differential equation produces the largest change in the expected value of a chosen observable, and proves the optimal perturbation is unique. It gives a numerical recipe for computing this perturbation and tests it on a double-well SDE.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Lemma 14's spectral-gap proof is the load-bearing gap: B is not compactly embedded in L1, and the strong-norm strict contraction used for simplicity is unjustified.","rationale":"The reader's weakest_assumption pointed to Lemma 14, and I agree this is the most load-bearing step: every later result, including Theorem 25 and Proposition 28, inherits the spectral gap of L0 on B. I part company on the precise diagnosis. The standard positivity argument does give a pointwise strict L1 contraction for zero-mean functions, which is already enough to rule out other unit-circle eigenvalues and non-unique fixed points; the problem is that the proof states this contraction in the strong norm ||.||_s, where it is not justified. More seriously, the invocation of Hennion's theorem is not backed by a compact embedding B -> L1; the bounded sequence f_n = n 1_{[x0, x0+1/n]} with x0 outside D is a concrete counterexample to such an embedding. Proposition 12 provides a weaker compactness property, and it is possible that a suitable spectral theorem accepts that property, but the paper does not say so. These are proof gaps rather than demonstrated falsehoods: the L1 contraction and the dissipativity of the SDE are strong hints that the conclusions are true and repairable. The numerical section is suggestive but not conclusive, and the false implication in Proposition 23 is an additional blemish. For these reasons the appropriate verdict remains CONDITIONAL, unchanged from the reader's assessment.","tokens_in":23262,"tokens_out":12997,"duration_ms":139712,"concrete_test":"Check the hypotheses of the version of Hennion's theorem cited as [10, Theorem B.14]. If it requires a compact embedding B -> L1, the proof of Lemma 14 fails immediately, since the sequence f_n = n 1_{[x0, x0+1/n]} with x0 outside D is bounded in B and has no L1-convergent subsequence; if the theorem instead only requires relative compactness of L(B) in L1, verify that this condition is exactly Proposition 12 and supply the missing argument. Independently, for the Ornstein-Uhlenbeck example of Remark 18, compute the ratio ||L w_n||_s / ||w_n||_s for w_n = Gaussian(x0+n, sigma) - Gaussian(x0-n, sigma) as n -> infinity; if the ratio does not stay strictly below 1, the strong-norm strict contraction asserted in Lemma 14 is false, although the spectral-gap conclusion may still be repairable by a different argument.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central results (Theorem 25 and Proposition 28) require L0 to have a spectral gap on the strong space B with 1 a simple isolated eigenvalue. Lemma 14 is the only place this is established, and its proof has two unsupported steps. First, Hennion's theorem is invoked via Lemma 13, but the standard hypothesis is a compact embedding of B into the weak space L1; this embedding fails. For x0 outside D, the sequence f_n = n 1_{[x0, x0+1/n]} satisfies ||f_n||_s = ||f_n||_{L1_2} + ||1_D f_n||_2 = O(rho_2(x0)) < infinity, yet f_n has no convergent subsequence in L1 (it concentrates to a Dirac mass). Proposition 12 only proves that L(B) is relatively compact in L1, not that B -> L1 is compact, and no argument shows that the cited theorem accepts this weaker condition. Second, the proof that 1 is the only eigenvalue on the unit circle and is simple asserts ||Lw||_s < ||w||_s for nonzero zero-mean w. Positivity of the kernel gives pointwise strict inequality and hence strict L1 contraction, but it does not give strict contraction in the weighted norm ||w||_s = ||w||_{L1_2} + ||1_D w||_2; L|w| may have larger weighted norm than |w|, and the L2-term on D is unrelated to positivity. The uniqueness conclusion may be recoverable from the L1 contraction together with the compactness of L(B), but as written Lemma 14 does not prove the spectral gap on which the linear response and optimization arguments rest. Proposition 23 also contains a false implication (zero integral does not imply zero L1 norm), but this is separate and more easily repaired.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies a dissipative stochastic differential equation on R^d and its annealed transfer operator L at time one, viewed as a kernel operator. The authors construct a 'strong' Banach space B of densities with a weighted L1 norm and an L2 norm on a compact set D, and cite results from prior work for Lasota-Yorke inequalities and BV-type compactness. They then consider infinitesimal perturbations of the kernel supported in D, prove a linear-response formula of the form (Id - L0)^{-1} \\dot L f0 in L1 (Theorem 25), and show that the problem of maximizing the response of an observable over a closed bounded strictly convex set P of perturbations has a unique solution (Proposition 28). A Fourier-based numerical scheme for approximating the optimal perturbation is presented and illustrated on a double-well SDE.","tokens_in":23573,"tokens_out":12793,"duration_ms":139771,"significance":"If fully justified, the paper extends the optimal-response framework from compact phase spaces to noncompact R^d for SDEs, which is a genuinely useful step for climate and control applications. The explicit linear-response formula and the constructive numerical algorithm, with code available on Zenodo, are notable strengths. The mathematical architecture is appropriate: Lasota-Yorke inequalities, resolvent bounds, and convex optimization are the right tools, and most of Sections 4 and 5 are standard once the spectral gap is granted. The main caveat is that two load-bearing proofs, the spectral-gap proof in Lemma 14 and the uniqueness proof in Proposition 23, are not rigorous as written. These should be fixed before the central claims can be considered established.","major_comments":[{"comment":"The simplicity argument for the eigenvalue 1 asserts the inequality ||L(u-v)||_s < ||u-v||_s from positivity of the kernel. This does not follow. Positivity gives pointwise |Lg| < L|g| for nonzero zero-mean g and hence strict contraction in L1, but it gives no such strict contraction in the strong norm ||g||_s = ||g||_{L1_2} + ||1_D g||_2; the L2(D)-component is unrelated to the positivity argument and can increase. In fact, if u-v is a fixed point, then ||L(u-v)||_s = ||u-v||_s, contradicting the asserted strict inequality. The simplicity of the eigenvalue 1 can likely be proved using the strict L1 contraction instead, but as written the proof is invalid.","section":"Lemma 14"},{"comment":"The proof that L has a spectral gap on B is incomplete. The cited Hennion theorem, in its standard form, requires a compact embedding of the strong space into the weak space; that embedding is false here. For x0 outside D, the sequence f_n = n 1_{[x0, x0+1/n]} satisfies ||f_n||_s = O(rho_2(x0)) < infinity yet has no subsequence convergent in L1, since it converges weakly to a Dirac mass. Proposition 12 proves only compactness of L(B) in L1, not compactness of the embedding B into L1, and its proof itself contains an invalid step: a Cauchy sequence in L1 is asserted to converge in BV2, and the closed unit ball of B is asserted to be sequentially compact in ||.||_s. If the authors intend a variant of Hennion's theorem that requires only relative compactness of L(B) in L1, that variant and its hypotheses should be stated explicitly and verified. As written, the spectral gap and hence Proposition 17, Theorem 25, and Proposition 28 rest on an unproved assumption.","section":"Lemma 14 / Proposition 12"},{"comment":"The uniqueness paragraph contains a false implication: from int(f_delta - g_delta) = 0 the text concludes ||f_delta - g_delta||_1 = 0. This is not valid; a nonzero function with zero integral exists in any nontrivial L1 space. Moreover, the perturbed kernel kappa_delta is not shown to be positive, so the positivity-based contraction argument for L0 cannot be applied to L_delta. Uniqueness of the invariant density f_delta may be recoverable from the simplicity of the eigenvalue 1 obtained from the Keller-Liverani stability statement immediately preceding, but that connection is not made, and the proof as written is unsupported.","section":"Proposition 23"},{"comment":"The text states that 'B is compactly immersed inside L1 (Lemma 12)', but Lemma 12 is the compactness of L(B) in L1, not the compactness of the embedding B into L1, which is false. This matters because the application of [25, Theorem 1] in Proposition 23 and Lemma 24 requires a correct compactness hypothesis for the perturbed operators. The authors should spell out the uniform conditions, for example a uniform Lasota-Yorke inequality together with relative compactness of L_delta(B) in L1, that justify the eigenvalue stability and the resolvent convergence. Without this, the perturbed spectral theory used in Theorem 25 is not rigorously grounded.","section":"Proposition 23 / Lemma 24"}],"minor_comments":[{"comment":"The statement contains a typo: the limit expression at the end repeats (Id - L_delta)^{-1} twice; the second occurrence should be (Id - L0)^{-1}.","section":"Lemma 24"},{"comment":"In the proof, 'By Theorem 13' should read 'By Lemma 13'.","section":"Lemma 22"},{"comment":"The last display writes a bound of the form C/delta ||1_D f0||_2 -> 0, which is inconsistent with the hypothesis ||r_delta||_2 = o(delta). The correct factor is ||r_delta/delta||_2 -> 0; please correct this line.","section":"Lemma 20"},{"comment":"The phrase 'compact inclusion of the strong space B in L1' in the introduction to Section 3.2 is misleading, since B is not compactly embedded in L1; what is needed is the compactness of L(B) in L1.","section":"Section 3.2"},{"comment":"The Zenodo code is reported to run with Delta x = 8e-3 while the figures use Delta x = 2e-3; the footnote explains how to modify the code, but for reproducibility it would be preferable to either provide the exact code used or explicitly state that the published figures require a parameter change.","section":"Section 6 / Code availability"}],"recommendation":"major_revision","confidential_remarks":"The paper is a good fit for a mathematical physics or nonlinear science journal. The main idea is promising and the numerical work is a genuine addition, but the spectral-gap proof in Lemma 14 and the uniqueness proof in Proposition 23 need to be repaired. I see no grounds for rejection: the errors appear fixable, but they are load-bearing because Theorem 25 and Proposition 28 inherit from them. I would encourage the authors to state the precise Hennion and Keller-Liverani theorems they use, including their compactness hypotheses, and to replace the false zero-integral implication in Proposition 23."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First, what's new: this is the first paper to set up and solve an optimal linear response problem for SDEs on non-compact R^d, with local kernel perturbations supported on a compact set. The main theorems (Theorem 25 and Proposition 28) give a clean formula for the response and existence/uniqueness of the optimal perturbation. The numerical pipeline is concrete, with code on Zenodo. The paper also does a good job of explaining the limiting assumptions (perturbations don't correspond to SDEs; drift perturbations left to future work) and relating to [2,12,13].\n\nThe soft spots are in the proofs. Lemma 14 is load-bearing: it claims a spectral gap on B via Hennion's theorem. The compact embedding B -> L1 is not proved and in fact fails — take x0 outside D and f_n = n on [x0,x0+1/n]; ||f_n||_s stays bounded but no subsequence converges in L1. Proposition 12 only gives compactness of L(B) in L1, not the embedding. Also, the proof of simplicity of the eigenvalue 1 uses ||L(u-v)||_s < ||u-v||_s, which does not follow from positivity of the kernel; positivity gives L1 contraction, not contraction in the weighted norm with its L2 term on D. So the spectral gap, and everything built on it, is not established as written.\n\nProposition 23's uniqueness proof contains a false inference: integral zero does not imply L1 norm zero. The subsequent L-infinity estimate goes through only if the L1 norm is zero, so uniqueness is also unproved. These gaps are probably repairable: the spectral stability theorem of Keller-Liverani [25] plus L1 contraction may do the job, and uniqueness can perhaps be obtained from spectral simplicity. But as written the paper doesn't close the argument.\n\nThe numerics look plausible but the section is thin: no convergence checks, no error diagnostics, and the normalization of the densities in Figure 2d is a bit odd. That's a minor concern relative to the proof gaps.\n\nBottom line: the paper deserves a serious referee. The idea is good, the setting is the right one, and the flaws look fixable. I wouldn't cite the current version in my own work, but I'd cite the corrected version. Send it to review, and expect a major revision.","headline":"Useful extension of optimal linear response to non-compact R^d with a concrete numerical scheme, but the key spectral gap proof has a hole and one uniqueness proof contains a false inference; likely repairable, but needs major revision.","tokens_in":24143,"tokens_out":3334,"would_cite":false,"duration_ms":34599,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["37H10","37C30","37M05","37N35","49N45","60H10"],"pacs":[],"model":"deepseek-v4-flash","headline":"For a dissipative SDE on R^d, the invariant density responds to local kernel perturbations with an explicit resolvent formula, and the optimal-perturbation problem has a unique solution.","keywords":["stochastic differential equations","transfer operators","linear response","optimal response","kernel perturbation","spectral gap","invariant density","Fokker-Planck equation"],"falsifier":"For the one-dimensional Ornstein-Uhlenbeck SDE $dX=-X\\,dt+dW$, compute or bound $\\|L|_{V_B}\\|_s$, the strong-norm operator norm of the transfer operator restricted to zero-average densities. If this norm is not strictly less than 1, the contraction premise behind Lemma 14 is false and the spectral gap, hence linear response and uniqueness of the optimal perturbation, must be established by another mechanism; the paper's own Remark 18 shows the analogous $L^2$ norm does not contract.","tokens_in":23033,"feed_emoji":"🎯","tokens_out":12255,"duration_ms":98281,"temperature":0.7,"pith_summary":"This paper studies an inverse control problem for stochastic differential equations: given a dissipative SDE on $\\mathbb{R}^d$ and an observable $\\varphi$, which small local change to the transition kernel produces the largest change in the expectation of $\\varphi$? The authors prove that the invariant density responds differentiably to such kernel perturbations, with the response given explicitly by $(\\mathrm{Id}-L_0)^{-1}\\dot{L}f_0$ in $L^1$, and that over a closed, bounded, strictly convex set of allowed perturbations the optimal perturbation exists and is unique. This extends optimal-response theory from compact phase spaces to the noncompact space $\\mathbb{R}^d$, the natural setting for many SDE models. The paper also gives a numerical recipe and demonstrates it on a double-well gradient SDE.","feed_headline":"A unique local kernel tweak optimizes SDE response","feed_subtitle":"Proof plus numerical recipe for steering an observable's expectation in dissipative SDEs on R^d.","key_machinery":"The load-bearing object is the annealed transfer operator $L_0 f(y)=\\int \\kappa(x,y)f(x)\\,dx$, whose kernel $\\kappa$ is the time-1 transition density of the SDE; two-sided Gaussian bounds and gradient estimates give it enough regularity to act on a strong space $B=\\{f\\in L^1_2(\\mathbb{R}^d): \\|f\\|_{L^1_2}+\\|\\mathbf{1}_D f\\|_2<\\infty\\}$. On this space a Lasota-Yorke inequality together with a compact embedding into $L^1$ yields a spectral gap, with the resolvent $(\\mathrm{Id}-L_0)^{-1}$ bounded on zero-average densities. The response operator $R(\\dot{\\kappa})=(\\mathrm{Id}-L_0)^{-1}\\dot{L}f_0$ is then a bounded linear map $L^2(D\\times D)\\to L^1$, so maximising $\\int \\varphi\\,R(\\dot{\\kappa})$ over a closed bounded strictly convex set is a strictly convex optimisation problem solved by the normalised Riesz representer $g/\\|g\\|_2$, whose Fourier coefficients are computable as $G_r = J(h_r)$ for an orthonormal basis $(h_r)$ of $L^2(D\\times D)$.","core_discovery":"The central claim is Theorem 25: for the time-1 annealed transfer operator $L_0$ of a dissipative SDE on $\\mathbb{R}^d$, a local kernel perturbation $\\kappa_\\delta = \\kappa_0 + \\delta\\dot{\\kappa} + r_\\delta$ with $r_\\delta = o(\\delta)$ supported on a compact set $D\\times D$ produces an invariant-density response satisfying $\\|(f_\\delta - f_0)/\\delta - (\\mathrm{Id}-L_0)^{-1}\\dot{L}f_0\\|_1 \\to 0$, where $\\dot{L}f(y)=\\int \\dot{\\kappa}(x,y)f(x)\\,dx$ and $f_0$ is the unperturbed invariant density. The second main claim, Proposition 28, states that for any $\\varphi\\in L^\\infty$ and any closed, bounded, strictly convex set $P\\subset L^2(D\\times D)$ containing zero in its relative interior, the problem of maximising the observable's response over $P$ has a unique solution, given by the normalised Riesz representer of the linear functional $J(\\dot{\\kappa})=\\int \\varphi\\,(\\mathrm{Id}-L_0)^{-1}\\dot{L}f_0$. These results rest on constructing a strong space $B$ of densities with weighted decay at infinity and local $L^2$ control on $D$, on which $L_0$ has a spectral gap. The paper derives a Fourier-approximation scheme for the optimal perturbation and illustrates it on symmetric and asymmetric double-well SDE examples with Gaussian observables.","pith_inferences":["The authors note that existing drift-perturbation estimates converge in $L^2$, which is insufficient for their spectral-gap strategy; if Frechet differentiability of the kernel in the drift were proved in the strong norm, the same optimisation framework would likely cover perturbations that genuinely correspond to SDEs.","Uniqueness requires strict convexity of $P$; for a feasible set like a box the maximum may be attained on a face, although the Riesz-representer formula still gives the steepest direction at the origin.","Computing each Fourier coefficient $G_r$ by solving the resolvent system separately is costly; a fast adjoint method that solves one dual equation per observable would produce all coefficients at once.","A numerical check of the contraction constant on the Ornstein-Uhlenbeck example, where the paper already shows the $L^2$ analogue fails, would directly test whether the spectral-gap mechanism is the right one."],"forward_implications":["The invariant-density derivative has an explicit, numerically evaluable form $(\\mathrm{Id}-L_0)^{-1}\\dot{L}f_0$, so linear response for dissipative SDEs on $\\mathbb{R}^d$ becomes a computational procedure rather than an existence statement.","The optimal kernel-perturbation problem is well-posed: for any closed, bounded, strictly convex set of admissible perturbations containing zero in its relative interior, a unique maximiser exists.","A truncated Fourier basis gives an implementable approximation of the optimal perturbation, with convergence as the basis grows; the experiments show the optimal kernel adjustment is localised where the observable is significant.","The results apply to any kernel transfer operator on $\\mathbb{R}^d$ with the same Gaussian bounds and strong-space properties, so they transfer to other dissipative SDEs beyond the double-well examples.","The perturbed kernel need not correspond to an SDE with a modified drift; the authors frame this local, model-independent perturbation as a first step toward treating drift perturbations in the same framework."],"supporting_citations":[{"why":"Supplies the two-sided Gaussian bounds and gradient estimates for the transition density that give the kernel its regularity and drive the contraction arguments.","marker":"[30]"},{"why":"Provides the weighted L1_alpha spaces, the compact embedding into L1, and the Lasota-Yorke estimates used to get the spectral gap on the noncompact phase space.","marker":"[11]"},{"why":"The Hennion-type theorem cited for quasicompactness, yielding the spectral gap from the Lasota-Yorke inequality.","marker":"[10]"},{"why":"Spectrum-stability theorem used to show the perturbed operators retain a simple leading eigenvalue and that resolvents converge uniformly.","marker":"[25]"},{"why":"Establishes the optimal-response framework for Markov Hilbert-Schmidt integral operators and contains the convex-optimisation propositions reused here.","marker":"[2]"},{"why":"Earlier linear-response framework for stochastic systems whose differentiability results the paper contrasts with its own strong-norm approach.","marker":"[21]"}],"fun_headline_variants":["Unique local kernel tweak optimizes SDE response","Optimal kernel perturbation exists uniquely","Steer SDE observables with optimal kernel shifts","Unique optimal response via kernel perturbation"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The spectral-gap mechanism rests on the assertion, stated without proof in Lemma 14, that the unperturbed transfer operator strictly contracts the strong norm on zero-average densities; the standard positivity argument gives strict contraction only in $L^1$, so if the strong-norm contraction fails, the linear-response and optimisation conclusions would not follow.","fun_headline_variants_meta":{"raw":{"variants":["Unique local kernel tweak optimizes SDE response","Optimal kernel perturbation exists uniquely","Steer SDE observables with optimal kernel shifts","Unique optimal response via kernel perturbation"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000509,"raw_usage":{"total_tokens":2494,"prompt_tokens":977,"completion_tokens":1517,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":593,"completion_tokens_details":{"reasoning_tokens":1474}},"tokens_in":593,"tokens_out":1517,"duration_ms":11489,"temperature":1.0,"reasoning_tokens":1474,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T22:00:06.065718+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"For the one-dimensional Ornstein-Uhlenbeck SDE $dX=-X\\,dt+dW$, compute or bound $\\|L|_{V_B}\\|_s$, the strong-norm operator norm of the transfer operator restricted to zero-average densities. If this norm is not strictly less than 1, the contraction premise behind Lemma 14 is false and the spectral gap, hence linear response and uniqueness of the optimal perturbation, must be established by another mechanism; the paper's own Remark 18 shows the analogous $L^2$ norm does not contract.","supporting_citations":[{"cited_title":"Menozzi, A","cited_arxiv_id":null,"evidence_quote":"Supplies the two-sided Gaussian bounds and gradient estimates for the transition density that give the kernel its regularity and drive the contraction arguments."},{"cited_title":"Extreme Value theory and Poisson statistics for discrete time samplings of stochastic differential equations","cited_arxiv_id":"2310.13972","evidence_quote":"Provides the weighted L1_alpha spaces, the compact embedding into L1, and the Lasota-Yorke estimates used to get the spectral gap on the noncompact phase space."},{"cited_title":"Demers, Niloofar Kiamari, and Carlangelo Liver ani","cited_arxiv_id":null,"evidence_quote":"The Hennion-type theorem cited for quasicompactness, yielding the spectral gap from the Lasota-Yorke inequality."},{"cited_title":"Stability of t he spectrum for transfer operators","cited_arxiv_id":null,"evidence_quote":"Spectrum-stability theorem used to show the perturbed operators retain a simple leading eigenvalue and that resolvents converge uniformly."},{"cited_title":"Optim al linear response for markov hilbert–schmidt integral operators and stochastic dynami cal systems","cited_arxiv_id":null,"evidence_quote":"Establishes the optimal-response framework for Markov Hilbert-Schmidt integral operators and contains the convex-optimisation propositions reused here."},{"cited_title":"A simple framework to j ustify linear response theory","cited_arxiv_id":null,"evidence_quote":"Earlier linear-response framework for stochastic systems whose differentiability results the paper contrasts with its own strong-norm approach."}],"review_version":1}