Full configuration interaction quantum Monte Carlo for accurate $\textit{ab initio}$ nuclear structure calculations: algorithms and calculation details

Resumo Técnico: Interação de Configuração Total por Monte Carlo Quântico para Cálculos Estruturais Nucleares Ab Initio Precisos

Definição do Problema
A teoria nuclear ab initio visa descrever núcleos finitos a partir de Hamiltonianos microscópicos derivados da teoria de campo efetivo quiral. Embora os métodos de expansão de muitos corpos (ex: Coupled Cluster, IMSRG, MBPT) estendam os cálculos para núcleos mais pesados, eles dependem de truncamentos que introduzem incertezas de muitos corpos residuais difíceis de quantificar. Por outro lado, o No-Core Shell Model (NCSM) fornece acesso direto à solução de interação de configuração total (FCI), mas é limitado pelo crescimento exponencial do espaço de Hilbert. Existe uma necessidade crítica de solvers estocásticos capazes de se aproximar do limite FCI em grandes espaços de modelo sem armazenar explicitamente a matriz do Hamiltoniano de muitos corpos, particularmente para sistemas com fortes correlações de ordem superior decorrentes de deformação, agrupamento (clustering) ou estruturas de camadas abertas.

Metodologia
Este trabalho adapta o algoritmo de Interação de Configuração Total por Monte Carlo Quântico (FCIQMC), originalmente desenvolvido para a química quântica, para cálculos de estrutura nuclear. O método emprega uma abordagem de projetor estocástico no tempo imaginário para resolver a equação de Schrödinger dentro de um setor de simetria escolhido (número de partículas, paridade e projeção de momento angular).

Componentes algorítmicos principais descritos incluem:

• Propagação no Tempo Imaginário: A função de onda é representada por uma população de walkers (caminhadores) com sinal, distribuídos através de determinantes de Slater. A propagação envolve etapas de geração (spawning), morte/clonagem e aniquilação para estabelecer estocasticamente a estrutura de sinal fermiônica.
• Geração de Excitações: Para melhorar a eficiência em relação à geração uniforme, os autores utilizam um gerador de heat-bath pré-computado que amostra preferencialmente excitações com grandes elementos de matriz do Hamiltoniano.
• Aproximação Iniciadora (i-FCIQMC): Para superar o problema de sinal em grandes espaços de Hilbert, a geração a partir de determinantes fracamente ocupados é restrita. Isso introduz um viés de walker finito que é sistematicamente reduzido conforme a população de walkers aumenta.
• Correção de Deslocamento Adaptativo (AS): Um deslocamento local é atribuído aos determinantes não-iniciadores para acelerar a convergência em direção ao limite FCI, preservando o limite exato conforme o número de walkers tende ao infinito.
• Estimadores de Observáveis:
  • Estimadores Projetados/de Teste: Utilizados para energias, utilizando determinantes de referência ou funções de onda de teste compactas.
  • Estimadores Puros: Essenciais para observáveis não comutativos, como raios. Estes são computados via amostragem da Matriz de Densidade Reduzida (RDM) usando duas réplicas independentes de walkers para eliminar o viés inerente aos estimadores quadráticos.
• Estados Excitados: O método é estendido para estados excitados via propagação multi-estado ortogonalizada, onde múltiplos conjuntos de walkers são propagados em paralelo e ortogonalizados em cada etapa.

Principais Contribuições

1. Adaptação do Formalismo: O artigo detalha a implementação específica de FCIQMC para Hamiltonianos nucleares, incluindo o tratamento de forças de três nucleons e os esquemas específicos de geração de excitação necessários para interações nucleares.
2. Validação contra FCI Determinístico: Em pequenos espaços de modelo ($e_{\text{max}}=2$) para $^4$He, $^8$Be, $^{12}$C e $^{16}$O, os resultados estocásticos de FCIQMC para energias e raios de protones pontuais ao quadrado ($R_p^2$) mostram reproduzir os valores exatos de FCI determinístico dentro das incertezas estatísticas.
3. Controle de Erro Sistemático: Os autores demonstram que o viés residual de walker finito nas aproximações iniciadoras pode ser controlado e quantificado. Eles introduzem uma estratégia usando extrapolações de lei de potência ($O(N_w) = O_\infty + \lambda N_w^{-\alpha}$) em grandes espaços de modelo ($e_{\text{max}}=10$). O expoente $\alpha$ é restringido pela exigência de consistência entre a energia projetada de convergência rápida e a energia pura de convergência mais lenta.
4. Extensão para Estados Excitados: O artigo demonstra a capacidade de calcular espectros de baixas energias (especificamente para $^6$Li) resolvendo múltiplos estados dentro do mesmo setor de simetria sem diagonalização explícita.

Resultados

• Benchmarking: No espaço $e_{\text{max}}=2$, as energias e raios de FCIQMC coincidem com os benchmarks de FCI. A correção de deslocamento adaptativo reduz significativamente o viés de walker finito e enfraquece a dependência do limiar do iniciador ($n_{\text{init}}$).
• Grandes Espaços de Modelo ($e_{\text{max}}=10$):
  • Para $^4$He, o cálculo converge diretamente com o aumento do número de walkers ($N_w$), coincidindo com os resultados de NCSM.
  • Para $^8$Be, $^{12}$C e $^{16}$O, a energia projetada converge rapidamente, enquanto os estimadores puros baseados em RDM (energia e raios) exibem convergência mais lenta. Ao aplicar a extrapolação de lei de potência restrita pela energia, os autores obtêm resultados que capturam correlações além dos métodos de expansão de muitos corpos truncados (MBPT, CC, IMSRG).
  • A estratégia de extrapolação é validada ao mostrar que os resultados são independentes do parâmetro de deslocamento do ajuste adaptativo ($\tilde{\Delta}$) e consistentes entre diferentes limiares de iniciador.
• Estados Excitados: Para $^6$Li, o método resolve com sucesso os cinco estados mais baixos ($1^+, 3^+, 0^+, 2^+_1, 2^+_2$). O espectro de excitação é corretamente ordenado mesmo em números menores de walkers, convergindo mais rápido que os resultados de NCSM em relação à truncagem de muitos corpos ($N_{\text{max}}$).

Significância
O artigo posiciona o FCIQMC como um solver estocástico sistematicamente aprimorável capaz de se aproximar do limite FCI para sistemas nucleares onde métodos determinísticos falham devido ao tamanho do espaço de Hilbert. Ele serve a dois papéis primários:

1. Solver Preciso: Fornece benchmarks de alta precisão para energias e raios de estado fundamental em grandes espaços de modelo, capturando altas ordens de correlação que métodos de expansão truncados perdem.
2. Ferramenta de Benchmark: Oferece um meio para avaliar as incertezas residuais de muitos corpos de outros métodos de expansão de muitos corpos truncados ao fornecer uma referência próxima à solução exata.

O trabalho estabelece a fundação algorítmica e numérica para os resultados físicos relatados em uma carta complementar [80], confirmando que o FCIQMC pode tratar núcleos de camadas fechadas e abertas dentro do mesmo framework e se estender à espectroscopia de estados excitados. Os autores observam que trabalhos futuros envolverão interações explícitas de três nucleons, bases de partícula única maiores e aplicações a observáveis coletivos e eletrofracos.