Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userReconstructing the Density of States by History-Dependent Metadynamics
52
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We present a novel method for the calculation of the energy density of states D(E) for systems described by classical statistical mechanics. The method builds on an extension of a recently proposed strategy that allows the free energy profile of a canonical system to be recovered within a pre-assigned accuracy,[A. Laio and M. Parrinello, PNAS 2002]. The method allows a good control over the error on the recovered system entropy. This fact is exploited to obtain D(E) more efficiently by combining measurements at different temperatures. The accuracy and efficiency of the method are tested for the two-dimensional Ising model (up to size 50x50) by comparison with both exact results and previous studies. This method is a general one and should be applicable to more realistic model systems.
rate research
Read More
Sampling complex potential energies is one of the most pressing challenges of contemporary computational science. Inspired by recent efforts that use quantum effects and discretized Feynmans path integrals to overcome large barriers we propose a replica exchange method. In each replica two copies of the same system with halved potential strengths interact via inelastic springs. The strength of the spring is varied in the different replicas so as to bridge the gap between the infinitely strong spring, that corresponds to the Boltzmann replica and the less tight ones. We enhance the spring length fluctuations using Metadynamics. We test the method on simple yet challenging problems.
It is shown that the algorithm introduced in [1] and conceived to deal with continuous degrees of freedom models is well suited to compute the density of states in models with a discrete energy spectrum too. The q=10 D=2 Potts model is considered as a test case, and it is shown that using the Maxwell construction the interface free energy can be obtained, in the thermodynamic limit, with a good degree of accuracy.
For a classical system of noninteracting particles we establish recursive integral equations for the density of states on the microcanonical ensemble. The recursion can be either on the number of particles or on the dimension of the system. The solution of the integral equations is particularly simple when the single-particle density of states in one dimension follows a power law. Otherwise it can be obtained using a Laplace transform method. Since the Laplace transform of the microcanonical density of states is the canonical partition function, it factorizes for a system of noninteracting particles and the solution of the problem is straightforward. The results are illustrated on several classical examples.
In this paper we propose a new formalism to map history-dependent metadynamics in a Markovian process. We apply this formalism to a model Langevin dynamics and determine the equilibrium distribution of a collection of simulations. We demonstrate that the reconstructed free energy is an unbiased estimate of the underlying free energy and analytically derive an expression for the error. The present results can be applied to other history-dependent stochastic processes such as Wang-Landau sampling.
Since the Time-Dependent Density Functional Theory is mathematically formulated through non-linear coupled time-dependent 3-dimensional partial differential equations it is natural to expect a strong sensitivity of its solutions to variations of the initial conditions, akin to the butterfly effect ubiquitous in classical dynamics. Since the Schrodinger equation for an interacting many-body system is however linear and (mathematically) the exact equations of the Density Functional Theory reproduce the corresponding one-body properties, it would follow that the Lyapunov exponents are also vanishing within a Density Functional Theory framework. Whether for realistic implementations of the Time-Dependent Density Functional Theory the question of absence of the butterfly effect and whether the dynamics provided is indeed a predictable theory was never discussed. At the same time, since the time-dependent density functional theory is a unique tool allowing us the study of non-equilibrium dynamics of strongly interacting many-fermion systems, the question of predictability of this theoretical framework is of paramount importance. Our analysis, for a number of quantum superfluid any-body systems (unitary Fermi gas, nuclear fission, and heavy-ion collisions) with a classical equivalent number of degrees of freedom ${cal O}(10^{10})$ and larger, suggests that its maximum Lyapunov are negligible for all practical purposes.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
