BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Benjamin Jourdain

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

ResearcherUnclaimed source record
math.PRq-fin.CPmath.APq-fin.PRmath.NAmath.OCmath.STMethodologyq-fin.MFq-fin.PMStatistics Theory

Catalog footprint

What is connected

35works
11topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

35 published item(s)

preprint2026arXiv

Linear independence properties of the signature components of time-augmented stochastic processes

The addition of the running time as a component of a path before computing its signature is a widespread approach to ensure the one-to-one property between them and leads to universal approximation theorems (Cuchiero, Primavera and Svaluto-Ferro, 2023). However, this also leads to the linear dependence of the components of the terminal value of the signature of the time-augmented path. More precisely, for a given natural number $N$, the signature components associated with words of length $N$ have the same linear span as the signature components associated with words of length not greater than $N$. We generalize this result by exhibiting other subfamilies of signature components with the same spanning properties. In particular we recover the result of Dupire and Tissot-Daguette which states that the spanning of the iterated integrals with the last integrator different from the time variable is the same as the spanning of all iterated integrals. We check that this choice leads to the minimal computation time when the terms of the signature are calculated using Chen's relation in a backward way. The same optimal computation time is symmetrically achieved in a forward way for the iterated integrals with the first integrator different from the time variable. Building on these results, we derive several results regarding the linear independence of the signature components of a time-augmented stochastic process. We show that if the stochastic process we consider is solution to some SDE with additive Brownian noise then any subfamily of components proposed previously is linearly independent. We also prove that the linear independence of these subfamilies of components is still true when we consider the discretization of the sample paths of this stochastic process on a grid with a sufficiently small discretization time step.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Central limit theorem over non-linear functionals of empirical measures: beyond the iid setting

The central limit theorem is, with the strong law of large numbers, one of the two fundamental limit theorems in probability theory. Benjamin Jourdain and Alvin Tse have extended to non-linear functionals of the empirical measure of independent and identically distributed random vectors the central limit theorem which is well known for linear functionals. The main tool permitting this extension is the linear functional derivative, one of the notions of derivation on the Wasserstein space of probability measures that have recently been developed. The purpose of this work is to generalize what has been done by Jourdain and Tse: provide a Central Limit Theorem for non-linear functionals of independent and non equidistributed random vectors such as the successive values of an ergodic Markov chain.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Lipschitz continuity of the Wasserstein projections in the convex order on the line

Wasserstein projections in the convex order were first considered in the framework of weak optimal transport, and found application in various problems such as concentration inequalities and martingale optimal transport. In dimension one, it is well-known that the set of probability measures with a given mean is a lattice w.r.t. the convex order. Our main result is that, contrary to the minimum and maximum in the convex order, the Wasserstein projections are Lipschitz continuity w.r.t. the Wasserstein distance in dimension one. Moreover, we provide examples that show sharpness of the obtained bounds for the 1-Wasserstein distance.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Stability of the Weak Martingale Optimal Transport Problem

While many questions in (robust) finance can be posed in the martingale optimal transport (MOT) framework, others require to consider also non-linear cost functionals. Following the terminology of Gozlan, Roberto, Samson and Tetali this corresponds to weak martingale optimal transport (WMOT). In this article we establish stability of WMOT which is important since financial data can give only imprecise information on the underlying marginals. As application, we deduce the stability of the superreplication bound for VIX futures as well as the stability of stretched Brownian motion and we derive a monotonicity principle for WMOT.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

One dimensional martingale rearrangement couplings

We are interested in martingale rearrangement couplings. As introduced by Wiesel [37] in order to prove the stability of Martingale Optimal Transport problems, these are projections in adapted Wasserstein distance of couplings between two probability measures on the real line in the convex order onto the set of martingale couplings between these two marginals. In reason of the lack of relative compactness of the set of couplings with given marginals for the adapted Wasserstein topology, the existence of such a projection is not clear at all. Under a barycentre dispersion assumption on the original coupling which is in particular satisfied by the Hoeffding-Fréchet or comonotone coupling, Wiesel gives a clear algorithmic construction of a martingale rearrangement when the marginals are finitely supported and then gets rid of the finite support assumption by relying on a rather messy limiting procedure to overcome the lack of relative compactness. Here, we give a direct general construction of a martingale rearrangement coupling under the barycentre dispersion assumption. This martingale rearrangement is obtained from the original coupling by an approach similar to the construction we gave in [24] of the inverse transform martingale coupling, a member of a family of martingale couplings close to the Hoeffding-Fréchet coupling, but for a slightly different injection in the set of extended couplings introduced by Beiglböck and Juillet [9] and which involve the uniform distribution on [0, 1] in addition to the two marginals. We last discuss the stability in adapted Wassertein distance of the inverse transform martingale coupling with respect to the marginal distributions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Strong solutions to a beta-Wishart particle system

