Source author record

Tony Lelièvre

Tony Lelièvre 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

31works
15topics
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

31 published item(s)

preprint2022arXiv

A probabilistic study of the kinetic Fokker-Planck equation in cylindrical domains

We consider classical solutions to the kinetic Fokker-Planck equation on a bounded domain $\mathcal O \subset~\mathbb{R}^d$ in position, and we obtain a probabilistic representation of the solutions using the Langevin diffusion process with absorbing boundary conditions on the boundary of the phase-space cylindrical domain $D = \mathcal O \times \mathbb{R}^d$. Furthermore, a Harnack inequality, as well as a maximum principle, is provided on $D$ for solutions to this kinetic Fokker-Planck equation, together with the existence of a smooth transition density for the associated absorbed Langevin process. This transition density is shown to satisfy an explicit Gaussian upper-bound. Finally, the continuity and positivity of this transition density at the boundary of $D$ is also studied. All these results are in particular crucial to study the behavior of the Langevin diffusion process when it is trapped in a metastable state defined in terms of positions.

preprint2022arXiv

Eyring-Kramers exit rates for the overdamped Langevin dynamics: the case with saddle points on the boundary

Let $(X_t)_{t\ge 0}$ be the stochastic process solution to the overdamped Langevin dynamics $$dX_t=-\nabla f(X_t) \, dt +\sqrt h \, dB_t$$ and let $Ω\subset \mathbb R^d $ be the basin of attraction of a local minimum of $f: \mathbb R^d \to \mathbb R$. Up to a small perturbation of $Ω$ to make it smooth, we prove that the exit rates of $(X_t)_{t\ge 0}$ from $Ω$ through each of the saddle points of $f$ on $\partial Ω$ can be parametrized by the celebrated Eyring-Kramers laws, in the limit $h \to 0$. This result provides firm mathematical grounds to jump Markov models which are used to model the evolution of molecular systems, as well as to some numerical methods which use these underlying jump Markov models to efficiently sample metastable trajectories of the overdamped Langevin dynamics.

preprint2021arXiv

The Adaptive Biasing Force algorithm with non-conservative forces and related topics

We propose a study of the Adaptive Biasing Force method's robustness under generic (possibly non-conservative) forces. We first ensure the flat histogram property is satisfied in all cases. We then introduce a fixed point problem yielding the existence of a stationary state for both the Adaptive Biasing Force and Projected Adapted Biasing Force algorithms, relying on generic bounds on the invariant probability measures of homogeneous diffusions. Using classical entropy techniques, we prove the exponential convergence of both biasing force and law as time goes to infinity, for both the Adaptive Biasing Force and the Projected Adaptive Biasing Force methods.

preprint2020arXiv

Adaptive force biasing algorithms: new convergence results and tensor approximations of the bias

A modification of the Adaptive Biasing Force method is introduced, in which the free energy is approximated by a sum of tensor products of one-dimensional functions. This enables to handle a larger number of reaction coordinates than the classical algorithm. We prove the algorithm is well-defined and prove the long-time convergence toward a regularized version of the free energy for an idealized version of the algorithm. Numerical experiments demonstrate that the method is able to capture correlations between reaction coordinates.

preprint2019arXiv

Local and Global Perspectives on Diffusion Maps in the Analysis of Molecular Systems

Diffusion maps approximate the generator of Langevin dynamics from simulation data. They afford a means of identifying the slowly-evolving principal modes of high-dimensional molecular systems. When combined with a biasing mechanism, diffusion maps can accelerate the sampling of the stationary Boltzmann-Gibbs distribution. In this work, we contrast the local and global perspectives on diffusion maps, based on whether or not the data distribution has been fully explored. In the global setting, we use diffusion maps to identify metastable sets and to approximate the corresponding committor functions of transitions between them. We also discuss the use of diffusion maps within the metastable sets, formalising the locality via the concept of the quasi-stationary distribution and justifying the convergence of diffusion maps within a local equilibrium. This perspective allows us to propose an enhanced sampling algorithm. We demonstrate the practical relevance of these approaches both for simple models and for molecular dynamics problems (alanine dipeptide and deca-alanine).

preprint2018arXiv

Effective dynamics for non-reversible stochastic differential equations: a quantitative study

Coarse-graining is central to reducing dimensionality in molecular dynamics, and is typically characterized by a mapping which projects the full state of the system to a smaller class of variables. While extensive literature has been devoted to coarse-graining starting from reversible systems, not much is known in the non-reversible setting. In this article, starting with a non-reversible dynamics, we introduce and study an effective dynamics which approximates the (non-closed) projected dynamics. Under fairly weak conditions on the system, we prove error bounds on the trajectorial error between the projected and the effective dynamics. In addition to extending existing results to the non-reversible setting, our error estimates also indicate that the notion of mean force motivated by this effective dynamics is a good one.

