---
title: Weak Variance Alpha-Gamma Process
url: https://www.emergentmind.com/topics/weak-variance-alpha-gamma-process
type: topic
---

# Weak Variance Alpha-Gamma Process

{"query":"Weak variance alpha-gamma process arXiv WVAG weak subordination selfdecomposability calibration", "max_results": 10}
The weak variance alpha-gamma process, usually abbreviated **WVAG** and also called the **weak variance-\(\alpha\)-gamma** or **weak VaG** process, is a multivariate Lévy process obtained by weakly subordinating an \(n\)-dimensional Brownian motion by an \(n\)-dimensional alpha-gamma subordinator. In the cited literature, it is presented as a multivariate generalization of the variance gamma framework that relaxes the classical requirement of independent Brownian components, thereby allowing an arbitrary covariance matrix \(\Sigma\) while retaining the Lévy property [1609.04481]. It is also a principal example of the broader class of weak variance generalized gamma convolution processes, \(WVGG^n\), obtained by weak subordination with multivariate Thorin subordinators [1712.03640].

## 1. Definition and origin

Weak subordination was introduced by Buchmann, Lu, and Madan as an extension of both univariate and multivariate classical subordination. In the classical multivariate setting, the pathwise operation \(X\circ T\) generally preserves the Lévy property only in two special situations: when the components of \(T\) are indistinguishable, or when the components of \(X\) are independent. Weak subordination replaces literal pathwise composition by a law-based construction \(X\odot T\) that still yields a Lévy process even when \(X\) has dependent components [1609.04481].

Within that framework, let \(n\ge 2\), let \(B\sim BM^n(\mu,\Sigma)\) be an \(n\)-dimensional Brownian motion with drift \(\mu\in\mathbb{R}^n\) and arbitrary covariance matrix \(\Sigma\in\mathbb{R}^{n\times n}\), and let \(T\sim AG_S^n(a,\alpha)\) be an \(n\)-dimensional alpha-gamma subordinator. Then
\[
X \sim WVAG^n(a,\mu,\Sigma,\alpha)
\]
is called a weak variance alpha-gamma process provided
\[
X = B \odot T.
\]
Another paper writes the parameter ordering as \(WVAG^n(a,\mu,\alpha,\Sigma)\); both notational conventions occur in the cited literature [1712.03640].

The alpha-gamma subordinator is built from independent gamma subordinators \(G_0,\dots,G_n\). With \(a>0\), \(\alpha=(\alpha_1,\dots,\alpha_n)\in(0,\infty)^n\), and \(a\alpha_k<1\), define
\[
\beta_k=\frac{1-a\alpha_k}{\alpha_k},
\]
take
\[
G_0\sim \Gamma_S(a,1), \qquad G_k\sim \Gamma_S(\beta_k,1/\alpha_k),\quad 1\le k\le n,
\]
independently, and set
\[
T = (T_1,\dots,T_n) \overset{d}{=} G_0\,\alpha + (G_1,\dots,G_n).
\]
This gives each component a **common gamma factor** \(G_0\) together with an **idiosyncratic gamma factor** \(G_k\) [1712.03640].

The direct predecessor is the strong variance-\(\alpha\)-gamma process of Luciano and Semeraro, obtained by traditional subordination with an alpha-gamma clock and a Brownian motion with independent components. The weak model removes the restriction that \(\Sigma\) be diagonal. When \(\Sigma\) is diagonal, the weak model reduces to the classical strong \(VAG\) model [1801.08852].

## 2. Characteristic exponent, decomposition, and dependence structure

For \(X\sim WVAG^n(a,\mu,\alpha,\Sigma)\), the cited calibration paper gives the characteristic exponent
\[
\Psi_X(\theta) = -a\ln\left\{ 1-i\langle \mu,\theta\rangle+\frac12\|\theta\|_\Sigma^2 \right\}
-\sum_{k=1}^n \beta_k \ln\left\{ 1-i\alpha_k\mu_k\theta_k+\frac12\alpha_k\Sigma_{kk}\theta_k^2 \right\}.
\]
This formula exhibits a decomposition into a common multivariate variance-gamma-type term and coordinatewise corrections [1801.08852].

