Polynomial slowdown in an angle-dependent 2d branching Brownian motion

Résumé technique : Ralentissement polynomial dans un mouvement brownien branchant bidimensionnel dépendant de l'angle

Énoncé du problème
L'article étudie le comportement asymptotique du déplacement maximal, $M_t$, dans un mouvement brownien branchant (MBB) spatialement inhomogène dans $\mathbb{R}^2$. Contrairement au MBB standard où les particules se branchent à un taux constant, ce modèle présente un taux de branchement $b(\theta)$ qui dépend uniquement de la coordonnée angulaire $\theta$ de la particule. La fonction $b$ est supposée être $2\pi$-périodique, atteignant son maximum unique en $\theta = 0$ (la direction privilégiée), et satisfaisant le développement local $b(\theta) = 1 - \beta|\theta|^\alpha + O(\theta^2)$ lorsque $\theta \to 0$, avec $\alpha \in (2/3, 2)$ et $\beta > 0$.

La question centrale est de savoir comment cette inhomogénéité angulaire affecte la vitesse et les fluctuations des particules de tête par rapport au cas homogène. Dans le MBB 1d standard, le maximum croît comme $\sqrt{2t} - \frac{3}{2\sqrt{2}}\log t$. Dans le MBB 2d homogène, la correction logarithmique devient $-\frac{1}{\sqrt{2}}\log t$. Les auteurs visent à déterminer le développement asymptotique précis de $M_t$ lorsque le taux de branchement est maximal uniquement dans une direction spécifique et décroît polynomialement avec l'angle.

Méthodologie
La stratégie de preuve combine des techniques probabilistes avec l'analyse d'équations aux dérivées partielles (EDP), en exploitant spécifiquement la formule de Feynman-Kac pour relier l'espérance du processus de branchement à la solution d'une équation de la chaleur avec un potentiel.

1. Lemmes Many-to-One et Many-to-Few : Les auteurs réduisent les espérances sur la population branchante à des espérances sur un (ou deux) mouvement(s) brownien(s) pondéré(s) par l'intégrale du taux de branchement le long du chemin.
   $$\mathbb{E}\left[\sum_{u \in N(t)} f(\text{chemin}_u)\right] = \mathbb{E}\left[f(\text{chemin}) \exp\left(\int_0^t b(\theta_s) ds\right)\right]$$
2. Analyse d'EDP via la théorie de Sturm-Liouville : Un composant critique est l'analyse du noyau $G(s, x; t, y)$ associé au mouvement brownien pondéré. En redimensionnant le temps et l'espace, le problème est mappé sur une EDP inhomogène en temps sur $[0, 1) \times \mathbb{R}$ :
   $$\partial_r u = \varrho \left( \partial_{xx} u - (1-r)^{-\alpha}|x|^\alpha u \right)$$
   Les auteurs utilisent la théorie de Sturm-Liouville pour analyser l'opérateur associé $L_q f = -f'' + q|x|^\alpha f$. Ils établissent des estimations précises pour la solution fondamentale de cette EDP, montrant que la décroissance est régie par la valeur propre principale $\lambda_0$ de l'opérateur $L = -f'' + |x|^\alpha f$.
3. Localisation des trajectoires : La preuve repose sur la localisation des trajectoires des particules extrémales. Il est démontré que les particules réalisant le déplacement maximal doivent rester dans un cône angulaire étroit autour de $\theta=0$. Plus précisément, la coordonnée $Y$ (transverse à la direction privilégiée) est contrainte à se mettre à l'échelle comme $t^{\kappa/2}$ où $\kappa = \frac{2\alpha}{2+\alpha}$. Cette localisation permet l'approximation $\theta_s \approx Y_s / (\sqrt{2}s)$, transformant la dépendance angulaire en une dépendance de la coordonnée transverse.
4. Bornes supérieures et inférieures :
   • Borne supérieure : Établie en utilisant un calcul de premier moment (lemme many-to-one) combiné à une « pseudo-martingale dérivée » et des arguments de barrière. Les auteurs comparent le processus à un MBB 1d standard avec une barrière courbe pour contrôler les corrections logarithmiques.
   • Borne inférieure : Établie en utilisant un calcul de second moment (lemme many-to-two) sur un ensemble spécifique de particules $\Gamma_{s_0}(t)$ satisfaisant des contraintes de trajectoire strictes (rester dans le cône angulaire et près du maximum prédit). Il est démontré que le second moment est borné par le carré du premier moment, assurant l'existence de telles particules avec une probabilité élevée.

