The Duality Principle for Multidimensional Optional Semimartingales

1. Introduction

In this study, our objective is to build upon the exploration of the duality principle in option pricing for multidimensional optional semimartingales, as introduced in Eberlein et al. (2009). The options under consideration have payoffs dependent on various factors, such as foreign exchange rates, interest rates, and multiple stocks, including examples like swap options and quanto options.

The pricing of these complex options often necessitates a comprehensive understanding of joint probabilities and involves intricate integration procedures. A suggested alternative approach, as outlined in Eberlein et al. (2009), involves transforming the original problem into its dual option pricing problem, rather than directly solving it. The subsequent resolution of this dual problem provides valuable insights into the valuation of options characterized by intricate dependencies and multiple influencing factors.

As stated earlier, the model we employ to describe the evolution of asset price processes is in the context of multidimensional optional semimartingale models. The family of optional semimartingales constitutes a diverse class of stochastic processes that includes cadlag semimartingales as a subset. Generally, such processes lack cadlag modifications, i.e., they are not inherently right-continuous with finite left-limits. Recently, numerous articles have explored optional processes and their relevance in financial and energy markets. For a comprehensive discussion on their application in the energy market, particularly addressing spike-related issues, please refer to Abdelghani et al. (2022).

Applying the outcomes derived from our study, we can establish duality relationships among swap options, quanto options, and standard call and put options within the framework of optional semimartingale models. This revelation carries notable implications, primarily manifesting in a substantial reduction in the computational complexity associated with the valuation of these financial instruments.

By leveraging our findings, we unveil a practical and efficient means of interrelating swap and quanto options with their standard call and put counterparts in optional semimartingale models. This not only enhances the understanding of the intricate connections between these financial instruments but also provides a tangible advantage in terms of computational efficiency when determining their respective valuations.

Extensive literature has delved into the duality principle of option pricing for multidimensional models. Notable contributions include studies by Margrabe (1978); Geman et al. (1995), and Gerber and Shiu (1996) for the Black-Scholes model, Eberlein and Papapantoleon (2005) for examining time-inhomogeneous Levy processes, Fajardo and Mordecki (2006) for Levy processes, and Schröder (2015) is focusing on cadlag semimartingales. In Eberlein et al. (2009) work, Eberlein, Papapantoleon, and Shiryaev establish the predictable characteristics triplet for one-dimensional semimartingales under the dual measure, employing the Esscher change of measure. On the contrary, limited attention has been given to optional semimartingales.

The paper is organized as follows. In Section 2, we will establish crucial definitions and notations that serve as the foundation for understanding the subsequent sections. Additionally, we will review pertinent results regarding the canonical decomposition of optional semimartingales. This includes an exploration of their quadruplet of predictable characteristics, and their the Laplace cumulant process. Section 3 contains details and lemmas tailored for Exponentially Optional Semimartingale models. We included the main result of this work in Section 4, where we present results concerning multidimensional optional semimartingales and the treatment of linear transformations applied to multidimensional semimartingales. In Section 5, we delve into the duality relation for optional processes, specifically examining the European, Margrabe, and quanto options. Alongside this, in Section 6, we provide explicit examples elucidating the application of the duality principles examined in the previous section.

2. Canonical Decomposition of Optional Semimartingales

Let’s begin by defining key terms and notations crucial for understanding the upcoming sections. Specifically, we need to clarify the concept of an unusual probability space. An unusual probability space is represented as $(\Omega ,\mathcal{F},\mathbf{F}={\left({\mathcal{F}}_t\right)}_{t\ge 0},\mathbf{P})$, wherein the natural filtration ${\mathcal{F}}_t$ is characterized by the absence of both right- and left-continuity. The completeness of the space arises from the fact that $\mathcal{F}$ incorporates all its P-null sets, and P serves as a completed probability measure. We define ${\mathbf{F}}_{+}=\left({\mathcal{F}}_{t^{+}}\right)$, with ${\mathcal{F}}_{t^{+}}={\bigcap }_{t\le s}{\mathcal{F}}_s$. Now, considering ${\mathbf{F}}^{*}:={\mathbf{F}}_{+}^P$ as the standard filtration—both right-continuous and complete—the construction of ${\mathbf{F}}_{+}^P$ involves enriching ${\mathbf{F}}_{+}$ with P-null sets.

Let ${\mathbb{R}}^d$ represent the d-dimensional Euclidean space. In this space, the scalar product between two vectors u and v, both elements of ${\mathbb{R}}^d$, is denoted as $\langle u,v\rangle =u^{\prime }v$, where $u^{\prime }$ denotes the transpose of the vector (or matrix) u.

Moving forward, the next step involves defining optional processes, which serve as the foundational processes under consideration in this paper. All definitions in this regard are drawn from Abdelghani and Melnikov (2020). To formulate the definition of optional processes, it is essential to establish the concept of an optional $\sigma$-algebra. In this context, we designate the $\sigma$-algebra $\mathcal{O}\left(F\right)$ on the interval $[0,\infty )\times \Omega$ as optional if it is generated by all right-continuous ${\mathcal{F}}_t$-adapted processes with left-hand limits.

Definition 1.

A random process $X=\left(X_t\right),t\in [0,\infty )$, is said to be optional if it is $\mathcal{O}\left(F\right)$-measurable. Optional processes are progressively measurable, and thus clearly measurable. In general, optional processes have right and left limits but are not necessarily continuous on the right, left, or otherwise.

Additionally, when it comes to optional stochastic processes, we can introduce the following definitions. The left process, denoted as $X_{-}={\left(X_{t_{-}}\right)}_{t\ge 0}$, where $X_{0_{-}}=X_0$ and $X_{t_{-}}={\lim }_{s\uparrow t}{X}_s$. Similarly, we have the right processes $X_{+}={\left(X_{t_{+}}\right)}_{t\ge 0}$. Also, the regular differential process $\Delta X={(\Delta X_t)}_{t\ge 0}$ where $\Delta X_t=X_t-X_{t-}$ and the forward differential process ${\Delta }^{+}X={\left({\Delta }^{+}{X}_t\right)}_{t\ge 0}$ where ${\Delta }^{+}{X}_t=X_{t+}-X_t$.

It is also crucial to understand the classification of stopping times, a topic that is thoroughly discussed in Chapter 2 of Abdelghani and Melnikov (2020). A solid understanding of these definitions is paramount for delving into the subsequent sections. Here, we would like to briefly define stopping (or Markov) times. Then we will restate two important categories of this random time which is important in the rest of the paper.

Definition 2.

Let $\left({\mathcal{F}}_t\right)$ be a filtration on $(\Omega ,\mathcal{F},P)$. A random variable T on Ω with values in $R_{+}\cup \{+\infty \}$ is called a stopping time or optional time of the filtration $\left({\mathcal{F}}_t\right)$ if, for all $t\in R_{+}$, the event $\{T\le t\}$ belongs to ${\mathcal{F}}_t$ and belongs to the trivial filtration. The set of all stopping times is denoted by $\mathcal{T}$.

Furthermore, Stopping times are also called Markov time. The simplest examples of them are $1_{T\le t}$ and every positive constant.

Definition 3.

T is called a wide sense stopping time if for all t,

$$\{T<t\}\in F_t$$

(1)

Also, we can say that T as a wide sense stopping time is optional in the wide sense if and only if it is optional to the family ${\mathit{F}}_{+}$. Another characterization of wide sense stopping times is, for all $t>0$, $T\wedge t$ is $F_t$-measurable.

We also require another classification of these stopping times, specifically, the concept of totally inaccessible stopping times, which will prove beneficial later on.

Definition 4.

$T>0$ is totally inaccessible if for every sequence ${\left(T_n\right)}_{n=1,2,\cdots }$ of stopping times that increases to T the event $\{{\lim }_{n\uparrow \infty }{T}_n=T,T<\infty \}$ has zero probability

Understanding the definition of a predictable process and $\sigma$-algebra is also valuable. We designate the $\sigma$-algebra $\mathcal{P}\left(\mathcal{F}\right)$ on $\Omega \times [0,\infty )$ as predictable if it is generated by all left-continuous ${\mathcal{F}}_t$-adapted processes with right-hand limits or by sets $\left\{\right(\omega ,t):S(\omega )<t\le T(\omega \left)\right\}$, where $S,T$ vary across all Markov times.

Definition 5.

A random process $Y=\left(Y_t\right)$, $t\in [0,\infty )$, whose trajectories have limits from the right is said to be predictable if $P\left(\mathcal{F}\right)$-measurable. Furthermore, it is termed strongly predictable if $Y_t\in P\left(\mathcal{F}\right)$ and $Y_t^{+}\in \mathcal{O}\left(\mathcal{F}\right)$ holds for all t.

Predictable processes are also optional, i.e., $\mathcal{P}\left(\mathcal{F}\right)\subseteq \mathcal{O}\left(\mathcal{F}\right)$. Moreover, the definition of strongly predictable processes implies that for every stopping time T, the random variables $X_{T^{+}}{1}_{T<\infty }$ and $X_{T^{-}}{1}_{T<\infty }$ are ${\mathcal{F}}_{T^{+}}$-measurable and ${\mathcal{F}}_{T^{-}}$-measurable, respectively. The set of strongly predictable processes is denoted by $P_s\left(\mathcal{F}\right)$ or $P_s$.

Let’s revisit the definition of optional martingales, highlighting the difference from the definition of regular martingales on ${\mathbf{F}}^{*}$ satisfying the usual conditions. For an adapted process $X=\left(X_t\right)$ with $t\in {\mathbb{R}}_{+}$, it qualifies as a martingale (or alternatively, a supermartingale or submartingale) concerning ${\mathbf{F}}^{*}$ if, for every $t\in {\mathbb{R}}_{+}$, $X_t$ is integrable, and the conditional expectation $X_s=\mathbb{E}\left[X_t\right|{\mathbf{F}}_s^{*}]$ (or $X_s\ge \mathbb{E}\left[X_t\right|{\mathbf{F}}_s^{*}]$, $X_s\le \mathbb{E}\left[X_t\right|{\mathbf{F}}_s^{*}]$) holds almost surely for all $s\le t$. Referring to the existence and uniqueness of optional martingales in Chapter 5 of Abdelghani and Melnikov (2020), we outline the definitions of optional martingale and optional local martingale here.

Definition 6.

We define $M={\left(M_t\right)}_{t\in {\mathbb{R}}_{+}}$, as an optional martingale (or alternatively, an optional supermartingale or an optional submartingale) if;

(a)

M is an optional process (i.e., $M\in O\left(\mathcal{F}\right)$),

(b)

The random variable $M_T·1_{\{T<\infty \}}$ is integrable for any $T\in \mathcal{T}$,

(c)

