Contrôle en chimie quantique : conception et analyse de schémas d'optimisation

Gabriel Turinici

La thèse concerne l'analyse numérique des méthodes de simulation en chimie et physique quantique.

Michel Caffarel, Anthony Scemama

cbenmoha

Nadia Jbili

Cette thèse porte sur des techniques de contrôle optimal pour des systèmes issus de la mécanique quantique et de la résonance magnétique nucléaire. Le travail présenté dans ce mémoire est divisé en quatre parties.Dans la première partie, nous nous sommes intéressés au contrôle optimal simultané de l’équation de Schrödinger dépendante du temps via un champ laser qui représente le contrôle et que nous supposons soumis à une famille de perturbations. Ceci nous conduit à considérer un problème d’optimisation multi-critère via l’introduction d’un ensemble de fonctionnelles de coût à minimiser (au sens de Pareto).Dans la deuxième partie, nous étudions le cadre mathématique de l’équation de Bloch périodique. Les conditions d’optimalité nécessaires du premier ordre sont étudiées. Plus précisément, nous prouvons l’existence d’une solution périodique, ainsi que l’existence d’un optimum.Dans la troisième partie, nous présentons un nouvel algorithme d’optimisation pour les dynamiques périodiques. Cet algorithme est appliqué à la maximisation du signal sur bruit en résonance magnétique nucléaire. Le travail réalisé est ici avant tout numérique et algorithmique. Il s’agit à notre connaissance du premier algorithme de contrôle quantique permettant de considérer des dynamiques périodiques en temps. Nous avons montré l’efficacité de cette méthode pour le cas d’un système de spins homogènes et inhomogènes.La dernier partie permet de présenter l’algorithme de Shinnar-Le-Roux (SLR) qui est une méthode d’optimisation analytique. Des résultats numériques ont été réalisés en comparant cette méthode avec une méthode itérative de type GRAPE introduite dans les chapitres précédents. Le résultat de cette comparaison donne un avantage à l’algorithme SLR.

César Feniou

L'équation de Schrödinger polyélectronique régit l'évolution non relativiste des systèmes quantiques tels que les électrons, atomes et molécules. La chimie quantique ab initio consiste à trouver des solutions approximées de cette équation afin de prédire des propriétés moléculaires comme les longueurs de liaison, les spectres vibrationnels et les mécanismes réactionnels. Malgré les immenses progrès de la puissance de calcul, atteindre la précision chimique (une erreur d'énergie inférieure à 1 kcal/mol) reste exponentiellement coûteux pour les algorithmes classiques lorsque la taille du système augmente. Cette limite a conduit à envisager l'informatique quantique comme alternative pour la simulation moléculaire, les ordinateurs quantiques obéissant aux mêmes principes physiques que les systèmes simulés. Il a été établi que, sous certaines conditions, l'équation de Schrödinger moléculaire pouvait être résolue en temps polynomial sur un ordinateur quantique, permettant en principe une accélération exponentielle par rapport aux méthodes classiques. La concrétisation de cet avantage reste un problème ouvert et requiert des avancées tant matérielles (stabilité des processeurs) que algorithmiques (consommation de ressources). Ce dernier aspect constitue l'objectif principal de cette thèse. Trois défis fondamentaux sont identifiés: établir une correspondance et une discrétisation efficace entre le problème de structure électronique et les qubits; concevoir des algorithmes pour préparer de bonnes approximations de l'état fondamental ; décomposer les opérations quantiques en circuits implémentables à coût minimal. La première partie traite du court terme, où le faible nombre de qubits et le bruit limitent les opérations exécutables. Dans ce contexte, le Variational Quantum Eigensolver (VQE) et son extension adaptative ADAPT-VQE sont des références établies. Nous proposons plusieurs améliorations pour la construction adaptative d'ansätze. L'algorithme Greedy Gradient-Free Adaptive VQE simplifie le processus d'optimisation, réduisant considérablement la consommation de ressources et permettant même une démonstration sur un processeur de 25 qubits. La méthode Overlap-ADAPT-VQE introduit un schéma de compactification variationnelle qui reconstruit des approximations d'états quantiques cibles, permettant en chimie quantique d'accélérer significativement la convergence du VQE pour divers niveaux de corrélation électronique. Nous développons ensuite une stratégie d'initialisation des registres quantiques à partir de fonctions d'onde corrélées issues de calculs classiques, pour combiner calcul haute performance et post-traitement quantique. La seconde partie adopte une perspective de long terme centrée sur le comportement asymptotique des algorithmes. Nous présentons une décomposition de la porte NOT multi-contrôlée avec une profondeur de circuit polylogarithmique en nombre de qubits de contrôle, offrant une amélioration exponentielle sur les approches linéaires existantes. Enfin, nous proposons une méthode de chimie quantique en espace réel, compatible avec les algorithmes quantiques d'extraction de valeurs propres. L'objectif étant d'améliorer la convergence spatiale en remplaçant les discrétisations uniformes, utilisées dans la plupart des approches en première quantification, par des grilles non uniformes, adaptées à la molécule, qui concentrent les points de discrétisation dans les zones de forte densité électronique. Nous étendons cette stratégie aux hamiltoniens transcorrélés : des opérateurs qui régularisent les singularités du potentiel coulombien, et dont les fonctions propres sont lisses dans l'espace. Notre approche est plurielle, explorant tous les niveaux d'abstraction, de la théorie à l'implémentation, et couvrant des algorithmes exploitables sur des machines actuelles, et d'autres destinés aux architectures de long terme. Ces contributions renforcent les fondements des algorithmes quantiques pour la simulation moléculaire.

Jean Hladik, M. Chrysos, Pierre-Emmanuel Hladik +1 more