Source author record

Aurélien Alfonsi

Aurélien Alfonsi appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

18works
11topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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.

Published work

18 published item(s)

preprint2022arXiv

Approximation of Stochastic Volterra Equations with kernels of completely monotone type

In this work, we develop a multifactor approximation for $d$-dimensional Stochastic Volterra Equations (SVE) with Lipschitz coefficients and kernels of completely monotone type that may be singular. First, we prove an $L^2$-estimation between two SVEs with different kernels, which provides a quantification of the error between the SVE and any multifactor Stochastic Differential Equation (SDE) approximation. For the particular rough kernel case with Hurst parameter lying in $(0,1/2)$, we propose various approximating multifactor kernels, state their rates of convergence and illustrate their efficiency for the rough Bergomi model. Second, we study a Euler discretization of the multifactor SDE and establish a convergence result towards the SVE that is uniform with respect to the approximating multifactor kernels. These obtained results lead us to build a new multifactor Euler scheme that reduces significantly the computational cost in an asymptotic way compared to the Euler scheme for SVEs. Finally, we show that our multifactor Euler scheme outperforms the Euler scheme for SVEs for option pricing in the rough Heston model.

preprint2021arXiv

Constrained overdamped Langevin dynamics for symmetric multimarginal optimal transportation

The Strictly Correlated Electrons (SCE) limit of the Levy-Lieb functional in Density Functional Theory (DFT) gives rise to a symmetric multi-marginal optimal transport problem with Coulomb cost, where the number of marginal laws is equal to the number of electrons in the system, which can be very large in relevant applications. In this work, we design a numerical method, built upon constrained overdamped Langevin processes to solve Moment Constrained Optimal Transport (MCOT) relaxations (introduced in A. Alfonsi, R. Coyaud, V. Ehrlacher and D. Lombardi, Math. Comp. 90, 2021, 689--737) of symmetric multi-marginal optimal transport problems with Coulomb cost. Some minimizers of such relaxations can be written as discrete measures charging a low number of points belonging to a space whose dimension, in the symmetrical case, scales linearly with the number of marginal laws. We leverage the sparsity of those minimizers in the design of the numerical method and prove that any local minimizer to the resulting problem is actually a \emph{global} one. We illustrate the performance of the proposed method by numerical examples which solves MCOT relaxations of 3D systems with up to 100 electrons.

preprint2016arXiv

Maximum Likelihood Estimation for Wishart processes

In the last decade, there has been a growing interest to use Wishart processes for modelling, especially for financial applications. However, there are still few studies on the estimation of its parameters. Here, we study the Maximum Likelihood Estimator (MLE) in order to estimate the drift parameters of a Wishart process. We obtain precise convergence rates and limits for this estimator in the ergodic case and in some nonergodic cases. We check that the MLE achieves the optimal convergence rate in each case. Motivated by this study, we also present new results on the Laplace transform that extend the recent findings of Gnoatto and Grasselli and are of independent interest.

preprint2015arXiv

Dynamic optimal execution in a mixed-market-impact Hawkes price model

We study a linear price impact model including other liquidity takers, whose flow of orders either follows a Poisson or a Hawkes process. The optimal execution problem is solved explicitly in this context, and the closed-formula optimal strategy describes in particular how one should react to the orders of other traders. This result enables us to discuss the viability of the market. It is shown that Poissonian arrivals of orders lead to quite robust Price Manipulation Strategies in the sense of Huberman and Stanzl. Instead, a particular set of conditions on the Hawkes model balances the self-excitation of the order flow with the resilience of the price, excludes Price Manipulation Strategies and gives some market stability.

preprint2015arXiv

Extension and calibration of a Hawkes-based optimal execution model

We provide some theoretical extensions and a calibration protocol for our former dynamic optimal execution model. The Hawkes parameters and the propagator are estimated independently on financial data from stocks of the CAC40. Interestingly, the propagator exhibits a smoothly decaying form with one or two dominant time scales, but only so after a few seconds that the market needs to adjust after a large trade. Motivated by our estimation results, we derive the optimal execution strategy for a multi-exponential Hawkes kernel and backtest it on the data for round trips. We find that the strategy is profitable on average when trading at the midprice, which is in accordance with violated martingale conditions. However, in most cases, these profits vanish when we take bid-ask costs into account.

preprint2015arXiv