There exists an integrable random variable $\widehat{M}$ such that $M_T=\mathbb{E}\left[\widehat{M}\right|{\mathcal{F}}_T]$ (or, $M_T\ge \mathbb{E}\left[\widehat{M}\right|{\mathcal{F}}_T]$, $M_T\le \mathbb{E}\left[\widehat{M}\right|{\mathcal{F}}_T]$) almost surely on $\{T<\infty \}$ for any $T\in \mathcal{T}$.

Definition 7.

A process $M={\left(M_t\right)}_{t\ge 0}$ is called an optional local martingale if there exists a sequence $(R_n,M^{\left(n\right)}),n\in \mathbb{N},$ where $R_n$ are wide sense stopping times, $R_n\uparrow \infty$ a.s. and $M^{\left(n\right)}$ is an optional martingale, such that $M{1}_{[0,R_n]}=M^{\left(n\right)}{1}_{[0,R_n]}$ and the random variable $M_{R_{n+}}$ is integrable for any $n\in \mathbb{N}$.

Consider ${\mathcal{M}}^O\left(P\right)$ (${\mathcal{M}}_{loc}^O\left(P\right)$), representing the set of optional (local) martingales in relation to the probability measure P. Let $\mathcal{V}$ denote the collection of all $\mathbf{F}$-adapted processes with a finite variation where the variation for an optional process $A=\left(A_t\right)$, is given by

$$\mathrm{Var}{\left(A\right)}_t=\sum _{0\le s\le t}|A_{s^{+}}-A_s|+{\int }_{0^{+}}^t\left|d{A}^r\right|<\infty$$

(2)

where $A^r$ is a right-continuous finite-variation process, and the series is absolutely convergent. This discussion lays the groundwork for the exploration of optional semimartingales and their canonical decomposition, which is the next crucial step in this study.

Definition 8.

The stochastic process X is called an optional semimartingale denoted by the set $\mathcal{S}$, if

$$X=X_0+M+A,$$

where $M\in {\mathcal{M}}_{loc}^O,A\in \mathcal{V},A_0=M_0=0$ and $X_0$ is an ${\mathcal{F}}_0$-measurable finite random variable.

If this decomposition holds with a strongly predictable process A of locally integrable variation, then the optional semimartingale X attains the special designation of a special optional semimartingale, denoted by the set ${\mathcal{S}}_s$.

The canonical and component representation of semimartingales is also essential to our analysis of optional semimartingales. The canonical and component representation of optional semimartingales can be seen as a natural consequence of the decomposition

$$X=X_0+X^c+X^d+X^g$$

where $X^c$ is a continuous optional semimartingale with the decomposition $X^c=a+m$, where a is continuous and strongly predictable with locally integrable variation, and m is a continuous local martingale. The discrete optional semimartingale parts, $X^d=a^d+m^d$ and $X^g=a^g+m^g$, are expressible in terms of some underlying measures of right and left jumps, respectively.

Let’s proceed to define the characteristics of an optional semimartingale. Consider a d-dimensional optional semimartingale denoted as $Y={\left(Y_t\right)}_{0\le t\le T}$, where $Y={(Y_1,\dots ,Y_d)}^{\prime }$ with the decomposition $Y_i=A_i+M_i$, where $M_i\in {\mathcal{M}}_{\mathrm{loc}}^O$, $A_i\in \mathcal{V}$, $M_{i,0}=0$, $i=1,\dots ,d$. Let ${\left(T_n\right)}_{n\ge 1}$ and ${\left(U_n\right)}_{n\ge 1}$ be sequences of totally inaccessible stopping times and totally inaccessible wide sense stopping times, respectively. For $(E,\mathcal{E})$, consider the Lusin space where $E=({\mathbb{R}}^d\setminus \left\{0\right\})\cup \left\{{\delta }^d\right\}\cup \left\{{\delta }^g\right\}$ such that ${\delta }^d$, and ${\delta }^g$ are some supplementary points or are the set of processes with finite variation on any segment $[0,t]$, P-a.s.; $\mathcal{E}=\mathcal{B}\left(E\right)$. On the $\sigma$-algebra $\tilde{\mathcal{E}}=\mathcal{B}\left(E\right)\times \mathcal{E}$, define integer random measures on $({\mathbb{R}}_{+}\times E,\tilde{\mathcal{E}})$

$$\begin{array}{cccc}\hfill {\mu }^r(\Gamma )& =\sum _{n\ge 1}{\mathbf{1}}_{\Gamma }(T_n,{\beta }_{T_n}^r),\hfill & \hfill {\mu }^g(\Gamma )& =\sum _{n\ge 1}{\mathbf{1}}_{\Gamma }(U_n,{\beta }_{U_n}^g)\hfill \end{array}$$

where ${\mathbf{1}}_{\Gamma }(·)$ is an indicator function of a set $\Gamma \in \tilde{\mathcal{E}},$${\beta }_t^r=\Delta Y_t$ if $\Delta Y_t\ne 0$ and ${\beta }_t^r=\delta$ if $\Delta Y_t=0,{\beta }_t^g={\Delta }^{+}{Y}_t$ if ${\Delta }^{+}{Y}_t\ne 0,{\beta }_t^g=\delta$ if ${\Delta }^{+}{Y}_t=0,t>0.$ Under the unusual conditions on probability space $(\Omega ,\mathcal{F},\mathbf{F}={\left({\mathcal{F}}_t\right)}_{t\ge 0},\mathbf{P})$, any optional semimartingale Y can be decomposed as follows:

$$\begin{array}{c}{Y}_t=Y_0+B_t+C_t+\underset{]0,t]\times E}{\int }h\left(x\right)d({\mu }^r-{\nu }^r)+\underset{[0,t[\times E}{\int }h\left(x\right)d({\mu }^g-{\nu }^g)\\ +\underset{]0,t]\times E}{\int }(x-h\left(x\right))d{\mu }^r+\underset{[0,t[\times E}{\int }(x-h\left(x\right))d{\mu }^g\end{array}$$

or in short notation

$$Y=Y_0+B+C+\sum _{j=r,g}\left[h\left(x\right)*({\mu }^j-{\nu }^j)+(x-h\left(x\right))*{\mu }^j\right],$$

(3)

Here, $Y_0$ is ${\mathcal{F}}_0$-measurable random variable, $B\in \mathcal{V},B_0=0,$ and $C\in {\mathcal{M}}_{loc}^O,C_0=0,$ and it is continuous. The ${\nu }^j$ represent the compensators of ${\mu }^j$ (refer to chapter 7 of Abdelghani and Melnikov (2020) for a detailed definition of compensators). The function $h\left(x\right)$ is a truncation function, meaning it is a bounded function with compact support that behaves as $h\left(x\right)=x$ around the origin. As an example, one could choose $h\left(x\right)=x{1}_{\left\{\right|x|\le 1\}}$. Additionally, $W*\mu$ denotes the integral process, and $W*(\mu -\nu )$ denotes the stochastic integral with respect to the compensated random measure $(\mu -\nu )$.

The collection of the processes $B,C$ and the measures ${\nu }^r,{\nu }^g$ is called the (local) characteristics of the semimartingale Y with respect to the probability measure P. Let us call them a quadruplet and denote them by

$$Q\left(Y\right|P)=(B,C,{\nu }^r,{\nu }^g).$$

(4)

In addition, there exists an increasing predictable process A, predictable processes b, c and two transition kernels $F^r$ and $F^g$ from $(\Omega \times [0,T],\mathbf{P})$ into $(E,\mathcal{E})$ such that

$$\begin{array}{cc}\hfill B& =b·A,B_t={\int }_0^t{b}_{s}d{A}_s\hfill \\ \hfill C& =c·A,C_t={\int }_0^t{b}_{s}d{A}_s\hfill \\ \hfill {\nu }^j& =F^j\otimes A,{\nu }^j([0,t]\times E)={\int }_0^t{\int }_E{F}_s^j\left(dx\right)d{A}_s,j=r,g\hfill \end{array}$$

Every optional semimartingale Y with characteristics $Q\left(Y\right|P)=(B,C,{\nu }^r,{\nu }^g)$ can be associated with an optional Laplace cumulant process defined by

$$G\left(u\right)=\langle u,B\rangle +\frac{1}{2}\langle u,Cu\rangle +(e^{\langle u,x\rangle }-1-\langle u,h\left(x\right)\rangle )*{\nu }^r+(e^{\langle u,x\rangle }-1-\langle u,h\left(x\right)\rangle )*{\nu }^g$$

(5)

Furthermore, expressing the optional Laplace cumulant process can be achieved using the process A, represented as $G\left(u\right)=g\left(u\right)·A$ with

$$g\left(u\right)=\langle u,b\rangle +\frac{1}{2}\langle u,cu\rangle +(e^{\langle u,x\rangle }-1-\langle u,h\left(x\right)\rangle )*F^r+(e^{\langle u,x\rangle }-1-\langle u,h\left(x\right)\rangle )*F^g$$

(6)

Another significant distinction of the calculus of optional semimartingales on the unusual probability space, in comparison to regular semimartingales on the usual probability space, is their integration. To conclude this section and emphasize this difference, let us briefly discuss the results presented in Chapter 7 of Abdelghani and Melnikov (2020) concerning integration with respect to optional semimartingales.

The optional stochastic integral with respect to optional semimartingales X is defined in terms of the stochastic integrals with respect to its components A and m as

$$\begin{array}{cc}\hfill Y_t& ={\int }_0^t{h}_{s}d{X}_s=h\circ X_t={\int }_0^t{h}_{s}d{A}_s+{\int }_0^t{h}_{s}d{m}_s\hfill \\ \hfill & =h\circ A_t+h\circ m_t=h\circ A_t^r+h\circ A_t^g+h\circ m_t^r+h\circ m_t^g,\hfill \end{array}$$

(7)

Where

$$h\circ A_t^r={\int }_{0^{+}}^t{h}_{s^{-}}d{A}_s^r,h\circ A_t^g={\int }_0^{t^{-}}{h}_{s}d{A}_{s^{+}}^g,$$

$$h\circ m_t^r={\int }_{0^{+}}^t{h}_{s^{-}}d{m}_s^r,h\circ m_t^g={\int }_0^{t^{-}}{h}_{s}d{m}_{s^{+}}^g,$$

Moreover, given that $X=X^r+X^g=X^c+X^d+X^g$, we can express the integral with respect to X as

$$h\circ X=h\circ X^r+h\circ X^g=h\circ X^c+h\circ X^d+h\circ X^g,$$

where $X^r=A^r+m^r=A^c+m^c+A^d+m^d$ and $X^g=A^g+m^g$.

3. Exponentially Optional semimartingale models

In this section, we delve into specific details and lemmas tailored for Exponentially Optional Semimartingale models. The forthcoming discussions will cover various results, including equivalent statements for exponentially special optional semimartingales, a martingale version of the Lévy–Khintchine formula applicable to optional semimartingales, multiplicative decomposition for optional semimartingale, and a corollary addressing the uniqueness of the exponential optional compensator.

To start this section, we present definitions for a special optional semimartingale and an exponential optional compensator. This initial step aims to show differences between the findings in this section and those in the references Jacod and Shiryaev (1987) and Eberlein et al. (2009).

Definition 9.

An optional semimartingale X is called exponentially special optional semimartingales if $\exp (X-X_0)$ is a special optional semimartingale.

Definition 10.

Let X be an optional semimartingale. A strongly predictable process $V\in \mathcal{V}$ is called exponential optional compensator of X if $\exp (X-X_0-V)\in {\mathcal{M}}_{loc}^O\left(P\right).$

In the following lemma, our objective is to reformulate the outcome presented in Corollary II.2.42 of Jacod and Shiryaev (1987) for exponentially special optional semimartingales. This restatement serves as an extension of the original result by Jacob and Shiryaev.

Lemma 1.

Suppose X is an exponentially special optional semimartingale. The following statements are equivalent:

a) X is an optional semimartingale, and it admits the local quadruplet of characteristics $Q\left(X\right|P)=(B,C,{\nu }^r,{\nu }^g)$.

b) For each $z\in \mathbb{R}$, the process $e^{zX}-\left(e^{z{X}_{-}}\right)\circ G\left(z\right)\in {\mathcal{M}}_{loc}^O.$

