First passage time distribution in underdamped harmonic oscillators

Technische Samenvatting: Eerste Passagetijdverdeling in Ondergedempte Harmonische Oscillatoren

Probleemstelling
Het artikel behandelt de berekening van de eerste passagetijdverdeling (FPTd), aangeduid als $P(t_{fp})$, voor de positie $x$ van een ondergedempte harmonische oscillator om een drempelwaarde $x_B$ te overschrijden. Hoewel de FPTd goed gevestigd is voor overdempte systemen, blijft het complex voor ondergedempte systemen (gekarakteriseerd door een hoge kwaliteitsfactor $Q$) vanwege de door traagheid veroorzaakte langetermijngeheugeneffecten. De vorm van de verdeling is sterk afhankelijk van $Q$, wat verschillende theoretische benaderingen vereist voor verschillende dempingsregimes. De auteurs streven ernaar een uitgebreide afleiding van $P(t_{fp})$ te bieden die het volledige bereik van $Q$ dekt, van intermediaire waarden ($Q \sim 10$) tot zeer grote waarden ($Q \gg 100$), en om de specifieke ruispatronen te analyseren die ontsnappingsevenementen aansturen.

Methodologie
De auteurs modelleren het systeem met behulp van de Langevin-vergelijking voor een harmonische oscillator met massa $m$, stijfheid $k$ en viskeuze demping $\gamma$, gekoppeld aan een thermisch bad bij temperatuur $T$. De dynamica wordt geanalyseerd in de faseruimte $(x, v)$ met behulp van dimensieloze eenheden. De afleiding van $P(t_{fp})$ wordt opgedeeld in drie afzonderlijke bijdragen op basis van tijdschalen en beginvoorwaarden:

1. Instantane bijdrage ($P_I$): Houdt rekening met trajecten die boven de drempelwaarde beginnen ($x(0) \ge x_B$).
2. Kortetermijnbijdrage ($P_{II}$): Houdt rekening met trajecten met een hoge initiële energie ($E(0) \ge B$) die de drempelwaarde binnen de eerste oscillatieperiode overschrijden.
3. Langetermijnbijdrage ($P_{III}$): Houdt rekening met trajecten die beginnen met een lage energie ($E(0) < B$) en die thermische activatie vereisen om de barrière te overschrijden.

Om het langetermijngedrag af te leiden, maken de auteurs gebruik van twee verschillende formalismen, afhankelijk van de kwaliteitsfactor:

• Voor intermediaire $Q$ ($Q \sim 10$): De auteurs maken gebruik van de Kramers lineaire differentiaaloperator met absorberende randvoorwaarden. Ze lossen de eigenwaarden en eigenfuncties van deze niet-zelfadjointe operator op. De langetermijn-FPTd wordt gedomineerd door de traagst vervallende eigenwaarde ($\lambda_1$), die overeenkomt met de ontsnappingssnelheid.
• Voor zeer grote $Q$ ($Q \gg 100$): Een model voor energiediffusie wordt afgeleid. Door te transformeren naar poolcoördinaten (energie $E$ en fase $\theta$) en te middelen over de snelle fasevariabele, wordt het systeem gereduceerd tot een Fokker-Planck-vergelijking voor energie alleen. De gemiddelde eerste passagetijd wordt berekend door de adjoint-vergelijking voor het energiediffusieproces op te lossen, gevolgd door een middeling over de initiële evenwichtsverdeling.

Daarnaast analyseren de auteurs de gemiddelde trajecten van positie, snelheid en thermische ruis die leiden tot $t_{fp}$ met behulp van instanton-theorie en numerieke simulaties om specifieke ruispatronen te identificeren.

Belangrijkste Resultaten