Multivariate transient price impact and matrix-valued positive definite functions

We consider a model for linear transient price impact for multiple assets that takes cross-asset impact into account. Our main goal is to single out properties that need to be imposed on the decay kernel so that the model admits well-behaved optimal trade execution strategies. We first show that the existence of such strategies is guaranteed by assuming that the decay kernel corresponds to a matrix-valued positive definite function. An example illustrates, however, that positive definiteness alone does not guarantee that optimal strategies are well-behaved. Building on previous results from the one-dimensional case, we investigate a class of nonincreasing, nonnegative and convex decay kernels with values in the symmetric $K\times K$ matrices. We show that these decay kernels are always positive definite and characterize when they are even strictly positive definite, a result that may be of independent interest. Optimal strategies for kernels from this class are well-behaved when one requires that the decay kernel is also commuting. We show how such decay kernels can be constructed by means of matrix functions and provide a number of examples. In particular we completely solve the case of matrix exponential decay.

preprint2015arXiv

Optimal transport bounds between the time-marginals of a multidimensional diffusion and its Euler scheme

In this paper, we prove that the time supremum of the Wasserstein distance between the time-marginals of a uniformly elliptic multidimensional diffusion with coefficients bounded together with their derivatives up to the order $2$ in the spatial variables and H{ö}lder continuous with exponent $γ$ with respect to the time variable and its Euler scheme with $N$ uniform time-steps is smaller than $C \left(1+\mathbf{1}\_{γ=1} \sqrt{\ln(N)}\right)N^{-γ}$. To do so, we use the theory of optimal transport. More precisely, we investigate how to apply the theory by Ambrosio, Gigli and Savar{é} to compute the time derivative of the Wasserstein distance between the time-marginals. We deduce a stability inequality for the Wasserstein distance which finally leads to the desired estimation.

preprint2015arXiv

Smile with the Gaussian term structure model

We propose an affine extension of the Linear Gaussian term structure Model (LGM) such that the instantaneous covariation of the factors is given by an affine process on semidefinite positive matrices. First, we set up the model and present some important properties concerning the Laplace transform of the factors and the ergodicity of the model. Then, we present two main numerical tools to implement the model in practice. First, we obtain an expansion of caplets and swaptions prices around the LGM. Such a fast and accurate approximation is useful for assessing the model behavior on the implied volatility smile. Second, we provide a second order scheme for the weak error, which enables to calculate exotic options by a Monte-Carlo algorithm.

preprint2013arXiv

A remark on the optimal transport between two probability measures sharing the same copula

We are interested in the Wasserstein distance between two probability measures on $\R^n$ sharing the same copula $C$. The image of the probability measure $dC$ by the vectors of pseudo-inverses of marginal distributions is a natural generalization of the coupling known to be optimal in dimension $n=1$. It turns out that for cost functions $c(x,y)$ equal to the $p$-th power of the $L^q$ norm of $x-y$ in $\R^n$, this coupling is optimal only when $p=q$ i.e. when $c(x,y)$ may be decomposed as the sum of coordinate-wise costs.

preprint2013arXiv

Capacitary measures for completely monotone kernels via singular control

We give a singular control approach to the problem of minimizing an energy functional for measures with given total mass on a compact real interval, when energy is defined in terms of a completely monotone kernel. This problem occurs both in potential theory and when looking for optimal financial order execution strategies under transient price impact. In our setup, measures or order execution strategies are interpreted as singular controls, and the capacitary measure is the unique optimal control. The minimal energy, or equivalently the capacity of the underlying interval, is characterized by means of a nonstandard infinite-dimensional Riccati differential equation, which is analyzed in some detail. We then show that the capacitary measure has two Dirac components at the endpoints of the interval and a continuous Lebesgue density in between. This density can be obtained as the solution of a certain Volterra integral equation of the second kind.

preprint2013arXiv

Exact and high order discretization schemes for Wishart processes and their affine extensions

This work deals with the simulation of Wishart processes and affine diffusions on positive semidefinite matrices. To do so, we focus on the splitting of the infinitesimal generator, in order to use composition techniques as Ninomiya and Victoir or Alfonsi. Doing so, we have found a remarkable splitting for Wishart processes that enables us to sample exactly Wishart distributions, without any restriction on the parameters. It is related but extends existing exact simulation methods based on Bartlett's decomposition. Moreover, we can construct high-order discretization schemes for Wishart processes and second-order schemes for general affine diffusions. These schemes are in practice faster than the exact simulation to sample entire paths. Numerical results on their convergence are given.