Proof. $\left(a\right)\Rightarrow \left(b\right):$ We consider the following decomposition of $X:$

$$X=M+B+\sum _{s\le t}\left[\Delta X_s-h(\Delta X_s)\right]+\sum _{s<t}\left[{\Delta }^{+}{X}_s-h\left({\Delta }^{+}{X}_s\right)\right].$$

Applying the change of variables formula to $f\left(x\right)=e^{zx}$ and X, we get

$$\begin{array}{cc}\hfill e^{zX}-e^{z{X}_0}=& u{e}^{z{X}_{-}}\circ B+z{e}^{z{X}_{-}}\circ M\hfill \\ \hfill & -\frac{1}{2}{z}^2{e}^{z{X}_{-}}·C\hfill \\ \hfill & +e^{z{X}_{-}}\left(e^{zx}-1-zh\left(x\right)\right)*{\mu }^r\hfill \\ \hfill & +e^{zX}\left(e^{zx}-1-zh\left(x\right)\right)*{\mu }^g.\hfill \end{array}$$

(8)

Let us consider the right-hand side of (8); the second term is an optional local martingale; other terms are optional processes with a finite variation. However, since $e^{zX}$ is a special optional semimartingale, $e^{z{X}_{-}}\left(e^{zx}-1-zh\left(x\right)\right)*{\mu }^r+e^{zX}\left(e^{zx}-1-zh\left(x\right)\right)*{\mu }^g$ is actually a process of locally integrable variation. Thus,

$$e^{z{X}_{-}}\left(e^{zx}-1-zh\left(x\right)\right)*({\mu }^r-{\nu }^r)+e^{zX}\left(e^{zx}-1-zh\left(x\right)\right)*({\mu }^g-{\nu }^g)\in {\mathcal{M}}_{loc}^O,$$

and the result follows.

$\left(b\right)\Rightarrow \left(a\right):$ By hypothesis, $e^{zX}$ is an optional semimartingale for each $z\in \mathbb{R}.$ Then, since $\log \left(x\right)\in C^2,X$ is an optional semimartingale.

Let $(B^{*},C^{*},{\nu }^{r*},{\nu }^{g*})$ be a good version of the characteristics of $X.$ For each $z\in \mathbb{R}$ we associate to $(B^{*},C^{*},{\nu }^{r*},{\nu }^{g*})$ a process $G^{*}\left(z\right).$ We have proved the implication $\left(a\right)\Rightarrow \left(b\right),$ thus $e^{zX}-e^{zX}\circ G^{*}\left(z\right)\in {\mathcal{M}}_{loc}^O$. Then the hypothesis and the uniqueness of the canonical decomposition of the special semimartingale $e^{zX}$ show that $e^{zX}\circ G^{*}\left(z\right)=e^{zX}\circ G\left(z\right)$ up to an evanescent set. Using the orthogonality of $G^r$ and $G^g$ and integrating the processes $e^{-z{X}_{-}}$ and $e^{-zX},$ we obtain that $G\left(z\right)$ and $G^{*}\left(z\right)$ are indistinguishable. Therefore, the set N of all $\omega$ for which there exists $z\in \mathbb{Q}$ and $t\in {\mathbb{Q}}_{+}$ with $G{\left(z\right)}_{t\left(\omega \right)}\ne G^{*}{\left(z\right)}_{t\left(\omega \right)}$ is P-null.

Now, we observe that the optional cumulant function is continuous and, as a result, completely characterized by its values on $\mathbb{Q}.$ Consequently, outside of the set $N,$ we have $B_t^{*}:=B_t^{r*}+B_t^{g*},B_t:=B^r+B^g,C_t^{*}=C_t$ and ${\nu }^{j*}([0,t]\times ·)={\nu }^j([0,t]\times ·),j=r,g,$ for all $t\in {\mathbb{Q}}_{+}$ (see Gnedenko and Kolmogorov (1954)). This relationship also holds for all $t\in \mathbb{R}+$ due to the right-continuity of $B^{r*},B^r,C^{*},C,{\nu }^{r*},{\nu }^r,$ and the left-continuity of $B^{g*},B^g,{\nu }^{g*},{\nu }^g.$ Therefore, $(B,C,{\nu }^r,{\nu }^g)$ also serves as a version of the local characteristics of X. □

A significant finding in Jacod and Shiryaev (1987) (II.2.48 Corollary) and Eberlein et al. (2009) presents a martingale version of the Lévy–Khintchine formula specifically tailored for regular semimartingales. In this context, our goal is to broaden this result to encompass optional semimartingales.

Corollary 1.

Suppose X is an exponentially special optional semimartingale and $\mathcal{E}\left(G\right(z\left)\right)\ne 0$ for all $z\in \mathbb{R}$. If for each $z\in \mathbb{R},\frac{e^{zX}}{\mathcal{E}\left(G\right(z\left)\right)}\in {\mathcal{M}}_{loc}^O$. Then X is an optional semimartingale with characteristics $(B,C,{\nu }^r,{\nu }^g).$

Proof. 

Denote $M\left(z\right)=\frac{e^{zX}}{\mathcal{E}\left(G\right(z\left)\right)}\in {\mathcal{M}}_{loc}^O.$ Then by Lemma 5.2.3 in Pak (2021) we get

$$\begin{array}{cc}\hfill e^{zX}=& \mathcal{E}\left(G\right(z\left)\right)M\left(z\right)\hfill \\ \hfill =& e^{z{X}_0}+\mathcal{E}\left(G\left(z\right)\right)·M_t^r\left(z\right)+M_{-}\left(z\right)·{\mathcal{E}}_t^r\left(G\left(z\right)\right)\hfill \\ \hfill & +{\mathcal{E}}_{+}\left(G\left(z\right)\right)\odot M_{t+}^g\left(z\right)+M\left(z\right)\odot {\mathcal{E}}_{t+}^g\left(G\left(z\right)\right)\hfill \\ \hfill =& e^{z{X}_0}+\mathcal{E}\left(G\left(z\right)\right)·M_t^r\left(z\right)+M_{-}\left(z\right){\mathcal{E}}_{-}^r\left(G\left(z\right)\right)·G_t^r\left(z\right)\hfill \\ \hfill & +{\mathcal{E}}_{+}\left(G\left(z\right)\right)\odot M_{t+}^g\left(z\right)+M\left(z\right){\mathcal{E}}^r\left(G\left(z\right)\right)\odot G_{t+}^g\left(z\right)\hfill \\ \hfill =& e^{z{X}_0}+\mathcal{E}\left(G\left(z\right)\right)·M_t^r\left(z\right)+{\mathcal{E}}_{+}\left(G\left(z\right)\right)\odot M_{t+}^g\left(z\right)\hfill \\ \hfill & +e^{zX}\circ G_t\left(z\right)\hfill \end{array}$$

Note that we have used the definition of $\mathcal{E}\left(G\right(z\left)\right)$ above and the notation ⊙ for the optional integral with respect to the optional martingale process. It follows from the above that $e^{zX}-e^{z{X}_{-}}\circ G\left(z\right)\in {\mathcal{M}}_{loc}^O.$ The result follows from $\left(b\right)\Rightarrow \left(a\right)$ of Lemma 1. □

In addition to the canonical decomposition for an optional semimartingale, we want to establish that every optional semimartingale also allows for a multiplicative decomposition, as demonstrated in the following theorem.

Theorem 1.

Let X be an optional semimartingale with $X_0=1,$ such that $X,X_{-}$ and $X_{+}$ take their values in $(0,\infty ).$ Then X admits a multiplicative decomposition $X=LD,$ where $L\in {\mathcal{M}}_{loc}^O$, $L>0$ and $D\in \mathcal{V}$ is a positive strongly predictable process and $L_0=D_0=1,$ if and only if X is a special semimartingale.

In this case, the multiplicative decomposition is unique (up to evanescence), and is given as follows, where $X=1+M+A$ is the canonical (additive) decomposition of $X,$ and $H^r:=\frac{1}{X_{-}+\Delta A}$ and $H^g:=\frac{1}{X+{\Delta }^{+}A}$ are necessarily locally bounded and positive:

$$L=\mathcal{E}\left(H\circ M\right),D=\frac{1}{\mathcal{E}\left(-H\circ A\right)}$$

(9)

Proof. 

If a multiplicative decomposition $X=LD$ exists, by Lemma 5.2.3 in Pak (2021) we have

$$\begin{array}{c}\hfill X=LD=1+{\int }_{]0,t]}{D}_{s}d{L}_s^r+{\int }_{]0,t]}{L}_{s-}d{D}_s^r+{\int }_{[0,t[}{D}_{s+}d{L}_{s+}^g+{\int }_{[0,t[}{L}_{s}d{D}_{s+}^g.\end{array}$$

(10)

But ${\int }_{]0,t]}{D}_{s}d{L}_s^r+{\int }_{[0,t[}{L}_{s}d{D}_{s+}^g$ is an optional local martingale and ${\int }_{]0,t]}{L}_{s-}d{D}_s^r+{\int }_{[0,t[}{L}_{s}d{D}_{s+}^g$ is a strongly predictable process with finite variation, thus X is an optional special semimartingale.

Suppose that we have two multiplicative decompositions $X=LD=\widehat{L}\widehat{D},$ for which we can write (10). The uniqueness of the canonical decomposition yields

$${\int }_{]0,t]}{L}_{s-}d{D}_s^r+{\int }_{[0,t[}{L}_{s}d{D}_{s+}^g={\int }_{]0,t]}{\widehat{L}}_{s-}d{\widehat{D}}_s^r+{\int }_{[0,t[}{\widehat{L}}_{s}d{\widehat{D}}_{s+}^g(a.s.).$$

Since $D^r$ and $D^g$ are orthogonal, we have

$${\int }_{]0,t]}{L}_{s-}d{D}_s^r={\int }_{]0,t]}{\widehat{L}}_{s-}d{\widehat{D}}_s^r,{\int }_{[0,t[}{L}_{s}d{D}_{s+}^g={\int }_{[0,t[}{\widehat{L}}_{s}d{\widehat{D}}_{s+}^g(a.s.).$$

Since $L,\widehat{L},D,\widehat{D}$ are positive and $L/\widehat{L}=\widehat{D}/D,$ we deduce that

$${\int }_{]0,t]}1/D_{s-}d{D}_s^r={\int }_{]0,t]}1/{\widehat{D}}_{s-}d{\widehat{D}}_s^r$$