preprint2015arXiv

Accelerated dynamics: Mathematical foundations and algorithmic improvements

We present a review of recent works on the mathematical analysis of algorithms which have been proposed by A.F. Voter and co-workers in the late nineties in order to efficiently generate long trajectories of metastable processes. These techniques have been successfully applied in many contexts, in particular in the field of materials science. The mathematical analysis we propose relies on the notion of quasi stationary distribution.

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.

preprint2015arXiv

Local and global solution for a nonlocal Fokker-Planck equation related to the adaptive biaising force processes

We prove global existence, uniqueness and regularity of the mild, Lp and classical solution of a non-linear Fokker-Planck equation arising in an adaptive importance sampling method for molecular dynamics calculations. The non- linear term is related to a conditional expectation, and is thus non-local. The proof uses tools from the theory of semigroups of linear operators for the local existence result, and an a priori estimate based on a supersolution for the global existence result.

preprint2015arXiv

Long-time convergence of an adaptive biasing force method: Variance reduction by Helmholtz projection

In this paper, we propose an improvement of the adaptive biasing force (ABF) method, by projecting the estimated mean force onto a gradient. The associated stochastic process satisfies a non linear stochastic differential equation. Using entropy techniques, we prove exponential convergence to the stationary state of this stochastic process. We finally show on some numerical examples that the variance of the approximated mean force is reduced using this technique, which makes the algorithm more efficient than the standard ABF method.

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]).

preprint2014arXiv

A Generalized Parallel Replica Dynamics

Metastability is a common obstacle to performing long molecular dynamics simulations. Many numerical methods have been proposed to overcome it. One method is parallel replica dynamics, which relies on the rapid convergence of the underlying stochastic process to a quasi-stationary distribution. Two requirements for applying parallel replica dynamics are knowledge of the time scale on which the process converges to the quasi-stationary distribution and a mechanism for generating samples from this distribution. By combining a Fleming-Viot particle system with convergence diagnostics to simultaneously identify when the process converges while also generating samples, we can address both points. This variation on the algorithm is illustrated with various numerical examples, including those with entropic barriers and the 2D Lennard-Jones cluster of seven atoms.

preprint2014arXiv

Free energy calculations from adaptive molecular dynamics simulations with adiabatic reweighting

We propose an adiabatic reweighting algorithm for computing the free energy along an external parameter from adaptive molecular dynamics simulations. The adaptive bias is estimated using Bayes identity and information from all the sampled configurations. We apply the algorithm to a structural transition in a cluster and to the migration of a crystalline defect along a reaction coordinate. Compared to standard adaptive molecular dynamics, we observe an acceleration of convergence. With the aid of the algorithm, it is also possible to iteratively construct the free energy along the reaction coordinate without having to differentiate the gradient of the reaction coordinate or any biasing potential.

preprint2014arXiv

Macroscopic limit of a one-dimensional model for aging fluids

We study a one-dimensional equation arising in the multiscale modeling of some non-Newtonian fluids. At a given shear rate, the equation provides the instantaneous mesoscopic response of the fluid, allowing to compute the corresponding stress. In a simple setting, we study the well-posedness of the equation and next the long-time behavior of its solution. In the limit of a response of the fluid much faster than the time variations of the ambient shear rate, we derive some equivalent macroscopic differential equations that relate the shear rate and the stress. Our analytical conclusions are confronted to some numerical experiments. The latter quantitatively confirm our derivations.

preprint2014arXiv

Mathematical Analysis of Temperature Accelerated Dynamics

We give a mathematical framework for temperature accelerated dynamics (TAD), an algorithm proposed by M.R. Sørensen and A.F. Voter to efficiently generate metastable stochastic dynamics. Using the notion of quasistationary distributions, we propose some modifications to TAD. Then considering the modified algorithm in an idealized setting, we show how TAD can be made mathematically rigorous.

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.

preprint2014arXiv

The parallel replica method for simulating long trajectories of Markov chains

The parallel replica dynamics, originally developed by A.F. Voter, efficiently simulates very long trajectories of metastable Langevin dynamics. We present an analogous algorithm for discrete time Markov processes. Such Markov processes naturally arise, for example, from the time discretization of a continuous time stochastic dynamics. Appealing to properties of quasistationary distributions, we show that our algorithm reproduces exactly (in some limiting regime) the law of the original trajectory, coarsened over the metastable states.

preprint2013arXiv

