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

Gabriel Stoltz

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

ResearcherUnclaimed source record
cond-mat.stat-mechmath-phmath.MPcond-mat.mtrl-scimath.PRmath.FAmath.NAcond-mat.othercond-mat.str-elmath.STphysics.comp-phStatistics TheoryBiomoleculesComputationcond-mat.softhep-thmath.APmath.SPNumerical Analysis

Catalog footprint

What is connected

35works
19topics
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)

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
preprint2016arXiv

An accurate scheme to solve cluster dynamics equations using a Fokker-Planck approach

We present a numerical method to accurately simulate particle size distributions within the formalism of rate equation cluster dynamics. This method is based on a discretization of the associated Fokker-Planck equation. We show that particular care has to be taken to discretize the advection part of the Fokker-Planck equation, in order to avoid distortions of the distribution due to numerical diffusion. For this purpose we use the Kurganov-Noelle-Petrova scheme coupled with the monotonicity-preserving reconstruction MP5, which leads to very accurate results. The interest of the method is highlighted on the case of loop coarsening in aluminum. We show that the choice of the models to describe the energetics of loops does not significantly change the normalized loop distribution, while the choice of the models for the absorption coefficients seems to have a significant impact on it.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Error Analysis of Modified Langevin Dynamics

We consider Langevin dynamics associated with a modified kinetic energy vanishing for small momenta. This allows us to freeze slow particles, and hence avoid the re-computation of inter-particle forces, which leads to computational gains. On the other hand, the statistical error may increase since there are a priori more correlations in time. The aim of this work is first to prove the ergodicity of the modified Langevin dynamics (which fails to be hypoelliptic), and next to analyze how the asymptotic variance on ergodic averages depends on the parameters of the modified kinetic energy. Numerical results illustrate the approach, both for low-dimensional systems where we resort to a Galerkin approximation of the generator, and for more realistic systems using Monte Carlo simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Improving dynamical properties of metropolized discretizations of overdamped Langevin dynamics

The discretization of overdamped Langevin dynamics, through schemes such as the Euler-Maruyama method, can be corrected by some acceptance/rejection rule, based on a Metropolis-Hastings criterion for instance. In this case, the invariant measure sampled by the Markov chain is exactly the Boltzmann-Gibbs measure. However, rejections perturb the dynamical consistency of the resulting numerical method with the reference dynamics. We present in this work some modifications of the standard correction of discretizations of overdamped Langevin dynamics on compact spaces by a Metropolis-Hastings procedure, which allow us to either improve the strong order of the numerical method, or to decrease the bias in the estimation of transport coefficients characterizing the effective dynamical behavior of the dynamics. For the latter approach, we rely on modified numerical schemes together with a Barker rule for the acceptance/rejection criterion.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Coupling a reactive potential with a harmonic approximation for atomistic simulations of material failure

Molecular dynamics (MD) simulations involving reactive potentials can be used to model material failure. The empirical potentials which are used in such simulations are able to adapt to the atomic environment, at the expense of a significantly higher computational cost than non-reactive potentials. However, during a simulation of failure, the reactive ability is needed only in some limited parts of the system, where bonds break or form and the atomic environment changes. Therefore, simpler non-reactive potentials can be used in the remainder of the system, provided that such potentials reproduce correctly the behavior of the reactive potentials in this region, and that seamless coupling is ensured at the interface between the reactive and non-reactive regions. In this article, we propose a methodology to combine a reactive potential with a non-reactive approximation thereof, made of a set of harmonic pair and angle interactions and whose parameters are adjusted to predict the same energy, geometry and Hessian in the ground state of the potential. We present a methodology to construct the non-reactive approximation of the reactive potential, and a way to couple these two potentials. We also propose a criterion for on-the-fly substitution of the reactive potential by its non-reactive approximation during a simulation. We illustrate the correctness of this hybrid technique for the case of MD simulation of failure in two-dimensional graphene originally modeled with REBO potential.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Permutation-invariant distance between atomic configurations

We present a permutation-invariant distance between atomic configurations, defined through a functional representation of atomic positions. This distance enables to directly compare different atomic environments with an arbitrary number of particles, without going through a space of reduced dimensionality (i.e. fingerprints) as an intermediate step. Moreover, this distance is naturally invariant through permutations of atoms, avoiding the time consuming associated minimization required by other common criteria (like the Root Mean Square Distance). Finally, the invariance through global rotations is accounted for by a minimization procedure in the space of rotations solved by Monte Carlo simulated annealing. A formal framework is also introduced, showing that the distance we propose verifies the property of a metric on the space of atomic configurations. Two examples of applications are proposed. The first one consists in evaluating faithfulness of some fingerprints (or descriptors), i.e. their capacity to represent the structural information of a configuration. The second application concerns structural analysis, where our distance proves to be efficient in discriminating different local structures and even classifying their degree of similarity.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The computation of averages from equilibrium and nonequilibrium Langevin molecular dynamics

