Skip to main content

Path-dependent backward stochastic Volterra integral equations with jumps, differentiability and duality principle

Abstract

We study the existence and uniqueness of a solution to path-dependent backward stochastic Volterra integral equations (BSVIEs) with jumps, where path-dependence means the dependence of the free term and generator of a path of a càdlàg process. Furthermore, we prove path-differentiability of such a solution and establish the duality principle between a linear path-dependent forward stochastic Volterra integral equation (FSVIE) with jumps and a linear path-dependent BSVIE with jumps. As a result of the duality principle we get a comparison theorem and derive a class of dynamic coherent risk measures based on path-dependent BSVIEs with jumps.

1 Introduction

Since Pardoux and Peng provided in Pardoux and Peng (1990) the existence and uniqueness to a non-linear backward stochastic differential equation (BSDE) with Lipschitz continuous generator and measurable terminal condition, this type of equation has attracted enormous attention in probability theory and its applications. Their research has been extended in terms of BSDEs with Poisson jumps by Tang and Li (1994) and Barles et al. (1997) and for BSDEs driven by random measures by Becherer (2006). Moreover, path-dependent BSDEs with jumps, where the generator and the terminal condition both depend on a path of a càdlàg process, were studied by Wang (2015) and Kromer et al. (2017). BSDEs are important because of their connections to financial mathematics as studied in El Karoui et al. (1997), Rong (1997), Carmona (2008), Delong (2013), Crépey (2013), and many others.

Furthermore, BSDEs are special cases of backward stochastic Volterra integral equations (BSVIEs). Stochastic Volterra integral equations were first investigated by Berger and Mizel (1980) and extended by Protter (1985) and Pardoux and Protter (1990). Lin (2002) studied these equations in an L2-framework under global Lipschitz conditions on the coefficients the existence and uniqueness of a solution to backward stochastic nonlinear Volterra integral equations of the form

$$Y(t)=\xi+\int_{t}^{T}f\left(t,s,Y(s),Z(t,s)\right)ds+\int_{t}^{T}[g(t,s,Y(s))+Z(t,s)]dW(s),\:t\in[0,T]. $$

This result was extended by Aman and N’Zi (2005) by weakening the global Lipschitz condition to a local one. Wang and Zhang (2007) proved existence and uniqueness of a solution to BSVIEs with jumps and non-Lipschitz coefficients. Moreover, they studied the regularity of this solution. A preliminary theory for BSVIEs was established by Yong (2006) considering a BSVIE of the form

$$Y(t)=f(t)-\int_{t}^{T}h\left(t,s,Y(s),Z(t,s),Z(s,t)\right)ds-\int_{t}^{T}Z(t,s)dW(s),\quad t\in[0,T]. $$

These results were extended by Anh and Yong (2006) and Yong (2008), where they considered the notion of adapted M-solutions and established the well-posedness of adapted M-solutions to BSVIEs. Furthermore, Yong (2006) and Wang and Yong (2015) introduced the duality principle between linear forward stochastic Volterra integral equations (FSVIEs) and linear BSVIEs to prove the comparison theorem. Their study was motivated firstly by the use of BSVIEs in the context of dynamic risk measures, as proposed by Yong (2007), and secondly by the connection between BSVIEs and stochastic differential utilities as introduced and studied in Duffie and Epstein (1992) and Lazrak and Quenez (2003).

Further, Ren (2010) studied the existence and uniqueness of a solution to BSVIEs driven by a cylindrical Brownian motion on a separable Hilbert space and a Poisson random measure with a non-Lipschitz coefficient. Moreover, he proved a duality principle between linear FSVIEs with jumps and linear BSVIEs with jumps.

The aim of this work is to study path-dependent BSVIEs with jumps, where the path-dependence means the dependence of the free term and the generator on a path of a càdlàg process. These types of path-dependent stochastic equations came into the focus of stochastic analysis, particularly since the development of a path-dependent Itô formula as in Dupire (2009), Cont and Fournié (2010), and Levental et al. (2013). It was first considered for backward equations by Peng and Wang in Peng and Wang (2016) and developed in Wang (2015). They also showed differentiability of path-dependent BSDEs in the Dupire sense and applied it to path-dependent partial differential equations. We extend their differentiability result to the BSVIE case. The Differentiability of BSDEs and BSVIEs was considered in a different context, e.g., Ankirchner et al. (2007) and Kromer and Overbeck (2014); Kromer and Overbeck (2017); Kromer et al. (2017). One application of BSVIEs is the risk measures associated with them. Path-dependency can then be incorporated to exhibit the dependency of the risk measures and the underlying portfolio from a driving exogenous factor process or asset price process X. Differentiability can then be used to measure the sensitivity of the risk with respect to this driving factor process. In Kromer and Overbeck (2014); Kromer and Overbeck (2017) a different kind of differentiability was used for the purpose of capital allocation for risk measures.

We derive an existence and uniqueness result for adapted M-solutions to path-dependent BSVIEs with jumps. Furthermore, we prove path-differentiability of this solution, where we use the functional Itô formula introduced by Dupire (2009) and extended by Cont and Fournié (2010) and Levental et al. (2013). Moreover, we establish a duality principle between linear path-dependent FSVIEs with jumps and linear path-dependent BSVIEs with jumps. The main challenge of our work is to handle the dependence of the generator and terminal condition of the path of a càdlàg process in the differentiability result, duality principle, and comparison theorem for path-dependent BSVIEs with jumps.

This work is organized as follows. In Section 2, we present necessary notations, spaces, and results to path-dependent BSDEs with jumps. Furthermore, we introduce the considered path-dependent BSVIE with jumps. In Section 3, we provide the existence and uniqueness result for path-dependent BSVIEs with jumps and give a stability estimate to such a solution, extending the results of Peng and Wang (2016) on path-dependent BSDE without jumps and Yong (2008). In Section 4, we study the path-differentiability of solutions to path-dependent BSVIEs with jumps, as done for BSDEs in Kromer et al. (2017). In Section 5, we will establish a duality principle between a linear path-dependent FSVIE with jumps and a linear path-dependent BSVIE with jumps as well as a comparison theorem for adapted M-solutions to BSVIEs, again generalizing the results in Wang and Yong (2015) by introducing path-dependency and jumps. Finally, a class of dynamic coherent risk measures based on path-dependent BSVIEs with jumps are derived in Section 6.

2 Preliminaries

Let t≥0 and T>t be fixed. Throughout this paper, we consider a probability space \((\Omega,\mathbb {F},\mathbb {P})\) with the filtration \(\mathbb {F}=\left (\mathcal {F}^{t}_{s}\right)_{s\in [t,T]}\) defined by , where (W(s))s[0,T] is a d-dimensional Brownian motion, N a random measure, and \(\mathcal {N}\) the \(\mathbb {P}\)-null set.

The space of all càdlàg \(\mathbb {R}^{d}\)-valued functions on [0,T] endowed with the supremum metric is denoted by \(D\left ([0,T],\mathbb {R}^{d}\right)\). For each \(\gamma \in D\left ([0,T],\mathbb {R}^{d}\right)\) the value of γ at s[0,T] is γ(s) and we define the path of γ stopped at t[0,T], denoted by γ t , as for all u[0,T]. The path shifted by \(h\in \mathbb {R}^{d}\) in t is denoted by , for u[0,T]. Furthermore, for \(t,\bar {t}\in [0,T]\) it is

$$ \Vert\gamma_{t}-\bar{\gamma}_{\bar{t}}\Vert_{D}:= \Vert \gamma_{t}-\bar{\gamma}_{\bar{t}}\Vert_{D(t,\bar{t})} =\Vert \gamma_{t}-\bar{\gamma}_{\bar{t}}\Vert+|t-\bar{t}| =\sup_{s\in[0,T]}|\gamma_{t}(s)-\bar{\gamma}_{\bar{t}}(s)| +|t-\bar{t}|. $$

For any \(\gamma \in D\left ([0,T],\mathbb {R}^{d}\right)\) and a fixed t[0,T], we are looking for a unique solution \((Y_{\gamma _{t}}(\cdot),Z_{\gamma _{t}}(\cdot,\cdot), U_{\gamma _{t}}(\cdot,\cdot,\cdot))\) to the following path-dependent backward stochastic Volterra integral equation (BSVIE) with jumps

$$ \begin{aligned} Y_{\gamma_{t}}(s)&=\Phi(X_{\gamma_{t},s}) +\int_{s}^{T}f(s,r,X_{\gamma_{t},r-},Y_{\gamma_{t}}(r-), Z_{\gamma_{t}}(s,r),Z_{\gamma_{t}}(r,s),U_{\gamma_{t}}(s,r), U_{\gamma_{t}}(r,s))dr \\ &\quad-\int_{s}^{T}Z_{\gamma_{t}}(s,r)dW(r) -\int_{s}^{T}\int_{\mathbb{R}^{l}} U_{\gamma_{t}}(s,r,z)\tilde{N}(dr,dz), \end{aligned} $$
(1)

