Stochastic Volterra Integral System

• Stochastic Volterra integral systems are defined by integrals with memory kernels, where the current state depends on the past trajectory.
• They incorporate various noises—such as Lévy, fractional Brownian, and Poisson—and are analyzed using SPDE-lifting, fixed-point, and operator-theoretic techniques.
• They have practical applications in optimal control, numerical approximation, and game theory, addressing challenges like singular kernels and path dependence.

Searching arXiv for recent and foundational papers on stochastic Volterra integral systems to ground the article. A stochastic Volterra integral system is a stochastic dynamical system in which the state at time $t$ depends on its past through Volterra-type integral terms, typically with coefficients or kernels depending on both the observation time and the integration time. In the literature summarized here, this class includes forward stochastic Volterra integral equations (SVIEs), backward stochastic Volterra integral equations (BSVIEs), forward-backward systems (FBSVIEs), delay and path-dependent variants, fuzzy-valued systems, and formulations under generalized noises such as Lévy noise, compensated Poisson random measures, fractional Brownian motion, cylindrical Wiener noise, and $G$-Brownian motion. The common structural feature is memory: unlike standard SDEs, the present state depends on the entire past trajectory through deterministic or stochastic kernels, so the resulting dynamics are typically non-Markovian (Benth et al., 2019).

1. Defining structure and principal model classes

A basic forward SVIE studied on a separable Hilbert space $U$ has the form

$$X(t)=x_0(t)+\int_0^t \mu(t,s,X(s))\,ds+\int_0^t \sigma(t,s,X(s-))\,dL(s), \qquad t\ge 0,$$

where $X$ is an adapted càdlàg $U$-valued process, $x_0$ is an $\mathcal F_0$-measurable initial process, $\mu$ is a drift kernel, $\sigma$ is a diffusion kernel, and $G$0 is a square-integrable Hilbert-space-valued Lévy process (Benth et al., 2019). In scalar or finite-dimensional Brownian settings, representative models take the form

$G$1

with weakly singular kernels, or

$G$2

when jumps are included (Wang et al., 2023, Khalaf et al., 2020).

The term “system” is used in several non-equivalent but related senses. It may denote a vector-valued forward equation; a coupled forward-backward configuration such as

$G$3

or a path-dependent state-control system with memory in both state and control channels (Wang et al., 2016, Li et al., 2023). It may also denote a delayed equation, where the coefficients depend on $G$4, or a symmetric implicit system in fuzzy stochastic analysis in which fuzzy stochastic integrals appear on both sides of the equation (Konan et al., 26 Mar 2026, Malinowski, 2024).

Backward formulations introduce additional asymmetry. A generic BSVIE on $G$5 is written as

$G$6

and the dependence on $G$7 leads to the adapted M-solution concept rather than a purely adapted solution (Wang et al., 2012). Anticipated BSVIEs extend this further by allowing future-point and future-average terms such as $G$8 and $G$9 in the generator (Yang et al., 24 Jan 2025).

This range of models suggests that the phrase “stochastic Volterra integral system” is best understood as a family of non-Markovian stochastic integral systems unified by Volterra memory kernels rather than as a single canonical equation.

2. Functional-analytic formulations and well-posedness theory

A central analytic development is the SPDE-lifting method for Hilbert-space-valued SVIEs. The state equation is embedded into a larger function space $U$0 of $U$1-valued functions on $U$2, equipped with the right-shift semigroup

$U$3

whose generator is denoted $U$4. The lifted coefficients are defined by

$U$5

and the associated first-order SPDE is

$U$6

Under Lipschitz and linear growth conditions on $U$7 and $U$8, the SPDE has a unique adapted càdlàg mild solution, and the original Volterra solution is recovered as the boundary trace $U$9. This yields existence, uniqueness, and finite second moments on finite horizons for the SVIE (Benth et al., 2019).

A concrete realization of this lifting is the weighted Filipović-type space $$X(t)=x_0(t)+\int_0^t \mu(t,s,X(s))\,ds+\int_0^t \sigma(t,s,X(s-))\,dL(s), \qquad t\ge 0,$$0, consisting of absolutely continuous $$X(t)=x_0(t)+\int_0^t \mu(t,s,X(s))\,ds+\int_0^t \sigma(t,s,X(s-))\,dL(s), \qquad t\ge 0,$$1-valued functions with norm