We consider numerical methods for thermodynamic sampling, i.e. computing sequences of points distributed according to the Gibbs-Boltzmann distribution, using Langevin dynamics and overdamped Langevin dynamics (Brownian dynamics). A wide variety of numerical methods for Langevin dynamics may be constructed based on splitting the stochastic differential equations into various component parts, each of which may be propagated exactly in the sense of distributions. Each such method may be viewed as generating samples according to an associated invariant measure that differs from the exact canonical invariant measure by a stepsize-dependent perturbation. We provide error estimates a la Talay-Tubaro on the invariant distribution for small stepsize, and compare the sampling bias obtained for various choices of splitting method. We further investigate the overdamped limit and apply the methods in the context of driven systems where the goal is sampling with respect to a nonequilibrium steady state. Our analyses are illustrated by numerical experiments.

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Nonlinear Fluctuating Hydrodynamics in One Dimension: the Case of Two Conserved Fields

We study the BS model, which is a one-dimensional lattice field theory taking real values. Its dynamics is governed by coupled differential equations plus random nearest neighbor exchanges. The BS model has exactly two locally conserved fields. Through numerical simulations the peak structure of the steady state space-time correlations is determined and compared with nonlinear fluctuating hydrodynamics, which predicts a traveling peak with KPZ scaling function and a standing peak with a scaling function given by the completely asymmetric Levy distribution with parameter $α= 5/3$. As a by-product, we completely classify the universality classes for two coupled stochastic Burgers equations with arbitrary coupling coefficients.

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

Local density dependent potential for compressible mesoparticles

We focus on finding a coarse grained description able to reproduce the thermodynamic behavior of a molecular system by using mesoparticles representing several molecules. Interactions between mesoparticles are modelled by an interparticle potential, and an additional internal equation of state is used to account for the thermic contribution of coarse grained internal degrees of freedom. Moreover, as strong non-equilibrium situations over a wide range of pressure and density are targeted, the internal compressibility of these mesoparticles has to be considered. This is done by introducing a dependence of the potential on the local environment of the mesoparticles, either by defining a spherical local density or by means of a Voronoi tessellation. As an example, a local density dependent potential is fitted to reproduce the Hugoniot curve of a model of nitromethane, where each mesoparticle represents one thousand molecules.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A mathematical formulation of the random phase approximation for crystals

This works extends the recent study on the dielectric permittivity of crystals within the Hartree model [E. Cances and M. Lewin, Arch. Rational Mech. Anal., 197 (2010) 139--177] to the time-dependent setting. In particular, we prove the existence and uniqueness of the nonlinear Hartree dynamics (also called the random phase approximation in the physics literature), in a suitable functional space allowing to describe a local defect embedded in a perfect crystal. We also give a rigorous mathematical definition of the microscopic frequency-dependent polarization matrix, and derive the macroscopic Maxwell-Gauss equation for insulating and semiconducting crystals, from a first order approximation of the nonlinear Hartree model, by means of homogenization arguments.

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Nonequilibrium shear viscosity computations with Langevin dynamics

We study the mathematical properties of a nonequilibrium Langevin dynamics which can be used to estimate the shear viscosity of a system. More precisely, we prove a linear response result which allows to relate averages over the nonequilibrium stationary state of the system to equilibrium canonical expectations. We then write a local conservation law for the average longitudinal velocity of the fluid, and show how, under some closure approximation, the viscosity can be extracted from this profile. We finally characterize the asymptotic behavior of the velocity profile, in the limit where either the transverse or the longitudinal friction go to infinity. Some numerical illustrations of the theoretical results are also presented.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Anomalous diffusion for a class of systems with two conserved quantities

We introduce a class of one dimensional deterministic models of energy-volume conserving interfaces. Numerical simulations show that these dynamics are genuinely super-diffusive. We then modify the dynamics by adding a conservative stochastic noise so that it becomes ergodic. System of conservation laws are derived as hydrodynamic limits of the modified dynamics. Numerical evidence shows these models are still super-diffusive. This is proven rigorously for harmonic potentials.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Free Energy Methods for Bayesian Inference: Efficient Exploration of Univariate Gaussian Mixture Posteriors

