Come fanno le lucciole a lampeggiare insieme

Nessun direttore, e funziona lo stesso

Livello Base — linguaggio semplice, senza matematica

Su certi fiumi del Sud-est asiatico, migliaia di lucciole sugli alberi lungo la riva lampeggiano nello stesso istante, una volta dopo l'altra, per ore. Interi alberi si accendono e si spengono insieme. Ai primi viaggiatori europei che lo riferirono non si credette per un secolo, perché le spiegazioni ovvie - un capo, un segnale, un battito di palpebre dell'osservatore - erano tutte sbagliate.

Non c'è nessun capo. Ogni lucciola segue una sola minuscola regola: ha il suo ritmo interno, e quando vede lampeggiare una vicina sposta un pochino avanti il proprio timer. È tutta qui la regola. Nessuno comanda, nessuno sa che cosa stia facendo il gruppo, e ogni insetto reagisce soltanto a quelli che riesce a vedere.

Da quella piccola regola l'ordine compare da solo. Falle partire tutte a caso e all'inizio è solo un caos di scintillii. Poi piccoli gruppi si allineano, quei gruppi tirano dentro i vicini, e nel giro di un minuto o due tutto il prato pulsa insieme. Non è stato organizzato niente: si è montato da sé.

La parte interessante è che non succede sempre. Se le lucciole si badano appena l'una con l'altra, o se i loro ritmi naturali sono troppo diversi, la spintarella non ce la fa e restano sparpagliate per sempre. C'è una soglia, ed è netta: sotto, caos; sopra, ordine improvviso. Muovi lentamente il cursore dell'accoppiamento e guarda l'istante in cui scatta.

Il modello di Kuramoto e una vera transizione di fase

Livello Studente — le equazioni principali

Riduci una lucciola all'essenziale e hai un oscillatore di fase: una sola variabile \(\theta\) che scorre costante da 0 a \(2\pi\) e lampeggia quando si riavvolge. Dai a ognuno la sua frequenza naturale \(\omega_i\) pescata da una distribuzione, e accoppia ognuno a tutti gli altri con un termine che tira la sua fase verso quella dei vicini. Ne esce il modello di Kuramoto, \(\dot\theta_i = \omega_i + \frac{K}{N}\sum_j \sin(\theta_j - \theta_i)\), una delle equazioni più studiate della scienza non lineare.

Il trucco che lo rende risolvibile è il parametro d'ordine. Definisci \(re^{i\psi}\) come la media di \(e^{i\theta_j}\) su tutti gli oscillatori: \(r\) misura quanto le fasi siano raggruppate, da 0 per sparpagliate a 1 per identiche, e \(\psi\) è la fase media. Il termine di accoppiamento allora collassa in \(Kr\sin(\psi - \theta_i)\): ogni oscillatore viene tirato verso la media del gruppo con una forza proporzionale a quanto il gruppo sia già sincronizzato. La sincronia si autorinforza, ed è esattamente per questo che la transizione è netta invece che graduale.

Sotto un accoppiamento critico \(K_c\) l'unico stato stabile ha \(r = 0\) e gli oscillatori vanno alla deriva indipendenti. Sopra, una frazione di loro si aggancia alla fase media mentre gli estremi continuano a derivare, e \(r\) cresce come \(\sqrt{K - K_c}\): la crescita a radice quadrata è la firma di una transizione di fase del secondo ordine, la stessa forma matematica di un magnete che si ordina raffreddandosi. Per una dispersione lorentziana di frequenze naturali di semilarghezza \(\gamma\), la soglia è esattamente \(K_c = 2\gamma\).

Le lucciole vere non sono proprio così. Interagiscono attraverso impulsi discreti invece che con un accoppiamento continuo, e per lo più non correggono il ciclo in corso ma reimpostano la fase del successivo. I modelli a impulsi accoppiati di tipo Mirollo-Strogatz lo catturano meglio, e danno un risultato più forte: per oscillatori identici con accoppiamento eccitatorio a impulsi, la sincronia si raggiunge da quasi ogni condizione iniziale. La versione della natura è più robusta del modello liscio, non meno.

