alphaFeedback

WorksMathematics

Markovian Lifts and Multifactor Approximations for Weighted Sub-fractional Brownian Motion

AI contributions
Not recorded for this paper
View details

Abstract

We construct a Hilbert-space Markovian lift of weighted sub-fractional Brownian motion for 0<b<20<b<2, including the logarithmic case b=1b=1. A Laplace factorization represents the process as the boundary value of a forward curve. In the weighted Sobolev space with derivative weight rγ/(1+r)2r^γ/(1+r)^2, 0<γ<10<γ<1, right translation is strongly continuous and the noise operator is Hilbert–Schmidt if and only if b+γ>1b+γ>1. The state norm bounds the total variation of initial curves on compact intervals. For locally bounded deterministic temporal weights, this bound allows the driven scalar equation to be restarted from every lifted state. Under spatial Lipschitz assumptions, with additional time and derivative regularity in the Young case, the solution and forward curve form a time-inhomogeneous Markov process with CbC_b-Feller transitions; both coefficients may have linear growth. Discretization of the Laplace variable gives independent Ornstein–Uhlenbeck factors coupled to the original process through the same Gaussian noise. The nonlinear approximations converge uniformly in time in every finite LpL^p. For uniform spectral partitions with balanced cutoffs, the solution rate is M−(b−1)/3M^{-(b-1)/3} for 1<b<21<b<2, while the uniform driver rate is M−b/4log⁡(2+M)M^{-b/4}\sqrt{\log(2+M)} for 0<b≤10<b\le1. Under stronger coefficient assumptions, we derive error bounds for Young solutions and for the fully discrete scheme, retaining the dependence on the spectral cutoff and specifying joint mesh conditions.

keywords: MSC 2020: 60G15; 60H10; 60H15; 60H35; 65C30

Introduction and main results

Weighted sub-fractional Brownian motion is a Gaussian process whose covariance combines a temporal weight with a memory kernel. Its current value generally does not determine the conditional law of future increments. We record the contribution of past Brownian noise to future values in a forward curve, which provides a Markovian state for the driver. Discretizing the same covariance representation gives finitely many Ornstein–Uhlenbeck factors. We study the resulting approximations of both the driver and the scalar equation it drives.

Fix 0<b<20 < b < 2 and a non-negative Borel function f∈Lloc1(R+)f \in L_{\mathrm{loc}}^{1}(\mathbb{R}_{+}). We write ζf,b\zeta_{f,b} for the centered Gaussian process with covariance

Rf,b(s,t)=∫0s∧tf(u)Qb(s−u,t−u) du,s,t≥0,(1)R_{f,b}(s,t) = \int_{0}^{s \wedge t} f(u)Q_b(s-u,t-u)\,\mathrm{d}u,\qquad s,t \ge0, \tag*{(1)}

where

