Variational non-gaussian approach to interacting spin-boson models

Résumé Technique : Approche Variationnelle Non-Gaussienne pour les Modèles Spin-Boson Interactifs

Énoncé du Problème
Les modèles spin-boson interactifs, où des degrés de liberté quantiques discrets (spins) se couplent à un ou plusieurs modes bosoniques (photons ou phonons), sont au cœur de la science quantique contemporaine, incluant l'électrodynamique quantique de cavité (QED de cavité), la QED de circuit et les systèmes spin-phonon. Un défi primaire de la simulation de ces systèmes provient de la nature de dimension infinie de l'espace de Hilbert bosonique combinée à de fortes corrélations à plusieurs corps dans le secteur de spin. Les méthodes numériques conventionnelles font face à un compromis : la troncature de la base d'occupation bosonique introduit un biais, tandis que l'assouplissement du secteur de spin vers des descriptions de champ moyen échoue à capturer l'intrication essentielle spin-boson. De plus, les approches directes par réseaux de tenseurs (comme la DMRG) appliquées au système composite complet souffrent souvent d'une croissance rapide de la dimension de liaison en raison de la nécessité de représenter explicitement le champ bosonique aux côtés des spins.

Méthodologie
Les auteurs proposent un cadre variationnel hybride qui découple l'optimisation du secteur bosonique du secteur de spin. L'état fondamental est représenté au sein d'un variational manifold non-gaussien compact défini par l'ansatz :
$$|\psi_{\text{NGS}}\rangle = U_\lambda (U_{\text{GS}} |0_b\rangle \otimes |\phi_s\rangle)$$
où :

1. Couche Bosonique Gaussienne ($U_{\text{GS}}$) : Une transformation gaussienne comprenant des opérateurs de déplacement ($\Delta_R$) et de compression ($\xi$) agissant sur le vide bosonique $|0_b\rangle$. Cela capture les fluctuations bosoniques et les décalages de champ moyen.
2. Transformation de Revêtement ($U_\lambda$) : Un opérateur unitaire qui introduit des corrélations spin-boson au-delà du niveau de champ moyen. Il agit comme un déplacement dépendant du spin (transformation de type polaron) contrôé par un paramètre variationnel $\lambda$. Cette couche transforme le couplage linéaire spin-boson en interactions spin-spin effectives.
3. Secteur de Spin Effectif ($|\phi_s\rangle$) : L'état de spin n'est pas un simple état produit mais l'état fondamental exact d'un Hamiltonien de spin effectif ($H_{\text{eff}}$). Cet Hamiltonien est obtenu en moyennant l'Hamiltonien complet sur l'état gaussien bosonique et la transformation de revêtement.

La méthode emploie une boucle d'optimisation auto-cohérente :

• Pour des paramètres bosoniques fixes ($\Delta_R, \xi, \lambda$), l'Hamiltonien de spin effectif est résolu via la DMRG (Density Matrix Renormalization Group) pour obtenir $|\phi_s\rangle$ et les valeurs d'attente du spin.
• Ces observables de spin sont utilisées pour mettre à jour les paramètres variationnels bosoniques en minimisant la fonctionnelle d'énergie variationnelle totale.
• La boucle itère jusqu'à convergence de l'énergie et des observables sélectionnées.

Cette approche évite la troncature explicite du champ photonique. Au lieu de cela, les degrés de liberté bosoniques sont intégrés variationnellement, projetant le problème sur un modèle de spin auto-cohérent.

Contributions Clés et Résultats
L'article évalue ce cadre NGS-DMRG par rapport aux modèles Dicke et Dicke-Ising, en comparant les résultats avec des simulations DMRG de référence convergées (utilisant des coupures explicites de photons).

1. Précision et Amélioration Variationnelle :

   • L'approche non-gaussienne produit systématiquement des énergies de l'état fondamental plus basses que la limite gaussienne (champ moyen), confirmant l'amélioration variationnelle.
   • Pour le modèle de Dicke, la méthode reproduit avec précision la transition de phase quantique (QPT) superradiante et le scaling de taille finie, correspondant aux résultats de la DMRG de référence avec des ressources de calcul nettement réduites.
   • Dans le modèle Dicke-Ising, le cadre cartographie avec succès un diagramme de phases complexe présentant :
     • Des transitions continues (second ordre) entre les phases normales et superradiantes.
     • Des transitions discontinues (premier ordre) et des points tricritiques dans le régime ferromagnétique.
     • Une phase intermédiaire étroite de coexistence de l'ordre antiferromagnétique et de la superradiance (AFM-SP) dans le régime antiferromagnétique, une caractéristique prédite par la théorie du champ moyen mais difficile à résoudre numériquement.
   • La méthode capture les fluctuations de taille finie et converge vers la limite thermodynamique, où l'intrication spin-photon s'annule, conformément aux limites analytiques établies.

2. Efficacité Computationnelle :

   • Réduction de la Dimension de Liaison : En traitant les bosons de manière variationnelle, le solveur DMRG opère sur un Hamiltonien effectif composé uniquement de spins. Cela résulte en une dimension de liaison MPS ($\chi$) nettement plus faible que celle de la DMRG de référence, qui doit gérer explicitement la chaîne composite spin-boson.
   • Temps d'Exécution : Malgré le surcoût de la boucle auto-cohérente, la réduction de la complexité du réseau de tenseurs conduit à une diminution substantielle du temps d'exécution total. La boucle d'optimisation externe démontre une convergence exponentielle.
   • L'avantage est plus marqué dans les régimes de fortes corrélations et de frontières de phases complexes (par exemple, la phase de coexistence), où la DMRG de référence nécessite de très grandes dimensions de liaison et des coupures de photons élevées.

Signification
L'article affirme que ce cadre hybride offre une voie algorithmique distincte par rapport aux méthodes de réseaux de tenseurs d'optimisation de base locale (LBO). En capturant les modes photoniques via des paramètres variationnels (déplacement, compression, revêtement) plutôt que par une reconfiguration itérative de la base, la méthode permet au solveur DMRG de se concentrer exclusivement sur les corrélations de spin.

La signification réside dans la fourniture d'une méthode numérique robuste et non-perturbative capable d'accéder à la physique spin-boson fortement corrélée sans éliminer ni tronquer arbitrairement le champ bosonique. Les auteurs démontrent que cette approche peut reproduire avec précision des diagrammes de phases complexes et des comportements critiques dans des modèles tels que la chaîne Dicke-Ising, atteignant une haute fidélité avec un coût de calcul réduit. Les auteurs fournissent une implémentation open-source qui supporte les couplages spin-spin de portée arbitraire et le désordre (disponible sur : https://github.com/Jpedromend/NGS), avec des extensions potentielles vers des systèmes multi-modes et des géométries de plus haute dimension via des solveurs alternatifs (ex: PEPS, QMC).