First passage time for an underdamped harmonic oscillator and application to the power of an information engine

Technische Samenvatting: Eerste Passage Tijd voor een Ondergedempde Harmonische Oscillator en Toepassing op Informatie-engines

Probleemstelling
De schatting van de verdeling van de Eerste Passage Tijd (FPT) is een fundamenteel probleem in de statistische fysica met toepassingen die variëren van chemische reactiesnelheden tot informatie-engines. Hoewel FPT-verdelingen goed gevestigd zijn voor overgedempte systemen die worden beschreven door Langevin-vergelijkingen, blijft het ondergedempde regime een uitdaging. In ondergedempte systemen introduceert traagheid een tweede afgeleide van de verplaatsing in de stochastische bewegingsvergelijkingen, waardoor het proces niet-Markoviaans wordt. Bijgevolg is het berekenen van de FPT-waarschijnlijkheid historisch gezien een "ontmoedigende taak" geweest, waarbij gesloten vorm oplossingen over het algemeen alleen beschikbaar zijn voor specifieke limieten (bijv. het vrije deeltje in het lange-tijd limiet). Dit artikel behandelt de berekening van de FPT-verdeling, $P(t_{fp})$, voor een Brownse harmonische oscillator met een eindige kwaliteitsfactor $Q$ die een positiedrempel $x_B$ overschrijdt.

Methodologie
De auteurs stellen een verenigd theoretisch kader voor dat drie verschillende benaderingen combineert om het volledige bereik van tijdschalen en kwaliteitsfactoren te dekken:

1. Decompositie van de PDF: De totale waarschijnlijkheidsdichtheidsfunctie $P(t_{fp})$ wordt gedecomposed in drie bijdragen:

   • $P_I \delta(t_{fp})$: Vertegenwoordigt instantane triggering waarbij de initiële positie $x(0) \geq x_B$.
   • $P_{II}(t_{fp})$: Beschrijft de dynamica binnen de eerste oscillatieperiode ($0 < t_{fp} \leq 2\pi$).
   • $P_{III}(t_{fp})$: Vangt het gedrag op de lange termijn op ($t_{fp} > 2\pi$).

2. Korte-tijd regime ($t_{fp} \leq 2\pi$): Voor lage demping ($Q \gg 1$) gebruiken de auteurs een Hamiltoniaanse benadering. Ze nemen aan dat de totale energie $E$ nagenoeg constant is over één periode.

   • Instantane Trigger ($P_I$): Berekend met behulp van de Boltzmann-verdeling voor de initiële positie, wat resulteert in $P_I = \frac{1}{2}\text{erfc}(\sqrt{B})$, waarbij $B = \frac{1}{2}x_B^2$ de energiebarrière is.
   • Eerste Periode ($P_{II}$): De dynamica wordt behandeld als deterministische sinusoïdale beweging $x(t) = x_E \cos(\theta_0 + t)$. De waarschijnlijkheid wordt afgeleid door te integreren over de gezamenlijke verdeling van de initiële energie $E$ en fase $\theta_0$. Dit levert een plateau op gevolgd door een snelle afname, bepaald door de maximale toegankelijke energie $E^\dagger(t)$ op tijdstip $t$.

3. Lange-tijd regime ($t_{fp} > 2\pi$): Voor tijden die groter zijn dan één oscillatieperiode, wordt het gedrag van het systeem beheerst door energiediffusie door thermische ruis.

   • Hoge $Q$ Limiet: Het probleem wordt gemapt naar een Kramers-ontsnappingsprobleem over een energiebarrière. De ontsnappingssnelheid $\Gamma(t)$ wordt berekend op basis van de tijdsafhankelijke energiesdrempel $E^\dagger(t)$, om uiteindelijk te convergeren naar een constante snelheid $\Gamma_B$.
   • Eindige $Q$: De auteurs bepalen de ontsnappingssnelheid $\Gamma$ door het langzaamste eigengetal $-\lambda_1(x_B)$ van de Kramers-differentiaaloperator te berekenen die geassocieerd is met de stochastische vergelijking, onderhevig aan een absorberende randvoorwaarde bij $x_B$.
   • De resulterende staart $P_{III}$ volgt een exponentiële verval $e^{-\Gamma_B t_{fp}}$.

Belangrijkste Resultaten
Het artikel presenteert een analytische gesloten vorm expressie voor de FPT-verdeling (Eq. 2) die uitsluitend geparametriseerd is door de energiebarrière $B$ en de kwaliteitsfactor $Q$.

• Experimentele Validatie: De theoretische voorspellingen werden getest tegen een experiment met een ondergedempde micro-cantilever ($Q \approx 7$) bij kamertemperatuur. De positie werd gemeten met picometer resolutie bij een bemonsteringsfrequentie van 2 MHz.
• Overeenkomst: De experimentele data vertoonden "uitstekende overeenstemming" met het theoretische model over drie verschillende regimes:
  1. De waarschijnlijkheid van instantane triggers ($P_I$).
  2. De plateau-waarde van de verdeling voor $0 < t < \pi$ ($P_{II}$).
  3. De exponentiële vervalsnelheid van de lange-tijd staart ($P_{III}$).
• Toepassing op Informatie-engines: De afgeleide FPT-verdeling werd toegepast om het vermogen van een informatie-engine te schatten. In deze opstelling wordt arbeid geëxtraheerd wanneer de oscillator een drempel $x_B$ overschrijdt, wat de potentiaalput verschuift. Het gemiddelde vermogen $\langle P \rangle$ werd berekend als de ratio van de gemiddelde geëxtraheerde arbeid tot de som van de gemiddelde FPT en een verblijftijd. De theoretische vermogenscurve kwam overeen met de experimentele data, waarbij een optimaal werkpunt werd geïdentificeerd waar de drempel $x_B \approx 2$ en de verschuiving $L \approx 1$.

Betekenis en Claims
De auteurs beweren dat dit werk "nieuw inzicht" biedt in het FPT-probleem voor ondergedempte systemen door een oplossing te bieden die het volledige bereik van kwaliteitsfactoren dekt, waarbij de kloof tussen korte-tijd Hamiltonianse dynamica en lange-tijd diffusieve ontsnapping wordt overbrugd.

• Theoretische Vooruitgang: Het artikel stelt dat de afgeleide gesloten vorm oplossing (Eq. 2) een vooruitgang vertegenwoordigt in de FPT-statistiek voor ondergedempte systemen. Er wordt opgemerkt dat voorheen de enige bekende gesloten vorm oplossing voor een onderdempt systeem die van een vrij deeltje (random acceleration process) was in het lange-tijd limiet.
• Onderscheid van Kramers Theorie: De auteurs benadrukken dat hun beschrijving verschilt van de standaard Kramers-theorie, die doorgaans van toepassing is op de diffusieve oversteek van grote barrières. In hun verdeling verschijnt de exponentiële staart pas voor tijden die aanzienlijk groter zijn dan één oscillatieperiode, terwijl het korte-tijd gedrag wordt gedomineerd door inertiële oscillaties.
• Praktisch Nut: De methode maakt de precieze schatting van het vermogen in informatie-engines met traagheid mogelijk, een parameter die voorheen moeilijk te kwantificeren was zonder volledige kennis van de verdeling. De auteurs stellen dat dit de weg vrijmaakt voor het toepassen van FPT-statistiek op andere gebieden die betrokken zijn bij ondergedempte dynamica, zoals de stochastische thermodynamica.

Het artikel concludeert door te vermelden dat de volledige details van de berekeningen en numerieke simulaties met betrekking tot het effect van de kwaliteitsfactor worden verstrekt in een begeleidend artikel [22].