Study of Liapunov Exponents and the Reversibility of Molecular Dynamics Algorithms
نویسندگان
چکیده
Molecular dynamics algorithms like the Hybrid Monte Carlo [1] and the Kramers equation [2,3] algorithm have been playing a major role in numerical simulations of QCD on a Euclidean lattice. In order for these algorithms to fulfill the detailed balance condition, the classical motion governed by the set of Hamilton’s equations for the system should be reversible. However, as noticed some time ago [4], this reversibility condition is violated due to the round-off errors in the numerical integration of the equations of motion. Recently, it was pointed out [3] that, due to the chaotic nature of the equations of motion, the round-off errors in the integration which violate the reversibility condition get magnified exponentially with a positive Liapunov exponent ν. We will focus our discussion on the equations that arise in the simulations of Wilson QCD with two flavors of quarks with degenerate masses. Periodic boundary conditions have been taken for all the fields. The full partition function for Wilson QCD is given by,
منابع مشابه
Instabilities and Non-Reversibility of Molecular Dynamics Trajectories
The theoretical justification of the Hybrid Monte Carlo algorithm depends upon the molecular dynamics trajectories within it being exactly reversible. If computations were carried out with exact arithmetic then it would be easy to ensure such reversibility, but the use of approximate floating point arithmetic inevitably introduces violations of reversibility. In the absence of evidence to the c...
متن کاملEnergy study at different solvents for potassium Channel Protein by Monte Carlo, Molecular and Langevin Dynamics Simulations
Potassium Channels allow potassium flux and are essential for the generation of electric current acrossexcitable membranes. Potassium Channels are also the targets of various intracellular controlmechanisms; such that the suboptimal regulation of channel function might be related to pathologicalconditions. Realistic studies of ion current in biologic channels present a major challenge for compu...
متن کاملLiapunov Spectra for Infinite Chains of Nonlinear Oscillators
We argue that the spectrum of Liapunov exponents for long chains of nonlinear oscillators, at large energy per mode, may be well approximated by the Liapunov exponents of products of independent random matrices. If, in addition, statistical mechanics applies to the system, the elements of these random matrices have a distribution which may be calculated from the potential and the energy alone. ...
متن کاملAsymptotic Liapunov Exponents Spectrum For An Extended Chaotic Coupled Map Lattice
The scaling hypothesis for the coupled chaotic map lattices (CML) is formulated. Scaling properties of the CML in the regime of extensive chaos observed numerically before is justified analytically. The asymptotic Liapunov exponents spectrum for coupled piece-wise linear chaotic maps is found.
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
عنوان ژورنال:
دوره شماره
صفحات -
تاریخ انتشار 1996