Programme de Reconstruction Modulaire

I. Le problème

Conjecture. Soit $|\psi\rangle$ un état quantique dans un espace de Hilbert $\mathcal{H}$, et soit $\{A\}$ le treillis des sous-algèbres de von Neumann d'observables localisées. À chaque $A$ est associé un opérateur modulaire $\Delta_A$ de Tomita-Takesaki relatif à $|\psi\rangle$. Alors la collection $\{\Delta_A\}_{A}$ détermine uniquement, à isométrie près, la géométrie lorentzienne (si elle émerge) associée à $(\mathcal{H}, |\psi\rangle, \{A\})$.

Cette conjecture est le pari central du cadre Bath-TT relationnel, posé dans l'entrée 042, §V. Elle n'est pas établie. Elle n'a pas même, à ma connaissance, été posée dans ces termes-là par la littérature. Des fragments existent, aucun n'est global.

Ce document inaugure le programme de recherche correspondant. Il se met à jour à mesure que les jalons sont atteints.

II. Pourquoi l'averaging ne suffit pas

L'approche thermodynamique — Jacobson 1995, Clausius local, équation d'état — travaille exclusivement au niveau du scalaire d'entropie $S(A)$. Elle dérive les équations d'Einstein sur un espace-temps déjà donné. Elle ne construit pas l'espace-temps ; elle lit ses lois effectives.

Ryu-Takayanagi 2006 dépasse ce cadre dans AdS/CFT en donnant un algorithme constructif (aire minimale) pour relier $S(A)$ à une métrique holographique. Mais cet algorithme dépend crucièrement de la structure asymptotique anti-de-Sitter. Il ne se généralise pas.

Van Raamsdonk 2010 est qualitatif : couper l'intrication coupe l'espace. Slogan éclairant, pas algorithme.

Le problème commun : $S(A)$ est un scalaire par sous-algèbre. Une géométrie lorentzienne requiert infiniment plus d'information. L'écriture

cache derrière elle un opérateur $\rho_A$ (infiniment-dimensionnel) et son logarithme $K_A = -\log \rho_A$ (le hamiltonien modulaire). Le spectre de $\rho_A$ est une fonction réelle ; $S(A)$ en est un moment intégré. Reconstruire la géométrie demande le spectre entier, pas son moment.

III. Les phases

IV. Premier jalon : CFT libre 1+1D

Le cas minimal où la conjecture peut être testée concrètement. Ce qui suit est une reformulation rigoureuse de résultats existants, rassemblés en un seul énoncé.

Setup

Soit $\phi$ un champ scalaire réel sans masse en dimension 1+1 sur l'espace de Minkowski. Son vide $|0\rangle$ est l'état fondamental conforme. À chaque intervalle ouvert $I = (a, b) \subset \mathbb{R}$ sur la tranche $t = 0$, on associe l'algèbre de von Neumann $\mathcal{A}(I)$ engendrée par les opérateurs localisés dans le domaine de dépendance causale $\mathcal{D}(I) = \{(t, x) : |t - (a+b)/2| < (b-a)/2 - |x - (a+b)/2|\}$.

Le modulaire explicite

Théorème (Hislop-Longo 1982, spécialisé en 1+1D). Le hamiltonien modulaire de $\mathcal{A}(I)$ relatif au vide est :

Le flot modulaire $\sigma_t^I = e^{iK_I t}$ agit géométriquement comme la transformation de Möbius qui préserve les points $a$ et $b$ et fixe la règle de paramétrisation $e^{-2\pi t}$.

Extraction de l'information géométrique

(i) Longueur. L'entropie d'intrication est bien connue (Calabrese-Cardy 2004) :

La longueur $|I|$ est extractible du spectre de $\Delta_I$ via cette entropie.

(ii) Position relative. Pour deux intervalles disjoints $I_1 = (a_1, b_1)$ et $I_2 = (a_2, b_2)$, la structure modulaire conjointe de $\mathcal{A}(I_1) \vee \mathcal{A}(I_2)$ dépend du cross-ratio

La collection $\{\eta(I_i, I_j)\}$ sur toutes les paires détermine les positions relatives jusqu'à translation globale et réflexion.

(iii) Orientation temporelle. Le flot $\sigma_t^I$ a une direction privilégiée : pour $t > 0$, il contracte vers $b$ (avec $b > a$). Ceci fixe l'orientation future du cône de lumière via Bisognano-Wichmann inversé : le flot ressemblant à un boost vers la droite identifie la direction spatiale positive.

(iv) Cutoff conforme. La composante $c/3$ du coefficient de $\log |I|$ dans $S(I)$ identifie la théorie comme un scalaire libre ($c = 1$) parmi les CFT 1+1D.

Théorème de reconstruction (Phase 1)

Ce résultat n'est pas nouveau dans ses pièces (Bisognano-Wichmann, Hislop-Longo, Calabrese-Cardy, Casini-Huerta sont tous antérieurs). Ce qui est nouveau, c'est l'énoncé unifié comme théorème de reconstruction : la donnée modulaire $\{\Delta_I\}$ seule suffit. Pas d'averaging. Pas de RT-type recipe holographique. Juste le flot modulaire sur chaque sous-algèbre.

C'est la base. La phase 1 est accomplie dans son cas minimal.

V. Deuxième jalon : vérification numérique sur chaîne harmonique

Le Théorème 1.1 est un énoncé analytique. La question pratique : sur un modèle lattice concret, les ingrédients géométriques sont-ils numériquement extractibles de la donnée modulaire ? Ce test est fait.

Méthode

Chaîne harmonique discrète à $N = 400$ sites, conditions aux bords ouvertes, masse infrarouge $m = 0.001$ pour régulariser le zéro-mode. Méthode de Peschel (J. Phys. A 36, L205, 2003) : le spectre modulaire d'un état gaussien se calcule depuis les matrices de corrélation $\langle q_i q_j\rangle$ et $\langle p_i p_j\rangle$ restreintes à une sous-région. Script modular-chain.sage, exécuté avec SageMath 10.8.

Résultats qualitatifs

Test 1 — Longueur depuis l'entropie. $S(|I|)$ suit un scaling logarithmique avec $\log(\text{chord})$, confirmant Calabrese-Cardy. Coefficient numérique ajusté : 0.206. Prédiction CFT naïve $c/6 = 0.167$ pour $c = 1$. Écart 24%.

Test 2 — Longueur depuis le spectre modulaire. Le gap entre les plus basses valeurs propres $\epsilon_k$ du hamiltonien modulaire varie inversement avec $\log|I|$ :

Accord qualitatif fort avec $\epsilon_k \sim 2\pi k/\log|I|$, quantitatif à facteur $\sim 2$.

Test 3 — Position relative depuis cross-ratios. L'information mutuelle $I(I_1:I_2) = S(I_1) + S(I_2) - S(I_1 \cup I_2)$ dépend manifestement du cross-ratio $\eta$ : 0.34 à $\eta = 0.69$, 0.06 à $\eta = 0.15$. La position relative des deux intervalles est bien encodée dans la donnée modulaire.

Test 4 — Structure universelle. Le mode fondamental $\epsilon_1 \cdot \log|I|/(2\pi) \approx 0.6$, compatible avec la valeur universelle O(1) prédite pour les CFT 1+1D libres.

Désagréments quantitatifs (honnêtes)

Les écarts de 20-100% avec les formules CFT simplifiées ont trois causes identifiées :

(a) Les formules analytiques utilisées (ex. $I_\text{mut} = -(c/6)\log(1-\eta)$) sont des approximations dans des régimes restreints. La vraie formule CFT libre pour l'information mutuelle (Calabrese-Cardy-Tonni 2011) implique plusieurs cross-ratios et des thêta-fonctions de Jacobi. Elle n'est pas implémentée dans le script actuel.

