Truncated Conjugate Gradient (TCG): an optimal strategy for the analytical evaluation of the many-body polarization energy and forces in molecular simulations

Abstract : We introduce a new class of methods, denoted as Truncated Conjugate Gradient (TCG) methods, to solve the many-body polarization energy and its associated forces in molecular simulations encountered in molecular dynamics (MD) and Monte-Carlo techniques. The method consists of a fixed number of Conjugate Gradient (CG) iterations. TCG approaches provide a scalable solution to the polarization problem at a user-chosen cost and a corresponding optimal accuracy and complexity. The optimality of the CG-method guarantees that the number of the required matrix-vector products are reduced to a minimum compared to other iterative methods. This family of methods is non empirical, fully adaptive and provides analytical gradients, avoiding therefore any energy drift in MD as compared to popular iterative solvers. Besides speed, one great advantage of this class of approximate methods is that their accuracy is systematically improvable. Indeed, as the CG-method is a Krylov subspace method, the associated error is monotonically reduced at each iteration. On top of that, two improvements can be proposed at virtually no cost: (i) the use of preconditioners can be employed, which leads to the Truncated Preconditioned Conjugate Gradient (TPCG); (ii) since the residual of the final step of the CG-method is available, one additional Picard fi xed point iteration ("peek"), equivalent to one step of Jacobi Over Relaxation (JOR) with relaxation parameter omega, can be made at almost no cost. This method is denoted by TCG-n(omega). Black box adaptive methods to find omega are provided and discussed. Results show that TPCG-3(omega) is converged to high accuracy for various types of systems including proteins and highly charged systems at the fixed cost of 4 matrix-vector products: (3 CG iterations+the initial CG descent direction) whereas T(P)CG-2(omega) provides robust results at a reduced cost (3 matrix-vector products) and off ers new perspectives for long polarizable MD as a production algorithm. The T(P)CG-1(omega) level provides less accurate solutions for inhomogeneous systems, but its applicability to well-conditioned problems such as water is remarkable, with only two matrix-vector product evaluations.
Type de document :
Article dans une revue
Journal of Chemical Theory and Computation, American Chemical Society, 2017, 13 (1), pp.180-190. 〈http://pubs.acs.org/doi/abs/10.1021/acs.jctc.6b00981〉. 〈10.1021/acs.jctc.6b00981〉
Liste complète des métadonnées

Littérature citée [2 références]  Voir  Masquer  Télécharger

https://hal.archives-ouvertes.fr/hal-01395833
Contributeur : Jean-Philip Piquemal <>
Soumis le : jeudi 8 décembre 2016 - 18:43:42
Dernière modification le : jeudi 6 septembre 2018 - 14:58:12
Document(s) archivé(s) le : jeudi 23 mars 2017 - 08:50:58

Fichier

acs.jctc.6b00981.pdf
Publication financée par une institution

Identifiants

Collections

Citation

Félix Aviat, Antoine Levitt, Benjamin Stamm, Yvon Maday, Pengyu Ren, et al.. Truncated Conjugate Gradient (TCG): an optimal strategy for the analytical evaluation of the many-body polarization energy and forces in molecular simulations . Journal of Chemical Theory and Computation, American Chemical Society, 2017, 13 (1), pp.180-190. 〈http://pubs.acs.org/doi/abs/10.1021/acs.jctc.6b00981〉. 〈10.1021/acs.jctc.6b00981〉. 〈hal-01395833v4〉

Partager

Métriques

Consultations de la notice

8109

Téléchargements de fichiers

78