preprint2012arXiv

A Mean-Reverting SDE on Correlation matrices

We introduce a mean-reverting SDE whose solution is naturally defined on the space of correlation matrices. This SDE can be seen as an extension of the well-known Wright-Fisher diffusion. We provide conditions that ensure weak and strong uniqueness of the SDE, and describe its ergodic limit. We also shed light on a useful connection with Wishart processes that makes understand how we get the full SDE. Then, we focus on the simulation of this diffusion and present discretization schemes that achieve a second-order weak convergence. Last, we explain how these correlation processes could be used to model the dependence between financial assets.

preprint2012arXiv

Optimal execution and price manipulations in time-varying limit order books

This paper focuses on an extension of the Limit Order Book (LOB) model with general shape introduced by Alfonsi, Fruth and Schied. Here, the additional feature allows a time-varying LOB depth. We solve the optimal execution problem in this framework for both discrete and continuous time strategies. This gives in particular sufficient conditions to exclude Price Manipulations in the sense of Huberman and Stanzl or Transaction-Triggered Price Manipulations (see Alfonsi, Schied and Slynko). These conditions give interesting qualitative insights on how market makers may create or not price manipulations.

preprint2012arXiv

Strong convergence of some drift implicit Euler scheme. Application to the CIR process

We study the convergence of a drift implicit scheme for one-dimensional SDEs that was considered by Alfonsi for the Cox-Ingersoll-Ross (CIR) process. Under general conditions, we obtain a strong convergence of order 1. In the CIR case, Dereich, Neuenkirch and Szpruch have shown recently a strong convergence of order 1/2 for this scheme. Here, we obtain a strong convergence of order 1 under more restrictive assumptions on the CIR parameters.

preprint2010arXiv

Optimal execution strategies in limit order books with general shape functions

We consider optimal execution strategies for block market orders placed in a limit order book (LOB). We build on the resilience model proposed by Obizhaeva and Wang (2005) but allow for a general shape of the LOB defined via a given density function. Thus, we can allow for empirically observed LOB shapes and obtain a nonlinear price impact of market orders. We distinguish two possibilities for modeling the resilience of the LOB after a large market order: the exponential recovery of the number of limit orders, i.e., of the volume of the LOB, or the exponential recovery of the bid-ask spread. We consider both of these resilience modes and, in each case, derive explicit optimal execution strategies in discrete time. Applying our results to a block-shaped LOB, we obtain a new closed-form representation for the optimal strategy, which explicitly solves the recursive scheme given in Obizhaeva and Wang (2005). We also provide some evidence for the robustness of optimal strategies with respect to the choice of the shape function and the resilience-type.

preprint2006arXiv

A Call-Put Duality for Perpetual American Options

It is well known that in models with time-homogeneous local volatility functions and constant interest and dividend rates, the European Put prices are transformed into European Call prices by the simultaneous exchanges of the interest and dividend rates and of the strike and spot price of the underlying. This paper investigates such a Call Put duality for perpetual American options. It turns out that the perpetual American Put price is equal to the perpetual American Call price in a model where, in addition to the previous exchanges between the spot price and the strike and between the interest and dividend rates, the local volatility function is modified. We prove that equality of the dual volatility functions only holds in the standard Black-Scholes model with constant volatility. Thanks to these duality results, we design a theoretical calibration procedure of the local volatility function from the perpetual Call and Put prices for a fixed spot price $x_0$. The knowledge of the Put (resp. Call) prices for all strikes enables to recover the local volatility function on the interval $(0,x_0)$ (resp. $(x_0,+\infty)$).

preprint2006arXiv

General Duality for Perpetual American Options

In this paper, we investigate the generalization of the Call-Put duality equality obtained in [1] for perpetual American options when the Call-Put payoff $(y-x)^+$ is replaced by $ϕ(x,y)$. It turns out that the duality still holds under monotonicity and concavity assumptions on $ϕ$. The specific analytical form of the Call-Put payoff only makes calculations easier but is not crucial unlike in the derivation of the Call-Put duality equality for European options. Last, we give some examples for which the optimal strategy is known explicitly.