WGQED1D: a microscopic N-atom waveguide-QED toolkit for single-photon scattering and weak-measurement experiments

Resumo Técnico de WGQED1D

Definição do Problema
A propagação de luz fraca através de um conjunto atômico é um problema fundamental na óptica quântica, relevante para fenômenos como transparência induzida eletromagneticamente, propagação de luz lenta e tempos de residência negativos pós-selecionados. Embora a teoria de resposta linear macroscópica (por exemplo, absorção de Beer–Lambert e dispersão de Kramers–Kronig) seja bem estabelecida, estabelecer a ponte entre a dinâmica microscópica de $N$ átomos individuais e a resposta óptica integrada de uma nuvem requer um tratamento cuidadoso de entrada–saída. Os autores observam que implementações diretas de fórmulas padrão frequentemente resultam em desacordos quantitativos com o limite de Beer–Lambert por fatores de ordem unitária. Esses desacordos são atribuídos a escolhas sutis de convenção relativas à estrutura bidirecional do acoplamento de guia de onda e à calibração entre a profundidade óptica e a eficiência de acoplamento por átomo.

Metodologia
O artigo apresenta o WGQED1D, um kit de ferramentas Python leve e validado, projetado para a simulação microscópica de espalhamento de fóton único por $N$ átomos de dois níveis em um guia de onda bidirecional de uma dimensão. A metodologia baseia-se nos seguintes componentes principais:

1. Modelo Físico: O sistema é modelado no subespaço de excitação única, justificado por entradas de fóton único e baixa saturação. A dinâmica é governada por um Hamiltoniano efetivo não-Hermitiano padrão ($H_{\text{eff}}$) atuando sobre $N$ amplitudes complexas ($b_j$). Este Hamiltoniano inclui:
   • Interações dipolo-dipolo coletivas mediadas por troca de fótons retardados ($e^{ik_0|z_j-z_k|}$).
   • Perda diagonal para o espaço livre.
   • Um termo de condução representando o pacote de ondas gaussiano de entrada.
2. Relação Entrada–Saída: O kit de ferramentas reconstrói o campo de saída ($\phi_{\text{out}}$) a partir das amplitudes atômicas usando uma relação entrada–saída específica que contabiliza a acumulação de fase dos átomos em relação ao limite da nuvem.
3. Calibração Crítica: Uma contribuição metodológica central é a identificação e verificação de uma calibração específica entre a eficiência de acoplamento de guia de onda por átomo ($\beta = \Gamma_{1D}/\Gamma$) e a profundidade óptica alvo ($OD$). Os autores estabelecem que o acordo quantitativo com o limite de Beer–Lambert em espaço livre requer $\beta = OD/(2N)$. Este fator de dois surge da natureza bidirecional da guia de onda, onde os átomos emitem para frente e para trás com igual probabilidade, contrastando com idealizações unidirecionais ("quirais").
4. Implementação Numérica: O kit de ferramentas utiliza evolução determinística via um método explícito de Runge–Kutta 4(5). Opera em $O(N^2)$ de memória e operações por passo de tempo, utilizando álgebra linear densa. Inclui ferramentas de análise espectral para extrair o tempo de excitação atômica pós-selecionado ($\tau_T$) usando a fórmula de valor fraco de Wiseman–Steinberg.

Principais Contribuições

• Toolkit WGQED1D: Um módulo Python reprodutível e de código aberto (~400 linhas) que implementa o modelo microscópico, a reconstrução entrada–saída e a análise de valor fraco espectral. Depende apenas de bibliotecas científicas padrão (NumPy, SciPy, Matplotlib).
• Calibração Microscópica: A derivação explícita e a verificação numérica da relação $\beta = OD/(2N)$, que resolve desacordos quantitativos sistemáticos entre simulações microscópicas e previsões macroscópicas de Beer–Lambert.
• Infraestrutura de Validação: O artigo documenta três validações distintas:
  1. Acordo com a transmissão analítica de pulso médio de Beer–Lambert para profundidades ópticas $OD \in [0.25, 3]$ com erros $<0,1\%$.
  2. Verificação da fórmula de valor fraco de Wiseman–Steinberg para tempos de excitação pós-selecionados.
  3. Comparação direta com dados experimentais de Angulo et al. (2026) referentes a tempos de residência atômica negativos.

Resultados

• Acordo com Beer–Lambert: Simulações com $N=500$ átomos reproduzem a transmissão analítica de Beer–Lambert com precisão superior a 0,1% em toda a faixa de profundidades ópticas testadas, confirmando a validade da calibração $\beta = OD/(2N)$ no limite diluído.
• Extração de Valor Fraco: O kit de ferramentas extrai com sucesso tempos de excitação pós-selecionados ($\tau_T$) consistentes com a fórmula analítica de valor fraco. As discrepâncias estão dentro de 10% e são atribuídas à resolução espectral finita e correções coletivas de $N$ finito em altas profundidades ópticas.
• Comparação Experimental: Quando comparado aos dados de Angulo et al. sobre tempos de residência pós-selecionados negativos em nuvens de rubídio frias, o kit de ferramentas reproduz três de seis pontos de dados dentro de seus erros de $1\sigma$. Para os pontos restantes (principalmente em $OD=4$), a simulação segue a tendência qualitativa e, em alguns casos, alinha-se mais proximamente com as medições experimentais do que a fórmula analítica de resposta linear.
• Desempenho: O kit de ferramentas roda eficientemente em hardware padrão, produzindo resultados convergidos para conjuntos realistas ($N \sim 200–1000$) em segundos ou minutos.

Significância e Alegações
O artigo posiciona o WGQED1D como uma "plataforma numérica reprodutível" para estudar medição fraca, dispersão anômala, atraso de grupo e observáveis pós-selecionados na transmissão de fóton único através de nuvens atômicas opticamente espessas. Os autores enfatizam que o kit de ferramentas prioriza a confiabilidade e a verificabilidade numérica em vez do desempenho bruto, oferecendo uma correspondência direta linha por linha com as equações governantes.

A significância do trabalho reside em fornecer uma ferramenta validada que lida corretamente com as sutilezas do acoplamento de guia de onda bidirecional, permitindo assim que pesquisadores realizem simulações microscópicas que correspondam quantitativamente aos limites macroscópicos e observações experimentais recentes. Os autores afirmam explicitamente que, para usuários que exigem tamanhos de sistema maiores ($N \gtrsim 10^4$), a estrutura algorítmica pode ser estendida para álgebra linear esparsa ou backends de GPU, embora estes não estejam implementados atualmente. O kit de ferramentas é lançado sob a licença MIT para facilitar a investigação científica aberta.