Expected perimeter of the convex hull of planar Brownian motion stopped upon exiting the unit disk

Ce papier calcule la valeur exacte du périmètre attendu de l'enveloppe convexe d'un mouvement brownien plan arrêté à la sortie du disque unité, en réduisant le problème à l'évaluation d'une déviation horizontale maximale via la mesure harmonique, tout en établissant des bornes pour l'aire attendue dont le calcul exact reste plus complexe.

L'essentiel

🌊 Le Voyage d'une Particule Perdue

Imaginez une petite bille (une particule de mouvement brownien) lâchée au centre d'une grande piscine ronde (le disque unité). Cette bille ne suit pas une ligne droite ; elle bouge de manière totalement aléatoire, comme si elle était poussée par des milliers de micro-vents invisibles. Elle dérive, tourne, fait des zigzags, jusqu'à ce qu'elle touche le bord de la piscine. À ce moment précis, elle s'arrête.

Les mathématiciens Hugo Panzo et Stjepan Šebek se sont posé une question fascinante sur le chemin parcouru par cette bille : Quelle est la taille de la "bulle" imaginaire qui engloberait tout son trajet ?

🎈 La "Bulle" de l'Histoire (L'Enveloppe Convexe)

Pour visualiser cela, imaginez que vous prenez un élastique et que vous l'étirez autour de tous les points où la bille est passée. L'élastique va se tendre pour former la forme la plus petite possible qui contient tout le chemin. En mathématiques, on appelle cela l'enveloppe convexe.

Le papier cherche à répondre à deux questions principales sur cette "bulle" :