The same paper states that \(X\) has the decomposition in law
\[
X \stackrel{D}{=} V_0+\sum_{k=1}^n V_k e_k,
\]
where \(V_0\) is multivariate \(VG\) and the \(V_k\) are univariate \(VG\) terms. In that sense, the model is a sum of one common multivariate variance-gamma component and several idiosyncratic univariate variance-gamma components. A later paper on energy spread options makes this representation explicit in a closely related notation, writing the WVAG process as a sum of independent \(VG\) processes plus deterministic drift [2507.11480].

Each marginal is univariate \(VG\):
\[
X_k\sim VG^1(1/\alpha_k,\mu_k,\Sigma_{kk}), \qquad k=1,\dots,n.
\]
Thus the model preserves variance-gamma marginals while allowing multivariate dependence beyond the strong \(VAG\) setup [1801.08852].

A central feature is the covariance structure. For \(k\neq \ell\),
\[
\operatorname{Cov}(X_k(1),X_\ell(1))
=
a(\alpha_k\wedge \alpha_\ell)\Sigma_{k\ell}
+
a\alpha_k\alpha_\ell \mu_k\mu_\ell.
\]
The cited paper emphasizes that the term \(a(\alpha_k\wedge \alpha_\ell)\Sigma_{k\ell}\) gives \(WVAG\) a broader dependence range than \(VAG\). Dependence can enter through Brownian correlation \(\Sigma_{k\ell}\), drift interaction \(\mu_k\mu_\ell\), and the shared gamma-clock effect via \(a\) and the \(\alpha_k\) [1801.08852].

The model is therefore characterized in the literature by **common and idiosyncratic time changes**, **VG marginals with possibly different kurtoses**, **full support jump measure**, and a **wider range of dependence** than \(VAG\) [1609.04481].

## 3. Thorin representation and the \(WVGG^n\) framework

The weak variance alpha-gamma process is a special case of a weak variance generalized gamma convolution process. In the general \(WVGG^n\) construction, if \(B\sim BM^n(\mu,\Sigma)\) and \(T\sim GGC_S^n(d,\mathcal{T})\), then
\[
X=B\odot T\sim WVGG^n(d,\mu,\Sigma,\mathcal{T}).
\]
For the alpha-gamma case, the Thorin measure is particularly simple [1712.03640].

Specifically, the Thorin measure of a driftless alpha-gamma subordinator is
\[
\mathcal{T} = a\,\delta_{\,\alpha/\|\alpha\|^2} +\sum_{k=1}^n \beta_k\,\delta_{\,e_k/\alpha_k}.
\]
Thus the weak variance alpha-gamma process is a weakly subordinated Brownian motion by a Thorin subordinator whose Thorin measure is supported on finitely many rays or points [1712.03640].

In the general \(WVGG^n\) setting, the characteristic exponent is
\[
\Psi_X(\theta) = i\langle \mu,\theta\rangle -\frac12 \|\theta\|_\Sigma^2 -\int_{[0,\infty)^n_*} \ln\!\left( \frac{\|\xi\|^2 - i\langle \mu,\xi\rangle + \tfrac12\|\xi\|_\Sigma^2}{\|\xi\|^2} \right)\,\mathcal{T}(d\xi), \qquad \theta\in\mathbb{R}^n,
\]
and in the driftless case \(\mu=0\) this simplifies to
\[
\Psi_X(\theta) = -\frac12\|\theta\|_\Sigma^2 -\int_{[0,\infty)^n_*} \ln\!\left( \frac{\|\xi\|^2+\tfrac12\|\theta\|_\Sigma^2}{\|\xi\|^2} \right)\,\mathcal{T}(d\xi).
\]
For \(WVAG^n\), this specialization becomes a finite linear combination of logarithmic terms corresponding to the common factor \(G_0\) and the idiosyncratic gamma components \(G_k\) [1712.03640].

The Lévy-measure interpretation developed in the weak-subordination paper is consistent with this structure. Jumps occur in two ways: a **common jump** driven by \(G_0\) affecting all components together, and **idiosyncratic jumps** driven by \(G_k\) affecting component \(k\) individually. This decomposition is the hallmark of the weak VaG construction [1609.04481].