A nonintrusive method to approximate linear systems with nonlinear parameter dependence

We consider a family of linear systems $A_μα=C$ with system matrix $A_μ$ depending on a parameter $μ$ and for simplicity parameter-independent right-hand side $C$. These linear systems typically result from the finite-dimensional approximation of a parameter-dependent boundary-value problem. We derive a procedure based on the Empirical Interpolation Method to obtain a separated representation of the system matrix in the form $A_μ\approx\sum_{m}β_m(μ)A_{μ_m}$ for some selected values of the parameter. Such a separated representation is in particular useful in the Reduced Basis Method. The procedure is called nonintrusive since it only requires to access the matrices $A_{μ_m}$. As such, it offers a crucial advantage over existing approaches that instead derive separated representations requiring to enter the code at the level of assembly. Numerical examples illustrate the performance of our new procedure on a simple one-dimensional boundary-value problem and on three-dimensional acoustic scattering problems solved by a boundary element method.

preprint2013arXiv

Accurate and efficient evaluation of the a posteriori error estimator in the reduced basis method

The reduced basis method is a model reduction technique yielding substantial savings of computational time when a solution to a parametrized equation has to be computed for many values of the parameter. Certification of the approximation is possible by means of an a posteriori error bound. Under appropriate assumptions, this error bound is computed with an algorithm of complexity independent of the size of the full problem. In practice, the evaluation of the error bound can become very sensitive to round-off errors. We propose herein an explanation of this fact. A first remedy has been proposed in [F. Casenave, Accurate \textit{a posteriori} error evaluation in the reduced basis method. \textit{C. R. Math. Acad. Sci. Paris} \textbf{350} (2012) 539--542.]. Herein, we improve this remedy by proposing a new approximation of the error bound using the Empirical Interpolation Method (EIM). This method achieves higher levels of accuracy and requires potentially less precomputations than the usual formula. A version of the EIM stabilized with respect to round-off errors is also derived. The method is illustrated on a simple one-dimensional diffusion problem and a three-dimensional acoustic scattering problem solved by a boundary element method.

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.

preprint2013arXiv

Greedy algorithms for high-dimensional eigenvalue problems

In this article, we present two new greedy algorithms for the computation of the lowest eigenvalue (and an associated eigenvector) of a high-dimensional eigenvalue problem, and prove some convergence results for these algorithms and their orthogonalized versions. The performance of our algorithms is illustrated on numerical test cases (including the computation of the buckling modes of a microstructured plate), and compared with that of another greedy algorithm for eigenvalue problems introduced by Ammar and Chinesta.

preprint2013arXiv

Low temperature asymptotics for Quasi-Stationary Distributions in a bounded domain

We analyze the low temperature asymptotics of the quasi-stationary distribution associated with the overdamped Langevin dynamics (a.k.a. the Einstein-Smoluchowski diffusion equation) in a bounded domain. This analysis is useful to rigorously prove the consistency of an algorithm used in molecular dynamics (the hyperdynamics), in the small temperature regime. More precisely, we show that the algorithm is exact in terms of state-to-state dynamics up to exponentially small factor in the limit of small temperature. The proof is based on the asymptotic spectral analysis of associated Dirichlet and Neumann realizations of Witten Laplacians. In order to cover a reasonably large range of applications, the usual assumptions that the energy landscape is a Morse function has been relaxed as much as possible.

preprint2012arXiv

Derivation of Langevin Dynamics in a Nonzero Background Flow Field

We propose a derivation of a nonequilibrium Langevin dynamics for a large particle immersed in a background flow field. A single large particle is placed in an ideal gas heat bath composed of point particles that are distributed consistently with the background flow field and that interact with the large particle through elastic collisions. In the limit of small bath atom mass, the large particle dynamics converges in law to a stochastic dynamics. This derivation follows the ideas of [D. Dürr, S. Goldstein, and J. L. Lebowitz, 1981 and 1983; P. Calderoni, D. Dürr, and S. Kusuoka, 1989] and provides extensions to handle the nonzero background flow. The derived nonequilibrium Langevin dynamics is similar to the dynamics in [M. McPhie, et al., 2001]. Some numerical experiments illustrate the use of the obtained dynamic to simulate homogeneous liquid materials under flow.

preprint2012arXiv

On the length of one-dimensional reactive paths

Motivated by some numerical observations on molecular dynamics simulations, we analyze metastable trajectories in a very simplecsetting, namely paths generated by a one-dimensional overdamped Langevin equation for a double well potential. More precisely, we are interested in so-called reactive paths, namely trajectories which leave definitely one well and reach the other one. The aim of this paper is to precisely analyze the distribution of the lengths of reactive paths in the limit of small temperature, and to compare the theoretical results to numerical results obtained by a Monte Carlo method, namely the multi-level splitting approach.

