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

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

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

A generic construction for high order approximation schemes of semigroups using random grids

Our aim is to construct high order approximation schemes for general semigroups of linear operators $P_{t},t\geq 0$. In order to do it, we fix a time horizon $T$ and the discretization steps $h_{l}=\frac{T}{n^{l}},l\in \mathbb{N}$ and we suppose that we have at hand some short time approximation operators $Q_{l}$ such that $P_{h_{l}}=Q_{l}+O(h_{l}^{1+\alpha })$ for some $\alpha >0$. Then, we consider random time grids $\Pi (\omega )=\{t_0(\omega )=0<t_{1}(\omega )<...<t_{m}(\omega )=T\}$ such that for all $1\le k\le m$, $t_{k}(\omega )-t_{k-1}(\omega )=h_{l_{k}}$ for some $l_{k}\in \mathbb{N}$, and we associate the approximation discrete semigroup $P_{T}^{\Pi (\omega )}=Q_{l_{n}}...Q_{l_{1}}.$ Our main result is the following: for any approximation order $\nu$, we can construct random grids $\Pi_{i}(\omega )$ and coefficients $c_{i}$, with $i=1,...,r$ such that $$P_{t}f=\sum_{i=1}^{r}c_{i}\mathbb{E}(P_{t}^{\Pi _{i}(\omega )}f(x))+O(n^{-\nu})$$% with the expectation concerning the random grids $\Pi _{i}(\omega ).$ Besides, $\text{Card}(\Pi _{i}(\omega ))=O(n)$ and the complexity of the algorithm is of order $n$, for any order of approximation $\nu$. The standard example concerns diffusion processes, using the Euler approximation for~$Q_l$. In this particular case and under suitable conditions, we are able to gather the terms in order to produce an estimator of $P_tf$ with finite variance. However, an important feature of our approach is its universality in the sense that it works for every general semigroup $P_{t}$ and approximations. Besides, approximation schemes sharing the same $\alpha$ lead to the same random grids $\Pi_{i}$ and coefficients $c_{i}$. Numerical illustrations are given for ordinary differential equations, piecewise deterministic Markov processes and diffusions

Optimal execution and absence of price manipulations in limit order book models

We continue the analysis of optimal execution strategies in the model for a limit order book with nonlinear price impact and exponential resilience that was considered in our earlier paper with A. Fruth. We now allow for non-homogeneous resilience rates and arbitrary trading dates and consider the extended problem of optimizing jointly over trading dates and sizes. Our main results show that, under general conditions on the shape function of the limit order book, placing the deterministic trade sizes identified in our earlier paper at trading dates that are homogeneously spaced is optimal also within the large class of adaptive strategies with arbitrary trading dates. Perhaps even more importantly, our analysis yields as a corollary that our model does not admit price manipulation strategies in the sense of Huberman and Stanzl. This latter result contrasts the recent findings of Gatheral, where, in a related but different model, exponential resilience was found to give rise to price manipulation strategies when price impact is nonlinear.

Optimal trade execution and absence of price manipulations in limit order book models

We analyze the existence of price manipulation and optimal trade execution strategies in a model for an electronic limit order book with nonlinear price impact and exponential resilience. Our main results show that, under general conditions on the shape function of the limit order book, placing deterministic trade sizes at trading dates that are homogeneously spaced is optimal within a large class of adaptive strategies with arbitrary trading dates. This extends results from our earlier work with A. Fruth. Perhaps even more importantly, our analysis yields as a corollary that our model does not admit price manipulation strategies. This latter result contrasts the recent findings of Gatheral [12], where, in a related but different model, exponential resilience was found to give rise to price manipulation strategies when price impact is nonlinear.

Nonnegativity preserving convolution kernels. Application to Stochastic Volterra Equations in closed convex domains and their approximation

This work defines and studies convolution kernels that preserve nonnegativity. When the past dynamics of a process is integrated with a convolution kernel like in Stochastic Volterra Equations or in the jump intensity of Hawkes processes, this property allows to get the nonnegativity of the integral. We give characterizations of these kernels and show in particular that completely monotone kernels preserve nonnegativity. We then apply these results to analyze the stochastic invariance of a closed convex set by Stochastic Volterra Equations. We also get a comparison result in dimension one. Last, when the kernel is a positive linear combination of decaying exponential functions, we present a second order approximation scheme for the weak error that stays in the closed convex domain under suitable assumptions. We apply these results to the rough Heston model and give numerical illustrations