and

$${\int }_{[0,t[}1/D_{s}d{D}_{s+}^g={\int }_{[0,t[}1/{\widehat{D}}_{s}d{\widehat{D}}_{s+}^g$$

then if

$$U={\int }_{]0,t]}1/D_{s-}d{D}_s^r+{\int }_{[0,t[}1/D_{s}d{D}_{s+}^g,$$

we also have

$${\int }_{]0,t]}1/{\widehat{D}}_{s-}d{\widehat{D}}_s^r+{\int }_{[0,t[}1/{\widehat{D}}_{s}d{\widehat{D}}_{s+}^g=U,$$

which means, we have $D=\widehat{D}=\mathcal{E}\left(U\right),$ and this proves the uniqueness.

It remains to prove the following: if $X=1+M+A$ is the canonical decomposition of the optional special semimartingale X (with $X>0,X_{-}>0$ and $X_{+}>0$), and if $N=H\circ M$ and $B\circ A,$ then the strongly predictable process H is both locally bounded and positive. Moreover, $X=LD,$ where $L=\mathcal{E}\left(N\right)$ and $D=1/\mathcal{E}(-B)$ (as defined in (9)). Note that N, and L are optional local martingale, while B, and D are strongly predictable with finite variation. First, we observe that

$${}^{p}X={}^{p}A+{}^{p}M=A_{-}+\Delta A+M_{-}=X_{-}+\Delta A,$$

where the notation ${ }^p$ denotes a predictable projection of the process. Then for any predictable stopping time $T,$ and since $X>0,$ we get $X_{T-}+\Delta A_T=\mathbf{E}\left(X_T\right|{\mathcal{F}}_T)>0$ a.s. on $\{T<\infty \}.$ Then Theorem 2.4.52 in Abdelghani and Melnikov (2020) yields that $X_{-}+\Delta A>0$ outside P-null set, and we deduce $0<H^r<\infty .$ Next, notice that

$${}^{o}X_{+}={}^{o}A_{+}+{}^{o}M_{+}=A+{\Delta }^{+}A+M=X+{\Delta }^{+}A,$$

where ^o denotes an optional projection of the process. Then for any stopping time $T,$ and since $X>0$, we get $X_T+{\Delta }^{+}{A}_T=\mathbf{E}\left(X_{T+}\right|{\mathcal{F}}_T)>0$ a.s. on $\{T<\infty \}.$ Then Theorem 2.4.52 in Abdelghani and Melnikov (2020) yields that $X+{\Delta }^{+}A>0$ outside P-null set, and we have we deduce $0<H^g<\infty .$ Furthermore, Lemma $H^r$ and $H^g$ are locally bounded because they both are strongly predictable. Since we know that H is locally bounded we can define L and D as above, and it remains to prove that $X=LD=\mathcal{E}\left(N\right)/\mathcal{E}(-B).$

Observe that $\Delta B=\frac{\Delta A}{X_{-}+\Delta A}<1$ and ${\Delta }^{+}B=\frac{{\Delta }^{+}A}{X+{\Delta }^{+}A}<1$ identically, hence $\mathcal{E}(-B)>0$ and $\mathcal{E}{(-B)}_{-}>0,\mathcal{E}{(-B)}_{+}>0.$ Then we apply the change of variables formula to the function $f(x,y)=x/y,$ or rather to a $C^2$ function coinciding with f for y outside an arbitrary small neighborhood of $0,$ to obtain with $D^{*}=1/D$ and $\widehat{X}=LD:$

$$\begin{array}{cc}\hfill {\widehat{X}}_t=& 1+{\widehat{X}}_{-}·N_t^r+{\widehat{X}}_{-}·B_t^r+\sum _{s\le t}\left(\frac{1+\Delta N_s}{1-\Delta B_s}-1-\Delta N_s-\Delta B_s\right){\widehat{X}}_{s-}\hfill \\ \hfill & +\widehat{X}\odot N_{t+}^g+\widehat{X}\odot B_{t+}^g+\sum _{s<t}\left(\frac{1+{\Delta }^{+}{N}_s}{1-{\Delta }^{+}{B}_s}-1-{\Delta }^{+}{N}_s-{\Delta }^{+}{B}_s\right){\widehat{X}}_s\hfill \\ \hfill =& 1+{\widehat{X}}_{-}·N_t^r+{\widehat{X}}_{-}·B_t^r+\sum _{s\le t}{\widehat{X}}_{s-}\frac{\Delta N_s+\Delta B_s}{1-\Delta B_s}\Delta B_s\hfill \\ \hfill & +\widehat{X}\odot N_{t+}^g+\widehat{X}\odot B_{t+}^g+\sum _{s<t}{\widehat{X}}_{s-}\frac{{\Delta }^{+}{N}_s+{\Delta }^{+}{B}_s}{1-{\Delta }^{+}{B}_s}{\Delta }^{+}{B}_s\hfill \\ \hfill =& 1+{\widehat{X}}_{-}\left(1+\frac{\Delta B_s}{1-\Delta B_s}\right)·N_t^r+{\widehat{X}}_{-}\left(1+\frac{\Delta B_s}{1-\Delta B_s}\right)·B_t^r\hfill \\ \hfill & +\widehat{X}\left(1+\frac{{\Delta }^{+}{B}_s}{1-{\Delta }^{+}{B}_s}\right)\odot N_{t+}^g+\widehat{X}\left(1+\frac{{\Delta }^{+}{B}_s}{1-{\Delta }^{+}{B}_s}\right)\odot B_{t+}^g\hfill \\ \hfill =& 1+\frac{{\widehat{X}}_{-}{H}^r}{1-\Delta B}·X_t^r+\frac{\widehat{X}{H}^g}{1-{\Delta }^{+}B}\odot X_{t+}^g\hfill \\ \hfill =& 1+\frac{{\widehat{X}}_{-}}{X_{-}}·X_t^r+\frac{\widehat{X}}{X}\odot X_{t+}^g.\hfill \end{array}$$

Note that $\Delta B=H^r\Delta A,{\Delta }^{+}B=H^g{\Delta }^{+}A$ and $1-H^r\Delta A=H^r{X}_{-},1-H^g{\Delta }^{+}A=H^{g}X,$ hence $\frac{H^r}{1-\Delta B}=\frac{1}{X_{-}}$ and $\frac{H^g}{1-{\Delta }^{+}B}=\frac{1}{X}.$ Then

$$\frac{1}{{\widehat{X}}_{-}}·{\widehat{X}}^r+\frac{1}{\widehat{X}}\odot {\widehat{X}}^g=\frac{1}{X_{-}}·X^r+\frac{1}{X}\odot X^g$$

Thus, similarly to how we proved in the beginning, we deduce that $\widehat{X}=X.$ □

To conclude this section, we present a corollary addressing the uniqueness of the exponential optional compensator.

Corollary 2.

A real-valued optional semimartingale X has a unique (up to indistinguishably) exponential optional compensator V if X is exponentially special optional semimartingales.

Proof. 

Suppose that X is exponentially special optional semimartingales. By Theorem 1, there exists a unique positive process $D\in \mathcal{V}$ such that $D_0=1$ and $L:=\frac{\exp (X-X_0)}{D}\in {\mathcal{M}}_{loc}^O.$ Since $\exp (X-X_0)>0$ and $\frac{\exp (X_{-}-X_0)}{D_{-}}=L_{-}<\infty ,$ we have $D_{-}>0.$ Therefore, $V:=\log \left(D\right)\in \mathcal{V}.$ □

4. Multidimensional optional dual measures

In this section, our objective is to extend the findings presented in Eberlein et al. (2009) to the case of optional semimartingales. Eberlein, Papapantoleon, and Shiryaev, in their study, provided a characterization for one-dimensional semimartingales. These semimartingales were defined as scalar products of the driving semimartingale and d-dimensional vectors, under a suitable equivalent probability measure. Our goal is to generalize this result for optional semimartingales H, emphasizing its relevance to interest rate modeling through the optional stochastic integral in the Esscher transform.

Let $L\left(H\right)$ represent the space of predictable integrable processes concerning the semimartingale H. As previously defined in Section 2, $v^{\prime }\circ H$ signifies the optional stochastic integral of $v^{\prime }$ with respect to the optional semimartingale H. Additionally, $\tilde{G}\left(v\right)$ is used to denote the logarithmic transform of the cumulant process $G\left(v\right)$, expressed as $\mathcal{E}\left(G\left(v\right)\right)=\exp \left(\tilde{G}\left(v\right)\right)$.

Theorem 2.

Let H be an ${\mathbb{R}}^d$-valued optional semimartingale with characteristics $Q\left(H\right|P)=(B,C,{\nu }^r,{\nu }^g).$ Let u be a vector in ${\mathbb{R}}^d.$ Consider an ${\mathbb{R}}^d$-valued strongly predictable process $v,$ such that $v\in L\left(H\right)$ and $v^{\prime }\circ H$ is exponentially special optional semimartingale. Define the measure $P_v$ via the Radon-Nikodym derivative

$$\frac{d{P}_v}{dP}=\exp \left(v^{\prime }\circ H_T-\tilde{G}{\left(v\right)}_T\right),$$

assuming that $e^{v^{\prime }\circ H-\tilde{G}\left(v\right)}\in {\mathcal{M}}^O\left(P\right)$. Then the process $H^u$ with $H^u=\langle u,H\rangle =u^{\prime }H$ is a 1-dimensional optional semimartingale with characteristics $Q\left(H^u\right|P_v)=(B^u,C^u,{\nu }^{ur},{\nu }^{ug})$ of the form

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill B^u=& u^{\prime }B+u^{\prime }dCv+\sum _{j=r,g}\left(h\left(u^{\prime }x\right)\frac{e^{v^{\prime }x}}{1+W^j\left(v\right)}-u^{\prime }h\left(x\right)\right)*{\nu }^j,\hfill \\ \hfill C^u=& u^{\prime }Cu,\hfill \\ \hfill {\nu }^{ur}(\omega & ,t,E)=1_E\left(u^{\prime }x\right)\frac{e^{v^{\prime }x}}{1+W^r\left(v\right)}*{\nu }^r,\hfill \\ \hfill {\nu }^{ug}(\omega & ,t,E)=1_E\left(u^{\prime }x\right)\frac{e^{v^{\prime }x}}{1+W^g\left(v\right)}*{\nu }^g,\hfill \end{array}\end{array}$$

(11)

where $E\in \mathcal{B}(\mathbb{R}\setminus \{0\left\}\right)$, and $W^j{\left(v\right)}_t=\int (e^{v^{\prime }w}-1){\nu }^j(\left\{t\right\}\times dw),j=r,g$.

Proof. 

Using Corollary 1 in Section 3, it is sufficient to show that

$$\frac{e^{z{H}^u}}{\mathcal{E}\left(G^u\left(z\right)\right)}\in {\mathcal{M}}_{loc}^O\left(P_v\right)$$

(12)

for all $z\in \mathcal{Z},$ where $\mathcal{Z}\subseteq \mathbb{R}$ open, such that

$$1_{\left\{\right|y|>1\}}{e}^{zy}*{\nu }^{ur}\in \mathcal{V},1_{\left\{\right|y|>1\}}{e}^{zy}*{\nu }^{ug}\in \mathcal{V}.$$

(13)

Here $G^u$ denotes the Laplace cumulant process associated with the characteristics

$$(B^u,C^u,{\nu }^{ur},{\nu }^{ug})$$

Using Lemma 3.1 in Gasparyan (1988), (12) is equivalent to

$$Z\frac{e^{z{H}^u}}{\mathcal{E}\left(G^u\left(z\right)\right)}=\frac{e^{v^{\prime }\circ H}}{\mathcal{E}\left(G\right(v\left)\right)}\frac{e^{z{H}^u}}{\mathcal{E}\left(G^u\left(z\right)\right)}\in {\mathcal{M}}_{loc}^O\left(P\right)$$

(14)

Furthermore, using (11), condition (13) transforms to

$$1_{\left\{\right|u^{\prime }x|>1\}}{e}^{{(zu+v)}^{\prime }x}*{\nu }^r\in \mathcal{V},1_{\left\{\right|u^{\prime }x|>1\}}{e}^{{(zu+v)}^{\prime }x}*{\nu }^g\in \mathcal{V},$$

(15)

because $W^r\left(v\right)$ and $W^g\left(v\right)$ do not depend on $x.$ Now, the exponent of the numerator in (14) is

$$v^{\prime }\circ H+z{H}^u=v^{\prime }\circ H+z{u}^{\prime }H={(v+zu)}^{\prime }\circ H,$$

which has the unique exponential optional compensator by Corollary 2 in Section 3

$$G(v+zu)$$

under the measure $P,$ for the cumulant defined in (5); that is,

$$\frac{e^{{(v+zu)}^{\prime }\circ H}}{\mathcal{E}\left(G\right(v+zu\left)\right)}\in {\mathcal{M}}_{loc}^O\left(P\right).$$

(16)

Therefore, to complete the proof, it suffices to show

$$\mathcal{E}\left(G\left(v\right)\right)\mathcal{E}\left(G^u\left(z\right)\right)=\mathcal{E}\left(G(v+zu)\right).$$

(17)

Now, using the multiplication rule for optional stochastic exponential (see Lemma 8.1.2, Abdelghani and Melnikov (2020)) we get

$$\mathcal{E}\left(G\left(v\right)\right)\mathcal{E}\left(G^u\left(z\right)\right)=\mathcal{E}(G\left(v\right)+G^u\left(z\right)+[G\left(v\right),G^u\left(z\right)]).$$

Hence, (17) reduces to showing that

$$G\left(v\right)+G^u\left(z\right)+[G\left(v\right),G^u\left(z\right)]=G(v+zu).$$

(18)

On the right-hand side of (18), we have

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill G(v+zu)=& {(v+zu)}^{\prime }·B+\frac{1}{2}\int {(v+zu)}^{\prime }dC(v+zu)\hfill \\ \hfill & +\sum _{j=r,g}\left(e^{{(v+zu)}^{\prime }x}-1-{(v+zu)}^{\prime }h\left(x\right)\right)*{\nu }^j\hfill \end{array}\end{array}$$

(19)

Where operator · is the lebesgue integral. On the left-hand side, we have similarly

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill G\left(v\right)=& v^{\prime }·B+\frac{1}{2}\int v^{\prime }dCv+\sum _{j=r,g}\left(e^{v^{\prime }x}-1-v^{\prime }h\left(x\right)\right)*{\nu }^j.\hfill \end{array}\end{array}$$

(20)

The second term on the left-hand side of (18), using (11), is

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill G^u\left(z\right)=& z{B}^u+\frac{1}{2}{z}^2{C}^u+\sum _{j=r,g}(e^{zy}-1-zh\left(y\right))*{\nu }^{uj}\hfill \\ \hfill =& z{u}^{\prime }B+\int z{v}^{\prime }dCu+\sum _{j=r,g}z\left(h\left(u^{\prime }x\right)\frac{e^{v^{\prime }x}}{1+W^j\left(v\right)}-u^{\prime }h\left(x\right)\right)*{\nu }^j\hfill \\ \hfill & +\frac{1}{2}{z}^2{u}^{\prime }Cu+\sum _{j=r,g}\left(e^{z{u}^{\prime }x}-1-zh\left(u^{\prime }x\right)\right)\frac{e^{v^{\prime }x}}{1+W^j\left(v\right)}*{\nu }^j\hfill \\ \hfill =& z{u}^{\prime }B+\int z{v}^{\prime }dCu+\frac{1}{2}{z}^2{u}^{\prime }Cu\hfill \\ \hfill & +\sum _{j=r,g}\left(\frac{e^{{(zu+v)}^{\prime }x}-e^{v^{\prime }x}}{1+W^j\left(v\right)}-z{u}^{\prime }h\left(x\right)\right)*{\nu }^j.\hfill \end{array}\end{array}$$

(21)

The last term on the left-hand side of (18) is

$$[G\left(v\right),G^u\left(z\right)]=\sum _{t\le ·}\Delta G{\left(v\right)}_t\Delta G^u{\left(z\right)}_t+\sum _{t<·}{\Delta }^{+}G{\left(v\right)}_t{\Delta }^{+}{G}^u{\left(z\right)}_t$$

(22)

since $G,G^u$ are optional processes of finite variation. Now, we proceed as follows. First, we show that the drift terms on the left-hand side and right-hand side of (18) are equal. Then we show that the diffusive terms are equal. Finally, we prove that the jump terms on the left-hand side and right-hand side are equal. The drift term and the diffusive term are rather easy to handle; indeed from (20) and (21), we have that the drift term of the left-hand side of (18) is

$$v^{\prime }·B+z{u}^{\prime }B=(v+z{u}^{\prime }B)·B$$

(23)

Similarly, the diffusive term of the left-hand side of (18) is

$$\frac{1}{2}\left(\int v^{\prime }dCv+\int z{v}^{\prime }dCu+\int z{u}^{\prime }dCv+z^2{u}^{\prime }Cu\right)=\frac{1}{2}\int {(v+zu)}^{\prime }dC(v+zu),$$

(24)

since the matrix C is symmetric. Hence, both terms agree with the right-hand side.

The jump terms are more difficult to manipulate due to the presence of the fixed times of discontinuity for the optional semimartingale $H,$ which entails that the Laplace cumulant process is discontinuous. Regarding the fixed times of discontinuity, we have

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill \Delta G{\left(v\right)}_t=& \int (e^{v^{\prime }x}-1){\nu }^r(\left\{t\right\}\times dx)=W^r{\left(v\right)}_t,\hfill \\ \hfill {\Delta }^{+}G{\left(v\right)}_t=& \int (e^{v^{\prime }x}-1){\nu }^g(\left\{t\right\}\times dx)=W^g{\left(v\right)}_t,\hfill \end{array}\end{array}$$

(25)

and

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill \Delta G^u{\left(z\right)}_t=& \int (e^{zy}-1){\nu }^{ur}(\left\{t\right\}\times dy)\hfill \\ \hfill =& \int (e^{z{u}^{\prime }x}-1)\frac{e^{v^{\prime }x}}{1+W^r{\left(v\right)}_t}{\nu }^r(\left\{t\right\}\times dx)\hfill \\ \hfill =& \int (e^{{(v+zu)}^{\prime }x}-e^{v^{\prime }x})\frac{1}{1+W^r{\left(v\right)}_t}{\nu }^r(\left\{t\right\}\times dx).\hfill \end{array}\end{array}$$

(26)

Similarly,

$$\begin{array}{c}\hfill {\Delta }^{+}{G}^u{\left(z\right)}_t=\int (e^{{(v+zu)}^{\prime }x}-e^{v^{\prime }x})\frac{1}{1+W^g{\left(v\right)}_t}{\nu }^g(\left\{t\right\}\times dx).\end{array}$$

(27)

An important observation here is that $W\left(v\right)$ does not depend on the integrating variable $x.$ Hence, we can pull it out of the integration, and get

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill \Delta G^u{\left(z\right)}_t=& \frac{1}{1+W^r{\left(v\right)}_t}\int (e^{{(v+zu)}^{\prime }x}-e^{v^{\prime }x}){\nu }^r(\left\{t\right\}\times dx),\hfill \\ \hfill {\Delta }^{+}{G}^u{\left(z\right)}_t=& \frac{1}{1+W^g{\left(v\right)}_t}\int (e^{{(v+zu)}^{\prime }x}-e^{v^{\prime }x}){\nu }^g(\left\{t\right\}\times dx).\hfill \end{array}\end{array}$$

(28)

Hence, for the last term in (18), we can calculate further using (28)

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill [G\left(v\right),G^u\left(z\right)]=& \sum _{t\le ·}\Delta G{\left(v\right)}_t\Delta G^u{\left(z\right)}_t+\sum _{t<·}{\Delta }^{+}G{\left(v\right)}_t{\Delta }^{+}{G}^u{\left(z\right)}_t\hfill \\ \hfill =& \sum _{t\le ·}\frac{W^r{\left(v\right)}_t}{1+W^r{\left(v\right)}_t}\int (e^{{(v+zu)}^{\prime }x}-e^{v^{\prime }x}){\nu }^r(\left\{t\right\}\times dx)\hfill \\ \hfill & +\sum _{t<·}\frac{W^g{\left(v\right)}_t}{1+W^g{\left(v\right)}_t}\int (e^{{(v+zu)}^{\prime }x}-e^{v^{\prime }x}){\nu }^g(\left\{t\right\}\times dx)\hfill \\ \hfill =& \sum _{t\le ·}\frac{1+W^r{\left(v\right)}_t-1}{1+W^r{\left(v\right)}_t}\int (e^{{(v+zu)}^{\prime }x}-e^{v^{\prime }x}){\nu }^r(\left\{t\right\}\times dx)\hfill \end{array}\end{array}$$

(29)

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill & +\sum _{t<·}\frac{1+W^g{\left(v\right)}_t-1}{1+W^g{\left(v\right)}_t}\int (e^{{(v+zu)}^{\prime }x}-e^{v^{\prime }x}){\nu }^g(\left\{t\right\}\times dx)\hfill \\ \hfill =& \sum _{t\le ·}(\int (e^{{(v+zu)}^{\prime }x}-e^{v^{\prime }x}){\nu }^r(\left\{t\right\}\times dx)\hfill \\ & -\frac{1}{1+W^r{\left(v\right)}_t}\int (e^{{(v+zu)}^{\prime }x}-e^{v^{\prime }x}){\nu }^r(\left\{t\right\}\times dx))\hfill \\ \hfill & +\sum _{t<·}(\int (e^{{(v+zu)}^{\prime }x}-e^{v^{\prime }x}){\nu }^g(\left\{t\right\}\times dx)\hfill \\ & -\frac{1}{1+W^g{\left(v\right)}_t}\int (e^{{(v+zu)}^{\prime }x}-e^{v^{\prime }x}){\nu }^g(\left\{t\right\}\times dx)).\hfill \end{array}\end{array}$$

Now, using the above observations, we can show that the jump terms on the left-hand side and right-hand side of (18) are equal. Indeed, we can express the integrals with respect to the compensators ${\nu }^r$ and ${\nu }^g$ as sums in (20) and (21). Then we have that the jump term on the left-hand side of (18) is

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill & \sum _{t\le ·}\int (e^{v^{\prime }x}-1-v^{\prime }h\left(x\right)+\frac{e^{{(zu+v)}^{\prime }x}-e^{v^{\prime }x}}{1+W^r{\left(v\right)}_t}\hfill \\ \hfill & -z{u}^{\prime }h\left(x\right)+e^{{(v+zu)}^{\prime }x}-e^{v^{\prime }x}\hfill \\ \hfill & -\frac{1}{1+W^r{\left(v\right)}_t}\int (e^{{(v+zu)}^{\prime }x}-e^{v^{\prime }x})){\nu }^r(\left\{t\right\}\times dx)\hfill \\ \hfill & +\sum _{t\le ·}\int (e^{v^{\prime }x}-1-v^{\prime }h\left(x\right)+\frac{e^{{(zu+v)}^{\prime }x}-e^{v^{\prime }x}}{1+W^g{\left(v\right)}_t}\hfill \\ \hfill & -z{u}^{\prime }h\left(x\right)+e^{{(v+zu)}^{\prime }x}-e^{v^{\prime }x}\hfill \\ \hfill & -\frac{1}{1+W^g{\left(v\right)}_t}\int (e^{{(v+zu)}^{\prime }x}-e^{v^{\prime }x})){\nu }^g(\left\{t\right\}\times dx)\hfill \\ \hfill =& \sum _{t\le ·}\int \left(e^{{(v+zu)}^{\prime }x}-1-{(v+zu)}^{\prime }h\left(x\right)\right){\nu }^r(\left\{t\right\}\times dx)\hfill \\ \hfill & +\sum _{t<·}\int \left(e^{{(v+zu)}^{\prime }x}-1-{(v+zu)}^{\prime }h\left(x\right)\right){\nu }^g(\left\{t\right\}\times dx)\hfill \\ \hfill =& \sum _{t\le ·}\left(e^{{(v+zu)}^{\prime }x}-1-{(v+zu)}^{\prime }h\left(x\right)\right)*{\nu }^r\hfill \\ \hfill & +\left(e^{{(v+zu)}^{\prime }x}-1-{(v+zu)}^{\prime }h\left(x\right)\right)*{\nu }^g\hfill \end{array}\end{array}$$

which equals the corresponding quantity on the right-hand side of (18). The proof is finished. □

5. Duality relations

In this section, we establish an equivalent relation, namely the duality relation, between the prices of Margrabe options, and quanto options with European call and put options in a market characterized by an unusual probability space. This analysis takes into account a portfolio composed of optional processes.

To initiate, we aim to establish the duality relationship between the value of a swap option, specifically a Margrabe option or Spread option, with a payoff of $M_T={(S_T^1-S_T^2)}^{+}$, and the payoff of European call and put option. It is important to note that we represent the payoff of a European call option at maturity T by $C_T={(S_T-K)}^{+}$, where K is the strike price, and for the put option, it is denoted by $P_T={(K-S_T)}^{+}$.

Theorem 3.

Assume that the asset price processes $S^1$ and $S^2$ are exponential special optional semimartingales and $e^{H^i}\in {\mathcal{M}}^O\left(P\right),i=1,2.$ Then we can relate the value of a swap option and a European option via the following equality:

$$M_T=P_T^v(1,S_T^u)=C_T^{\theta }(S_T^{\zeta },1),$$

(30)

Or

$$\mathbf{E}{(S_T^1-S_T^2)}^{+}={\mathbf{E}}_v{(1-S_T^u)}^{+}={\mathbf{E}}_{\theta }{(S_T^{\zeta }-1)}^{+},$$

(31)

where the characteristics $(C^u,{\nu }^{ur},{\nu }^{ug})$ and $(C^{\zeta },{\nu }^{\zeta r},{\nu }^{\zeta g})$ of $H^u=\log S^u$ and $H^{\zeta }=\log S^{\zeta },$ respectively, are given by Theorem 2 for $v={(1,0)}^{\prime },u={(-1,1)}^{\prime },$ and $\theta ={(0,1)}^{\prime },\zeta ={(1,-1)}^{\prime }.$

Proof. 

We will use asset $S^1$ as the numeraire asset; if we use asset $S^2$ instead, then we get the duality relationship with a call option. The value of the swap option is

$$\begin{array}{cc}\hfill M_T=\mathbf{E}{(S_T^1-S_T^2)}^{+}=& \mathbf{E}\left[e^{H_T^1}{\left(1-\frac{S_T^2}{S_T^1}\right)}^{+}\right]\hfill \\ \hfill =& \mathbf{E}\left[e^{v^{\prime }{H}_T}{\left(1-\frac{S_T^2}{S_T^1}\right)}^{+}\right]\hfill \end{array}$$

(32)

where $v={(1,0)}^{\prime }.$ Moreover, $e^{v^{\prime }H}\in {\mathcal{M}}^O\left(P\right)$ by assumption. Define a new measure $P_v$ via the Radon-Nikodym derivative

$$\frac{d{P}_v}{dP}=e^{v^{\prime }{H}_T},$$

then the pricing problem (32) becomes

$$M_T=\mathbf{E}{(S_T^1-S_T^2)}^{+}={\mathbf{E}}_v\left[{\left(1-\frac{S_T^2}{S_T^1}\right)}^{+}\right]$$

where we define the process $S^u={\left(S_t^u\right)}_{0\le t\le T}$ via

$$S_t^u=\frac{S_t^2}{S_t^1}=\frac{e^{H_t^2}}{e^{H_t^1}}=e^{u^{\prime }{H}_t}=e^{H_t^u},0\le t\le T,$$

for $u={(-1,1)}^{\prime }.$ The triplet of local characteristics of the semimartingale $H^u$ is given by theorem 2 for $v={(1,0)}^{\prime }$ and $u={(-1,1)}^{\prime }.$ Now, applying Lemma 3.1 in Gasparyan (1988), we obtain that

$$e^{u^{\prime }H}\in {\mathcal{M}}^O\left(P_v\right)$$

since

$$e^{u^{\prime }H}{e}^{v^{\prime }H}=e^{H^2}\in {\mathcal{M}}^O\left(P\right)$$

Therefore, we conclude that

$$M_T=\mathbf{E}{(S_T^1-S_T^2)}^{+}={\mathbf{E}}_v{(1-S_T^u)}^{+}=P_T^v(1,S_T^u).$$

Using the same methodology, we derive the proof of equivalency for the Margrabe option and the call option as well. □

Additionally, we demonstrate another duality, this time between a quanto call option, with the payoff $Q_T^C=S_T^1{(S_T^2-K)}^{+}$, and a European call option with payoff as defined before.

Theorem 4.

Assume that the asset price processes $S^1$ and $S^2$ are exponential special optional semimartingales and $e^{H^i}\in {\mathcal{M}}^O\left(P\right),i=1,2.$ Then we can relate the value of a quanto call option and a European call option via the following duality:

$$Q_T^C=C_T^v(S_T^u,K)$$

(33)

Or

$$\mathbf{E}\left[S_T^1{(S_T^2-K)}^{+}\right]={\mathbf{E}}_v{(S_T^u-K)}^{+},$$

(34)

where the characteristics $(C^u,{\nu }^{ur},{\nu }^{ug})$ of $H^u=\log S^u$ are given by Theorem 2 for $v={(1,0)}^{\prime },u={(0,1)}^{\prime }.$ An analogous duality result relates the quanto put option and the European put option.

Proof. 

The value of the quanto call option is

$$\begin{array}{cc}\hfill Q_T^C=\mathbf{E}\left[S_T^1{(S_T^2-K)}^{+}\right]=& \mathbf{E}\left[e^{H_T^1}{(S_T^2-K)}^{+}\right]\hfill \\ \hfill =& \mathbf{E}\left[e^{v^{\prime }H}{(e^{H_T^2}-K)}^{+}\right]\hfill \\ \hfill =& {\mathbf{E}}_v\left[{(e^{H_T^u}-K)}^{+}\right],\hfill \end{array}$$

where $\frac{d{P}_v}{dP}=e^{v^{\prime }{H}_T}$ for $v={(1,0)}^{\prime }$ and $H^u=u^{\prime }H$ for $u={(0,1)}^{\prime }.$ Hence, the statement follows. □

Additionally, we delve into the Call-Put duality within this market, exploring two relations for these two contracts. The proof for one of them (a) requires utilizing theorem (2) to establish this duality, while the proof for the other relation (b) involves portfolio value and the Black-Scholes model.

Theorem 5.

Assume that the asset price process S is exponential special optional semimartingales and $e^H\in {\mathcal{M}}^O\left(P\right)$. Then we can relate the value of a European call option and a European put option via the following equivalencies:

$$a)P_T(K,S_T)=K{C}_T^v(S_T^u,1/K)$$

(35)

where the characteristics $(C^u,{\nu }^{ur},{\nu }^{ug})$ of $H^u=\log S^u$ is given by Theorem 2 for $v=1,u=-1$. Furthermore, we can relate the value of a standard put option with volatility σ to a standard call option via the following equality with volatility $-\sigma$, as follows:

$$b)P_T(K,S_T)=C_T(-S_T,-k),\mathrm{Or}{P}_T(K,S_T,\sigma )=C_T(-S_T,-k,-\sigma )$$

(36)

 

a. The value of a put option is

$$\begin{array}{cc}\hfill P_T(K,S_T)=\mathbf{E}\left[{(K-S_T)}^{+}\right]=\mathbf{E}& \left[e^{H_T}{(K{e}^{-H_T}-1)}^{+}\right]\hfill \\ \hfill =K\mathbf{E}& \left[e^{H_T}{(e^{-H_T}-1/K)}^{+}\right]\hfill \\ \hfill =K\mathbf{E}& \left[e^{v{H}_T}{(e^{-H_T}-1/K)}^{+}\right]\hfill \\ \hfill =K{\mathbf{E}}_v& \left[{(e^{-H_T}-1/K)}^{+}\right]\hfill \end{array}$$

where $\frac{d{p}_v}{dp}=e^{v{H}_T}$ for $v=1$ and $H^u=uH$ for $u=-1$. Thus, the statement (a) follows. □

 

b. To prove this duality, we employ a different approach than the previous, specifically utilizing the Black-Scholes model. In the following, we provide detailed proof of the pricing of a European call option and subsequently investigate the pricing of a European put option. We then establish the relationship between these two.

Assume that the market consists of two types of securities B and S and a portfolio $\pi =(\eta ,\xi )$ which is composed of the optional processes $\eta$ and $\eta .\xi$ is the volume of the reference asset B, while $\xi$ is the volume of the security S. Suppose $B_t>0$ and $S_t\ge 0$ for all $t\ge 0$ and write the ratio process $R_t=S_t/B_t$. Then, the discounted value of the portfolio is $Y_t={\eta }_t+{\xi }_t{R}_t$ which is a real-valued optional semimartingale that had right and left limits. Furthermore, we restrict the portfolio, $\pi$, to be self-financing that is $d{Y}_t=\xi d{R}_t$ where the interest rate is zero. Now let’s consider the augmented Balck-Scholes model with left and right jumps as described in Abdelghani and Melnikov (2020),

$$\begin{array}{cc}\hfill B_t=& B_0+{\int }_{]0,t]}r{B}_{s}ds\hfill \end{array}$$

