Pierre Henry-Labordere

Pierre Henry-Labordere, Nadia Oudjane, Xiaolu Tan et al.

We provide a representation result of parabolic semi-linear PD-Es, with polynomial nonlinearity, by branching diffusion processes. We extend the classical representation for KPP equations, introduced by Skorokhod (1964), Watanabe (1965) and McKean (1975), by allowing for polynomial nonlinearity in the pair $(u, Du)$, where $u$ is the solution of the PDE with space gradient $Du$. Similar to the previous literature, our result requires a non-explosion condition which restrict to "small maturity" or "small nonlinearity" of the PDE. Our main ingredient is the automatic differentiation technique as in Henry Labordere, Tan and Touzi (2015), based on the Malliavin integration by parts, which allows to account for the nonlinearities in the gradient. As a consequence, the particles of our branching diffusion are marked by the nature of the nonlinearity. This new representation has very important numerical implications as it is suitable for Monte Carlo simulation. Indeed, this provides the first numerical method for high dimensional nonlinear PDEs with error estimate induced by the dimension-free Central limit theorem. The complexity is also easily seen to be of the order of the squared dimension. The final section of this paper illustrates the efficiency of the algorithm by some high dimensional numerical experiments.

Pierre Henry-Labordere, Xiaolu Tan, Nizar Touzi

We propose an unbiased Monte-Carlo estimator for $\mathbb{E}[g(X_{t_1}, \cdots, X_{t_n})]$, where $X$ is a diffusion process defined by a multi-dimensional stochastic differential equation (SDE). The main idea is to start instead from a well-chosen simulatable SDE whose coefficients are updated at independent exponential times. Such a simulatable process can be viewed as a regime-switching SDE, or as a branching diffusion process with one single living particle at all times. In order to compensate for the change of the coefficients of the SDE, our main representation result relies on the automatic differentiation technique induced by Bismu-Elworthy-Li formula from Malliavin calculus, as exploited by Fournié et al.(1999) for the simulation of the Greeks in financial applications. In particular, this algorithm can be considered as a variation of the (infinite variance) estimator obtained in Bally and Kohatsu-Higa [Section 6.1](2014) as an application of the parametrix method.

Pierre Henry-Labordere, Nizar Touzi

We provide a probabilistic representations of the solution of some semilinear hyperbolicand high-order PDEs based on branching diffusions. These representations pave theway for a Monte-Carlo approximation of the solution, thus bypassing the curse ofdimensionality. We illustrate the numerical implications in the context of some popularPDEs in physics such as nonlinear Klein-Gordon equation, a simplied scalar versionof the Yang-Mills equation, a fourth-order nonlinear beam equation and the Gross-Pitaevskii PDEas an example of nonlinear Schrodinger equations.

Pierre Henry-Labordere, Xiaolu Tan, Nizar Touzi

We generalize the algorithm for semi-linear parabolic PDEs in Henry-Labordère (2012) to the non-Markovian case for a class of Backward SDEs (BSDEs). By simulating the branching process, the algorithm does not need any backward regression. To prove that the numerical algorithm converges to the solution of BSDEs, we use the notion of viscosity solution of path dependent PDEs introduced by Ekren, Keller, Touzi and Zhang (2012) and extended in Ekren, Touzi and Zhang (2013).