Soluzioni esatte, topologia di rete e la sincronia come principio di progetto

Livello Esperto — profondità matematica completa

01L'ansatz di Ott-Antonsen

Per decenni il modello di Kuramoto è stato trattabile solo allo stato stazionario. Ott e Antonsen hanno mostrato nel 2008 che per una distribuzione lorentziana di frequenze la dinamica a dimensione infinita collassa esattamente su una varietà bidimensionale, riducendo l'intero sistema a un'equazione differenziale ordinaria per il parametro d'ordine complesso. Questo ha reso risolvibile il comportamento transitorio, non solo i punti fissi, e ha aperto all'analisi esatta stati chimera, bistabilità e sincronizzazione forzata. La riduzione non è un'approssimazione: è una varietà invariante che attrae quasi tutte le condizioni iniziali.

02Stati chimera

Nel 2002 Kuramoto e Battogtokh hanno trovato una cosa che si riteneva impossibile: in un anello di oscillatori identici e accoppiati in modo identico, uno stato stabile in cui una regione si aggancia in perfetta sincronia mentre un'altra deriva in modo incoerente. Niente distingue le due regioni: la simmetria si rompe spontaneamente. Le chimere sono state poi prodotte in oscillatori chimici, metronomi meccanici e schiere di laser, e si sostiene siano rilevanti per il sonno uniemisferico di delfini e uccelli, dove metà cervello dorme mentre l'altra metà resta sveglia.

03Topologia, e la funzione di stabilità principale

L'accoppiamento tutti-con-tutti è una comoda finzione. Su una rete, che la sincronia sia stabile è determinato dallo spettro laplaciano del grafo: la funzione di stabilità principale di Pecora e Carroll separa la dinamica del singolo oscillatore dalla topologia, quindi la sincronizzabilità si riduce al rapporto fra autovalori \(\lambda_N/\lambda_2\). Il ricablaggio a piccolo mondo lo migliora drasticamente, il che è una ragione plausibile per topologie del genere nel tessuto neurale. Il risultato controintuitivo in quest'area è che aggiungere collegamenti può distruggere la sincronia, l'analogo su rete del paradosso di Braess.

04Che cosa fanno davvero le lucciole, biologicamente

In Pteroptyx malaccae la sincronia si ottiene aggiustando il periodo invece che reimpostando la fase: l'insetto cambia la propria frequenza libera per adattarla allo stimolo, e può raggiungere sfasamento zero, cosa che il semplice reset non può fare. La funzione probabile è un argomento di rapporto segnale-rumore: un coro sincronizzato è visibile alle femmine molto più lontano attraverso la vegetazione di un singolo individuo, e al suo interno il lampo di un maschio non viene mascherato dai vicini. Nella specie americana Photinus carolinus lo schema è ancora diverso - raffiche di lampi separate da buio, con la sincronia che emerge dentro ogni raffica e lavori recenti che mostrano come sia mediata da un piccolo numero di vicini visibili invece che da una media globale.

05Dove salta fuori la stessa equazione

Il modello non è una curiosità sugli insetti. Le reti elettriche sono ben descritte da oscillatori di Kuramoto del secondo ordine, coi generatori come fasi e la perdita di sincronia come blackout; lo stesso parametro d'ordine misura la coerenza della rete. Le cellule pacemaker circadiane del nucleo soprachiasmatico si sincronizzano attraverso un termine di accoppiamento esattamente di questa forma. Le schiere di giunzioni Josephson sono matematicamente identiche. E la sincronia patologica nel nucleo subtalamico è il bersaglio della stimolazione cerebrale profonda nella malattia di Parkinson, dove i protocolli di stimolazione desincronizzante sono progettati direttamente da questa teoria.