The purpose of this paper is to study the existence and uniqueness of solutions to a Stochastic Differential Equation (SDE) coming from the eigenvalues of Wishart processes. The coordinates are non-negative, evolve as Cox-Ingersoll-Ross (CIR) processes and repulse each other according to a Coulombian like interaction force. We show the existence of strong and pathwise unique solutions to the system until the first multiple collision, and give a necessary and sufficient condition on the parameters of the SDE for this multiple collision not to occur in finite time.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Asymptotic error distribution for the Ninomiya-Victoir scheme in the commutative case

In a previous work, we proved strong convergence with order $1$ of the Ninomiya-Victoir scheme $X^{NV}$ with time step $T/N$ to the solution $X$ of the limiting SDE when the Brownian vector fields commute. In this paper, we prove that the normalized error process $N \left(X - X^{NV}\right)$ converges to an affine SDE with source terms involving the Lie brackets between the Brownian vector fields and the drift vector field. This result ensures that the strong convergence rate is actually $1$ when the Brownian vector fields commute, but at least one of them does not commute with the drift vector field. When all the vector fields commute the limit vanishes. Our result is consistent with the fact that the Ninomiya-Victoir scheme solves the SDE in this case.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Asymptotics for the normalized error of the Ninomiya-Victoir scheme

In a previous work, we proved strong convergence with order $1/2$ of the Ninomiya-Victoir scheme $X^{NV,η}$ with time step $T/N$ to the solution $X$ of the limiting SDE. In this paper we check that the normalized error defined by $\sqrt{N}\left(X - X^{NV,η}\right)$ converges to an affine SDE with source terms involving the Lie brackets between the Brownian vector fields. The limit does not depend on the Rademacher random variables $η$. This result can be seen as a first step to adapt to the Ninomiya-Victoir scheme the central limit theorem of Lindeberg Feller type, derived by M. Ben Alaya and A. Kebaier for the multilevel Monte Carlo estimator based on the Euler scheme. When the Brownian vector fields commute, the limit vanishes. This suggests that the rate of convergence is greater than $1/2$ in this case and we actually prove strong convergence with order $1$.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Ninomiya-Victoir scheme : strong convergence properties and discretization of the involved Ordinary Differential Equations

In this paper, we summarize the results about the strong convergence rate of the Ninomiya-Victoir scheme and the stable convergence in law of its normalized error that we obtained in previous papers. We then recall the properties of the multilevel Monte Carlo estimators involving this scheme that we introduced and studied before. Last, we are interested in the error introduced by discretizing the ordinary differential equations involved in the Ninomiya-Victoir scheme. We prove that this error converges with strong order 2 when an explicit Runge-Kutta method with order 4 (resp. 2) is used for the ODEs corresponding to the Brownian (resp. Stratonovich drift) vector fields. We thus relax the order 5 for the Brownian ODEs needed by Ninomiya and Ninomiya (2009) to obtain the same order of strong convergence. Moreover, the properties of our multilevel Monte-Carlo estimators are preserved when these Runge-Kutta methods are used.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A multitype sticky particle construction of Wasserstein stable semigroups solving one-dimensional diagonal hyperbolic systems with large monotonic data

This article is dedicated to the study of diagonal hyperbolic systems in one space dimension, with cumulative distribution functions, or more generally nonconstant monotonic bounded functions, as initial data. Under a uniform strict hyperbolicity assumption on the characteristic fields, we construct a multitype version of the sticky particle dynamics and obtain existence of global weak solutions by compactness. We then derive a $L^p$ stability estimate on the particle system uniform in the number of particles. This allows to construct nonlinear semigroups solving the system in the sense of Bianchini and Bressan [Ann. of Math. (2), 2005]. We also obtain that these semigroup solutions satisfy a stability estimate in Wasserstein distances of all orders, which encompasses the classical $L^1$ estimate and generalises to diagonal systems the results by Bolley, Brenier and Loeper [J. Hyperbolic Differ. Equ., 2005] in the scalar case. Our results are obtained without any smallness assumption on the variation of the data, and only require the characteristic fields to be Lipschitz continuous and the system to be uniformly strictly hyperbolic.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A trajectorial interpretation of the dissipations of entropy and Fisher information for stochastic differential equations