Because of their multimodality, mixture posterior distributions are difficult to sample with standard Markov chain Monte Carlo (MCMC) methods. We propose a strategy to enhance the sampling of MCMC in this context, using a biasing procedure which originates from computational Statistical Physics. The principle is first to choose a "reaction coordinate", that is, a "direction" in which the target distribution is multimodal. In a second step, the marginal log-density of the reaction coordinate with respect to the posterior distribution is estimated; minus this quantity is called "free energy" in the computational Statistical Physics literature. To this end, we use adaptive biasing Markov chain algorithms which adapt their targeted invariant distribution on the fly, in order to overcome sampling barriers along the chosen reaction coordinate. Finally, we perform an importance sampling step in order to remove the bias and recover the true posterior. The efficiency factor of the importance sampling step can easily be estimated \emph{a priori} once the bias is known, and appears to be rather large for the test cases we considered. A crucial point is the choice of the reaction coordinate. One standard choice (used for example in the classical Wang-Landau algorithm) is minus the log-posterior density. We discuss other choices. We show in particular that the hyper-parameter that determines the order of magnitude of the variance of each component is both a convenient and an efficient reaction coordinate. We also show how to adapt the method to compute the evidence (marginal likelihood) of a mixture model. We illustrate our approach by analyzing two real data sets.

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

Microscopic calculations of Hugoniot curves of neat TATB and of its detonation products

We compute the Hugoniot curves of both neat TATB and its detonation products mixture using atomistic simulation tools. To compute the Hugoniot states, we adapted our "Sampling Constraints in Average" (SCA) method (Maillet et al., Applied Math. Research eXpress 2008, 2009) to Monte-Carlo simulations. For neat TATB, we show that the potential proposed by Rai (Rai et al., J. Chem. Phys. 129, 2008) is not accurate enough to predict the Hugoniot curve and requires some optimization of its parameters. Concerning detonation products, thermodynamic properties at chemical equilibrium are computed using a specific RxMC method (Bourasseau et al., Phys. Chem. Chem. Phys. 13, 2011) taking into account the presence of carbon clusters in the fluid mixture. We show that this explicit description of the solid phase immersed in the fluid phase modifies the chemical equilibrium.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Negative thermal conductivity of chains of rotors with mechanical forcing

We consider chains of rotors subjected to both thermal and mechanical forcings,in a nonequilibrium steady-state. Unusual nonlinear profiles of temperature and velocities are observed in the system. In particular, the temperature is maximal in the center, which is an indication of the nonlocal behavior of the system. In spite of that, local equilibrium holds for long enough chains. Our numerical results also show that, when the mechanical forcing is strong enough, the energy current can be increased by an inverse temperature gradient. This counterintuitive result again reveals the complexity of nonequilibrium states.

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
preprint2010arXiv

The Microscopic Origin of the Macroscopic Dielectric Permittivity of Crystals: A Mathematical Viewpoint

The purpose of this paper is to provide a mathematical analysis of the Adler-Wiser formula relating the macroscopic relative permittivity tensor to the microscopic structure of the crystal at the atomic level. The technical level of the presentation is kept at its minimum to emphasize the mathematical structure of the results. We also briefly review some models describing the electronic structure of finite systems, focusing on density operator based formulations, as well as the Hartree model for perfect crystals or crystals with a defect.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Constant entropy sampling and release waves of shock compressions

We present several equilibrium methods that allow to compute isentropic processes, either during the compression or the release of the material. These methods are applied to compute the isentropic release of a shocked monoatomic liquid at high pressure and temperature. Moreover, equilibrium results of isentropic release are compared to the direct nonequilibrium simulation of the same process. We show that due to the viscosity of the liquid but also to nonequilibrium effects, the release of the system is not strictly isentropic.

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

Reducing the thermal conductivity of carbon nanotubes below the random isotope limit

We find that introducing segmented isotopic disorder patterns may considerably reduce the thermal conductivity of pristine carbon nanotubes below the uncorrelated disorder value. This is a result of the interplay between different length scales in the phonon scattering process. We use ab-initio atomistic Green's function calculations to quantify the effect of various types of segmentation similar to that experimentally produced by coalescence of isotope-engineered fullerenes.

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
preprint2009arXiv

Thermal conductivity of the Toda lattice with conservative noise

We study the thermal conductivity of the one dimensional Toda lattice perturbed by a stochastic dynamics preserving energy and momentum. The strength of the stochastic noise is controlled by a parameter $γ$. We show that heat transport is anomalous, and that the thermal conductivity diverges with the length $n$ of the chain according to $κ(n) \sim n^α$, with $0 < α\leq 1/2$. In particular, the ballistic heat conduction of the unperturbed Toda chain is destroyed. Besides, the exponent $α$ of the divergence depends on $γ$.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Adiabatic approximation, Gell-Mann and Low theorem and degeneracies: A pedagogical example

We study a simple system described by a 2x2 Hamiltonian and the evolution of the quantum states under the influence of a perturbation. More precisely, when the initial Hamiltonian is not degenerate,we check analytically the validity of the adiabatic approximation and verify that, even if the evolution operator has no limit for adiabatic switchings, the Gell-Mann and Low formula allows to follow the evolution of eigenstates. In the degenerate case, for generic initial eigenstates, the adiabatic approximation (obtained by two different limiting procedures) is either useless or wrong, and the Gell-Mann and Low formula does not hold. We show how to select initial states in order to avoid such failures.

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

0 citations0 reviewsNo code linkedOpen access