1. Decompositie van FPTd: De totale verdeling wordt uitgedrukt als de som $P(t_{fp}) = P_I \delta(t_{fp}) + P_{II}(t_{fp}) + P_{III}(t_{fp})$.
   • $P_I$ en $P_{II}$ worden afgeleid uit Hamiltoniaanse dynamica en zijn onafhankelijk van $Q$.
   • $P_{III}$ is afhankelijk van $Q$ en beschrijft de ontsnapping uit de potentiaalput.
2. Regime-afhankelijke ontsnappingssnelheden:
   • Bij intermediaire $Q$ wordt de ontsnappingssnelheid $\Gamma_B$ bepaald door de dominante eigenwaarde $\lambda_1$ van de Kramers-operator. Numerieke diagonalisatie van deze operator laat een uitstekende overeenstemming zien met directe Langevin-simulaties.
   • Bij zeer grote $Q$ levert de energiediffusie-benadering een gemiddelde eerste passagetijd $\tau_Z(B)$ en een bijbehorende snelheid $\Gamma_Z$ op. De auteurs stellen vast dat voor grote $Q$, de tijd om de drempelwaarde te bereiken na het overschrijden van de energiebarrière ongeveer $\pi$ (een halve periode) is, wat leidt tot een totale snelheid $\Gamma_Z = 1/(\pi + \tau_Z)$.
   • Het artikel merkt op dat voor intermediaire $Q$ de energiediffusie-benadering de gemiddelde FPT onderschat, omdat de oscillator de energiebarrière meerdere malen kan overschrijden voordat hij de ruimtelijke drempelwaarde $x_B$ bereikt.
3. Transient gedrag: Voor grote $Q$ vertoont de FPTd een "trapsgewijze" afname voordat deze het asymptotische exponentiële regime bereikt. Deze plateaus, met een periodiciteit van $2\pi$, ontstaan door de bijdrage van hogere orde eigenwaarden ($\lambda_{n>1}$) die imaginaire delen hebben die gerelateerd zijn aan de oscillatiefrequentie.
4. Ruispatronen en gemiddelde trajecten: De analyse van de gemiddelde trajecten die eindigen bij $t_{fp}$ onthult een specifieke vorm: een exponentiële groei van een sinusoïdaal signaal gevolgd door een exponentieel verval. De gemiddelde ruis $\langle \eta(t) \rangle$ vertoont een duidelijk patroon dat resonant is met de frequentie van de oscillator. De snelheidsverdeling op het moment van overschrijding wordt gevonden als $P(v) = v e^{-v^2/2}$, onafhankelijk van $Q$.

Betekenis en Claims
De auteurs beweren dat hun theoretische kader een volledige beschrijving biedt van de FPTd voor ondergedempte harmonische oscillatoren over alle kwaliteitsfactoren. De primaire betekenis ligt in:

• Verenigende benaderingen: Het aantonen hoe verschillende wiskundige instrumenten (Kramers-operator eigenwaarden versus energiediffusie) vereist zijn om het systeem accuraat te beschrijven, afhankelijk van het dempingsregime.
• Kwantitatieve overeenstemming: De theoretische voorspellingen vertonen een uitstekende overeenstemming met zowel directe numerieke simulaties van de Langevin-vergelijking als met experimentele gegevens die eerder zijn gemeten op een oscillator met $Q=7$ (verwezen in een begeleidende Letter [17]).
• Fysiek inzicht: De identificatie van specifieke ruispatronen die ontsnappingsgebeurtenissen aansturen, suggereert dat de gemiddelde FPT in ondergedempte systemen niet louter een statistisch gemiddelde is, maar wordt gedreven door een reproduceerbaar dynamisch pad (instanton) dat wordt gekenmerkt door een specifiek ruissignatuur.

Het artikel concludeert dat hoewel de exponentiële afname van de FPTd goed begrepen is, het transiënte gedrag (plateaus) en de specifieke ruispatronen die de ontsnapping aansturen nog gebieden zijn waar de analytische determinatie van hogere orde effecten uitdagend blijft, hoewel hun bestaan duidelijk wordt ondersteund door simulaties.