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

Resumo Técnico: Tempo de Primeira Passagem para um Oscilador Harmônico Subamortecido e Aplicação a Motores de Informação

Enunciado do Problema
A estimativa da distribuição do Tempo de Primeira Passagem (FPT, do inglês First Passage Time) é um problema fundamental na física estatística, com aplicações que variam desde taxas de reação química até motores de informação. Embora as distribuições de FPT sejam bem estabelecidas para sistemas superamortecidos descritos por equações de Langevin, o regime subamortecido permanece desafiador. Em sistemas subamortecidos, a inércia introduz uma segunda derivada do deslocamento nas equações de movimento estocásticas, tornando o processo não-markoviano. Consequentemente, calcular a probabilidade de FPT tem sido historicamente uma "tarefa assombrosa", com soluções em forma fechada geralmente disponíveis apenas para limites específicos (por exemplo, a partícula livre no limite de longo tempo). Este artigo aborda o cálculo da distribuição de FPT, $P(t_{fp})$, para um oscilador harmônico browniano com um fator de qualidade $Q$ finito cruzando um limiar de posição $x_B$.

Metodologia
Os autores propõem um arcabouço teórico unificado que combina três abordagens distintas para cobrir toda a gama de escalas de tempo e fatores de qualidade:

1. Decomposição da PDF: A função de densidade de probabilidade total $P(t_{fp})$ é decomposta em três contribuições:

   • $P_I \delta(t_{fp})$: Representa o disparo instantâneo onde a posição inicial $x(0) \geq x_B$.
   • $P_{II}(t_{fp})$: Descreve a dinâmica dentro do primeiro período de oscilação ($0 < t_{fp} \leq 2\pi$).
   • $P_{III}(t_{fp})$: Captura o comportamento de longo tempo ($t_{fp} > 2\pi$).

2. Regime de Curto Tempo ($t_{fp} \leq 2\pi$): Para baixo amortecimento ($Q \gg 1$), os autores empregam uma aproximação Hamiltoniana. Eles assumem que a energia total $E$ é aproximadamente constante ao longo de um período.

   • Disparo Instantâneo ($P_I$): Calculado usando a distribuição de Boltzmann para a posição inicial, resultando em $P_I = \frac{1}{2}\text{erfc}(\sqrt{B})$, onde $B = \frac{1}{2}x_B^2$ é a barreira de energia.
   • Primeiro Período ($P_{II}$): A dinâmica é tratada como um movimento sinusoidal determinístico $x(t) = x_E \cos(\theta_0 + t)$. A probabilidade é derivada integrando sobre a distribuição conjunta da energia inicial $E$ e da fase $\theta_0$. Isso resulta em um patamar seguido por um decaimento rápido, determinado pela energia máxima $E^\dagger(t)$ acessível no tempo $t$.

3. Regime de Longo Tempo ($t_{fp} > 2\pi$): Para tempos que excedem um período de oscilação, o comportamento do sistema é governado pela difusão de energia devido ao ruído térmico.

   • Limite de $Q$ Alto: O problema é mapeado para um problema de escape de Kramers através de uma barreira de energia. A taxa de escape $\Gamma(t)$ é computada com base no limiar de energia dependente do tempo $E^\dagger(t)$, convergindo eventualmente para uma taxa constante $\Gamma_B$.
   • Q Finito: Os autores determinam a taxa de escape $\Gamma$ calculando o autovalor mais lento $-\lambda_1(x_B)$ do operador diferencial de Kramers associado à equação estocástica, sujeito a uma condição de contorno absorvente em $x_B$.
   • A cauda resultante $P_{III}$ segue um decaimento exponencial $e^{-\Gamma_B t_{fp}}$.

Principais Resultados
O artigo apresenta uma expressão analítica em forma fechada para a distribuição de FPT (Eq. 2), parametrizada unicamente pelo limiar de energia $B$ e pelo fator de qualidade $Q$.

• Validação Experimental: As previsões teóricas foram testadas contra um experimento utilizando um micro-cantilever subamortecido ($Q \approx 7$) à temperatura ambiente. A posição foi medida com resolução de picômetros a uma taxa de amostragem de 2 MHz.
• Concordância: Os dados experimentais mostraram "excelente concordância" com o modelo teórico em três regimes distintos:
  1. A probabilidade de disparos instantâneos ($P_I$).
  2. O valor do patamar da distribuição para $0 < t < \pi$ ($P_{II}$).
  3. A taxa de decaimento exponencial da cauda de longo tempo ($P_{III}$).
• Aplicação em Motores de Informação: A distribuição de FPT derivada foi aplicada para estimar a potência de um motor de informação. Nesta configuração, o trabalho é extraído quando o oscilador cruza um limiar $x_B$, deslocando o poço de potencial. A potência média $\langle P \rangle$ foi calculada como a razão entre o trabalho médio extraído e a soma do FPT médio e um tempo de permanência (dwell time). A curva de potência teórica coincidiu com os dados experimentais, identificando um ponto de operação ideal onde o limiar $x_B \approx 2$ e o deslocamento $L \approx 1$.

Significância e Alegações
Os autores alegam que este trabalho fornece uma "nova visão" sobre o problema de FPT para sistemas subamortecidos ao oferecer uma solução que cobre toda a gama de fatores de qualidade, unindo a dinâmica Hamiltoniana de curto tempo e a fuga difusiva de longo tempo.

• Avanço Teórico: O artigo afirma que a solução em forma fechada derivada (Eq. 2) representa um avanço nas estatísticas de FPT para sistemas subamortecidos. Observa-se que a única solução em forma fechada conhecida anteriormente para um sistema subamortecido era para a partícula livre (processo de aceleração aleatória) no limite de longo tempo.
• Distinção da Teoria de Kramers: Os autores enfatizam que sua descrição difere da teoria de taxa de Kramers padrão, que tipicamente se aplica à travessia difusiva de grandes barreiras. Em sua distribuição, a cauda exponencial só emerge para tempos significativamente maiores que um período de oscilação, enquanto o comportamento de curto tempo é dominado por oscilações inerciais.
• Utilidade Prática: O método permite a estimativa precisa da potência em motores de informação com inércia, um parâmetro anteriormente difícil de quantificar sem o conhecimento da distribuição completa. Os autores afirmam que isso abre caminho para a aplicação das estatísticas de FPT a outros campos envolvendo dinâmica subamortecida, como a termodinâmica estocástica.

O artigo conclui observando que os detalhes completos dos cálculos e simulações numéricas referentes ao efeito do fator de qualidade são fornecidos em um artigo complementar [22].