The dissipation of general convex entropies for continuous time Markov processes can be described in terms of backward martingales with respect to the tail filtration. The relative entropy is the expected value of a backward submartingale. In the case of (non necessarily reversible) Markov diffusion processes, we use Girsanov theory to explicit the Doob-Meyer decomposition of this submartingale. We deduce a stochastic analogue of the well known entropy dissipation formula, which is valid for general convex entropies, including the total variation distance. Under additional regularity assumptions, and using Itô's calculus and ideas of Arnold, Carlen and Ju \cite{Arnoldcarlenju}, we obtain moreover a new Bakry Emery criterion which ensures exponential convergence of the entropy to $0$. This criterion is non-intrisic since it depends on the square root of the diffusion matrix, and cannot be written only in terms of the diffusion matrix itself. We provide examples where the classic Bakry Emery criterion fails, but our non-intrisic criterion applies without modifying the law of the diffusion process.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Computation of sensitivities for the invariant measure of a parameter dependent diffusion

We consider the solution to a stochastic differential equation with a drift function which depends smoothly on some real parameter $λ$, and admitting a unique invariant measure for any value of $λ$ around $λ$ = 0. Our aim is to compute the derivative with respect to $λ$ of averages with respect to the invariant measure, at $λ$ = 0. We analyze a numerical method which consists in simulating the process at $λ$ = 0 together with its derivative with respect to $λ$ on long time horizon. We give sufficient conditions implying uniform-in-time square integrability of this derivative. This allows in particular to compute efficiently the derivative with respect to $λ$ of the mean of an observable through Monte Carlo simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Ninomiya-Victoir scheme: strong convergence, antithetic version and application to multilevel estimators

In this paper, we are interested in the strong convergence properties of the Ninomiya-Victoir scheme which is known to exhibit weak convergence with order 2. We prove strong convergence with order $1/2$. This study is aimed at analysing the use of this scheme either at each level or only at the finest level of a multilevel Monte Carlo estimator: indeed, the variance of a multilevel Monte Carlo estimator is related to the strong error between the two schemes used on the coarse and fine grids at each level. Recently, Giles and Szpruch proposed a scheme permitting to construct a multilevel Monte Carlo estimator achieving the optimal complexity $O\left(ε^{-2}\right)$ for the precision $ε$. In the same spirit, we propose a modified Ninomiya-Victoir scheme, which may be strongly coupled with order $1$ to the Giles-Szpruch scheme at the finest level of a multilevel Monte Carlo estimator. Numerical experiments show that this choice improves the efficiency, since the order $2$ of weak convergence of the Ninomiya-Victoir scheme permits to reduce the number of discretization levels.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On the long time behavior of stochastic vortices systems

In this paper, we are interested in the long-time behaviour of stochastic systems of n interacting vortices: the position in R2 of each vortex evolves according to a Brownian motion and a drift summing the influences of the other vortices computed through the Biot and Savart kernel and multiplied by their respective vorticities. For fixed n, we perform the rescalings of time and space used successfully by Gallay and Wayne [5] to study the long-time behaviour of the vorticity formulation of the two dimensional incompressible Navier-Stokes equation, which is the limit as n $\rightarrow$ $\infty$ of the weighted empirical measure of the system under mean-field interaction. When all the vorticities share the same sign, the 2n-dimensional process of the rescaled positions of the vortices is shown to converge exponentially fast as time goes to infinity to some invariant measure which turns out to be Gaussian if all the vorticities are equal. In the particular case n = 2 of two vortices, we prove exponential convergence in law of the 4-dimensional process to an explicit random variable, whatever the choice of the two vorticities. We show that this limit law is not Gaussian when the two vorticities are not equal.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Optimal convergence rate of the multitype sticky particle approximation of one-dimensional diagonal hyperbolic systems with monotonic initial data