## 4. Self-decomposability

A major structural question for \(WVAG\) is self-decomposability. The general result for \(WVGG^n\) states that if the Brownian motion subordinate is driftless, then the process is self-decomposable:
\[
\mu = 0 \quad \Longrightarrow \quad X\sim SD^n.
\]
More precisely, for \(n\ge 2\),
\[
X\sim WVGG^n_{\mathbf 0} \quad \Longrightarrow \quad X\sim SD^n.
\]
This extends an earlier strong-subordination result to weak subordination [1712.03640].

For the weak variance alpha-gamma process, the specialization is especially sharp. Corollary \(\ref{corWValphaG}\) in the cited paper states that if
\[
X\sim WVAG^n(a,\mu,\Sigma,\alpha),
\]
then
\[
\mu=0 \quad \Longrightarrow \quad X\sim SD^n,
\]
while if
\[
|\Sigma|\neq 0 \quad \text{and}\quad \mu\neq 0,
\]
then
\[
X\not\sim SD^n.
\]
Accordingly, for the weak variance alpha-gamma process, self-decomposability is completely characterized by driftlessness, provided the Brownian covariance is invertible [1712.03640].

The negative direction in the general \(WVGG^n\) theory depends on moment conditions for the Thorin measure. A principal sufficient condition is
\[
0<\int_{(0,\infty)^n}\left(1+\|\xi\|^{1/2}\right)
\left(\frac{\|\xi\|^n}{\prod_{k=1}^n \xi_k}\right)^{1/2}\,\mathcal{T}(d\xi) <\infty.
\]
If \(n\ge 2\), \(|\Sigma|\neq 0\), \(\mu\neq 0\), and this condition holds, then
\[
X\not\sim SD^n.
\]
For measures of the form
\[
\mathcal{T}=\sum_{k=1}^m \int_{(0,\infty)} \delta_{v_k u}\,\mathcal{T}_k(du),
\]
the key integrability condition is equivalent to
\[
\int_{(1,\infty)} u^{1/2}\,\mathcal{T}_k(du)<\infty, \qquad 1\le k\le m.
\]
In the alpha-gamma case, because the Thorin measure is discrete on rays, these conditions become transparent and can be checked explicitly [1712.03640].

The cited paper also emphasizes an important subtlety. The moment conditions are sufficient for non-self-decomposability, but not always necessary in a naive form. An example is constructed of a \(WVGG^n\) process that is still self-decomposable even though the Brownian motion subordinate has nonzero drift. This example is not the alpha-gamma process itself. For \(WVAG\), however, the result is complete under invertibility of \(\Sigma\): driftlessness is exactly the criterion [1712.03640].

## 5. Calibration and statistical inference

A dedicated calibration study compares three estimation methods for the \(WVAG\) parameters \((a,\mu,\alpha,\Sigma)\): **method of moments (MOM)**, **maximum likelihood estimation (MLE)**, and **digital moment estimation (DME)**. MOM is based on least-squares matching of theoretical and empirical moments. MLE uses a numerically computed likelihood via Fourier inversion. DME fits the model by matching empirical and theoretical quantiles or probabilities [1801.08852].

Because the density is not known in closed form, MLE requires Fourier inversion of the characteristic function:
\[
f_{Y(t)}(x) = (2\pi)^{-n} \int_{\mathbb{R}^n} e^{-i\langle x,\theta\rangle}\Phi_{X(t)}(\theta)\,d\theta, \qquad \Phi_{X(t)}(\theta)=\exp\big(t\Psi_X(\theta)\big),
\]
provided \(\Phi_{X(t)}\in L^1\). For \(X\sim WVAG^n(a,\mu,\alpha,\Sigma)\) and \(Y=I+X\), assuming \(\Sigma\) is invertible, the cited paper derives the sufficient condition
\[
\left(\frac{a}{n}+\min_{1\le k\le n}\beta_k\right)t>\frac12,\qquad t>0,
\]
which implies
\[
\Phi_{X(t)},\Phi_{Y(t)}\in L^1.
\]
This is the Fourier invertibility condition used to justify numerical likelihood evaluation [1801.08852].

