Search

This article proposes a so-called Quasi Time-Reversible (QTR G-Ext) scheme based on Grassmann extrapolation of density matrices for an accurate calculation of initial guesses in Born-Oppenheimer Molecular Dynamics simulations. The method shows excellent results on four large molecular systems, ranging from 21 to 94 atoms simulated with Kohn-Sham density functional theory surrounded with a classical environment with 6k to 16k atoms. Namely, it clearly reduces the number of self-consistent field iterations, while keeping a similar energy drift as in the extended Lagrangian Born-Oppenheimer method.

The Kohn-Sham method uses a single model system, and corrects it by a density functional the exact user friendly expression of which is not known and is replaced by an approximated, usable, model. We propose to use instead more than one model system, and use a greedy extrapolation method to correct the results of the model systems. Evidently, there is a higher price to pay for it. However, there are also gains: within the same paradigm, e.g., excited states and physical properties can be obtained.

Born-Oppenheimer Molecular Dynamics (BOMD) is a powerful but expensive technique. The main bottleneck in a density functional theory bomd calculation is the solution to the Kohn-Sham (KS) equations, that requires an iterative procedure that starts from a guess for the density matrix. Converged densities from previous points in the trajectory can be used to extrapolate a new guess, however, the non-linear constraint that an idempotent density needs to satisfy make the direct use of standard linear extrapolation techniques not possible. In this contribution, we introduce a locally bijective map between the manifold where the density is defined and its tangent space, so that linear extrapolation can be performed in a vector space while, at the same time, retaining the correct physical properties of the extrapolated density using molecular descriptors. We apply the method to real-life, multiscale polarizable QM/MM.

In this work, we provide the mathematical elements we think essential for a proper understanding of the calculus of the electrostatic energy of point-multipoles of arbitrary order under periodic boundary conditions. The emphasis is put on the expressions of the so-called self parts of the \Ewald\, summation where different expressions can be found in literature. Indeed, such expressions are of prime importance in the context of new generation polarizable force field where the self field appears in the polarization equations. We provide a general framework, where the idea of the \Ewald\ splitting is applied to the electric potential and subsequently, all other quantities such as the electric field, the energy and the forces are derived consistently thereof. Mathematical well-posedness is shown for all these contributions for any order of multipolar distribution.