Search

Molecular dynamics (MD) simulations are used in biochemistry, physics, and other fields to study the motions, thermodynamic properties, and the interactions between molecules. Computational limitations and the complexity of these problems, however, create the need for approximations to the standard MD methods and for uncertainty quantification and reliability assessment of those approximations. In this paper, we exploit the intrinsic two-scale nature ofMDto construct a class of large-scale dynamics approximations. The reliability of these methods is evaluated here by measuring the differences between full, classical MD simulations and those based on these large-scale approximations. Molecular dynamics evolutions are non-linear and chaotic, so the complete details of molecular evolutions cannot be accurately predicted even using full, classical MD simulations. This paper provides numerical results that demonstrate the existence of computationally efficient large-scale MD approximations which accurately model certain large-scale properties of the molecules: the energy, the linear and angular momenta, and other macroscopic features of molecular motions. [ABSTRACT FROM AUTHOR]

Most of today's molecular-dynamics simulations of materials are based on the Born-Oppenheimer approximation. There are many cases, however, in which the coupling of the electrons and nuclei is important and it is necessary to go beyond the Born-Oppenheimer approximation. In these methods, the non-adiabatic coupling vectors are fundamental since they represent the link between the classical atomic motion of the nuclei and the time evolution of the quantum electronic state. In this paper we analyze the calculation of non-adiabatic coupling vectors in a basis set of local orbitals and derive an expression to calculate them in a practical and computationally efficient way. Some examples of the application of this expression using a local-orbital density functional theory approach are presented for a few simple molecules: H3, formaldimine, and azobenzene. These results show that the approach presented here, using the Slater transition-state density, is a very promising way for the practical calculation of non-adiabatic coupling vectors for large systems. [ABSTRACT FROM AUTHOR]

The multiscale coarse-graining (MS-CG) method is a method for constructing a coarse-grained (CG) model of a system using data obtained from molecular dynamics simulations of the corresponding atomically detailed model. The formal statistical mechanical derivation of the method shows that the potential energy function extracted from an MS-CG calculation is a variational approximation for the true potential of mean force of the CG sites, one that becomes exact in the limit that a complete basis set is used in the variational calculation if enough data are obtained from the atomistic simulations. Most applications of the MS-CG method have employed a representation for the nonbonded part of the CG potential that is a sum of all possible pair interactions. This approach, despite being quite successful for some CG models, is inadequate for some others. Here we propose a systematic method for including three body terms as well as two body terms in the nonbonded part of the CG potential energy. The current method is more general than a previous version presented in a recent paper of this series [L. Larini, L. Lu, and G. A. Voth, J. Chem. Phys. 132, 164107 (2010)], in the sense that it does not make any restrictive choices for the functional form of the three body potential. We use hierarchical multiresolution functions that are similar to wavelets to develop very flexible basis function expansions with both two and three body basis functions. The variational problem is solved by a numerical technique that is capable of automatically selecting an appropriate subset of basis functions from a large initial set. We apply the method to two very different coarse-grained models: a solvent free model of a two component solution made of identical Lennard-Jones particles and a one site model of SPC/E water where a site is placed at the center of mass of each water molecule. These calculations show that the inclusion of three body terms in the nonbonded CG potential can lead to significant improvement in the accuracy of CG potentials and hence of CG simulations. [ABSTRACT FROM AUTHOR]

Stochastic simulation of coupled chemical reactions is often computationally intensive, especially if a chemical system contains reactions occurring on different time scales. In this paper, we introduce a multiscale methodology suitable to address this problem, assuming that the evolution of the slow species in the system is well approximated by a Langevin process. It is based on the conditional stochastic simulation algorithm (CSSA) which samples from the conditional distribution of the suitably defined fast variables, given values for the slow variables. In the constrained multiscale algorithm (CMA) a single realization of the CSSA is then used for each value of the slow variable to approximate the effective drift and diffusion terms, in a similar manner to the constrained mean-force computations in other applications such as molecular dynamics. We then show how using the ensuing Fokker-Planck equation approximation, we can in turn approximate average switching times in stochastic chemical systems. [ABSTRACT FROM AUTHOR]

An efficient time-dependent importance sampling method is developed for the Monte Carlo calculation of time correlation functions via the initial value representation (IVR) of semiclassical (SC) theory. A prefactor-free time-dependent sampling function weights the importance of a trajectory based on the magnitude of its contribution to the time correlation function, and global trial moves are used to facilitate the efficient sampling the phase space of initial conditions. The method can be generally applied to sampling rare events efficiently while avoiding being trapped in a local region of the phase space. Results presented in the paper for two system-bath models demonstrate the efficiency of this new importance sampling method for full SC-IVR calculations. [ABSTRACT FROM AUTHOR]

Trotter's theorem forms the theoretical basis of most modern molecular dynamics. In essence this theorem states that a time displacement operator (a Lie operator) constructed by exponentiating a sum of noncommuting operators can be approximated by a product of single operators provided the time interval is 'very small.' In theory 'very small' implies infinitesimally small (at which point the approximate product becomes exact), while in practical analysis a finite time interval is divided into several small subintervals or steps. It follows, therefore, that the larger the number of steps the better the approximation to the exact time displacement operator. The question therefore arises: How many steps are sufficient? For bounded operators, standard theorems are available to provide the answer. In this paper we show that a very simple combinatorial formula can be derived which allows the computation of the global differences (as a function of the number of steps) between the Taylor coefficients of the exact time displacement operator and an approximate one constructed by using a finite number of steps. The formula holds for both bounded and nonbounded operators and shows, quantitatively, what is qualitatively expected-that the error decreases with increasing number of steps. Furthermore, the formula applies irrespective of the complexity of the system, boundary conditions, or the thermodynamic ensemble employed for averaging the initial conditions. The analysis yields explicit expressions for the Taylor coefficients which are then used to compute the errors. In the case of the algorithmically based practical numerical simulations in which fixed, albeit small, steps are repeatedly applied, the rise in the number of steps does not reduce the size of the steps but increases the total time of interest. The combinatorial formula shows that, here, the errors diverge. Furthermore, this work can be used to supplement other efforts such as the use of shadow Hamiltonians where the truncation of the series expansion of the latter will produce errors in the higher order propagator moments. [ABSTRACT FROM AUTHOR]