Source author record

Piergiacomo Sabino

Piergiacomo Sabino appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

math.PR q-fin.CP q-fin.PR q-fin.MF

Catalog footprint

What is connected

7works

4topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

Efficient Simulation of $p$-Tempered $α$-Stable OU Processes

We develop efficient methods for simulating processes of Ornstein-Uhlenbeck type related to the class of $p$-tempered $α$-stable ($\ts$) distributions. Our results hold for both the univariate and multivariate cases and we consider both the case where the $\ts$ distribution is the stationary law and where it is the distribution of the background driving Lévy process (BDLP). In the latter case, we also derive an explicit representation for the transition law as this was previous known only in certain special cases and only for $p=1$ and $α\in[0,1)$. Simulation results suggest that our methods work well in practice.

preprint2022arXiv

Exchange option pricing under variance gamma-like models

In this article we focus on the pricing of exchange options when the dynamic of logprices follows either the well-known variance gamma or the recent variance gamma++ process introduced in Gardini et al [19]. In particular, for the former model we can derive a Margrabe's type formula whereas, for the latter one we can write an "integral free" formula. Furthermore, we show how to construct a general multidimensional versions of the variance gamma++ processes preserving both the mathematical and numerical tractability. Finally we apply the derived models to German and French energy power markets: we calibrate their parameters using real market data and we accordingly evaluate exchange options with the derived closed formulas, Fourier based methods and Monte Carlo techniques.

preprint2020arXiv

Correlating Lévy processes with Self-Decomposability: Applications to Energy Markets

Based on the concept of self-decomposability, we extend some recent multivariate Lévy models built using multivariate subordination with the aim of capturing situations in which a sudden event in one market is propagated onto related markets after a certain stochastic time delay. Consequently, we study the properties of such processes, derive closed form expressions for the characteristic function and detail how a Monte Carlo scheme can be easily implemented. We illustrate the applicability of our approach in the context of gas and power Energy markets focusing on the calibration and on the pricing of spread options written on different underlying assets using simulations techniques.

preprint2016arXiv

Cointegrating Jumps: an Application to Energy Facilities

Based on the concept of self-decomposable random variables we discuss the application of a model for a pair of dependent Poisson processes to energy facilities. Due to the resulting structure of the jump events we can see the self-decomposability as a form of cointegration among jumps. In the context of energy facilities, the application of our approach to model power or gas dynamics and to evaluate transportation assets seen as spread options is straightforward. We study the applicability of our methodology first assuming a Merton market model with two underlying assets; in a second step we consider price dynamics driven by an exponential mean-reverting Geometric Ornstein-Uhlenbeck plus compound Poisson that are commonly used in the energy field. In this specific case we propose a price spot dynamics for each underlying that has the advantage of being treatable to find non-arbitrage conditions. In particular we can find close-form formulas for vanilla options so that the price and the Greeks of spread options can be calculated in close form using the Margrabe formula (if the strike is zero) or some other well known approximation.

preprint2011arXiv

Multidimensional Quasi-Monte Carlo Malliavin Greeks

We investigate the use of Malliavin calculus in order to calculate the Greeks of multidimensional complex path-dependent options by simulation. For this purpose, we extend the formulas employed by Montero and Kohatsu-Higa to the multidimensional case. The multidimensional setting shows the convenience of the Malliavin Calculus approach over different techniques that have been previously proposed. Indeed, these techniques may be computationally expensive and do not provide flexibility for variance reduction. In contrast, the Malliavin approach exhibits a higher flexibility by providing a class of functions that return the same expected value (the Greek) with different accuracies. This versatility for variance reduction is not possible without the use of the generalized integral by part formula of Malliavin Calculus. In the multidimensional context, we find convenient formulas that permit to improve the localization technique, introduced in Fournié et al and reduce both the computational cost and the variance. Moreover, we show that the parameters employed for variance reduction can be obtained \textit{on the flight} in the simulation. We illustrate the efficiency of the proposed procedures, coupled with the enhanced version of Quasi-Monte Carlo simulations as discussed in Sabino, for the numerical estimation of the Deltas of call, digital Asian-style and Exotic basket options with a fixed and a floating strike price in a multidimensional Black-Scholes market.

preprint2010arXiv

Convenient Multiple Directions of Stratification

This paper investigates the use of multiple directions of stratification as a variance reduction technique for Monte Carlo simulations of path-dependent options driven by Gaussian vectors. The precision of the method depends on the choice of the directions of stratification and the allocation rule within each strata. Several choices have been proposed but, even if they provide variance reduction, their implementation is computationally intensive and not applicable to realistic payoffs, in particular not to Asian options with barrier. Moreover, all these previously published methods employ orthogonal directions for multiple stratification. In this work we investigate the use of algorithms producing convenient directions, generally non-orthogonal, combining a lower computational cost with a comparable variance reduction. In addition, we study the accuracy of optimal allocation in terms of variance reduction compared to the Latin Hypercube Sampling. We consider the directions obtained by the Linear Transformation and the Principal Component Analysis. We introduce a new procedure based on the Linear Approximation of the explained variance of the payoff using the law of total variance. In addition, we exhibit a novel algorithm that permits to correctly generate normal vectors stratified along non-orthogonal directions. Finally, we illustrate the efficiency of these algorithms in the computation of the price of different path-dependent options with and without barriers in the Black-Scholes and in the Cox-Ingersoll-Ross markets.

preprint2009arXiv

Pricing and Hedging Asian Basket Options with Quasi-Monte Carlo Simulations

In this article we consider the problem of pricing and hedging high-dimensional Asian basket options by Quasi-Monte Carlo simulation. We assume a Black-Scholes market with time-dependent volatilities and show how to compute the deltas by the aid of the Malliavin Calculus, extending the procedure employed by Montero and Kohatsu-Higa (2003). Efficient path-generation algorithms, such as Linear Transformation and Principal Component Analysis, exhibit a high computational cost in a market with time-dependent volatilities. We present a new and fast Cholesky algorithm for block matrices that makes the Linear Transformation even more convenient. Moreover, we propose a new-path generation technique based on a Kronecker Product Approximation. This construction returns the same accuracy of the Linear Transformation used for the computation of the deltas and the prices in the case of correlated asset returns while requiring a lower computational time. All these techniques can be easily employed for stochastic volatility models based on the mixture of multi-dimensional dynamics introduced by Brigo et al. (2004).