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

Résumé technique : Temps de premier passage pour un oscillateur harmonique sous-amorti et application aux moteurs d'information

Énoncé du problème
L'estimation de la distribution du temps de premier passage (FPT - First Passage Time) est un problème fondamental en physique statistique, avec des applications allant des taux de réaction chimique aux moteurs d'information. Bien que les distributions FPT soient bien établies pour les systèmes suramortis décrits par des équations de Langevin, le régime sous-amorti reste complexe. Dans les systèmes sous-amortis, l'inertie introduit une dérivée seconde du déplacement dans les équations du mouvement stochastique, rendant le processus non markovien. Par conséquent, le calcul de la probabilité de FPT a historiquement été une « tâche colossale », des solutions analytiques n'étant généralement disponibles que pour des cas spécifiques (par exemple, la particule libre dans la limite de temps long). Cet article traite du calcul de la distribution du FPT, $P(t_{fp})$, pour un oscillateur harmonique brownien présentant un facteur de qualité $Q$ fini et franchissant un seuil de position $x_B$.

Méthodologie
Les auteurs proposent un cadre théorique unifié combinant trois approпредеches distinctes pour couvrir toute la gamme des échelles de temps et des facteurs de qualité :

1. Décomposition de la PDF : La fonction de densité de probabilité totale $P(t_{fp})$ est décomposée en trois contributions :

   • $P_I \delta(t_{fp})$ : Représente le déclenchement instantané où la position initiale $x(0) \geq x_B$.
   • $P_{II}(t_{fp})$ : Décrit la dynamique au cours de la première période d'oscillation ($0 < t_{fp} \leq 2\pi$).
   • $P_{III}(t_{fp})$ : Capture le comportement à long terme ($t_{fp} > 2\pi$).

2. Régime de court terme ($t_{fp} \leq 2\pi$) : Pour un faible amortissement ($Q \gg 1$), les auteurs utilisent une approximation hamiltonienne. Ils supposent que l'énergie totale $E$ est approximativement constante sur une période.

   • Déclenchement instantané ($P_I$) : Calculé à l'aide de la distribution de Boltzmann pour la position initiale, ce qui donne $P_I = \frac{1}{2}\text{erfc}(\sqrt{B})$, où $B = \frac{1}{2}x_B^2$ est la barrière d'énergie.
   • Première période ($P_{II}$) : La dynamique est traitée comme un mouvement sinusoïdal déterministe $x(t) = x_E \cos(\theta_0 + t)$. La probabilité est dérivée par intégration sur la distribution jointe de l'énergie initiale $E$ et de la phase $\theta_0$. Cela produit un plateau suivi d'une décroissance rapide, déterminée par l'énergie maximale $E^\dagger(t)$ accessible au temps $t$.

3. Régime de long terme ($t_{fp} > 2\pi$) : Pour des temps dépassant une période d'oscillation, le comportement du système est régi par la diffusion de l'énergie due au bruit thermique.

   • Limite de haut $Q$ : Le problème est mis en correspondance avec un problème d'échappement de Kramers à travers une barrière d'énergie. Le taux d'échappement $\Gamma(t)$ est calculé sur la base du seuil d'énergie dépendant du temps $E^\dagger(t)$, convergeant finalement vers un taux constant $\Gamma_B$.
   • Q fini : Les auteurs déterminent le taux d'échappement $\Gamma$ en calculant la plus petite valeur propre $-\lambda_1(x_B)$ de l'opérateur différentiel de Kramers associé à l'équation stochastique, soumis à une condition limite absorbante à $x_B$.
   • La queue résultante $P_{III}$ suit une décroissance exponentielle $e^{-\Gamma_B t_{fp}}$.

Résultats clés
L'article présente une expression analytique sous forme fermée pour la distribution FPT (Éq. 2), paramétrée uniquement par la barrière d'énergie $B$ et le facteur de qualité $Q$.

• Validation expérimentale : Les prédictions théoriques ont été testées par rapport à une expérience utilisant un micro-levier (micro-cantilever) sous-amorti ($Q \approx 7$) à température ambiante. La position a été mesurée avec une résolution de l'ordre du picomètre à une fréquence d'échantillonnage de 2 MHz.
• Accord : Les données expérimentales ont montré un « excellent accord » avec le modèle théorique à travers trois régimes distincts :
  1. La probabilité de déclenchements instantanés ($P_I$).
  2. La valeur du plateau de la distribution pour $0 < t < \pi$ ($P_{II}$).
  3. Le taux de décroissance exponentielle de la queue à long terme ($P_{III}$).
• Application aux moteurs d'information : La distribution FPT dérivée a été appliquée pour estimer la puissance d'un moteur d'information. Dans cette configuration, un travail est extrait lorsque l'oscillateur franchit un seuil $x_B$, déplaçant le puits de potentiel. La puissance moyenne $\langle P \rangle$ a été calculée comme le rapport entre le travail moyen extrait et la somme du FPT moyen et d'un temps de séjour. La courbe de puissance théorique correspondait aux données expérimentales, identifiant un point de fonctionnement optimal où le seuil $x_B \approx 2$ et le déplacement $L \approx 1$.

Signification et revendications
Les auteurs affirment que ce travail apporte une « nouvelle perspective » sur le problème du FPT pour les systèmes sous-amortis en offrant une solution couvrant toute la gamme des facteurs de qualité, comblant ainsi le fossé entre la dynamique hamiltonienne à court terme et l'échappement diffusif à long terme.

• Avancée théorique : L'article affirme que la solution analytique dérivée (Éq. 2) représente une avancée dans les statistiques FPT pour les systèmes sous-amortis. Il note que la seule solution connue auparavant pour un système sous-amorti était celle de la particule libre (processus d'accélération aléatoire) dans la limite de temps long.
• Distinction de la théorie de Kramers : Les auteurs soulignent que leur description diffère de la théorie de taux de Kramers standard, qui s'applique généralement au franchissement diffusif de barrières importantes. Dans leur distribution, la queue exponentielle n'émerge que pour des temps significativement supérieurs à une période d'oscillation, tandis que le comportement à court terme est dominé par les oscillations inertielles.
• Utilité pratique : La méthode permet d'estimer précisément la puissance dans les moteurs d'information possédant de l'inertie, un paramètre auparavant difficile à quantifier sans la connaissance complète de la distribution. Ils affirment que cela ouvre la voie à l'application des statistiques FPT à d'autres domaines impliquant une dynamique sous-amortie, tels que la thermodynamique stochastique.

L'article conclut en précisant que les détails complets des calculs et des simulations numériques concernant l'effet du facteur de qualité sont fournis dans un article accompagnant [22].