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

Questo studio determina un'espressione esatta per il perimetro atteso dell'involucro convesso di un moto browniano planare arrestato all'uscita dal disco unitario, riducendo il problema al calcolo dello spostamento orizzontale massimo e fornendo anche limiti non banali per l'area attesa, il cui calcolo esatto risulta invece molto più complesso.

L'essenziale

🌊 Il Viaggio Casuale e la sua "Ombra"

Immaginate di avere una goccia d'inchiostro (o un piccolo esploratore) che si trova esattamente al centro di un grande lago circolare perfetto (il "disco unitario").

Questo esploratore non sa dove andare: cammina in modo completamente casuale, come se fosse ubriaco o guidato dal vento. In matematica, questo movimento si chiama moto browniano.

L'esploratore cammina finché non tocca la riva del lago. Appena tocca la riva, il viaggio finisce.

Ora, immagina di tracciare con un pennarello verde tutto il percorso fatto dall'esploratore. Se prendessi un elastico e lo tirassi attorno a tutto quel percorso disordinato, l'elastico formerebbe una forma chiusa. Questa forma è chiamata inviluppo convesso (o "convex hull"). È come se l'elastico definisse il "contorno" minimo che racchiude tutto il viaggio.

📏 La Grande Domanda: Quanto è lungo il contorno?

Gli autori di questo studio (Hugo Panzo e Stjepan Šebek) si sono chiesti una cosa molto specifica:

"Se facciamo questo esperimento milioni di volte, qual è la lunghezza media di questo elastico verde?"

In termini matematici, vogliono calcolare il perimetro atteso di questa forma.

La Soluzione Geniale: Un Trucco di Geometria

Calcolare la lunghezza di una forma così irregolare e casuale sembra un incubo. Ma gli autori hanno usato un trucco intelligente, un po' come se avessero usato una lente magica.

1. Il Trucco della Rotazione: Hanno notato che, poiché il lago è perfettamente circolare e il movimento è casuale in tutte le direzioni, non importa da che parte guardiamo. La "massima distanza" che l'esploratore ha raggiunto verso destra è statisticamente uguale alla massima distanza verso l'alto, o verso qualsiasi altra direzione.
2. La Misura Semplice: Invece di misurare l'intero contorno complicato, hanno scoperto che basta guardare quanto l'esploratore è andato avanti (verso destra) prima di fermarsi. Se chiamiamo questa distanza massima $M$, il perimetro totale è semplicemente $2\pi$ volte questa distanza. È come dire: "Se so quanto è alto il tuo salto massimo, posso calcolare la lunghezza della tua ombra totale".
3. La Mappa Magica (Conformale): Per calcolare esattamente questa distanza massima $M$, hanno usato una tecnica di "piegatura dello spazio". Hanno trasformato il loro lago circolare (con un pezzo tagliato via) in un mezzo piano (come un foglio infinito tagliato a metà). In questo nuovo mondo "stirato", i calcoli diventano molto più facili, come risolvere un puzzle su un tavolo piatto invece che su una palla.

Il Risultato Numerico

Dopo tutti questi calcoli complessi, hanno trovato una formula precisa.
La lunghezza media dell'elastico verde è circa 3,21.
Per fare un confronto: il lago ha un raggio di 1, quindi la sua circonferenza è circa 6,28. L'elastico che racchiude il viaggio casuale è circa la metà della circonferenza del lago.

📐 E l'Area? (Il Mistero irrisolto)

Gli autori hanno anche provato a calcolare quanto spazio occupa l'area verde (l'area dentro l'elastico).
Qui la storia cambia. Mentre il perimetro è stato risolto con una formula elegante, l'area è molto più ostica.

Immagina di dover calcolare l'area di una nuvola di forma imprevedibile. È molto più difficile.

