Fourier Acceleration of Langevin Molecular Dynamics

Fourier acceleration has been successfully applied to the simulation of lattice field theories for more than a decade. In this paper, we extend the method to the dynamics of discrete particles moving in continuum. Although our method is based on a mapping of the particles' dynamics to a regular grid so that discrete Fourier transforms may be taken, it should be emphasized that the introduction of the grid is a purely algorithmic device and that no smoothing, coarse-graining or mean-field approximations are made. The method thus can be applied to the equations of motion of molecular dynamics (MD), or its Langevin or Brownian variants. For example, in Langevin MD simulations our acceleration technique permits a straightforward spectral decomposition of forces so that the long-wavelength modes are integrated with a longer time step, thereby reducing the time required to reach equilibrium or to decorrelate the system in equilibrium. Speedup factors of up to 30 are observed relative to pure (unaccelerated) Langevin MD. As with acceleration of critical lattice models, even further gains relative to the unaccelerated method are expected for larger systems. Preliminary results for Fourier-accelerated molecular dynamics are presented in order to illustrate the basic concepts. Possible extensions of the method and further lines of research are discussed.Comment: 11 pages, two illustrations included using graphic

Stochastic resonance-free multiple time-step algorithm for molecular dynamics with very large time steps

Molecular dynamics is one of the most commonly used approaches for studying the dynamics and statistical distributions of many physical, chemical, and biological systems using atomistic or coarse-grained models. It is often the case, however, that the interparticle forces drive motion on many time scales, and the efficiency of a calculation is limited by the choice of time step, which must be sufficiently small that the fastest force components are accurately integrated. Multiple time-stepping algorithms partially alleviate this inefficiency by assigning to each time scale an appropriately chosen step-size. However, such approaches are limited by resonance phenomena, wherein motion on the fastest time scales limits the step sizes associated with slower time scales. In atomistic models of biomolecular systems, for example, resonances limit the largest time step to around 5-6 fs. In this paper, we introduce a set of stochastic isokinetic equations of motion that are shown to be rigorously ergodic and that can be integrated using a multiple time-stepping algorithm that can be easily implemented in existing molecular dynamics codes. The technique is applied to a simple, illustrative problem and then to a more realistic system, namely, a flexible water model. Using this approach outer time steps as large as 100 fs are shown to be possible

Comparison of two non-primitive methods for path integral simulations: Higher-order corrections vs. an effective propagator approach

Two methods are compared that are used in path integral simulations. Both methods aim to achieve faster convergence to the quantum limit than the so-called primitive algorithm (PA). One method, originally proposed by Takahashi and Imada, is based on a higher-order approximation (HOA) of the quantum mechanical density operator. The other method is based upon an effective propagator (EPr). This propagator is constructed such that it produces correctly one and two-particle imaginary time correlation functions in the limit of small densities even for finite Trotter numbers P. We discuss the conceptual differences between both methods and compare the convergence rate of both approaches. While the HOA method converges faster than the EPr approach, EPr gives surprisingly good estimates of thermal quantities already for P = 1. Despite a significant improvement with respect to PA, neither HOA nor EPr overcomes the need to increase P linearly with inverse temperature. We also derive the proper estimator for radial distribution functions for HOA based path integral simulations.Comment: 17 pages, latex, 6 postscript figure

Network development in biological gels: role in lymphatic vessel development

In this paper, we present a model that explains the prepatterning of lymphatic vessel morphology in collagen gels. This model is derived using the theory of two phase rubber material due to Flory and coworkers and it consists of two coupled fourth order partial differential equations describing the evolution of the collagen volume fraction, and the evolution of the proton concentration in a collagen implant; as described in experiments of Boardman and Swartz (Circ. Res. 92, 801–808, 2003). Using linear stability analysis, we find that above a critical level of proton concentration, spatial patterns form due to small perturbations in the initially uniform steady state. Using a long wavelength reduction, we can reduce the two coupled partial differential equations to one fourth order equation that is very similar to the Cahn–Hilliard equation; however, it has more complex nonlinearities and degeneracies. We present the results of numerical simulations and discuss the biological implications of our model

Some aspects of the Liouville equation in mathematical physics and statistical mechanics

This paper presents some mathematical aspects of Classical Liouville theorem and we have noted some mathematical theorems about its initial value problem. Furthermore, we have implied on the formal frame work of Stochastic Liouville equation (SLE)

Gibbs' Paradox according to Gibbs and slightly beyond

The so-called Gibbs paradox is a paradigmatic narrative illustrating the necessity to account for the N! ways of permuting N identical particles when summing over microstates. Yet, there exist some mixing scenarios for which the expected thermodynamic outcome depends on the viewpoint one chooses to justify this combinatorial term. After a brief summary on Gibbs' paradox and what is the standard rationale used to justify its resolution, we will allow ourself to question from a historical standpoint whether the Gibbs paradox has actually anything to do with Gibbs' work. In so doing, we also aim at shedding a new light with regards to some of the theoretical claims surrounding its resolution. We will then turn to the statistical thermodynamics of discrete and continuous mixtures and introduce the notion of composition entropy to characterise these systems. This will enable us to address, in a certain sense, a "curiosity" pointed out by Gibbs in a paper published in 1876. Finally, we will �nish by proposing a connexion between the results we propose and a recent extension of the Landauer bound regarding the minimum amount of heat to be dissipated to reset one bit of memory

Constrained molecular dynamics in the isothermal-isobaric ensemble and its adaptation for adiabatic free energy dynamics

The implementation of holonomic constraints within measure-preserving integrators for molecular dynamics simulations in the isothermal-isobaric ensemble is considered. We review the basic methodology of generating measure-preserving integrators for the microcanonical, canonical, and isothermal-isobaric ensembles and proceed to show how the standard SHAKE and RATTLE algorithms must be modified for the isothermal-isobaric ensemble. Comparison is made between constrained and unconstrained simulations employing multiple time scale integration techniques. Finally, we describe a temperature accelerated version of the isothermal-isobaric molecular dynamics approach, in which the cell matrix is adiabatically decoupled from the particles and maintained at a high temperature as a means of exploring polymorphism in molecular crystals. We demonstrate that constraints can be easily adapted for this new approach and, again, we compare the performace of this temperature-accelerated scheme with and without bond constraints

Endpoint-restricted adiabatic free energy dynamics approach for the exploration of biomolecular conformational equilibria.

A method for calculating the free energy difference between two structurally defined conformational states of a chemical system is developed. A path is defined using a previously reported collective variable that interpolates between two or more conformations, and a restraint is introduced in order to keep the system close to the path. The evolution of the system along the path, which typically presents a high free energy barrier, is generated using enhanced sampling schemes. Although the formulation of the method in terms of a path is quite general, an important advance in this work is the demonstration that prior knowledge of the path is, in fact, not needed and that the free energy difference can be obtained using a simplified definition of the path collective variable that <i>only</i> involves the endpoints. We first validate this method on cyclohexane isomerization. The method is then tested for an extensive conformational change in a realistic molecular system by calculating the free energy difference between the <i>α</i> -helix and <i>β</i> -hairpin conformations of deca-alanine in solution. Finally, the method is applied to a biologically relevant system to calculate the free energy difference of an observed and a hypothetical conformation of an antigenic peptide bound to a major histocompatibility complex