Simulation numérique de systèmes quantiques ouverts en grande dimension

José Lages

Julien Roussel

Resoudre exactement l'equation de Lindblad, qui modelise l'evolution des systemes quantiques ouverts, n'est possible que pour des tailles petites. C'est pourquoi on souhaite disposer d'une gamme de methodes permettant de calculer des solutions approchees. On presentera une methode deterministe dite de rang reduit, qui repose sur une forme reduite approchee de la matrice densite, puis une methode stochastique de Monte-Carlo qui donne la solution a une erreur statistique pres. On exposera d'une part une extension de la methode de rang reduit qui permet d'utiliser un rang adaptatif en majorant l'erreur commise, et d'autre part une nouvelle methode de reduction de variance. Cette nouvelle methode permet d'accelerer considerablement les simulations de Monte-Carlo en utilisant la solution approchee fournie par la methode rang reduit comme variable de controle.

Fatma Hamzaoui, Besma Othmani, Roselyne Chotin-Avot +3 more

Mazyar Mirrahimi

Dans cette thèse, nous étudions trois classes de modèles utilisés dans la littérature pour représenter les systèmes quantiques: 1 -L'équation de Schrödinger où le contrôle agit sur le système de façon bilinéaire; 2 -L'équation de Lindblad; 3 -Les filtres quantiques (modèles stochastiques). Les contributions de la thèse concernant l'équation de Schrödinger se répartissent en trois parties. Dans le premier chapitre, nous étudions la contrôlabilité d'un tel système. Le cas de dimension finie étant déjà bien exploré, nous traitons l'exemple d'un oscillateur harmonique quantique comme un cas typique des problèmes de dimension infinie. Parmi les résultats obtenus nous retrouvons transposée dans les termes de la théorie du contrôle, l'assertion bien connue des physiciens: ``les sources classiques de contrôle ne peuvent générer que de la lumière classique''. La question de la génération des trajectoires est abordée dans le Chapitre 2. Le contrôle en boucle ouverte du système est alors traité à l'aide des méthodes de stabilisation de Lyapounov. Ces méthodes de contrôle par feedback sont utilisées en simulation et le contrôle retrouvé est ensuite inséré en boucle ouverte dans le système physique. La convergence est étudiée dans différentes configurations et des exemples numériques tirés de la chimie quantique sont testés. Enfin dans le chapitre 3, nous étudions le problème inverse d'identification de l'Hamiltonien. Malgré le grand intérêt pratique que présente ce problème, peu de contributions ont été apportées jusqu'à maintenant. Nous étudions d'abord le problème mathématique d'identifiabilité. Une première réponse positive à cette question est apportée. Ensuite nous considérons le problème d'identification. A l'aide de méthodes numériques d'optimisation, nous proposons une première approche qui permet de résoudre ce problème inverse. Au sujet de l'équation de Lindblad, la contribution de cette thése se résume à la réduction du modèle lorsque certaines hypothèses sur les durées de vie atomiques sont vérifiées. Cette étude peut être considérée comme une première étape vers le contrôle en boucle fermée d'un ensemble statistique de systèmes quantiques. Finalement dans le chapitre 5, nous considérons les filtres quantiques. Certaines méthodes issues de la théorie des probabilités ainsi que les techniques de Lyapounov stochastiques nous permettent d'étudier la stabilisation globale de ces modéles.

Kaelan Donatella

Dynamique hors-équilibre des systèmes quantiques ouverts et fermés Cette thèse est consacrée à l'étude de la dynamique de hors-équilibre dans les systèmes quantiques fermés et ouverts. La question de l'obtention d'équations dynamiques numériquement efficaces pour simuler de tels systèmes, qui souffrent d'une malédiction de dimensionnalité lorsqu'ils sont de grande taille est au centre de la thèse, ainsi que les divers nouveaux résultats physiques qui peuvent être obtenus grâce à ces équations efficaces. La question de l'obtention de représentations efficaces des états quantiques est également centrale, à travers la construction d'Ansatz. En particulier, nous présentons la méthode du sous-espace dynamique, une nouvelle méthode numérique pour la simulation de l'évolution temporelle de systèmes quantiques ouverts avec une entropie modérée. En représentant la matrice densité par un sous-espace constitué de M états purs, la complexité de l'intégration de l'équation maîtresse de Lindblad est considérablement réduite, et l'algorithme obtenu est entièrement déterministe. Nous appliquons ensuite cette méthode pour modéliser des algorithmes quantiques bruités, en simulant un algorithme quantique bruité de transformation de Fourier qui est au centre de l'algorithme de Shor pour la factorisation de nombre premiers. Nous trouvons une loi d'échelle pour l'erreur en fonction du taux de dissipation, et le comportement de l'infidélité par rapport au circuit idéal est modélisé en fonction des états initiaux qui sont fournis en entrée à l'algorithme. Ensuite, l'application des méthodes d'apprentissage automatique à la mécanique quantique est présentée, à travers la classe d'Ansatz des états quantiques neuronaux (NQS). La construction générale des NQS est introduite, ainsi que sa relation avec d'autres Ansatz tels que les réseaux tensoriels. Ensuite, le problème de la simulation de la dynamique des systèmes fermés avec NQS est étudié, ainsi que les problèmes courants liés à la géométrie de l'espace des paramètres qui peuvent survenir. Un schéma alternatif est introduit et appliqué aux quenches dépendantes du temps dans les systèmes critiques, où le mécanisme de Kibble-Zurek peut être étudié. Enfin, différentes constructions d'Ansatz pour les systèmes ouverts sont présentées, notamment la combinaison de la méthode du sous-espace dynamique avec un Ansatz neuronal. Enfin, nous étudions théoriquement la propagation de l'intrication dans les systèmes photoniques dissipatifs à forte interaction, motivés par la réalisation expérimentale récente de telles plateformes. Nos résultats montrent qu'en dépit du caractère dissipatif du système étudié, la propagation de l'intrication est ballistique avec une vitesse en accord avec l'image de quasiparticules impliquées dans la dynamique, des doublons et des holons photoniques. Notre analyse révèle que la dissipation est fortement asymétrique dans les deux cas étudiés de l'injection d'un doublon ou de son extraction, avec une influence plus forte dans le cas de l'extraction.

Michel Le Bellac