$$X(t)=x_0(t)+\int_0^t \mu(t,s,X(s))\,ds+\int_0^t \sigma(t,s,X(s-))\,dL(s), \qquad t\ge 0,$$2

where $$X(t)=x_0(t)+\int_0^t \mu(t,s,X(s))\,ds+\int_0^t \sigma(t,s,X(s-))\,dL(s), \qquad t\ge 0,$$3 is nondecreasing, $$X(t)=x_0(t)+\int_0^t \mu(t,s,X(s))\,ds+\int_0^t \sigma(t,s,X(s-))\,dL(s), \qquad t\ge 0,$$4, and $$X(t)=x_0(t)+\int_0^t \mu(t,s,X(s))\,ds+\int_0^t \sigma(t,s,X(s-))\,dL(s), \qquad t\ge 0,$$5. In this space the shift semigroup is strongly continuous and satisfies

$$X(t)=x_0(t)+\int_0^t \mu(t,s,X(s))\,ds+\int_0^t \sigma(t,s,X(s-))\,dL(s), \qquad t\ge 0,$$6

The same framework gives directly checkable assumptions on the kernels $$X(t)=x_0(t)+\int_0^t \mu(t,s,X(s))\,ds+\int_0^t \sigma(t,s,X(s-))\,dL(s), \qquad t\ge 0,$$7 and $$X(t)=x_0(t)+\int_0^t \mu(t,s,X(s))\,ds+\int_0^t \sigma(t,s,X(s-))\,dL(s), \qquad t\ge 0,$$8 (Benth et al., 2019).

Well-posedness in other settings is typically obtained by fixed-point or successive approximation arguments. For weakly singular Brownian SVIEs with $$X(t)=x_0(t)+\int_0^t \mu(t,s,X(s))\,ds+\int_0^t \sigma(t,s,X(s-))\,dL(s), \qquad t\ge 0,$$9, Lipschitz drift and diffusion imply existence and uniqueness, together with

$X$0

where $X$1 (Wang et al., 2023). For jump-driven SVIEs under a non-Lipschitz modulus condition governed by a concave $X$2, existence and uniqueness follow from Picard iteration and Bihari’s inequality rather than Grönwall’s inequality (Khalaf et al., 2020). For equations driven by fractional Brownian motion with Hurst parameter $X$3, the stochastic integral is interpreted pathwise as a Riemann–Stieltjes integral, and existence and uniqueness are proved in $X$4 via deterministic Volterra estimates and a weighted fixed-point argument (Besalú et al., 2010).

Path-dependent SVIEs push the dependence from $X$5 to the stopped trajectory $X$6. In that setting, strong well-posedness is established in an $X$7 Bochner-space framework, first for $X$8 and then locally in time for all $X$9, with Hölder regularity estimates and continuity of the flow in the initial condition (Gnabeyeu et al., 22 Mar 2026). This suggests that the well-posedness theory has evolved from kernel-specific finite-dimensional arguments to a broad operator-theoretic framework capable of handling path dependence, singularity, and nonstandard state spaces.

3. Backward, forward-backward, and comparison structures

The backward theory is structurally more delicate than the forward theory. For BSVIEs with generators depending on $U$0 but not $U$1, comparison can be established for adapted solutions under monotonicity in $U$2 and sign conditions on $U$3. By contrast, when the generator depends on $U$4, the adapted M-solution framework is required, and the comparison conclusion weakens to an inequality for conditional expectations of time-integrated solutions rather than a pointwise order for $U$5 (Wang et al., 2012).

The infinite-dimensional singular BSVIE theory substantially enlarges the admissible kernel class. The generator satisfies

$U$6

with $U$7-type spaces and $U$8 in a refined local singularity class $U$9. This framework covers fractional kernels, Volterra Heston kernels, and completely monotone kernels, and yields unique adapted M-solutions in natural Hilbert-valued spaces (Wang et al., 2023).

Forward-backward stochastic Volterra systems arise naturally in control. For optimal control problems of FBSVIEs with closed, possibly non-convex control regions, the state system couples a forward SVIE with a backward SVIE, and first-order necessary conditions are obtained through a duality principle between linear backward stochastic Volterra equations and linear stochastic Fredholm-Volterra equations with conditional expectation. The resulting maximum principle uses one adjoint system and only first-order differentiability of the coefficients, with the adjacent cone replacing spike variation (Wang et al., 2016).

