Giovanni Bussi

Canonical sampling through velocity rescaling

The authors 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. The authors 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 of the thermostat parameter also with regard to the dynamic properties.

PLUMED: a portable plugin for free-energy calculations with molecular dynamics

Here we present a program aimed at free-energy calculations in molecular systems. It consists of a series of routines that can be interfaced with the most popular classical molecular dynamics (MD) codes through a simple patching procedure. This leaves the possibility for the user to exploit many different MD engines depending on the system simulated and on the computational resources available. Free-energy calculations can be performed as a function of many collective variables, with a particular focus on biological problems, and using state-of-the-art methods such as metadynamics, umbrella sampling and Jarzynski-equation based steered MD. The present software, written in ANSI-C language, can be easily interfaced with both Fortran and C/C++ codes.

Well-Tempered Metadynamics : A Smoothly Converging and Tunable Free-Energy Method

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. Convergence 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 an alanine dipeptide free-energy landscape.

PLUMED 2: New feathers for an old bird

Enhancing sampling and analyzing simulations are central issues in molecular simulation. Recently, we introduced PLUMED, an open-source plug-in that provides some of the most popular molecular dynamics (MD) codes with implementations of a variety of different enhanced sampling algorithms and collective variables (CVs). The rapid changes in this field, in particular new directions in enhanced sampling and dimensionality reduction together with new hardware, require a code that is more flexible and more efficient. We therefore present PLUMED 2 here—a complete rewrite of the code in an object-oriented programming language (C++). This new version introduces greater flexibility and greater modularity, which both extends its core capabilities and makes it far easier to add new methods and CVs. It also has a simpler interface with the MD engines and provides a single software library containing both tools and core facilities. Ultimately, the new code better serves the ever-growing community of users and contributors in coping with the new challenges arising in the field.

Equilibrium Free Energies from Nonequilibrium Metadynamics

In this Letter we propose a new formalism to map history-dependent metadynamics in a Markovian process. We apply this formalism to model Langevin dynamics and determine the equilibrium distribution of a collection of simulations. We demonstrate that the reconstructed free energy is an unbiased estimate of the underlying free energy and analytically derive an expression for the error. The present results can be applied to other history-dependent stochastic processes, such as Wang-Landau sampling.

Excitons in Carbon Nanotubes: An Ab Initio Symmetry-Based Approach

The optical absorption spectrum of the carbon (4,2) nanotube is computed using an ab initio many-body approach which takes into account excitonic effects. We develop a new method involving a local basis set which is symmetric with respect to the screw-symmetry of the tube. Such a method has the advantages of scaling faster than plane-wave methods and allowing for a precise determination of the symmetry character of the single-particle states, two-particle excitations, and selection rules. The binding energy of the lowest, optically active states is approximately 0.8 eV. The corresponding exciton wave functions are delocalized along the circumference of the tube and localized in the direction of the tube axis.

Enhanced Sampling in Molecular Dynamics Using Metadynamics, Replica-Exchange, and Temperature-Acceleration

We review a selection of methods for performing enhanced sampling in 1 molecular dynamics simulations. We consider methods based on collective variable biasing 2

Isothermal-isobaric molecular dynamics using stochastic velocity rescaling

The authors present a new molecular dynamics algorithm for sampling the isothermal-isobaric ensemble. In this approach the velocities of all particles and volume degrees of freedom are rescaled by a properly chosen random factor. The technical aspects concerning the derivation of the integration scheme and the conservation laws are discussed in detail. The efficiency of the barostat is examined in Lennard-Jones solid and liquid near the triple point and compared to the deterministic Nose-Hoover and the stochastic Langevin methods. In particular, the dependence of the sampling efficiency on the choice of the thermostat and barostat relaxation times is systematically analyzed.

Accurate sampling using Langevin dynamics.

We show how to derive a simple integrator for the Langevin equation and illustrate how it is possible to check the accuracy of the obtained distribution on the fly, using the concept of effective energy introduced in a recent paper [J. Chem. Phys. 126, 014101 (2007)]. Our integrator leads to correct sampling also in the difficult high-friction limit. We also show how these ideas can be applied in practical simulations, using a Lennard-Jones crystal as a paradigmatic case.

Solid State Effects on Exciton States and Optical Properties of PPV

We perform ab initio calculations of optical properties for a typical semiconductor conjugated polymer, poly-para-phenylenevinylene, in both isolated chain and crystalline packing. In order to obtain results for excitonic energies and real-space wave functions we explicitly include electron-hole interaction within the density-matrix formalism. We find that the details of crystalline arrangement crucially affect the optical properties, leading to a richer exciton structure and opening nonradiative decay channels. This has implications for the optical activity and optoelectronic applications of polymer films.