preprint2012arXiv

Two mathematical tools to analyze metastable stochastic processes

We present how entropy estimates and logarithmic Sobolev inequalities on the one hand, and the notion of quasi-stationary distribution on the other hand, are useful tools to analyze metastable overdamped Langevin dynamics, in particular to quantify the degree of metastability. We discuss the interest of these approaches to estimate the efficiency of some classical algorithms used to speed up the sampling, and to evaluate the error introduced by some coarse-graining procedures. This paper is a summary of a plenary talk given by the author at the ENUMATH 2011 conference.

preprint2011arXiv

Periodic long-time behaviour for an approximate model of nematic polymers

We study the long-time behaviour of a nonlinear Fokker-Planck equation, which models the evolution of rigid polymers in a given flow, after a closure approximation. The aim of this work is twofold: first, we propose a microscopic derivation of the classical Doi closure, at the level of the kinetic equation ; second, we prove the convergence of the solution to the Fokker-Planck equation to periodic solutions in the long-time limit.

preprint2010arXiv

A multiple replica approach to simulate reactive trajectories

A method to generate reactive trajectories, namely equilibrium trajectories leaving a metastable state and ending in another one is proposed. The algorithm is based on simulating in parallel many copies of the system, and selecting the replicas which have reached the highest values along a chosen one-dimensional reaction coordinate. This reaction coordinate does not need to precisely describe all the metastabilities of the system for the method to give reliable results. An extension of the algorithm to compute transition times from one metastable state to another one is also presented. We demonstrate the interest of the method on two simple cases: a one-dimensional two-well potential and a two-dimensional potential exhibiting two channels to pass from one metastable state to another one.

preprint2010arXiv

Enhanced sampling of multidimensional free-energy landscapes using adaptive biasing forces

We propose an adaptive biasing algorithm aimed at enhancing the sampling of multimodal measures by Langevin dynamics. The underlying idea consists in generalizing the standard adaptive biasing force method commonly used in conjunction with molecular dynamics to handle in a more effective fashion multidimensional reaction coordinates. The proposed approach is anticipated to be particularly useful for reaction coordinates, the components of which are weakly coupled, as illuminated in a mathematical analysis of the long-time convergence of the algorithm. The strength as well as the intrinsic limitation of the method are discussed and illustrated in two realistic test cases.

preprint2010arXiv

Reduced basis techniques for stochastic problems

We report here on the recent application of a now classical general reduction technique, the Reduced-Basis approach initiated in [C. Prud'homme, D. Rovas, K. Veroy, Y. Maday, A. T. Patera, and G. Turinici. Reliable real-time solution of parametrized partial differential equations: Reduced-basis output bounds methods. Journal of Fluids Engineering, 124(1):7080, 2002.], to the specific context of differential equations with random coefficients. After an elementary presentation of the approach, we review two contributions of the authors: [S. Boyaval, C. Le Bris, Y. Maday, N.C. Nguyen, and A.T. Patera. A reduced basis approach for variational problems with stochastic parameters: Application to heat conduction with variable Robin co-efficient. Computer Methods in Applied Mechanics and Engineering, 198(4144):3187-3206, 2009.], which presents the application of the RB approach for the discretization of a simple second order elliptic equation supplied with a random boundary condition, and [S. Boyaval and T. Lelièvre, A variance reduction method for parametrized stochastic differential equations using the reduced basis paradigm with T. Lelièvre, Commun. Math. Sci. 8, special Issue "Mathematical Issue on Complex Fluids" P. Zhang ed., to appear, 2010, ARXIV preprint arXiv:0906.3600], which uses a RB type approach to reduce the variance in the Monte-Carlo simulation of a stochastic differential equation. We conclude the review with some general comments and also discuss possible tracks for further research in the direction.

preprint2006arXiv

An efficient sampling algorithm for Variational Monte Carlo

We propose a new algorithm for sampling the $N$-body density $|Ψ({\bf R})|^2/\int_{\mathbb{R}^{3N}} |Ψ|^2$ in the Variational Monte Carlo (VMC) framework. This algorithm is based upon a modified Ricci-Ciccotti discretization of the Langevin dynamics in the phase space $({\bf R},{\bf P})$ improved by a Metropolis acceptation/rejection step. We show through some representative numerical examples (Lithium, Fluorine and Copper atoms, and phenol molecule), that this algorithm is superior to the standard sampling algorithm based on the biased random walk (importance sampling).