A related but distinct forward-backward structure appears in linear-quadratic control of forward SVIEs. There the optimality system is a coupled FSVIE and Type-II BSVIE, and its decoupling requires a causal auxiliary state, a Type-III BSVIE, and a path-dependent Riccati equation on path space. The resulting open-loop optimal control admits a causal state feedback representation, and when the control appears only in the diffusion term, the feedback reduces to a Markovian state feedback despite the non-Markovian state equation (Wang et al., 2022).

These developments show that stochastic Volterra integral systems do not merely extend SDEs by memory terms; they require distinct notions of solution, comparison, duality, and decoupling.

4. Optimal control, maximum principles, and games

The control theory of stochastic Volterra systems is dominated by Pontryagin-type principles rather than dynamic programming, because the Volterra memory destroys the usual semigroup flow property. In the linear-quadratic setting for forward SVIEs, the state equation

$x_0$0

is paired with a quadratic cost, and the open-loop optimal control is represented as

$x_0$1

where $x_0$2 is a causal auxiliary state and $x_0$3 is derived from a path-dependent Riccati equation. Under $x_0$4, $x_0$5, and $x_0$6, the Riccati equation has a unique strongly regular solution (Wang et al., 2022).

For generalized Volterra control systems formulated as stochastic integral-differential equations with memory in both state and control,

$x_0$7

the adjoint equation is a generalized anticipated BSDE whose drift contains future conditional expectations of integrals against the kernels $x_0$8 and $x_0$9. The stochastic maximum principle then yields the first-order condition

$\mathcal F_0$0

with a corresponding explicit optimal control in the linear-quadratic example (Li et al., 2023).

Delay changes the structure again. For delayed SVIEs,

$\mathcal F_0$1

Hida–Malliavin calculus is used to derive an adjoint anticipated backward stochastic Volterra integral equation (ABSVIE). The representation

$\mathcal F_0$2

is central in both necessary and sufficient maximum principles. The sufficient principle is expressed through conditional maximization of a Volterra Hamiltonian, and the necessary principle is the conditional stationarity of $\mathcal F_0$3 (Konan et al., 26 Mar 2026).

Game-theoretic extensions use ABSVIEs as adjoint equations for nonzero-sum stochastic differential games of delayed Volterra systems. In that setting, the generator may contain both pointwise time-advanced and average time-advanced terms, and well-posedness, comparison, and Malliavin regularity are established before deriving the Nash equilibrium conditions. In the linear-quadratic SDVIE game, the equilibrium controls take the explicit form

$\mathcal F_0$4

(Yang et al., 24 Jan 2025).

A consistent pattern across these results is that memory terms migrate into the adjoint equation rather than disappearing under augmentation. This suggests that, for stochastic Volterra integral systems, the adjoint object is often itself a Volterra-type backward equation.

5. Numerical approximation and computational schemes

Numerical analysis of stochastic Volterra systems is shaped by kernel singularity and non-Markovian dependence. For weakly singular SVIEs with non-differentiable drift, the randomized Milstein scheme avoids Taylor expansion of the drift by introducing i.i.d. auxiliary variables $\mathcal F_0$5 into the drift quadrature. The principal strong error estimate is

$\mathcal F_0$6

and the numerical experiments recover the predicted slopes $\mathcal F_0$7 for $\mathcal F_0$8 and $\mathcal F_0$9 for $\mu$0 (Wang et al., 2023). The same paper addresses the practical simulation of multiple singular stochastic integrals via a Riemann–Stieltjes discretization of nested integrals.

For path-dependent SVIEs, an interpolated $\mu$1-integrated Euler–Maruyama scheme is designed to respect the Volterra structure by integrating the kernels over each mesh cell rather than using point evaluations. Under Hölder-Lipschitz assumptions on the coefficients and regularity assumptions on the kernels, the fixed-time strong error satisfies

$\mu$2

and the uniform-in-time error is controlled by

$\mu$3

(Gnabeyeu et al., 22 Mar 2026).

Other numerical frameworks exploit the specific Volterra structure differently. The Walsh-function method approximates the linear SVIE

$\mu$4

by expanding the solution and kernels in a Walsh basis and converting the equation into a finite algebraic system

$\mu$5

Under Lipschitz assumptions on $\mu$6, $\mu$7, and $\mu$8, the basis approximation error is $\mu$9 and the full solution error satisfies

