BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Gabriel Stoltz

Gabriel Stoltz contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math-phmath.MPmath.PRcond-mat.stat-mechmath.FABiomoleculescond-mat.othercond-mat.str-elmath.APmath.NAmath.SPNumerical Analysisphysics.comp-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
9works
0followers
13topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

9 published item(s)

preprint2022arXiv

Error estimates and variance reduction for nonequilibrium stochastic dynamics

Equilibrium properties in statistical physics are obtained by computing averages with respect to Boltzmann-Gibbs measures, sampled in practice using ergodic dynamics such as the Langevin dynamics. Some quantities however cannot be computed by simply sampling the Boltzmann-Gibbs measure, in particular transport coefficients, which relate the current of some physical quantity of interest to the forcing needed to induce it. For instance, a temperature difference induces an energy current, the proportionality factor between these two quantities being the thermal conductivity. From an abstract point of view, transport coefficients can also be considered as some form of sensitivity analysis with respect to an added forcing to the baseline dynamics. There are various numerical techniques to estimate transport coefficients, which all suffer from large errors, in particular large statistical errors. This contribution reviews the most popular methods, namely the Green-Kubo approach where the transport coefficient is expressed as some time-integrated correlation function, and the approach based on longtime averages of the stochastic dynamics perturbed by an external driving (so-called nonequilibrium molecular dynamics). In each case, the various sources of errors are made precise, in particular the bias related to the time discretization of the underlying continuous dynamics, and the variance of the associated Monte Carlo estimators. Some recent alternative techniques to estimate transport coefficients are also discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Weighted $L^2$-contractivity of Langevin dynamics with singular potentials

Convergence to equilibrium of underdamped Langevin dynamics is studied under general assumptions on the potential $U$ allowing for singularities. By modifying the direct approach to convergence in $L^2$ pioneered by F. Hérau and developped by Dolbeault, Mouhot and Schmeiser, we show that the dynamics converges exponentially fast to equilibrium in the topologies $L^2(dμ)$ and $L^2(W^* dμ)$, where $μ$ denotes the invariant probability measure and $W^*$ is a suitable Lyapunov weight. In both norms, we make precise how the exponential convergence rate depends on the friction parameter $γ$ in Langevin dynamics, by providing a lower bound scaling as $\min(γ, γ^{-1})$. The results hold for usual polynomial-type potentials as well as potentials with singularities such as those arising from pairwise Lennard-Jones interactions between particles.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Coulomb gases under constraint: some theoretical and numerical results

We consider Coulomb gas models for which the empirical measure typically concentrates, when the number of particles becomes large, on an equilibrium measure minimizing an electrostatic energy. We study the behavior when the gas is conditioned on a rare event. We first show that the special case of quadratic confinement and linear constraint is exactly solvable due to a remarkable factorization, and that the conditioning has then the simple effect of shifting the cloud of particles without deformation. To address more general cases, we perform a theoretical asymptotic analysis relying on a large deviations technique known as the Gibbs conditioning principle. The technical part amounts to establishing that the conditioning ensemble is an I-continuity set of the energy. This leads to characterizing the conditioned equilibrium measure as the solution of a modified variational problem. For simplicity, we focus on linear statistics and on quadratic statistics constraints. Finally, we numerically illustrate our predictions and explore cases in which no explicit solution is known. For this, we use a Generalized Hybrid Monte Carlo algorithm for sampling from the conditioned distribution for a finite but large system.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Large deviations of empirical measures of diffusions in weighted topologies

We consider large deviations of empirical measures of diffusion processes. In a first part, we present conditions to obtain a large deviations principle (LDP) for a precise class of unbounded functions. This provides an analogue to the standard Cramér condition in the context of diffusion processes, which turns out to be related to a spectral gap condition for a Witten-Schrödinger operator. Secondly, we study more precisely the properties of the Donsker-Varadhan rate functional associated with the LDP. We revisit and generalize some standard duality results as well as a more original decomposition of the rate functional with respect to the symmetric and antisymmetric parts of the dynamics. Finally, we apply our results to overdamped and underdamped Langevin dynamics, showing the applicability of our framework for degenerate diffusions in unbounded configuration spaces.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Machine learning force fields and coarse-grained variables in molecular dynamics: application to materials and biological systems

Machine learning encompasses a set of tools and algorithms which are now becoming popular in almost all scientific and technological fields. This is true for molecular dynamics as well, where machine learning offers promises of extracting valuable information from the enormous amounts of data generated by simulation of complex systems. We provide here a review of our current understanding of goals, benefits, and limitations of machine learning techniques for computational studies on atomistic systems, focusing on the construction of empirical force fields from ab-initio databases and the determination of reaction coordinates for free energy computation and enhanced sampling.

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
preprint2010arXiv

Free energy calculations: An efficient adaptive biasing potential method

We develop an efficient sampling and free energy calculation technique within the adaptive biasing potential (ABP) framework. By mollifying the density of states we obtain an approximate free energy and an adaptive bias potential that is computed directly from the population along the coordinates of the free energy. Because of the mollifier, the bias potential is "nonlocal" and its gradient admits a simple analytic expression. A single observation of the reaction coordinate can thus be used to update the approximate free energy at every point within a neighborhood of the observation. This greatly reduces the equilibration time of the adaptive bias potential. This approximation introduces two parameters: strength of mollification and the zero of energy of the bias potential. While we observe that the approximate free energy is a very good estimate of the actual free energy for a large range of mollification strength, we demonstrate that the errors associated with the mollification may be removed via deconvolution. The zero of energy of the bias potential, which is easy to choose, influences the speed of convergence but not the limiting accuracy. This method is simple to apply to free energy or mean force computation in multiple dimensions and does not involve second derivatives of the reaction coordinates, matrix manipulations nor on-the-fly adaptation of parameters. For the alanine dipeptide test case, the new method is found to gain as much as a factor of ten in efficiency as compared to two common adaptive biasing force formulations and it is shown to be as efficient as well-tempered metadynamics with the post-process deconvolution giving a clear advantage to the mollified density of states method.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Gell-Mann and Low formula for degenerate unperturbed states

The Gell-Mann and Low switching allows to transform eigenstates of an unperturbed Hamiltonian $H_0$ into eigenstates of the modified Hamiltonian $H_0 + V$. This switching can be performed when the initial eigenstate is not degenerate, under some gap conditions with the remainder of the spectrum. We show here how to extend this approach to the case when the ground state of the unperturbed Hamiltonian is degenerate. More precisely, we prove that the switching procedure can still be performed when the initial states are eigenstates of the finite rank self-adjoint operator $\cP_0 V \cP_0$, where $\cP_0$ is the projection onto a degenerate eigenspace of $H_0$.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The many-body Green function of degenerate systems

A rigorous non perturbative adiabatic approximation of the evolution operator in the many-body physics of degenerate systems is derived. This approximation is used to solve the long-standing problem of the choice of the initial states of H0 leading to eigenstates of H0+V for degenerate systems. These initial states are eigenstates of P0 V P0, where P0 is the projection onto a degenerate eigenspace of H0. This result is used to give the proper definition of the Green function, the statistical Green function and the non-equilibrium Green function of degenerate systems. The convergence of these Green functions is established.

0 citations0 reviewsNo code linkedOpen access