2025-11-21T02:34:15.165429

Maximal Entropy Random Walks in Z: Random and non-random environments

Thibaut, Gerin, Offret
The Maximal Entropy Random Walk (MERW) is a natural process on a finite graph, introduced a few years ago with motivations from theoretical physics. The construction of this process relies on Perron-Frobenius theory for adjacency matrices. Generalizing to infinite graphs is rather delicate, and in this article, we treat in a fairly exhaustive manner the case of the MERW on Z with loops, for both random and nonrandom loops. Thanks to an explicit combinatorial representation of the corresponding Perron-Frobenius eigenvectors, we are able to precisely determine the asymptotic behavior of these walks. We show, in particular, that essentially all MERWs on Z with loops have positive speed.
academic

Marches Aléatoires d'Entropie Maximale dans Z : Environnements aléatoires et non-aléatoires

Informations Fondamentales

  • ID de l'article: 2503.15957
  • Titre: Maximal Entropy Random Walks in Z: Random and non-random environments
  • Auteurs: Thibaut Duboux (Université Bourgogne Europe), Lucas Gerin (École Polytechnique), Yoann Offret (Université Bourgogne Europe)
  • Classification: math.CO (Mathématiques Combinatoires), math.PR (Théorie des Probabilités)
  • Date de soumission: 20 novembre 2025
  • Lien de l'article: https://arxiv.org/abs/2503.15957

Résumé

Les marches aléatoires d'entropie maximale (MERW) constituent un processus naturel sur les graphes finis, introduit il y a quelques années par des motivations issues de la physique théorique. La construction de ce processus repose sur la théorie de Perron-Frobenius de la matrice d'adjacence. La généralisation aux graphes infinis s'avère subtile. Cet article étudie en détail le modèle MERW sur la grille entière Z avec boucles, couvrant les environnements aléatoires et non-aléatoires. Grâce à une représentation combinatoire explicite des vecteurs propres de Perron-Frobenius correspondants, les auteurs parviennent à déterminer précisément le comportement asymptotique de ces marches. En particulier, il est démontré que presque toutes les MERW avec boucles sur Z possèdent une vitesse positive.

Contexte et Motivation de la Recherche

Contexte du Problème

  1. Limitations des marches aléatoires standards: Les marches aléatoires standards constituent un processus stochastique fondamental en théorie des probabilités, physique statistique et analyse de réseaux, avec probabilités de transition données par pi,j=ai,j/ai,p_{i,j} = a_{i,j}/\sum_\ell a_{i,\ell}. Cependant, ces marches ne maximisent pas nécessairement l'entropie des trajectoires.
  2. Motivation physique des MERW: Les MERW proviennent de la formulation par intégrales de chemin de la mécanique quantique, avec probabilités de transition pi,j=ai,jψj/(λψi)p_{i,j} = a_{i,j}\psi_j/(\lambda\psi_i), où ψ\psi est le vecteur propre positif correspondant au rayon spectral λ\lambda. Cette construction permet aux MERW de maximiser le taux d'entropie des trajectoires sur les graphes finis.
  3. Défis pour les graphes infinis:
    • Sur les graphes infinis, la théorie de Perron-Frobenius ne s'applique plus directement
    • Il est nécessaire de distinguer les cas R-récurrent et R-transient
    • Dans le cas R-transient, le vecteur propre positif n'est pas unique, et il existe un ensemble convexe de solutions extrémales
  4. Lacune de recherche: Bien que les MERW sur les graphes finis aient été largement étudiées, il existe un manque d'études systématiques pour les graphes infinis, en particulier dans les cas avec perturbations aléatoires.

Importance de la Recherche

  1. Signification théorique: Établit des connexions entre la probabilité combinatoire, les processus stochastiques et la théorie de la localisation d'Anderson
  2. Valeur applicative: Perspectives d'application dans les algorithmes de détection de communautés dans les réseaux complexes
  3. Contribution méthodologique: Fournit de nouveaux outils pour analyser les MERW sur les graphes infinis

Contributions Principales

