An analytical solution of a quantum system with non-Markovian behavior: The Bixon-Jortner system in time domain

Resumo Técnico: Uma Solução Analítica do Sistema de Bixon-Jortner no Domínio do Tempo

Enunciado do Problema
O artigo aborda o desafio de obter uma descrição exata e analítica da dinâmica não-markoviana em sistemas quânticos. Embora o comportamento não-markoviano seja tipicamente estudado em sistemas bipartidos (sistema mais ambiente) onde o traço sobre o ambiente leva a efeitos de memória, as abordagens padrão frequentemente dependem de aproximações como a equação mestra de Lindblad (limite markoviano) ou correções perturbativas (ex: formalismo de Nakajima-Zwanzig). Os autores focam no modelo de Bixon-Jortner (BJ), um sistema consistindo em um único estado discreto $|\phi\rangle$ acoplado a uma "escada" infinita de estados $|n\rangle$. Diferente dos modelos padrão de sistema-ambiente, o espaço de Hilbert do BJ é uma soma direta em vez de um produto tensorial. Embora o modelo BJ seja conhecido por exibir decaimento exponencial (limite de Wigner-Weisskopf) e subsequente dinâmica de decaimento-recorrência, tratamentos analíticos anteriores baseavam-se em transformadas de Laplace que obscureciam a origem física dos "chutes" (transições não suaves) observados em intervalos de tempo inteiros. O objetivo é derivar a dinâmica exata no domínio do tempo e demonstrar explicitamente a natureza não-markoviana do sistema sem recorrer a aproximações.

Metodologia
Os autores atacam o problema da evolução temporal diretamente no domínio do tempo através das seguintes etapas:

1. Formulação: Partindo da equação de Schrödinger dependente do tempo na imagem de interação, eles derivam uma equação íntegro-diferencial para o coeficiente $b_\phi(t)$ do estado discreto. Esta equação elimina os coeficientes dos estados da escada em favor de uma integral sobre o histórico passado de $b_\phi$.
2. Transformação para uma Equação Diferencial com Atraso (DDE): Ao especializar para os parâmetros do modelo BJ (espaçamento uniforme da escada $\Delta$ e $E_\phi=0$) e aplicar a fórmula de somação de Poisson, a equação íntegro-diferencial é transformada em uma equação diferencial com atraso. Nesta forma, a derivada temporal de $b_\phi$ no tempo $\tau$ depende dos valores de $b_\phi$ em tempos anteriores $\tau - m$ (onde $m$ são inteiros), tornando explícitos os efeitos de memória não-markovianos.
3. Solução Analítica via Ansatz: Para resolver a DDE resultante, os autores empregam um ansatz de decomposição aditiva. Eles definem uma função transformada $f(\tau) = b_\phi(\tau)e^{\beta\tau}$ e assumem que a solução pode ser escrita como uma soma de peças aditivas $f_n(\tau-n)$ que "entram em ação" em cada tempo inteiro $\tau = n$.
4. Recorrência e Funções Geradoras: Ao analisar a equação dentro de intervalos de tempo específicos $(k, k+1)$, eles derivam uma relação de recorrência para os componentes aditivos. Esta relação é resolvida usando uma abordagem de função geradora, que vincula a solução aos polinômios de Laguerre $L_n(x)$.

Principais Resultados

• Equação Não-Markoviana Exata: O artigo deriva uma equação diferencial com atraso exata (Eq. 19) para a amplitude do estado discreto. A equação mostra explicitamente que a evolução do sistema é governada por um termo de decaimento exponencial markoviano mais uma soma de termos atrasados. Esses termos atrasados introduzem "chutes" a cada unidade de tempo inteiro $\tau = n$, causando desvios do decaimento exponencial simples e levando a recorrências parciais.
• Solução em Forma Fechada: Os autores obtêm uma solução analítica exata para a amplitude de probabilidade $b_\phi(\tau)$ (Eq. 56). A solução é expressa como uma soma envolvendo polinômios de Laguerre:
  $$b_\phi(\tau) = e^{-\beta\tau} \left( 1 + \sum_{n=1}^{\infty} H(\tau-n) e^{\beta n} [L_n(2\beta(\tau-n)) - L_{n-1}(2\beta(\tau-n))] \right)$$
  onde $H$ é a função degrau de Heaviside e $\beta$ é um parâmetro de acoplamento adimensional.
• Dinâmica de Decaimento e Recorrência: A solução confirma que, para $0 \le \tau < 1$, o sistema passa por um decaimento exponencial puro com uma taxa $\Gamma = 2\pi v^2 / \Delta$ (consistente com a regra de ouro de Fermi). Em $\tau = 1, 2, 3, \dots$, os termos atrasados são ativados, modificando a taxa de decaimento e fazendo com que a probabilidade exiba transições não suaves e recorrências parciais, conforme visualizado nos gráficos fornecidos para vários valores de $\beta$.

Significância e Alegações
O artigo alega que sua principal contribuição é fornecer um estudo de caso explícito, acessível e exato de não-markovianidade em um sistema quântico.

• Transparência: Diferente dos métodos anteriores de transformada de Laplace que produzem fórmulas analíticas mas obscurecem o mecanismo físico dos "chutes", a abordagem de equação diferencial com atraso no domínio do tempo torna a natureza não-markoviana do sistema completamente transparente. Os efeitos de memória são diretamente visíveis como dependências de valores passados em intervalos inteiros.
• Acessibilidade: As técnicas de solução (transformação de íntegro-diferencial para DDE, ansatz simples e funções geradoras) são apresentadas como intuitivas e acessíveis, oferecendo um exemplo pedagógico claro de como a dinâmica não-markoviana surge e pode ser resolvida analiticamente.
• Utilidade do Modelo: Os autores posicionam o sistema BJ como um modelo "prístino" valioso para entender fenômenos de decaimento e recorrência, removendo as complicações matemáticas de interações luz-matéria mais complexas enquanto retém a física essencial.

O artigo conclui modestamente, sugerindo que, embora a solução exata seja derivada para este sistema específico de soma direta, a abordagem poderia potencialmente ser generalizada para outros sistemas, como escadas com acoplamento não constante ou sistemas bipartidos tradicionais de produto tensorial, para comparação com formulações padrão de equações mestras.