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

Resumen Técnico: Tiempo de Primer Paso para un Oscilador Armónico Subamortiguado y Aplicación a Motores de Información

Planteamiento del Problema
La estimación de la distribución del Tiempo de Primer Paso (FPT, por sus siglas en la inglés) es un problema fundamental en la física estadística con aplicaciones que van desde las tasas de reacción química hasta los motores de información. Aunque las distribuciones de FPT están bien establecidas para sistemas sobreamortiguados descritos por ecuaciones de Langevin, el régimen subamortiguado sigue siendo un desafío. En los sistemas subamortiguados, la inercia introduce una segunda derivada del desplazamiento en las ecuaciones de movimiento estocásticas, lo que convierte al proceso en no markoviano. En consecuencia, computar la probabilidad de FPT ha sido históricamente una "tarea desalentadora", con soluciones analíticas disponibles generalmente solo para límites específicos (por ejemplo, la partícula libre en el límite de tiempo largo). Este artículo aborda el cálculo de la distribución de FPT, $P(t_{fp})$, para un oscilador armónico browniano con un factor de calidad $Q$ finito que cruza un umbral de posición $x_B$.

Metodología
Los autores proponen un marco teórico unificado que combina tres enfoques distintos para cubrir todo el rango de escalas de tiempo y factores de calidad:

1. Descomposición de la PDF: La función de densidad de probabilidad total $P(t_{fp})$ se descompone en tres contribuciones:

   • $P_I \delta(t_{fp})$: Representa el disparo instantáneo donde la posición inicial $x(0) \geq x_B$.
   • $P_{II}(t_{fp})$: Describe la dinámica dentro del primer periodo de oscilación ($0 < t_{fp} \leq 2\pi$).
   • $P_{III}(t_{fp})$: Captura el comportamiento de largo tiempo ($t_{fp} > 2\pi$).

2. Régimen de Corto Tiempo ($t_{fp} \leq 2\pi$): Para un bajo amortiguamiento ($Q \gg 1$), los autores emplean una aproximación hamiltoniana. Asumen que la energía total $E$ es aproximadamente constante durante un periodo.

   • Disparo Instantáneo ($P_I$): Calculado mediante la distribución de Boltzmann para la posición inicial, resultando en $P_I = \frac{1}{2}\text{erfc}(\sqrt{B})$, donde $B = \frac{1}{2}x_B^2$ es la barrera de energía.
   • Primer Periodo ($P_{II}$): La dinámica se trata como un movimiento sinusoidal determinista $x(t) = x_E \cos(\theta_0 + t)$. La probabilidad se deriva integrando sobre la distribución conjunta de la energía inicial $E$ y la fase $\theta_0$. Esto produce una meseta seguida de un rápido decaimiento, determinado por la energía máxima $E^\dagger(t)$ accesible en el tiempo $t$.

3. Régimen de Largo Tiempo ($t_{fp} > 2\pi$): Para tiempos que exceden un periodo de oscilación, el comportamiento del sistema está gobernado por la difusión de energía debido al ruido térmico.

   • Límite de $Q$ Alto: El problema se mapea a un problema de escape de Kramers a través de una barrera de energía. La tasa de escape $\Gamma(t)$ se calcula basándose en el umbral de energía dependiente del tiempo $E^\dagger(t)$, convergiendo eventualmente a una tasa constante $\Gamma_B$.
   • Q Finito: Los autores determinan la tasa de escape $\Gamma$ calculando el autovalor más lento $-\lambda_1(x_B)$ del operador diferencial de Kramers asociado a la ecuación estocástica, sujeto a una condición de contorno absorbente en $x_B$.
   • La cola resultante $P_{III}$ sigue un decaimiento exponencial $e^{-\Gamma_B t_{fp}}$.

Resultos Clave
El artículo presenta una expresión analítica de forma cerrada para la distribución FPT (Ec. 2) parametrizada únicamente por la barrera de energía $B$ y el factor de calidad $Q$.

• Validación Experimental: Las predicciones teóricas fueron contrastadas con un experimento utilizando un micro-cantiléver subamortiguado ($Q \approx 7$) a temperatura ambiente. La posición se midió con resolución de picómetros a una tasa de muestreo de 2 MHz.
• Acuerdo: Los datos experimentales mostraron un "excelente acuerdo" con el modelo teórico a través de tres regímenes distintos:
  1. La probabilidad de disparos instantáneos ($P_I$).
  2. El valor de la meseta de la distribución para $0 < t < \pi$ ($P_{II}$).
  3. La tasa de decaimiento exponencial de la cola de largo tiempo ($P_{III}$).
• Aplicación a Motores de Información: La distribución FPT derivada se aplicó para estimar la potencia de un motor de información. En esta configuración, el trabajo se extrae cuando el oscilador cruza un umbral $x_B$, desplazando el pozo de potencial. La potencia media $\langle P \rangle$ se calculó como la relación entre el trabajo extraído medio y la suma del FPT medio y un tiempo de permanencia. La curva de potencia teórica coincidió con los datos experimentales, identificando un punto de operación óptimo donde el umbral $x_B \approx 2$ y el desplazamiento $L \approx 1$.

Significancia y Reivindicaciones
Los autores afirman que este trabajo proporciona una "nueva visión" al problema del FPT para sistemas subamortiguados al ofrecer una solución que cubre todo el rango de factores de calidad, tendiendo un puente entre la dinámica hamiltoniana de corto tiempo y la escape difusivo de largo tiempo.

• Avance Teórico: El artículo afirma que la solución analítica derivada (Ec. 2) representa un avance en las estadísticas de FPT para sistemas subamortiguados. Señalan que la única solución de forma cerrada conocida anteriormente para un sistema subamortiguado era para la partícula libre (proceso de aceleración aleatoria) en el límite de tiempo largo.
• Distinción de la Teoría de Kramers: Los autores enfatizan que su descripción difiere de la teoría de tasa de Kramers estándar, que típicamente se aplica al cruce difusivo de grandes barreras. En su distribución, la cola exponencial solo emerge para tiempos significativamente mayores que un periodo de oscilación, mientras que el comportamiento de corto tiempo está dominado por oscilaciones inerciales.
• Utilidad Práctica: El método permite la estimación precisa de la potencia en motores de información con inercia, un parámetro previamente difícil de cuantificar sin el conocimiento de la distribución completa. Los autores afirman que esto abre el camino para aplicar las estadísticas de FPT a otros campos que involucran dinámica subamortiguada, como la termodinámica estocástica.

El artículo concluye señalando que los detalles completos de los cálculos y las simulaciones numéricas respecto al efecto del factor de calidad se proporcionan en un artículo adjunto [22].