{"id":"58fa236e-8304-451c-8cec-b9bfb05aabb9","arxiv_id":"2505.20217","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"For one-dimensional drift-diffusion equations with small noise, the long-time limit of the solution at each critical time scale equals a Markov chain on a hierarchical set of metastable wells with explicit jump rates.","lead":"This paper proves a precise rule for how solutions of a one-dimensional diffusion equation with small noise change over very long times, moving among stable states through a hierarchy of slow time scales. It gives explicit limiting formulas at every scale and could help predict rare transitions in systems modeled by one-dimensional random motions.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The derivation of the higher-level jump rates (2.22) rests on an unproved Lemma 6.4, and the multi-saddle case is not a routine analogue of Lemma 5.5.","rationale":"The paper's central theorem is plausible and the resolvent strategy is coherent, but the proof of the higher-level resolvent convergence contains an articulated gap: Lemma 6.4 is stated without proof, with the explanation that it is similar to Lemma 5.5. The reader's verdict already identified Lemma 6.4 as a missing piece, but the formal weakest_assumption in the reader's report is the nondegeneracy of the drift zeros. Nondegeneracy is an explicit hypothesis of the theorem and not a correctness risk, so the most load-bearing unresolved point is the omitted proof of Lemma 6.4 in the multi-saddle situation. My agreement with the reader is therefore partial: same overall verdict, but a different emphasis on where the main risk lies. If Lemma 6.4 fails, Lemma 6.3 fails, and the effective rates (2.22) for the Markov chain at every level p >= 2 would not follow, breaking the hierarchy underlying Theorems 2.5 and 3.2. If the missing estimate is supplied with a o(1) bound, the central claim appears sound. No evidence suggests the claim is false; the issue is completeness of the proof, consistent with a conditional acceptance.","tokens_in":60569,"tokens_out":21241,"duration_ms":215662,"concrete_test":"Supply the missing proof of Lemma 6.4 for the case W^{(p)}_{k,k+1} contains two distinct saddles sigma^- < sigma^+. Specifically, estimate the integral over B_r(sigma^-) by an integration by parts and use Lemma 6.6 to bound the oscillation of phi_{p,epsilon} on that binding set in terms of |f(M_p(k+1)) - f(M_p(k))| and eta/delta = 1/(2 log epsilon^{-1}). If the resulting bound is o(1), Lemma 6.4 holds unchanged; if it leaves a non-vanishing term, the jump rates in (2.22) must be corrected and Theorem 2.5(a) for p >= 2 fails as stated.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Lemma 6.4 is needed in Lemma 6.3 to show that the binding-set contribution to the resolvent identity vanishes for every p > 1. The text says the proof is omitted because it is similar to Lemma 5.5. That similarity is not automatic. Lemma 5.5 uses Proposition 5.2, which gives local constancy of the resolvent on the entire interval containing the binding set. For p > 1, when the level-p barrier between M_p(k) and M_p(k+1) has two distinct saddles, sigma^- < sigma^+, the right binding set B_r(sigma^-) and the left binding set B_l(sigma^+) lie in the region between the saddles. There the resolvent is not constant; Lemma 6.6 shows it interpolates between f(M_p(k)) and f(M_p(k+1)) through the weight varpi^+_{p,k}. Controlling the integral over these binding sets therefore requires an estimate on the oscillation of the resolvent on intervals of length eta = epsilon/delta where the function changes by an O(1) amount over a distance comparable to delta. This is a genuinely different argument from Lemma 5.5. Until that argument is written and the integral is proved to vanish, the derivation of the rates (2.22) and hence Theorem 2.5(a) for p >= 2 is incomplete. The nondegeneracy of the zeros of b is an explicit assumption, not a hidden gap; the missing lemma is the real soft spot.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the one-dimensional parabolic equation ∂_t u_ε = b(x) ∂_x u_ε + ε a(x) ∂_x^2 u_ε on R, where b and a are 1-periodic C^1 functions with Lipschitz derivatives, a is bounded below by a positive constant, and every zero of b is simple. The authors construct, for each level p = 1,...,q, a time scale θ_ε^(p), a partition of the stable equilibria into sets M_p(k), and an explicit continuous-time Markov chain X_p with rates given by (2.21)–(2.22). Their main theorem, Theorem 2.5(a), states that u_ε(θ_ε^(p) t, x) converges to a convolution of an explicit hitting distribution h_p(x,·), the semigroup of X_p, and the π-weighted average of u_0 over M_p(ℓ), with no fitted constants. Theorems 2.5(b) and 2.5(c) characterize intermediate time scales, and Theorem 2.8 treats the final layer for periodic initial data. The proof uses the stochastic representation of solutions and the resolvent approach to metastability: Theorem 3.3 shows that resolvent solutions are asymptotically constant on wells and that all limit points solve the reduced resolvent equation; this is then converted into convergence of finite-dimensional distributions (Theorem 3.2) and, via Lemma 8.2, into the parabolic limit.","tokens_in":60852,"tokens_out":7762,"duration_ms":90393,"significance":"If the proof is completed, this is a complete and explicit description of metastability for a general one-dimensional non-reversible periodic diffusion at every critical time scale. It extends earlier results that assumed constant a or reversible Gibbs structure, and it provides parameter-free formulas for the renormalized jump rates and for the hitting and invariant weights. The resolvent framework is applied systematically, and the hierarchical construction is detailed and checkable. The main weakness is that one load-bearing estimate for p ≥ 2 is asserted without proof, which leaves the higher-level rates and hence Theorem 2.5(a) for p ≥ 2 incomplete as it stands.","major_comments":[{"comment":"Lemma 6.4 is load-bearing: it is used in Lemma 6.3 to show that the binding-set contribution to the resolvent identity (6.9) vanishes, and Lemma 6.3 is what produces the reduced generator L_p in the p > 1 case. The proof is omitted as 'similar to Lemma 5.5', but that analogy is not automatic. In the p = 1 proof, Proposition 5.2 gives local constancy of φ_{1,ε} on the whole interval containing the binding set B_ε(k), and the estimate in Lemma 5.5 uses this constancy on the entire binding interval. For p ≥ 2, when the level-p barrier between M_p(k) and M_p(k+1) has two distinct saddles σ^- < σ^+, the right binding set B_r(σ^-) and the left binding set B_l(σ^+) lie between the two saddles. There the resolvent is not constant: Lemma 6.6 shows that φ_{p,ε} interpolates between f_ε(M_p(k)) and f_ε(M_p(k+1)) through the weight ϖ^+_{p,k} over a distance of order δ = sqrt(ε log ε^{-1}), while the binding intervals have length η = ε/δ. Controlling the integral therefore requires an oscillation estimate on intervals of length η in a region where the resolvent changes by O(1) over distance δ, which is genuinely different from the argument of Lemma 5.5. Until this estimate is supplied, the derivation of the rates (2.22) and the conclusion of Theorem 2.5(a) for p ≥ 2 are incomplete.","section":"§6, Lemma 6.4"}],"minor_comments":[{"comment":"Part (c) is a statement in the main theorem, but its proof ends with 'The details are left to the reader'. The argument may indeed be a repetition of the proof of part (b) with θ^{(p)}_ε replaced by ε^{-1}, but the details should be written out or the statement should be presented as a corollary with a precise derivation.","section":"§8, proof of Theorem 2.5(c)"},{"comment":"The estimate E^ε_x[τ(a,b)]/θ^{(p)}_ε = o_ε(1) is asserted by reference to Lemmas 4.2 and A.11. Since the interval [a,b] contains two distinct saddles and possibly many local minima, the reader would benefit from a displayed verification that the maximum of S(x)-S(y) over the relevant interval is strictly below h_p by a uniform (ε-independent) margin.","section":"§6, proof of Lemma 6.6"},{"comment":"The superscript placement in the quantities h_{p+1,+}^k and h_{p+1,k}^+ is inconsistent between the displayed rates (2.22) and the surrounding text; the notation should be unified.","section":"§2, equations (2.21)–(2.22)"},{"comment":"The sentence 'By Corollary 4.9, h_p(x,M_p(j)) = P_x[τ(M_p(j))=τ(M_p)] + o_ε(1)' says the error is uniform in the parameters on which it depends; it would be clearer to state explicitly over which parameters (x, j, p, k) the uniformity holds.","section":"§3, Remark 3.4"}],"recommendation":"major_revision","confidential_remarks":"The missing proof of Lemma 6.4 is the only substantive obstacle I see to the paper's central claims. It is likely repairable within the manuscript's scope, but it cannot be waved away as routine in view of the multi-saddle interpolation described by Lemma 6.6. I would ask the authors to provide the omitted argument in full. The paper is otherwise in scope for a probability journal and makes a significant contribution."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read it. The main theorem is new and the proof strategy is coherent: the authors extend the resolvent approach to general periodic non-reversible drifts, where no stationary measure is known, and obtain the entire hierarchy of critical and intermediate scales with explicit weights and no fitted constants. That is a real step beyond the reversible-gradient cases in [18,19] and the constant-diffusion case in [23]. The stochastic representation linking the parabolic solution to the diffusion, the explicit formulas for h_p and the jump rates, and the inductive construction of the hierarchy are all clearly laid out. The paper earns its place in the metastability literature.\n\nThe soft spots are real but localized. Lemma 6.4 is asserted without proof, and the stress-test note is right that it is not a routine analogue of Lemma 5.5. For p > 1 with two distinct saddles, the binding sets lie in the region where the resolvent interpolates between the two adjacent wells, so controlling the integral requires an oscillation estimate on intervals of length eta = epsilon/delta, not the local-constancy argument used at level one. Until that argument is supplied, the derivation of the higher-level rates (2.22) and hence Theorem 2.5(a) for p >= 2 rests on an unproved step. The proof of Theorem 2.5(c) also leaves details to the reader, but that is minor. The nondegeneracy assumption b'(x_j) != 0 is explicit and standard in this context; it is not a hidden flaw.\n\nMy overall assessment: the central argument is plausible and likely fillable, but the paper as posted is conditional. I would send it to a serious referee rather than desk reject it, and the referee should ask for a complete proof of Lemma 6.4 before accepting. The intended audience is researchers in metastability and small-noise stochastic dynamics; for them the paper is worth reading now.","headline":"Full multi-scale metastable asymptotics for a natural class of 1D parabolic equations, with one genuinely missing lemma in the higher-level rate derivation.","tokens_in":61365,"tokens_out":1383,"would_cite":true,"duration_ms":21264,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["60J60","82C31"],"pacs":[],"model":"deepseek-v4-flash","headline":"A complete Markov-chain hierarchy describes all long-time asymptotics of a one-dimensional diffusion PDE.","keywords":["metastability","parabolic differential equations","small-noise diffusion","Markov chain hierarchy","time scales","resolvent approach","Eyring-Kramers rates","periodic drift"],"falsifier":"Compute, numerically or analytically, the first critical-scale limit for a periodic drift with a simple saddle, for instance $b(x)=\\sin(2\\pi x)$, $a(x)\\equiv 1$, where the theorem predicts $\\theta_\\epsilon^{(1)}=e^{1/(\\pi\\epsilon)}$ and a nearest-neighbour symmetric random walk with jump rate 1; a mismatch in the prefactor or rate would disprove the formula.","tokens_in":60364,"feed_emoji":"⏳","tokens_out":11131,"duration_ms":120707,"temperature":0.7,"pith_summary":"Consider the parabolic equation $\\partial_t u_\\epsilon = b(x)\\partial_x u_\\epsilon + \\epsilon a(x)\\partial_x^2 u_\\epsilon$ with periodic drift $b$ and strictly positive diffusion $a$, started from a bounded continuous $u_0$. The paper proves that as $\\epsilon\\to 0$ the solution evolves in a finite cascade of exponentially separated time scales $\\theta_\\epsilon^{(1)}\\ll\\cdots\\ll\\theta_\\epsilon^{(q)}$, and that on each of these scales the limit of $u_\\epsilon(t\\theta_\\epsilon^{(p)},x)$ is the expectation of $u_0$ under an explicit continuous-time Markov chain $X_p$ on aggregates of stable equilibria of $\\dot x=b(x)$. All data of the limit — the time scales, the chain's jump rates, the entrance probabilities $h_p(x,\\cdot)$, and the local weights $\\pi$ — are computed from Gaussian integrals around the critical points of the potential $S(x)=-\\int b/a$, with no fitted constants. This matters because it turns a singularly perturbed PDE into a finite, explicit hierarchy of Markov-chain reductions for arbitrary periodic one-dimensional coefficients with simple zeros of $b$.","feed_headline":"Markov chains predict every long-time scale of a 1D diffusion PDE","feed_subtitle":"Tiny-noise heat equations relax in layers, each governed by explicit jump rates computed from the drift and diffusion.","key_machinery":"The carrying object is the resolvent approach to metastability. For each level $p$ and $\\lambda>0$ one solves the resolvent equation $(\\lambda-\\theta_\\epsilon^{(p)} L_\\epsilon)\\phi_{p,\\epsilon}=G$ for a function $G$ constant on each well $E(M_p(k))$; Theorem 3.3 shows $\\phi_{p,\\epsilon}$ becomes constant on each well as $\\epsilon\\to0$ and that its limit solves the reduced resolvent equation $(\\lambda-L_p)f=g$ for the Markov chain $X_p$. This step is powered by local ergodicity estimates, hitting-time bounds, and Gaussian saddle-point expansions near unstable equilibria, encoded in test functions built from the linearized generator around each saddle. From the resolvent asymptotics one deduces convergence of the finite-dimensional distributions of the sped-up diffusion $X_\\epsilon(\\theta_\\epsilon^{(p)}\\cdot)$ to $X_p$ (Theorem 3.2), and the stochastic representation $u_\\epsilon(t,x)=E_x^\\epsilon[u_0(X_\\epsilon(t))]$ transfers that convergence to the solution of the parabolic equation.","core_discovery":"The central assertion is Theorem 2.5(a): for each level $1\\le p\\le q$, with $\\theta_\\epsilon^{(p)}=e^{h_p/\\epsilon}$, $$\\lim_{\\epsilon\\to0} u_\\epsilon(\\theta_\\$epsilon^{{(p)}}$ t, x) = \\sum_{k\\in\\mathbb{Z}} h_p(x,M_p(k)) \\sum_{\\ell\\in\\mathbb{Z}} $p_t^{{(p)}}$(M_p(k),M_p(\\ell)) \\sum_{m'\\in M_p(\\ell)} \\frac{\\pi(m')}{\\pi(M_p(\\ell))} u_0(m').$$ Here the sets $M_p(k)$ form a partition of the stable equilibria $\\mathcal{M}$ of $\\dot x=b(x)$, obtained by merging the closed irreducible classes of the previous Markov chain $X_{p-1}$; $h_p(x,\\cdot)$ is the probability that the diffusion started at $x$ first hits the aggregate $M_p(k)$; and $\\pi(m)=\\sqrt{-2\\pi/(b'(m)a(m))}$ at a stable equilibrium $m$, with $\\pi(M)=\\sum_{m\\in M}\\pi(m)$. The jump rates of $X_p$ are $R_p(M_p(k),M_p(k\\pm1))=\\pi_p(k)^{-1}\\sigma_p(k,k\\pm1)\\mathbf{1}\\{h_k^{p,\\pm}=h_p\\}$, where $\\sigma_p$ sums Gaussian saddle weights $\\sqrt{2\\pi a(\\sigma)/b'(\\sigma)}$ over the highest points between adjacent aggregates. Consequently every quantity in the limit — the exponential time scales, the Markov-chain rates, the entrance laws, and the local stationary weights — is computed directly from the coefficients $a$ and $b$, with no parameter fitting.","pith_inferences":["Although the paper works on the line with periodic coefficients, its estimates are local, so the same ratio-of-saddle-weights structure should describe non-periodic one-dimensional landscapes with finitely or countably many wells, with periodicity mainly simplifying the final layer and normalisation.","A direct finite-$\\epsilon$ numerical comparison of the PDE solution against the Markov-chain prediction — including the prefactors — would give quantitative evidence for the rate of convergence to the limit, a question the paper does not address.","The degenerate case $b'(x_j)=0$ is the natural next target: one expects modified saddle-point weights, possibly of Gamma-function type for a quartic saddle, and possibly slower-than-exponential corrections to the escape scale; the paper's own remarks point in this direction."],"forward_implications":["At every critical time scale $\\theta_\\epsilon^{(p)}$, the parabolic solution's limit is the expectation of $u_0$ under the Markov chain $X_p$ with explicit rates, so no fitted constants or unknown stationary measures enter.","Between two consecutive critical scales, the solution relaxes to the $\\pi$-weighted average over the next aggregate $M_{p+1}(k)$; in particular, starting exactly at a saddle point gives the one-half average of the two adjacent wells, as stated in Theorem 2.5(c).","The hierarchy is finite — the number of aggregates strictly decreases at each layer — so the full long-time behaviour is described by at most $N$ Markov-chain reductions, and the final-layer limit on the torus is reversible with the explicit weights when $S(0)=S(1)$.","The diffusion's own trajectory at time $\\theta_\\epsilon^{(p)} t$ converges in finite-dimensional law to $X_p$ started from the entrance distribution $h_p(x,\\cdot)$, so the PDE limit is pointwise sharp for all $x$ and $t$, including points at saddle points."],"supporting_citations":[{"why":"Supplies the resolvent approach: resolvent solutions become asymptotically constant on wells, yielding convergence of the order process to the reduced Markov chain.","marker":"[21]"},{"why":"Provides the route from convergence of the trace to convergence of the finite-dimensional distributions of the speeded-up process.","marker":"[20]"},{"why":"Gives the stochastic representation of parabolic solutions used to transfer the diffusion convergence to $u_\\epsilon$.","marker":"[12]"},{"why":"Provides existence, uniqueness, and elliptic regularity for the weak solutions of the resolvent equations.","marker":"[7]"}],"fun_headline_variants":["Metastable cascade in 1D heat equation: all time scales via Markov chains","Small-noise heat PDEs relax in layers, each controlled by explicit jump rates","From drift and diffusion to full metastable hierarchy of a 1D diffusion","All long-time asymptotics of a 1D diffusion, computed from coefficients","Multiscale limit of parabolic equations: explicit Markov chain structure"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that every zero of the periodic drift $b$ is simple, $b'(x_j)\\neq 0$, with $a>0$ of class $C^1$ and Lipschitz derivative; if a critical point of the potential $S(x)=-\\int b/a$ were degenerate, the Gaussian prefactors and exponential time scales would no longer be the correct asymptotics.","fun_headline_variants_meta":{"raw":{"variants":["Metastable cascade in 1D heat equation: all time scales via Markov chains","Small-noise heat PDEs relax in layers, each controlled by explicit jump rates","From drift and diffusion to full metastable hierarchy of a 1D diffusion","All long-time asymptotics of a 1D diffusion, computed from coefficients","Multiscale limit of parabolic equations: explicit Markov chain structure"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000319,"raw_usage":{"total_tokens":2034,"prompt_tokens":1411,"completion_tokens":623,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":1027,"completion_tokens_details":{"reasoning_tokens":523}},"tokens_in":1027,"tokens_out":623,"duration_ms":7623,"temperature":1.0,"reasoning_tokens":523,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T13:57:17.906187+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute, numerically or analytically, the first critical-scale limit for a periodic drift with a simple saddle, for instance $b(x)=\\sin(2\\pi x)$, $a(x)\\equiv 1$, where the theorem predicts $\\theta_\\epsilon^{(1)}=e^{1/(\\pi\\epsilon)}$ and a nearest-neighbour symmetric random walk with jump rate 1; a mismatch in the prefactor or rate would disprove the formula.","supporting_citations":[{"cited_title":"Landim and D","cited_arxiv_id":null,"evidence_quote":"Supplies the resolvent approach: resolvent solutions become asymptotically constant on wells, yielding convergence of the order process to the reduced Markov chain."},{"cited_title":"Landim, M","cited_arxiv_id":null,"evidence_quote":"Provides the route from convergence of the trace to convergence of the finite-dimensional distributions of the speeded-up process."},{"cited_title":"Friedman:Stochastic differential equations and applications","cited_arxiv_id":null,"evidence_quote":"Gives the stochastic representation of parabolic solutions used to transfer the diffusion convergence to $u_\\epsilon$."},{"cited_title":"Evans: Partial Differential Equations","cited_arxiv_id":null,"evidence_quote":"Provides existence, uniqueness, and elliptic regularity for the weak solutions of the resolvent equations."}],"review_version":1}