{"id":"b8a44d7e-21ed-43d7-9411-4db07540f6b5","arxiv_id":"2412.11559","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"For Riemann-Liouville fractional Brownian motion, the solution of the FLEFE, the time-averaged MSD converges to the mean-squared increment rather than the MSD for 1/2<α<3/2, yielding spurious nonergodicity, while strong aging restores ergodicity.","lead":"This paper derives exact formulas for how a fractional Langevin model far from equilibrium behaves in single-trajectory and ensemble measurements. It shows that apparent nonergodicity in a wide parameter range is an artifact of comparing the wrong averages, and that waiting before measuring can restore ergodicity.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The universal claim fails at α=1: ordinary Brownian motion lies inside 1/2<α<3/2, and there Eqs. (27) and (43) both give 2K₁Δ, so TAMSD=MSD and no spurious nonergodicity occurs; all 'any α>1/2' statements need the caveat α≠1.","rationale":"I read the paper as a largely self-contained derivation of MSD, MSI, TAMSD, higher-order increments, and aging for RL-FBM/FLEFE. The main asymptotic results for fractional α appear internally consistent: for 1/2<α<3/2 the mean TAMSD is the time average of the MSI, and once the MSI becomes asymptotically stationary the TAMSD must share its prefactor, so Eq. (43) follows from Eq. (28) plus the standard beta-integral identity in Eq. (32). For α>3/2 and α=3/2 the stated power-law and logarithmic forms also follow by direct integration of the MSI. I therefore do not think the H-function expansion in Appendix B is the most load-bearing point: it is one route to the same asymptotic result, and the key cancellations are already visible in the MSI integration. The more concrete and falsifiable problem is the universal claim itself. α=1 lies exactly in the claimed spurious-nonergodicity window, and there the formulas coincide, so the headline assertion is false without a caveat. This is a genuine correctness issue in the central claim, but it is local and easily repaired by excluding α=1. The absence of simulation error bars and code is a reproducibility weakness, but it does not affect the analytical counterexample at α=1. Hence I recommend keeping a conditional verdict, with the condition that all universal statements be corrected to exclude the Brownian limit.","tokens_in":29968,"tokens_out":23985,"duration_ms":218055,"concrete_test":"Plot or evaluate the ratio R(α)=⟨δ²(Δ)⟩/⟨x²(Δ)⟩ from Eqs. (27) and (43) for α∈(1/2,3/2), e.g. α=0.6, 0.8, 1.0, 1.2, 1.4, with any common Kα. At α=1 the ratio is exactly 1, while at the other listed values it differs from 1. Equivalently, re-check the Appendix B calculation at α=1, which already yields I₁(Δ,T)=2K₁(T−Δ) and hence TAMSD=MSD=2K₁Δ. If this is confirmed, revise the abstract and Sec. IV.C claims 'for any α>1/2' and 'in the entire domain' to 'for α>1/2, α≠1' or add an explicit exception for the Brownian limit.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The paper's stated central claim is that the mean TAMSD and MSD do not coincide for any α>1/2, and that spurious nonergodicity occurs in the whole interval 1/2<α<3/2 (abstract and Sec. IV.C). This is false at α=1, which is inside that interval. At α=1, RL-FBM reduces to ordinary Brownian motion. Substituting α=1 into Eq. (27) gives ⟨x²(Δ)⟩=2K₁Δ, and substituting into Eq. (43) gives ⟨δ²(Δ)⟩=2K₁Δ, because Γ(2)=1 and |cos π|=1. The paper's own Appendix B treats α=1 as a special integer case and obtains I₁(Δ,T)=2K₁(T−Δ), which yields exactly the Brownian TAMSD 2K₁Δ. Thus the TAMSD/MSD ratio is 1, not spurious nonergodicity. The prefactor ratio R(α)=(2α−1)Γ(α)²/[Γ(2α)|cos πα|] equals 1 only at α=1 in this regime. This does not invalidate the fractional-order results, but the universal quantifier 'any α>1/2' and the phrase 'in the entire domain' are incorrect as written; the correct statement must exclude α=1 or explicitly identify it as the ordinary ergodic Brownian limit.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper analyzes the fractional Langevin equation far from equilibrium (FLEFE), whose zero-initial-condition solution is Riemann-Liouville fractional Brownian motion (RL-FBM). The authors derive exact and asymptotic expressions for the MSD, the mean-squared increment (MSI), the time-averaged MSD (TAMSD), higher-order increments, and aging effects. The central claim is that for 1/2<α<3/2 the mean TAMSD does not converge to the MSD but instead to the stationary MSI, leading to ``spurious nonergodicity,'' while for α≥3/2 the first increments are nonergodic and higher-order increments restore stationarity in appropriate intervals of α. Strong aging is claimed to restore ergodicity in 1/2<α<3/2. The main results are derived via hypergeometric and H-function expansions and are compared with numerical simulations.","tokens_in":30338,"tokens_out":11174,"duration_ms":105059,"significance":"If the results are properly qualified, the paper is a valuable contribution: it provides closed-form expressions for central observables of RL-FBM with no fitted parameters, it draws attention to the distinction between MSD and MSI for processes with nonstationary increments, and it offers testable predictions for single-particle-tracking data. The higher-order increment stationarity and aging analysis extend the known properties of RL-FBM. The paper is self-contained, and the simulation results support the main asymptotic formulas. However, the headline universal claim fails at α=1, which is inside the stated spurious-nonergodicity interval, and the higher-order increment derivation relies on an unproved derivative approximation; both points need to be addressed before the claims can be accepted as stated.","major_comments":[{"comment":"The statement that the TAMSD and MSD differ for every α>1/2, and that spurious nonergodicity occurs throughout 1/2<α<3/2, is false at α=1. At α=1 the FLEFE reduces to ordinary Brownian motion: Eq. (27) gives ⟨x²(t)⟩=2K₁t, while Eq. (33) and Eq. (43) both give 2K₁Δ, so the TAMSD equals the MSD and the ratio is 1. This is confirmed by the special case in Appendix B: Eq. (B22) yields I₁=2K₁T(1−Δ/T), which together with Eq. (B1) gives ⟨δ²(Δ)⟩=2K₁Δ. The universal quantifiers ``any α>1/2'' and ``the entire domain'' in the abstract, Sec. IV.C, and the conclusions must therefore be qualified by excluding α=1, or α=1 must be explicitly identified as the ordinary ergodic Brownian limit.","section":"Abstract; Sec. IV.C; Eq. (46)"},{"comment":"The derivation of higher-order increment stationarity rests on the approximation Δ^{(1)}x(t;τ₁) ≈ τ₁ dx(t)/dt for t≫τ₁, stated in Eq. (50), without an error estimate. Since Eq. (53), the general higher-order increment result, is a central claim, the authors should justify this replacement in the mean-square sense, for example by showing that the neglected terms are subdominant in the appropriate limit, or by deriving the second-order increment result directly from the stochastic integral representation. The numerical agreement in Fig. 3 is encouraging, but it does not replace a bound or a rigorous statement of the approximation's validity regime.","section":"Sec. V; Eq. (50)"}],"minor_comments":[{"comment":"The title contains a typo: ``Rieman n-Liouville'' should be ``Riemann-Liouville''.","section":"Title"},{"comment":"The word ``revels'' should be ``reveals'' in the sentence describing the stationarity of the second-order increment.","section":"Sec. V, text after Eq. (52)"},{"comment":"The text contains the typo ``TMASD'' where ``TAMSD'' is meant; this should be corrected.","section":"Appendix B"},{"comment":"The caption says the mean TAMSD converges to the MSI ``in the strong aging limit,'' but the surrounding text concerns the long-time limit T/Δ≫1; please clarify whether aging is involved in this figure.","section":"Fig. 2 caption"},{"comment":"The sentence ``the relation δ ≫ δ is not perfectly fulfilled'' appears to contain a typo; presumably it should read ``δ ≪ Δ'' or ``δ ≪ Δ is not perfectly fulfilled.'","section":"Appendix D"}],"recommendation":"major_revision","confidential_remarks":"The α=1 counterexample is a genuine flaw in the headline claim, but it is localized and fixable by an explicit caveat; I do not see grounds for rejection. The higher-order increment derivation should also be strengthened as requested. No concerns about citation practices or novelty disclosure came up in my reading."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a useful, largely self-contained derivation of the TAMSD for Riemann-Liouville FBM (the FLEFE solution), covering all fractional orders and adding a clean aging analysis. The core result—that for 1/2<α<3/2 the mean TAMSD converges to the MSI rather than the MSD—is new and the simulations back it up. That said, the universal statements 'any α>1/2' and 'in the entire domain' are literally false at α=1, where the process is ordinary Brownian motion and Eq. (27) and Eq. (43) both give 2K₁Δ, so TAMSD=MSD and there is no spurious nonergodicity. The paper's own Appendix B treats α=1 separately and produces exactly the Brownian TAMSD, so this is an oversight in the wording, not a flaw in the fractional-order math. It should be fixed by excluding α=1 (or explicitly labeling it the ergodic Brownian limit) throughout the abstract and conclusions.\n\nWhat's new: the TAMSD for all α>1/2 (Eqs. 43-46), the higher-order increment stationarity (Eqs. 52-53), and the strong-aging ergodicity restoration. The exact expressions for MSI and MSD were already in Lim and Eab-Lim, but the TAMSD is new, as is the distinction between MSD and MSI for nonstationary increments. The derivations are careful and the H-function machinery is used competently, even if the asymptotic expansions in App. B lack explicit uniformity bounds. That's a minor concern given the simulation agreement.\n\nSoft spots: (i) the α=1 exception, which should be in the abstract; (ii) the higher-order increment argument relies on the approximation Δx(t;τ)/τ ≈ dx/dt without error control—this is heuristic but the result seems right and is verified numerically; (iii) no error bars or code are provided for the simulations, which is a reproducibility nit, not a correctness issue.\n\nBottom line: worth a serious referee. The fractional-order results are likely correct, and the paper fills a real gap in the RL-FBM literature. Send it to review with a request to fix the α=1 caveat and tighten the higher-order justification. I'd take it to reading group and would cite it for the TAMSD formulas, with the α=1 caveat noted.","headline":"Solid fractional-α results, but the 'any α>1/2' claim misses the ordinary Brownian exception at α=1.","tokens_in":30837,"tokens_out":3193,"would_cite":true,"duration_ms":24958,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["60G22","82C31"],"pacs":["05.40.Fb","05.40.-a"],"model":"deepseek-v4-flash","headline":"For the fractional Langevin equation far from equilibrium, the time-averaged mean-squared displacement converges to the mean-squared increment rather than the mean-squared displacement in the regime $1/2<\\alpha<3/2$—spurious nonergodicity…","keywords":["fractional Langevin equation","Riemann-Liouville fractional Brownian motion","FBM II","spurious nonergodicity","time-averaged mean squared displacement","mean squared increment","aging","anomalous diffusion"],"falsifier":"Simulate RL-FBM trajectories via the discretized integral in Appendix C for $\\alpha=0.8$ and $\\alpha=1.2$; fix a small lag $\\Delta$ and measure the ratio $\\langle\\delta^2(\\Delta)\\rangle/\\langle x^2(\\Delta)\\rangle$ as $T/\\Delta$ grows. The paper's claim requires this ratio to tend to $R(\\alpha)=\\frac{(2\\alpha-1)\\Gamma(\\alpha)^2}{\\Gamma(2\\alpha)|\\cos(\\pi\\alpha)|}$ (about 1.12 for $\\alpha=0.8$) rather than to 1. Equivalently, for $\\alpha>3/2$ the mean TAMSD must scale as $T^{2\\alpha-3}\\Delta^2$ and not converge to the MSD; a direct evaluation of the exact H-function expression (41) that shows convergence to the MSD for large $T$ would falsify the central claim.","tokens_in":29786,"feed_emoji":"📉","tokens_out":14330,"duration_ms":99809,"temperature":0.7,"pith_summary":"Fractional Langevin equations are used to model diffusion with memory. This paper studies the version driven by white noise rather than by the equilibrium noise required by the fluctuation-dissipation theorem—the 'fractional Langevin equation far from equilibrium'—whose solution is Riemann-Liouville fractional Brownian motion. It establishes that for fractional-derivative order $1/2<\\alpha<3/2$, the time-averaged mean-squared displacement (TAMSD) converges to the mean-squared increment (MSI, or structure function) rather than to the ensemble mean-squared displacement (MSD). The MSD and MSI differ by a constant prefactor, so a standard ergodicity test based on TAMSD-versus-MSD would incorrectly declare the process nonergodic even though its increments become stationary. For $\\alpha\\geq 3/2$ the first increments are genuinely nonergodic, but higher-order increments recover stationarity in successive $\\alpha$-bands, and strong aging restores ergodicity in $1/2<\\alpha<3/2$.","feed_headline":"Fractional Langevin motion fakes nonergodicity in time averages","feed_subtitle":"For 1/2 < α < 3/2 the time average tracks the mean-squared increment, not the mean-squared displacement.","key_machinery":"The object carrying the argument is Riemann-Liouville fractional Brownian motion (RL-FBM, also called FBM II), defined as the fractional integral of white noise, $x(t)=\\sqrt{2K_\\alpha}\\int_0^t \\frac{(t-t')^{\\alpha-1}}{\\Gamma(\\alpha)}\\xi(t')dt'$, which is the zero-initial-condition solution of the FLEFE. The derivational machinery is the asymptotic expansion of the Fox H-function, reached through the hypergeometric form $_2F_1$, in Appendices A and B. The load-bearing identity is the long-time stationary limit of the MSI integral, $\\int_0^\\infty [(1+s)^{\\alpha-1}-s^{\\alpha-1}]^2\\,ds+\\frac{1}{2\\alpha-1}=\\frac{\\Gamma(\\alpha)^2}{\\Gamma(2\\alpha)|\\cos(\\pi\\alpha)|}$ for $1/2<\\alpha<3/2$, which produces the spurious-nonergodicity prefactor; the power-logarithmic variant of the H-function expansion handles the borderline case $\\alpha=3/2$.","core_discovery":"The central discovery is that RL-FBM, the solution of the FLEFE, displays spurious nonergodicity. For all $\\alpha>1/2$ the mean TAMSD and the MSD do not coincide even in the limit of long trajectories; specifically, for $1/2<\\alpha<3/2$ the mean TAMSD converges to the stationary MSI, $\\langle\\delta^2(\\Delta)\\rangle\\sim \\langle x^2_\\Delta(t)\\rangle\\sim \\frac{2K_\\alpha}{\\Gamma(2\\alpha)|\\cos(\\pi\\alpha)|}\\Delta^{2\\alpha-1}$, which differs from the MSD $\\langle x^2(t)\\rangle = \\frac{2K_\\alpha}{(2\\alpha-1)\\Gamma(\\alpha)^2} t^{2\\alpha-1}$ by the factor $R(\\alpha)=\\frac{(2\\alpha-1)\\Gamma(\\alpha)^2}{\\Gamma(2\\alpha)|\\cos(\\pi\\alpha)|}$. Since the increments become asymptotically stationary in this regime, the TAMSD-to-MSI convergence is the physically appropriate ergodic criterion, and the TAMSD-to-MSD mismatch is therefore spurious. For $\\alpha\\geq 3/2$, the first-order increments are not ergodic, but the $(n+1)$th-order increments restore stationarity in the bands $(2n+1)/2<\\alpha<(2n+3)/2$; under strong aging ($t_a\\gg T$) the aged MSD and aged mean TAMSD coincide for $1/2<\\alpha<3/2$, restoring ergodicity.","pith_inferences":["Because the ratio R(α)=⟨δ²⟩/⟨x²⟩ = (2α−1)Γ(α)²/[Γ(2α)|cos(πα)|] is a closed-form function of α alone, an experimenter could use an observed TAMSD/MSD ratio to estimate α and to test whether a system is better described by RL-FBM than by FBM or FLE (which have R=1); the paper derives R(α) but does not propose this estimator.","The paper does not analyze confining potentials, but its closed-form MSI and TAMSD for the free process provide the ingredients to derive corresponding results for a harmonically trapped FLEFE, which would make contact with optical-tweezer experiments.","The strong-aging result implies that in systems where measurements start well after preparation (e.g., long climate or financial records), apparent ergodicity may reflect the aging time rather than equilibrium; distinguishing these requires varying the delay between preparation and measurement, an experimental protocol the paper does not discuss."],"forward_implications":["For single-particle-tracking data, a TAMSD that does not match the MSD is not by itself evidence of nonergodicity: for processes with asymptotically stationary increments, the structure function (MSI) is the correct reference, and RL-FBM is ergodic in that sense for 1/2 < α < 3/2.","In the regime α ≥ 3/2, ergodicity is recovered by moving to higher-order increments: the (n+1)th-order MSI becomes stationary in each band (2n+1)/2 < α < (2n+3)/2, so differentiating (or incrementing) trajectories more times restores standard ergodicity tests.","Strong aging (t_a ≫ T) makes the aged MSD and aged TAMSD coincide and the increment autocovariance function become lag-time-only dependent for 1/2 < α < 3/2; experiments that wait long enough after preparing a system will see apparent ergodicity even though the unaged process fails the usual test.","The FLEFE is a minimal Langevin-type model for systems that violate the fluctuation-dissipation theorem, remaining well-defined for all α > 1/2 and thus covering subdiffusive, superdiffusive, and ballistic-like regimes within one equation."],"supporting_citations":[{"why":"It introduced the FLEFE and supplied the exact solution in the form of RL-FBM that this paper analyzes.","marker":"[39]"},{"why":"It derived the MSI of RL-FBM and its stationary long-time asymptotics in the domain 1/2 < α < 3/2, which the paper extends to TAMSD and aging.","marker":"[41]"},{"why":"It defined fractional Brownian motion and its Riemann-Liouville form (FBM II), the process studied here.","marker":"[17]"},{"why":"It introduced the concept of spurious ergodicity breaking in Ornstein-Uhlenbeck processes, which this paper generalizes to RL-FBM.","marker":"[42]"},{"why":"It provided the H-function expansion formulas used to derive the long-time TAMSD and MSI asymptotics.","marker":"[66]"},{"why":"It supplied the hypergeometric identities (Pfaff transformation and Gauss summation) used in the derivations.","marker":"[64]"},{"why":"It defined the ergodicity criterion comparing TAMSD with MSD in single-particle tracking, the convention the paper argues should be replaced by the MSI comparison.","marker":"[15]"}],"fun_headline_variants":["Riemann-Liouville FBM shows spurious nonergodicity","Aging restores ergodicity in fractional Langevin motion","For 1/2<α<3/2, time averages track MSI, not MSD","Spurious nonergodicity cured by aging in fractional Langevin","Fractional motion fakes nonergodicity in time averages"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the Fox H-function asymptotic expansions in Appendix B (Eqs. B9 and B14) remain valid uniformly in $\\alpha$, including the special values $\\alpha=1$, $\\alpha=3/2$, and integer $\\alpha$, with the stated cancellations of leading terms; the paper gives no error bounds or uniformity proof, so if these expansions fail in any regime the claimed prefactor in Eq. (43), and hence the spurious-nonergodicity ratio, would be affected.","fun_headline_variants_meta":{"raw":{"variants":["Riemann-Liouville FBM shows spurious nonergodicity","Aging restores ergodicity in fractional Langevin motion","For 1/2<α<3/2, time averages track MSI, not MSD","Spurious nonergodicity cured by aging in fractional Langevin","Fractional motion fakes nonergodicity in time averages"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000482,"raw_usage":{"total_tokens":2495,"prompt_tokens":1171,"completion_tokens":1324,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":787,"completion_tokens_details":{"reasoning_tokens":1226}},"tokens_in":787,"tokens_out":1324,"duration_ms":11293,"temperature":1.0,"reasoning_tokens":1226,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T14:49:38.192322+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Simulate RL-FBM trajectories via the discretized integral in Appendix C for $\\alpha=0.8$ and $\\alpha=1.2$; fix a small lag $\\Delta$ and measure the ratio $\\langle\\delta^2(\\Delta)\\rangle/\\langle x^2(\\Delta)\\rangle$ as $T/\\Delta$ grows. The paper's claim requires this ratio to tend to $R(\\alpha)=\\frac{(2\\alpha-1)\\Gamma(\\alpha)^2}{\\Gamma(2\\alpha)|\\cos(\\pi\\alpha)|}$ (about 1.12 for $\\alpha=0.8$) rather than to 1. Equivalently, for $\\alpha>3/2$ the mean TAMSD must scale as $T^{2\\alpha-3}\\Delta^2$ and not converge to the MSD; a direct evaluation of the exact H-function expression (41) that shows convergence to the MSD for large $T$ would falsify the central claim.","supporting_citations":[{"cited_title":"Panja, Generalized Langevin equation formulation f or anomalous polymer dynamics, J","cited_arxiv_id":null,"evidence_quote":"It introduced the FLEFE and supplied the exact solution in the form of RL-FBM that this paper analyzes."},{"cited_title":"Lutz, Fractional Langevin equation, Phys","cited_arxiv_id":null,"evidence_quote":"It derived the MSI of RL-FBM and its stationary long-time asymptotics in the domain 1/2 < α < 3/2, which the paper extends to TAMSD and aging."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It introduced the concept of spurious ergodicity breaking in Ornstein-Uhlenbeck processes, which this paper generalizes to RL-FBM."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It provided the H-function expansion formulas used to derive the long-time TAMSD and MSI asymptotics."},{"cited_title":"Magdziarz, A","cited_arxiv_id":null,"evidence_quote":"It supplied the hypergeometric identities (Pfaff transformation and Gauss summation) used in the derivations."}],"review_version":1}