L'article apporte les contributions majeures suivantes :

  1. Représentation explicite des vecteurs propres (Théorème 2.9): Pour les environnements agréables, une description combinatoire complète et explicite des deux vecteurs propres extrémaux ψ+\psi^+ et ψ\psi^- est fournie :\beta_{-1}\cdots\beta_i, & i < 0\\ 1, & i = 0\\ (\beta_0\beta_1\cdots\beta_{i-1})^{-1}, & i > 0 \end{cases}$$ où $\beta_i$ et $\alpha_i$ peuvent être exprimés par des fonctions génératrices ou des fractions continues.
  2. Caractérisation complète pour les environnements déterministes (Propositions 2.3, 2.5):
    • Démonstration que la matrice AA est R-transiente
    • Description précise de la monotonie et du comportement limite des vecteurs propres extrémaux
    • Preuve que toutes les MERW dans les environnements agréables sont transitoires
  3. Vitesse linéaire en environnement aléatoire (Théorème 3.1): Pour les environnements aléatoires i.i.d., les MERW extrémales possèdent une vitesse constante positive : limnXn+n=v>0(presque suˆrement)\lim_{n\to\infty} \frac{X^+_n}{n} = v > 0 \quad \text{(presque sûrement)}v=1/Eμ[S]v = 1/\mathbb{E}_\mu[S], avec SS ayant une expression explicite.
  4. Analyse de couplage pour les MERW non-extrémales (Théorème 3.4): Pour les MERW correspondant aux vecteurs propres mixtes ψ(κ)=κψ++(1κ)ψ\psi^{(\kappa)} = \kappa\psi^+ + (1-\kappa)\psi^-, il est démontré qu'une vitesse linéaire vv est également atteinte sous la condition Xn(κ)+X^{(\kappa)}_n \to +\infty.
  5. Calculs exacts pour l'environnement de Bernoulli (Proposition 3.6): Lorsque l'environnement satisfait P(wk=M)=pP(w_k=M)=p, il est démontré que : limp0vp,M=14(2+M)2>0\lim_{p\to 0} v_{p,M} = \sqrt{1-\frac{4}{(2+M)^2}} > 0 Ceci montre que la vitesse présente un saut de discontinuité lorsque p0p\to 0.
  6. Contraste avec les environnements périodiques (Proposition A.1): Il est démontré que pour les environnements déterministes périodiques, les MERW sont nullement récurrentes et hautement localisées, formant un contraste frappant avec les environnements aléatoires.

Détails Méthodologiques

Définition de la Tâche

Entrées:

  • Environnement w=(wi)iZw = (w_i)_{i\in\mathbb{Z}}, où wi0w_i \geq 0 est le poids de la boucle au sommet ii
  • Matrice d'adjacence AA définie par : poids unitaire entre sommets adjacents, poids de boucle wiw_i au sommet ii

Sorties:

  • Rayon spectral combinatoire λ\lambda
  • Vecteur propre positif ψ\psi satisfaisant Aψ=λψA\psi = \lambda\psi
  • Comportement asymptotique de la marche aléatoire MERW correspondante (Xn)n0(X_n)_{n\geq 0}

Contraintes:

  • L'environnement est M-agréable : borné, non identiquement égal à MM, et pour tout ε>0\varepsilon>0 et r0r\geq 0, il existe ii tel que wi,,wi+rw_i,\ldots,w_{i+r} soient tous Mε\geq M-\varepsilon

Architecture de la Méthode Principale

1. Calcul du Rayon Spectral Combinatoire (Lemme 2.2)

La preuve que λ=2+M\lambda = 2+M est donnée par comptage de chemins :

Stratégie de preuve de la borne inférieure:

  • Considération du nombre de chemins de ii à ii restant au-dessus de ii : uriiu^{i\circlearrowleft i}_r
  • Utilisation de la fonction génératrice Hi,i[i],w(z)=r0dr00zrH^{[\geq i],w}_{i,i}(z) = \sum_{r\geq 0} d^{0\circlearrowleft 0}_r z^r
  • Décomposition en arches (arch-decomposition) donnant : Hi,i[i],w(z)=11(Mε)zz2Hi,i[i],w(z)H^{[\geq i],w}_{i,i}(z) = \frac{1}{1-(M-\varepsilon)z - z^2H^{[\geq i],w}_{i,i}(z)}
  • Résolution donnant la singularité principale en z=(2+Mε)1z^* = (2+M-\varepsilon)^{-1}
  • Utilisation du théorème de transfert pour obtenir dr00c(2+Mε)rr3/2d^{0\circlearrowleft 0}_r \geq c(2+M-\varepsilon)^r r^{-3/2}

2. Construction des Vecteurs Propres Extrémaux (Proposition 2.3)

Méthode d'approximation par troncature:

  • Pour k0k\leq 0, définition de ψ(k,ε)\psi^{(k,\varepsilon)} satisfaisant les conditions aux limites : 0, & n\leq k-2\\ \varepsilon, & n=k-1\\ 2\varepsilon, & n=k \end{cases}$$ et la relation de récurrence $\psi_{n+1} + w_n\psi_n + \psi_{n-1} = \lambda\psi_n$
  • Ajustement de εk\varepsilon_k de sorte que ψ0(k,εk)=1\psi^{(k,\varepsilon_k)}_0 = 1
  • Utilisation d'un argument diagonal pour obtenir une sous-suite convergente
  • Exploitation de l'estimée de convexité pour prouver limnψn+=\lim_{n\to\infty}\psi^+_n = \infty

3. Dérivation de la Représentation Combinatoire (Théorème 2.9)

Étapes clés:

Définition des fonctions génératrices : βi=1λHi,i[i],w(1λ),αi=1λHi,i[i],w(1λ)\beta_i = \frac{1}{\lambda}H^{[\leq i],w}_{i,i}\left(\frac{1}{\lambda}\right), \quad \alpha_i = \frac{1}{\lambda}H^{[\geq i],w}_{i,i}\left(\frac{1}{\lambda}\right)

