اشترك بالحزمة الذهبية واحصل على وصول غير محدود شمرا أكاديميا
تسجيل مستخدم جديدCanonical sampling through velocity-rescaling
72
0
0.0
(
0
)
اسأل ChatGPT حول البحث
ﻻ يوجد ملخص باللغة العربية
We present a new molecular dynamics algorithm for sampling the canonical distribution. In this approach the velocities of all the particles are rescaled by a properly chosen random factor. The algorithm is formally justified and it is shown that, in spite of its stochastic nature, a quantity can still be defined that remains constant during the evolution. In numerical applications this quantity can be used to measure the accuracy of the sampling. We illustrate the properties of this new method on Lennard-Jones and TIP4P water models in the solid and liquid phases. Its performance is excellent and largely independent on the thermostat parameter also with regard to the dynamic properties.
قيم البحث
اقرأ أيضاً
We present a method for determining the free energy dependence on a selected number of collective variables using an adaptive bias. The formalism provides a unified description which has metadynamics and canonical sampling as limiting cases. Converge
nce and errors can be rigorously and easily controlled. The parameters of the simulation can be tuned so as to focus the computational effort only on the physically relevant regions of the order parameter space. The algorithm is tested on the reconstruction of alanine dipeptide free energy landscape.
We discuss the use of a Langevin equation with a colored (correlated) noise to perform constant-temperature molecular dynamics simulations. Since the equations of motion are linear in nature, it is easy to predict the response of a Hamiltonian system
to such a thermostat and to tune at will the relaxation time of modes of different frequency. This allows one to optimize the time needed to thermalize the system and generate independent configurations. We show how this frequency-dependent response can be exploited to control the temperature of Car-Parrinello-like dynamics, keeping at low temperature the electronic degrees of freedom, without affecting the adiabatic separation from the vibrations of the ions.
In this thesis we have used Quantum Monte Carlo techniques to study two systems that can be regarded as the archetype for neutral strongly interacting systems: 4He, and its fermionic counterpart 3He.More specifically, we have used the Path Integral G
round State and the Path Integral Monte Carlo methods to study a system of two dimensional 3He (2d-3He) and a system of 4He adsorbed on Graphene-Fluoride (GF) and Graphane (GH) at both zero and finite temperature. The purpose of the study of 4He on GF (GH) was the research of new physical phenomena, whereas in the case of 2d-3He it was the application of novel methodologies for the ab-initio study of static and dynamic properties of Fermi systems. In the case of 2d-3He we have computed the spin susceptibility as function of density which turned out to be in very good agreement with experimental data; we have also obtained the first ab-initio evaluation of the zero-sound mode and the dynamic structure factor of 2d-3He that is in remarkably good agreement with experiments. In the case of 4He adsorbed on GF (GH), we determined the zero temperature equilibrium density of the first monolayer of 4He showing also that the commensurate sqrt(3) x sqrt(3) R30 phase is unstable on both substrates; at equilibrium density we found that 4He on GF (GH) is a modulated superfluid with an anisotropic phono-rotonic spectrum; at high coverages we found an incommensurate triangular solid and, on both GF and GH, a commensurate phase at filling factor x= 2/7 that is locally stable or at least metastable. Remarkably, in this commensurate solid phase and for both GF and GH, our computations show preliminary evidence of the presence of a superfluid fraction.
Numerical simulations and finite-size scaling analysis have been carried out to study the percolation behavior of straight rigid rods of length $k$ ($k$-mers) on two-dimensional square lattices. The $k$-mers, containing $k$ identical units (each one
occupying a lattice site), were adsorbed at equilibrium on the lattice. The process was monitored by following the probability $R_{L,k}(theta)$ that a lattice composed of $L times L$ sites percolates at a concentration $theta$ of sites occupied by particles of size $k$. A nonmonotonic size dependence was observed for the percolation threshold, which decreases for small particles sizes, goes through a minimum, and finally asymptotically converges towards a definite value for large segments. This striking behavior has been interpreted as a consequence of the isotropic-nematic phase transition occurring in the system for large values of $k$. Finally, the universality class of the model was found to be the same as for the random percolation model.
Classical molecular dynamics simulations have recently become a standard tool for the study of electrochemical systems. State-of-the-art approaches represent the electrodes as perfect conductors, modelling their responses to the charge distribution o
f electrolytes via the so-called fluctuating charge model. These fluctuating charges are additional degrees of freedom that, in a Born-Oppenheimer spirit, adapt instantaneously to changes in the environment to keep each electrode at a constant potential. Here we show that this model can be treated in the framework of constrained molecular dynamics, leading to a symplectic and time-reversible algorithm for the evolution of all the degrees of freedom of the system. The computational cost and the accuracy of the new method are similar to current alternative implementations of the model. The advantage lies in the accuracy and long term stability guaranteed by the formal properties of the algorithm and in the possibility to systematically introduce additional kinematic conditions of arbitrary number and form. We illustrate the performance of the constrained dynamics approach by enforcing the electroneutrality of the electrodes in a simple capacitor consisting of two graphite electrodes separated by a slab of liquid water.
سجل دخول لتتمكن من نشر تعليقات
التعليقات
جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
