Dynamique quantique non-markovienne

Onset non-markovien a g*tau_c > 0.1

Distance-trace > 1% a g*tau_c > 0.1

Perte de cohérence différentielle par régime

MET: dimension d'embedding = rang de Hankel

Borne tight sur N_G (Edge B)

Noyau de mémoire K(tau) = sum a_k exp(-mu_k tau). Le régime markovien correspond à une limite delta (pointillés). Les modes lents (faible Re(mu)) maintiennent le flux d'information inverse plus longtemps.

Théorème MET (2.1, PROUVÉ) : le système original S + environnement E avec noyau de mémoire K(tau) rationnel est exactement équivalent à une dynamique Lindblad dans l'espace élargi S+anc. de dimension = rang de Hankel du noyau.

Cohérence rho_01(t) pour les trois régimes. En dessous du seuil P_c = 0.098 (g*tau_c < 0.1) : décroissance exponentielle (markovien). Au-dessus : oscillations et flux inverse caractéristiques de la dynamique non-markovienne.

N_G géométrique par régime (barres) et valeur propre de la matrice de Choi à fréquence nulle (bande). Vert : CP-divisible. Bleu foncé : flux inverse fort (mode lent, N_G moyen = 1.67 vs 0.19 pour mode rapide).

Seuil H2 dépassé : NM détectable (SNR > 3)

MODÈLE ILLUSTRATIF - Démonstration déterministe fondée sur les seuils reportés dans la publication. Pas une simulation quantique en direct.

Faiblement NM (0.1 < g·τ_c < 1)

Physique quantique · Systèmes ouverts · Dynamique non-markovienne

Dynamique quantique non-markovienne via Nakajima-Zwanzig

Formalisme exact de projection pour les systèmes quantiques ouverts sans hypothèse de structure Lindblad présupposée

Les équations de Lindblad-GKSL (Gorini-Kossakowski-Sudarshan 1976, Lindblad 1976) constituent la pierre angulaire de la dynamique quantique ouverte sous l'hypothèse de Markov. Elles supposent que le couplage système-bain est si faible et la mémoire du bain si courte que l'évolution du système ne dépend pas de son histoire passée. Cette hypothèse, conceptuellement simple, est structurellement incorrecte pour les environnements à corrélation temporelle non-négligeable.

Ce travail part de l'équation de Nakajima-Zwanzig (NZ) exacte et montre que toute dynamique à noyau de mémoire rationnel admet une représentation markovienne exacte dans un espace élargi. Ce résultat le Théorème d'Embedding Markovien (MET) unifie la non-markovianité avec la théorie des systèmes ouverts sous une seule structure géométrique à quatre cônes.

Par la méthode de projection (opérateurs P et Q = 1 - P), on isole la dynamique du sous-système S de celle du bain E. L'équation exacte pour rho_S(t) est :

K_NZ(t,s) : super-opérateur de mémoire exact

H_LS : Hamiltonien de Lamb shift

Le noyau K_NZ(t,s) encode toute l'histoire des corrélations système-bain. Sous la forme de Prony (noyau rationnel) : K(tau) = sum_k a_k exp(-mu_k tau), avec Re(mu_k) > 0, ce noyau admet une représentation finie exacte. La limite de Markov correspond à K(tau) to mathcalL_GKSL * delta(tau) quand tau_c to 0 (Théorème 2.2).

Si K(tau) est un noyau rationnel de rang de Hankel r,

il existe un système élargi (S + r ancilles) tel que :

Dimension minimale de l'espace élargi = rang de Hankel du noyau

Le rang de Hankel du noyau détermine la dimension minimale de l'espace élargi nécessaire. Pour un noyau à un seul mode (K = a exp(-mu tau)), une seule ancille suffit. Pour r modes, r ancilles sont nécessaires et suffisantes. Ce résultat généralise le théorème de dilatation de Stinespring à la dynamique non-markovienne.

Corollaire 2.2 (vérifié numériquement) : l'algorithme d'extraction du rang de Hankel identifie la dimension minimale d'embedding avec 396/400 essais concordants dans la validation numérique (taux de succès 99%).

Condition d'admissibilité et taux markovien asymptotique

L'opérateur K(0) = integral_0^inf K(tau) dtau est le taux markovien asymptotique. La condition d'admissibilité exige que K(0) soit un générateur valide de semigroupe GKSL. Lorsque cette condition est satisfaite, la dynamique NZ converge bien vers la limite de Markov quand tau_c to 0.

Définition 3.2 : mesure géométrique N_G

