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

Mathias Rousset

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

ResearcherUnclaimed source record
math.PRmath.NAcond-mat.stat-mechComputationmath-phmath.APmath.FAmath.MPNumerical Analysis

Catalog footprint

What is connected

11works
9topics
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

11 published item(s)

preprint2020arXiv

Exact targeting of Gibbs distributions using velocity-jump processes

This work introduces and studies a new family of velocity jump Markov processes directly amenable to exact simulation with the following two properties: i) trajectories converge in law when a time-step parameter vanishes towards a given Langevin or Hamil-tonian dynamics; ii) the stationary distribution of the process is always exactly given by the product of a Gaussian (for velocities) by any target log-density whose gradient is pointwise computabe together with some additional explicit appropriate upper bound. The process does not exhibit any velocity reflections (jump sizes can be controlled) and is suitable for the 'factorization method'. We provide a rigorous mathematical proof of: i) the small time-step convergence towards Hamiltonian/Langevin dynamics, as well as ii) the exponentially fast convergence towards the target distribution when suitable noise on velocity is present. Numerical implementation is detailed and illustrated.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

More on the long time stability of Feynman-Kac semigroups

Feynman-Kac semigroups appear in various areas of mathematics: non-linear filtering, large deviations theory, spectral analysis of Schrodinger operators among others. Their long time behavior provides important information, for example in terms of ground state energy of Schrodinger operators, or scaled cumulant generating function in large deviations theory. In this paper, we propose a simple and natural extension of the stability of Markov chains for these non-linear evolutions. As other classical ergodicity results, it relies on two assumptions: a Lyapunov condition that induces some compactness, and a minorization condition ensuring some mixing. Illustrative examples are provided, where the stability of the non-linear semigroup arises either from the underlying dynamics or from the Feynman-Kac weight function. We also use our technique to provide uniform in the time step convergence estimates for discretizations of stochastic differential equations

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A $N$-uniform quantitative Tanaka's theorem for the conservative Kac's $N$-particle system with Maxwell molecules

This paper considers the space homogenous Boltzmann equation with Maxwell molecules and arbitrary angular distribution. Following Kac's program, emphasis is laid on the the associated conservative Kac's stochastic $N$-particle system, a Markov process with binary collisions conserving energy and total momentum. An explicit Markov coupling (a probabilistic, Markovian coupling of two copies of the process) is constructed, using simultaneous collisions, and parallel coupling of each binary random collision on the sphere of collisional directions. The euclidean distance between the two coupled systems is almost surely decreasing with respect to time, and the associated quadratic coupling creation (the time variation of the averaged squared coupling distance) is computed explicitly. Then, a family (indexed by $δ> 0$) of $N$-uniform ''weak'' coupling / coupling creation inequalities are proven, that leads to a $N$-uniform power law trend to equilibrium of order ${\sim}_{ t \to + \infty} t^{-δ} $, with constants depending on moments of the velocity distributions strictly greater than $2(1 + δ)$. The case of order $4$ moment is treated explicitly, achieving Kac's program without any chaos propagation analysis. Finally, two counter-examples are suggested indicating that the method: (i) requires the dependance on $>2$-moments, and (ii) cannot provide contractivity in quadratic Wasserstein distance in any case.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Analysis of Adaptive Multilevel Splitting algorithms in an idealized case

The Adaptive Multilevel Splitting algorithm is a very powerful and versatile method to estimate rare events probabilities. It is an iterative procedure on an interacting particle system, where at each step, the $k$ less well-adapted particles among $n$ are killed while $k$ new better adapted particles are resampled according to a conditional law. We analyze the algorithm in the idealized setting of an exact resampling and prove that the estimator of the rare event probability is unbiased whatever $k$. We also obtain a precise asymptotic expansion for the variance of the estimator and the cost of the algorithm in the large $n$ limit, for a fixed $k$.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Scalable and Quasi-Contractive Markov Coupling of Maxwell Collision

