Two of the most successful methods that are presently available for simulating the quantum dynamics of condensed phase systems are centroid molecular dynamics (CMD) and ring polymer molecular dynamics (RPMD). Despite their conceptual differences, practical implementations of these methods differ in just two respects: the choice of the Parrinello-Rahman mass matrix and whether or not a thermostat is applied to the internal modes of the ring polymer during the dynamics. Here we explore a method which is halfway between the two approximations: we keep the path integral bead masses equal to the physical particle masses but attach a Langevin thermostat to the internal modes of the ring polymer during the dynamics. We justify this by showing analytically that the inclusion of an internal mode thermostat does not affect any of the desirable features of RPMD: thermostatted RPMD (TRPMD) is equally valid with respect to everything that has actually been proven about the method as RPMD itself. In particular, because of the choice of bead masses, the resulting method is still optimum in the short-time limit, and the transition state approximation to its reaction rate theory remains closely related to the semiclassical instanton approximation in the deep quantum tunneling regime. In effect, there is a continuous family of methods with these properties, parameterised by the strength of the Langevin friction. Here we explore numerically how the approximation to quantum dynamics depends on this friction, with a particular emphasis on vibrational spectroscopy. We find that a broad range of frictions approaching optimal damping give similar results, and that these results are immune to both the resonance problem of RPMD and the curvature problem of CMD.

We present a comparative assessment of the accuracy of two different approaches for evaluating dispersion interactions: interatomic pairwise corrections and semiempirical meta-generalized-gradient-approximation (meta-GGA)-based functionals. This is achieved by employing conventional (semi)local and (screened-)hybrid functionals, as well as semiempirical hybrid and nonhybrid meta-GGA functionals of the M06 family, with and without interatomic pairwise Tkatchenko-Scheffler corrections. All of those are tested against the benchmark S22 set of weakly bound systems, a representative larger molecular complex (dimer of NiPc molecules), and a representative dispersively bound solid (hexagonal boron nitride). For the S22 database, we also compare our results with those obtained from the pairwise correction of Grimme (DFT-D3) and nonlocal Langreth-Lundqvist functionals (vdW-DF1 and vdW-DF2). We find that the semiempirical kinetic-energy-density dependence introduced in the M06 functionals mimics some of the nonlocal correlation needed to describe dispersion. However, long-range contributions are still missing. Pair-wise interatomic corrections, applied to conventional semilocal or hybrid functionals, or to M06 functionals, provide for a satisfactory level of accuracy irrespectively of the underlying functional. Specifically, screened-hybrid functionals such as the Heyd-Scuseria-Ernzerhof (HSE) approach reduce self-interaction errors in systems possessing both localized and delocalized orbitals and can be applied to both finite and extended systems. Therefore, they serve as a useful underlying functional for dispersion corrections.

Progress in the atomic-scale modelling of matter over the past decade has been tremendous. This progress has been brought about by improvements in methods for evaluating interatomic forces that work by either solving the electronic structure problem explicitly, or by computing accurate approximations of the solution and by the development of techniques that use the Born-Oppenheimer (BO) forces to move the atoms on the BO potential energy surface. As a consequence of these developments it is now possible to identify stable or metastable states, to sample configurations consistent with the appropriate thermodynamic ensemble, and to estimate the kinetics of reactions and phase transitions. All too often, however, progress is slowed down by the bottleneck associated with implementing new optimization algorithms and/or sampling techniques into the many existing electronic-structure and empirical-potential codes. To address this problem, we are thus releasing a new version of the i-PI software. This piece of software is an easily extensible framework for implementing advanced atomistic simulation techniques using interatomic potentials and forces calculated by an external driver code. While the original version of the code[1] was developed with a focus on path integral molecular dynamics techniques, this second release of i-PI not only includes several new advanced path integral methods, but also offers other classes of algorithms. In other words, i-PI is moving towards becoming a universal force engine that is both modular and tightly coupled to the driver codes that evaluate the potential energy surface and its derivatives.

The properties of hydrogen ions in aqueous solution are governed by the ability of water to incorporate ions in a dynamical hydrogen bond network, characterized by a structural variability that has complicated the development of a consistent molecular level description of H(+)(aq). Isolated protonated water clusters, H(+)(H2O)n, serve as finite model systems for H(+)(aq), which are amenable to highly sensitive and selective gas phase spectroscopic techniques. Here, we isolate and assign the infrared (IR) signatures of the Zundel-type and Eigen-type isomers of H(+)(H2O)6, the smallest protonated water cluster for which both of these characteristic binding motifs coexist, down into the terahertz spectral region. We use isomer-selective double-resonance population labeling spectroscopy on messenger-tagged H(+)(H2O)6·H2 complexes from 260 to 3900 cm(-1). Ab initio molecular dynamics calculations qualitatively recover the IR spectra of the two isomers and allow attributing the increased width of IR bands associated with H-bonded moieties to anharmonicities rather than excited state lifetime broadening. Characteristic hydrogen-bond stretching bands are observed below 400 cm(-1).

The polyalanine-based peptide series Ac-Alan-LysH + (n=5-20) is a prime example that a secondary structure motif which is well-known from the solution phase (here: helices) can be formed in vacuo. We here revisit this conclusion for n=5,10,15, using density-functional theory (van der Waals corrected generalized gradient approximation), and gas-phase infrared vibrational spectroscopy. For the longer molecules (n=10,15) α-helical models provide good qualitative agreement (theory vs. experiment) already in the harmonic approximation. For n=5, the lowest energy conformer is not a simple helix, but competes closely with α-helical motifs at 300 K. Close agreement between infrared spectra from experiment and ab initio molecular dynamics (including anharmonic effects) supports our findings.