La mise en œuvre pratique d'algorithmes numériques pour intégrer des équations maîtresses pour matrices de densité, incluant choix de discrétisation, intégrateurs, représentations creuses et en tenseurs, exponentiation par sous‑espace de Krylov, et démêlages stochastiques, avec attention à la préservation numérique de la trace et de la positivité et au coût de calcul.

Principe

La simulation numérique doit équilibrer précision, stabilité, efficacité et physicalité : choisir des pas temporels adaptés aux taux les plus rapides pertinents, exploiter la parcimonie et les quantités conservées, utiliser des intégrateurs respectant la structure du générateur ou appliquer un post‑traitement pour restaurer la positivité, et estimer erreurs et mise à l'échelle des ressources.

Démonstration

Simuler un oscillateur harmonique amorti et drivé dans une base de Fock tronquée peut utiliser une représentation creuse des superopérateurs de Lindblad et des méthodes de sous‑espace de Krylov pour l'action de exp(tL), ou employer Monte‑Carlo des sauts quantiques pour échantillonner des trajectoires lorsque la dimension de l'espace de Hilbert rend la propagation complète de la matrice de densité infaisable.

Mauvaise application

Employer des intégrateurs explicites avec des pas temporels trop grands sur des générateurs dissipatifs raides conduit à des populations négatives ou à la perte de trace ; une troncature négligente de l'espace de Hilbert sans vérifications de convergence introduit des canaux de relaxation artificiels ; ne pas moyenner suffisamment de trajectoires donne d'importantes erreurs statistiques.

Conséquence

Des simulations numériques bien conçues produisent des séries temporelles et des états stationnaires précis dans des tolérances quantifiées et rendent des problèmes à grande échelle tractables ; de mauvaises pratiques numériques conduisent à des résultats non physiques, au gaspillage de ressources et à des conclusions trompeuses sur la dynamique et les propriétés stationnaires.

Inversion

Une alternative à l'intégration numérique est la solution analytique lorsque disponible ou la diagonalisation directe du Liouvillien complet pour des systèmes très petits ; elles évitent les erreurs de pas temporel mais ne sont pas extensibles aux grands espaces de Hilbert.

Limite

Cible des espaces de Hilbert finis ou tronqués et des équations maîtresses locales en temps ou convenablement converties depuis des formes non locales ; exclut les traitements de bain continu sans discrétisation et les méthodes négligeant la stabilité et les contraintes de conservation intrinsèques au problème quantique.

Tension sémantique

Tension entre l'utilisation de démêlages stochastiques (efficaces en mémoire, erreurs statistiques) et de solveurs déterministes de la matrice de densité (sans bruit d'échantillonnage, mémoire plus élevée) ; autre tension entre imposer strictement la positivité complète à chaque étape et minimiser la charge de calcul.

Synthèse

La simulation numérique d'équations maîtresses est la sélection et l'application disciplinées de discrétisations, algorithmes de résolution, troncations de base et contrôles d'erreur qui réalisent une évolution temporelle et un calcul d'état stationnaire précis, stable et physiquement cohérent pour systèmes quantiques ouverts dans les ressources informatiques disponibles.