Brenier and Grenier [SIAM J. Numer. Anal., 1998] proved that sticky particle dynamics with a large number of particles allow to approximate the entropy solution to scalar one-dimensional conservation laws with monotonic initial data. In [arXiv:1501.01498], we introduced a multitype version of this dynamics and proved that the associated empirical cumulative distribution functions converge to the viscosity solution, in the sense of Bianchini and Bres-san [Ann. of Math. (2), 2005], of one-dimensional diagonal hyperbolic systems with monotonic initial data of arbitrary finite variation. In the present paper, we analyse the L 1 error of this approximation procedure, by splitting it into the discretisation error of the initial data and the non-entropicity error induced by the evolution of the particle system. We prove that the error at time t is bounded from above by a term of order (1 + t)/n, where n denotes the number of particles, and give an example showing that this rate is optimal. We last analyse the additional error introduced when replacing the multitype sticky particle dynamics by an iterative scheme based on the typewise sticky particle dynamics, and illustrate the convergence of this scheme by numerical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Optimal scaling for the transient phase of the random walk Metropolis algorithm: The mean-field limit

We consider the random walk Metropolis algorithm on $\mathbb{R}^n$ with Gaussian proposals, and when the target probability measure is the $n$-fold product of a one-dimensional law. In the limit $n\to\infty$, it is well known (see [Ann. Appl. Probab. 7 (1997) 110-120]) that, when the variance of the proposal scales inversely proportional to the dimension $n$ whereas time is accelerated by the factor $n$, a diffusive limit is obtained for each component of the Markov chain if this chain starts at equilibrium. This paper extends this result when the initial distribution is not the target probability measure. Remarking that the interaction between the components of the chain due to the common acceptance/rejection of the proposed moves is of mean-field type, we obtain a propagation of chaos result under the same scaling as in the stationary case. This proves that, in terms of the dimension $n$, the same scaling holds for the transient phase of the Metropolis-Hastings algorithm as near stationarity. The diffusive and mean-field limit of each component is a diffusion process nonlinear in the sense of McKean. This opens the route to new investigations of the optimal choice for the variance of the proposal distribution in order to accelerate convergence to equilibrium (see [Optimal scaling for the transient phase of Metropolis-Hastings algorithms: The longtime behavior Bernoulli (2014) To appear]).

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Reducing the debt : is it optimal to outsource an investment?

We deal with the problem of outsourcing the debt for a big investment, according two situations: either the firm outsources both the investment (and the associated debt) and the exploitation to a private consortium, or the firm supports the debt and the investment but outsources the exploitation. We prove the existence of Stackelberg and Nash equilibria between the firm and the private consortium, in both situations. We compare the benefits of these contracts. We conclude with a study of what happens in case of incomplete information, in the sense that the risk aversion coefficient of each partner may be unknown by the other partner.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Stochastic particle approximation of the Keller-Segel equation and two-dimensional generalization of Bessel processes

The Keller-Segel partial differential equation is a two-dimensional model for chemotaxis. When the total mass of the initial density is one, it is known to exhibit blow-up in finite time as soon as the sensitivity $χ$ of bacteria to the chemo-attractant is larger than $8π$. We investigate its approximation by a system of $N$ two-dimensional Brownian particles interacting through a singular attractive kernel in the drift term. In the very subcritical case $χ\textless{}2π$, the diffusion strongly dominates this singular drift: we obtain existence for the particle system and prove that its flow of empirical measures converges, as $N\to\infty$ and up to extraction of a subsequence, to a weak solution of the Keller-Segel equation. We also show that for any $N\ge 2$ and any value of $χ\textgreater{}0$, pairs of particles do collide with positive probability: the singularity of the drift is indeed visited. Nevertheless, when $χ\textless{}2πN$, it is possible to control the drift and obtain existence of the particle system until the first time when at least three particles collide. We check that this time is a.s. infinite, so that global existence holds for the particle system, if and only if $χ\leq 8π(N-2)/(N-1)$. Finally, we remark that in the system with $N=2$ particles, the difference between the two positions provides a natural two-dimensional generalization of Bessel processes, which we study in details.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Capital distribution and portfolio performance in the mean-field Atlas model

We study a mean-field version of rank-based models of equity markets such as the Atlas model introduced by Fernholz in the framework of Stochastic Portfolio Theory. We obtain an asymptotic description of the market when the number of companies grows to infinity. Then, we discuss the long-term capital distribution. We recover the Pareto-like shape of capital distribution curves usually derived from empirical studies, and provide a new description of the phase transition phenomenon observed by Chatterjee and Pal. Finally, we address the performance of simple portfolio rules and highlight the influence of the volatility structure on the growth of portfolios.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Optimal scaling for the transient phase of Metropolis Hastings algorithms: The longtime behavior