(b) Les corrections finies-$N$ et les contributions des modes de bord (OBC avec intervalle central) ne sont pas modélisées.

(c) La masse régulatrice $m = 0.001$ introduit une légère déviation du scaling conforme exact.

Ce que ce jalon établit

Acquis. Les ingrédients géométriques (longueur, position relative, structure universelle) sont présents dans la donnée modulaire $\{\Delta_I\}$ numérique. La conjecture de reconstruction n'est pas réfutée par ces tests lattice ; elle est numériquement plausible au niveau qualitatif.

Reste à faire — Jalon 2.1-bis. L'extraction quantitative exacte du continuum limit. Cela demande : (i) implémenter les formules CFT complètes (Calabrese-Cardy-Tonni) ; (ii) faire le scaling propre en $N \to \infty$, $m \to 0$ avec rapport $m\xi$ fixé ; (iii) isoler les contributions UV vs IR dans l'entropie. Modérément difficile : travail de codage et soin numérique, sans obstacle conceptuel.

VI. Troisième jalon : cas massif en 1+1D

Le Théorème 1.1 est établi pour le cas conforme ($m = 0$). La question suivante : que se passe-t-il lorsque l'invariance conforme est brisée par une masse non-nulle ? Si la conjecture de reconstruction modulaire est vraie génériquement, elle doit tenir au-delà des CFT. Test fait sur chaîne harmonique massive, $N = 400$ sites, scan sur six ordres de grandeur : $m \in \{10^{-3}, 10^{-2}, 0.05, 0.1, 0.3, 1.0\}$. Script modular-mass-scan.sage.

Scan 1 — Transition conforme → massif dans l'entropie

L'entropie $S(|I|)$ mesurée sur intervalles centrés de tailles $|I| \in \{8, 16, 32, 64, 128\}$, et ajustement logarithmique $S = c_\text{log} \log|I| + \text{const}$ :