The same paper reports two simulation regimes. For a bivariate \(WVAG\) model with \(c=1\), the condition is satisfied and **MLE gives the best overall fit**. For \(c=0.1\), the condition is violated and **DME gives the best fit**, while **MLE still works reasonably well** and MOM again performs worst [1801.08852].

In an empirical application to daily log-returns of the S&P500 and FTSE100 over five years, the paper finds that the **WVAG model fits better than VAG** across fit diagnostics, that for WVAG **DME is best** among the three estimation methods, and that the likelihood-ratio test strongly rejects the VAG restriction \(\Sigma_{12}=0\):
\[
D=514.03,\qquad p<10^{-4}.
\]
The estimated drift vector is reported to be very close to zero, and the authors note that, using the self-decomposability result for \(WVAG\), they do not reject self-decomposability at the \(5\%\) level [1801.08852].

## 6. Ornstein–Uhlenbeck embeddings, Esscher transform, and applications

The \(WVAG\) process has also been embedded into multivariate Lévy-driven Ornstein–Uhlenbeck dynamics in two distinct ways. In **WVAG-OU**, the stationary distribution of the OU process is \(WVAG\). In **OU-WVAG**, the background driving Lévy process (BDLP) itself is \(WVAG\). For the OU recursion at equally spaced times \(t_k=k\Delta\),
\[
X(t_k)=e^{-\lambda \Delta}X(t_{k-1}) + e^{-\lambda\Delta} Z^*(\Delta)^{(k)},
\]
where
\[
Z^*(\Delta) := \int_0^{\lambda\Delta} e^s\,dZ(s).
\]
This innovation term is the basis of the likelihood and simulation theory in both models [2011.14542].

For the **WVAG-OU** model, the cited paper derives an explicit BDLP representation and shows that the BDLP is a compound Poisson process with drift. The corresponding innovation law is a discrete-continuous mixture, which yields exact simulation and a modified likelihood based on a dominating mixture measure. For the **OU-WVAG** model, the paper proves that the innovation term \(Z^*(\Delta)\) is absolutely continuous, so the standard Lebesgue likelihood applies, although the exponent lacks the closed form available in the \(WVAG\)-OU case [2011.14542].

A later application uses \(WVAG\) as the BDLP in an OU model for energy prices, with
\[
\log S(t)=\boldsymbol{\Lambda}(t)+X(t).
\]
That paper emphasizes that \(WVAG\) is a flexible multivariate jump model with a natural dependence structure, and studies forward pricing, Carr–Madan call pricing, and Hurd–Zhou FFT pricing for spread options [2507.11480].

A central structural result in that setting is that \(WVAG\) is **not closed under the Esscher transform**. Although the transformed process remains a sum of independent \(VG\) components, it is generally not a \(WVAG\) process of the original form, because the transformed marginal components are no longer jointly parameterized in the WVAG structure. The paper identifies this non-closure as a key difference between \(WVAG\) and the multivariate \(VG\) class [2507.11480].

The literature situates \(WVAG\) in several finance-facing contexts. The self-decomposability paper notes applications in **instantaneous portfolio theory**, **multivariate stock return modeling**, and **calibration work in finance**, and stresses that weak subordination is attractive because it allows general covariance matrices \(\Sigma\), unlike strong subordination models that often require diagonality or independent components [1712.03640]. The energy paper extends that line to **energy spread options** under variance gamma-driven OU dynamics [2507.11480].

A related but indirect development is the study of gamma-driven stochastic differential equations of the form
\[
dX_t=\sigma(X_{t-})\,dL_t,\qquad X_0=0,
\]
where \(L\) is a gamma process. That work does **not** introduce a process called the weak variance alpha-gamma process, but it is described as a mathematically relevant analogue for gamma-driven, monotone, volatility-modulated models with weak solution theory and likelihood-ratio formulas [2108.11891]. A plausible implication is that the broader analytical toolkit around gamma subordinators, absolute continuity of laws, and likelihood-based inference continues to reinforce the mathematical setting in which \(WVAG\) is studied.

Source: https://www.emergentmind.com/topics/weak-variance-alpha-gamma-process