Search

Specialized computational chemistry packages have permanently reshaped the landscape of chemical and materials science by providing tools to support and guide experimental efforts and for the prediction of atomistic and electronic properties. In this regard, electronic structure packages have played a special role by using first-principle-driven methodologies to model complex chemical and materials processes. Over the past few decades, the rapid development of computing technologies and the tremendous increase in computational power have offered a unique chance to study complex transformations using sophisticated and predictive many-body techniques that describe correlated behavior of electrons in molecular and condensed phase systems at different levels of theory. In enabling these simulations, novel parallel algorithms have been able to take advantage of computational resources to address the polynomial scaling of electronic structure methods. In this paper, we briefly review the NWChem computational chemistry suite, including its history, design principles, parallel tools, current capabilities, outreach, and outlook.

In the milestoning framework, and more generally in related transition interface sampling schemes, one significantly enhances the calculation of relaxation rates for complex equilibrium kinetics from molecular dynamics simulations between the milestones or interfaces. The goal of the present paper is to advance milestoning applications into the realm of non-equilibrium statistical mechanics, in particular, to calculate entire time correlation functions. In order to accomplish this, we introduce a novel methodology for obtaining the flux through a given milestone configuration as a function of both time and initial configuration and build upon it with a novel formalism describing autocorrelation for Langevin motion in a discrete configuration space. The method is then applied to three different test systems: a harmonic oscillator, which we solve analytically, a two-well potential, which is solved numerically, and an atomistic molecular dynamics simulation of alanine dipeptide.