Le coefficient logarithmique décroît continuement de 0.206 (quasi-CFT) vers 0.000 (loi d'aire stricte). La transition s'effectue autour de $|I| \approx \xi$. L'échelle de saturation trahit la masse : elle est extractible de la seule donnée $S(|I|)$.

Scan 2 — Fermeture du gap modulaire

Le gap $\varepsilon_1$ entre les deux plus basses valeurs propres du hamiltonien modulaire :

Deux comportements distincts : en régime conforme, $\varepsilon_1 \sim 1/\log|I|$ (décroissance lente attendue) ; en régime massif, $\varepsilon_1$ s'effondre dès que $|I|$ dépasse $\xi$. Le spectre modulaire lui-même encode la structure massive.

Scan 3 — La signature la plus remarquable

Entropie d'un intervalle de taille fixe $|I| = 32$, déplacé à cinq positions le long d'une chaîne de 400 sites :

À $m = 0.001$ (quasi-CFT), l'entropie dépend de la position — la proximité des bords ouverts brise la translation-invariance, exactement comme en CFT avec bords. À $m \geq 0.05$, l'entropie devient strictement translation-invariante, à quatre décimales près. Les corrélations s'épuisent dans une enveloppe de taille $\xi$ trop petite pour voir les bords, et l'observable devient bulk-only.

Ce que ce jalon établit

Acquis. La reconstruction modulaire n'est pas un artefact de l'invariance conforme. En présence de masse, la donnée $\{\Delta_A\}$ :

(i) continue d'encoder la longueur des intervalles via le spectre modulaire ;
(ii) encode la masse via l'échelle de saturation de l'entropie ;
(iii) distingue les régimes (local/bulk vs boundary-sensitif) selon $|I|/\xi$.

Conjecture renforcée. En théorie massive 1+1D, $\{\Delta_A\}$ détermine conjointement (i) la géométrie minkowskienne émergente, (ii) la masse du champ. La signature modulaire de $m$ est l'échelle de bulk-ification de l'entropie.

Reste à faire — Jalon 1.1-bis. Formaliser analytiquement la conjecture ci-dessus. Ingrédients disponibles : Borchers (1992) pour le flot modulaire non-géométrique, Longo-Xu pour les théories avec spectre massif, Hollands-Sanders pour la propriété d'Haag en contexte massif. Le saut analytique : prouver que l'échelle lue dans le scan $|I| \to \infty$ coïncide avec la longueur de Compton. Non-trivial ; ni bloqué.

VII. Quatrième jalon : passage en 3+1D

Le saut qui compte. Tous les résultats précédents vivent en dimension 1+1 — dimension suffisamment basse pour qu'on puisse raisonnablement soupçonner des artefacts de pauvreté géométrique. En $d = 1$ spatial, il n'y a qu'une direction ; intervalles et cross-ratios suffisent. En $d = 3$, la géométrie devient strictement plus riche : orientation, volume vs surface, fonctions à distance en $1/r$ plutôt que logarithmiques.

Test fait : réseau cubique $N = 20$ soit $8000$ sites, conditions aux bords périodiques, champ scalaire libre régularisé IR par $m = 0.01$. Calcul FFT des corrélations du vide, boules centrées de rayons $R \in \{2, 3, 4, 5, 6\}$, puis méthode de Peschel sur les sous-matrices. Script modular-3d-scan.sage.

Scan 1 — Loi d'aire, signature de la dimension

Ajustement linéaire $S = \alpha R^2 + \beta$ :

La qualité du fit est remarquable : $S/\text{aire}$ reste entre 0.031 et 0.034 sur toute la plage, alors que $S/|B|$ décroît de 0.051 à 0.015. La dimension spatiale est lisible dans la donnée modulaire : scaling en $R^2$ contre $\log R$ en $d = 1$ ou $R^1$ en $d = 2$.

Scan 2 — Invariance par translation à la précision machine

Boule de rayon $R = 4$ placée à quatre centres différents du réseau PBC, dont le coin $(0,0,0)$ :

Scan 3 — Information mutuelle et structure globale

Deux boules de rayon $r = 2$, distance $d$ variable. L'information mutuelle $I(B_1 : B_2) = S(B_1) + S(B_2) - S(B_1 \cup B_2)$ encode la corrélation entre les régions :

Décroissance monotone pour $d \leq 7$, puis un uptick à $d = 10$. Cet uptick n'est pas du bruit : c'est la signature topologique du tore — pour $N = 20$, deux points à distance nominale 10 sont également à distance 10 par le chemin inverse (périodicité). La géométrie du tore discret est lisible dans la décroissance de l'information mutuelle.

Ce que ce jalon établit

Acquis. En 3+1D avec tranche spatiale $\mathbb{T}^3$ (tore PBC), la donnée modulaire $\{\Delta_B\}$ encode :

(i) la dimension spatiale $d = 3$ via le scaling en $R^2$ de l'entropie ;
(ii) l'homogénéité via l'invariance du spectre sous translation (précision machine) ;
(iii) la distance via la décroissance de l'information mutuelle ;
(iv) la topologie ($\mathbb{T}^3$ vs $\mathbb{R}^3$) via les corrections périodiques visibles dans $I(d)$.

Conjecture confirmée génériquement. Le résultat 1+1D n'était pas un artefact dimensionnel. La reconstruction modulaire fonctionne en $d = 3$ avec une précision remarquable, sans aucun ajustement de paramètres ni hypothèse supplémentaire.

Reste à faire — Jalons 1.2-bis et 1.2-ter.

1.2-bis : reconstruction de l'orientation (groupe $SO(3)$) via des régions non sphériques (ellipsoïdes). Si la conjecture tient pleinement, le spectre modulaire d'un ellipsoïde doit distinguer ses axes principaux. Test calculable, en cours de conception.

1.2-ter : structure causale lorentzienne via diamants 3+1D (pas juste la tranche spatiale). Demande de sortir du vide-au-temps-fixé pour échantillonner des régions qui s'étendent dans le temps. Plus délicat numériquement.

VIII. Cinquième jalon : orientation $SO(3)$

La dimension et l'homogénéité sont reconstruites (§VII). Il reste la troisième structure de $\mathbb{R}^3$ — les rotations. Si la donnée modulaire encode la géométrie euclidienne complète, elle doit distinguer les formes (sphère vs ellipsoïde) et refléter correctement les symétries rotationnelles.

Test : ellipsoïdes de semi-axes $(a, b, c)$ à divers orientations, même réseau que §VII ($N=20$, PBC, $m=0.01$). Script modular-ellipsoid-scan.sage.

Scan 1 — La forme est lisible

L'entropie varie de 4.3 à 7.6 selon la forme — un facteur presque 2 à volume quasi-constant. La forme est pleinement lisible.

Scan 2 — Symétrie cubique exacte

Ellipsoïde prolate de semi-axes $(2, 2, 5)$ aligné le long des trois axes du réseau. Distance maximale entre les spectres (triés) :

Test complémentaire : ellipsoïde tri-axial $(3, 4, 5)$ avec les 4 permutations testables des axes. Toutes les permutations donnent des spectres identiques à $10^{-15}$ près. La symétrie cubique du réseau (sous-groupe de $SO(3)$) est préservée de manière absolument exacte.

Scan 3 — Discrétisation des rotations hors-réseau

Même ellipsoïde $(2, 2, 5)$ roté par des angles non-alignés sur le réseau :

Les rotations non-cristallographiques donnent 0.5–6% d'écart, avec $|B|$ qui varie parce que la boule digitalisée change d'un nombre entier de sites. Ce sont des artefacts de discrétisation purs : en continuum limit ($N \to \infty$ avec $R \cdot a \to $ const), ils doivent tendre vers zéro selon $O(1/R)$ (boundary term). L'erreur grandit quand la rotation s'éloigne des symétries cubiques ; c'est exactement la signature attendue.

Ce que ce jalon établit (bilan complet $\mathbb{R}^3$)

Ensemble avec §VII, le bilan pour la tranche spatiale est remarquable :

(i) dimension $d = 3$ via loi d'aire $S \propto R^2$ ;
(ii) homogénéité via invariance par translation à $10^{-15}$ ;
(iii) topologie ($\mathbb{T}^3$ vs $\mathbb{R}^3$) via structure périodique de $I(d)$ ;
(iv) isotropie $SO(3)$ via variation contrôlée du spectre sous forme/orientation.

Le groupe euclidien complet $E(3) = \mathbb{R}^3 \rtimes O(3)$ est donc lu dans la donnée $\{\Delta_B\}$ des boules et ellipsoïdes sur le réseau. La tranche spatiale de la conjecture est pleinement reconstructible.

Ce qui reste. Le temps et le caractère lorentzien. La tranche $t = 0$ n'est pas le plein espace-temps : l'information temporelle vit dans les diamants causaux qui s'étendent dans le temps. C'est l'objet du Jalon 1.2-ter.

IX. Sixième jalon : structure temporelle via profil de Casini

Jusqu'ici, toute la structure extraite vit sur la tranche spatiale $t = 0$. Le temps n'apparaît pas explicitement. Mais la géométrie de Minkowski $\mathbb{R}^{3,1}$ n'est pas qu'un produit $\mathbb{R} \times \mathbb{R}^3$ ; elle est l'espace-temps tout entier avec sa structure causale. Comment extraire la direction temporelle de données collectées à un seul instant ?

Réponse : par le théorème de Casini-Hislop-Longo. Pour une boule $B$ de rayon $R$ centrée à l'origine au temps $t=0$ dans une CFT, le hamiltonien modulaire est explicitement

Le poids $f(x) = (R^2 - |x|^2)/(2R)$ n'est pas arbitraire : c'est la composante $\xi^0$ du vecteur de Killing conforme qui engendre le flot modulaire de la diamant causal $\mathcal{D}(B) = \{(t,x) : |t| + |x| < R\}$. Ce vecteur définit la direction « temporelle » propre à la diamant, et le flot modulaire $\sigma_t = e^{iK_B t}$ l'incarne comme un boost conforme. Extraire $f(x)$ des données modulaires revient à lire la direction temporelle lorentzienne.

Méthode : retrait de site et première loi d'intrication

Bhattacharya-Nozaki-Takayanagi (2013) : à l'ordre dominant, $\delta S = 2\pi \langle \delta K_B \rangle$. Pour le vide avec $\langle T_{00} \rangle = \varepsilon_\text{vac}$ uniforme, une perturbation qui retire un site $x$ de $B$ donne

Le signe : retirer un site intérieur (haut $f$) crée une cavité dans $B$, augmente l'aire totale, donc augmente $S$. $\delta S < 0$. Retirer un site au bord (bas $f$) rétrécit l'aire extérieure. $\delta S > 0$. La corrélation prédite est négative.

Résultat : le profil conforme est extrait

Script modular-time-scan.sage. Réseau 20³ PBC, même setup que §VII-VIII.

La corrélation entre $\delta S(x)$ et le profil de Casini est de $0.90$ à $0.99$ à travers quatre rayons. Le signe est exactement négatif comme prédit. Le coefficient universel $|B| \cdot R \approx 0.27$ est quasi-stable sur la plage testée, reflet de la nature UV (indépendante de $R$) de la densité de vide $\varepsilon_\text{vac}$.

Isotropie du profil

Vérification : pour $R = 4$, $\delta S(x)$ doit dépendre uniquement de $|x|$, pas de la direction (le vecteur de Killing est radial).

Les trois axes donnent $-0.06978$ à la précision d'affichage ; les trois diagonales donnent $-0.07753$. L'isotropie résiduelle est exacte sur les symétries cubiques du réseau. Le profil est radial.

Bilan dimensionnel complet

Avec ce résultat, le programme clot sa séquence dimensionnelle :

(i) dimension spatiale $d = 3$ — loi d'aire $S \propto R^2$ (§VII)
(ii) homogénéité $\mathbb{R}^3$ — translation à $10^{-15}$ (§VII)
(iii) topologie $\mathbb{T}^3$ vs $\mathbb{R}^3$ — $I(d)$ avec wrap-around (§VII)
(iv) isotropie $SO(3)$ — formes et rotations (§VIII)
(v) direction temporelle — profil de Casini / CKV (§IX, ici)

Le groupe de Poincaré complet $\mathbb{R}^{3,1} \rtimes O(3,1)$ — ou plutôt sa structure isométrique locale — est extractible de $\{\Delta_B\}$ sur les seules régions convexes d'une tranche spatiale. La conjecture tient en dimension physique.

X. Septième jalon : cas massif en 3+1D

Le carré logique compte quatre coins : massless/massif croise 1+1D/3+1D. Trois sont testés. Reste le dernier : que devient le profil de Casini quand la masse brise la conformité en dimension physique ?

Hypothèse naturelle : le profil doit se dégrader à l'échelle $\xi = 1/m$, comme en 1+1D (Jalon 1.1). Pour $R = 4$, la transition devrait être visible vers $m \sim 1/R = 0.25$. Test : scan en masse du test précédent. Script modular-3d-mass-scan.sage.

Résultat : la corrélation est remarquablement stable

La corrélation reste entre $-0.92$ et $-0.94$ alors que $m$ varie sur trois ordres de grandeur et que $\xi = 1/m$ passe de $1000$ à $1$ (traverse $R = 4$ sans transition visible dans la corrélation). Cette stabilité est plus forte que ce qu'on attendait du théorème de Casini-Hislop-Longo, qui suppose la conformité.

Interprétation : le profil de Casini est une signature géométrique pure

La lecture qui émerge : le profil $f(x) = (R^2 - |x|^2)/(2R)$ encode simplement la distance de $x$ au bord de $B$. C'est une propriété de la région elle-même, pas de la théorie qu'elle accueille. Dans l'état fondamental de n'importe quelle théorie raisonnable, retirer un site affecte l'intrication proportionnellement à sa centralité dans la région. Ceci vaut que la théorie soit conforme ou non.

Ce qui change avec la masse : le coefficient global $B$. Il décroît de $-0.071$ à bas $m$ vers $-0.047$ à $m = 1$ (factor 1.5). La masse est donc lue non pas dans la forme du profil, mais dans son amplitude.

Vérification indirecte : transition vers loi d'aire pure

Mesure complémentaire : $S(B)$ pour $R \in \{2, 3, 4, 5\}$ aux différentes masses. Le ratio $S(R=4)/S(R=2)$ doit tendre vers $(4/2)^2 = 4$ en régime loi-d'aire pure, et être inférieur en régime conforme (corrections log) :

Le ratio passe de $3.02$ (corrections log CFT) à $3.9$ (loi d'aire pure, ratio idéal $4$). La transition vers le régime massif est réelle dans le global, même si le profil local reste robuste. Cohérent.

Bilan du carré logique

Les quatre coins :

La conjecture de reconstruction modulaire tient dans les quatre coins accessibles numériquement. Le pattern est clair : {$\Delta_A$} encode conjointement géométrie et masse, avec la masse signalée différemment en 1+1D (dimension où la forme s'aplatit) et 3+1D (dimension où le coefficient change).

XI. Le pont vers la dynamique : première loi d'intrication

Sept jalons accomplis — la partie numérique-exploratoire ferme le carré logique (massless/massif × 1+1D/3+1D) au niveau kinematique. Pour faire émerger la dynamique lorentzienne (les équations d'Einstein via Jacobson 1995), il faut la première loi d'intrication,

Avec $K_B = -\log \rho_B$ le hamiltonien modulaire du vide réduit à $B$. Injectée dans Clausius local le long des horizons de Rindler, cette identité donne les équations d'Einstein.

Deux tentatives incorrectes, une réussie

Tentative 1 (inadaptée) : perturbation de la masse $m \to m + \delta m$ (modular-first-law.sage). Résultat : ratio stable $\approx -0.07$. Diagnostic : changer la masse n'est pas une perturbation d'état au sens du théorème — c'est un changement d'hamiltonien.

Tentative 2 (approximation Casini) : perturbation thermique $\rho_T = e^{-\beta H}/Z$, $K_B$ évalué via la formule de Casini $K_B^\text{HL} = 2\pi \int_B f(x) T_{00}(x) d^3x$ (modular-first-law-thermal.sage). Résultat : ratios divergent. Diagnostic : Casini est valable en CFT exacte ; sur réseau à $m = 0.01 \neq 0$, $K_B$ contient des corrections que Casini ne capture pas.

Tentative 3 (correcte) : perturbation thermique + hamiltonien modulaire $h_0$ extrait numériquement du vide du réseau via la méthode de Peschel généralisée (modular-first-law-exact.sage). Pour un état gaussien bosonique $\rho$ de covariance $\Gamma = \text{diag}(X, P)$ avec $\langle qp\rangle = 0$, on peut expliciter

puis $\langle K_0 \rangle_\rho = \tfrac{1}{2}[\text{Tr}(h_{qq} X_\rho) + \text{Tr}(h_{pp} P_\rho)]$ (à une constante près qui se cancelle dans les différences).

Résultat : la première loi est validée

À $T = 0.001$ (soit $\beta\omega_\text{min} \sim 10$, régime genuinement basse-température), ratio $= 1.0015$, à 0.15 % près de la valeur théorique. La convergence vers 1 est monotone depuis les températures plus hautes. L'entropie relative $S_\text{rel} = \delta\langle K_0\rangle - \delta S$ est partout positive (Klein's inequality) et décroît vers zéro plus rapidement que tout polynôme en $T$ dans le régime propre.

Extension au régime massif

Test supplémentaire à $m = 0.3$, soit $\xi = 3.3 < R = 4$ — régime où la symmétrie conforme est cassée et où la formule de Casini ne devrait pas s'appliquer. Résultat : le ratio converge encore plus vite vers 1 (ratio = 1.008 à $T = 0.01$, intercept du fit linéaire = 0.938). L'entropie relative $S_\text{rel}$ décroît plus rapidement encore (scaling approximatif $\alpha \approx 9$ en log-log) — cohérent avec la décroissance exponentielle $e^{-\beta m}$ dans le régime de gap.

La première loi tient donc indépendamment du régime conforme, dès lors que $h_0$ est extrait numériquement du vide correspondant. C'est la version forte attendue : le théorème n'est pas conditionné à la classe de théories.

Conséquence pour le programme

Avec la première loi validée, la chaîne kinematique $\to$ dynamique est complète :

(i) donnée modulaire $\{\Delta_B\}$ encode la géométrie (Jalons 1.2–1.2-quater)
(ii) première loi relie variations de $S_B$ à variations de $\langle T_{00} \rangle$ (§XI)
(iii) Clausius local injecte (ii) dans la thermodynamique des horizons de Rindler
(iv) identité de Bianchi + conservation donne $R_{\mu\nu} - \tfrac{1}{2}g_{\mu\nu}R = 8\pi G T_{\mu\nu}$ (Jacobson 1995)

La reconstruction de l'espace-temps lorentzien avec sa dynamique est ainsi, sur le cadre numérique teste, conséquence directe de la structure modulaire du vide. Ce qui reste : extension analytique rigoureuse (Jalons 1.1-bis, 3.1) et identification précise de $G$ par la structure microscopique du bain thermique (pont avec Bath-TT).

XII. Huitième jalon : premier test interagissant

Toute la structure établie ci-dessus repose sur des états gaussiens — vide du champ scalaire libre, définis par leur covariance. La question : la reconstruction modulaire est-elle un artefact de la gaussianité ? Si oui, le programme s'arrête aux théories libres. Sinon, il s'étend aux théories interagissantes — le vrai domaine physique.

Test : chaîne d'Ising 1+1D avec champs transverse et longitudinal,

Pour $h_z = 0$, le modèle se réduit à des fermions libres via Jordan-Wigner (intégrable). Pour $h_z \neq 0$, il devient non-intégrable (Kim-Huse, Banuls et al.) et sa fonction d'onde fondamentale n'est pas gaussienne. Diagonalisation exacte sur $L = 12$ sites via Lanczos, traces partielles directes sans hypothèse gaussienne. Script modular-interacting.py.

Profil de Casini 1+1D en régime non-intégrable

Pour un intervalle $I = [a, b)$ en 1+1D, le profil de Casini attendu est $f(x) = (x-a)(b-x)/(b-a)$ (l'analogue 1D de $(R^2 - |x|^2)/(2R)$). Régime critique $h_x = 1$ avec champ longitudinal $h_z = 0.1$, intervalle $|I| = 8$ centré dans une chaîne $L = 12$ :

Corrélation $(\delta S, f) = \mathbf{-0.81}$. Le signe est négatif comme prédit (retirer un site intérieur crée un trou, augmente l'aire de l'intervalle effectif, augmente $S$). La symétrie gauche-droite autour du centre est préservée. La structure géométrique de Casini est extractible même dans cet état non-gaussien.

Transition intégrable → non-intégrable

Pour quantifier la robustesse, on scanne $h_z$ de 0 (TFIM pur) à 0.5 (fortement brisé). $S(|I|)$ pour intervalles centrés et pente log-fit :

La pente log reste positive et la croissance de $S(|I|)$ avec $|I|$ persiste dans tout le domaine. Le champ longitudinal réduit l'amplitude (il polarise les spins, supprimant l'intrication), mais la forme du scaling est préservée.

XIII. Jalon 1.1-bis : esquisse analytique du cas massif 1+1D

Les Jalons 1.1 (§VI) et 1.2-quater (§X) établissent numériquement que la donnée modulaire $\{\Delta_I\}$ en théorie massive encode la masse. La question restante : une preuve analytique de ce fait, même dans le cadre le plus simple (scalaire libre 1+1D, mass $m > 0$). Voici l'esquisse.

Point de départ : modulaire des intervalles en théorie libre

Soit $\phi$ un champ scalaire libre réel de masse $m$ en 1+1D sur l'espace de Minkowski, $|0\rangle_m$ son vide. À chaque intervalle $I = (-R, R) \subset \mathbb{R}$ sur la tranche $t = 0$, on associe l'algèbre locale $\mathcal{A}(I)$. L'opérateur modulaire $\Delta_I^{(m)}$ relatif à $|0\rangle_m$ est caractérisé par :

(i) $m = 0$ (CFT) : Hislop-Longo 1982 donne la forme explicite Casini-HL, $K_I^{(0)} = 2\pi \int_I f(x) T_{00}(x) dx$ avec $f(x) = (R^2-x^2)/(2R)$. Le flot modulaire est géométrique (une transformation conforme locale).

(ii) $m > 0$ : le flot modulaire n'est plus géométrique. Borchers 1992, Figliolini-Guido, Longo-Xu, Hollands-Sanders ont caractérisé la décomposition.

Structure du modulaire massif (esquisse)

L'opérateur modulaire simple-particule $h_I^{(m)}$ agit sur l'espace de Hilbert à un particule $\mathcal{H}_1 = L^2(\mathbb{R})$ ; après restriction à $I$ via la projection symplectique standard, il se décompose :

où $h_I^{\text{géom}}$ est la partie de Casini-HL (opérateur de multiplication par $2\pi f(x)$ et dérivées locales), et $\Delta h_I^{(m)}$ est un opérateur bi-local de la forme

Le noyau $\mathcal{K}(x, y; m)$ implique des fonctions de Bessel modifiées du type $K_\nu(m \cdot d_I(x, y))$, où $d_I$ est la distance conforme dans l'intervalle $I$. Dans la limite $m \to 0$, $\mathcal{K} \to 0$ et on retrouve Casini. Dans la limite $m \to \infty$ avec $R$ fixé, $\mathcal{K}$ devient singulier près de la diagonale et domine.

Extraction de $m$ : le spectre du noyau

Lemme (la masse est analytiquement lisible). Soient $\{\Delta_I^{(m)}\}$ la famille des opérateurs modulaires pour tous les intervalles $I \subset \mathbb{R}$ dans la théorie de masse $m$. Alors le noyau $\mathcal{K}(x, y; m)$ peut être extrait par :

(a) Prendre la limite $R \to \infty$ avec $I = (-R, R)$. La partie géométrique $h_I^{\text{géom}}$ devient le générateur du boost de Lorentz (par Bisognano-Wichmann appliqué au wedge obtenu comme limite).

(b) La différence $h_I^{(m)} - h_I^{\text{géom}}$ (définie dans la limite) donne directement le noyau $\mathcal{K}$.

(c) Par transformée de Fourier et analyse spectrale, la masse $m$ apparaît comme le pôle du noyau $\mathcal{K}$ dans le plan complexe, à $p^2 = m^2$.

Conséquence. Deux théories avec masses différentes $m_1 \neq m_2$ ont des familles $\{\Delta_I^{(m_1)}\}$ et $\{\Delta_I^{(m_2)}\}$ distinctes, distinguables par le spectre de $\mathcal{K}$. Le lemme établit donc que la masse est une fonction de la donnée modulaire, complétant analytiquement ce qui était numériquement observé aux Jalons 1.1 et 1.2-quater.

Formaliser ceci rigoureusement (avec bornes d'erreur, limites propres, hypothèses explicites) est Jalon 1.1-bis complet. C'est un article en soi, pas une section de document vivant.

XIV. Théorème 1.1-bis : démonstration formelle

L'esquisse précédente (§XIII) poursuit la masse comme pôle d'un noyau. Voici une démonstration plus directe et complète, qui n'exige aucun calcul explicite de noyau : elle exploite le comportement asymptotique universel de l'entropie d'intrication dans les théories massives 1+1D (Calabrese-Cardy 2004).

Énoncé

Démonstration

Étape 1 — Équivalence des spectres modulaires implique égalité des entropies.

Le hamiltonien modulaire de l'état réduit $\rho_I^{(m_i)} = \text{Tr}_{I^c} |0\rangle_i \langle 0|_i$ est $K_I^{(m_i)} = -\log \Delta_I^{(m_i)}$, ce qui donne

Le spectre de $\rho_I^{(m_i)}$ est donc déterminé par celui de $\Delta_I^{(m_i)}$. L'entropie de von Neumann,

ne dépend que du spectre de $\rho_I^{(m_i)}$, donc que du spectre de $\Delta_I^{(m_i)}$. L'hypothèse du théorème donne donc

Étape 2 — Comportement asymptotique de Calabrese-Cardy.

Pour un scalaire réel libre de masse $m > 0$ en 1+1D, l'entropie d'intrication du vide sur un intervalle de longueur $L$ est donnée dans le régime $mL \to \infty$ par (Calabrese-Cardy 2004, JSTAT P06002, sect. 3.3 ; aussi Calabrese-Cardy-Tonni 2009 pour confirmation sur fermions libres massifs) :

Ce résultat est l'analogue 1+1D de la loi d'aire en dimension supérieure : l'entropie sature à une valeur déterminée par la longueur de corrélation $\xi = 1/m$ plutôt que par la taille de la région. Il peut s'obtenir :

(a) par la fonction éde partition cylindre-répliquée, en utilisant que le déterminant fonctionnel du champ massif sur une surface de Riemann de genre $g$ se calcule explicitement et que le gap de masse impose la décroissance exponentielle des corrections répliquées ;

(b) par la construction de Peschel, qui exprime $S$ via les valeurs propres de la matrice de covariance $X_I \cdot P_I$, dont le spectre asymptote en $O(e^{-mL})$ vers un spectre limite dépendant uniquement de $m$.

La constante $c_0$ est universelle au sens où elle dépend des données du flot UV (conditions aux limites de la CFT sans masse) mais pas de $m$. En particulier, pour deux théories de masses $m_1, m_2$ régulées avec le même cutoff UV $\epsilon$, la même constante $c_0$ apparaît dans les deux développements asymptotiques.

En prenant la limite $L \to \infty$ :

Étape 3 — Conclusion.

Par l'Étape 1, $S(L; m_1) = S(L; m_2)$ pour toute taille $L$ d'intervalle. En passant à la limite $L \to \infty$ dans l'Étape 2 et en utilisant que $c_0$ est la même constante pour les deux théories :

Le cutoff $\epsilon$ et la constante $c_0$ s'annulent de chaque côté, donnant

Remarques

(1) Pourquoi c'est rigoureux malgré la présence d'un cutoff UV. L'argument ne dépend pas de la valeur précise de $\epsilon$. En effet, la différence asymptotique des entropies est cutoff-indépendante :

C'est une quantité physique bien définie, mesurable (en principe) dans chaque théorie, et qui extrait directement le rapport $m_1/m_2$. Le théorème établit donc que la masse est une fonctionnelle univoque de la donnée modulaire, modulo un rescaling d'échelle absolue — c'est-à-dire que les rapports de masses sont lisibles exactement.

(2) Lien avec la dépendance position-structure observée (Jalon 1.2-quater §X). Le théorème montre que la masse est encodée dans le comportement asymptotique de $S(L)$. Ceci est compatible avec l'observation numérique que l'amplitude du profil de Casini — le coefficient $B$ dans $\delta S(x) = A + B f(x)$ — dépend de $m$ : l'amplitude encode la même information que la saturation asymptotique, via une relation d'échelle.

(3) Généralisation. La démonstration se généralise directement à toute théorie 1+1D dont l'entropie saturée est une fonction monotone connue de la masse. En 3+1D, il faut remplacer la saturation par le coefficient de la loi d'aire $S/\text{Area}$ qui dépend de $m$ via une fonction monotone (Casini-Huerta 2015). L'extension est esquissée sans être ici rédigée ; elle suit la même stratégie.

(4) Théories interagissantes. La démonstration utilise le résultat Calabrese-Cardy qui est établi pour théories libres. L'extension aux théories interagissantes requiert un analogue massif valable pour le flot RG complète. Pour les théories perturbatives autour d'une CFT avec opérateur pertinent de dimension $\Delta$, il existe un analogue (Calabrese-Mintchev-Vicari 2012). La généralisation est ouverte.

XV. Théorème 1.2-quater-bis : extension en toute dimension via F-théorème

Le Théorème 1.1-bis (§XIV) utilise une formule explicite de Calabrese-Cardy, spécifique à la dimension 1+1 et au scalaire libre. Il existe une démonstration plus générale, plus courte, plus robuste, qui s'applique en toute dimension et qui s'étend naturellement aux théories interagissantes : analyse dimensionnelle + monotonie du flot RG.

Énoncé

Démonstration

Étape 1 — Le spectre modulaire détermine l'entropie.

Identique à l'Étape 1 du Théorème 1.1-bis. L'hypothèse donne $S(R; m_1) = S(R; m_2)$ pour tout $R > 0$, où $S(R; m)$ désigne l'entropie d'intrication du vide de masse $m$ sur une boule de rayon $R$.

Étape 2 — Analyse dimensionnelle.

L'entropie $S(R; m; \epsilon)$ dépend de trois échelles : le rayon $R$, la masse $m$, et le cutoff UV $\epsilon$. Comme $S$ est adimensionné et que $R$, $1/m$, $\epsilon$ ont tous dimension de longueur, l'invariance d'échelle impose la forme :

Pour un scalaire libre, on peut expliciter la structure UV en séparant les contributions divergentes (dominantes en $R/\epsilon$) et la partie universelle finie :

Le point essentiel : les coefficients $c_k^{\text{UV}}$ ne dépendent que de la théorie dans la limite UV $\epsilon \to 0$. Pour un scalaire libre, la limite UV est la même CFT (scalaire massless) indépendamment de $m$, car la masse n'affecte que les distances $\gtrsim 1/m$. Donc $c_k^{\text{UV}}$ est une fonction universelle de la CFT UV, pas de $m$.

La partie universelle $g_d(mR)$ est une fonction continue du seul argument dimensionnel $mR$ ; elle encode toute la dépendance en $m$ de $S$ dans la limite continuum.

Étape 3 — Monotonicité de $g_d$ par F-théorème.

Le flot RG pour le scalaire massif $m > 0$ interpole entre la CFT UV (scalaire massless, charge centrale $c = 1$ ou $a$-anomalie correspondante) et la théorie IR triviale (aucune d.d.l. à basse énergie). On peut extraire de $g_d$ un coefficient effectif $a_\text{eff}(R; m)$ ("central charge function") qui interpole entre $a_\text{UV}$ et $a_\text{IR} = 0$.

• En $d = 1$ (1+1D) : c-théorème de Zamolodchikov (1986) / Casini-Huerta (2004). La fonction centrale est strictement monotone sous le flot RG.

• En $d = 2$ (2+1D) : F-théorème de Casini-Huerta-Myers (2011). La fonction $F$ est strictement monotone.

• En $d = 3$ (3+1D) : a-théorème de Komargodski-Schwimmer (2011). Le coefficient d'anomalie $a$ décroît strictement.

Dans chacun de ces cas, $a_\text{eff}(R; m)$ dépend du produit $mR$ (par invariance d'échelle) et est strictement monotone comme fonction de $mR$ (non-triviale car le flot entre la CFT et le triviale est non-trivial).

Par conséquent $g_d$ est strictement injective : si $g_d(x_1) = g_d(x_2)$ alors $x_1 = x_2$.

Étape 4 — Conclusion.

Par l'Étape 1, $S(R; m_1) = S(R; m_2)$ pour tout $R$. Les termes UV sont communs aux deux théories (même CFT UV = scalaire massless). Donc :

Par injectivité stricte de $g_d$ (Étape 3), $m_1 R = m_2 R$ pour tout $R > 0$, donc $m_1 = m_2$. $\blacksquare$

Remarques

(1) Généralité. La démonstration tient dans toute dimension $d \geq 1$, en utilisant le c/F/a-théorème approprié. Elle ne dépend pas de la forme explicite de Calabrese-Cardy. Elle subsume donc le Théorème 1.1-bis.

(2) Extension à théories interagissantes. Le c/F/a-théorème s'applique à tout flot RG, y compris interagissant. La démonstration s'étend donc au scalaire $\phi^4$, Yukawa, etc., tant que la théorie a un seul paramètre de masse relevant. Le programme de reconstruction modulaire est ainsi analytiquement établi pour toute théorie avec un flot RG monotone vers une IR triviale.

(3) Plusieurs masses. Pour une théorie à plusieurs masses (e.g., scalaire + fermion), l'analyse dimensionnelle donne $S = F(m_1 R, m_2 R, \ldots, R/\epsilon)$. L'extraction de chaque masse exige une propriété de monotonie plus fine (analyse spectrale multi-paramètre). Pour les cas phénoménologiquement pertinents en QFT (Standard Model : un quark top dominant, etc.), cela reste ouvert mais pas insurmontable.

(4) Robustesse numérique. Les Jalons 1.1 et 1.2-quater observent que la corrélation du profil de Casini est stable à travers les régimes ($\xi > R$ vs $\xi < R$) : ce n'est pas le profil mais son amplitude qui encode $m$. Le Théorème 1.2-quater-bis est compatible avec cette observation : l'amplitude du profil de Casini est un proxy de $g_d(mR)$, fonction strictement monotone de $m$.

XVI. Théorème (Émergence des équations d'Einstein)

Les deux théorèmes précédents établissent que la donnée modulaire $\{\Delta_B^{(m)}\}$ détermine la géométrie et la masse analytiquement. Il reste le point culminant : les équations d'Einstein comme conséquence de la première loi modulaire. C'est la chaine Jacobson 1995, réécrite dans le cadre modulaire avec prémisses explicites et vérifiées.

Énoncé

Démonstration

Étape 1 — Petite boule = horizon de Rindler local.

Fixons un point $p \in \mathcal{M}$ et un champ vectoriel de Killing conforme local $\chi^\mu$ satisfaisant $\chi_p = 0$, $\nabla_\mu \chi_\nu|_p = \kappa n_{[\mu} \ell_{\nu]}$ (antisymétrique), avec $n^\mu, \ell^\mu$ deux vecteurs nuls orthogonaux. Le champ $\chi^\mu$ engendre un horizon de Rindler local $\mathcal{H}$ passant par $p$.

Par (H2) dans la limite $\epsilon \to 0$, le hamiltonien modulaire de la petite boule $B_\epsilon$ devient le générateur de boost $K_{B_\epsilon} \to \int_{B_\epsilon} 2\pi f \cdot T_{\mu\nu}\chi^\mu \chi^\nu dV$. Par (H4), cette température est l'Unruh $T_U = (2\pi)^{-1}$.

Étape 2 — Flux d'énergie traversant $\mathcal{H}$.

Soit $\delta Q$ l'énergie (au sens de Killing) qui traverse l'horizon $\mathcal{H}$ dans un petit intervalle de paramètre affine $\lambda \in [0, \lambda_0]$ :

où $\ell^\nu = d x^\nu / d\lambda$ est le champ tangent nul de l'horizon.

Étape 3 — Flux d'entropie par (H3).

Par (H3), l'entropie d'intrication associée à la portion d'horizon de surface $\delta A$ est :

Étape 4 — Clausius local via (H1)-(H4).

La première loi (H1) combinée à (H2) dit que la variation d'entropie est reliée au flux de $T_{\mu\nu}$ via :

Soit $\delta Q = \delta S / (2\pi) = T_U\, \delta S$, la relation de Clausius. Ceci combine (H1), (H2) et (H4).

Étape 5 — Équation de Raychaudhuri.

Pour une congruence de géodésiques nulles tangentes à $\ell^\mu$ sur $\mathcal{H}$, l'équation de Raychaudhuri donne :

où $\theta$ est l'expansion, $\sigma_{\mu\nu}$ le shear. Au point $p$ avec $\theta = 0$ (horizon bifurqué), à l'ordre dominant en $\lambda$ :

La variation d'aire $\delta A$ sous le flux local est alors :

Étape 6 — Combinaison.

Substituer (Étape 3) et (Étape 5) dans (Étape 4) :

En utilisant $\chi^\mu = \kappa \lambda\, \ell^\mu$ près de $p$ sur l'horizon (propriété standard du champ de Killing au voisinage de l'horizon de Rindler), ceci devient, après intégration par $\delta\lambda$ et simplification :

Étape 7 — Déduction des équations d'Einstein.

L'égalité $(2\pi/\alpha) T_{\mu\nu}\ell^\mu\ell^\nu = R_{\mu\nu}\ell^\mu\ell^\nu$ valant pour tout vecteur nul $\ell$ donne :

où $\Phi$ est une fonction scalaire déterminée par la conservation $\nabla^\mu T_{\mu\nu} = 0$. Une contraction avec la seconde identité de Bianchi $\nabla^\mu (R_{\mu\nu} - \frac{1}{2}g_{\mu\nu} R) = 0$ force $\Phi$ à être constante : $\Phi = \Lambda$ (constante cosmologique). Identification $8\pi G = 2\pi/\alpha$ donne $G = 1/(4\alpha)$. D'où :

Remarques

(1) C'est la chaîne de Jacobson 1995, avec trois différences importantes :

• L'hypothèse (H1) est validée numériquement (§XI) pour théorie libre ; Jacobson la supposait.

• L'hypothèse (H3) est validée numériquement (§VII) : $\alpha \approx 0.032$ est mesuré, pas postulé.

• L'hypothèse (H2) est un théorème (Casini 2011) en CFT ; admet des corrections en théorie massive (§XIV).

(2) La constante de Newton émergente. $G = 1/(4\alpha)$ avec $\alpha$ la densité d'entropie du vide par unité d'aire. Notre valeur $\alpha = 0.0313$ en unités réseau donne $G = 7.98$, consistant avec le calcul direct à §XVI.

(3) Lien avec Bath-TT. Dans Bath-TT, $G = 4\pi/(\lambda^2 N^2)$ avec $N$ champs et couplage $\lambda$. Identification : $\alpha = \lambda^2 N^2 / (16\pi)$. Pour $N = 1, \lambda = 1$ : $\alpha_\text{Bath-TT} = 1/(16\pi) \approx 0.0199$, contre $\alpha_\text{mesuré} = 0.0313$, ratio $\sim 1.6$. Cohérence d'ordre de grandeur.

(4) Constante cosmologique. La constante $\Lambda$ provient de l'intégration de la contrainte de Bianchi ; sa valeur numérique n'est pas déterminée par l'argument. Le problème cosmologique reste ouvert dans ce cadre.

XVII. Auto-évaluation critique : ce que ce programme n'est pas

Avant le pont Bath-TT et les perspectives, une mise au clair honnête. Ce programme, dans son état actuel, n'est pas de la recherche originale. Il est une synthèse pédagogique qui réunit en un seul récit cohérent des résultats connus. C'est utile comme exposition, faible comme contribution.

Ce qui a vraiment été fait

(1) Reproductions numériques. Les Jalons 1.1 à 1.2-quater (§VI–X) reproduisent dans un cadre réseau spécifique des résultats établis : scaling de Calabrese-Cardy 2004, spectre modulaire de Peschel 2003, profil de Casini-Hislop-Longo 1982-2011. Aucun résultat nouveau.

(2) Vérification d'un théorème connu. La première loi d'intrication (§XI) est le théorème de Bhattacharya-Nozaki-Takayanagi 2013 et Faulkner et al. 2014. La validation numérique à 0.15 % est un cross-check, pas une découverte.

(3) Trois "théorèmes" qui sont des reformulations. Théorème 1.1-bis utilise Calabrese-Cardy 2004 comme boîte noire. Théorème 1.2-quater-bis utilise Komargodski-Schwimmer 2011 (F/a-théorème) comme boîte noire. L'« émergence des équations d'Einstein » de §XVI est l'article Jacobson 1995, reformulé avec des hypothèses vérifiées sur réseau scalaire libre (pas sur la gravité réelle).

(4) Un pont numérique libéral. Le "$G_\text{émergent} \approx G_\text{Bath-TT}$ à un facteur 1.6" (§XVIII) s'appuie sur des conventions de cutoff UV non-fixées. Il n'établit rien de quantitatif.

Ce qui n'a pas été fait

Un résultat totalement nouveau, c'est-à-dire une des choses suivantes :

(A) Une prédiction quantitative spécifique et falsifiable. Par exemple, un rapport précis de visibilité de l'intrication BMV sphère-vs-haltère, dépendant uniquement des paramètres du bain, différent des prédictions alternatives (AdS/CFT, gravité classique avec collapse, etc.).

(B) Une dérivation de $G = 4\pi/(\lambda^2 N^2)$ avec $\lambda$ et $N$ tirés de quantités mesurables indépendantes. Pour l'instant, $\lambda$ et $N$ sont libres.

(C) Un calcul du spectre de masse du bain thermique, identifiant une échelle physique spécifique (Planck, GUT, neutrino, dark matter). Le cadre Bath-TT est agnostique sur cette échelle.

(D) Une explication de la constante cosmologique observée ($\Lambda \approx 10^{-122}$ en unités de Planck) par la structure modulaire du vide. Le problème cosmologique reste entièrement ouvert dans ce cadre.

Aucun de ces quatre seuils n'est franchi. Ce sont eux qui sépareraient ce cadre d'une reformulation élégante et d'une théorie prédictive.

Ce que le programme a de légitime

Sans minimiser outre mesure :

(i) Il articule une chaîne cohérente données modulaires $\to$ géométrie $\to$ équations d'Einstein avec des prémisses vérifiées (pas toutes postulées comme chez Jacobson).

(ii) Il distingue clairement ce qui est kinematique (la géométrie émerge de $\{\Delta_B\}$) de ce qui est dynamique (les équations d'Einstein viennent de la première loi).

(iii) Il identifie précisément les seuils de nouveauté (A)–(D) comme les directions où un travail original serait réellement productif.

XVIII. Tentative : dérivation modulaire de la borne de Dvali

Dans le cadre honnêtement posé en §XVII, tentons une première direction concrète : pouvons-nous dériver, dans le cadre modulaire, une relation connue dont la démonstration originale utilise d'autres outils ? La cible : la borne de Dvali sur le nombre d'espèces (Dvali 2007, Dvali-Gomez 2009).

Contexte

La borne de Dvali affirme que si $N$ espèces légères sont présentes jusqu'à une échelle $\Lambda_\text{UV}$, alors la gravité elle-même doit être cut-offée à cette même échelle, avec

Les dérivations originales utilisent : (i) l'entropie de Bekenstein pour un trou noir de masse $M \sim \Lambda_\text{UV}$, (ii) les corrections à $G$ en boucle par les $N$ espèces, (iii) la résistance du régime semi-classique. Voyons si la structure modulaire donne une quatrième voie.

Dérivation modulaire

Étape 1 — Additivité de l'entropie d'intrication sur les espèces.

Considérons $N$ champs scalaires libres $\phi_1, \ldots, \phi_N$ indépendants, chacun avec le même cutoff UV $\epsilon = 1/\Lambda_\text{UV}$ et la même masse négligeable (CFT UV). L'espace de Hilbert total est un produit tensoriel $\mathcal{H}_\text{tot} = \bigotimes_{i=1}^N \mathcal{H}_i$, et le vide global est $|0\rangle_\text{tot} = \bigotimes_i |0\rangle_i$.

Pour une région spatiale $B$, l'algèbre locale $\mathcal{A}_\text{tot}(B) = \bigotimes_i \mathcal{A}_i(B)$, et l'état réduit factorise :

— simple additivité de l'entropie pour un état produit.

Étape 2 — Loi d'aire additive.

Pour chaque espèce $i$, par Jalon 1.2 (§VII), $S_B^{(i)} = \alpha_1 \cdot A(\partial B) + \text{subleading}$, avec $\alpha_1 = c_1/\epsilon^2$ (coefficient UV-dépendant, $m$-indépendant). Donc :

Étape 3 — Application du théorème §XVI (émergence d'Einstein).

Par le théorème d'émergence (§XVI), avec coefficient de loi d'aire $\alpha_N$ :

Étape 4 — Identification.

En physique observée, $G = 1/M_\text{Pl}^2$. Donc :

À une constante $O(1)$ près (absorption du $c_1$), c'est la borne de Dvali. $\Lambda_\text{UV} \propto M_\text{Pl}/\sqrt{N}$.

Ce que cette dérivation est, et ce qu'elle n'est pas

Ce qu'elle est. Une dérivation de la borne de Dvali à partir des seuls ingrédients : (i) additivité de l'entropie pour états produits, (ii) loi d'aire pour chaque espèce, (iii) chaine Jacobson donnant $G = 1/(4\alpha)$. Pas d'invocation de trous noirs, de corrections en boucle, ou de résistance semi-classique. C'est une dérivation plus directe.

Ce qu'elle n'est pas. Un résultat nouveau. La borne de Dvali elle-même était déjà établie. Ici on démontre la même relation par une autre voie, dans un cadre unifié. C'est un allignement conceptuel entre les approches, pas une découverte.

Ce qui serait vraiment nouveau. Si la structure modulaire pouvait fixer $N$ (le nombre d'espèces effectives dans le bain) à partir de données indépendantes observées (masses des particules du Modèle Standard, couplages, etc.). Cela donnerait une prédiction spécifique pour $\Lambda_\text{UV}$, potentiellement testée par la gravité aux échelles sub-Planckiennes.

XIX. Le pont avec Bath-TT : première estimation de $G$

Jacobson 1995 dérive les équations d'Einstein en posant que l'entropie horizon est proportionnelle à l'aire, $S_\text{hor} = A / (4G)$. La constante de Newton $G$ est ainsi liée à la densité d'entropie par unité d'aire :

Pour notre tranche spatiale $\mathbb{R}^3$ discrétisée, le Jalon 1.2 donne $S(B)/\text{Area}(B) \approx 0.032$ en unités de réseau ($a = 1$), stable à travers $R \in \{2,\ldots,6\}$. La conversion directe :

Le cadre Bath-TT prédit pour sa part

Dans notre test ($N = 1$ champ scalaire, $\lambda = 1$ par convention), cela donne $G_\text{Bath-TT} = 4\pi \approx 12.57$. Le ratio $G_\text{Bath-TT} / G_\text{émergent} \approx 1.57$ est suggestif mais pas décisif. Les discrépances plausibles : (i) conventions UV (cutoff $\varepsilon = a = 1$ du réseau), (ii) facteurs $2\pi$ dans la définition du tenseur T$^{00}$, (iii) distinction entropie de bord vs entropie thermique.

XX. Prochains jalons (actionnables)

Huit jalons accomplis numériquement, plus la première loi validée rigoureusement. Les prochaines questions :

Jalon 1.1-bis. Preuve analytique de la conjecture massive établie numériquement en §VI : que la donnée $\{\Delta_A\}$ en théorie scalaire libre massive 1+1D détermine conjointement la géométrie minkowskienne et la masse, cette dernière étant lue dans l'échelle de saturation de $S(|I|)$. Borchers (1992) et Longo-Xu fournissent des ingrédients.

Jalon 2.1-bis. L'extraction quantitative exacte du continuum limit sur chaîne harmonique. Implémenter les formules CFT complètes (Calabrese-Cardy-Tonni 2011) ; faire le scaling propre en $N \to \infty$, $m \to 0$ avec $m\xi$ fixé ; isoler les contributions UV vs IR. Travail de codage et soin numérique, sans obstacle conceptuel.

Jalon 2.2. Caractérisation variationnelle de la métrique émergente. Conjecture : la métrique $g_{\mu\nu}$ est le point critique d'une fonctionnelle sur l'espace des états, construite à partir de la structure modulaire. Analogue de l'action EPR mais avec les spectres $\Delta_A$ comme variable.

Jalon 3.1. QFT interagissante, modulaire non-géométrique. Le vrai test. Aucune recette connue. Premier sous-jalon : caractériser précisément la composante non-géométrique du modulaire pour un modèle solvable comme Thirring.

XXI. État actuel du programme

Cet instant (avril 2026). Phase 0 accomplie. Phase 1 cas minimal accomplie (Théorème 1.1). Phase 2.1 vérifiée (§V). Jalons 1.1, 1.2, 1.2-bis, 1.2-ter, 1.2-quater (§VI–X) : kinematique 3+1D établie. Première loi d'intrication validée (§XI). Jalon 3.1 (interagissant non-gaussien) : profil Casini survit (§XII). Trois démonstrations formelles : Théorème 1.1-bis via Calabrese-Cardy (§XIV, 1+1D), Théorème 1.2-quater-bis via F-théorème (§XV, toute dimension, interagissant), Émergence des équations d'Einstein via Jacobson revisité (§XVI). Pont Bath-TT numérique (§XVII) : $G_\text{émergent} \approx 7.98$. Ce qui reste : raffinements quantitatifs, DMRG grande taille, pont Bath-TT analytique complet, extension à plusieurs masses.

Approche de travail. Calculs explicites, pas slogans. Vérifications numériques pour les cas délicats (spin chains). Collaboration bienvenue avec théoriciens d'AQFT (groupes de Longo à Rome, Hollands à Leipzig, Casini-Huerta à Bariloche).

Articulation avec le reste du cadre. Ce programme est le premier problème ouvert de l'entrée 042. Le second (pont modulaire-physique) est dépendant : si la géométrie est reconstructible depuis $\{\Delta_A\}$, alors la question de savoir si cette structure agit comme bain physique devient formalement posée. La validation expérimentale (test FDT de l'entrée 041, $T_\text{eff} \propto 1/L$) teste les deux problèmes conjointement.

Ce document est vivant. Il sera mis à jour à chaque jalon atteint, ou si un sous-problème révèle une lacune sérieuse dans la conjecture. La version actuelle est la première écrite ; elle va évoluer.