Key research themes
1. How can realized volatility and high-frequency data improve estimation in stochastic volatility models?
This research focuses on leveraging high-frequency intraday data to construct realized volatility measures, facilitating more precise estimation and inference in stochastic volatility (SV) models. Realized volatility, computed as sums of squared intraday returns, offers consistent, although biased in finite samples, estimates of integrated volatility. The central inquiry is how to quantify and correct for errors between realized volatility and the latent integrated volatility in SV models, improving estimation without resorting to simulation-heavy methods.
2. What are the theoretical and empirical advancements in continuous-time Gaussian and self-similar stochastic volatility models, especially regarding option pricing and small-time asymptotics?
This theme investigates continuous-time Gaussian SV models with self-similar volatility processes, focusing on their small-time asymptotic behavior, option price and implied volatility sensitivities, and the mathematical characterization of volatility dynamics via fractional and Gaussian processes. The role of self-similarity parameter (Hurst index) in controlling the short maturity implied volatility surface and providing model-based estimators is emphasized.
3. How do high-dimensional and multivariate stochastic volatility models address computational challenges and improve covariance modeling in financial time series?
Investigates methodological advancements to scale stochastic volatility models to high-dimensional multivariate settings while overcoming computational bottlenecks of Bayesian MCMC and Monte Carlo likelihood estimation. Focus is on penalized ordinary least squares (OLS) frameworks incorporating sparsity-inducing penalties, factor structures, and state-space representations to enable efficient estimation and forecasting of large covariance matrices relevant to portfolio risk and asset dependence structures.


















![As mentioned before, there exist many different formulas, for example, by Kahl and Jackel [28], Lewis[29] or Zhylyevskyy [36]. We will use the approach by Lewis [29] because it is well suited for more complex models with jumps. It is also well behaved compared with the formula by Albrecher et al. [27] but has the numerical advantage that we only have to calculate one integral for each call option price. where X = In(S/K) + rr and](https://figures.academia-assets.com/118628658/figure_001.jpg)






![New pricing formula for the Yan Hanson model is (31) with otk) defined in (23), or (33) in terms of @(k). 4.2. Example: a new fractional SVJD model We want to modify the Bates model from Section 3 by using an approximate fractional Brownian motion in the volatility process. This process has a long memory for Hurst parameter H > 0.5, and for H = 0.5, it turns into a standard Wiener process. Thao [39] defined approximative fractional process as an It6 integral,](https://figures.academia-assets.com/118628658/figure_014.jpg)


![where (W,),59 is a standard Wiener process, p < 0, A, vy 2 0 and (Z,),50 is a Lévy subordinator (independent on W,)!!. The model price of a call option with strike K can be retrieved using a standard Fourier transform [17]. Option price (A3) can be rewritten using the notation C(u) = C(u, tT) = exp{A(u, tT) + B(u, r)}. Also, note that C(—i) = 1 for any positive tT.](https://figures.academia-assets.com/118628658/figure_019.jpg)






![where we used the Fubini’s theorem (the integrand is measurable and integrable) and the fact that ie g(y)dy = 1. W substitute t = T — t and define h(k, v, t) by to obtain from (25) the following equation for h: which is equal to equation (2.7) on p. 38 in Lewis [29] and has a fundamental solution F with initial value F (k,v,0) = 1 (in Lewis [29], it is called fundamental transform). It is regular as a function of k = k, + ik; within a strip k,; < k; < kp. From this fundamental solution, we can derive the explicit formula for the option price](https://figures.academia-assets.com/118628658/figure_013.jpg)

![and @(k) is defined in (32) with @(k) as in Example |. Similarly, we could use the formula (31). Equation (41) is a Ricatti equation and can be solved explicitly (for example, PospiSil and Sobotka [26], Proposition 2.1), and then we obtain C; by integrating (42). The formula (33) for this model with approximate fractional Brownian motion is](https://figures.academia-assets.com/118628658/figure_015.jpg)



![Extending the work of Cerny [1] to the case where the price of a risky asset price is represented by a lognormal geometric brownian motion with constant parameters, as in (2.2), we obtain explicit expressions for both the mean-value process, H(r), of the option to be hedged, and the optimal control, u(r), to be applied at the rebalancing instant 7. The main results are given by Theorems 2.1 and 2.2 stated below. Full proofs can be found in Maiali [6]. In what follows we use the following notation:](https://figures.academia-assets.com/105324473/figure_001.jpg)





