We consider the Random Walk Metropolis algorithm on $\mathbb{R}^n$ with Gaussian proposals, and when the target probability measure is the $n$-fold product of a one-dimensional law. It is well known (see Roberts et al. (Ann. Appl. Probab. 7 (1997) 110-120)) that, in the limit $n\to\infty$, starting at equilibrium and for an appropriate scaling of the variance and of the timescale as a function of the dimension $n$, a diffusive limit is obtained for each component of the Markov chain. In Jourdain et al. (Optimal scaling for the transient phase of the random walk Metropolis algorithm: The mean-field limit (2012) Preprint), we generalize this result when the initial distribution is not the target probability measure. The obtained diffusive limit is the solution to a stochastic differential equation nonlinear in the sense of McKean. In the present paper, we prove convergence to equilibrium for this equation. We discuss practical counterparts in order to optimize the variance of the proposal distribution to accelerate convergence to equilibrium. Our analysis confirms the interest of the constant acceptance rate strategy (with acceptance rate between $1/4$ and $1/3$) first suggested in Roberts et al. (Ann. Appl. Probab. 7 (1997) 110-120). We also address scaling of the Metropolis-Adjusted Langevin Algorithm. When starting at equilibrium, a diffusive limit for an optimal scaling of the variance is obtained in Roberts and Rosenthal (J. R. Stat. Soc. Ser. B. Stat. Methodol. 60 (1998) 255-268). In the transient case, we obtain formally that the optimal variance scales very differently in $n$ depending on the sign of a moment of the distribution, which vanishes at equilibrium. This suggest that it is difficult to derive practical recommendations for MALA from such asymptotic results.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The small noise limit of order-based diffusion processes

We introduce order-based diffusion processes as the solutions to multidimensional stochastic differential equations, with drift coefficient depending only on the ordering of the coordinates of the process and diffusion matrix proportional to the identity. These processes describe the evolution of a system of Brownian particles moving on the real line with piecewise constant drifts, and are the natural generalization of the rank-based diffusion processes introduced in stochastic portfolio theory or in the probabilistic interpretation of nonlinear evolution equations. Owing to the discontinuity of the drift coefficient, the corresponding ordinary differential equations are ill-posed. Therefore, the small noise limit of order-based diffusion processes is not covered by the classical Freidlin-Wentzell theory. The description of this limit is the purpose of this article. We first give a complete analysis of the two-particle case. Despite its apparent simplicity, the small noise limit of such a system already exhibits various behaviours. In particular, depending on the drift coefficient, the particles can either stick into a cluster, the velocity of which is determined by elementary computations, or drift away from each other at constant velocity, in a random ordering. The persistence of randomness in the zero noise limit is of the very same nature as in the pioneering works by Veretennikov (Mat. Zametki, 1983) and Bafico and Baldi (Stochastics, 1981) concerning the so-called Peano phenomenon. In the case of rank-based processes, we use a simple convexity argument to prove that the small noise limit is described by the sticky particle dynamics introduced by Brenier and Grenier (SIAM J. Numer. Anal., 1998), where particles travel at constant velocity between collisions, at which they stick together. In the general case of order-based processes, we give a sufficient condition on the drift for all the particles to aggregate into a single cluster, and compute the velocity of this cluster. Our argument consists in turning the study of the small noise limit into the study of the long time behaviour of a suitably rescaled process, and then exhibiting a Lyapunov functional for this rescaled process.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Convergence of the Wang-Landau algorithm

We analyze the convergence properties of the Wang-Landau algorithm. This sampling method belongs to the general class of adaptive importance sampling strategies which use the free energy along a chosen reaction coordinate as a bias. Such algorithms are very helpful to enhance the sampling properties of Markov Chain Monte Carlo algorithms, when the dynamics is metastable. We prove the convergence of the Wang-Landau algorithm and an associated central limit theorem.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Propagation of chaos for rank-based interacting diffusions and long time behaviour of a scalar quasilinear parabolic equation

We study a quasilinear parabolic Cauchy problem with a cumulative distribution function on the real line as an initial condition. We call 'probabilistic solution' a weak solution which remains a cumulative distribution function at all times. We prove the uniqueness of such a solution and we deduce the existence from a propagation of chaos result on a system of scalar diffusion processes, the interactions of which only depend on their ranking. We then investigate the long time behaviour of the solution. Using a probabilistic argument and under weak assumptions, we show that the flow of the Wasserstein distance between two solutions is contractive. Under more stringent conditions ensuring the regularity of the probabilistic solutions, we finally derive an explicit formula for the time derivative of the flow and we deduce the convergence of solutions to equilibrium.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