• Non sono riusciti a trovare una formula "pulita" e semplice per l'area esatta.
• Tuttavia, hanno creato dei limiti: hanno detto "l'area è sicuramente più grande di X e più piccola di Y".
• Hanno anche fatto delle simulazioni al computer (come se avessero fatto correre l'esploratore 100.000 volte virtualmente) e hanno scoperto che l'area media è circa 0,66.

🌟 Perché è importante?

Questo lavoro è importante perché:

1. Completa il quadro: Prima, si studiava spesso come si comporta un esploratore che rimbalza contro i muri (moto browniano riflesso). Qui, invece, studiano l'esploratore che viene "ucciso" (fermato) appena tocca il bordo. È una situazione diversa, più simile a un viaggio di sola andata.
2. Mostra la bellezza della matematica: Dimostra come concetti astratti (come le mappe conformi e le misure armoniche) possano essere usati per risolvere problemi pratici su forme casuali.
3. Metodi nuovi: Hanno mostrato che per alcune cose (come il perimetro) possiamo trovare risposte esatte, mentre per altre (come l'area) dobbiamo accontentarci di stime molto precise o limiti rigorosi.

In sintesi

Immaginate un esploratore che vaga a caso in un cerchio finché non tocca il bordo. Gli autori hanno scoperto che, nonostante il caos del suo cammino, la lunghezza media del suo "contorno" è un numero preciso e calcolabile (3,21), grazie a un geniale trucco geometrico. L'area, invece, rimane un mistero più profondo, per il quale abbiamo solo stime molto accurate.

Sommario tecnico

1. Problema e Contesto

Il lavoro si concentra sullo studio del convesso inviluppo (convex hull) generato da un moto browniano planare standard $W_t = (X_t, Y_t)$, iniziato dall'origine, che viene "ucciso" (fermato) al momento del primo uscita dal disco unitario aperto $D = \{x \in \mathbb{R}^2 : \|x\| < 1\}$.

L'obiettivo principale è calcolare il valore atteso del perimetro di questo inviluppo convesso, denotato come $E[P_{\tau_D}]$, dove $\tau_D$ è il tempo di uscita dal disco.
Il problema si inserisce in un contesto di ricerca esistente che ha già trattato:

• Il perimetro atteso per un moto browniano fissato a tempo deterministico (risolto da Takács nel 1980).
• Inviluppi convessi di moti browniani riflessi in geometrie confinate.
  Tuttavia, questo studio è originale nel considerare un moto browniano con condizioni al bordo di Dirichlet (assorbimento al bordo) e un tempo di arresto casuale $\tau_D$, rappresentando un'"escursione singola" dal centro al bordo.

2. Metodologia

Gli autori adottano un approccio che combina strumenti di geometria convessa, teoria della probabilità e analisi complessa (conformal invariance). La strategia si articola nei seguenti passaggi:

A. Riduzione tramite la Formula di Cauchy

Utilizzando la formula di Cauchy per l'area superficiale (Cauchy's surface area formula), il perimetro $P_{\tau_D}$ viene espresso come l'integrale della funzione di supporto $h(\theta)$ dell'inviluppo convesso lungo l'angolo $\theta \in [0, 2\pi)$:
$$P_{\tau_D} = \int_0^{2\pi} h(\theta) \, d\theta$$
Grazie all'invarianza rotazionale del moto browniano e del disco unitario, il valore atteso della funzione di supporto $E[h(\theta)]$ è costante e uguale al valore atteso del massimo della coordinata orizzontale $M = \sup_{0 \le t \le \tau_D} X_t$.
Ciò riduce il problema al calcolo di:
$$E[P_{\tau_D}] = 2\pi E[M]$$

B. Caratterizzazione della Distribuzione di $M$ tramite Misura Armonica

Per calcolare $E[M]$, gli autori determinano la funzione di distribuzione cumulativa (CDF) di $M$.

1. Dominio Troncato: Viene introdotto un dominio geometrico chiamato "disco troncato" $D_a = D \cap \{z : \Re(z) < a\}$. L'evento $\{M \ge a\}$ è equivalente all'evento in cui il moto browniano esce dal dominio $D_a$ attraverso la corda verticale $V_a$ (il segmento rettilineo che taglia il disco a $x=a$).
2. Mappatura Conforme: Viene costruita una mappatura conforme $f_a$ che trasforma il dominio $D_a$ nel semipiano superiore $\mathbb{H}$. Questa mappatura è composta da:
   • Una trasformazione lineare frazionaria $l_a$ che mappa $D_a$ in un settore angolare (wedge) $W_a$.
   • Una mappa di potenza $p_a$ che "apre" il settore angolare nel semipiano superiore $\mathbb{H}$.
3. Invarianza Conforme: Sfruttando l'invarianza conforme del moto browniano, la misura armonica di $V_a$ vista dall'origine in $D_a$ è uguale alla misura armonica dell'immagine di $V_a$ nel semipiano superiore, vista dall'immagine dell'origine.
4. Nucleo di Poisson: La probabilità $P(M \ge a)$ viene calcolata esplicitamente utilizzando il nucleo di Poisson del semipiano superiore, ottenendo una formula chiusa in termini di funzioni inverse trigonometriche.

C. Stima dell'Area

Per l'area dell'inviluppo convesso $E[A_{\tau_D}]$, gli autori notano che il problema è significativamente più complesso a causa della dipendenza dal tempo di massimizzazione casuale. Utilizzano la formula dell'area di Blaschke (o di Cauchy) e forniscono limiti rigorosi:

• Limite Superiore: Derivato direttamente dalla formula dell'area e dal calcolo del secondo momento di $M$.
• Limite Inferiore: Ottenuto confrontando l'inviluppo convesso con lo "star hull" (inviluppo stellato) del percorso, per il quale è possibile calcolare un'espressione esatta.

3. Risultati Chiave

Perimetro Atteso

Il risultato principale è l'espressione esatta per la funzione di distribuzione di $M$ e il conseguente valore atteso del perimetro.

• CDF di $M$: Per $0 < a < 1$:
  $$P(M < a) = \frac{2 \arcsin a}{\pi - \arccos a}$$
• Valore Atteso di $M$:
  $$E[M] = \int_0^{\pi/2} \left( 1 - \frac{2t}{\frac{\pi}{2} + t} \right) \cos t \, dt \approx 0.511655$$
• Perimetro Atteso:
  $$E[P_{\tau_D}] = 2\pi E[M] \approx 3.214826$$

Area Attesa

Non è stato possibile trovare una forma chiusa per il valore atteso dell'area $E[A_{\tau_D}]$. Tuttavia, sono stati stabiliti i seguenti limiti rigorosi:

• Limite Superiore: $E[A_{\tau_D}] \le \pi E[M^2] \approx 1.139699$.
• Limite Inferiore: Basato sull'area dello star hull, $E[A_{\tau_D}] \ge \pi - \frac{8}{3} \approx 0.474925$.
• Simulazioni Monte Carlo: Le simulazioni numeriche (su $10^5$ traiettorie) hanno fornito una stima dell'area di circa 0.6612, confermando che il valore reale si colloca all'interno dei limiti teorici.

4. Contributi e Significato

1. Soluzione Esatta per un Problema Aperto: Il paper fornisce la prima espressione esatta per il perimetro atteso dell'inviluppo convesso di un moto browniano planare arrestato al primo contatto con il bordo di un dominio (condizioni di Dirichlet), colmando un vuoto nella letteratura rispetto ai casi a tempo fisso o riflessi.
2. Metodologia Innovativa: L'uso combinato della formula di Cauchy, della mappatura conforme su domini non standard (disco troncato) e della misura armonica dimostra una potente sinergia tra geometria complessa e probabilità.
3. Complessità dell'Area: Il lavoro evidenzia chiaramente la difficoltà intrinseca nel calcolare l'area attesa rispetto al perimetro. Mentre il perimetro dipende solo dal massimo di una coordinata, l'area richiede la conoscenza congiunta della posizione e del tempo di massimizzazione, che perde le proprietà di simmetria e indipendenza tipiche dei tempi fissi.
4. Validazione Numerica: L'accordo tra il valore teorico calcolato ($\approx 3.2148$) e la stima Monte Carlo ($\approx 3.2136$) valida sia la derivazione matematica che l'implementazione numerica.

In sintesi, questo articolo stabilisce un nuovo risultato fondamentale nella teoria dei moti browniani confinati, offrendo strumenti analitici precisi per il perimetro e limiti rigorosi per l'area, aprendo la strada a futuri studi sulle proprietà geometriche di tali processi stocastici.