Contributions et résultats clés

• Ralentissement polynomial : Le résultat principal (Théorème 1.1) est que le déplacement maximal $M_t$ présente un ralentissement polynomial par rapport à la vitesse standard $\sqrt{2t}$. Le développement asymptotique est donné par :
  $$M_t = \sqrt{2t} - \vartheta_1 \sqrt{2} t^{(2-\alpha)/(2+\alpha)} - \left( \frac{3}{2\sqrt{2}} - \frac{\alpha}{2\sqrt{2}(2+\alpha)} \right) \log t + O_P(1)$$
  Ici, $\vartheta_1$ est une constante explicite dépendant de $\beta$ et de la valeur propre principale $\lambda_0$ de l'opérateur $-f'' + |x|^\alpha f$.
• Optimalité des exposants : L'article prouve que l'exposant du terme polynomial, $\frac{2-\alpha}{2+\alpha}$, est optimal. De plus, la Proposition 1.2 établit que les particules à la distance $m(t)$ ont une coordonnée $Y$ d'ordre $t^{\alpha/(2+\alpha)}$ et un angle d'ordre $t^{-2/(2+\alpha)}$.
• Tightness : Les auteurs prouvent que la suite $(M_t - m(t))_{t \ge 1}$ est tendue (bornée en probabilité), bien qu'ils ne prouvent pas la convergence en distribution. Ils notent que la variable aléatoire limite serait probablement la limite d'une martingale approchée plutôt que la martingale dérivée trouvée dans les cas homogènes.
• Transition de phase : Les résultats sont spécifiques au régime $\alpha \in (2/3, 2)$. Les auteurs discutent du fait que pour $\alpha \ge 2$, la correction polynomiale devrait disparaître (ne laissant qu'une correction logarithmique), et pour $\alpha \le 2/3$, le comportement devient plus complexe avec potentiellement des termes polynomiaux supplémentaires.

Importance et affirmations
L'article prétend révéler une « variété plus riche de comportements » dans les processus de branchement en rendant le taux de branchement dépendant de l'angle plutôt que d'une coordonnée linéaire (ce qui ressemblerait à un modèle de bande). L'importance réside dans :

1. Correction polynomiale universelle : Le modèle fournit une nouvelle occurrence de ralentissement polynomial pour le déplacement maximal, distincte de la correction $T^{1/3}$ trouvée dans le MBB inhomogène en temps. Crucialement, ce modèle permet une plage continue d'exposants polynomiaux dans $(0, 1/2)$ en variant $\alpha$, whereas les modèles inhomogènes en temps produisent généralement un exposant fixe.
2. Interaction entre géométrie et branchement : Le travail démontre comment la géométrie de la direction privilégiée (l'angle) interagit avec la diffusion et le branchement pour créer un régime d'échelle spécifique pour la coordonnée transverse ($t^{\kappa/2}$).
3. Nouveauté technique : La dérivation d'estimations précises du noyau de chaleur pour une EDP inhomogène en temps avec un potentiel singulier $(1-t)^{-\alpha}|x|^\alpha$ en utilisant la théorie de Sturm-Liouville est une contribution technique significative, étendant les méthodes utilisées dans des travaux antérieurs sur le MBB inhomogène en temps (par exemple, Maillard et Zeitouni).

Les auteurs restent modestes concernant la convergence en distribution, reconnaissant que l'absence d'une structure de martingale exacte rend la preuve de la distribution limite une question ouverte. Ils laissent également les cas $\alpha \le 2/3$ et $\alpha \ge 2$ comme problèmes ouverts pour des recherches futures.