(37)

$$\begin{array}{cc}\hfill S_t=& S_0+{\int }_{]0,t]}{S}_{s^{-}}(\mu ds+\sigma d{W}_s+ad{L}_s^r)+{\int }_{[0,t[}b{S}_{s}d{L}_{s^{+}}^g\hfill \end{array}$$

(38)

where $L_t^r=L_t-\lambda t$, $L_t^g=-{\overline{L}}_{t^{-}}+\gamma t$, and r, $\mu$, $\sigma$, a, and b are constants. W is diffusion term and L and $\overline{L}$ are independent Poisson with constant intensity $\lambda$ and $\gamma$, respectively. Let the initial wealth account for $B_0$ and the initial price be $S_0$. We can write S as $S_t=S_0\mathcal{E}\left(H\right)$, where $H_t=\mu t+\sigma W_t+a(L_t-\lambda t)+b(\gamma t-{\overline{L}}_{t^{-}})$, with $H_0=0$, and $B_t=B_0\exp \left(h_t\right)$ where $h_t=rt$. Furthermore, we are gonna use the stochastic exponential form mentioned in Abdelghani and Melnikov (2020) as follows:

$$\mathcal{E}{\left(H\right)}_t=\exp [H_t-\frac{1}{2}\langle H^c,H^c\rangle ]\prod _{0<s\le t}(1+\Delta S_s)e^{-\Delta S_s}\prod _{0\le s<t}(1+{\Delta }^{+}{S}_s)e^{-{\Delta }^{+}{S}_s}$$

(39)

In this case, the Ratio process is:

$$\begin{array}{c}\hfill R_t=S_0\exp \left\{\left(\mu -r-\frac{1}{2}({\sigma }^2-\lambda a^2-\gamma b^2)\right)t+\sigma W_t\right\}\times \\ \hfill \prod _{0<s\le t}\left[(1+a\Delta L_t)e^{-a\Delta L_t}\right]\prod _{0\le s<t}\left[(1-b{\Delta }^{+}{\overline{L}}_{s^{-}})e^{-b{\Delta }^{+}{\overline{L}}_{s^{-}}}\right]\end{array}$$

and is not a local optional martingale. So we need to find a suitable derivative to transfer this ratio under a martingale measure Q such that $Y_t$ is a local optional martingale with respect to Q (see Melnikov and Shiryaev (1996) and for further explanation, you can refer to Abdelghani and Melnikov (2020) chapter 9). In the following, we repeat some of the key concepts. So we want the derivative $Z=\mathcal{E}\left(N\right)$ for which we have to find a local martingale N such that $\Psi (h,H,N)$,

$$\begin{array}{cc}\hfill \Psi (h,H,N)=& N_t+H_t-h_t+{\langle h^c-N^c,h^c-H^c\rangle }_t\hfill \\ \hfill +& \sum _{0<s\le t}\frac{(\Delta h_s-\Delta N_s)(\Delta h_s-\Delta H_s)}{1+\Delta h_s}\hfill \\ \hfill +& \sum _{0<s\le t}\frac{({\Delta }^{+}{h}_s-{\Delta }^{+}{N}_s)({\Delta }^{+}{h}_s-{\Delta }^{+}{H}_s)}{1+{\Delta }^{+}{h}_s}\hfill \end{array}$$

is a local martingale. Thus, for the case of European options, chapter 9 of Abdelghani and Melnikov (2020) takes $N_t=\varsigma W_t+c(L_t-\lambda t)+d(\gamma t-{\overline{L}}_{t^{-}})$ which is an optional local martingale that will render Z as an optional scaling factor. By this choice of N we get the following:

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill \Psi (h,H,N)=& (\varsigma +\sigma )W_t+(a+c+ac)(L_t-\lambda t)\hfill \\ \hfill +& (b+d-bd)(\gamma t-{\overline{L}}_{t^{-}})\hfill \\ \hfill +& (\mu -r+\varsigma \sigma +2ac\lambda +2bd\gamma )t\hfill \end{array}\end{array}$$

(40)

is local martingale if

$$\mu -r+\varsigma \sigma +2ac\lambda +2bd\gamma =0$$

(41)

Thus we have to find $(\varsigma ,\sigma ,d)$ such that (41) is satisfied. We find the martingale measure Q by solving this equation which has infinitely many solutions that mean the market of Black-Scholes with left and right jumps is incomplete.

Now let’s turn our attention to the problem at hand, pricing a European call option. As outlined in chapter 9, section 6 of Abdelghani and Melnikov (2020) with the following choice of parameters for equation (41):

$$\left\{c=\frac{-a}{1+a},d=\frac{-b}{1-b},\varsigma =\frac{r-\mu +2\lambda \frac{a^2}{1+a}+2\gamma \frac{b^2}{1+b}}{\sigma }\right\}$$

(42)

leads to the normalized price of R under Q that is:

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill R_t=& R_0\mathcal{E}\left(\Psi (h,H,N)\right)\hfill \\ \hfill =& R_0\mathcal{E}\left((\varsigma +\sigma )W_t\right)\hfill \\ \hfill =& S_0\exp \{(\varsigma +\sigma )W_t-\frac{1}{2}{(\varsigma +\sigma )}^{2}t\},\hfill \\ \hfill B_0=& 1\hfill \end{array}\end{array}$$

(43)

It is just a function of the Wiener process and, all the left jumps are absorbed. Thus, $R_t$ is $Q-$martingale. Thus,

$$E_Q\left[R_T\right|F_t]=R_t=e^{-rt}{B}_t$$

(44)

As calculated in Abdelghani and Melnikov (2020), with this choice of parameters, the price of a call option at arbitrary time $0\le t<T$ is:

$$C_t=e^{rt}{R}_t\Phi \left({\tilde{K}}_t+(\varsigma +\sigma )\sqrt{T-t}\right)+e^{-r(T-t)}K\Phi \left({\tilde{K}}_t\right)$$

(45)

where:

$$\begin{array}{c}\hfill \tilde{K}=\frac{1}{(\varsigma +\sigma )\sqrt{T}}\left[\ln \left(\frac{S_0}{K}\right)+(r-\frac{1}{2}{(\varsigma +\sigma )}^2)T\right]\end{array}$$

(46)

Now, we proceed with the pricing of a European Put Option for optional processes. The payoff function of a European put option can be written as:

$$P_T={(K-S_T)}^{+}={(S_T-K)}^{+}+K-S_T$$

(47)

Using the payoff of the call option, we can write:

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill P_t=& E_Q\left[P_T\right|F_t]\hfill \\ \hfill =& E_Q\left[{(S_T-K)}^{+}\right|F_t]+E_Q[k-S_T|F_t]\hfill \\ \hfill =& C_t+E_Q[k-S_T|F_t]\hfill \end{array}\end{array}$$

(48)

Where $C_t$ is the price of the call option. Now if we denote:

$${\tilde{P}}_T=\frac{P_T}{B_T}$$

(49)

We have:

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill {\tilde{P}}_t=& E_Q\left[\frac{P_T}{B_T}\right|F_t]\hfill \\ \hfill =& E_Q\left[\frac{C_T}{B_T}\right|F_t]+E_Q\left[\frac{K-S_T}{B_T}\right|F_t]\hfill \\ \hfill =& {\tilde{C}}_t+K{e}^{-rT}-E_Q\left[R_T\right|F_t]\hfill \end{array}\end{array}$$

(50)

To compute $E_Q\left[R_T\right|F_t]$, we use the same martingale measure as it is used for pricing the European call option. Thus, by (44) We can rewrite (50) using (49) as follows:

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill {\tilde{P}}_t=& e^{-rt}{P}_t\hfill \\ \hfill =& e^{-rt}{C}_t+K{e}^{-rT}-e^{-rt}{S}_t\hfill \end{array}\end{array}$$

(51)

By simplifying we get the following relation which is called Call-Put Parity.

$$P_t=C_t+K{e}^{-r(T-t)}-S_t$$

(52)

By using this Call-Put Parity relation, we can calculate the fair price of the put option and also establish the duality relation. Using the fact that $\Phi (-x)=1-\Phi \left(x\right)$, we have:

$$\begin{array}{cc}\hfill P_0=& C_0+K{e}^{-rT}-S_0\hfill \\ \hfill =& S_0\Phi (\tilde{K}+(\varsigma +\sigma )\sqrt{T})-e^{-rT}K\Phi \left(\tilde{K}\right)+K{e}^{-rT}-S_0\hfill \\ \hfill =& S_0\left[\Phi (\tilde{K}+(\varsigma +\sigma )\sqrt{T})-1\right]-e^{-rT}(-K)(\Phi \left(\tilde{K}\right)-1)\hfill \\ \hfill =& (-S_0)\Phi (-\tilde{K}-(\varsigma +\sigma )\sqrt{T})-e^{-rT}(-K)\Phi (-\tilde{K})\hfill \end{array}$$

Therefore, using (46) definition, we have:

$$\begin{array}{c}\hfill -\tilde{K}=\frac{1}{-(\varsigma +\sigma )\sqrt{T}}\left[\ln \left(\frac{-S_0}{-K}\right)+(r-\frac{1}{2}{(-\varsigma -\sigma )}^2)T\right]\end{array}$$

(53)

Now given the definition of $\varsigma$ in (42), we can see that:

$$-\varsigma =\frac{r-\mu +2\lambda \frac{a^2}{1+a}+2\gamma \frac{b^2}{1+b}}{-\sigma }$$

is achieved simply by using $-\sigma$ thus we established the following relation which is called Call-Put Dulaity.

$$P_0(S_0,K,\sigma )=C_0(-S_0,-K,-\sigma )$$

(54)

And for arbitrary time $0\le t<T$, we have:

$$\begin{array}{cc}\hfill P_t=& C_t+K{e}^{-r(T-t)}-S_t\hfill \\ \hfill =& r^{rt}{R}_t\Phi ({\tilde{K}}_t+(\varsigma +\sigma )\sqrt{T-t})-e^{-r(T-t)}K\Phi \left({\tilde{K}}_t\right)+K{e}^{-r(T-t)}-S_t\hfill \\ \hfill =& (-S_t)\Phi (-{\tilde{K}}_t-(\varsigma +\sigma )\sqrt{T-t})-e^{-r(T-t)}(-K)\Phi (-{\tilde{K}}_t)\hfill \end{array}$$

where:

$$\begin{array}{c}\hfill {\tilde{K}}_t=\frac{1}{(\varsigma +\sigma )\sqrt{T-t}}\left[\ln \left(\frac{S_t}{K}\right)+(r-\frac{1}{2}{(\varsigma +\sigma )}^2)(T-t)\right]\end{array}$$

(55)

Thus, the Call-Put Duality for arbitrary time $0\le t<T$ is:

$$P_t(S_t,K,\sigma )=C_t(-S_t,-K,-\sigma )$$

(56)

Upon examining the duality relation in chapter 4 of the Melnikov (2011) and drawing a parallel, it becomes apparent that the (56) relation bears similarity to the duality relation involving cadlag processes. □

6. Applications and Examples

In this section, we proceed with concrete examples and explicit calculations for the previously obtained results in Section 5. The first example delves into an optional version of the well-known Merton’s jump-diffusion model and calculate the price of a quanto call option.

Example 1.

Let’s consider an optional version of the well-known Merton’s jump-diffusion model. The two assets $S^1$ and $S^2$ are modeled as optional exponential jump-diffusion processes:

$$S_t^i=\exp \left(b_{i}t+W_t^i+\sum _{k=0}^{N_t^{i,r}}{U}_k^i+\sum _{k=0}^{N_t^{i,g}}{Z}_k^i\right),i=1,2,$$

(57)

where the drift terms $b^i$ are determined by martingale condition; a Wiener process $W_t^i$ have a variance ${\sigma }_i\ge 0$ and a correlation coefficient $-1\le \rho \le 1,$ i.e. $\langle W^1,W^2\rangle =\rho$ and

$$c=\left(\begin{array}{cc}{\sigma }_1^2& \rho {\sigma }_1{\sigma }_2\\ \rho {\sigma }_1{\sigma }_2& {\sigma }_2^2\end{array}\right);$$

the process $N_t^{i,r},i=1,2$ is a right-continuous Poisson process with intensity ${\lambda }_i^r$ while the process $N_t^{i,r},i=1,2$ is a left-continuous version of a Poisson process with intensity ${\lambda }_i^g,$ the random variables $U_k^i\sim Normal(s_i^r,0)$ and $Z_k^i\sim Normal(s_i^g,0).$ The Poisson processes $N_t^{i,r},i=1,2$ and $N_t^{i,g},i=1,2$ are independent. Thus, the local characteristics of $H=(H^1,H^2)$ are

$$b=\left(\begin{array}{c}{b}_1\\ b_2\end{array}\right);c=\left(\begin{array}{cc}{\sigma }_1^2& \rho {\sigma }_1{\sigma }_2\\ \rho {\sigma }_1{\sigma }_2& {\sigma }_2^2\end{array}\right);$$

(58)

and $F^r$ and $F^g$ have the Lebesgue densities $f^r$ and $f^g$ where

$$f^j(x_1,x_2)=\frac{1}{s_1^j{s}_2^{j}2\pi }\exp \left(-\frac{x_1^2}{2{\left(s_1^j\right)}^2}-\frac{x_2^2}{2{\left(s_2^j\right)}^2}\right){\lambda }_1^j{\lambda }_2^j,j=r,g.$$

(59)

Now, the price of a quanto call option with strike price K, according to Theorem 4, is equal to the price of a call option with the same strike $K,$ on an asset $S^u$ with characteristics $(B^u,C^u,{\nu }^{ur},{\nu }^{ug})$ provided by Theorem 2 for $v={(1,0)}^{\prime }$ and $u={(0,1)}^{\prime }.$ Hence, we can calculate

$$\begin{array}{cc}\hfill b^u=& b_2+\rho {\sigma }_1{\sigma }_2,\hfill \end{array}$$

(60)

$$\begin{array}{cc}\hfill c^u=& {\sigma }_2^2;\hfill \end{array}$$

(61)

for the Levy measures, using the independence of the normal variables and completing the square, we have for $y\in \mathbb{R},E\in \mathcal{B}(\mathbb{R}\setminus 0):$

$$\begin{array}{cc}\hfill 1_E\left(y\right)*F^{uj}=& 1_E\left(x_2\right)e^{x_1}*F\hfill \\ \hfill =& {\int }_{{\mathbb{R}}^2}{1}_E\left(x_2\right)e^{x_1}\frac{1}{s_1^j{s}_2^{j}2\pi }\exp \left(-\frac{x_1^2}{2{\left(s_1^j\right)}^2}-\frac{x_2^2}{2{\left(s_2^j\right)}^2}\right){\lambda }_1^j{\lambda }_2^{j}d{x}_{1}d{x}_2\hfill \\ \hfill =& {\lambda }_2^j{\int }_E\frac{1}{s_2^j\sqrt{2\pi }}\exp \left(-\frac{x_2^2}{2{\left(s_2^j\right)}^2}\right)d{x}_2\hfill \\ \hfill & \times {\lambda }_1^j{e}^{{\left(s_1^j\right)}^2/2}{\int }_{\mathbb{R}}\frac{1}{s_1^j\sqrt{2\pi }}\exp \left(-\frac{{(x_1-{\left(s_1^j\right)}^2)}^2}{2{\left(s_1^j\right)}^2}\right)d{x}_1\hfill \\ \hfill =& {\lambda }_2^j{\lambda }_1^j{e}^{{\left(s_1^j\right)}^2/2}{\int }_E\frac{1}{s_2^j\sqrt{2\pi }}\exp \left(-\frac{x_2^2}{2{\left(s_2^j\right)}^2}\right)d{x}_2,j=r,g.\hfill \end{array}$$

(62)

Therefore, pricing a quanto option in case the two assets follow the optional version of Merton’s jump-diffusion model, is equivalent to pricing a call option in a univariate jump-diffusion model of the same class, with parameters given by (60)-(62). In particular, left jumps occur according to a compound Poisson process with intensity ${\lambda }^{ur}={\lambda }_2^r{\lambda }_1^r{e}^{{\left(s_1^r\right)}^2/2}$, and jump heights are normally distributed with jump variance $s_2^r$ and zero mean, while right jumps occur according to a left-continuous compound Poisson process with intensity ${\lambda }^{ug}={\lambda }_2^g{\lambda }_1^g{e}^{{\left(s_1^g\right)}^2/2},$ and jump heights are normally distributed with jump variance $s_2^g$ and zero mean.

Moreover, the second example extends the previous model by incorporating exponential jumps and proceeds to calculate the price of a quanto call option.

Example 2.

Let us now consider the previous example in a case when the random variables $U_k^i$ and $Z_k^i$ are mutually independent and exponentially distributed with parameters ${\alpha }^i$ and ${\beta }^i$ respectively. Thus, the density functions $f^r$ and $f^g$ have the following forms

$$\begin{array}{c}\hfill f^r(x_1,x_2)={\alpha }^1{\alpha }^2{e}^{-{\alpha }^1{x}_1-{\alpha }^2{x}_2},\\ \hfill f^g(x_1,x_2)={\beta }^1{\beta }^2{e}^{-{\beta }^1{x}_1-{\beta }^2{x}_2}.\end{array}$$

Now, the price of a quanto call option with strike $K,$ according to Theorem 4, is equal to the price of a call option with the same strike $K,$ on an asset $S^u$ with characteristics $(B^u,C^u,{\nu }^{ur},{\nu }^{ug})$ provided by Theorem 2 for $v={(1,0)}^{\prime }$ and $u={(0,1)}^{\prime }.$ Consequently, in this case $b^u$ and $c^u$ are the same as in Example 1, and we obtain for $y\in \mathbb{R},E\in \mathcal{B}(\mathbb{R}\setminus \{0\left\}\right):$

$$\begin{array}{cc}\hfill 1_E\left(y\right)*F^{ur}=& {\lambda }_1^r{\lambda }_2^r{\int }_{{\mathbb{R}}^2}{1}_E\left(x_2\right)e^{x_1}{\alpha }_1{\alpha }_2{e}^{-{\alpha }_1{x}_1-{\alpha }_2{x}_{2}d{x}_{1}d{x}_2}\hfill \\ \hfill =& {\lambda }_1^r{\lambda }_2^r{\alpha }_1{\alpha }_2{\int }_{\mathbb{R}}{e}^{x_1(1-{\alpha }_1)}d{x}_1{\int }_E{e}^{-{\alpha }_2{x}_2}d{x}_2\hfill \\ \hfill =& {\lambda }_1^r{\lambda }_2^r\frac{{\alpha }_1{\alpha }_2}{{\alpha }_1-1}{\int }_E{e}^{-{\alpha }_2{x}_2}d{x}_2.\hfill \end{array}$$

The levy measure for left jumps is calculated in a similar way. As a result, pricing a quanto option in case of two assets are modeled as optional jump-diffusion’s with exponential jumps, is equivalent to pricing a call option in a univariate optional jump-diffusion model of the same class. In particular, left jumps occur according to a compound Poisson process with intensity ${\lambda }^{ur}={\lambda }_2^r{\lambda }_1^r\frac{{\alpha }_1}{{\alpha }_1-1}$, and jump heights are exponentially distributed with a parameter ${\alpha }_2$, while right jumps occur according to a left-continuous compound Poisson process with intensity ${\lambda }^{ug}={\lambda }_2^g{\lambda }_1^g\frac{{\beta }_1}{{\beta }_1-1},$ and jump heights are exponentially distributed with a parameter ${\beta }_2$.