This paper considers space homogenous Boltzmann kinetic equations in dimension $d$ with Maxwell collisions (and without Grad's cut-off). An explicit Markov coupling of the associated conservative (Nanbu) stochastic $N$-particle system is constructed, using plain parallel coupling of isotropic random walks on the sphere of two-body collisional directions. The resulting coupling is almost surely decreasing, and the $L_2$-coupling creation is computed explicitly. Some quasi-contractive and uniform in $N$ coupling / coupling creation inequalities are then proved, relying on $2+α$-moments ($α>0$) of velocity distributions; upon $N$-uniform propagation of moments of the particle system, it yields a $N$-scalable $α$-power law trend to equilibrium. The latter are based on an original sharp inequality, which bounds from above the coupling distance of two centered and normalized random variables $(U,V)$ in $\R^d$, with the average square parallelogram area spanned by $(U-U_\ast,V-V_\ast)$, $(U_\ast,V_\ast)$ denoting an independent copy. Two counter-examples proving the necessity of the dependance on $>2$-moments and the impossibility of strict contractivity are provided. The paper, (mostly) self-contained, does not require any propagation of chaos property and uses only elementary tools.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

On a probabilistic interpretation of shape derivatives of Dirichlet groundstates with application to Fermion nodes

This paper considers Schrödinger operators, and presents a probabilistic interpretation of the variation (or shape derivative) of the Dirichlet groundstate energy when the associated domain is perturbed. This interpretation relies on the distribution on the boundary of a stopped random process with Feynman-Kac weights. Practical computations require in addition the explicit approximation of the normal derivative of the groundstate on the boundary. We then propose to use this formulation in the case of the so-called fixed node approximation of fermion groundstates, defined by the bottom eigenelements of the Schrödinger operator of a fermionic system with Dirichlet conditions on the nodes (the set of zeros) of an initially guessed skew-symmetric function. We show that the shape derivative of the fixed node energy vanishes if and only if either (i) the distribution on the nodes of the stopped random process is symmetric; or (ii) the nodes are exactly the zeros of a skew-symmetric eigenfunction of the operator. We propose an approximation of the shape derivative of the fixed node energy that can be computed with a Monte-Carlo algorithm, which can be referred to as Nodal Monte-Carlo (NMC). The latter approximation of the shape derivative also vanishes if and only if either (i) or (ii) holds.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Individual-based models for bacterial chemotaxis in the diffusion asymptotics

We discuss velocity-jump models for chemotaxis of bacteria with an internal state that allows the velocity jump rate to depend on the memory of the chemoattractant concentration along their path of motion. Using probabilistic techniques, we provide a pathwise result that shows that the considered process converges to an advection-diffusion process in the (long-time) diffusion limit. We also (re-)prove using the same approach that the same limiting equation arises for a related, simpler process with direct sensing of the chemoattractant gradient. Additionally, we propose a time discretization technique that retains these diffusion limits exactly, i.e., without error that depends on the time discretization. In the companion paper \cite{variance}, these results are used to construct a coupling technique that allows numerical simulation of the process with internal state with asymptotic variance reduction, in the sense that the variance vanishes in the diffusion limit.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Langevin dynamics with constraints and computation of free energy differences

In this paper, we consider Langevin processes with mechanical constraints. The latter are a fundamental tool in molecular dynamics simulation for sampling purposes and for the computation of free energy differences. The results of this paper can be divided into three parts. (i) We propose a simple discretization of the constrained Langevin process based on a standard splitting strategy. We show how to correct the scheme so that it samples {\em exactly} the canonical measure restricted on a submanifold, using a Metropolis rule in the spirit of the Generalized Hybrid Monte Carlo (GHMC) algorithm. Moreover, we obtain, in some limiting regime, a consistent discretization of the overdamped Langevin (Brownian) dynamics on a submanifold, also sampling exactly the correct canonical measure with constraints. The corresponding numerical methods can be used to sample (without any bias) a probability measure supported by a submanifold. (ii) For free energy computation using thermodynamic integration, we rigorously prove that the longtime average of the Lagrange multipliers of the constrained Langevin dynamics yields the gradient of a rigid version of the free energy associated with the constraints. A second order time discretization using the Lagrange multipliers is proposed. (iii) The Jarzynski-Crooks fluctuation relation is proved for Langevin processes with mechanical constraints evolving in time. An original numerical discretization without time-step error is proposed. Numerical illustrations are provided for (ii) and (iii).

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Simulating individual-based models of bacterial chemotaxis with asymptotic variance reduction

We discuss variance reduced simulations for an individual-based model of chemotaxis of bacteria with internal dynamics. The variance reduction is achieved via a coupling of this model with a simpler process in which the internal dynamics has been replaced by a direct gradient sensing of the chemoattractants concentrations. In the companion paper \cite{limits}, we have rigorously shown, using a pathwise probabilistic technique, that both processes converge towards the same advection-diffusion process in the diffusive asymptotics. In this work, a direct coupling is achieved between paths of individual bacteria simulated by both models, by using the same sets of random numbers in both simulations. This coupling is used to construct a hybrid scheme with reduced variance. We first compute a deterministic solution of the kinetic density description of the direct gradient sensing model; the deviations due to the presence of internal dynamics are then evaluated via the coupled individual-based simulations. We show that the resulting variance reduction is \emph{asymptotic}, in the sense that, in the diffusive asymptotics, the difference between the two processes has a variance which vanishes according to the small parameter.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Computation of free energy differences through nonequilibrium stochastic dynamics: the reaction coordinate case

The computation of free energy differences through an exponential weighting of out of equilibrium paths (known as the Jarzynski equality) is often used for transitions between states described by an external parameter $λ$ in the Hamiltonian. We present here an extension to transitions between states defined by different values of some reaction coordinate, using a projected Brownian dynamics. In contrast with other approaches, we use a projection rather than a constraining potential to let the constraints associated with the reaction coordinate evolve. We show how to use the Lagrange multipliers associated with these constraints to compute the work associated with a given trajectory. Appropriate discretizations are proposed. Some numerical results demonstrate the applicability of the method for the computation of free energy difference profiles.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Computation of free energy profiles with parallel adaptive dynamics

We propose a formulation of adaptive computation of free energy differences, in the ABF or nonequilibrium metadynamics spirit, using conditional distributions of samples of configurations which evolve in time. This allows to present a truly unifying framework for these methods, and to prove convergence results for certain classes of algorithms. From a numerical viewpoint, a parallel implementation of these methods is very natural, the replicas interacting through the reconstructed free energy. We show how to improve this parallel implementation by resorting to some selection mechanism on the replicas. This is illustrated by computations on a model system of conformational changes.

0 citations0 reviewsNo code linkedOpen access