1. Quelle est la longueur de son périmètre ? (La taille de l'élastique).
2. Quelle est sa surface ? (La quantité d'eau à l'intérieur de la bulle).

📏 Le Périmètre : Une Surprise Facile à Calculer

Pour trouver la longueur du périmètre, les auteurs ont utilisé une astuce de génie basée sur la symétrie.

• L'analogie du compas : Imaginez que vous regardez le trajet de la bille sous tous les angles possibles (comme si vous tourniez autour de la piscine). À chaque angle, la bille s'est éloignée d'un certain maximum du centre.
• La découverte clé : Grâce à une propriété magique des mouvements aléatoires (l'invariance conforme), ils ont découvert que le périmètre total est simplement lié à la plus grande distance horizontale que la bille a jamais atteinte avant de toucher le bord.

C'est comme si, pour connaître la taille de tout le cercle, il suffisait de savoir jusqu'où la bille est allée vers la droite.

Le résultat exact :
Les auteurs ont trouvé une formule précise. La longueur moyenne attendue de ce périmètre est d'environ 3,21.
Pour vous donner une idée : le périmètre de la piscine elle-même (un cercle de rayon 1) est de $2\pi \approx 6,28$. La "bulle" formée par le trajet aléatoire est donc beaucoup plus petite que la piscine, ce qui est logique car la bille n'a pas eu le temps de tout explorer avant de sortir.

📐 La Surface : Le Mystère Non Résolu

C'est ici que l'histoire devient plus compliquée.

Calculer la surface (l'aire) de cette bulle est un défi bien plus grand que de calculer le périmètre.

• Pourquoi ? Pour connaître la surface, il ne suffit pas de savoir jusqu'où la bille est allée. Il faut aussi savoir quand elle y est arrivée et quelle était sa position verticale à ce moment précis. C'est comme essayer de dessiner la forme exacte d'une tache d'encre qui s'étale de manière chaotique : c'est beaucoup plus dur à prédire que juste mesurer la largeur de la tache.

Les auteurs n'ont pas réussi à trouver une formule simple et exacte pour la surface moyenne (comme ils l'ont fait pour le périmètre). C'est un problème mathématique très difficile qui résiste encore à une solution élégante.

Cependant, ils ne sont pas restés les bras croisés :

1. Ils ont créé des bornes (une estimation minimale et maximale) pour cette surface. Ils savent que la surface moyenne est comprise entre 0,47 et 1,14.
2. Ils ont utilisé des simulations informatiques (des millions de voyages virtuels de la bille) pour estimer la valeur réelle. Leurs calculs suggèrent que la surface moyenne est d'environ 0,66.

🧠 En Résumé

Ce papier est une belle histoire de deux parties :

1. Le succès : Ils ont résolu l'énigme du périmètre en trouvant une formule exacte et élégante, en reliant le mouvement chaotique à une simple mesure de distance maximale.
2. Le défi : Ils ont montré que l'énigme de la surface est beaucoup plus tenace. Même si nous ne pouvons pas encore écrire la formule exacte, nous savons désormais dans quelles limites elle se situe grâce à des méthodes ingénieuses et des simulations.

C'est un exemple magnifique de la façon dont les mathématiques peuvent transformer un problème de "désordre" (le mouvement aléatoire) en une structure prévisible, tout en nous rappelant qu'il reste encore des mystères à résoudre !

Résumé technique

1. Problématique et Contexte

L'article étudie les propriétés géométriques stochastiques du mouvement brownien plan standard $W = (X_t, Y_t)_{t \ge 0}$ démarrant à l'origine, arrêté au moment de sa première sortie du disque unité ouvert $D = \{x \in \mathbb{R}^2 : \|x\| < 1\}$. Ce temps d'arrêt est noté $\tau_D = \inf\{t \ge 0 : W_t \notin D\}$.

L'objet principal d'étude est l'enveloppe convexe $H_{\tau_D}$ de la trajectoire du mouvement brownien sur l'intervalle de temps $[0, \tau_D]$. Les auteurs visent à calculer l'espérance mathématique du périmètre de cette enveloppe convexe, notée $E[P_{\tau_D}]$.

Ce travail se distingue de la littérature existante qui traite principalement :

• Du mouvement brownien arrêté à un temps fixe déterministe (résultats connus pour $E[P_1]$).
• Du mouvement brownien réfléchi à l'intérieur d'un domaine (conditions aux limites de Neumann).
  Ici, les auteurs considèrent un mouvement « tué » (conditions de Dirichlet) avec un temps de vie aléatoire fini presque sûrement, représentant une excursion unique du centre vers la frontière.

2. Méthodologie

La démonstration repose sur une combinaison ingénieuse d'outils probabilistes et d'analyse complexe :

A. Réduction via la Formule de Cauchy

Les auteurs utilisent la formule de la surface de Cauchy pour relier le périmètre d'un ensemble convexe compact à sa fonction de support $h(\theta)$.
$$P_{\tau_D} = \int_0^{2\pi} h(\theta) \, d\theta$$
où $h(\theta) = \sup_{0 \le t \le \tau_D} \langle W_t, e_\theta \rangle$.
Grâce à l'invariance par rotation du mouvement brownien et du disque, la loi de $h(\theta)$ est identique à celle du maximum horizontal $M = \sup_{0 \le t \le \tau_D} X_t$. Cela permet de réduire le problème au calcul de l'espérance de $M$ :
$$E[P_{\tau_D}] = 2\pi E[M]$$

B. Transformation Conformes et Mesure Harmonique

Pour déterminer la loi de $M$, les auteurs introduisent un domaine géométrique appelé disque tronqué $D_a = D \cap \{z \in \mathbb{C} : \Re(z) < a\}$.
L'événement $\{M \ge a\}$ équivaut à l'événement où le mouvement brownien sort du domaine $D_a$ par la corde verticale $V_a$ (et non par l'arc circulaire).
La probabilité de cet événement est égale à la mesure harmonique de $V_a$ vue de l'origine dans le domaine $D_a$.

Pour calculer cette mesure harmonique, les auteurs construisent une application conforme $f_a$ qui transforme le disque tronqué $D_a$ en le demi-plan supérieur $\mathbb{H}$ :

1. Étape 1 : Une transformation de Möbius $l_a$ envoie $D_a$ sur un secteur angulaire (un « coin ») $W_a$.
2. Étape 2 : Une application puissance $p_a$ ouvre ce coin pour le transformer en demi-plan supérieur $\mathbb{H}$.

En utilisant l'invariance conforme du mouvement brownien, la probabilité recherchée est égale à la mesure harmonique de l'image de la corde $V_a$ (qui devient le demi-axe réel positif $[0, \infty)$) vue de l'image de l'origine $f_a(0)$ dans $\mathbb{H}$.

C. Noyau de Poisson

Une fois dans le demi-plan supérieur, la mesure harmonique est calculée explicitement à l'aide du noyau de Poisson. Cela permet d'obtenir une expression analytique exacte pour la fonction de répartition de $M$.

3. Résultats Principaux

A. Distribution et Espérance du Maximum Horizontal

Les auteurs établissent la fonction de répartition cumulative de $M$ :
$$P(M < a) = \frac{2 \arcsin a}{\pi - \arccos a}, \quad 0 < a < 1$$
L'espérance $E[M]$ est ensuite déduite par intégration :
$$E[M] = \int_0^{\pi/2} \left( 1 - \frac{2t}{\frac{\pi}{2} + t} \right) \cos t \, dt \approx 0.511655$$

B. Périmètre Attendu (Théorème Principal)

En combinant les résultats précédents, l'article fournit une expression exacte pour le périmètre attendu de l'enveloppe convexe :
$$E[P_{\tau_D}] = 2\pi E[M] \approx 3.214826$$
Cette valeur est obtenue de manière rigoureuse et exacte, contrairement aux approximations asymptotiques souvent utilisées dans la littérature pour les temps longs.

C. Aire de l'Enveloppe Convexe

Pour l'aire $A_{\tau_D}$, la situation est plus complexe. La formule de Blaschke relie l'aire à $E[h(0)^2]$ et $E[h'(0)^2]$.

• Le terme $E[h(0)^2]$ (lié au carré du maximum) est calculable.
• Le terme $E[h'(0)^2]$ (lié à la coordonnée $Y$ au moment où le maximum en $X$ est atteint) est difficile car le temps d'arrêt $\tau_D$ brise la symétrie temporelle nécessaire pour appliquer la loi de l'arcsinus classique.

Les auteurs ne parviennent pas à trouver une forme close pour $E[A_{\tau_D}]$, mais ils fournissent :

• Une borne supérieure : $E[A_{\tau_D}] \le \pi E[M^2] \approx 1.1397$.
• Une borne inférieure non triviale : En utilisant l'enveloppe étoilée (star hull) $H^\star_{\tau_D}$, ils prouvent que $E[A_{\tau_D}] \ge \pi - 8/3 \approx 0.4749$.
• Simulations Monte Carlo : Des simulations numériques donnent une estimation de l'aire autour de $0.6612$.

4. Signification et Contributions

• Résultat Exact : C'est l'un des premiers résultats donnant une expression exacte (sous forme d'intégrale définie) pour une fonctionnelle géométrique (périmètre) d'une enveloppe convexe de mouvement brownien arrêtée à un temps d'arrêt aléatoire (temps de sortie).
• Méthodologie Élégante : L'approche combinant la formule de Cauchy, les transformations conformes spécifiques (disque tronqué vers demi-plan) et la mesure harmonique offre un cadre puissant pour résoudre des problèmes similaires de géométrie stochastique.
• Complémentarité : Ce travail complète les études sur le mouvement brownien réfléchi (conditions de Neumann) en traitant le cas du mouvement tué (conditions de Dirichlet), offrant ainsi une vue plus complète du comportement des enveloppes convexes dans des géométries confinées.
• Difficulté de l'Aire : L'article met en lumière la difficulté fondamentale de calculer l'aire attendue dans ce contexte, soulignant la complexité introduite par le couplage entre les coordonnées et le temps d'arrêt aléatoire, contrairement au cas du périmètre qui se factorise grâce à l'invariance rotationnelle.

En résumé, cet article apporte une solution analytique précise au problème du périmètre attendu et pose des bornes rigoureuses pour l'aire, tout en illustrant la puissance des méthodes de l'analyse complexe appliquées à la théorie des probabilités.