Variational non-gaussian approach to interacting spin-boson models

Resumo Técnico: Abordagem Variacional Não-Gaussiana para Modelos de Interação Spin-Bóson

Enunciado do Problema
Modelos de interação spin-bóson, onde graus de liberdade quânticos discretos (spins) acoplam-se a um ou mais modos bosônicos (fótons ou fônons), são centrais para a ciência quântica contemporânea, incluindo QED de cavidade, QED de circuito e sistemas spin-fônon. Um desafio primário na simulação desses sistemas surge da natureza de dimensão infinita do espaço de Hilbert bosônico combinada com correlações de muitos corpos fortes no setor de spin. Métodos numéricos convencionais enfrentam um compromisso (trade-off): a truncagem da base de ocupação bosônica introduz viés, enquanto o relaxamento do setor de spin para descrições de campo médio falha em capturar o emaranhamento essencial spin-bóson. Além disso, abordagens diretas de redes de tensores (como DMRG) aplicadas ao sistema composto completo frequentemente sofrem com o crescimento rápido da dimensão de ligação devido à necessidade de representar explicitamente o campo bosônico junto aos spins.

Metodologia
Os autores propõem um arcabouço variacional híbrido que desacopla a otimização do setor bosônico do setor de spin. O estado fundamental é representado dentro de um manifold variacional não-gaussiano compacto definido pelo ansatz:
$$|\psi_{\text{NGS}}\rangle = U_\lambda (U_{\text{GS}} |0_b\rangle \otimes |\phi_s\rangle)$$
onde:

1. Camada Bosônica Gaussiana ($U_{\text{GS}}$): Uma transformação gaussiana compreendendo operadores de deslocamento ($\Delta_R$) e compressão (squeezing, $\xi$) atuando sobre o vácuo bosônico $|0_b\rangle$. Isso captura as flutuações bosônicas e deslocamentos de campo médio.
2. Transformação de Vestimenta (Dressing Transformation, $U_\lambda$): Um operador unitário que introduz correlações spin-bósons além do nível de campo médio. Ele atua como um deslocamento dependente do spin (transformação do tipo polaron) controlado por um parâmetro variacional $\lambda$. Esta camada mapeia o acoplamento linear spin-bóson em interações efetivas spin-spin.
3. Setor de Spin Efetivo ($|\phi_s\rangle$): O estado de spin não é um simples estado produto, mas o estado fundamental exato de um Hamiltoniano de spin efetivo ($H_{\text{eff}}$). Este Hamiltoniano é obtido através da média do Hamiltoniano completo sobre o estado gaussiano bosônico e a transformação de vestimenta.

O método emprega um ciclo de otimização autoconsistente:

• Para parâmetros bosônicos fixos ($\Delta_R, \xi, \lambda$), o Hamiltoniano de spin efetivo é resolvido usando DMRG (Grupo de Renormalização de Matriz de Densidade) para obter $|\phi_s\rangle$ e valores de expectativa de spin.
• Esses observáveis de spin são usados para atualizar os parâmetros variacionais bosônicos através da minimização do funcional de energia variacional total.
• O ciclo itera até a convergência da energia e de observáveis selecionados.

Esta abordagem evita a truncagem explícita do campo fotônico. Em vez disso, os graus de liberdade bosônicos são integrados variacionalmente, mapeando o problema em um modelo de spin autoconsistente.

Principais Contribuições e Resultados
O artigo compara este arcabouço NGS-DMRG com os modelos Dicke e Dicke-Ising, confrontando os resultados com simulações DMRG de referência (que utilizam truncagens explícitas de fótons).

1. Precisão e Melhoria Variacional:

   • A abordagem não-gaussiana produz consistentemente energias de estado fundamental menores do que o limite gaussiano (campo médio), confirmando a melhoria variacional.
   • Para o modelo Dicke, o método reproduz com precisão a transição de fase quântica (QPT) superradiante e o escalonamento de tamanho finito, coincidindo com os resultados do DMRG de referência com recursos computacionais significativamente reduzidos.
   • No modelo Dicke-Ising, o arcabouço mapeia com sucesso um diagrama de fases complexo apresentando:
     • Transições contínuas (segunda ordem) entre fases normal e superradiante.
     • Transições descontínuas (primeira ordem) e pontos tricríticos no regime ferromagnético.
     • Uma fase intermediária estreita de coexistência de ordem antiferromagnética e superradiância (AFM-SP) no regime antiferromagnético, uma característica prevista pela teoria de campo médio, mas difícil de resolver numericamente.
   • O método captura flutuações de tamanho finito e converge para o limite termodinâmico, onde o emaranhamento spin-fóton desaparece, sendo consistente com limites analíticos estabelecidos.

2. Eficiência Computacional:

   • Redução da Dimensão de Ligação: Ao tratar os bósons variacionalmente, o solver DMRG opera em um Hamiltoniano efetivo apenas de spin. Isso resulta em uma dimensão de ligação MPS ($\chi$) significamente menor do que o DMRG de referência, que deve lidar explicitamente da cadeia composta spin-bóson.
   • Tempo de Execução: Apesar do overhead do ciclo autoconsistente, a redução da complexidade da rede de tensores leva a uma diminuição substancial no tempo total de execução. O ciclo de otimização externo demonstra convergência exponencial.
   • A vantagem é mais pronunciada em regimes com fortes correlações e fronteiras de fase complexas (ex: a fase de coexistência), onde o DMRG de referência exige dimensões de ligação e cutoffs de fótons muito grandes.

Significância
O artigo afirma que este arcabouço híbrido oferece um caminho algorítmico distinto em comparação com métodos de redes de tensores de otimização de base local (LBO). Ao capturar os modos fotônicos através de parâmetros variacionais (deslocamento, compressão, vestimenta) em vez de reconfiguração iterativa de base, o método permite que o solver DMRG foque exclusivamente nas correlações de spin.

A significância reside em fornecer um método numérico robusto e não-perturbativo capaz de acessar a física de spin-bóson fortemente correlacionada sem eliminar ou truncar arbitrariamente o campo bosônico. Os autores demonstram que esta abordagem pode reproduzir com precisão diagramas de fases complexos e comportamentos críticos em modelos como a cadeia Dicke-Ising, alcançando alta fidelidade com custo computacional reduzido. Os autores fornecem uma implementação de código aberto disponível em: https://github.com/Jpedromend/NGS, que suporta acoplamentos spin-spin de alcance arbitrário e desordem, com potenciais extensões para sistemas de múltiplos modos e geometrias de maior dimensão via outros solvers (ex: PEPS, QMC).