Development of High Performance Molecular Dynamics : Fast Evaluation of Polarization Forces / Félix Aviat ; sous la direction de Jean-Philip Piquemal et de Alessandra Carbone

Date :

Type : Livre / Book

Type : Thèse / Thesis

Langue / Language : anglais / English

Dynamique moléculaire

Systèmes biologiques

Polarisation (physique nucléaire)

Gradient conjugué, Méthode du

Classification Dewey : 541.22

Piquemal, Jean-Philip (19..-....) (Directeur de thèse / thesis advisor)

Carbone, Alessandra (19..-.... ; mathématicienne) (Directeur de thèse / thesis advisor)

Vuilleumier, Rodolphe (19..-....) (Président du jury de soutenance / praeses)

Dumont, Elise (1980-....) (Rapporteur de la thèse / thesis reporter)

Caffarel, Michel (1961-....) (Rapporteur de la thèse / thesis reporter)

Montes, Matthieu (1980-....) (Membre du jury / opponent)

Sorbonne université (Paris ; 2018-....) (Organisme de soutenance / degree-grantor)

École doctorale Chimie physique et chimie analytique de Paris Centre (Paris ; 2000-....) (Ecole doctorale associée à la thèse / doctoral school)

Laboratoire de chimie théorique (Paris ; 1997-....) (Laboratoire associé à la thèse / thesis associated laboratory)

Résumé / Abstract : La dynamique moléculaire classique est un outil précieux pour explorer le monde de l’infiniment petit, notamment lorsque l’on parle de systèmes biologiques (protéines). Ces simulations sont basées sur des modèles plus ou moins avancés, où la prise en compte des électrons n’est pas évidente. La polarisabilité permet la prise en compte de la mobilité des nuages électronique tout en conservant un cadre classique. Elle peut-être décrite par des dipôles induits, dont le calcul demande la résolution d’un schéma auto-cohérent (SCF), coûteux en terme de temps de calcul et source d’instabilité. 
Dans cette thèse, nous proposons un nouvel algorithme permettant un traitement plus rapide et plus stable des dipôles induits, basé sur la troncature du Gradient Conjugué. La précision, la versatilité, l’efficacité du Gradient Conjugué Tronqué (TCG) sont évalués sur plusieurs système. Son application à des calculs plus sensibles d’énergie libre est également testée. Le TCG est enfin utilisé afin de mettre au point de nouveaux intégrateurs de dynamique moléculaire permettant d’importantes accélérations. On démontre finalement que TCG est un outil polyvalent, adaptable, efficace, qui permet une accélération considérable des dynamiques polarisables longues.

Résumé / Abstract : Classical molecular dynamics is a precious tools to explore the infinitely small world, e.g. when considering biological systems (such as proteins). These simulations are based on physical models of various precision and complexity, where taking electrons into account is not easy. Polarizability allows one to take into account the mobility of the electronic density, while keeping the classical description framework. It can be described using induced dipoles, whose computation is done through a Self-Consistent procedure, which is costly in terms of computational time and can also cause instability. In this thesis, we introduce a new algorithm allowing a faster and more stable treatment of the induced dipoles, based on the truncation of the Conjugate Gradient. Accuracy, versatility, efficiency of the so-called Truncated Conjugate Gradient (TCG) are evaluated on various systems. Its applicability to free energy calculations is also tested. TCG is finally used in order to derive new molecular dynamics integrators allowing for considerable accelerations. To sum this up, TCG proves to be a polyvalent, adaptable, efficient tool, which allows for substantial acceleration of long polarizable dynamics.