First passage time distribution in underdamped harmonic oscillators

Sintesi Tecnica: Distribuzione del Tempo di Primo Passaggio in Oscillatori Armonici Sottosmorzati

Enunciato del Problema
Il documento affronta il calcolo della distribuzione del tempo di primo passaggio (FPTd), denotata come $P(t_{fp})$, per la posizione $x$ di un oscillatore armonico sottosmorzato che supera una soglia $x_B$. Sebbene la FPTd sia ben stabilita per i sistemi sovrasmorzati, essa rimane complessa per i sistemi sottosmorzati (caratterizzati da un alto fattore di qualità $Q$) a causa degli effetti di memoria a lungo termine indotti dall'inerzia. La forma della distribuzione dipende fortemente da $Q$, richiedendo approcci teorici differenti per diversi regimi di smorzamento. Gli autori mirano a fornire una derivazione completa di $P(t_{fp})$ coprendo l'intero intervallo di $Q$, da valori intermedi ($Q \sim 10$) a valori molto elevati ($Q \gg 100$), e ad analizzare i pattern di rumore specifici che guidano gli eventi di fuga.

Metodologia
Gli autori modellano il sistema utilizzando l'equazione di Langevin per un oscillatore armonico con massa $m$, rigidezza $k$ e smorzamento viscoso $\gamma$, accoppiato a un bagno termico alla temperatura $T$. La dinamica è analizzata nello spazio delle fasi $(x, v)$ utilizzando unità adimensionali. La derivazione di $P(t_{fp})$ è segmentata in tre contributi distinti basati su scale temporali e condizioni iniziali:

1. Contributo Istantaneo ($P_I$): Tiene conto delle traiettorie che partono sopra la soglia ($x(0) \ge x_B$).
2. Contributo a Breve Termine ($P_{II}$): Tiene conto delle traiettorie con alta energia iniziale ($E(0) \ge B$) che attraversano la soglia entro il primo periodo di oscillazione.
3. Contributo a Lungo Termine ($P_{III}$): Tiene conto delle traiettorie che partono con bassa energia ($E(0) < B$) e che richiedono l'attivazione termica per superare la barriera.

Per derivare il comportamento a lungo termine, gli autori impiegano due formalismi distinti a seconda del fattore di qualità:

• Per $Q$ Intermedio ($Q \sim 10$): Gli autori utilizzano l'operatore differenziale lineare di Kramers con condizioni al contorno assorbenti. Risolvono per gli autovalori e le autofunzioni di questo operatore non autoadiabatico. La FPTd a lungo termine è dominata dallo slow-decaying eigenvalue più lento ($\lambda_1$), che corrisponde al tasso di fuga.
• Per $Q$ Molto Elevato ($Q \gg 100$): Viene derivato un modello di diffusione dell'energia. Trasformando in coordinate polari (energia $E$ e fase $\theta$) e mediando sulla variabile di fase veloce, il sistema viene ridotto a un'equazione di Fokker-Planck per la sola energia. Il tempo medio di primo passaggio viene calcolato risolvendo l'equazione adiabatica per il processo di diffusione dell'energia, seguito da una media sulla distribuzione di equilibrio iniziale.

Inoltre, gli autori analizzano le traiettorie medie di posizione, velocità e rumore termico che precedono $t_{fp}$ utilizzando la teoria dell'instanton e simulazioni numeriche per identificare specifici pattern di rumore.

Risultati Chiave

1. Decomposizione della FPTd: La distribuzione totale è espressa come la somma $P(t_{fp}) = P_I \delta(t_{fp}) + P_{II}(t_{fp}) + P_{III}(t_{fp})$.
   • $P_I$ e $P_{II}$ sono derivate dalla dinamica hamiltoniana e sono indipendenti da $Q$.
   • $P_{III}$ dipende da $Q$ e descrive la fuga dal pozzo di potenziale.
2. Tassi di Fuga Dipendenti dal Regime:
   • In regime di $Q$ intermedio, il tasso di fuga $\Gamma_B$ è determinato dall'autovalore dominante $\lambda_1$ dell'operatore di Kramers. La diagonalizzazione numerica di questo operatore mostra un eccellente accordo con le simulazioni dirette di Langevin.
   • Per $Q$ molto elevato, l'approccio di diffusione dell'energia fornisce un tempo medio di primo passaggio $\tau_Z(B)$ e un corrispondente tasso $\Gamma_Z$. Gli autori riscontrano che, per $Q$ elevato, il tempo per raggiungere la soglia dopo aver attraversato la barriera di energia è approssimativamente $\pi$ (mezzo periodo), portando a un tasso totale $\Gamma_Z = 1/(\pi + \tau_Z)$.
   • Il documento nota che per $Q$ intermedio, l'approssimazione di diffusione dell'energia sottostima il tempo medio di FPT perché l'oscillatore può attraversare la barriera di energia più volte prima di raggiungere la soglia spaziale $x_B$.
3. Comportamento Transitorio: Per $Q$ elevato, la FPTd presenta un decadimento "a gradini" prima di raggiungere il regime esponenziale asintotico. Questi plateau, con una periodicità di $2\pi$, derivano dal contributo degli autovalori di ordine superiore ($\lambda_{n>1}$) che possiedono parti immaginarie legate alla frequenza di oscillazione.
4. Pattern di Rumore e Traiettorie Medie: L'analisi delle traiettorie medie che terminano in $t_{fp}$ rivela una forma specifica: una crescita esponenziale di un segnale sinusoidale seguita da un decadimento esponenziale. Il rumore medio $\langle \eta(t) \rangle$ mostra un pattern distinto in risonanza con la frequenza dell'oscillatore. Si scopre che la distribuzione di velocità al momento del superamento è $P(v) = v e^{-v^2/2}$, indipendente da $Q$.

Significato e Rivendicazioni
Gli autori affermano che il loro quadro teorico fornisce una descrizione completa della FPTd per oscillatori armonici sottosmorzati attraverso tutti i fattori di qualità. La portata principale risiede nel:

• Unificare gli Approcci: Dimostrare come siano necessari diversi strumenti matematici (autovalori dell'operatore di Kramers rispetto alla diffusione dell'energia) per descrivere accuratamente il sistema a seconda del regime di smorzamento.
• Accordo Quantitativo: Le previsioni teoriche mostrano un eccellente accordo sia con le simulazioni numeriche dirette dell'equazione di Langevin che con i dati sperimentali precedentemente misurati su un oscillatore con $Q=7$ (riferito in una lettera allegata [17]).
• Intuizione Fisica: L'identificazione di specifici pattern di rumore che guidano gli eventi di fuga suggerisce che la FPT media nei sistemi sottosmorzati non è solo una media statistica, ma è guidata da un percorso dinamico riproducibile (instanton) caratterizzato da una specifica firma del rumore.

Il documento conclude che, sebbene il decadimento esponenziale della FPTd sia ben compreso, il comportamento transitorio (plateau) e i pattern di rumore specifici che guidano la fuga rimangono aree in cui la determinazione analitica degli effetti di ordine superiore è impegnativa, sebbene la loro esistenza sia chiaramente supportata dalle simulazioni.