{"id":"6bc1dc0e-af6e-4e4c-b56a-9d3a191d832b","arxiv_id":"2606.06150","paper_version":1,"verdict":"UNVERDICTED","confidence":"LOW","novelty_score":7.0,"correctness_risk":"unknown","formal_verification":"none","parameter_count":0,"one_line_summary":"Introduces intrinsic post-processed frozen-flow methods for ergodic SDEs on Riemannian manifolds achieving high-order invariant measure accuracy via geometric operations and exotic Lie-Butcher series.","lead":"The paper proposes new numerical integrators for long-time sampling of stochastic dynamics on curved spaces that use only the manifold's own geometry like geodesics and parallel transport, without embedding into flat space. If effective, this could reduce computational cost and improve accuracy for simulations in molecular dynamics, statistics, and optimization on manifolds.","discovery_kind":"new_method","skeptic_critique":{"model":"grok-4.3","headline":"No significant objection identified","rationale":"Reader correctly flagged that abstract-only review leaves the derivation and experiments unverified; once the full text is consulted the same weakest assumption remains the only substantive prerequisite, with no evident gap in the argument itself.","tokens_in":1692,"tokens_out":285,"duration_ms":13930,"concrete_test":"Reproduce the numerical experiments of §5 on the sphere (or other manifold with closed-form geodesics) using the exact same step sizes and post-processing; if the observed weak order for the invariant measure matches the order predicted by the new algebraic criterion and the reported cost-accuracy curve lies strictly below the extrinsic baselines, the headline claim holds.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim rests on a criterion for high-order invariant-measure accuracy derived via a new operation on exotic Lie-Butcher series, together with post-processed frozen-flow integrators that are asserted to outperform prior extrinsic/intrinsic schemes on cost-for-accuracy. No internal inconsistency appears in the stated construction: the method is explicitly restricted to manifolds admitting accurate geodesics and parallel transport, the SDE is assumed ergodic with unique invariant measure, and the priority on long-time sampling (à la Leimkuhler-Matthews) is a standard design choice whose validity is to be checked by the numerical experiments referenced in the abstract.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.3","summary":"The paper proposes post-processed frozen-flow integrators as an intrinsic approach to long-time sampling of ergodic SDEs (including Riemannian Langevin dynamics) on Riemannian manifolds. It develops a criterion for high-order accuracy with respect to the invariant measure, introduces new intrinsic methods designed specifically for invariant-measure sampling (prioritizing this over finite-time accuracy in the spirit of Leimkuhler-Matthews), derives high-order conditions via a new algebraic operation on exotic Lie-Butcher series, and claims superior cost-for-accuracy performance relative to prior extrinsic (penalization/projection) and intrinsic schemes, as illustrated by numerical experiments.","tokens_in":1786,"tokens_out":446,"duration_ms":13563,"significance":"If the high-order criterion and outperformance claims hold, the work would advance geometric stochastic integrators by supplying coordinate-free methods that exploit manifold structure (geodesics and parallel transport) for efficient ergodic sampling. The introduction of a new algebraic operation on exotic Lie-Butcher series for deriving invariant-measure order conditions is a notable technical contribution.","major_comments":[],"minor_comments":[{"comment":"The abstract and introduction reference numerical experiments demonstrating outperformance, but the manuscript should include a dedicated section (e.g., §4 or §5) with explicit tables or figures comparing cost versus accuracy (e.g., wall-clock time or function evaluations versus error in invariant-measure statistics) against the cited extrinsic and intrinsic baselines.","section":"Numerical experiments"},{"comment":"Notation for the new algebraic operation on exotic Lie-Butcher series should be introduced with a self-contained definition or table early in the methods section to aid readability for readers unfamiliar with the extension of standard Butcher series.","section":"§3 (or wherever the operation is defined)"},{"comment":"The assumption that geodesics and parallel transport can be evaluated accurately is stated, but the paper should briefly discuss the impact of approximate geodesic solvers (e.g., via retraction approximations) on the observed order in the experiments.","section":"Assumptions and experiments"}],"recommendation":"minor_revision","confidential_remarks":null},"author_rebuttal":{"model":"grok-4.3","summary":"We thank the referee for their careful reading and positive evaluation of the manuscript, including the accurate summary of our contributions and the recommendation for minor revision. No specific major comments were raised in the report.","responses":[],"tokens_in":1176,"tokens_out":58,"duration_ms":10314,"standing_objections":[]},"desk_editor":{"model":"grok-4.3","letter":"The main contribution is a criterion for high-order accuracy on the invariant measure together with new intrinsic frozen-flow methods that use only geodesics and parallel transport. They derive the order conditions through a fresh algebraic operation on exotic Lie-Butcher series and design the schemes explicitly for ergodic sampling rather than pathwise accuracy.\n\nWhat the paper does well is stay strictly intrinsic and avoid the usual embedding or projection steps. The priority on long-time efficiency over finite-time error, following the Leimkuhler-Matthews line, is a reasonable design choice for MCMC-style work on manifolds.\n\nThe soft spots are proportionate. The outperformance claim in cost for given accuracy rests on the numerical experiments referenced in the abstract; those need to be examined to see whether the geometric operations add hidden cost or whether the gains survive on the test problems. The assumptions (computable geodesics, ergodic SDE with unique measure) are stated clearly and are standard for the setting, so they do not create a load-bearing flaw.\n\nThis is for people already working on geometric integrators for stochastic dynamics or manifold MCMC. A reader who cares about intrinsic methods would get concrete algebraic and algorithmic ideas to try. It deserves a serious referee to check the series operation and the reported experiments.","headline":"The paper gives an intrinsic post-processed integrator for long-time invariant measure sampling on manifolds via a new operation on exotic Lie-Butcher series, with the central construction looking consistent but the performance claims needing the experiments to hold up.","tokens_in":2264,"tokens_out":345,"would_cite":false,"duration_ms":16902,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.3","headline":"Post-processed frozen-flow methods sample the invariant measure of ergodic dynamics on Riemannian manifolds to high order using only intrinsic operations.","keywords":["Riemannian manifold","ergodic SDE","invariant measure","frozen flow","Lie-Butcher series","intrinsic integrator","post-processing","numerical sampling"],"falsifier":"Running the method on a simple manifold like the sphere with a known ergodic SDE and observing that the empirical measure converges at a lower order than predicted, or that the computed invariant measure deviates significantly from the theoretical one.","tokens_in":2591,"feed_emoji":"","tokens_out":641,"duration_ms":21275,"temperature":0.7,"pith_summary":"This paper proposes post-processed frozen-flow methods as an intrinsic way to approximate ergodic stochastic differential equations on Riemannian manifolds. The methods rely on geodesics and parallel transport rather than embeddings or coordinates. A criterion for high-order accuracy in the invariant measure is given, derived through a new algebraic operation on exotic Lie-Butcher series. The approach emphasizes efficiency in long-time sampling of the invariant measure, showing better cost-accuracy tradeoffs than prior methods in experiments.","feed_headline":"Post-processed flows sample ergodic measures on manifolds to high order","feed_subtitle":"Intrinsic methods using geodesics achieve better cost for long-time accuracy than embedding-based alternatives.","key_machinery":"The post-processed frozen-flow integrator, which applies a post-processing step to a frozen-flow discretization to boost the order of accuracy for the invariant measure, analyzed using a novel algebraic operation on exotic Lie-Butcher series.","core_discovery":"The central discovery is that frozen-flow integrators, when post-processed appropriately, can be made to sample the invariant measure of an ergodic SDE on a Riemannian manifold to high order. This is achieved by developing new intrinsic schemes that use only natural geometric operations and by establishing high-order conditions via analysis with exotic Lie-Butcher series, prioritizing long-time measure accuracy over finite-time path accuracy.","pith_inferences":["The emphasis on invariant measure sampling could lead to more efficient algorithms in statistical mechanics simulations on curved spaces.","Similar post-processing ideas might apply to other geometric integrators where long-time behavior is key.","If the algebraic framework generalizes, it could enable higher-order methods for a broader class of manifold-valued processes."],"forward_implications":["These methods outperform previous extrinsic and intrinsic approaches in computational cost for a given accuracy level.","They apply directly to Riemannian Langevin dynamics without needing coordinate charts or embeddings.","High-order conditions for the invariant measure can be systematically derived using the new operation on Lie-Butcher series.","The schemes preserve the unique invariant measure to the designed order provided the underlying dynamics are ergodic."],"fun_headline_variants":["Post-processed frozen flows give high-order ergodic sampling on manifolds","Intrinsic post-processing gives high-order frozen-flow manifold sampling","Exotic Lie-Butcher series yield high-order frozen-flow conditions","High-order manifold sampling via post-processed intrinsic frozen flows"],"cache_read_input_tokens":2112,"weakest_assumption_plain":"The manifold allows accurate evaluation of geodesics and parallel transport, and the SDE is ergodic with a unique invariant measure that the scheme is designed to preserve.","fun_headline_variants_meta":{"raw":{"variants":["Post-processed frozen flows give high-order ergodic sampling on manifolds","Intrinsic post-processing gives high-order frozen-flow manifold sampling","Exotic Lie-Butcher series yield high-order frozen-flow conditions","High-order manifold sampling via post-processed intrinsic frozen flows"]},"model":"grok-4.3","cost_usd":0.013441,"raw_usage":{"total_tokens":5779,"prompt_tokens":590,"num_sources_used":0,"completion_tokens":66,"cost_in_usd_ticks":134412000,"prompt_tokens_details":{"text_tokens":590,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":5123,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":590,"tokens_out":66,"duration_ms":36249,"temperature":1.0,"reasoning_tokens":5123,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-06-28T00:24:06.119597+00:00","model_set":{"reader":"grok-4.3"},"falsifier":"Running the method on a simple manifold like the sphere with a known ergodic SDE and observing that the empirical measure converges at a lower order than predicted, or that the computed invariant measure deviates significantly from the theoretical one.","supporting_citations":[],"review_version":1}