Qb(x,y)={xb+yb−(x+y)b1−b,b≠1,(x+y)log⁡(x+y)−xlog⁡x−ylog⁡y,b=1,x,y≥0,(2)Q_b(x,y) = \begin{cases} \dfrac{x^b+y^b-(x+y)^b}{1-b}, & b \ne1,\\ (x+y)\log(x+y)-x\log x-y\log y, & b=1, \end{cases} \qquad x,y \ge0, \tag*{(2)}

with 0log⁡0:=00\log0 := 0. For b≠1b \ne1, this is the weighted sub-fractional Brownian covariance in [18 Section 2.1, Eq. (4)]. The power-weight family originates in occupation-time fluctuation limits of age-dependent branching systems [16 Theorem 2.1(i), Eq. (2.1)]. The logarithmic case b=1b = 1 appears for power weights in [16 Theorem 2.1(ii), p. 596]. Besides its origin in occupation-time fluctuation limits, this Gaussian family has been used in animal movement modeling. Ramírez-González, Murillo-Salas, and Sun [20] fitted its logarithmic member to bat trajectories using longitude and latitude observations. The same family is included among the Gaussian models considered in the multidimensional telemetry analysis of [19]. Inference for stochastic differential equations driven by weighted sub-fractional Brownian motion is studied in [17]. The Euler approximation is used to construct a trajectory likelihood that retains the dependence between the Gaussian driving increments. Neural-network and radial-basis representations are used to estimate the drift and diffusion coefficients together with the normalized temporal weight and covariance exponent.

The regularity of the boundary process changes at the critical value b=1b = 1. If f∈Lloc∞(R+)f \in L_{\mathrm{loc}}^{\infty}(\mathbb{R}_{+}) and 0<b<10 < b < 1, then ζf,b\zeta_{f,b} admits a modification whose paths are Hölder continuous of every order

0<η<1+b2.0 < \eta< \frac{1+b}{2}.

At b=1b = 1 the increment variance acquires a logarithmic correction, but there is still a modification with Hölder continuous paths of every order below one. For later statements it is convenient to write

Hbδ:={(1+b)/2,0<b<1,1,b=1,H_b^\delta:= \begin{cases} (1+b)/2, & 0<b<1,\\ 1, & b=1, \end{cases}

where the superscript indicates that the endpoint is not asserted. Hence the nonlinear equation is a Young equation for 0<b≤10<b\leq1. When 1<b<21<b<2, the boundary admits an absolutely continuous realization and the same equation is interpreted in the Lebesgue–Stieltjes sense.

Relation to the literature.

Laplace representations of fractional Gaussian processes lead to infinite-dimensional Markovian states and finite-factor approximations; see [7, 8]. Affine representations of fractional processes are developed in [14], and strong approximation rates in [13]. For stochastic Volterra equations, boundary representations through first-order SPDEs are studied in [5]; weighted Sobolev realizations and polynomially weighted shift spaces are treated in [6, 15]. Related results on affine Volterra processes and their Markovian lifts appear in [1, 10, 12].

The nonlinear equations in [6, 15] are formulated as stochastic Volterra equations with Brownian Itô integrals. In our construction, the Gaussian forward curve is also an Itô stochastic evolution, while its boundary drives the scalar equation in the Young or Lebesgue–Stieltjes sense. The main issue for the nonlinear equation is to restart it from an arbitrary curve in the state space: such a curve need not have the Hölder regularity of the Gaussian process. The total-variation bound below provides the required control. For the finite-factor approximations, we use the same Gaussian noise and estimate the driver error in the path norms used by the solution map. This gives convergence under the well-posedness assumptions; quantitative bounds require the additional hypotheses specified below.

The state space and the Hilbert–Schmidt condition.

For

μb(dλ):=1Γ(2−b)λ1−b dλ,λ>0,(3)\mu_b(\mathrm{d}\lambda) := \frac{1}{\Gamma(2-b)}\lambda^{1-b}\,\mathrm{d}\lambda,\qquad\lambda>0, \tag*{(3)}

and

ϕr(λ):=b1−e−λrλ,r≥0,(4)\phi_r(\lambda) := \sqrt{b}\frac{1-e^{-\lambda r}}{\lambda},\qquad r\geq0, \tag*{(4)}

let Ub:=L2((0,∞),μb)U_b := L^2((0,\infty),\mu_b). Here λ\lambda is the positive relaxation rate in the Laplace representation; the term spectral below refers to this variable. Proposition 2.1 proves Qb(x,y)=⟨ϕx,ϕy⟩UbQ_b(x,y)=\langle\phi_x,\phi_y\rangle_{U_b}.

For 0<γ<10<\gamma<1, we use the polynomial weight

wγ(r):=rγ(1+r)2,r>0,(5)w_\gamma(r) := \frac{r^\gamma}{(1+r)^2},\qquad r>0, \tag*{(5)}

and the Hilbert space

Eγ:={θ∈ACloc(R+):∥θ∥γ2:=∣θ(0)∣2+∫0∞rγ(1+r)2∣θ′(r)∣2 dr<∞}.(6)\mathcal{E}_{\gamma}:=\left\{\theta\in AC_{\mathrm{loc}}(\mathbb{R}_{+}):\|\theta\|_{\gamma}^{2}:=|\theta(0)|^{2}+\int_{0}^{\infty}\frac{r^{\gamma}}{(1+r)^{2}}|\theta'(r)|^{2}\,\mathrm{d}r<\infty\right\}. \tag*{(6)}

The value assigned to wγw_{\gamma} at the single point r=0r=0 is immaterial. rγr^{\gamma} controls the singularity of the feature derivative at the origin, while (1+r)−2(1+r)^{-2} provides sufficient decay at infinity and is stable under right translation. Consider the feature rule

(Gbh)(r):=⟨ϕr,h⟩Ub,h∈Ub, r≥0.(7)(G_bh)(r):=\langle\phi_r,h\rangle_{\mathcal{U}_b},\qquad h\in\mathcal{U}_b,\ r\ge0. \tag*{(7)}

This rule defines a Hilbert–Schmidt operator from Ub\mathcal{U}_b to Eγ\mathcal{E}_{\gamma} precisely when

b+γ>1.(8)b+\gamma>1. \tag*{(8)}
∥Gb∥HS(Ub,Eγ)2=b2b−2Γ(b+γ−1)Γ(3−b−γ).\|G_b\|_{\mathrm{HS}(\mathcal{U}_b,\mathcal{E}_{\gamma})}^{2}=b2^{b-2}\Gamma(b+\gamma-1)\Gamma(3-b-\gamma).

The parameters of the lifted equation satisfy

0<γ<1,b+γ>1.(9)0<\gamma<1,\qquad b+\gamma>1. \tag*{(9)}

Thus 1−b<γ<11-b<\gamma<1 for 0<b<10<b<1 and 0<γ<10<\gamma<1 for 1≤b<21\le b<2. Lemma 2.4 also covers γ=0\gamma=0 when 1<b<21<b<2; we use positive γ\gamma to state the construction in one parameter range.

Write (Stθ)(r)=θ(t+r)(S_t\theta)(r)=\theta(t+r) for right translation and let AA be its generator. Lemma 2.3 proves that (St)t≥0(S_t)_{t\ge0} is strongly continuous on Eγ\mathcal{E}_{\gamma} and that boundary evaluation δ0θ=θ(0)\delta_0\theta=\theta(0) is continuous. Theorem 2.7 then constructs the mild solution of

dΘt=AΘt dt+f(t)Gb dWt,Θ0=0,\mathrm{d}\Theta_t=A\Theta_t\,\mathrm{d}t+\sqrt{f(t)}G_b\,\mathrm{d}W_t,\qquad\Theta_0=0,

where WW is cylindrical Brownian motion on Ub\mathcal{U}_b. Its boundary is Θt(0)=ζt,f,b\Theta_t(0)=\zeta_{t,f,b}, while

Θt(r)=E[ζt+r,f,b∣Ft].\Theta_t(r)=\mathbb{E}[\zeta_{t+r,f,b}\mid\mathcal{F}_t].

Here Ft\mathcal{F}_t is the filtration carrying the driving Brownian motions. The curve records the contribution of past noise to future values; the remaining contribution uses independent future Brownian increments. This decomposition gives the Markov property.

The nonlinear equation.

Fix T>0T>0 and consider

Xt=x0+∫0ta(Xs,s) ds+∫0tσ(Xs,s) dζs,f,b,0≤t≤T,(10)X_t=x_0+\int_0^t a(X_s,s)\,\mathrm{d}s+\int_0^t \sigma(X_s,s)\,\mathrm{d}\zeta_{s,f,b},\qquad0\le t\le T, \tag*{(10)}

where x0∈Rx_0\in\mathbb{R} is deterministic and a,σ:R×R+→Ra,\sigma:\mathbb{R}\times\mathbb{R}_{+}\to\mathbb{R} are the drift and diffusion coefficients.

Assume additionally that f∈Lloc∞(R+)f \in L^\infty_{\mathrm{loc}}(\mathbb{R}_+). Under Assumption 3.5, Theorem 3.10 gives the Markov property for (X,Θ)(X,\Theta) and CbC_b-Feller continuity of its time-dependent transitions. In the Young regime this assumption allows linear growth of both coefficients. The drift and diffusion are continuous and globally Lipschitz in space; the diffusion is Lipschitz in time and has a bounded, uniformly α\alpha-Hölder spatial derivative, with α>(1−b)/(1+b)\alpha> (1-b)/(1+b) for b<1b < 1 and α>0\alpha> 0 for b=1b = 1. For b>1b > 1, only continuity and uniform global spatial Lipschitz continuity of both coefficients are required. Table 1 summarizes the parameter ranges for this Hilbert–Schmidt realization and the corresponding pathwise integrals.

Range of bbState exponent γ\gammaPathwise integral
0<b<10 < b < 11−b<γ<11-b < \gamma< 1Young
b=1b = 10<γ<10 < \gamma< 1Young
1<b<21 < b < 20<γ<10 < \gamma< 1Lebesgue–Stieltjes

Table 1. State exponents and pathwise interpretations used in the lifted theory.

Approximation and coefficient regularity.

References

  1. [1]E. Abi Jaber, M. Larsson, and S. Pulido. Affine Volterra processes. Annals of Applied Probability, 29(5):3155–3200, 2019. https://doi.org/10.1214/19-AAP1477.
  2. [2]R. J. Adler and J. E. Taylor. Random Fields and Geometry. Springer, 2007. https://doi.org/10.1007/978-0-387-48116-6.
  3. [3]A. Alfonsi and A. Kebajer. Approximation of stochastic Volterra equations with kernels of completely monotone type. Mathematics of Computation, 93(346):643–677, 2024. https://doi.org/10.1090/mcom/3911.
  4. [4]C. Bayer and S. Breneis. Markovian approximations of stochastic Volterra equations with the fractional kernel. Quantitative Finance, 23(1):53–70, 2023. https://doi.org/10.1080/14697688.2022.2139193.
  5. [5]F. E. Benth, N. Detering, and P. Krühner. Stochastic Volterra integral equations and a class of first-order stochastic partial differential equations. Stochastics, 94(7):1054–1076, 2022. https://doi.org/10.1080/17442508.2021.2019738.
  6. [6]L. A. Bianchi, S. Bonaccorsi, O. Cañadas, and M. Friesen. Limit theorems for stochastic Volterra processes. Stochastics and Partial Differential Equations: Analysis and Computations, 14:2450–2521, 2026. https://doi.org/10.1007/s40072-026-00428-w.
  7. [7]P. Carmona and L. Coutin. Fractional Brownian motion and the Markov property. Electronic Communications in Probability, 3:95–107, 1998. https://doi.org/10.1214/ECP.v3-998.
  8. [8]P. Carmona, L. Coutin, and G. Montseny. Approximation of some Gaussian processes. Statistical Inference for Stochastic Processes, 3:161–171, 2000. https://doi.org/10.1023/A:1009999518898.
  9. [9]N. D. Cong, L. H. Duc, and P. T. Hong. Nonautonomous Young differential equations revisited. Journal of Dynamics and Differential Equations, 30:1921–1943, 2018. https://doi.org/10.1007/s10884-017-9634-y.
  10. [10]C. Cuchiero and J. Teichmann. Generalized Feller processes and Markovian lifts of stochastic Volterra processes: the affine case. Journal of Evolution Equations, 20:1301–1348, 2020. https://doi.org/10.1007/s00028-020-00557-2.
  11. [11]G. Da Prato and J. Zabczyk. Stochastic Equations in Infinite Dimensions. Cambridge University Press, second edition, 2014. https://doi.org/10.1017/CBO9781107295513.
  12. [12]Y. Hamaguchi. Markovian lifting and asymptotic log-Harnack inequality for stochastic Volterra integral equations. Stochastic Processes and their Applications, 178:104482, 2024. https://doi.org/10.1016/j.spa.2024.104482.
  13. [13]P. Harms. Strong convergence rates for Markovian representations of fractional processes. Discrete and Continuous Dynamical Systems–B, 26(10):5567–5579, 2021. https://doi.org/10.3934/dcdsb.2020367.
  14. [14]P. Harms and D. Stefanovits. Affine representations of fractional processes with applications in mathematical finance. Stochastic Processes and their Applications, 129(4):1185–1228, 2019. https://doi.org/10.1016/j.spa.2018.04.010.
  15. [15]F. Huber. Markovian lifts of stochastic Volterra equations in Sobolev spaces: solution theory, an Itô formula and invariant measures. Preprint, arXiv:2406.10352v4 [math.PR], 11 July 2026.arxiv.org/abs/2406.10352
  16. [16]J. A. López-Mimbela, A. Murillo-Salas, and J. H. Ramírez-González. Occupation time fluctuations of an age-dependent critical binary branching particle system. ALEA, Latin American Journal of Probability and Mathematical Statistics, 21:593–625, 2024. https://doi.org/10.30757/ALEA.v21-24.
  17. [17]J. H. Ramírez-González. Inference for stochastic differential equations driven by weighted sub-fractional Brownian motion using neural networks and the Euler approximation. Preprint, arXiv:2610.00793 [stat.ML], 2026.
  18. [18]J. H. Ramírez-González. Weighted Sub-fractional Brownian Motion: Covariance Structure, Path Properties and Euler Approximation. Preprint, arXiv:2409.04798v2, 28 August 2026.
  19. [19]J. H. Ramírez-González, E. A. Chacón-Montalván, P. Moraga, and Y. Sun. Multidimensional integral fractional Ornstein–Uhlenbeck process with an application to animal movement. Environmetrics, 37(6):e70143, 2026. https://doi.org/10.1002/env.70143.
  20. [20]J. H. Ramírez-González, A. Murillo-Salas, and Y. Sun. A new family of Gaussian processes for modeling animal movement: Application to bat telemetry data. Environmetrics, 37(7):e70155, 2026. https://doi.org/10.1002/env.70155.

Paper details

vlnc:W_d68ct5yg
v1 · latestHistory (1)

published · Fresh — added to the Hub within the last 30 days

added ·

Usage

views · 0

downloads · 0

A view counts once per browser per day, a download each time the PDF is requested.

CC BY 4.0

Original source

  • Mathematics

Not recorded for this paper · View statement →

Looking for related papers…
Contents