where \(X_{\gamma _{t}}\) is a \(\mathbb {R}^{d}\)-valued adapted càdlàg process on [t,T] with initial path \(\gamma \in D\left ([0,T],\mathbb {R}^{d}\right)\). N(ds,dx) is an integer-valued random measure and \(\tilde {N}(ds,dz)=N(ds,dz)-\nu (dz)ds\) is the compensated Poisson measure with a Lévy measure ν. For the definition of the spaces that we need for the later sections, we use the notation of Delong (2013), Levental et al. (2013), Peng and Wang (2016) and Yong (2008).

  • Let \(\mathcal {M}\left (\mathbb {R}^{l};\mathbb {R}^{m}\right)\) be the space of all measurable functions from \(\mathbb {R}^{l}\) to \(\mathbb {R}^{m}\).

  • \(L^{2}_{\mu }\left (\mathbb {R}^{l};\mathbb {R}^{m}\right)\) denotes the space of all \(\varphi \in \mathcal {M}\left (\mathbb {R}^{l};\mathbb {R}^{m}\right)\) which satisfy

    $$ \int_{\mathbb{R}^{l}}|\varphi(x)|^{2}\mu(dx)<\infty. $$
  • Let

    $$ \begin{aligned} L^{2}_{\mathcal{F}^{t}_{T}}(D([0,T],\mathbb{R}^{d});\mathbb{R}^{m}) &=\{\left.\psi:\Omega\times D\left([0,T],\mathbb{R}^{d}\right) \rightarrow\mathbb{R}^{m}\right|\psi(\cdot)\:\text{ is}\: \mathcal{F}^{t}_{T}\text{-measurable} \\& \quad\text{and}\: \mathbb{E}[|\psi(\cdot)|^{2}] <\infty\}. \end{aligned} $$
  • Similarly, we define

    $$ L^{2}_{\mathbb{F}}(t,T)=\left\{ \left.h:\Omega\times[t,T] \rightarrow\mathbb{R}^{m}\right|h(s)\:\text{is}\:\mathcal{F}^{t}_{s} \text{-measurable},\:\mathbb{E}\left[\int_{t}^{T}|h(s)|^{2}dr\right] <\infty\right\} $$

    and

    $$\begin{aligned} &L^{2}_{\mathcal{F}^{t}_{S}}(R,T) \\& =\left\{ \left. \psi:[R,T]\times\Omega \rightarrow\mathbb{R}^{m}\right| \psi(s)\:\text{ is }\:\mathcal{F}_{S}^{t}\text{-measurable}, \mathbb{E}\left[\int_{R}^{T}|\psi(u)|^{2}du\right]<\infty \right\}. \end{aligned} $$
  • We denote by S2[t,T] the space for the process Y with

    $$\begin{aligned} S^{2}[t,T]& =\left\{\left.Y:\Omega\times[t,T]\rightarrow\mathbb{R}^{m} \right| (Y(s))_{s\in[t,T]}\;\text{is a c\`{a}dl\`{a}g}\;\mathbb{F}\text{-adapted process and}\right. \\& \quad\left. \mathbb{E}\left[\sup_{s\in[t,T]}|Y(s)|^{2}\right] <\infty\right\}. \end{aligned} $$
  • For the process Z, we introduce the space

    $$\begin{aligned} M^{2}[t,T] &=\left\{ \left.Z:\Omega\times[t,T]\rightarrow\mathbb{R}^{m\times d} \right|(Z(s))_{s\in[t,T]}\;\text{is}\;\text{predictable and}\right. \\&\left. \quad\mathbb{E}\left[\int_{t}^{T}\Vert Z(s)\Vert^{2}ds\right]<\infty\right\}, \end{aligned} $$

    where we denote \(\Vert z \Vert =\sqrt {tr(zz^{*})}\) for \(z\in \mathbb {R}^{m\times d}\).

  • The space for the process U is denoted by

    $$\begin{aligned} L^{2}_{N}(t,T)= & \left\{U:\Omega\times[t,T]\times \mathbb{R}^{l}\rightarrow\mathbb{R}^{m} \left|\; (U(s,z))_{s\in[t,T],z\in\mathbb{R}^{l}} \text{ is predictable and}\right.\right. \\ &\left. \mathbb{E}\left[\int_{t}^{T}\int_{\mathbb{R}^{l}}|U(s,z)|^{2}\nu(dz)ds\right]<\infty \right\}. \end{aligned} $$
  • We define the space \(\mathbb {H}^{2}[t,T]:=S^{2}[t,T]\times M^{2}[t,T]\times L^{2}_{N}(t,T)\) and equip it with the norm

    $$\begin{aligned} & \Vert (Y,Z,U) \Vert_{\mathbb{H}^{2}[t,T]} \\& =\left(\mathbb{E}\left[\sup_{s\in[t,T]}|Y(s)|^{2} +\int_{t}^{T}\Vert Z(r)\Vert^{2}dr +\int_{t}^{T}\int_{\mathbb{R}^{l}} |U(r,z)|^{2}\nu(dz)dr\right]\right)^{\frac{1}{2}}. \end{aligned} $$
  • We denote the space for the process Z(·,·) of a solution to a BSVIE by

    $$\begin{aligned} L^{2}(S,T;M^{2}[R,T]) =& \left\{\left.Z:\Omega\times[S,T]\times[R,T]\rightarrow \mathbb{R}^{m\times d} \right|\ \text{for almost all}\:u\in[S,T] \right.\\&\left. Z(u,\cdot)\in M^{2}[R,T], \:\mathbb{E}\left[\int_{S}^{T}\int_{R}^{T} \Vert Z(u,r)\Vert^{2}dr\:du\right]< \infty\right\}. \end{aligned} $$
  • For the process U(·,·,·) of BSVIEs, we define the space

    $$\begin{aligned} &L^{2}(S,T;L^{2}_{N}(R,T)) \\& =\left\{\left.U:\Omega\times[S,T]\times[R,T]\times\mathbb{R}^{l} \rightarrow\mathbb{R}^{m}\right|\:\text{for almost all}\: u\in[S,T]\right. \\& \quad\left. U(u,\cdot,\cdot)\in L^{2}_{N}(R,T), \mathbb{E}\left[\int_{S}^{T}\int_{R}^{T}\int_{\mathbb{R}^{l}} |U(u,r,z)|^{2}\nu(dz)dr\:du\right]<\infty\right\}. \end{aligned} $$
  • We define the space \(\mathcal {H}^{2}[t,T]=S^{2}[t,T]\times L^{2}(t,T;M^{2}[t,T])\times L^{2}\left (t,T;L^{2}_{N}(t,T)\right)\) equipped with the norm

    $$\begin{aligned} &\Vert(Y(\cdot),Z(\cdot,\cdot),U(\cdot,\cdot,\cdot)) \Vert_{\mathcal{H}^{2}[t,T]} \\& =\left(\mathbb{E}\left[\sup_{s\in[t,T]}|Y(s)|^{2} +\int_{t}^{T}\int_{t}^{T}\Vert Z(s,r)\Vert^{2} dr\:ds \right.\right.\\& \quad\left.\left.+\int_{t}^{T}\int_{t}^{T}\int_{\mathbb{R}^{l}} |U(s,r,z)|^{2}\nu(dz)dr\:ds\right]\right)^{\frac{1}{2}}. \end{aligned} $$

Furthermore, we assume that the weak property of predictable representation (PR) holds. (PR) Any square integrable \(\mathbb {F}\)-martingale M has the unique representation

$$ M(s)=M(t)+\int_{t}^{s}Z(r)dW(r) +\int_{t}^{s}\int_{\mathbb{R}^{l}}U(r,x)\tilde{N}(dr,dx),\quad s\in[t,T], $$

where \((Z,U)\in M^{2}[t,T]\times L^{2}_{N}(t,T)\). For more details, see Section III.4 in Jacod and Shiryaev (2003) or Section XIII.2 in He et al. (1992).

For the reader’s convenience, we introduce some necessary definitions and results. From Levental et al. (2013), we use the following definition about path-derivatives. These derivatives are consistent with the derivatives from Dupire (2009).

Definition 1

Let \(F:D\left (\left [0,T\right ],\mathbb {R}^{d}\right)\rightarrow \mathbb {R}^{m}\) and {e i }i=1,…,d the canonical basis of \(\mathbb {R}^{d}\). For \(\gamma \in D\left ([0,T],\mathbb {R}^{d}\right)\), we denote by \(D_{i}^{1}F(\gamma ;[t,T])\), i=1,…,d, the derivative of F at γin the direction of the \(\mathbb {R}^{d}\)-valued process , which means

By

we denote the second-order derivative in the direction \(1_{\left [t,T\right ]}e_{i}\) and \(1_{\left [t,T\right ]}e_{j}\), for i,j=1,…,d. Note that the derivative \(D^{1}_{i}F(\gamma ;[t,T])\) is in \(\mathbb {R}^{m}\) and \(D^{1}F(\gamma ;[t,T])\in \mathbb {R}^{m\times d}\). Furthermore, the second-order derivative \(D_{ij}^{2}F(\gamma ;[t,T])\) is in \(\mathbb {R}^{m}\) and D2F(γ;[t,T]) is a tensor. F is said to be of class \(C^{1}\left (D\left ([0,T],\mathbb {R}^{d}\right)\right)\) if the derivatives \(D^{1}_{i}F(\gamma ;[t,T])\)exist and are continuous. The class \(C^{2}\left (D\left ([0,T],\mathbb {R}^{d}\right)\right)\) can be defined similarly.

In order to prove our main result on BSVIE in Section 3, we introduce an auxiliary path-dependent BSDE similar to that in Wang and Yong (2015) where the simpler case without jumps and no path-dependency is considered.

First, for a given path \(\gamma \in D\left ([0,T],\mathbb {R}^{d}\right)\) and t[0,T] we consider the following path-dependent BSDE with jumps

$$ \begin{aligned} Y_{\gamma_{t}}(s)=&\Phi(X_{\gamma_{t},T})+\int_{s}^{T} f(r,X_{\gamma_{t},r-},Y_{\gamma_{t}}(r-),Z_{\gamma_{t}}(r), U_{\gamma_{t}}(r))dr -\int_{s}^{T}Z_{\gamma_{t}}(r)dW(r) \\&-\int_{s}^{T}\int_{\mathbb{R}^{l}}U_{\gamma_{t}}(r,z)\tilde{N}(dr,dz), \quad s\in[t,T]. \end{aligned} $$
(2)

We assume that (A1) \(\Phi \in L^{2}_{\mathcal {F}^{t}_{T}}\left (D\left ([0,T],\mathbb {R}^{d}\right);\mathbb {R}^{m}\right)\), (A2) the generator \(f:\Omega \times [t,T]\times D\left ([0,T],\mathbb {R}^{d}\right)\times \mathbb {R}^{m}\times \mathbb {R}^{m\times d}\times \mathcal {M}(\mathbb {R}^{l};\mathbb {R}^{m})\rightarrow \mathbb {R}^{m}\) is \(\mathbb {F}\)-progressively measurable and Lipschitz continuous in the sense that

$$ \begin{aligned} & |f(\omega,s,\gamma,y,z,u) -f(\omega,s,\bar{\gamma},\bar{y},\bar{z},\bar{u})|\\ & \leq K(s)\left(\Vert\gamma-\bar{\gamma}\Vert_{D} +|y-\bar{y}|+\Vert z-\bar{z}\Vert+ \int_{\mathbb{R}^{l}}|u(z)-\bar{u}(z)| \nu(dz)\right), \end{aligned} $$

a.e. (ω,s)Ω×[t,T] for all \((\gamma,y,z,u),(\bar {\gamma },\bar {y},\bar {z},\bar {u})\in D\left ([0,T],\mathbb {R}^{d}\right)\times \mathbb {R}^{m}\times \mathbb {R}^{m\times d}\times L^{2}_{\nu }(\mathbb {R}^{l};\mathbb {R}^{m})\) and with \(K:=\sup _{s\in [t,T]}|K(s)|^{2}<\infty \), (A3) \(\mathbb {E}\left [\int _{t}^{T}|f_{0}(r,\gamma _{r-})|^{2}dr\right ]<\infty \), where f0(r,γ):=f(r,γ,0,0,0).

The following results can be found in Wang (2015) and Kromer et al. (2017). Since we do not need a comparison result for path-dependent BSDEs with jumps, it is sufficient to assume that the generator is progressively measurable and not predictable.

Theorem 1

Assume that (A1)–(A3) hold. For each \(\gamma \in D\left ([0,T],\mathbb {R}^{d}\right)\), t[0,T], BSDE (2) admits a unique adapted solution \((Y_{\gamma _{t}},Z_{\gamma _{t}},U_{\gamma _{t}})\in \mathbb {H}^{2}[t,T]\) with

$$ Y_{\gamma_{t}}(s) =\mathbb{E}\left[\left. \Phi(X_{\gamma_{t},T}) +\int_{s}^{T}f(r,X_{\gamma_{t},r-},Y_{\gamma_{t}}(r-), Z_{\gamma_{t}}(r),U_{\gamma_{t}}(r))dr\right|\mathcal{F}^{t}_{s} \right] $$
(3)

and the following stability estimate holds

$$ \left\Vert (Y_{\gamma_{t}},Z_{\gamma_{t}},U_{\gamma_{t}}) \right\Vert^{2}_{\mathbb{H}^{2}[t,T]} \leq C\left(\mathbb{E}\left[|\Phi(X_{\gamma_{t},T})|^{2}\right] + \mathbb{E}\left[\int_{t}^{T} |f_{0}(r,X_{\gamma_{t},r-})|^{2}dr\right]\right). $$
(4)

Let \(\left (\bar {f},\bar {\Phi }\right)\) also satisfy (A1)–(A3) and for any \(\bar {\gamma }\in D\left ([0,T],\mathbb {R}^{d}\right)\) let \(\left (\bar {Y}_{\bar {\gamma }_{t}}, \bar {Z}_{\bar {\gamma }_{t}}, \bar {U}_{\bar {\gamma }_{t}}\right)\in \mathbb {H}^{2}[t,T]\) be the unique adapted solution to BSDE (2) with respect to \((\bar {f},\bar {\Phi })\). Then,

$$ \begin{aligned} &\left\Vert \left(Y_{\gamma_{t}},Z_{\gamma_{t}},U_{\gamma_{t}}\right) -\left(\bar{Y}_{\bar{\gamma}_{t}}, \bar{Z}_{\bar{\gamma}_{t}}, \bar{U}_{\bar{\gamma}_{t}}\right) \right\Vert^{2}_{\mathbb{H}^{2}[t,T]} \leq C\mathbb{E}\left[|\Phi(X_{\gamma_{t},T}) -\bar{\Phi}(X_{\bar{\gamma}_{t},T})|^{2}\right.\\ &\left.+\int_{t}^{T}\left| f\left(r,X_{\gamma_{t},r-},Y_{\gamma_{t}}(r-),Z_{\gamma_{t}}(r), U_{\gamma_{t}}(r)\right) -\bar{f}\left(r,X_{\bar{\gamma}_{t},r-},Y_{\gamma_{t}}(r-),Z_{\gamma_{t}}(r), U_{\gamma_{t}}(r)\right)\right|^{2}dr\right]. \end{aligned} $$
(5)

3 Well-posedness of path-dependent BSVIEs with jumps

In this section, we look for a unique adapted solution to the path-dependent BSVIE (1). The notation in this section will follow the notation in Yong (2008) and Wang and Yong (2015).

For tR<ST we denote

$$\begin{aligned} &\Delta[R,S]=\{(r,s)\in[R,S]^{2}| R\leq s < r \leq S\}, \\&\Delta^{c}[R,S]=\{(r,s)\in[R,S]^{2}|R\leq r \leq s \leq S\}. \end{aligned} $$

In contrast to path-dependent BSDEs with jumps, the generator of a BSVIE depends on both Z(s,r), U(s,r) and Z(r,s), U(r,s). If the processes \((Y_{\gamma _{t}}(\cdot), Z_{\gamma _{t}}(\cdot,\cdot),U_{\gamma _{t}}(\cdot,\cdot,\cdot))\) satisfy a path-dependent BSVIE with jumps in the usual Itô sense, the processes \(Z_{\gamma _{t}}\) and \(U_{\gamma _{t}}\) of the solution are only unique determined on Δc[t,T], as shown in Example 1.1. in Yong (2008). There a more restrictive notion of solution is also introduced which we modify for our purpose.

Definition 2

Let t[0,T], S[t,T) and \(\gamma \in D\left ([0,T],\mathbb {R}^{d}\right)\). \((Y_{\gamma _{t}}(\cdot),Z_{\gamma _{t}}(\cdot,\cdot), U_{\gamma _{t}}(\cdot,\cdot,\cdot))\in \mathcal {H}^{2}[S,T]\) is called an adapted M-solution to BSVIE (1) on [S,T] if (1) holds in the usual Itô sense for almost all s[S,T] and, in addition, the following holds

$$ \begin{aligned} & Y_{\gamma_{t}}(u) \\ &=\mathbb{E}\left[\left. Y_{\gamma_{t}}(u)\right|\mathcal{F}_{S}^{t}\right] +\int_{S}^{u}Z_{\gamma_{t}}(u,r)dW(r) +\int_{S}^{u}\int_{\mathbb{R}^{l}} U_{\gamma_{t}}(u,r,z)\tilde{N}(dr,dz), \: \text{a.e.}\:u\in[S,T]. \end{aligned} $$
(6)

Remark 1

If \(\left (Y_{\gamma _{t}}(\cdot),Z_{\gamma _{t}}(\cdot,\cdot), U_{\gamma _{t}}(\cdot,\cdot,\cdot)\right)\in \mathcal {H}^{2}[S,T]\) is an adapted M-solution to BSVIE (1) on [S,T], it is also an adapted M-solution on \([\bar {S},T]\) for any \(\bar {S}\in (S,T)\).

For any \(\bar {S}\in (S,T)\)and almost all \(u\in [\bar {S},T]\), with Eq. (6) is

$$\mathbb{E}\left[\left.Y_{\gamma_{t}}(u)\right| \mathcal{F}_{\bar{S}}^{t}\right] =\mathbb{E}\left[\left.Y_{\gamma_{t}}(u)\right|\mathcal{F}_{S}^{t}\right] +\int_{S}^{\bar{S}} Z_{\gamma_{t}}(u,r)dW(r) +\int_{S}^{\bar{S}}\int_{\mathbb{R}^{l}} U_{\gamma_{t}}(u,r,z)\tilde{N}(dr,dz). $$

This yields for almost all \(u\in [\bar {S},T]\)

$$\begin{aligned} Y_{\gamma_{t}}(u) &=\mathbb{E}\left[\left.Y_{\gamma_{t}}(u)\right| \mathcal{F}_{\bar{S}}^{t}\right] -\int_{S}^{\bar{S}} Z_{\gamma_{t}}(u,r)dW(r) -\int_{S}^{\bar{S}}\int_{\mathbb{R}^{l}} U_{\gamma_{t}}(u,r,z)\tilde{N}(dr,dz) \\&\quad +\int_{S}^{u}Z_{\gamma_{t}}(u,r)dW(r) +\int_{S}^{u}\int_{\mathbb{R}^{l}} U_{\gamma_{t}}(u,r,z)\tilde{N}(dr,dz) \\&=\mathbb{E}\left[\left.Y_{\gamma_{t}}(u)\right| \mathcal{F}_{\bar{S}}^{t}\right] +\int_{\bar{S}}^{u} Z_{\gamma_{t}}(u,r)dW(r) +\int_{\bar{S}}^{u}\int_{\mathbb{R}^{l}} U_{\gamma_{t}}(u,r,z)\tilde{N}(dr,dz). \end{aligned} $$

As in Yong (2008), we prove existence and uniqueness of a solution to the path-dependent BSVIE with jumps (1) where the generator is independent of \(Y_{\gamma _{t}}(u)\), \(Z_{\gamma _{t}}(u,s)\) and \(U_{\gamma _{t}}(u,s,z)\) for \((s,u,z)\in \Delta [t,T]\times \mathbb {R}^{l}\). For the proof we use an auxiliary stochastic integral equation of the form

$$ \begin{aligned} \lambda_{\gamma_{t}}(s,r) =&\Phi(X_{\gamma_{t},s}) +\int_{r}^{T}h\left(s,u,X_{\gamma_{t},u-},\mu_{\gamma_{t}}(s,u), \varphi_{\gamma_{t}}(s,u)\right)du -\int_{r}^{T}\mu_{\gamma_{t}}(s,u)dW(u) \\&-\int_{r}^{T}\int_{\mathbb{R}^{l}}\varphi_{\gamma_{t}}(s,u,z) \tilde{N}(du,dz),\quad r\in\left[R,T\right],\:s\in\left[S,T\right], \end{aligned} $$
(7)

where R,S[t,T) are given, the process \(\lambda _{\gamma _{t}}(s,\cdot)\) is \(\mathbb {F}\)-adapted, and \(\mu _{\gamma _{t}}(s,\cdot),\varphi _{\gamma _{t}}(s,\cdot,\cdot)\) are predictable for all s[S,T]. For any \(s\in \left [S,T\right ]\) this stochastic integral equation (6) is a path-dependent BSDE with jumps on \(\left [R,T\right ]\). Similar to Yong (2008), we assume (B1) \(\Phi \in L^{2}_{\mathcal {F}^{t}_{T}}(D([0,T],\mathbb {R}^{d}); \mathbb {R}^{m})\). (B2)

$$ h:\Omega\times [S,T]\times[R,T]\times D([0,T],\mathbb{R}^{d}) \times \mathbb{R}^{m\times d} \times \mathcal{M}(\mathbb{R}^{l};\mathbb{R}^{m}) \rightarrow\mathbb{R}^{m} $$

is \(\mathcal {F}^{t}_{T}\otimes \mathcal {B}\left ([S,T]\times [R,T]\times D([0,T],\mathbb {R}^{d}) \times \mathbb {R}^{m\times d} \times \mathcal {M}(\mathbb {R}^{l};\mathbb {R}^{m})\right)\)-measurable such that \(r\mapsto h\left (s,r,\gamma,z,u\right)\) is \(\mathbb {F}\)-progressive measurable for all \((s,\gamma,z,u)\in [S,T]\times D\left ([0,T],\mathbb {R}^{d}\right)\times \mathbb {R}^{m\times d}\times L^{2}_{\nu }\left (\mathbb {R}^{l},\mathbb {R}^{m}\right)\). Furthermore, the generator h is Lipschitz continuous in the sense that

$$ \begin{aligned} & |h(\omega,s,r,\gamma,z,u)- h(\omega,s,r,\bar{\gamma},\bar{z},\bar{u})| \\& \leq K(s,r)\left(\Vert\gamma-\bar{\gamma}\Vert_{D} + \Vert z-\bar{z}\Vert + \int_{\mathbb{R}^{l}}\left|u(x)-\bar{u}(x)\right|\nu(dx)\right), \end{aligned} $$

a.e. (ω,s,r)Ω×[S,T]×[R,T] for all \((\gamma,z,u)\in D([0,T],\mathbb {R}^{d})\times \mathbb {R}^{m\times d}\times L^{2}_{\nu }\left (\mathbb {R}^{l},\mathbb {R}^{m}\right)\). \(K:[S,T]\times [R,T]\rightarrow [0,\infty)\) is a deterministic function such that for ε>0,

$$ \sup_{s\in[S,T]}\int_{R}^{T}|K(s,r)|^{2+\varepsilon}dr<\infty. $$

(B3) \(\mathbb {E}\left [\int _{S}^{T}\int _{R}^{T} |h_{0}(s,u,\gamma _{u-})|^{2}du \:ds\right ]<\infty \), where \(h_{0}\left (s,u,\gamma \right):=h(s,u,\gamma,0,0)\).

In the special case for R=S with given S[t,T), we define

$$\begin{aligned}Y_{\gamma_{t}}(s) & =\lambda_{\gamma_{t}}(s,s),\quad s\in[S,T],\\ Z_{\gamma_{t}}(s,u) & =\mu_{\gamma_{t}}(s,u),\quad (s,u)\in\Delta^{c}[S,T],\\ U_{\gamma_{t}}(s,u,z)&=\varphi_{\gamma_{t}}(s,u,z),\quad (s,u)\in\Delta^{c}[S,T],\:z\in\mathbb{R}^{l}. \end{aligned} $$

Similar to Yong (2008), we get for this special case of the auxiliary stochastic integral Eq. (7) the following result.

Corollary 1

Let (B1)–(B3) hold. For any \(\gamma \in D\left ([0,T],\mathbb {R}^{d}\right)\) and t[0,T] a path-dependent BSVIE with jumps of the form

$$ \begin{aligned} Y_{\gamma_{t}}(s) &= \Phi(X_{\gamma_{t},s}) +\int_{s}^{T}h\left(s,u,X_{\gamma_{t},u-},Z_{\gamma_{t}}(s,u), U_{\gamma_{t}}(s,u)\right)du -\int_{s}^{T}Z_{\gamma_{t}}(s,u)dW(u) \\&\quad-\int_{s}^{T}\int_{\mathbb{R}^{l}} U_{\gamma_{t}}(s,u,z)\tilde{N}(du,dz),\quad s\in[S,T], \end{aligned} $$
(8)

admits a unique adapted M-solution \((Y_{\gamma _{t}}(\cdot),Z_{\gamma _{t}}(\cdot,\cdot), U_{\gamma _{t}}(\cdot,\cdot,\cdot))\in \mathcal {H}^{2}[S,T]\) and the following estimates hold

$$ \begin{aligned} &\Vert (Y_{\gamma_{t}}(\cdot),Z_{\gamma_{t}}(\cdot,\cdot), U_{\gamma_{t}}(\cdot,\cdot,\cdot))\Vert^{2}_{\mathcal{H}^{2}[S,T]} \\&\leq C\left(\mathbb{E}\left[\int_{S}^{T}|\Phi(X_{\gamma_{t},s-})|^{2}ds\right] +\mathbb{E}\left[\int_{S}^{T}\int_{S}^{T} |h_{0}(s,u,X_{\gamma_{t},u-})|^{2} du\:ds\right]\right) \end{aligned} $$
(9)

and for TS small

$$ \begin{aligned} &\Vert (Y_{\gamma_{t}}(\cdot),Z_{\gamma_{t}}(\cdot,\cdot), U_{\gamma_{t}}(\cdot,\cdot,\cdot))\Vert^{2}_{\mathcal{H}^{2}[S,T]} \\&\leq C\left(\mathbb{E}\left[\int_{S}^{T}|\Phi(X_{\gamma_{t},s-})|^{2}ds\right] +\mathbb{E}\left[\int_{S}^{T}\left(\int_{S}^{T} |h_{0}(s,u,X_{\gamma_{t},u-})| du\right)^{2}ds\right]\right). \end{aligned} $$
(10)

Let \(\bar {h}\) and \(\bar {\Phi }\) also satisfy assumptions (B1)–(B3) and for a given path \(\bar {\gamma }\in D([0,T],\mathbb {R}^{d})\) let \((Y_{\gamma _{t}} (\cdot), Z_{\gamma _{t}}(\cdot,\cdot), U_{\gamma _{t}}(\cdot,\cdot,\cdot))\) be the unique adapted M-solution to BSVIE (8) with generator \(\bar {h}\) and terminal condition \(\bar {\Phi }\). For TS small enough, we get the estimate

$$ \begin{aligned} &\mathbb{E}\left[\sup_{s\in[S,T]}|Y_{\gamma_{t}}(s) -\bar{Y}_{\bar{\gamma}_{t}}(s)|^{2} +\int_{S}^{T} \Vert Z_{\gamma_{t}}(s,r)-\bar{Z}_{\bar{\gamma}_{t}}(s,r)\Vert^{2}dr \right.\\ &\left.\quad+\int_{S}^{T}\int_{\mathbb{R}^{l}} |U_{\gamma_{t}}(s,r,z) -\bar{U}_{\bar{\gamma}_{t}}(s,r,z)|^{2}\nu(dz)dr\right] \\ &\leq C\left(\mathbb{E}\left[|\Phi(X_{\gamma_{t},s})-\bar{\Phi}(X_{\bar{\gamma}_{t},s})|^{2}\right] +\mathbb{E}\left[\left(\int_{S}^{T} |h(s,u,X_{\gamma_{t},u-},Z_{\gamma_{t}}(s,u), U_{\gamma_{t}}(s,u))\right.\right.\right.\\&\quad \left.\left.\left.-\bar{h}(s,u,X_{\bar{\gamma}_{t},u-},Z_{\gamma_{t}}(s,u), U_{\gamma_{t}}(s,u))|du\right)^{2}\right]\right), \quad\text{a.e.}\;s\in[S,T]. \end{aligned} $$
(11)

Now, we show existence and uniqueness of a solution to a path-dependent BSVIE with jumps where the generator depends on processes \(Y_{\gamma _{t}}(s)\), \(Z_{\gamma _{t}}(u,s)\) and \(U_{\gamma _{t}}(u,s,z)\) for \((s,u)\in \Delta [t,T]\times \mathbb {R}^{l}\). The following assumptions are needed. (H1) \(\Phi \in L^{2}_{\mathcal {F}_{T}^{t}}\left (D\left ([0,T],\mathbb {R}^{d}\right); \mathbb {R}^{m}\right)\). (H2) The generator \(f:\Omega \times \Delta ^{c}[t,T]\times D\left ([0,T],\mathbb {R}^{d}\right)\times \mathbb {R}^{m} \times \mathbb {R}^{m\times d}\times \mathbb {R}^{m\times d} \times \mathcal {M}(\mathbb {R}^{l};\mathbb {R}^{m}) \times \mathcal {M}(\mathbb {R}^{l};\mathbb {R}^{m}) \rightarrow \mathbb {R}^{m}\) is \(\mathcal {F}_{T}^{t}\otimes \mathcal {B}(\Delta ^{c}[t,T]\times D\left ([0,T],\mathbb {R}^{d}\right) \times \mathbb {R}^{m}\times \mathbb {R}^{m\times d}\times \mathbb {R}^{m\times d} \times \mathcal {M}(\mathbb {R}^{l};\mathbb {R}^{m}) \times \mathcal {M}(\mathbb {R}^{l};\mathbb {R}^{m}))\)-measurable such that \(r\mapsto f\left (s,r,\gamma,y,z,\zeta,u,\xi \right)\) is \(\mathbb {F}\)-progressive measurable. Furthermore, generator f is Lipschitz continuous in the sense that

$$ \begin{aligned} & |f(\omega,s,r,\gamma,y,z,\zeta,u,\xi)- f(\omega,s,r,\bar{\gamma},\bar{y},\bar{z},\bar{\zeta}, \bar{u},\bar{\xi})| \\ & \leq K(s,r)\left(\Vert\gamma-\bar{\gamma}\Vert_{D}+|y-\bar{y}| + \Vert z-\bar{z}\Vert + \Vert \zeta-\bar{\zeta}\Vert +\int_{\mathbb{R}^{l}}\left|u(x)-\bar{u}(x)\right|\nu(dx) \right.\\ & \quad\left. +\int_{\mathbb{R}^{l}}\left|\xi(x)-\bar{\xi}(x)\right|\nu(dx) \right), \end{aligned} $$

a.e. (ω,s,r)Ω×Δc[t,T] for all \( (\gamma,y,z,\zeta,u,\xi),(\bar {\gamma },\bar {y},\bar {z},\bar {\zeta }, \bar {u},\bar {\xi })\in D\left ([0,T],\mathbb {R}^{d}\right)\times \mathbb {R}^{m}\times \mathbb {R}^{m\times d}\times \mathbb {R}^{m\times d}\times L^{2}_{\nu }\left (\mathbb {R}^{l};\mathbb {R}^{m}\right)\times L^{2}_{\nu }\left (\mathbb {R}^{l};\mathbb {R}^{m}\right)\). Moreover, \(K:\Delta ^{c}[t,T]\rightarrow [0,\infty)\) is a deterministic function such that for ε>0,

$$ \sup_{s\in[t,T]}\int_{s}^{T}|K(s,r)|^{2+\varepsilon}dr<\infty. $$

(H3) \(\mathbb {E}\left [\int _{t}^{T}\int _{t}^{T} |f_{0}(s,r,\gamma _{r-})|^{2}dr\:ds\right ]<\infty \), where f0(s,r,γr)=f(s,r,γr,0,0,0,0,0).

Theorem 2

Let (H1)–(H3) hold. Then, for any given path \(\gamma \in D([0,T],\mathbb {R}^{d})\), t[0,T], BSVIE (1) admits a unique adapted M-solution \((Y_{\gamma _{t}}(\cdot),Z_{\gamma _{t}}(\cdot,\cdot), U_{\gamma _{t}}(\cdot,\cdot,\cdot))\in \mathcal {H}^{2}[t,T]\). Moreover, the following estimate holds for all S[t,T]

$$ \begin{aligned} &\Vert (Y_{\gamma_{t}}(\cdot),Z_{\gamma_{t}}(\cdot,\cdot), U_{\gamma_{t}}(\cdot,\cdot,\cdot))\Vert^{2}_{\mathcal{H}^{2}[S,T]} \\&\leq C\left(\mathbb{E}\left[\int_{S}^{T}|\Phi(X_{\gamma_{t},s-})|^{2}ds\right] +\mathbb{E}\left[\int_{S}^{T}\int_{S}^{T} |f_{0}(s,r,X_{\gamma_{t},r-})|^{2} dr\:ds\right]\right). \end{aligned} $$
(12)

Proof

We use the techniques of Yong (2008) and split this proof into four steps.

First step: We show the existence and uniqueness of the adapted M-solution to BSVIE (1) on [S,T] for some S[t,T) and given \(\gamma \in D([0,T],\mathbb {R}^{d})\), t[0,T]. For any S[t,T) we define the space

$$\begin{aligned} \mathbb{M}^{2}[S,T] =&\left\{\left. (y(\cdot),z(\cdot,\cdot),u(\cdot,\cdot,\cdot)) \in\mathcal{H}^{2}[S,T]\right| y(s)=\mathbb{E}\left[\left.y(s)\right|\mathcal{F}^{t}_{S}\right]\right.\\ &\left.+\int_{S}^{s}z(s,r)dW(r)+\int_{S}^{s}\int_{\mathbb{R}^{l}}u(s,r,z)\tilde{N}(dr,dz),s\in[S,T] \right\}. \end{aligned} $$

For \((y(\cdot),z(\cdot,\cdot),u(\cdot,\cdot,\cdot))\in \mathbb {M}^{2}[S,T]\) we get

$$\begin{aligned} &\mathbb{E}\left[\sup_{s\in[S,T]}|y(s)|^{2} +\int_{S}^{T}\int_{s}^{T}\Vert z(s,r)\Vert^{2} dr\:ds +\int_{S}^{T}\int_{s}^{T}\int_{\mathbb{R}^{l}}|u(s,r,z)|^{2} \nu(dz)dr\:ds\right] \\& \leq\mathbb{E}\left[\sup_{s\in[S,T]}|y(s)|^{2} +\int_{S}^{T}\int_{s}^{T}\Vert z(s,r)\Vert^{2} dr\:ds +\int_{S}^{T}\int_{s}^{T}\int_{\mathbb{R}^{l}}|u(s,r,z)|^{2} \nu(dz)dr\:ds\right] \\&\quad +\mathbb{E}\left[ \int_{S}^{T}\int_{S}^{s}\Vert z(s,r)\Vert^{2} dr\:ds +\int_{S}^{T}\int_{S}^{s}\int_{\mathbb{R}^{l}}|u(s,r,z)|^{2} \nu(dz)dr\:ds\right] \\&\leq C\mathbb{E}\left[\sup_{s\in[S,T]}|y(s)|^{2} +\int_{S}^{T}\int_{s}^{T}\Vert z(s,r)\Vert^{2} dr\:ds +\int_{S}^{T}\int_{s}^{T}\int_{\mathbb{R}^{l}}|u(s,r,z)|^{2} \nu(dz)dr\:ds\right]. \end{aligned} $$

Thus,

$$\begin{aligned} &\Vert (y(\cdot),z(\cdot,\cdot),u(\cdot,\cdot,\cdot)) \Vert^{2}_{\mathbb{M}^{2}[S,T]} \\&=\mathbb{E}\left[\sup_{s\in[S,T]}|y(s)|^{2} +\int_{S}^{T}\int_{s}^{T}\Vert z(s,r)\Vert^{2} dr\:ds +\int_{S}^{T}\int_{s}^{T}\int_{\mathbb{R}^{l}}|u(s,r,z)|^{2} \nu(dz)dr\:ds\right] \end{aligned} $$

is an equivalent norm for \(\mathbb {M}^{2}[S,T]\). Now, we consider the following path-dependent BSVIE with jumps

$$ \begin{aligned} Y_{\gamma_{t}}(s)=&\Phi(X_{\gamma_{t},s}) +\int_{s}^{T}f\left(s,r,X_{\gamma_{t},r-}, y_{\gamma_{t}}(r-),Z_{\gamma_{t}}(s,r), z_{\gamma_{t}}(r,s), U_{\gamma_{t}}(s,r),u_{\gamma_{t}}(r,s)\right)dr\\ &-\int_{s}^{T}Z_{\gamma_{t}}(s,r)dW(r)-\int_{s}^{T}\int_{\mathbb{R}^{l}}U_{\gamma_{t}}(s,r,z) \tilde{N}(dr,dz),\quad s\in[S,T], \end{aligned} $$
(13)

where \(\left (y_{\gamma _{t}}(\cdot),z_{\gamma _{t}}(\cdot,\cdot), u_{\gamma _{t}}(\cdot,\cdot,\cdot)\right)\in \mathbb {M}^{2}[S,T]\) are given. Corollary 1 yields that this BSVIE admits a unique adapted M-solution \(\left (Y_{\gamma _{t}}(\cdot),Z_{\gamma _{t}}(\cdot,\cdot), U_{\gamma _{t}}(\cdot,\cdot,\cdot)\right)\in \mathcal {H}^{2}[S,T]\). With Eq. (10) the following estimate holds for TS small

$$\begin{aligned} &\mathbb{E}\left[\sup_{s\in[S,T]}|Y_{\gamma_{t}}(s)|^{2} +\int_{s}^{T}\Vert Z_{\gamma_{t}}(s,r)\Vert^{2}dr +\int_{s}^{T}\int_{\mathbb{R}^{l}} |U_{\gamma_{t}}(s,r,z)|^{2}\nu(dz)dr\right] \\&\leq C\left(\mathbb{E}\left[|\Phi(X_{\gamma_{t},s})|^{2}\right] +2\mathbb{E}\left[\left(\int_{S}^{T} |f_{0}(s,r,X_{\gamma_{t},r-})|dr\right)^{2}\right]\right. \\&\left.+2\mathbb{E}\left[\left(\int_{S}^{T}K(s,r) \left(|y_{\gamma_{t}}(r-)| +\Vert z_{\gamma_{t}}(r,s)\Vert +\int_{\mathbb{R}^{l}} |u_{\gamma_{t}}(r,s,z)|\nu(dz)\right)dr\right)^{2}\right]\right), \\&\text{a.e.}\:s\in[S,T]. \end{aligned} $$

Hölder’s inequality yields

$$ \begin{aligned} &\mathbb{E}\left[\sup_{s\in[S,T]}|Y_{\gamma_{t}}(s)|^{2} +\int_{s}^{T}\Vert Z_{\gamma_{t}}(s,r)\Vert^{2}dr +\int_{s}^{T}\int_{\mathbb{R}^{l}} |U_{\gamma_{t}}(s,r,z)|^{2}\nu(dz)dr\right]\\ &\leq 2C\left(\mathbb{E}\left[|\Phi(X_{\gamma_{t},s})|^{2}\right]+\mathbb{E}\left[\left(\int_{S}^{T} |f_{0}(s,r,X_{\gamma_{t},r-})|dr\right)^{2}\right] +\mathbb{E}\left[\left(\int_{S}^{T}|K(s,r)|^{2+\varepsilon}dr \right)^{\frac{2}{2+\varepsilon}} \right.\right.\\ &\left.\left.\cdot(T-S)^{\frac{\varepsilon}{2+\varepsilon}} \left(\sup_{s\in[S,T]}|y_{\gamma_{t}}(s)|^{2} +\int_{S}^{T}\Vert z_{\gamma_{t}}(r,s)\Vert^{2}dr +\int_{S}^{T}\int_{\mathbb{R}^{l}} |u_{\gamma_{t}}(r,s,z)|^{2}\nu(dz)dr\right)\right] \right). \end{aligned} $$
(14)

It is

$$ \begin{aligned} &\mathbb{E}\left[\sup_{s\in[S,T]}|y_{\gamma_{t}}(s)|^{2} +\int_{S}^{T}\int_{s}^{T} \Vert z_{\gamma_{t}}(r,s)\Vert^{2}dr\:ds +\int_{S}^{T}\int_{s}^{T}\int_{\mathbb{R}^{l}} |u_{\gamma_{t}}(r,s,z)|^{2}\nu(dz)dr\:ds\right] \\&\leq \mathbb{E}\left[\sup_{s\in[S,T]}|y_{\gamma_{t}}(s)|^{2} +\int_{S}^{T}\int_{S}^{T} \Vert z_{\gamma_{t}}(r,s)\Vert^{2}ds\: dr +\int_{S}^{T}\int_{S}^{T}\int_{\mathbb{R}^{l}} |u_{\gamma_{t}}(r,s,z)|^{2}\nu(dz)ds\:dr\right] \\&\leq C\mathbb{E}\left[\sup_{s\in[S,T]}|y_{\gamma_{t}}(s)|^{2} +\int_{S}^{T}\int_{r}^{T}\Vert z_{\gamma_{t}}(r,s)\Vert^{2} ds\:dr\right.\\ &\quad\left.+\int_{S}^{T}\int_{r}^{T}\int_{\mathbb{R}^{l}}|u_{\gamma_{t}}(r,s,z)|^{2} \nu(dz)ds\:dr\right] \end{aligned} $$
(15)

and this implies

$$\begin{aligned} &\Vert (Y_{\gamma_{t}}(\cdot),Z_{\gamma_{t}}(\cdot,\cdot), U_{\gamma_{t}}(\cdot,\cdot,\cdot))\Vert^{2}_{\mathbb{M}^{2}[S,T]} \\&\leq C\left\{\mathbb{E}\left[\int_{S}^{T} |\Phi(X_{\gamma_{t},s-})|^{2}ds\right] +\mathbb{E}\left[\int_{S}^{T}\int_{S}^{T}|f_{0}(s,r,X_{\gamma_{t},r-})|^{2}dr\:ds\right]\right.\\ &\left.\quad+\Vert (y_{\gamma_{t}}(\cdot),z_{\gamma_{t}}(\cdot,\cdot), u_{\gamma_{t}}(\cdot,\cdot,\cdot))\Vert_{\mathbb{M}^{2}[S,T]}^{2} \right\}. \end{aligned} $$

Therefore, \((Y_{\gamma _{t}}(\cdot),Z_{\gamma _{t}}(\cdot,\cdot), U_{\gamma _{t}}(\cdot,\cdot,\cdot))\) is in \(\mathbb {M}^{2}[S,T]\). Now, we define a map \(\Theta ^{t,\gamma }:\mathbb {M}^{2}[S,T]\rightarrow \mathbb {M}^{2}[S,T]\) by

$$\Theta^{t,\gamma}\left(y_{\gamma_{t}}(\cdot),z_{\gamma_{t}}(\cdot,\cdot), u_{\gamma_{t}}(\cdot,\cdot,\cdot)\right) =\left(Y_{\gamma_{t}}(\cdot),Z_{\gamma_{t}}(\cdot,\cdot), U_{\gamma_{t}}(\cdot,\cdot,\cdot)\right) $$

for all \((y_{\gamma _{t}}(\cdot),z_{\gamma _{t}}(\cdot,\cdot), u_{\gamma _{t}}(\cdot,\cdot,\cdot))\in \mathbb {M}^{2}[S,T]\). Furthermore, let \(\Theta ^{t,\gamma }(\bar {y}_{\gamma _{t}}(\cdot),\bar {z}_{\gamma _{t}}(\cdot,\cdot), \bar {u}_{\gamma _{t}}(\cdot,\cdot,\cdot)) =(\bar {Y}_{\gamma _{t}}(\cdot),\bar {Z}_{\gamma _{t}}(\cdot,\cdot), \bar {U}_{\gamma _{t}}(\cdot,\cdot,\cdot)) \) for \((\bar {y}_{\gamma _{t}}(\cdot), \bar {z}_{\gamma _{t}}(\cdot,\cdot), \bar {u}_{\gamma _{t}}(\cdot,\cdot,\cdot))\in \mathbb {M}^{2}[S,T]\). Corollary 1 yields for TS small

$$\begin{aligned} &\mathbb{E}\left[\sup_{s\in[S,T]}|Y_{\gamma_{t}}(s) -\bar{Y}_{\gamma_{t}}(s)|^{2} +\int_{s}^{T} \Vert Z_{\gamma_{t}}(s,r)-\bar{Z}_{\gamma_{t}}(s,r)\Vert^{2}dr\right.\\ &\left.+\int_{s}^{T}\int_{\mathbb{R}^{l}} |U_{\gamma_{t}}(s,r,z) -\bar{U}_{\gamma_{t}}(s,r,z)|^{2}\nu(dz)dr\right]\\ &\leq C\mathbb{E}\left[\left(\int_{S}^{T} K(s,r)\left(|y_{\gamma_{t}}(r-)-\bar{y}_{\gamma_{t}}(r-)| +\Vert z_{\gamma_{t}}(r,s)-\bar{z}_{\gamma_{t}}(r,s)\Vert\right.\right.\right.\\ &\quad\left.\left.\left.+\int_{\mathbb{R}^{l}} |u_{\gamma_{t}}(r,s,z)-\bar{u}_{\gamma_{t}}(r,s,z)|\nu(dz)\right)dr \right)^{2} \right] \\&\text{a.e.}\:s\in[S,T]. \end{aligned} $$

Similar to Eq. (15), we get with Hölder’s inequality and ε>0 for s[S,T]

$$ \begin{aligned} &\Vert \left(\Theta^{t,\gamma}\left(y_{\gamma_{t}}(\cdot), z_{\gamma_{t}}(\cdot,\cdot),u_{\gamma_{t}}(\cdot,\cdot,\cdot) \right) -\Theta^{t,\gamma}\left(\bar{y}_{\gamma_{t}}(\cdot), \bar{z}_{\gamma_{t}}(\cdot,\cdot), \bar{u}_{\gamma_{t}}(\cdot,\cdot,\cdot)\right) \right)\Vert^{2}_{\mathbb{M}^{2}[S,T]} \\ &\leq C(T-S)^{\frac{\varepsilon}{2+\varepsilon}}\mathbb{E}\left[(T-S)\sup_{s\in[S,T]} |y_{\gamma_{t}}(s)-\bar{y}_{\gamma_{t}}(s)|^{2} +\int_{S}^{T}\int_{S}^{T}\Vert z_{\gamma_{t}}(r,s)-\bar{z}_{\gamma_{t}}(r,s)\Vert^{2} dr\:ds \right. \\ &\left.\quad+\int_{S}^{T}\int_{S}^{T}\int_{\mathbb{R}^{l}} |u_{\gamma_{t}}(r,s,z)-\bar{u}_{\gamma_{t}}(r,s,z)|^{2}\nu(dz)dr\:ds \right.] \\&\leq C(T-S)^{\frac{\varepsilon}{2+\varepsilon}}\mathbb{E}\left[ (T-S)\sup_{s\in[S,T]}|y_{\gamma_{t}}(s)-\bar{y}_{\gamma_{t}}(s)|^{2} +\int_{S}^{T}\int_{r}^{T}\Vert z_{\gamma_{t}}(r,s)-\bar{z}_{\gamma_{t}}(r,s)\Vert^{2} ds\: dr \right.\\&\quad \left.+\int_{S}^{T}\int_{r}^{T}\int_{\mathbb{R}^{l}} |u_{\gamma_{t}}(r,s,z)-\bar{u}_{\gamma_{t}}(r,s,z)|^{2} \nu(dz)ds\:dr\right] \\&\leq C(T-S)^{\frac{\varepsilon}{2+\varepsilon}} \Vert \left(y_{\gamma_{t}}(\cdot),z_{\gamma_{t}}(\cdot,\cdot), u_{\gamma_{t}}(\cdot,\cdot,\cdot)\right) -\left(\bar{y}_{\gamma_{t}}(\cdot), \bar{z}_{\gamma_{t}}(\cdot,\cdot), \bar{u}_{\gamma_{t}}(\cdot,\cdot,\cdot)\right) \Vert^{2}_{\mathbb{M}^{2}[S,T]}. \end{aligned} $$
(16)

For TS>0 small enough, the map \(\Theta ^{t,\gamma }: \mathbb {M}^{2}[S,T] \rightarrow \mathbb {M}^{2}[S,T] \) is a contraction mapping and it admits a unique fixed point \(\left (Y_{\gamma _{t}}(\cdot), Z_{\gamma _{t}}(\cdot), U_{\gamma _{t}}(\cdot,\cdot,\cdot)\right)\in \mathbb {M}^{2}[S,T]\) which is the unique adapted M-solution of BSVIE (1) over [S,T]. Moreover, Eq. (14) yields the following estimate

$$ \begin{aligned} &\mathbb{E}\left[\sup_{s\in[S,T]}\left|Y_{\gamma_{t}}(s)\right|^{2} +\int_{S}^{T}\int_{S}^{T}\Vert Z_{\gamma_{t}}(s,r)\Vert^{2}dr\: ds +\int_{S}^{T}\int_{S}^{T}\int_{\mathbb{R}^{l}} \left|U_{\gamma_{t}}(s,r,z)\right|^{2}\nu(dz)dr\:ds\right] \\&\leq C\left(\mathbb{E}\left[\int_{S}^{T} |\Phi(X_{\gamma_{t},s-})|^{2}ds\right] +\mathbb{E}\left[\int_{S}^{T}\int_{S}^{T} |f_{0}(s,r,X_{\gamma_{t},r-})|^{2} dr\:ds\right]\right). \end{aligned} $$
(17)

Second step: Let R[t,S). Now, we want to determine \((Z_{\gamma _{t}}(s,r), U_{\gamma _{t}}(s,r,z))\) for (s,r)[S,T]×[R,S] and \(z\in \mathbb {R}^{l}\). Since \(Y_{\gamma _{t}}\in S^{2}[S,T]\), the process \(M(r)=\mathbb {E}\left [\left. Y_{\gamma _{t}}(s)\right |\mathcal {F}^{t}_{r}\right ]\), s[S,T], is a square integrable martingale and the martingale representation (PR) yields

$$\mathbb{E}\left[\left.Y_{\gamma_{t}}(s)\right| \mathcal{F}^{t}_{S}\right]= \mathbb{E}\left[\left.Y_{\gamma_{t}}(s)\right| \mathcal{F}^{t}_{R}\right]+\int_{R}^{S} Z_{\gamma_{t}}(s,r)dW(r) +\int_{R}^{S}\int_{\mathbb{R}^{l}} U_{\gamma_{t}}(s,r,z)\tilde{N} (dr,dz). $$

Furthermore, we obtain with the Burkholder–Davis–Gundy and Doob’s inequality

$$\begin{aligned} &\mathbb{E}\left[\int_{R}^{S} \Vert Z_{\gamma_{t}}(s,r)\Vert^{2}dr +\int_{R}^{S}\int_{\mathbb{R}^{l}} |U_{\gamma_{t}}(s,r,z)|^{2} \nu(dz)dr\right] \leq 4C \mathbb{E}\left[\left|Y_{\gamma_{t}}(s)\right|^{2}\right]. \end{aligned} $$

It follows with Eq. (17)

$$ \begin{aligned} &\mathbb{E}\left[ \sup_{s\in[S,T]}\left|Y_{\gamma_{t}}(s)\right|^{2} +\int_{S}^{T}\int_{R}^{S} \Vert Z_{\gamma_{t}}(s,r)\Vert^{2}dr\:ds +\int_{S}^{T}\int_{R}^{S}\int_{\mathbb{R}^{l}} |U_{\gamma_{t}}(s,r,z)|^{2} \nu(dz)dr\:ds\right] \\&\leq C\left(\mathbb{E}\left[ \int_{S}^{T}|\Phi(X_{\gamma_{t},s-})|^{2}ds\right] +\mathbb{E}\left[\int_{S}^{T}\int_{S}^{T} |f_{0}(s,r,X_{\gamma_{t},r-})|^{2} dr\:ds\right]\right). \end{aligned} $$
(18)

Combining this result and the result of the first step, we get the existence and uniqueness of the adapted M-solution \((Y_{\gamma _{t}}(\cdot),Z_{\gamma _{t}}(\cdot,\cdot), U_{\gamma _{t}}(\cdot,\cdot,\cdot))\) to BSVIE (1) on [S,T]×[R,T]. Estimates (17) and (18) provide the following estimate

$$\begin{aligned} &\mathbb{E}\left[ \sup_{s\in[S,T]}\left|Y_{\gamma_{t}}(s)\right|^{2} +\int_{S}^{T}\int_{R}^{T} \Vert Z_{\gamma_{t}}(s,r)\Vert^{2}dr\:ds +\int_{S}^{T}\int_{R}^{T}\int_{\mathbb{R}^{l}} |U_{\gamma_{t}}(s,r,z)|^{2} \nu(dz)dr\:ds\right] \\&\leq C\left(\mathbb{E}\left[\int_{S}^{T}|\Phi(X_{\gamma_{t},s-})|^{2}ds \right] +\mathbb{E}\left[\int_{S}^{T}\int_{S}^{T} |f_{0}(s,r,X_{\gamma_{t},r-})|^{2} dr\:ds\right]\right). \end{aligned} $$

Third step: Now, we consider a stochastic integral equation of the form

$$ \begin{aligned} \psi^{S}_{\gamma_{t}}(r)&=\Phi(X_{\gamma_{t},r})+\int_{S}^{T} f^{S}\left(r,s,X_{\gamma_{t},s-},Z_{\gamma_{t}}(r,s), U_{\gamma_{t}}(r,s)\right)ds \\&\quad-\int_{S}^{T}Z_{\gamma_{t}}(r,s)dW(s)-\int_{S}^{T}\int_{\mathbb{R}^{l}}U_{\gamma_{t}}(r,s,z) \tilde{N}(ds,dz), \quad r\in [R,S], \end{aligned} $$
(19)

where we define \(f^{S}:[R,S]\times [S,T]\times D([0,T],\mathbb {R}^{d})\times \mathbb {R}^{m\times d} \times \mathcal {M}(\mathbb {R}^{l},\mathbb {R}^{m})\rightarrow \mathbb {R}^{m}\) by

$$f^{S}\left(r,s,X_{\gamma_{t},s-},z,u\right) =f\left(r,s,X_{\gamma_{t},s-},Y_{\gamma_{t}}(s-),z,Z_{\gamma_{t}}(s,r), u,U_{\gamma_{t}}(s,r)\right). $$

Similar techniques to those in Yong (2008) can be used to prove an existence and uniqueness result to solutions

$$\left(\psi_{\gamma_{t}}^{S}(\cdot),Z_{\gamma_{t}}(\cdot,\cdot), U_{\gamma_{t}}(\cdot,\cdot,\cdot)\right)\in L^{2}_{\mathcal{F}_{S}^{t}}(R,S) \times L^{2}(R,S; M^{2}[S,T]) \times L^{2}\left(R,S;L^{2}_{N}(S,T)\right) $$

of the stochastic integral equation (19). Thus, \(Z_{\gamma _{t}}(r,s)\) and \(U_{\gamma _{t}}(r,s,x)\) are unique determined for (r,s)[R,S]×[S,T], \(x\in \mathbb {R}^{l}\). For TS small the following estimate holds

$$\begin{aligned} &\mathbb{E}\left[|\psi_{\gamma_{t}}^{S}(r)|^{2} +\int_{S}^{T}\Vert Z_{\gamma_{t}}(r,s)\Vert^{2}ds +\int_{S}^{T}\int_{\mathbb{R}^{l}} |U_{\gamma_{t}}(r,s,z)|^{2}\nu(dz)ds\right] \\&\leq C \left(\mathbb{E}\left[|\Phi(X_{\gamma_{t},r})|^{2}\right] +\mathbb{E}\left[\left(\int_{S}^{T} |f^{S}_{0}(r,s,X_{\gamma_{t},s-})|ds\right)^{2}\right]\right) \\&\leq C \left(\mathbb{E}\left[|\Phi(X_{\gamma_{t},r})|^{2}\right] +\mathbb{E}\left[\left(\int_{S}^{T} |f_{0}(r,s,X_{\gamma_{t},s-})|ds\right)^{2} \right]\right. \\&\quad +\left(T-S\right)^{\frac{\varepsilon}{2+\varepsilon}} K\mathbb{E}\left[ \sup_{s\in[S,T]}|Y_{\gamma_{t}}(s)|^{2}+\int_{S}^{T} \Vert Z_{\gamma_{t}}(s,r)\Vert^{2} ds\right. \\&\left.\left.\quad+\int_{S}^{T} \int_{\mathbb{R}^{l}}|U_{\gamma_{t}}(s,r,z)|^{2}\nu(dz)ds \right] \right) \\&\leq C \left(\mathbb{E}\left[|\Phi(X_{\gamma_{t},r})|^{2}\right] +\mathbb{E}\left[\left(\int_{S}^{T} |f_{0}(r,s,X_{\gamma_{t},s-})|ds\right)^{2} \right]\right.\\&\quad\left. +\mathbb{E}\left[\sup_{s\in[S,T]}|Y_{\gamma_{t}}(s)|^{2}+ \int_{S}^{T} \Vert Z_{\gamma_{t}}(s,r)\Vert^{2} ds +\int_{S}^{T}\int_{\mathbb{R}^{l}}|U_{\gamma_{t}}(s,r,z)|^{2}\nu(dz)ds \right] \right), \\&\text{a.e.}\;r\in[R,S]. \end{aligned} $$

Together with estimate (18), this yields

$$ \begin{aligned} &\mathbb{E}\left[\int_{R}^{S}|\psi_{\gamma_{t}}^{S}(r-)|^{2}dr +\int_{R}^{S}\int_{S}^{T} \Vert Z_{\gamma_{t}}(r,s)\Vert^{2}ds\:dr +\int_{R}^{S}\int_{S}^{T}\int_{\mathbb{R}^{l}} |U_{\gamma_{t}}(r,s,z)|^{2}\nu(dz)ds\:dr\right] \\&\leq C\left(\mathbb{E}\left[\int_{S}^{T} |\Phi(X_{\gamma_{t},r-})|^{2}dr\right] +\mathbb{E}\left[\int_{S}^{T}\int_{S}^{T} |f_{0}(r,s,X_{\gamma_{t},s-})|^{2} ds\:dr\right]\right) \end{aligned} $$
(20)

and for (r,s)[R,S]×[S,T], the processes \(Z_{\gamma _{t}}(r,s)\) and \(U_{\gamma _{t}}(r,s,x)\) are uniquely determined.

Fourth step: In the last step, we prove that the path-dependent BSVIE (1) admits a unique adapted M-solution on [R,S]. For r[R,S] we observe a path-dependent BSVIE with jumps of the form

$$ \begin{aligned} Y_{\gamma_{t}}(r)= &\psi^{S}_{\gamma_{t}}(r) +\int_{r}^{S}f\left(r,s,X_{\gamma_{t},s-},Y_{\gamma_{t}}(s-), Z_{\gamma_{t}}(r,s),Z_{\gamma_{t}}(s,r), U_{\gamma_{t}}(r,s),U_{\gamma_{t}}(s,r)\right)ds \\&-\int_{r}^{S}Z_{\gamma_{t}}(r,s)dW(s) -\int_{r}^{S}\int_{\mathbb{R}^{l}}U_{\gamma_{t}}(r,s,z) \tilde{N}(ds,dz),\quad r\in[R,S], \end{aligned} $$
(21)

where \(\psi _{\gamma _{t}}^{S}(\cdot)\in L^{2}_{\mathcal {F}^{t}_{S}}(R,S)\). Similarly to the first step, we can show that for SR small enough BSVIE (21) admits a unique adapted M-solution \((Y_{\gamma _{t}}(\cdot),Z_{\gamma _{t}}(\cdot,\cdot), U_{\gamma _{t}}(\cdot,\cdot,\cdot))\in \mathcal {H}^{2}[R,S]\) and it is

$$ \begin{aligned} &\mathbb{E}\left[\sup_{r\in[R,S]}\left|Y_{\gamma_{t}}(r)\right|^{2} +\int_{R}^{S}\int_{R}^{S}\Vert Z_{\gamma_{t}}(r,s)\Vert^{2}ds\: dr +\int_{R}^{S}\int_{R}^{S}\int_{\mathbb{R}^{l}} \left|U_{\gamma_{t}}(r,s,z)\right|^{2}\nu(dz)ds\:dr\right] \\&\leq C\left(\mathbb{E}\left[\int_{R}^{S} |\Phi(X_{\gamma_{t},r-})|^{2}dr\right] +\mathbb{E}\left[\int_{R}^{S}\int_{R}^{S} |f_{0}(r,s,X_{\gamma_{t},s-})|^{2} ds\:dr\right]\right). \end{aligned} $$
(22)

Insert Eq. (19) to get for r[R,S]

$$\begin{aligned} Y_{\gamma_{t}}(r) &= \psi^{S}_{\gamma_{t}}(r)+\int_{r}^{S} f\left(r,s,X_{\gamma_{t},s-}, Y_{\gamma_{t}}(s-),Z_{\gamma_{t}}(r,s),Z_{\gamma_{t}}(s,r), U_{\gamma_{t}}(r,s),U_{\gamma_{t}}(s,r)\right)ds \\&\quad -\int_{r}^{S} Z_{\gamma_{t}}(r,s)dW(s) -\int_{r}^{S}\int_{\mathbb{R}^{l}} U_{\gamma_{t}}(r,s,z)\tilde{N}(ds,dz) \\&=\Phi(X_{\gamma_{t},r})+\int_{r}^{T} f\left(r,s,X_{\gamma_{t},s-}, Y_{\gamma_{t}}(s-),Z_{\gamma_{t}}(r,s),Z_{\gamma_{t}}(s,r), U_{\gamma_{t}}(r,s),U_{\gamma_{t}}(s,r)\right)ds \\&\quad -\int_{r}^{T}Z_{\gamma_{t}}(r,s)dW(s) -\int_{r}^{T}\int_{\mathbb{R}^{l}}U_{\gamma_{t}}(r,s,z) \tilde{N}(ds,dz). \end{aligned} $$

Altogether, we determined the unique adapted M-solution to BSVIE (1) on [R,T] and the following estimate holds

$$\begin{aligned} &\mathbb{E}\left[\sup_{r\in[R,T]}|Y_{\gamma_{t}}(r)|^{2} +\int_{R}^{T}\int_{R}^{T} \Vert Z_{\gamma_{t}}(r,s)\Vert^{2} ds\:dr +\int_{R}^{T}\int_{R}^{T}\int_{\mathbb{R}^{l}} |U_{\gamma_{t}}(r,s,z)|^{2}\nu(dz)ds\:dr\right] \\&\leq C \mathbb{E}\left[\int_{R}^{T}|\Phi(X_{\gamma_{t},s-})|ds +\int_{R}^{T}\int_{R}^{T} \left|f_{0}\left(r,s,X_{\gamma_{t},s-}\right)\right|^{2}ds\:dr\right]. \end{aligned} $$

If we consider small segments of the rectangle [t,T]2 of the form [R k ,T]2, k=0,…,n, \(n\in \mathbb {N}\), R0=S, R1=R and R n =t, we get the existence and uniqueness of an adapted M-solution to BSVIE (1) over all rectangles [R k ,T]2. We use the same steps as before to show that BSVIE (1) has a unique adapted M-solution over [R k ,T] since we know that BSVIE (1) admits a unique adapted M-solution over [Rk−1,T]. □

The following result of the stability estimate will be useful for the next section. We consider a path-dependent BSVIE with jumps where the last components of the generator f enter as integrals with respect to Lévy measure ν

$$ \begin{aligned} & Y_{\gamma_{t}}(s) \\&=\Phi(X_{\gamma_{t},s}) +\int_{s}^{T} f\left(s,r,X_{\gamma_{t},r-},Y_{\gamma_{t}}(r-),Z_{\gamma_{t}}(s,r), Z_{\gamma_{t}}(r,s),\right. \\&\left.\int_{\mathbb{R}^{l}}U_{\gamma_{t}}(s,r,z)\delta(s,r,z)\nu(dz), \int_{\mathbb{R}^{l}}U_{\gamma_{t}}(r,s,z)\delta(r,s,z)\nu(dz) \right)dr \\&\quad-\int_{s}^{T}Z_{\gamma_{t}}(s,r)dW(r) -\int_{s}^{T}\int_{\mathbb{R}^{l}}U_{\gamma_{t}}(s,r,z)\tilde{N}(dr,dz), \quad s\in[t,T], \end{aligned} $$
(23)

where we assume (H’1) \(\Phi \in L^{2}_{\mathcal {F}_{T}^{t}}(D([0,T],\mathbb {R}^{d}); \mathbb {R}^{m})\). (H’2)

$$f:\Omega\times \Delta^{c}[t,T]\times D([0,T],\mathbb{R}^{d})\times\mathbb{R}^{m} \times \mathbb{R}^{m\times d}\times \mathbb{R}^{m\times d} \times \mathbb{R}^{m} \times \mathbb{R}^{m} \rightarrow\mathbb{R}^{m} $$

is \(\mathcal {F}_{T}^{t}\otimes \mathcal {B}\left (\Delta ^{c}[t,T]\times D([0,T],\mathbb {R}^{d})\times \mathbb {R}^{m} \times \mathbb {R}^{m\times d}\times \mathbb {R}^{m\times d} \times \mathbb {R}^{m} \times \mathbb {R}^{m}\right)\)-measurable such that

$$ r\mapsto f\left(s,r,\gamma,y,z,\zeta,u,\xi\right) $$

is \(\mathbb {F}\)-progressive measurable. Furthermore, generator f is Lipschitz continuous in the sense that

$$ \begin{aligned} & |f(\omega,s,r,\gamma,y,z,\zeta,u,\xi)- f(\omega,s,r,\bar{\gamma},\bar{y},\bar{z},\bar{\zeta}, \bar{u},\bar{\xi})| \\& \leq K(s,r)\left(\Vert\gamma-\bar{\gamma}\Vert_{D}+|y-\bar{y}| + \Vert z-\bar{z}\Vert + \Vert \zeta-\bar{\zeta}\Vert\right. \\& \quad\left. +\left|u(x)-\bar{u}(x)\right| +\left|\xi(x)-\bar{\xi}(x)\right| \right), \end{aligned} $$

a.e. (ω,s,r)Ω×Δc[t,T] for all \((\gamma,y,z,\zeta,u,\xi),(\bar {\gamma },\bar {y},\bar {z}, \bar {\zeta }, \bar {u},\bar {\xi })\in D([0,T],\mathbb {R}^{d})\times \mathbb {R}^{m}\times \mathbb {R}^{m\times d}\times \mathbb {R}^{m\times d}\times \mathbb {R}^{m}\times \mathbb {R}^{m}\). \(K:\Delta ^{c}\rightarrow [0,\infty)\) is a deterministic function such that for ε>0,

$$ \sup_{s\in[t,T]}\int_{s}^{T}|K(s,r)|^{2+\varepsilon}dr<\infty. $$

(H’3) \(\mathbb {E}\left [\int _{t}^{T}\int _{s}^{T} |f_{0}(s,r,\gamma _{r-})|^{2}dr\:ds\right ]<\infty \), where f0(s,r,γr)=f(s,r,γr,0,0,0,0,0). (H’4) \(\delta :\Omega \times [t,T]^{2}\times \mathbb {R}^{l}\rightarrow \mathbb {R}^{m}\) is measurable such that for any (s,r)[t,T]2

$$ \int_{\mathbb{R}^{l}} |\delta(s,r,z)|^{2}\nu(dz)<\infty. $$

Corollary 2

Let assumptions (H’1)–(H’4) hold and \((Y_{\gamma _{t}}(\cdot),Z_{\gamma _{t}}(\cdot,\cdot),U_{\gamma _{t}}(\cdot,\cdot,\cdot))\in \mathcal {H}^{2}{[t,T]}\) and \((\bar {Y}_{\gamma _{t}}(\cdot),\bar {Z}_{\gamma _{t}}(\cdot,\cdot), \bar {U}_{\gamma _{t}}(\cdot,\cdot,\cdot))\in \mathcal {H}^{2}{[t,T]}\) be the unique adapted M-solutions to BSVIE (23) with respect to (f,Φ) and \((\bar {f},\bar {\Phi })\), respectively. Then, for any \(\gamma \in D([0,T],\mathbb {R}^{d})\), the following stability estimate holds

$$ \begin{aligned} &\Vert (Y_{\gamma_{t}}(\cdot),Z_{\gamma_{t}}(\cdot,\cdot), U_{\gamma_{t}}(\cdot,\cdot,\cdot)) -(\bar{Y}_{\gamma_{t}}(\cdot),\bar{Z}_{\gamma_{t}}(\cdot,\cdot), \bar{U}_{\gamma_{t}}(\cdot,\cdot,\cdot)) \Vert^{2}_{\mathcal{H}^{2}[S,T]} \\[-2pt] &\leq C\mathbb{E}\left[\int_{S}^{T}|\Phi(X_{\gamma_{t},s-}) -\bar{\Phi}(X_{\gamma_{t},s-})|^{2}ds +\int_{S}^{T}\int_{S}^{T} \left|f\left(\vphantom{\int_{S}^{T}} s,r,X_{\gamma_{t},r-},Y_{\gamma_{t}}(r-),Z_{\gamma_{t}}(s,r),\right.\right.\right.\\[-2pt] &\quad\left. Z_{\gamma_{t}}(r,s), \int_{\mathbb{R}^{l}}U_{\gamma_{t}}(s,r,z)\delta(s,r,z)\nu(dz), \int_{\mathbb{R}^{l}}U_{\gamma_{t}}(r,s,z)\delta(r,s,z)\nu(dz)\right) \\[-2pt] &\quad -\bar{f}\left(s,r,X_{\gamma_{t},r-},Y_{\gamma_{t}}(r-),Z_{\gamma_{t}}(s,r), Z_{\gamma_{t}}(r,s), \int_{\mathbb{R}^{l}}U_{\gamma_{t}}(s,r,z)\delta(s,r,z)\nu(dz), \right.\\[-2pt] &\quad\left.\left.\left. \int_{\mathbb{R}^{l}}U_{\gamma_{t}}(r,s,z)\delta(r,s,z)\nu(dz) \right) \right|^{2} dr\:ds\right]. \end{aligned} $$
(24)

Proof

This proof follows the techniques in Yong (2008). We define for (s,r)Δc and \(z\in \mathbb {R}^{l}\)

$$\begin{aligned} \hat{\Phi}(X_{\gamma_{t},s}) =\Phi(X_{\gamma_{t},s}) -\bar{\Phi}(X_{\bar{\gamma}_{t},s}),\quad &\hat{Y}_{\gamma_{t}}(s) =Y_{\gamma_{t}}(s)-\bar{Y}_{\bar{\gamma}_{t}}(s),\quad &\\[-2pt] \hat{Z}_{\gamma_{t}}(s,r)=Z_{\gamma_{t}}(s,r) -\bar{Z}_{\bar{\gamma}_{t}}(s,r),\quad &\hat{U}_{\gamma_{t}}(s,r,z)=U_{\gamma_{t}}(s,r,z) -\bar{U}_{\bar{\gamma}_{t}}(s,r,z), \end{aligned} $$

and

$$\begin{aligned} &\hat{f}\left(s,r,X_{\gamma_{t},r-}\right)\\ &=f\left(\vphantom{\int_{\mathbb{R}^{l}}}s,r,X_{\gamma_{t},r-},Y_{\gamma_{t}}(r-), Z_{\gamma_{t}}(s,r),Z_{\gamma_{t}}(r,s),\right.\\[-2pt] &\quad\left.\int_{\mathbb{R}^{l}}U_{\gamma_{t}}(s,r,z)\delta(s,r,z)\nu(dz), \int_{\mathbb{R}^{l}}U_{\gamma_{t}}(r,s,z)\delta(r,s,z)\nu(dz) \right) \\[-2pt] &\quad -\bar{f}\left(\vphantom{\int_{\mathbb{R}^{l}}}s,r,X_{\gamma_{t},r-},Y_{\gamma_{t}}(r-), Z_{\gamma_{t}}(s,r),Z_{\gamma_{t}}(r,s),\right.\\[-2pt] &\quad\left.\int_{\mathbb{R}^{l}}U_{\gamma_{t}}(s,r,z)\delta(s,r,z)\nu(dz), \int_{\mathbb{R}^{l}}U_{\gamma_{t}}(r,s,z)\delta(r,s,z)\nu(dz) \right), \end{aligned} $$

where \(\bar {Z}^{i}_{\gamma _{t}}(s,r)= (\bar {Z}^{1}_{\gamma _{t}}(s,r),\ldots,\bar {Z}^{i}_{\gamma _{t}}(s,r), Z^{i+1}_{\gamma _{t}}(s,r),\ldots,Z^{d}_{\gamma _{t}}(s,r))\), i=1,…,d. \(\bar {f}_{y}(s,r,X_{\gamma _{t},r-})\), \(\bar {f}_{\zeta _{i}}(s,r,X_{\gamma _{t},r-})\), \(\bar {f}_{u}(s,r,X_{\gamma _{t},r-})\) and \(\bar {f}_{\xi }(s,r,X_{\gamma _{t}r-})\) are defined quite similarly.

Since \((Y_{\gamma _{t}}(\cdot),Z_{\gamma _{t}}(\cdot,\cdot), U_{\gamma _{t}}(\cdot,\cdot,\cdot)), (\bar {Y}_{\bar {\gamma }_{t}}(\cdot),\bar {Z}_{\bar {\gamma }_{t}}(\cdot,\cdot), \bar {U}_{\bar {\gamma }_{t}}(\cdot,\cdot,\cdot)) \in \mathcal {H}^{2}[t,T]\) are unique adapted M-solutions to BSVIE (1) with respect to (f,Φ) and \((\bar {f},\bar {\Phi })\), respectively, we get

$$ \begin{aligned} &\hat{Y}_{\gamma_{t}}(s) \\&=\hat{\Phi}(X_{\gamma_{t},s}) +\int_{s}^{T} \left[\hat{f}\left(s,r,X_{\gamma_{t},r-}\right)+\bar{f}_{y}(s,r,X_{\gamma_{t},r-}) \hat{Y}_{\gamma_{t}}(r-) +\sum_{i=1}^{d}\bar{f}_{z_{i}}(s,r,X_{\gamma_{t},r-}) \hat{Z}^{i}_{\gamma_{t}}(s,r)\right.\\ &\quad+\sum_{i=1}^{d}\bar{f}_{\zeta_{i}}(s,r,X_{\gamma_{t},r-}) \cdot\hat{Z}^{i}_{\gamma_{t}}(r,s) +\bar{f}_{u}(s,r,X_{\gamma_{t},r-})\cdot\int_{\mathbb{R}^{l}} U_{\gamma_{t}}(s,r,z)\delta(s,r,z)\nu(dz) \\&\quad\left. +\bar{f}_{\xi}(s,r,X_{\gamma_{t},r-})\cdot\int_{\mathbb{R}^{l}} U_{\gamma_{t}}(r,s,z)\delta(r,s,z)\nu(dz)\right]dr \\&\quad -\int_{s}^{T}\hat{Z}_{\gamma_{t}}(s,r)dW(r) -\int_{s}^{T}\int_{\mathbb{R}^{l}} \hat{U}_{\gamma_{t}}(s,r,z) \tilde{N}(dr,dz). \end{aligned} $$
(25)

Theorem 2 yields that \((\hat {Y}_{\gamma _{t}}(\cdot),\hat {Z}_{\gamma _{t}}(\cdot,\cdot), \hat {U}_{\gamma _{t}}(\cdot,\cdot,\cdot))\) is the unique adapted M-solution to BSVIE (25) and the following estimate holds

$$\begin{aligned} &\Vert(\hat{Y}_{\gamma_{t}}(\cdot),\hat{Z}_{\gamma_{t}}(\cdot,\cdot), \hat{U}_{\gamma_{t}}(\cdot,\cdot,\cdot))\Vert^{2}_{\mathcal{H}[t,T]} \\&\leq C\mathbb{E}\left[\int_{t}^{T}|\hat{\Phi}(X_{\gamma_{t},s-})|^{2}ds +\int_{t}^{T}\int_{t}^{T} \left|\hat{f}\left(s,r,X_{\gamma_{t},r-}\right)\right|^{2}dr\:ds \right], \end{aligned} $$

which provides the stability estimate. □

Corollary 3

Consider a BSVIE of the form

$$ \begin{aligned} & Y_{\gamma_{t}}(s) \\&=\Phi(X_{\gamma_{t},s}) +\int_{s}^{T}\left[ m\left(s,r,X_{\gamma_{t},r-},Y_{\gamma_{t}}(r-),Z_{\gamma_{t}}(s,r), Z_{\gamma_{t}}(r,s),\right.\right.\\ &\quad\left.\left.\int_{\mathbb{R}^{l}}U_{\gamma_{t}}(s,r,z)\delta(s,r,z)\nu(dz),\int_{\mathbb{R}^{l}}U_{\gamma_{t}}(r,s,z)\delta(r,s,z)\nu(dz)\right) +A(s,r,X_{\gamma_{t},r-})\right]dr\\ &\quad-\int_{s}^{T}Z_{\gamma_{t}}(s,r)dW(r)-\int_{s}^{T}\int_{\mathbb{R}^{l}}U_{\gamma_{t}} (s,r,z)\tilde{N}(dr,dz),\qquad s\in[t,T], \end{aligned} $$
(26)

where δsatisfies (H’4). Furthermore, let

$$m:\Omega\times\Delta^{c}[t,T]\times D([0,T],\mathbb{R}^{d})\times\mathbb{R}^{m}\times\mathbb{R}^{m\times d} \times\mathbb{R}^{m\times d}\times\mathbb{R}^{m} \times\mathbb{R}^{m}\rightarrow\mathbb{R}^{m} $$

be \(\mathcal {F}_{T}^{t}\otimes \mathcal {B}\left (\Delta ^{c}[t,T]\times D([0,T],\mathbb {R}^{d})\times \mathbb {R}^{m} \times \mathbb {R}^{m\times d} \times \mathbb {R}^{m\times d} \times \mathbb {R}^{m}\times \mathbb {R}^{m} \right)\)-measurable such that

$$ r\mapsto m(s,r,\gamma,y,z,\zeta,u,\xi) $$

is \(\mathbb {F}\)-progressive measurable.

Moreover, let m be linear in (y,z,ζ,u,ξ) such that m(s,r,γ,0,0,0,0,0)=0 a.s. for all \(\gamma \in D([0,T],\mathbb {R}^{d})\), (s,r)Δc[t,T], and

$$ |m(\omega,s,r,\gamma,y,z,\zeta,u,\xi)|\leq K(s,r) \left(|y|+\Vert z\Vert +\Vert \zeta\Vert +|u|+|\xi|\right), $$

a.e. (ω,s,r)Ω×Δc[t,T] for all \((\gamma,y,z,\zeta,u,\xi)\in D([0,T],\mathbb {R}^{d})\times \mathbb {R}^{m}\times \mathbb {R}^{m\times d}\times \mathbb {R}^{m\times d}\times \mathbb {R}^{m}\times \mathbb {R}^{m}\). \(K:\Delta ^{c}[t,T]\rightarrow [0,\infty)\) is a deterministic function such that for ε>0

$$ \sup_{s\in[t,T]}\int_{s}^{T}|K(s,r)|^{2+\varepsilon}dr<\infty. $$

\(A:\Omega \times \Delta ^{c}[t,T]\times D\left ([0,T],\mathbb {R}^{d}\right)\rightarrow \mathbb {R}^{m}\) is \(\mathcal {F}_{T}^{t}\otimes \mathcal {B}\left (\Delta ^{c}[t,T]\times D\left ([0,T],\mathbb {R}^{d}\right)\right)\)-measurable such that rA(s,r,γr) is \(\mathbb {F}\)-progressive measurable and \(\mathbb {E}\left [\int _{t}^{T}\int _{s}^{T} |A(s,r,\gamma _{r-})|^{2}dr\:ds\right ]<\infty \). Then, for all \(\gamma \in D\left ([0,T],\mathbb {R}^{d}\right)\), BSVIE (26) admits a unique adapted M-solution \((Y_{\gamma _{t}}(\cdot),Z_{\gamma _{t}}(\cdot,\cdot), U_{\gamma _{t}}(\cdot,\cdot,\cdot))\in \mathcal {H}^{2}[t,T]\) and the following estimate holds

$$ \begin{aligned} &\Vert(Y_{\gamma_{t}}(\cdot),Z_{\gamma_{t}}(\cdot,\cdot), U_{\gamma_{t}}(\cdot,\cdot,\cdot))\Vert^{2}_{\mathcal{H}^{2}[t,T]} \\&\leq C\left(\mathbb{E}\left[\int_{t}^{T}|\Phi(X_{\gamma_{t},s-})|^{2}ds \right] +\mathbb{E}\left[\int_{t}^{T}\int_{t}^{T} |A(s,r,X_{\gamma_{t},r-})|^{2}dr\:ds\right]\right). \end{aligned} $$
(27)

If \((Y_{\gamma _{t}}(\cdot),Z_{\gamma _{t}}(\cdot,\cdot), U_{\gamma _{t}}(\cdot,\cdot,\cdot)),\:(\bar {Y}_{\gamma _{t}}(\cdot),\bar {Z}_{\gamma _{t}}(\cdot,\cdot), \bar {U}_{\gamma _{t}}(\cdot,\cdot,\cdot))\in \mathcal {H}^{2}{[t,T]}\) are the unique adapted M-solution to BSVIE (23) with respect to (m,A,Φ) and \((\bar {m},\bar {A},\bar {\Phi })\), respectively, the following stability estimate holds

$$ \begin{aligned} &\Vert(Y_{\gamma_{t}}(\cdot),Z_{\gamma_{t}}(\cdot,\cdot), U_{\gamma_{t}}(\cdot,\cdot,\cdot)) -(\bar{Y}_{\gamma_{t}}(\cdot),\bar{Z}_{\gamma_{t}}(\cdot,\cdot), \bar{U}_{\gamma_{t}}(\cdot,\cdot,\cdot))\Vert^{2}_{\mathcal{H}^{2}[S,T]} \\&\leq C\left(\mathbb{E}\left[\int_{S}^{T} |\Phi(X_{\gamma_{t},s-})-\bar{\Phi}(X_{\gamma_{t},s-})|^{2}ds\right]\right.\\ &\quad\left. +\mathbb{E}\left[\int_{S}^{T}\int_{S}^{T} |A(s,r,X_{\gamma_{t},r-})-\bar{A}(s,r,X_{\gamma_{t},r-})|^{2} dr\:ds\right]\right) \end{aligned} $$
(28)

for any S[t,T).

Proof

This corollary is a consequence of the results in Theorem 2 and Corollary 2. □

4 Path-differentiability of path-dependent BSVIEs with jumps

In this section, we consider path-dependent BSVIEs with jumps of the form (23). To prove path-differentiability in the initial path of solutions to path-dependent BSVIEs, we need the following assumptions. (D1) \(\Phi \in L^{2}_{\mathcal {F}_{T}^{t}}(D([0,T],\mathbb {R}^{d}); \mathbb {R}^{m})\) is continuous and for \(i=1,\ldots,dD_{i}^{1}\Phi (\gamma ;[t,T])\in L^{2}_{\mathcal {F}_{T}^{t}}(D([0,T],\mathbb {R}^{d}); \mathbb {R}^{m})\) exist. (D2) The generator

$$ f:\Omega\times\Delta^{c}[t,T]\times D([0,T],\mathbb{R}^{d})\times\mathbb{R}^{m} \times\mathbb{R}^{m\times d} \times\mathbb{R}^{m\times d} \times\mathbb{R}^{m}\times\mathbb{R}^{m} \rightarrow \mathbb{R}^{m} $$

is \(\mathcal {F}_{T}^{t}\otimes \mathcal {B}\left (\Delta ^{c}[t,T]\times D([0,T],\mathbb {R}^{d})\times \mathbb {R}^{m} \times \mathbb {R}^{m\times d} \times \mathbb {R}^{m\times d} \times \mathbb {R}^{m}\times \mathbb {R}^{m} \right)\)-measurable such that

$$ r\mapsto f(s,r,\gamma,y,z,\zeta,u,\xi) $$

is \(\mathbb {F}\)-progressive measurable. Furthermore, (s,r)f(s,r,γ,y,z,ζ,u,ξ) is Lipschitz continuous in the sense that

$$ \begin{aligned} & |f(s,r,\gamma,y,z,\zeta,u,\xi) - f(s,r,\bar{\gamma},\bar{y},\bar{z},\bar{\zeta},\bar{u}, \bar{\xi})| \\& \leq K(s,r)\left(\Vert\gamma - \bar{\gamma}\Vert_{D}+|y-\bar{y}| +\Vert z-\bar{z}\Vert +\Vert \zeta-\bar{\zeta}\Vert +|u-\bar{u}| +|\xi-\bar{\xi}|\right), \end{aligned} $$

a.e. (ω,s,r)Ω×Δc[t,T], for all \((\gamma,y,z,\zeta,u,\xi),(\bar {\gamma },\bar {y},\bar {z}, \bar {\zeta }, \bar {u},\bar {\xi })\in D([0,T],\mathbb {R}^{d})\times \mathbb {R}^{m}\times \mathbb {R}^{m\times d}\times \mathbb {R}^{m\times d}\times \mathbb {R}^{m}\times \mathbb {R}^{m}\).

Moreover, \(K:\Delta ^{c}[t,T]\rightarrow [0,\infty)\) is a deterministic function such that for ε>0

$$ \sup_{s\in[t,T]}\int_{s}^{T}|K(s,r)|^{2+\varepsilon}dr<\infty. $$

(D3) \(\mathbb {E}\left [\int _{t}^{T}\int _{t}^{T} |f_{0}(s,r,\gamma _{r-})|^{2}dr\right ]<\infty \). (D4) For all \(\bar {\gamma }\in D([0,T],\mathbb {R}^{d})\), (s,r)Δc[t,T], the map

is of class \(C^{1}(\mathbb {R}^{d}\times \mathbb {R}^{m} \times \mathbb {R}^{m\times d} \times \mathbb {R}^{m\times d} \times \mathbb {R}^{m}\times \mathbb {R}^{m};\mathbb {R}^{m})\) and the first-order partial derivatives with respect to (x,y,z,ζ,u,ξ) are uniformly bounded. (D5) \(\delta :\Omega \times [t,T]^{2}\times \mathbb {R}^{l}\rightarrow \mathbb {R}^{m}\) is measurable such that for any (s,r)[t,T]2

$$ \int_{\mathbb{R}^{l}} |\delta(s,r,z)|^{2}\nu(dz)<\infty. $$

(D6) Let \(\gamma \mapsto U_{\gamma _{t}}(s,r,x)\) be differentiable for all \((s,r,x)\in [t,T]^{2}\times \mathbb {R}^{l}\) and assume that

$$ \int_{\mathbb{R}^{l}}|D^{1}_{i}U_{\gamma_{t}}(s,r,x;[t,T])|\nu(dx) <\infty, $$

for each (s,r)[t,T]2. Furthermore, assume that there exists a neighborhood (−ε,ε) and a process \(M_{\gamma _{t}}\) such that for each (s,r)[t,T]2\(\int _{\mathbb {R}^{l}} |M_{\gamma _{t}}(s,r,x)|\nu (dx)<\infty \) and such that for i=1,…,d

Theorem 3

Let (D1)–(D6) hold. Then, the function

$$D([0,T],\mathbb{R}^{d})\rightarrow\mathcal{H}^{2}[t,T], \qquad \gamma_{t}\mapsto (Y_{\gamma_{t}}(\cdot),Z_{\gamma_{t}}(\cdot,\cdot),U_{\gamma_{t}} (\cdot,\cdot,\cdot)) $$

is differentiable in the sense of Definition 1, where \(\mathbb {R}^{m}\) is replaced by the space \(\mathcal {H}^{2}[t,T]\). For i=1,…,d, s[t,T], the derivative is the unique adapted M-solution to a path-dependent BSVIE with jumps of the form

$$ \begin{aligned} & D^{1}_{i}Y_{\gamma_{t}}(s;[t,T]) \\&=D^{1}_{i}\Phi(X_{\gamma_{t},s};[t,T]) +\int_{s}^{T}D^{1}_{i}f\left(\vphantom{\int_{s}^{T}}s,r,X_{\gamma_{t},r-}, Y_{\gamma_{t}}(r-),Z_{\gamma_{t}}(s,r), Z_{\gamma_{t}}(r,s), \right.\\ &\quad\left. \int_{\mathbb{R}^{l}}U_{\gamma_{t}}(s,r,x)\delta(s,r,x)\nu(dx), \int_{\mathbb{R}^{l}}U_{\gamma_{t}}(r,s,x)\delta(r,s,x)\nu(dx); [t,T]\right) dr\\ &\quad+\int_{s}^{T} \left[ f_{y}\left(s,r,X_{\gamma_{t},r-},Y_{\gamma_{t}}(r-), Z_{\gamma_{t}}(s,r), Z_{\gamma_{t}}(r,s), \int_{\mathbb{R}^{l}}U_{\gamma_{t}}(s,r,x)\delta(s,r,x)\nu(dx),\right.\right. \\&\quad\left. \int_{\mathbb{R}^{l}}U_{\gamma_{t}}(r,s,x)\delta(r,s,x)\nu(dx) \right) \cdot D^{1}_{i}Y_{\gamma_{t}}(r;[t;T]) \\&\quad+\sum_{j=1}^{d} f_{z_{j}}\left(s,r,X_{\gamma_{t},r-},Y_{\gamma_{t}}(r-), Z_{\gamma_{t}}(s,r), Z_{\gamma_{t}}(r,s), \int_{\mathbb{R}^{l}}U_{\gamma_{t}}(s,r,x)\delta(s,r,x)\nu(dx),\right. \\&\quad\left. \int_{\mathbb{R}^{l}}U_{\gamma_{t}}(r,s,x)\delta(r,s,x)\nu(dx) \right)\cdot D^{1}_{i}Z^{j}_{\gamma_{t}}(s,r;[t,T]) \\&\quad+\sum_{j=1}^{d} f_{\zeta_{j}}\left(s,r,X_{\gamma_{t},r-},Y_{\gamma_{t}}(r-),Z_{\gamma_{t}}(s,r), Z_{\gamma_{t}}(r,s), \int_{\mathbb{R}^{l}}U_{\gamma_{t}}(s,r,x)\delta(s,r,x)\nu(dx),\right. \\&\quad\left. \int_{\mathbb{R}^{l}}U_{\gamma_{t}}(r,s,x)\delta(r,s,x)\nu(dx) \right) \cdot D^{1}_{i}Z^{j}_{\gamma_{t}}(r,s;[t,T]) \\&\quad+\int_{\mathbb{R}^{l}} f_{u} \left(s,r,X_{\gamma_{t},r-},Y_{\gamma_{t}}(r-),Z_{\gamma_{t}}(s,r), Z_{\gamma_{t}}(r,s), \int_{\mathbb{R}^{l}}U_{\gamma_{t}}(s,r,x)\delta(s,r,x)\nu(dx),\right. \\&\quad\left. \int_{\mathbb{R}^{l}}U_{\gamma_{t}}(r,s,x)\delta(r,s,x)\nu(dx) \right) \cdot D^{1}_{i}U_{\gamma_{t}}(s,r,x;[t,T]) \nu(dx)dr \\&\quad+\int_{\mathbb{R}^{l}} f_{\xi} \left(s,r,X_{\gamma_{t},r-},Y_{\gamma_{t}}(r-),Z_{\gamma_{t}}(s,r), Z_{\gamma_{t}}(r,s), \int_{\mathbb{R}^{l}}U_{\gamma_{t}}(s,r,x)\delta(s,r,x)\nu(dx),\right. \\&\quad\left.\left. \int_{\mathbb{R}^{l}}U_{\gamma_{t}}(r,s,x)\delta(r,s,x)\nu(dx) \right) \cdot D^{1}_{i}U_{\gamma_{t}}(r,s,x;[t,T]) \nu(dx)\right]dr \\&\quad -\int_{s}^{T}D^{1}_{i}Z_{\gamma_{t}}(s,r;[t,T])dW(r) -\int_{s}^{T}\int_{\mathbb{R}^{l}}D^{1}_{i} U_{\gamma_{t}}(s,r,x;[t,T]) \tilde{N}(dr,dx). \end{aligned} $$
(29)

Proof

To simplify the notation, we restrict ourselves to the one-dimensional case, m=d=l=1. The extension to the multidimensional case is straightforward. We define for (s,r)[t,T]2 and \(x\in \mathbb {R}\)

$$\begin{array}{@{}rcl@{}} &\mathcal{Y}^{h}(s) =\frac{1}{h}\left(Y_{\gamma_{t}^{h}}(s)-Y_{\gamma_{t}}(s)\right), &\mathcal{Z}^{h}(s,r) =\frac{1}{h}\left(Z_{\gamma_{t}^{h}}(s,r)-Z_{\gamma_{t}}(s,r)\right), \\&\mathcal{U}^{h}(s,r,x) =\frac{1}{h}\left(U_{\gamma_{t}^{h}}(s,r,x)-U_{\gamma_{t}}(s,r,x) \right), &\Delta\Phi^{h}(X_{\gamma_{t},s}) =\frac{1}{h}\left(\Phi(X_{\gamma_{t}^{h},s}) -\Phi(X_{\gamma_{t},s})\right). \end{array} $$

With the definition of \(X_{\gamma _{t},s}\), we get .

Define for any (s,r)[t,T]2, \(h\in \mathbb {R}\setminus \{0\}\), the mapping \(l^{s,r}_{h}:[0,1]\rightarrow D([0,T],\mathbb {R})\times \mathbb {R}\times \mathbb {R}\times \mathbb {R}\times \mathbb {R} \times \mathbb {R}\) by

$$\begin{aligned} & l^{s,r}_{h}(\theta) \\&= \left(X_{\gamma_{t},r-}+\theta (X_{\gamma_{t}^{h},r-}-X_{\gamma_{t},r-}), Y_{\gamma_{t}}(r-)+\theta h \mathcal{Y}^{h}(r-), Z_{\gamma_{t}}(s,r)+\theta h\mathcal{Z}^{h}(s,r), \right.\\ &\quad Z_{\gamma_{t}}(r,s)+\theta h\mathcal{Z}^{h}(r,s), \int_{\mathbb{R}}(U_{\gamma_{t}}(s,r,x)+\theta h \mathcal{U}^{h}(s,r,x))\delta(s,r,x)\nu(dx), \\&\quad\left. \int_{\mathbb{R}}(U_{\gamma_{t}}(r,s,x)+\theta h \mathcal{U}^{h}(r,s,x))\delta(r,s,x)\nu(dx)\right) \\&=\left(X_{\gamma_{t},r-}+\theta h, Y_{\gamma_{t}}(r-)+\theta h \mathcal{Y}^{h}(r-), Z_{\gamma_{t}}(s,r)+\theta h\mathcal{Z}^{h}(s,r),\right. \\&\quad Z_{\gamma_{t}}(r,s)+\theta h\mathcal{Z}^{h}(r,s), \int_{\mathbb{R}}(U_{\gamma_{t}}(s,r,x)+\theta h \mathcal{U}^{h}(s,r,x))\delta(s,r,x)\nu(dx), \\&\quad\left. \int_{\mathbb{R}}(U_{\gamma_{t}}(r,s,x)+\theta h \mathcal{U}^{h}(r,s,x))\delta(r,s,x)\nu(dx))\right). \end{aligned} $$

It follows for (s,r)[t,T]2

$$\begin{aligned} \frac{1}{h}\frac{\partial}{\partial \theta}l_{h}^{s,r}(\theta) =&\left(1,\mathcal{Y}^{h}(r-),\mathcal{Z}^{h}(s,r),\mathcal{Z}^{h}(r,s), \int_{\mathbb{R}}\mathcal{U}^{h}(s,r,x)\delta(s,r,x)\nu(dx),\right. \\&\left.\int_{\mathbb{R}}\mathcal{U}^{h}(r,s,x)\delta(r,s,x)\nu(dx) \right) \end{aligned} $$

and Corollary 3 provides that for \(h\in \mathbb {R}\setminus \{0\}(\mathcal {Y}^{h}(\cdot),\mathcal {Z}^{h}(\cdot,\cdot), \mathcal {U}^{h}(\cdot,\cdot,\cdot))\) is the unique adapted M-solution to the linearized path-dependent BSVIE with jumps

$$\begin{aligned} &\mathcal{Y}^{h}(s)\\ &= \Delta\Phi^{h}(X_{\gamma_{t},s})+\int_{s}^{T}\left[ m_{h}^{s,r}\left(\mathcal{Y}^{h}(r-),\mathcal{Z}^{h}(s,r), \mathcal{Z}^{h}(r,s), \int_{\mathbb{R}}\mathcal{U}^{h}(s,r,x)\delta(s,r,x)\nu(dx),\right.\right.\\ &\quad\left.\left.\int_{\mathbb{R}}\mathcal{U}^{h}(r,s,x)\delta(r,s,x)\nu(dx) \right) +A^{1}_{h}(s,r) \right]dr\\ &\quad -\int_{s}^{T}\mathcal{Z}^{h}(s,r)dW(r) -\int_{s}^{T}\int_{\mathbb{R}}\mathcal{U}^{h}(s,r,x) \tilde{N}(dr,dx),\qquad s\in[t,T], \end{aligned} $$

where \(m_{h}^{s,r}:\mathbb {R}\times \mathbb {R}\times \mathbb {R} \times \mathbb {R}\times \mathbb {R} \rightarrow \mathbb {R}\) is a linear function defined by

$$\begin{aligned} m_{h}^{s,r}(y,z,\zeta,u,\mu)=A^{2}_{h}(s,r)y +A^{3}_{h}(s,r)z +A^{4}_{h}(s,r)\zeta +A^{5}_{h}(s,r)u+ A^{6}_{h}(s,r)\mu, \end{aligned} $$

for

$$\begin{aligned} &A^{1}_{h}(s,r)=\int_{0}^{1}D^{1}f(l_{h}^{s,r}(\theta);[t,T]) d\theta, \qquad A^{2}_{h}(s,r) =\int_{0}^{1}\frac{\partial}{\partial y}f(l_{h}^{s,r}(\theta))d\theta,&\\& A^{3}_{h}(s,r) =\int_{0}^{1}\frac{\partial}{\partial z} f(l_{h}^{s,r}(\theta))d\theta, \qquad\qquad A^{4}_{h}(s,r) =\int_{0}^{1}\frac{\partial}{\partial \zeta}f(l_{h}^{s,r}(\theta))d\theta,&\\& A^{5}_{h}(s,r) =\int_{0}^{1}\frac{\partial}{\partial u} f(l_{h}^{s,r}(\theta))d\theta, \qquad\qquad A^{6}_{h}(s,r) =\int_{0}^{1}\frac{\partial}{\partial \mu} f(l_{h}^{s,r}(\theta))d\theta.& \end{aligned} $$

Additionally, we get with the inequality (28)

$$\begin{aligned} \mathbb{E}\left[\sup_{s\in[t,T]} |\mathcal{Y}^{h}(s)-\mathcal{Y}^{\bar{h}}(s)|^{2}\right] &\leq C\left(\mathbb{E}\left[\int_{t}^{T}|\Delta\Phi^{h}(X_{\gamma_{t},s-}) -\Delta\Phi^{\bar{h}}(X_{\gamma_{t},s-})|^{2}ds\right]\right.\\ &\left.+\mathbb{E}\left[\int_{t}^{T} \int_{t}^{T} |A^{1}_{h}(s,r)-A^{1}_{\bar{h}}(s,r)|^{2}dr\:ds\right]\right). \end{aligned} $$

We obtain for \(h,\bar {h}\in \mathbb {R}\setminus \{0\}\) with (D1)

$$\begin{aligned} &{\lim}_{h\rightarrow\bar{h}}\mathbb{E}\left[\int_{t}^{T} |\Delta\Phi^{h}(X_{\gamma_{t},s-}) -\Delta\Phi^{\bar{h}}(X_{\gamma_{t},s-})|^{2}ds\right] \\&\leq {\lim}_{h\rightarrow\bar{h}}\mathbb{E}\left[\int_{t}^{T} \left|\frac{1}{\bar{h}}\Phi(X_{\gamma_{t}^{\bar{h}},s-}) -\frac{1}{h}\Phi(X_{\gamma_{t}^{h},s-}) \right|^{2}ds\right]=0. \end{aligned} $$

(D4) provides

$$\begin{aligned} &{\lim}_{h\rightarrow\bar{h}} \mathbb{E}\left[\int_{t}^{T}\int_{t}^{T}\left| A^{1}_{h}(s,r)-A^{1}_{\bar{h}}(s,r)\right|^{2}dr\:ds\right] \\&\leq{\lim}_{h\rightarrow\bar{h}}\mathbb{E}\left[\int_{t}^{T} \int_{t}^{T} \int_{0}^{1}\left| D^{1}f(l_{h}^{s,r}(\theta);[t,T]) -D^{1}f(l_{\bar{h}}^{s,r}(\theta);[t,T])\right|^{2} d\theta \: dr\:ds\right]=0. \end{aligned} $$

Thus, it follows

$$\begin{aligned} {\lim}_{h\rightarrow\bar{h}}\mathbb{E}\left[\sup_{s\in[t,T]} |\mathcal{Y}^{h}(s)-\mathcal{Y}^{\bar{h}}(s)|^{2}\right]=0. \end{aligned} $$

With the same arguments, we get

$$\begin{aligned} {\lim}_{h\rightarrow\bar{h}}\mathbb{E}\left[\int_{t}^{T} \int_{t}^{T} |\mathcal{Z}^{h}(s,r)-\mathcal{Z}^{\bar{h}}(s,r)|^{2}dr\:ds\right]=0 \end{aligned} $$

and

$$\begin{aligned} {\lim}_{h\rightarrow\bar{h}}\mathbb{E}\left[\int_{t}^{T}\int_{t}^{T}\int_{\mathbb{R}} |\mathcal{U}^{h}(s,r,x) -\mathcal{U}^{\bar{h}}(s,r,x)|^{2}\nu(dx)dr\:ds\right]=0. \end{aligned} $$

Finally, this provides

$$\begin{aligned} {\lim}_{h\rightarrow\bar{h}}\Vert (\mathcal{Y}^{h},\mathcal{Z}^{h},\mathcal{U}^{h}) -(\mathcal{Y}^{\bar{h}},\mathcal{Z}^{\bar{h}}, \mathcal{U}^{\bar{h}}) \Vert^{2}_{\mathcal{H}^{2}[t,T]} =0. \end{aligned} $$

Since S2[t,T], L2(t,T;M2[t,T]) and \(L^{2}(t,T;L^{2}_{N}(t,T))\) are Banach spaces, the sequence \(\mathcal {Y}^{h_{n}}(s)\) converges to a process \(D^{1}Y_{\gamma _{t}}(s;[t,T])\), the sequence \(\mathcal {Z}^{h_{n}}(s,r)\) converges to a process \(D^{1}Z_{\gamma _{t}}(s,r;[t,T])\), and the process \(\mathcal {U}^{h_{n}}(s,r,x)\) converges to \(D^{1}U_{\gamma _{t}}(s,r,x;[t,T])\) for a sequence \((h_{n})_{n}\in \mathbb {R}\setminus \{0\}\) and with respect to the corresponding norms.

This implies that \((D^{1}Y_{\gamma _{t}}(\cdot ;[t,T]), D^{1}Z_{\gamma _{t}}(\cdot,\cdot ;[t,T]), D^{1}U_{\gamma _{t}}(\cdot,\cdot,\cdot ;[t,T]))\in \mathcal {H}^{2}[t,T]\) is the adapted M-solution to the path-dependent BSVIE with jumps (29). □

5 Duality Principle and Comparison Theorem

In this section, we prove the duality principle for linear path-dependent BSVIEs with jumps

$$ \begin{aligned} Y_{\gamma_{t}}(s) =&\Phi(X_{\gamma_{t},s}) +\int_{s}^{T}\left[ A^{T}(s,r,X_{\gamma_{t},r-})Y_{\gamma_{t}}(r-)+\sum_{j=1}^{d}B^{T}_{j}(s,X_{\gamma_{t},s-})Z^{j}_{\gamma_{t}}(r,s)\right. \\&\left.+\int_{\mathbb{R}^{l}} C^{T}(s,X_{\gamma_{t},s-},z)U_{\gamma_{t}}(r,s,z)\nu(dz)\right]dr\\ &-\int_{s}^{T}Z_{\gamma_{t}}(s,r)dW(r)-\int_{s}^{T}\int_{\mathbb{R}^{l}}U_{\gamma_{t}}(s,r,z)\tilde{N}(dr,dz),\qquad s\in[t,T], \end{aligned} $$
(30)

and linear path-dependent FSVIEs with jumps of the form

$$ \begin{aligned} X_{\gamma_{t}}(s)=&\varphi(X_{\gamma_{t},s}) +\int_{t}^{s}A(r,s,X_{\gamma_{t},s-})X_{\gamma_{t}}(r-)dr +\sum_{j=1}^{d}\int_{t}^{s}B_{j}(r, X_{\gamma_{t},r-}) X_{\gamma_{t}}(r-){dW}_{j}(r) \\&+\int_{t}^{s}\int_{\mathbb{R}^{l}} C(r,X_{\gamma_{t},r-},z)X_{\gamma_{t}}(r-) \tilde{N}(dr,dz), \qquad s\in[t,T], \\ X_{\gamma_{t}}(s)=&\gamma_{t}(s),\qquad s\in[0,t), \end{aligned} $$
(31)

where it is necessary that m=d for the dimensions. First, we show that such a FSVIE admits an adapted solution. We consider a more general path-dependent FSVIE with jumps of the form

$$ \begin{aligned} X_{\gamma_{t}}(s)=&\varphi(X_{\gamma_{t},s}) +\int_{t}^{s}\bar{A}(r,s,X_{\gamma_{t},s-})dr +\sum_{j=1}^{d}\int_{t}^{s}\bar{B}_{j}(r, X_{\gamma_{t},r-}){dW}_{j}(r) \\&+\int_{t}^{s}\int_{\mathbb{R}^{l}} \bar{C}(r,X_{\gamma_{t},r-},z) \tilde{N}(dr,dz), \qquad s\in[t,T], \\ X_{\gamma_{t}}(s)=&\gamma_{t}(s),\qquad s\in[0,t), \end{aligned} $$
(32)

and assume the following. (C1) Let \(\varphi (\cdot)\in L^{2}_{\mathbb {F}}(\Omega \times D([0,T],\mathbb {R}^{d});\mathbb {R}^{d})\) with

$$ |\varphi(\gamma)|^{2}\leq K\left(1+\sup_{u\in[0,T]}|\gamma(u)|^{2}\right) $$

and

$$ |\varphi(\gamma_{s})-\varphi(\bar{\gamma}_{\bar{s}})|^{2} \leq K\left(\sup_{u\in[0,T]}|\gamma(u)-\bar{\gamma}(u)|^{2} +|s-\bar{s}|^{2}\right), $$

for any \(\gamma,\bar {\gamma }\in D([0,T],\mathbb {R}^{d})\) and \(s,\bar {s}\in [t,T]\). (C2) The process \(\bar {A}:\Omega \times \Delta ^{c}[t,T]\times D([0,T],\mathbb {R}^{d})\rightarrow \mathbb {R}^{d}\) is \(\mathcal {F}^{t}_{T}\otimes \mathcal {B}(\Delta ^{c}[t,T]\times D([0,T],\mathbb {R}^{d}))\)-measurable such that \(r\mapsto \bar {A}(r,s,\gamma)\) is progressively measurable on [t,s] for any s[t,T], \(\gamma \in D([0,T],\mathbb {R}^{d})\). Moreover, map \(\bar {A}\) is Lipschitz continuous in the sense that

$$ \left| \bar{A}(r,s,\gamma) -\bar{A}(r,\bar{s},\bar{\gamma})\right|^{2} \leq K\left(\sup_{\tau\in[t,r]}|\gamma(\tau) -\bar{\gamma}(\tau)|^{2} +|s-\bar{s}|^{2}\right), $$

a.e. \((\omega,r)\in \Omega \times [t,s\wedge \bar {s}]\) for all \(s,\bar {s}\in [t,T]\), \(\gamma,\bar {\gamma }\in D([0,T],\mathbb {R}^{d})\), and satisfies the following growth condition

$$ \left| \bar{A}\left(r,s,\gamma\right)\right|^{2} \leq K\left(1+\sup_{u\in[0,s]}|\gamma(u)|^{2}\right). $$

(C3) The map \((r,x,z)\mapsto (\bar {B}_{1}(r,x),\ldots,\bar {B}_{d}(r,x),\bar {C}(r,x,z))\) defined on \([t,T]\times D([0,T],\mathbb {R}^{d})\times \mathbb {R}^{l}\) is Borel measurable and the following Lipschitz and growth conditions hold

$$ \sum_{j=1}^{d}| \bar{B}_{j}(s,\gamma) - \bar{B}_{j}(s,\bar{\gamma})|^{2} +\int_{\mathbb{R}^{l}} |\bar{C}(s,\gamma,z)-\bar{C}(s,\bar{\gamma},z)|^{2} \nu(dz) \leq K \sup_{u\in[0,s]}|\gamma(u)-\bar{\gamma}(u)|^{2} $$

and

$$ \sum_{j=1}^{d}| \bar{B}_{j}(s,\gamma)|^{2} +\int_{\mathbb{R}^{l}} |\bar{C}(s,\gamma,z)|^{2}\nu(dz) \leq K\left(1+\sup_{u\in[0,s]}|\gamma(u)|^{2}\right), $$

for each s[t,T] and for all \(\gamma,\bar {\gamma }\in D([0,T],\mathbb {R}^{d})\).

Theorem 4

Let (C1)–(C3) hold. Then, for any \(\gamma \in D\left ([0,T],\mathbb {R}^{d}\right)\) with \(\sup _{s\in [0,t]}|\gamma (s)|^{2}<\infty \), t[0,T], there exists a strong solution to FSVIE (32).

Proof

The proof is quite similar to Kromer et al. (2015). □

Now, let m=d and we consider the linear cases of path-dependent FSVIEs and BSVIEs with jumps (31) and (30), where we regard

$$\begin{array}{@{}rcl@{}} \bar{A}(r,s,\gamma)&=A(r,s,\gamma)\gamma(r-), \end{array} $$
(33)
$$\begin{array}{@{}rcl@{}} \bar{B}_{j}(r,\gamma)&=B_{j}(r,\gamma)\gamma(r-),& j=1,\cdots,d, \end{array} $$
(34)
$$\begin{array}{@{}rcl@{}} \bar{C}(r,\gamma,z)&=C(r,\gamma,z)\gamma(r-), \end{array} $$
(35)

for any \((r,s,\gamma)\in [t,T]^{2}\times \mathbb {R}^{l}\) and with \(A(\cdot,\cdot,\cdot),B_{j}(\cdot,\cdot),C(\cdot,\cdot) \in \mathbb {R}^{d\times d}\).

Theorem 5

Let \(\Phi \in L^{2}_{\mathcal {F}_{T}^{t}} (D([0,T],\mathbb {R}^{d});\mathbb {R})\), \(\gamma \in D([0,T],\mathbb {R}^{d})\), with \(\sup _{s\in [0,t]}|\gamma (s)|^{2}<\infty \), t[0,T] and \(\varphi (\cdot)\in L^{2}_{\mathbb {F}}(\Omega \times D([0,T],\mathbb {R}^{d});\mathbb {R}^{d})\) satisfying (C1).

Furthermore, \(A:\Omega \times [t,T]^{2}\times D([0,T],\mathbb {R}^{d})\rightarrow \mathbb {R}^{d\times d}\), \(B_{j}:\Omega \times [t,T]\times D([0,T],\mathbb {R}^{d})\rightarrow \mathbb {R}^{d\times d}\), \(C:\Omega \times [t,T]\times D([0,T],\mathbb {R}^{d})\times \mathbb {R}^{l}\rightarrow \mathbb {R}^{d\times d}\) are uniformly bounded and (C2)–(C3) hold for (33)–(35).

Let \(X_{\gamma _{t}}\) be the strong solution to FSVIE (31) and \((Y_{\gamma _{t}}(\cdot),Z_{\gamma _{t}}(\cdot,\cdot), U_{\gamma _{t}}(\cdot,\cdot,\cdot))\in \mathcal {H}^{2}[t,T]\) be the unique adapted M-solution to BSVIE (30). Then, we have the following duality principle

$$ \begin{aligned} \mathbb{E}\left[\int_{t}^{T} \langle \varphi(X_{\gamma_{t},r-}), Y_{\gamma_{t}}(r-)\rangle dr\right]= \mathbb{E}\left[\int_{t}^{T} \langle X_{\gamma_{t}}(r-), \Phi(X_{\gamma_{t},r-})\rangle dr\right]. \end{aligned} $$
(36)

Proof

Since \((Y_{\gamma _{t}}(\cdot),Z_{\gamma _{t}}(\cdot,\cdot), U_{\gamma _{t}}(\cdot,\cdot,\cdot))\in \mathcal {H}^{2}[t,T]\) is the unique adapted M-solution to BSVIE (30), we get

$$Y_{\gamma_{t}}(s)=\mathbb{E}\left[Y_{\gamma_{t}}(s)\right] +\int_{t}^{s}Z_{\gamma_{t}}(s,r)dW(r)+\int_{t}^{s}\int_{\mathbb{R}^{l}}U_{\gamma_{t}}(s,r,z) \tilde{N}(dr,dz). $$

With Fubini’s theorem, this implies

$$\begin{aligned} &\mathbb{E}\left[\int_{t}^{T} \langle \varphi(X_{\gamma_{t},r-}), Y_{\gamma_{t}}(r-)\rangle dr\right] \\&=\mathbb{E}\left[\int_{t}^{T} \langle X_{\gamma_{t}}(s-), Y_{\gamma_{t}}(s-)\rangle ds\right] -\mathbb{E}\left[\int_{t}^{T} \langle \int_{t}^{s}A(r,s,X_{\gamma_{t},s-})X_{\gamma_{t}}(r-)dr, Y_{\gamma_{t}}(s-)\rangle ds\right] \\&\quad-\mathbb{E}\left[\int_{t}^{T} \langle \int_{t}^{s}\sum_{j=1}^{d}B_{j}(r,s,X_{\gamma_{t},s-}) X_{\gamma_{t}}(r-){dW}_{j}(r)\right.\\ &\quad+\int_{t}^{s} \int_{\mathbb{R}^{l}} C(r,s,X_{\gamma_{t},s-},z)X_{\gamma_{t}}(r-)\tilde{N}(dr,dz), \\&\quad\left. \mathbb{E}\left[Y_{\gamma_{t}}(s-)\right] +\int_{t}^{s}Z_{\gamma_{t}}(s,r)dW(r)+\int_{t}^{s}\int_{\mathbb{R}^{l}}U_{\gamma_{t}}(s,r,z) \tilde{N}(dr,dz)\rangle ds\right] \\&=\mathbb{E}\left[\int_{t}^{T}\langle X_{\gamma_{t}}(s-), Y_{\gamma_{t}}(s-)-\int_{s}^{T} \left(A^{T}(s,r,X_{\gamma_{t},r-})Y_{\gamma_{t}}(r-)+\sum_{j=1}^{d}B^{T}_{j}(s,r,X_{\gamma_{t},r-})Z^{j}_{\gamma_{t}}(r,s)\right.\right.\\ &\quad\left.\left.+\int_{\mathbb{R}^{l}}C^{T}(s,r,X_{\gamma_{t},r-},z) U_{\gamma_{t}}(r,s,z)\rangle \nu(dz)\right)dr\rangle ds \right]\\ &=\mathbb{E}\left[\int_{t}^{T}\langle X_{\gamma_{t}}(s-), \Phi(X_{\gamma_{t},s-})\rangle ds\right]. \end{aligned} $$

Next, we consider a one-dimensional linear path-dependent FSVIE with jumps of the form

$$ \begin{aligned} X_{\gamma_{t}}(s) &=\varphi(X_{\gamma_{t},s}) +\int_{t}^{s}a(r,s,X_{\gamma_{t},s-})X_{\gamma_{t}}(r-)dr +\int_{t}^{s} b(r)X_{\gamma_{t}}(r-)dW(r) \\&\quad +\int_{t}^{s}\int_{\mathbb{R}^{l}} c(r,z)X_{\gamma_{t}}(r-)\tilde{N}(dr,dz), \quad s\in[t,T], \\ X_{\gamma_{t}}(s)&=\gamma_{t}(s),\quad s\in[0,t), \end{aligned} $$
(37)

where the coefficients b and c do not depend on the path of process \(X_{\gamma _{t}}\). For the proof of the comparison theorem to BSVIEs we need a comparison theorem for FSVIEs.

Theorem 6

Let \(a:\Omega \times [t,T]^{2}\times D([0,T],\mathbb {R})\rightarrow \mathbb {R}\) be \(\mathcal {F}^{t}_{T}\otimes \mathcal {B}([t,T]^{2} \times D([0,T],\mathbb {R}))\)-measurable such that for any s[t,T], \(\gamma \in D([0,T],\mathbb {R})\) the map \(r\rightarrow a(r,s,\gamma _{s-})\) is \(\mathbb {F}\)-progressively measurable on [t,s]. Moreover, a satisfies the following Lipschitz and growth conditions

$$ \begin{aligned} |a(r,s,\gamma_{s-})-a(r,\bar{s},\bar{\gamma}_{\bar{s}-})|^{2} \leq K\left(\sup_{s\in[t,r]}|\gamma(s)-\bar{\gamma}(s)|^{2}+|s-\bar{s}|^{2}\right), \end{aligned} $$
(38)

a.e. \((\omega,r)\in \Omega \times [t,s\wedge \bar {s}]\) for all \(s,\bar {s}\in [t,T]\), \(\gamma,\bar {\gamma }\in D([0,T],\mathbb {R})\), and

$$ \begin{aligned} |a(r,s,\gamma_{s-})|^{2}\leq K\left(1+\sup_{s\in[t,r]}|\gamma(s)|^{2}\right), \end{aligned} $$
(39)

a.e. (ω,r)Ω×[t,s] for all s[t,T], \(\gamma \in D([0,T],\mathbb {R})\). Furthermore, let

$$ \begin{aligned} a(r,\bar{s},\gamma_{\bar{s}-})- a(r,s,\gamma_{s-})\geq0, \end{aligned} $$
(40)

a.e. (ω,r)Ω×[t,s] for any \(\gamma \in D([0,T],\mathbb {R})\), \(s,\bar {s}\in [t,T]\) with \(\bar {s}\geq s\). Let \((r,z)\rightarrow \left (b(r),c(r,z)\right)\) be Borel measurable, where b is uniformly bounded and

$$\int_{\mathbb{R}^{l}}|c(r,z)|^{2}\nu(dz)<\infty, $$

for any r[t,T]. In addition, c(r,z)>−1 a.s. for all \((r,z)\in [t,T]\times \mathbb {R}^{l}\). Then, for any \(\varphi (\cdot)\in L^{2}_{\mathbb {F}}(\Omega \times D([0,T],\mathbb {R});\mathbb {R})\), satisfying (C1), with

$$ \varphi(\gamma_{\bar{s}})\geq \varphi(\gamma_{s})\geq 0, \quad s,\bar{s}\in[t,T],\:\bar{s}\geq s,\quad \text{a.s.}, $$
(41)

FSVIE (37) admits a unique solution \(X_{\gamma _{t}}(\cdot)\). Moreover, it is

$$ X_{\gamma_{t}}(s)\geq 0,\quad\:s \in[t,T],\:\text{a.s.} $$
(42)

Proof

We will use the techniques from Wang and Yong (2015) and Lu (2016). Let Π={t k ,0≤kN} be a finite sequence of real numbers with t=t0<t1<…<t N =T, \(\Vert \Pi \Vert =\max _{1\leq k\leq N}|t_{k}-t_{k-1}|\).

Define

where \(\varphi (X_{\gamma _{t},t_{k}})\leq \varphi (X_{\gamma _{t},t_{k+1}})\) a.s.

For any trsT, \(\gamma \in D([0,T],\mathbb {R}^{d})\), we note that with (C2)

$$\begin{aligned} &\left|a^{\Pi}(r,s,X_{\gamma_{t},s-}) - a(r,s,X_{\gamma_{t},s-})\right| \leq K\left(\sup_{\tau\in[s,\bar{s}],|s-\bar{s}|\leq \Vert\Pi\Vert}\left|X_{\gamma_{t},s}(\tau)-X_{\gamma_{t},\bar{s}} (\tau)\right| +\Vert\Pi\Vert\right). \end{aligned} $$

Furthermore, we obtain

$$\begin{aligned} &\mathbb{E}\left[\sup_{s\in[t,T]}|\varphi^{\Pi}(X_{\gamma_{t},s}) -\varphi(X_{\gamma_{t},s})|^{2}\right] \leq \mathbb{E}\left[ \sup_{|s-\bar{s}|\leq\Vert\Pi\Vert,\:s,\bar{s}\in[t,T]}\left| \varphi(X_{\gamma_{t},s})-\varphi(X_{\gamma_{t},\bar{s}}) \right|^{2}\right]. \end{aligned} $$

Now, let XΠ(·) be the solution to the following FSVIE with jumps

$$ \begin{aligned} X^{\Pi}_{\gamma_{t}}(s) &=\varphi^{\Pi}(X_{\gamma_{t},s}) +\int_{t}^{s}a^{\Pi}(r,s,X_{\gamma_{t},s-})X^{\Pi}_{\gamma_{t}}(r-)dr +\int_{t}^{s}b(r)X^{\Pi}_{\gamma_{t}}(r-)dW(r) \\&\quad+\int_{t}^{s}\int_{\mathbb{R}^{l}} c(r,z)X^{\Pi}_{\gamma_{t}}(r-)\tilde{N}(dr,dz), \; s\in[t,T], \\ X^{\Pi}_{\gamma_{t}}(s)&=\gamma(s),\quad s\in[0,t). \end{aligned} $$
(43)

We get

$$\begin{aligned} &\mathbb{E}\left[\sup_{s\in[t,T]}\left|X^{\Pi}_{\gamma_{t}}(s)-X_{\gamma_{t}}(s) \right|^{2}\right] \\&\leq C\mathbb{E}\left[\sup_{s\in[t,T]}|\varphi^{\Pi}(X_{\gamma_{t},s})-\varphi(X_{\gamma_{t},s})|^{2}+\sup_{s\in[t,T]}\left|\int_{t}^{s}b(r) (X^{\Pi}_{\gamma_{t}}(r-)-X_{\gamma_{t}}(r-))dW(r)\right|^{2}\right.\\ &\quad +\sup_{s\in[t,T]}\left(\int_{t}^{s}\left|a^{\Pi}(r,s,X_{\gamma_{t},s-}) X^{\Pi}_{\gamma_{t}}(r-)-a(r,s,X_{\gamma_{t},s-})X_{\gamma_{t}}(r-) \right|dr\right)^{2} \\&\quad\left.+\sup_{s\in[t,T]}\left|\int_{t}^{s}\int_{\mathbb{R}^{l}}c(r,z) (X^{\Pi}_{\gamma_{t}}(r-)-X_{\gamma_{t}}(r-))\tilde{N}(dr,dz)\right|^{2} \right]. \end{aligned} $$

It is

$$\begin{aligned} &\left|a^{\Pi}(r,s,X_{\gamma_{t},s-}) X^{\Pi}_{\gamma_{t}}(r-) -a(r,s,X_{\gamma_{t},s-})X_{\gamma_{t}}(r-)\right| \\&\leq \left|a^{\Pi}(r,s,X_{\gamma_{t},s-})\right| \left|X^{\Pi}_{\gamma_{t}}(r-)-X_{\gamma_{t}}(r-)\right| +\left|a^{\Pi}(r,s,X_{\gamma_{t},s-})-a(r,s,X_{\gamma_{t},s-})\right| \left|X_{\gamma_{t}}(r-)\right| \end{aligned} $$

and

$$\begin{aligned} &\mathbb{E}\left[\sup_{s\in[t,T]}\left(\int_{t}^{s} \left|a^{\Pi}(r,s,X_{\gamma_{t},s-})\right| \left|X^{\Pi}_{\gamma_{t}}(r-)-X_{\gamma_{t}}(r-)\right|dr\right)^{2}\right] \\&\leq \mathbb{E}\left[\sup_{s\in[t,T]}\int_{t}^{s} K\cdot\left(1+\sup_{\tau\in[t,r]}|X_{\gamma_{t}}(\tau)|^{2}\right) \left|X^{\Pi}_{\gamma_{t}}(r-)-X_{\gamma_{t}}(r-)\right|^{2}dr\right] \\& \leq K \mathbb{E}\left[\sup_{\tau\in[s,\bar{s}],s,\bar{s}\leq\Vert\Pi\Vert}\left|X^{\Pi}_{\gamma_{t}}(r-)-X_{\gamma_{t}}(r-)\right|^{2} \int_{t}^{T}\sup_{\tau\in[t,r]}|X_{\gamma_{t}}(\tau)|^{2}dr\right.\\ &\quad\left. +\int_{t}^{T}\left|X^{\Pi}_{\gamma_{t}}(r-)-X_{\gamma_{t}}(r-)\right|^{2}dr\right]. \end{aligned} $$

Since \(b:\omega \times [t,T]\rightarrow \mathbb {R}\), \(c:\omega \times [t,T]\times \mathbb {R}^{l}\rightarrow \mathbb {R}\) are Borel measurable and \(X_{\gamma _{t}}^{\Pi },X_{\gamma _{t}}\) are both càdlàg and adapted processes, which implies that \(X_{\gamma _{t}}^{\Pi }\) and \(X_{\gamma _{t}}\) are progressively measurable, the integrals

$$\int_{t}^{\cdot} b(r)(X_{\gamma_{t}}^{\Pi}(r-)-X_{\gamma_{t}}(r-))dW(r) \text{ and } \int_{t}^{\cdot}\int_{\mathbb{R}^{l}}c(r,z)(X_{\gamma_{t}}^{\Pi}(r-) -X_{\gamma_{t}}(r-))\tilde{N}(dr,dz) $$

are local martingales. Since \(b:\omega \times [t,T]\rightarrow \mathbb {R}\) is bounded, we get with Burkholder–Davis– Gundy inequality

$$\begin{aligned} &\mathbb{E}\left[\sup_{s\in[t,T]} \left|\int_{t}^{s}b(r)\left(X^{\Pi}_{\gamma_{t}}(r-) -X_{\gamma_{t}}(r-)\right)dW(r)\right|^{2}\right] \\&\leq C \mathbb{E}\left[\int_{t}^{T} \left|X^{\Pi}_{\gamma_{t}}(r-) -X_{\gamma_{t}}(r-)\right|^{2}dr\right]. \end{aligned} $$

Similarly, it is

$$\begin{aligned} &\mathbb{E}\left[\left| \int_{t}^{s}\int_{\mathbb{R}^{l}} c(r,z)\left(X^{\Pi}_{\gamma_{t}}(r-) -X_{\gamma_{t}}(r-)\right) \tilde{N}(dr,dz)\right|^{2}\right] \\&\leq C\mathbb{E}\left[\int_{t}^{T} \left|X^{\Pi}_{\gamma_{t}}(r-) -X_{\gamma_{t}}(r-)\right|^{2}dr\right] \end{aligned} $$

and with Gronwall’s lemma we get for any s[t,T]

$$\begin{aligned} &\mathbb{E}\left[\left|X_{\gamma_{t}}^{\Pi}(s) -X_{\gamma_{t}}(s)\right|^{2}\right] \\&\leq C\left\{ \mathbb{E}\left[ \sup_{s,\bar{s}\in[t,T],|s-\bar{s}|\leq\Vert\Pi\Vert}\left| \varphi(X_{\gamma_{t},s})-\varphi(X_{\gamma_{t},\bar{s}}) \right|^{2}\right]\right.\\ &\quad \left.+\mathbb{E}\left[\left(\sup_{\tau\in[s,\bar{s}],|s-\bar{s}|\leq \Vert\Pi\Vert}\left|X_{\gamma_{t},s-}(\tau)-X_{\gamma_{t},\bar{s}-} (\tau)\right|^{2} +\Vert\Pi\Vert^{2}\right) \int_{t}^{T} \left|X_{\gamma_{t}}(r-)\right|^{2}dr\right]\right\} e^{C(T-t)}. \end{aligned} $$

Hence, it follows for any s[t,T]

$$ {\lim}_{\Vert\Pi\Vert\rightarrow 0}\mathbb{E}\left[\left|X_{\gamma_{t}}^{\Pi}(s) -X_{\gamma_{t}}(s)\right|^{2}\right]=0. $$
(44)

For s[t,T1) we get with definition (43) a linear FSDE of the form

$$\begin{aligned} X^{\Pi}_{\gamma_{t}}(s) &=\varphi(X_{\gamma_{t},t}) +\int_{t}^{s}a(r,t,X_{\gamma_{t},t})X^{\Pi}_{\gamma_{t}}(r-) dr +\int_{t}^{s}b(r)X_{\gamma_{t}}^{\Pi}(r-)dW(r) \\&\quad+\int_{t}^{s}\int_{\mathbb{R}}c(r,z)X^{\Pi}_{\gamma_{t}}(r-) \tilde{N}(dr,dz), \end{aligned} $$

where the solution \(X^{\Pi }_{\gamma _{t}}(s)\), s[t,T1), is a Doléans–Dade exponential of the semimartingale \(X_{\gamma _{t}}\) which is positive for c(r,z)>−1 for all \((s,z)\in [t,t_{1})\times \mathbb {R}^{l}\), see Doléans-Dade (1970). In particular, it is

$$\begin{aligned} X^{\Pi}_{\gamma_{t}}(t_{1}-) &=\varphi(X_{\gamma_{t},t}) +\int_{t}^{t_{1}}a(r,t,X_{\gamma_{t},t})X^{\Pi}_{\gamma_{t}}(r-) dr +\int_{t}^{t_{1}}b(r)X_{\gamma_{t}}^{\Pi}(r-)dW(r) \\&\quad+\int_{t}^{t_{1}}\int_{\mathbb{R}}c(r,z)X^{\Pi}_{\gamma_{t}}(r-) \tilde{N}(dr,dz). \end{aligned} $$

Now, for s[t1,t2) it is

$$\begin{aligned} X^{\Pi}_{\gamma_{t}}(s) &=\hat{\varphi}(X_{\gamma_{t},t_{1}}) +\int_{t_{1}}^{s}a(r,t_{1},X_{\gamma_{t},t_{1}}) X^{\Pi}_{\gamma_{t}}(r-)dr +\int_{t_{1}}^{s}b(r)X^{\Pi}_{\gamma_{t}}(r-)dW(r) \\&\quad+\int_{t_{1}}^{s}\int_{\mathbb{R}}c(r,z)X^{\Pi}_{\gamma_{t}}(r-) \tilde{N}(dr,dz), \end{aligned} $$

where

$$\hat{\varphi}(X_{\gamma_{t},t_{1}}) =\varphi(X_{\gamma_{t},t_{1}})-\varphi(X_{\gamma_{t},t}) +X^{\Pi}_{\gamma_{t}}(t_{1}-) +\int_{t}^{t_{1}}\left(a(r,t_{1},X_{\gamma_{t},t_{1}}) -a(r,t,X_{\gamma_{t},t})\right)X^{\Pi}_{\gamma_{t}}(r-)dr. $$

Since (40) holds, \(\varphi (X_{\gamma _{t},t_{1}})\geq \varphi (X_{\gamma _{t},t})\) for tt1 and \(X_{\gamma _{t}}^{\Pi }(t_{1}-)\geq 0\), it is \(\hat {\varphi }(X_{\gamma _{t},t_{1}})\geq 0\). With the same arguments for s[t,T1), it is \(X^{\Pi }_{\gamma _{t}}(s)\geq 0\) a.s. for any s[t1,t2). By induction we get that \(X^{\Pi }_{\gamma _{t}}(s)\geq 0\) a.s. for all s[t,T]. Equation (44) yields that \(X_{\gamma _{t}}(s)\geq 0\) a.s. for any s[t,T]. □

For the comparison theorem we consider the following path-dependent BSVIE with jumps

$$ \begin{aligned} &Y_{\gamma_{t}}(s) \\&=\Phi(X_{\gamma_{t},s}) +\int_{s}^{T} f\left(s,r,X_{\gamma_{t},r-},Y_{\gamma_{t}}(r-)\right) +h\left(s,Z_{\gamma_{t}}(r,s), \int_{\mathbb{R}^{l}}U_{\gamma_{t}}(r,s,z)\nu(dz)\right)dr \\&\quad-\int_{s}^{T}Z_{\gamma_{t}}(s,r)dW(r) -\int_{s}^{T}\int_{\mathbb{R}^{l}} U_{\gamma_{t}}(s,r,z)\tilde{N}(dr,dz), \qquad s\in[t,T]. \end{aligned} $$
(45)

Theorem 7

Let \(f,\bar {f}:\Omega \times [t,T]^{2}\times D([0,T],\mathbb {R}^{d}) \times \mathbb {R}\rightarrow \mathbb {R}\) and \(h,\bar {h}:\Omega \times [t,T]\times D([0,T],\mathbb {R}^{d})\times \mathbb {R}\times \mathbb {R}\rightarrow \mathbb {R}\) satisfy assumptions (H2) and (H3) such that

$$ \begin{aligned} f\left(s,r,\gamma,y\right)&\geq \bar{f}\left(s,r,\gamma,y,\right), \quad\text{a.s.}, \\ h(s,z,u)&\geq\bar{h}(s,z,u),\quad\text{a.s.} \end{aligned} $$
(46)

for all \((s,r,\gamma,y,z,u)\in \Delta ^{c}[t,T]\times D([0,T],\mathbb {R}^{d})\times \mathbb {R}\times \mathbb {R}\times \mathbb {R}\). Moreover, \(y\mapsto f(s,r,\gamma,y)\in C^{1}(\mathbb {R};\mathbb {R})\) and \((z,u)\mapsto h(s,z,u)\in C^{1}(\mathbb {R}\times \mathbb {R};\mathbb {R})\) for any \(\gamma \in D([0,T],\mathbb {R})\) and (s,r)[t,T]2 with f y , h z , h u uniformly bounded and

$$|f_{y}(r,s,\gamma_{s-},y)\gamma(r-) -f_{y}(r,\bar{s},\bar{\gamma}_{\bar{s}-},y) \bar{\gamma}(r-)|^{2} \leq K\left(\sup_{s\in[t,r]}|\gamma(s)-\bar{\gamma}(s)|^{2}+|s-\bar{s}|^{2} \right). $$

Furthermore,

$$f_{y}(r,\bar{s},\gamma_{\bar{s}-},y)- f_{y}(r,s,\gamma_{s-},y)\geq0, $$

a.e. (ω,r)Ω×[t,s] for any \(\gamma \in D([0,T],\mathbb {R})\), \(s,\bar {s}\in [t,T]\) with \(\bar {s}\geq s\). For \(\Phi,\bar {\Phi }\in L^{2}_{\mathcal {F}_{T}^{t}}(D([0,T],\mathbb {R}^{d});\mathbb {R})\) we assume \(\Phi (\gamma)\geq \bar {\Phi }(\gamma)\) a.s. for all \(\gamma \in D([0,T],\mathbb {R}^{d})\).

Let \((Y_{\gamma _{t}}(\cdot),Z_{\gamma _{t}}(\cdot,\cdot), U_{\gamma _{t}}(\cdot,\cdot,\cdot))\) and \((\bar {Y}_{\gamma _{t}}(\cdot),\bar {Z}_{\gamma _{t}}(\cdot,\cdot), \bar {U}_{\gamma _{t}}(\cdot,\cdot,\cdot))\) be the adapted M-solution to BSVIE (45) corresponding to (f,h,Φ) and \((\bar {f},\bar {h},\bar {\Phi })\), respectively. Then, we get

$$ \begin{aligned} \mathbb{E}\left[\int_{t}^{T}Y_{\gamma_{t}}(u-)du\right] \geq \mathbb{E}\left[\int_{t}^{T}\bar{Y}_{\gamma_{t}}(u-)du\right]. \end{aligned} $$
(47)

Proof

Define for any (s,r)[t,T]2, \(h\in \mathbb {R}\setminus \{0\}\), the map \(l^{s,r}_{X_{\gamma _{t},r-}}:[0,1]\rightarrow \mathbb {R}\) by

$$l^{s,r}_{X_{\gamma_{t},r-}}(\theta) =\bar{Y}_{\gamma_{t}}(r-) +\theta\left(Y_{\gamma_{t}}(r-)-\bar{Y}_{\gamma_{t}}(r-)\right). $$

It follows that

$$\begin{aligned} & f\left(s,r,X_{\gamma_{t},r-},Y_{\gamma_{t}}(r-)\right)- f\left(s,r,X_{\gamma_{t},r-},\bar{Y}_{\gamma_{t}}(r-)\right) \\&=\int_{0}^{1} \frac{\partial}{\partial y}f\left(l^{s,r}_{X_{\gamma_{t},r-}}(\theta)\right)d\theta \cdot\left(Y_{\gamma_{t}}(r-)-\bar{Y}_{\gamma_{t}}(r-)\right). \end{aligned} $$

Similarly, we define for any s[t,T], \(h\in \mathbb {R}\setminus \{0\}\), a map \(g^{s}: [0,1]\rightarrow \mathbb {R}\times \mathbb {R}\) by

$$\begin{aligned} g^{s}(\theta)=&\left(\bar{Z}_{\gamma_{t}}(r,s)+\theta\left(Z_{\gamma_{t}}(r,s) -\bar{Z}_{\gamma_{t}}(r,s)\right),\right.\\ &\left.\int_{\mathbb{R}}\bar{U}_{\gamma_{t}}(r,s,z)+\theta\left(U_{\gamma_{t}}(r,s,z)-\bar{U}_{\gamma_{t}}(r,s,z)\right)\nu(dz)\right) \end{aligned} $$

and this yields

$$\begin{aligned} & h\left(s,Z_{\gamma_{t}}(r,s), \int_{\mathbb{R}}U_{\gamma_{t}}(r,s,z)\nu(dz)\right) -h\left(s,\bar{Z}_{\gamma_{t}}(r,s), \int_{\mathbb{R}}\bar{U}_{\gamma_{t}}(r,s,z)\nu(dz)\right) \\&=\int_{0}^{1}\frac{\partial}{\partial z}h\left(g^{s}(\theta)\right)d\theta\cdot \left(Z_{\gamma_{t}}(r,s)-\bar{Z}_{\gamma_{t},s-}(r,s)\right) \\&\quad+\int_{0}^{1}\frac{\partial}{\partial u}h\left(g^{s}(\theta)\right)d\theta\cdot \int_{\mathbb{R}}U_{\gamma_{t}}(r,s,z)-\bar{U}_{\gamma_{t}}(r,s,z) \nu(dz). \end{aligned} $$

Since \((Y_{\gamma _{t}}(\cdot),Z_{\gamma _{t}}(\cdot,\cdot), U_{\gamma _{t}}(\cdot,\cdot,\cdot))\), \((\bar {Y}_{\gamma _{t}}(\cdot),\bar {Z}_{\gamma _{t}}(\cdot,\cdot), \bar {U}_{\gamma _{t}}(\cdot,\cdot,\cdot))\) are adapted M-solutions to BSVIE (45) corresponding to (f,h,Φ) and \((\bar {f},\bar {h},\bar {\Phi })\), respectively, we get for s[t,T]

$$\begin{aligned} & Y_{\gamma_{t}}(s)-\bar{Y}_{\gamma_{t}}(s) \\&=\hat{\Phi}(X_{\gamma_{t},s}) +\int_{s}^{T}A_{0}(s,r,X_{\gamma_{t},r-}) \left(Y_{\gamma_{t}}(r-)-\bar{Y}_{\gamma_{t}}(r-)\right) +A_{1}(s) \left(Z_{\gamma_{t}}(r,s)-\bar{Z}_{\gamma_{t}}(r,s)\right)dr \\&\quad +\int_{s}^{T}\int_{\mathbb{R}}A_{2}(s) \left(U_{\gamma_{t}}(r,s,z)-\bar{U}_{\gamma_{t}}(r,s,z)\right) \nu(dz)dr -\int_{s}^{T}Z_{\gamma_{t}}(s,r)-\bar{Z}_{\gamma_{t}}(s,r)W(dr) \\&\quad -\int_{s}^{T}\int_{\mathbb{R}}U_{\gamma_{t}}(s,r,z) -\bar{U}_{\gamma_{t}}(s,r,z)\tilde{N}(dr,dz), \qquad s\in[t,T], \end{aligned} $$

where

$$\begin{aligned} &\hat{\Phi}(X_{\gamma_{t},s}) \\&= \Phi(X_{\gamma_{t},s})-\bar{\Phi}(X_{\gamma_{t},s}) +\int_{s}^{T} f\left(s,r,X_{\gamma_{t},r-},\bar{Y}_{\gamma_{t}}(r-)\right) -\bar{f}\left(s,r,X_{\gamma_{t},r-},\bar{Y}_{\gamma_{t}}(r-)\right)dr \\&+\int_{s}^{T}h\left(s,\bar{Z}_{\gamma_{t}}(r,s), \int_{\mathbb{R}^{l}}\bar{U}_{\gamma_{t}}(r,s,z) \nu(dz)\right) -\bar{h}\left(s,\bar{Z}_{\gamma_{t}}(r,s), \int_{\mathbb{R}^{l}}\bar{U}_{\gamma_{t}}(r,s,z) \nu(dz)\right)dr, \end{aligned} $$

and

$$\begin{aligned} &A_{0}(s,r,X_{\gamma_{t},r-})= \int_{0}^{1} \frac{\partial}{\partial y}f\left(l^{s,r}_{X_{\gamma_{t},r-}}(\theta)\right)d\theta, \quad A_{1}(s) = \int_{0}^{1}\frac{\partial}{\partial z}h\left(g^{s}(\theta)\right)d\theta, \\& A_{2}(s) = \int_{0}^{1}\frac{\partial}{\partial u}h\left(g^{s}(\theta)\right)d\theta. \end{aligned} $$

Pick any \(\varphi (\cdot)\in L^{2}_{\mathbb {F}}(\Omega \times D([0,T],\mathbb {R});\mathbb {R})\) with φ(γ)≥0, \(\gamma \in D([0,T];\mathbb {R})\), and consider a path-dependent linear FSVIE of the form

$$\begin{aligned} X_{\gamma_{t}}(s) &=\varphi(X_{\gamma_{t},s}) +\int_{t}^{s}A_{0}(r,s,X_{\gamma_{t},s-})X_{\gamma_{t}}(r-)dr +\int_{t}^{s}A_{1}(r)X_{\gamma_{t}}(r-)dW(r) \\&\quad+\int_{t}^{s}\int_{\mathbb{R}}A_{2}(r) X_{\gamma_{t}}(r-)\tilde{N}(dr,dz), \end{aligned} $$

for s[t,T]. Theorem 6 yields \(X_{\gamma _{t}}(s)\geq 0\) a.s. Moreover, it is \(\hat {\Phi }(X_{\gamma _{t},s})\geq 0\), a.s. The duality principle (Theorem 5) implies for s[t,T]

$$\begin{aligned} \mathbb{E}\left[\int_{t}^{T}\langle \varphi(X_{\gamma_{t},u-}), Y_{\gamma_{t}}(u-)-\bar{Y}_{\gamma_{t}}(u-)\rangle du\right] =\mathbb{E}\left[\int_{t}^{T}\langle X_{\gamma_{t}}(u-), \hat{\Phi}(X_{\gamma_{t},u-})\rangle du\right]\geq 0 \end{aligned} $$

and we obtain

$$\begin{aligned} \mathbb{E}\left[\int_{t}^{T}Y_{\gamma_{t}}(u-)du\right] \geq \mathbb{E}\left[\int_{t}^{T}\bar{Y}_{\gamma_{t}}(u-)du\right]. \end{aligned} $$

6 Path dependent BSVIEs and dynamic risk measures

As indicated in the introduction, BSVIEs are related to dynamic risk measures, see Yong (2007). If the risk measures are now based on path-dependent BSVIEs, the risk measures exhibit an additional dependence from an underlying process X. This process can be called a factor process or asset price process and may drive the value process Φ=Φ(X) as well as the generator of the BSVIE.

Then, differentiability enables the understanding of the sensitivity of the risk measure with respect to the underlying process X.

With the help of the comparison theorem, we are able to define a class of continuous-time dynamic risk measures by the solution of a path-dependent BSVIE with jumps. The following two definitions are taken from Yong (2007), which we modify for our purpose. Since the results of the previous sections can be applied, the proofs from Yong (2007) can also be modified.

Definition 3

A map \(\rho :L^{2}_{\mathcal {F}_{T}}(t,T)\rightarrow L^{2}_{\mathbb {F}}(t,T)\) is called a dynamic risk measure if the following hold

  1. 1.

    (Past Independence) For any \(\Phi (\cdot),\bar {\Phi }(\cdot) \in L^{2}_{\mathcal {F}_{T}^{t}}\left (D\left ([0,T],\mathbb {R}^{d}\right); \mathbb {R}\right)\), if

    $$ \Phi(\gamma_{s})=\bar{\Phi}(\gamma_{s}),\quad\text{a.s.}\ \omega\in\Omega,s\in[r,T], $$

    for some r[t,T], then

    $$ \rho(r;\Phi(\cdot))=\rho(r;\bar{\Phi}(\cdot)),\quad \text{a.s.}\ \omega\in\Omega. $$
  2. 2.

    (Monotonicity) For any \(\Phi (\cdot),\bar {\Phi }(\cdot)\in L^{2}_{\mathcal {F}_{T}^{t}}\left (D\left ([0,T],\mathbb {R}^{d}\right); \mathbb {R}\right)\), if

    $$ \Phi(\gamma)\leq\bar{\Phi}(\gamma), \quad\text{a.s.}\ \omega\in\Omega,\gamma\in D\left([0,T],\mathbb{R}^{d}\right), $$

    then

    $$ \rho(s;\Phi(\cdot))\geq\rho(s;\bar{\Phi}(\cdot)), \quad\text{a.s.}\ \omega\in\Omega,s\in[t,T]. $$

Definition 4

A dynamic risk measure \(\rho :L^{2}_{\mathcal {F}_{T}}(t,T)\rightarrow L^{2}_{\mathbb {F}}(t,T)\) is called a dynamic coherent risk measure if the following hold

  1. 1.

    (Positive homogeneity) For any \(\Phi (\cdot)\in L^{2}_{\mathcal {F}_{T}^{t}}\left (D\left ([0,T],\mathbb {R}^{d}\right); \mathbb {R}\right)\) and \(\lambda \in \mathbb {R}^{+}\),

    $$ \rho(s;\lambda\Phi(\cdot))=\lambda\cdot\rho(s;\Phi(\cdot)), \quad\text{a.s.}\ \omega\in\Omega,s\in[t,T]. $$
  2. 2.

    (Sub-additivity) For any \(\Phi (\cdot),\bar {\Phi }(\cdot)\in L^{2}_{\mathcal {F}_{T}^{t}}\left (D\left ([0,T],\mathbb {R}^{d}\right); \mathbb {R}\right)\)

    $$ \rho(s;\Phi(\cdot)+\bar{\Phi}(\cdot)) \leq \rho(s;\Phi(\cdot))+\rho(s;\bar{\Phi}(\cdot)), \quad\text{a.s.}\ \omega\in\Omega,s\in[t,T]. $$
  3. 3.

    (Translation invariance) There exists a deterministic integrable function r(·,·) such that for any \(\Phi (\cdot)\in L^{2}_{\mathcal {F}_{T}^{t}}\left (D\left ([0,T],\mathbb {R}^{d}\right); \mathbb {R}\right)\)

    $$ \rho(s;\Phi(\cdot)+c) =\rho(s;\Phi(\cdot))-ce^{-\int_{s}^{T}r(s,u)du}, \quad\text{a.s.}\ \omega\in\Omega,s\in[t,T]. $$

We denote

$$ \begin{aligned} \rho(s,\Phi(\cdot)) =\mathbb{E}\left[\int_{t}^{T}Y_{\gamma_{t}}(s-)ds\right], \end{aligned} $$
(48)

where \((Y_{\gamma _{t}}(\cdot),Z_{\gamma _{t}}(\cdot,\cdot), U_{\gamma _{t}}(\cdot,\cdot,\cdot))\) is the unique adapted M-solution to the following path-dependent BSVIE with jumps

$$\begin{aligned} &Y_{\gamma_{t}}(s) \\&=-\Phi(X_{\gamma_{t},s}) +\int_{s}^{T} f\left(s,r,X_{\gamma_{t},r-},Y_{\gamma_{t}}(r-)\right) +h\left(s,Z_{\gamma_{t}}(r,s), \int_{\mathbb{R}^{l}}U_{\gamma_{t}}(r,s,z)\nu(dz)\right)dr \\&\quad -\int_{s}^{T}Z_{\gamma_{t}}(s,r)dW(r) -\int_{s}^{T}\int_{\mathbb{R}^{l}} U_{\gamma_{t}}(s,r,z)\tilde{N}(dr,dz), \qquad s\in[t,T]. \end{aligned} $$

Proposition 1

  1. 1.

    Let f and h satisfy (H2) and (H3). Moreover, suppose that both are sub-additive, i.e.,

    $$\begin{aligned} f(s,r,\gamma,y_{1}+y_{2})&\leq f(s,r,\gamma,y_{1})+f(s,r,\gamma,y_{2}) \\ h(s,z_{1}+z_{2},u_{1}+u_{2})&\leq h(s,z_{1},u_{1})+h(s,z_{2},u_{2}) \end{aligned} $$

    for all \((s,r)\in [t,T]^{2},y_{1},y_{2}\in \mathbb {R}, z_{1},z_{2}\in \mathbb {R}^{d},u_{1},u_{2}\in \mathbb {R}\). Then, Φ(·)ρ(s;Φ(·)) is sub-additive, i.e.,

    $$ \rho(s;\Phi_{1}(\cdot)+\Phi_{2}(\cdot))\leq \rho(s;\Phi_{1}(\cdot))+\rho(s;\Phi_{2}(\cdot)),\:\text{a.e}. $$
  2. 2.

    If \(f:[t,T]^{2}\times D([0,T],\mathbb {R}^{d})\times \mathbb {R}\rightarrow \mathbb {R}\) and \(h:[t,T]\times \mathbb {R}^{d}\times \mathbb {R}\rightarrow \mathbb {R}\) are positively homogeneous, i.e.,

    $$\begin{aligned} f(s,r,\gamma,\lambda y)&=\lambda\cdot f(s,r,\gamma,y),\quad h(s,\lambda z,\lambda u)=\lambda\cdot h(s,z,u), \end{aligned} $$

    for any \(\lambda \in \mathbb {R}^{+}\), (s,r)[t,T]2, \(\gamma \in D([0,T],\mathbb {R}^{d})\), \(y\in \mathbb {R}\), \(z\in \mathbb {R}^{d}\) and \(u\in \mathbb {R}\) a.e., then Φ(·)ρ(s,Φ(·)) is also positively homogeneous.

  3. 3.

    If f is independent of \(\gamma \in D([0,T],\mathbb {R}^{d})\) and of the form f(s,r,y)=r(s,r)y with a deterministic integrable function r(·,·), then \(\Phi (\cdot)\rightarrow \rho (s;\Phi (\cdot))\) is transition invariant, i.e.,

    $$ \rho(s;\Phi(\cdot)+c) =\rho(s;\Phi(\cdot))-e^{\int_{s}^{T}r(s,u)du}, \quad\text{a.s.} \omega\in\Omega,s\in[t,T],\text{ for all} c\in\mathbb{R}. $$

    In particular, if f(s,r,γ,y)=0, then ρ(s,Φ(·)+c)=ρ(s,Φ(·))−c a.s. for any \(c\in \mathbb {R}\).

Theorem 8

Suppose f and h satisfy (H2) and (H3). Moreover, f(s,r,γ,y)=r(s,r)y, with r(·,·) being a bounded and deterministic integrable function. Then, ρ(·), defined by (48), is a dynamic coherent risk measure if f and h are positively homogeneous and sub-additive.

Proof

The proof is given by Proposition 1 and the comparison theorem for path-dependent BSVIEs with jumps (Theorem 7). □

References

Download references

Acknowledgments

The authors thank the editor and two anonymous referees for their helpful suggestions.

Author information

Authors and Affiliations

Authors

Contributions

Both authors read and approved the final manuscript.

Corresponding author

Correspondence to Jasmin A. L. Röder.

Ethics declarations

Competing interests

The authors declare that they have no competing interests.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License(http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Overbeck, L., Röder, J.L. Path-dependent backward stochastic Volterra integral equations with jumps, differentiability and duality principle. Probab Uncertain Quant Risk 3, 4 (2018). https://doi.org/10.1186/s41546-018-0030-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s41546-018-0030-2

Keywords