Back to search
2102.07760

Learning forecasts of rare stratospheric transitions from short simulations

Justin Finkel, Robert J. Webber, Edwin P. Gerber, Dorian S. Abbot, Jonathan Weare

correctmedium confidence
Category
Not specified
Journal tier
Strong Field
Processed
Sep 28, 2025, 12:56 AM

Audit review

The paper’s Appendix B derives (i) the Feynman–Kac boundary value problem for F defined in Eq. (B6) and its PDE/BC form (B13) via an Itô/martingale argument, (ii) the finite-time “stopped” operator identity (T^Δt_{D^c} − 1)F = I^Δt_{D^c}[VF − H] (B19), (iii) specializations to the committor and MFPT in the main text (13)–(15), and (iv) the homogenized Petrov–Galerkin system (B20)–(B22), together with Monte Carlo inner products (B23)–(B28) . The candidate solution reproduces items (i)–(iii) correctly and matches the operator/Galerkin framing, but it introduces a sign error in the RHS of the homogenized finite-time equation and hence the Galerkin right-hand side: the paper’s (B21)–(B22) gives −(T−1)ĝ − I[H − Vĝ], whereas the candidate writes −(T−I)ĝ + I[H − Vĝ], i.e., the integral term has the wrong sign. This propagates to their Monte Carlo estimator for b, where the integral should enter with a minus sign, not a plus .

Referee report (LaTeX)

\textbf{Recommendation:} minor revisions

\textbf{Journal Tier:} strong field

\textbf{Justification:}

Technically solid and well-motivated: the paper correctly derives the Feynman–Kac BVP (Appendix B, B6–B13) and the finite-time stopped operator identity with its Galerkin discretization (B19–B22), and connects these to committor and MFPT formulae in the main text (13)–(15). Minor improvements would be to make assumptions explicit and to maintain strict sign consistency when moving terms across equalities, which aids reproducibility in code implementations. The overall contribution is valuable for computing rare-event statistics from short-trajectory data in high-dimensional models.