Obtention par décomposition en arches des relations de récurrence : βi=1λwiβi1,αi=λwi1αi+1\beta_i = \frac{1}{\lambda - w_i - \beta_{i-1}}, \quad \alpha_i = \lambda - w_i - \frac{1}{\alpha_{i+1}}

Ceci est équivalent à l'expansion en fraction continue : αi=1λwi1λwi+11\alpha_i = \cfrac{1}{\lambda - w_i - \cfrac{1}{\lambda - w_{i+1} - \cfrac{1}{\ddots}}}

Vérification que ψi+=β1βi\psi^+_i = \beta_{-1}\cdots\beta_i (pour i<0i<0) satisfait effectivement l'équation aux valeurs propres.

4. Analyse de la Vitesse en Environnement Aléatoire

Identification de la chaîne de Markov:

  • Lorsque (wi)(w_i) est i.i.d., (βi)i(\beta_i)_i et (αi)i(\alpha_{-i})_i constituent des chaînes de Markov stationnaires ergodiques
  • Les MERW extrémales constituent des marches aléatoires en environnement aléatoire ergodique

Formule de vitesse (Équation 3.4): v=1Eμ[S],S=n0P0w(Xn+=0)v = \frac{1}{\mathbb{E}_\mu[S]}, \quad S = \sum_{n\geq 0} P^w_0(X^+_n = 0)

Preuve de finitude: Utilisation du contrôle : βi1{wiMδ}gMδ(βi1)+1{wi>Mδ}\beta_i \leq 1_{\{w_i\leq M-\delta\}} g_{M-\delta}(\beta_{i-1}) + 1_{\{w_i>M-\delta\}}

Obtention de : Eμ[β0β12βi2]Eμ[Z02]i<1\mathbb{E}_\mu[\beta_0\beta^2_{-1}\cdots\beta^2_{-i}] \leq \mathbb{E}_\mu[Z^2_0]^i < 1

Par conséquent, Eμ[S]<\mathbb{E}_\mu[S] < \infty, garantissant une vitesse positive.

Points d'Innovation Technique

  1. Méthode des fonctions génératrices de chemins: Transformation du problème de vecteur propre en problème de comptage de chemins, utilisant les outils de la combinatoire analytique.
  2. Technique des variables de Riccati: Les relations de récurrence satisfaites par αi\alpha_i et βi\beta_i sont liées aux variables de Riccati de la théorie de la localisation d'Anderson, établissant un lien entre les MERW et les opérateurs de Schrödinger aléatoires.
  3. Argument de couplage: Pour les MERW non-extrémales, construction d'un couplage avec les MERW extrémales, utilisant la finitude de l'ensemble des « mauvais temps » pour prouver la cohérence des vitesses.
  4. Analyse des frontières des fonctions génératrices: Contrôle précis du comportement des fonctions génératrices à la singularité principale, permettant l'obtention d'estimées asymptotiques du nombre de chemins.

Configuration Expérimentale

Paramètres de Simulation Numérique

Bien que cet article soit principalement théorique, il inclut des vérifications numériques :

  1. Figure 1: Simulation de MERW en environnement de Bernoulli
    • Paramètres : M=20M=20, p=0.02p=0.02 et p=0.05p=0.05
    • 200 trajectoires indépendantes, étapes temporelles n=600n=600
    • Vérification que la vitesse vp,Mv_{p,M} décroît en fonction de pp
  2. Figure 3: Variation de la vitesse vp,Mv_{p,M} en fonction de pp
    • Simulation Monte Carlo de la formule (3.21)
    • Différentes valeurs de MM : M=0.1,1,10M=0.1, 1, 10
    • Vérification du comportement asymptotique de la Proposition 3.6
  3. Figure 4: Variation de la vitesse vp,Mv_{p,M} en fonction de MM
    • Paramètres : p=0.3,0.5,0.8p=0.3, 0.5, 0.8
    • Démonstration que la vitesse n'est pas monotone en fonction de MM

Vérification Théorique

Calculs exacts pour l'environnement de Bernoulli:

  • Utilisation d'identités combinatoires liées aux nombres de Motzkin
  • Vérification que limp0vp,M=14/(2+M)2\lim_{p\to 0} v_{p,M} = \sqrt{1-4/(2+M)^2}
  • Preuve que vp,M3(1p)/(2+M)v_{p,M} \sim 3(1-p)/(2+M) lorsque p1p\to 1

Résultats Expérimentaux

Résultats Principaux

1. Environnement Déterministe (Section 2)

Vérification du Théorème 2.9:

  • Fourniture d'expressions complètement explicites pour ψ+\psi^+ et ψ\psi^-
  • Estimées de bornes : γαi,βi1\gamma \leq \alpha_i, \beta_i \leq 1, où γ=(λλ24)/2\gamma = (\lambda-\sqrt{\lambda^2-4})/2

Exemple jouet (fin de la section 2.3): Pour l'environnement fonction échelon (wi=Mw_i=M lorsque i0i\leq 0, wi=0w_i=0 lorsque i1i\geq 1) :

undefined