La mesure N_G quantifie le flux d'information inverse (backflow) depuis le bain vers le système. Formellement : N_G = sup_S_0 (1)/(V(0)) integral [dV/dt]_+ dt, où V(t) est le volume de trace d'un ensemble d'états initiaux S_0 et [x]_+ = max(0,x) désigne la partie positive.

N_G corrige la mesure de Lorenzo-Plastina-Paternostro (2013) qui confondait P-divisibilité et volume-monotonie. La structure à quatre cônes (Edge A du traité, prouvée dans Théorème 5.1) établit : semigroupe contenu dans CP-div contenu dans P-div contenu dans volume-monotone contenu dans toutes les dynamiques CPTP. N_G est fidèle à la CP-divisibilité, pas à la P-divisibilité.

Théorème 3.3 : CP-divisibilité via la matrice de Choi

La dynamique est CP-divisible si et seulement si :

C_K(tau) : matrice de Choi instantanée du noyau K(tau)

Condition nécessaire et suffisante (pas d'approximation)

La borne tight sur N_G (Edge B, 400 essais de validation) : N_G le C * sum_k | a_k| / Re(mu_k), avec C ~ 1.55 (borne empirique la plus serrée). Les modes lents (Re(mu_k) petit) génèrent un backflow plus important : N_G moyen = 1.67 pour les modes lents vs 0.19 pour les modes rapides. N_G décroît monotonement avec min_k Re(mu_k).

G_int = integral | rho_exact(t) - rho_Markov(t)|_1 dt quantifie l'erreur cumulée de l'approximation de Markov sur toute la trajectoire. G_int est petit quand gtau_c < 0.1 (régime markovien) et devient significatif au-dessus du seuil P_c = 0.098. Pour le dephasing pur (cas exact de Théorème 3.5), gamma(t) < 0 implique exactement non-markovien.

L'article formule 10 hypothèses explicitement falsifiables avec seuils numériques, testables sur les qubits supraconducteurs (IBM, Google, Rigetti). Chaque hypothèse spécifie la plateforme ciblée, le type d'expérience et le seuil de rejet. Les six hypothèses principales sont résumées ci-dessous ; H7 à H10 concernent les ensembles quantiques, la discrimination de modèles et la scalabilité.

H1 : onset non-markovien à gtau_c > 0.1 (P_c = 0.098 validé numériquement), SNR > 3, p < 0.01

H2 : distance-trace > 1% quand gtau_c > 0.1 (seuil observable expérimentalement)

H3 : N_G distingue CP-div de P-div avec DeltaBIC > 10 (confirmé DeltaBIC = 42)

H5 (MET) : dimension d'embedding = rang de Hankel ; Corollaire 2.2 : 396/400 essais concordants

H6 (Edge B) : borne tight N_G le C * sum | a_k| / Re(mu_k), C ~ 1.55, 0 violation / 400 essais

Théorème 5.1 : unification à quatre faces (PROUVÉ/VÉRIFIÉ-NUM)

Le Théorème 5.1 établit la structure à quatre cônes de la divisibilité quantique. Il est partiellement prouvé analytiquement et partiellement vérifié numériquement sur 597 essais aléatoires.

H1 confirmé : 98% des essais avec SNR = 8.3 montrent l'onset à P_c = 0.098

H2 confirmé : P_c = 0.098, distance-trace > 1% confirmée au-delà du seuil

H3 confirmé : DeltaBIC = 42 (largement au-dessus du seuil de 10)

H4 confirmé : DS = 0.15 bits (> 0.1 bits requis)

H5 partiel (85%) : R = 1.45 ; les cas non-concordants correspondent à des noyaux à proximité de dégénérescence de rang

H6 confirmé : d_L = 0.23 (distance de violation nulle), C ~ 1.55 sur 400 essais

Corollaire 2.2 : 396/400 essais concordants pour l'extraction du rang de Hankel

Théorème 5.1 face (iv) : 0 violation sur 597 essais aléatoires

Le code de validation numérique est publié sur GitHub : https://github.com/hopenmind/metcore. Il comprend les implémentations Python des algorithmes MET, de l'extraction du rang de Hankel, du calcul de N_G, de la matrice de Choi et des 10 protocoles d'hypothèses. La suite de tests reproduit tous les résultats de validation reportés.

Cet article est publié en accès ouvert sur Zenodo sous licence CC BY-NC-SA 4.0. Le code de validation est disponible sous licence MIT sur GitHub.

Acces libre · CC BY-NC-SA 4.0 · HOPE 'N MIND SASU · 2026