High order discretization schemes for stochastic volatility models

In usual stochastic volatility models, the process driving the volatility of the asset price evolves according to an autonomous one-dimensional stochastic differential equation. We assume that the coefficients of this equation are smooth. Using Itô's formula, we get rid, in the asset price dynamics, of the stochastic integral with respect to the Brownian motion driving this SDE. Taking advantage of this structure, we propose - a scheme, based on the Milstein discretization of this SDE, with order one of weak trajectorial convergence for the asset price, - a scheme, based on the Ninomiya-Victoir discretization of this SDE, with order two of weak convergence for the asset price. We also propose a specific scheme with improved convergence properties when the volatility of the asset price is driven by an Orstein-Uhlenbeck process. We confirm the theoretical rates of convergence by numerical experiments and show that our schemes are well adapted to the multilevel Monte Carlo method introduced by Giles [2008a, 2008b].

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Lévy flights in ecology

We are interested in modeling Darwinian evolution resulting from the interplay of phenotypic variation and natural selection through ecological interactions. The population is modeled as a stochastic point process whose generator captures the probabilistic dynamics over continuous time of birth, mutation, and death, as influenced by each individual's trait values, and interactions between individuals. An offspring usually inherits the trait values of her progenitor, except when a random mutation causes the offspring to take an instantaneous mutation step at birth to new trait values. In the case we are interested in, the probability distribution of mutations has a heavy tail and belongs to the domain of attraction of a stable law. We investigate the large-population limit with allometric demographies: larger populations made up of smaller individuals which reproduce and die faster, as is typical for micro-organisms. We show that depending on the allometry coefficient the limit behavior of the population process can be approximated by nonlinear Lévy flights of different nature: either deterministic, in the form of nonlocal fractional reaction-diffusion equations, or stochastic, as nonlinear super-processes with the underlying reaction and a fractional diffusion operator. These approximation results demonstrate the existence of such nontrivial fractional objects; their uniqueness is also proved.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Convergence of a stochastic particle approximation for fractional scalar conservation laws

We give a probabilistic numerical method for solving a partial differential equation with fractional diffusion and nonlinear drift. The probabilistic interpretation of this equation uses a system of particles driven by Lévy alpha-stable processes and interacting with their drift through their empirical cumulative distribution function. We show convergence to the solution for the associated Euler scheme.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Exact retrospective Monte Carlo computation of arithmetic average Asian options

Taking advantage of the recent litterature on exact simulation algorithms (Beskos, Papaspiliopoulos and Roberts) and unbiased estimation of the expectation of certain fonctional integrals (Wagner, Beskos et al. and Fearnhead et al.), we apply an exact simulation based technique for pricing continuous arithmetic average Asian options in the Black and Scholes framework. Unlike existing Monte Carlo methods, we are no longer prone to the discretization bias resulting from the approximation of continuous time processes through discrete sampling. Numerical results of simulation studies are presented and variance reduction problems are considered.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Regularity of the Exercise Boundary for American Put Options on Assets with Discrete Dividends

We analyze the regularity of the optimal exercise boundary for the American Put option when the underlying asset pays a discrete dividend at a known time $t_d$ during the lifetime of the option. The ex-dividend asset price process is assumed to follow Black-Scholes dynamics and the dividend amount is a deterministic function of the ex-dividend asset price just before the dividend date. The solution to the associated optimal stopping problem can be characterised in terms of an optimal exercise boundary which, in contrast to the case when there are no dividends, may no longer be monotone. In this paper we prove that when the dividend function is positive and concave, then the boundary is non-increasing in a left-hand neighbourhood of $t_d$, and tends to $0$ as time tends to $t_d^-$ with a speed that we can characterize. When the dividend function is linear in a neighbourhood of zero, then we show continuity of the exercise boundary and a high contact principle in the left-hand neighbourhood of $t_d$. When it is globally linear, then right-continuity of the boundary and the high contact principle are proved to hold globally. Finally, we show how all the previous results can be extended to multiple dividend payment dates in that case.

0 citations0 reviewsNo code linkedOpen access
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)$).

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access