$\sigma$0

(Paikaray et al., 2023).

A more global high-order method is the cubature method for Stratonovich SVIEs. Because the solution is generally not a semimartingale, the method first constructs a stochastic Taylor expansion via a functional Itô formula on an auxiliary two-time process $\sigma$1, then matches the relevant Volterra signatures by a discrete cubature measure. In the multi-block case the global approximation error is of order

$\sigma$2

up to explicit coefficient-dependent factors (Feng et al., 2021).

These schemes illustrate two recurring numerical themes: memory must be discretized rather than ignored, and kernel singularity usually determines the convergence rate.

6. Generalized noises, values, and integration frameworks

The concept of a stochastic Volterra integral system extends beyond standard Itô equations. In white-noise analysis, Brownian-driven Volterra processes with volatility modulation are treated on the Potthoff–Timpel distribution space $\sigma$3. For

$\sigma$4

integration with respect to $\sigma$5 is defined by

$\sigma$6

where

$\sigma$7

A generalized volatility modulation via Wick product is also introduced, and the Wick formulation coincides with the pointwise product under strong independence (Barndorff-Nielsen et al., 2013).

An infinite-dimensional analogue studies Hilbert-valued volatility-modulated Volterra processes driven by a cylindrical Wiener process. The stochastic integral

$\sigma$8

is defined as a Skorohod integral plus a correction term involving the Malliavin derivative of $\sigma$9, and the resulting calculus includes an Itô formula and a random-field formulation relevant to SPDE kernels (Benth et al., 2013).

Other generalizations change the state space rather than the noise. Symmetric fuzzy stochastic Volterra integral equations with constant retardation consider fuzzy-valued processes $G$00, with fuzzy stochastic Lebesgue–Aumann integrals and Hukuhara differences on both sides of the equation. Under hypotheses $G$01–$G$02, the initial value problem has a unique solution and depends continuously on the initial fuzzy path, kernels, and nonlinearities (Malinowski, 2024).

A different extension replaces classical probability by sublinear expectation. For $G$03-SVIEs of the form

$G$04

existence, uniqueness, mean-square continuity, and pathwise continuity of a modification are established in $G$05-type spaces. A comparison theorem is proved for the separable diffusion form $G$06, and the method avoids classical assumptions on partial derivatives of the coefficients by combining quasilinearization with a two-step approximation (Zhao et al., 30 Apr 2025).

These examples show that the Volterra structure is compatible with generalized state spaces, generalized stochastic integration, and uncertainty in the driving law.

7. Long-time behavior, stationarity, and asymptotics

Long-time analysis for stochastic Volterra systems differs sharply from the Markovian case. In the SPDE-lifting framework for homogeneous Hilbert-space-valued SVIEs, limiting distributions are obtained by importing invariant-measure criteria from the lifted SPDE. In the weighted space $G$07, if the coefficients vanish at infinity and the dissipativity condition

$G$08

holds on the vanishing-at-infinity subspace $G$09, then $G$10 converges weakly to a limit $G$11 on $G$12 (Benth et al., 2019). A second theorem handles the persistent constant component through a one-sided monotonicity condition on the drift at infinity (Benth et al., 2019).

A more recent stationarity theory emphasizes a negative conclusion. For forward SVIEs with convolution kernel $G$13, strong stationarity is essentially impossible unless the kernel is constant or the system is degenerate. In the affine mean-reverting case

$G$14

the paper shows that one can nevertheless induce a “fake stationary regime” by choosing a deterministic stabilizer $G$15 in $G$16. The stabilizer solves

$G$17

and forces constant first two moments and constant expected squared diffusion in type-I fake stationarity (Gnabeyeu et al., 5 Nov 2025).

The same work studies $G$18-confluence across initial values and long-run functional weak limits of the shifted process $G$19. Under suitable assumptions, the family is $G$20-tight in $G$21, and any limit process is weakly $G$22-stationary when the original process is in fake stationary regime of type I. The framework applies to fractional kernels

$G$23

and exponential-fractional kernels

$G$24

covering both rough behavior for $G$25 and long-memory or persistence for $G$26 (Gnabeyeu et al., 5 Nov 2025).

This suggests a useful correction to a common intuition: memory does not automatically produce stationary long-run dynamics. In the Volterra setting, stationarity usually has to be reformulated, weakened, or induced by additional deterministic structure rather than inherited from the kernel alone.