Full Configuration Interaction Quantum Monte Carlo for Accurate $\textit{Ab Initio}$ Nuclear Structure Calculations

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

Definição do Problema
O desenvolvimento de métodos ab initio de muitos corpos para núcleos fortemente interagentes é um desafio central na teoria nuclear moderna. Embora a teoria de campo efetivo quiral (EFT) forneça um arcabouço sistemático para derivar interações nucleares a partir da cromodinâmica quântica, resolver a equação de Schrödinger de muitos corpos para núcleos de massa média e pesada permanece computacionalmente proibitivo.

• Modelo de Casca Sem Núcleo (NCSM): Embora bem-sucedido para núcleos leves, o NCSM é limitado pelo crescimento exponencial da dimensão do espaço de modelo, que rapidamente excede a capacidade computacional dos supercomputadores atuais (atingindo $\sim 10^{12}$ dimensões). Consequentemente, o NCSM depende fortemente de extrapolações para sistemas maiores.
• Expansões de Muitos Corpos Truncadas: Métodos como a Teoria de Perturbação de Muitos Corpos (MBPT), o Grupo de Renormalização de Similaridade In-Medium (IMSRG) e a Teoria de Coupled-Cluster (CC) oferecem escalonamento polinomial, permitindo cálculos para núcleos mais pesados. No entanto, eles sofrem de erros de truncamento. Alcançar alta precisão exige levar essas expansões a ordens superiores, o que aumenta a complexidade matemática e o custo computacional. Além disso, fornecer estimativas confiáveis de incertezas de muitos corpos nesses esquemas truncados é difícil, particularmente para sistemas com forte agrupamento (clustering) ou correlações coletivas onde a convergência pode ser lenta.

Existe uma necessidade crítica de um método que retenha a precisão da Interação de Configuração Total (FCI) sem o armazenamento explícito do Hamiltoniano, servindo como um referencial (benchmark) para validar e guiar o desenvolvimento de métodos de expansão truncados.

Metodologia
Os autores introduzem o método de Interação de Configuração Total por Monte Carlo de Quantos (FCIQMC) ao domínio da estrutura nuclear ab initio. O FCIQMC é um método projetor estocástico que resolve a equação de Schrödinger no tempo imaginário:
$$-\frac{d}{d\tau} |\Psi(\tau)\rangle = (\hat{H} - E_0) |\Psi(\tau)\rangle$$
As principais características da implementação incluem:

• Amostragem Estocástica: A função de onda é representada por um conjunto de "caminhantes" (walkers — partículas com sinal) distribuídos através de determinantes de Slater em todo o espaço de configuração completo. O coeficiente de um determinante é determinado pela população líquida de sinais dos caminhantes naquele determinante.
• Etapas de Propagação: O conjunto de caminhantes evolui através de três etapas estocásticas por passo de tempo imaginário ($\Delta\tau$):
  1. Geração (Spawning): Amostra acoplamentos de Hamiltoniano fora da diagonal ($H_{ij}$) para criar novos caminhantes em determinantes conectados.
  2. Morte/Clonagem (Death/Cloning): Amostra o termo diagonal ($H_{ii}$) deslocado por um parâmetro $S$, ajustando a população de caminhantes no determinante atual.
  3. Aniquilação: Cancela caminhantes de sinais opostos residindo no mesmo determinante, o que é crucial para superar o problema do sinal de férmion sem exigir uma aproximação de nó fixo ou uma função de onda guia.
• Aproximação do Iniciador: Para controlar o problema do sinal em grandes espaços de modelo, os autores empregam a aproximação do iniciador combinada com uma correção de deslocamento adaptativo. Isso permite estimativas controláveis e aperfeiçoáveis do limite FCI conforme a população de caminhantes ($N_w$) aumenta.
• Configuração do Sistema: Os cálculos foram realizados para $^4$He, $^8$Be, $^{12}$C e $^{16}$O usando uma interação de dois núcleons otimizada (N2LOopt) dentro de uma base de oscilador harmônico esférico ($\hbar\omega = 20$ MeV). Nenhuma truncagem de configuração foi aplicada nos cálculos de FCIQMC.

Resultos Principais
O estudo apresenta energias de estado fundamental ($E$) e raios de carga ($R_{ch}$) para os núcleos selecionados com incertezas de muitos corpos de nível sub-percentual.

• Validação em Espaços Pequenos: Em um espaço de modelo pequeno ($e_{max}=2$), os resultados do FCIQMC coincidiram com os resultados exatos de FCI obtidos via diagonalização direta, validando a precisão da amostragem estocástica.
• Desempenho em Grandes Espaços de Modelo: Em um espaço de modelo grande ($e_{max}=10$), onde a dimensão do espaço de Hilbert para $^{16}$O atinge $\sim 10^{33}$ (tornando a diagonalização direta impossível), o FCIQMC forneceu resultados precisos.
  • Para $^4$He, os resultados do FCIQMC concordaram perfeitamente com as extrapolações quase exatas do NCSM.
  • Para núcleos mais pesados ($^8$Be, $^{12}$C, $^{16}$O), onde o NCSM falha em convergir, o FCIQMC entregou resultados precisos.
• Benchmarking de Métodos Truncados:
  • MBPT: Mostrou os maiores desvios em relação aos resultados do FCIQMC.
  • IMSRG: Geralmente razoável, mas exibiu desvios significativos para o $^8$Be, onde correções de ordem superior (IMSRG(3f2)+T) desviaram-se ainda mais dos valores exatos, indicando a ausência de correlações coletivas.
  • Coupled-Cluster (CC): Mostrou melhoria sistemática de CCSD para CCSDT-3. No entanto, para o raio de carga do $^{16}$O, a inclusão de triplos iterativos líderes levou a um desvio maior comparado ao FCIQMC, destacando a dificuldade de convergência em certos observáveis.
• Efeitos de Agrupamento (Clustering): Os resultados para o $^8$Be, que possui uma estrutura molecular/de agrupamento pronunciada, demonstraram que as expansões baseadas em referência (como o IMSRG de referência única) lutam para convergir devido à omissão de correlações coletivas de alta ordem. O FCIQMC, ao incluir correlações a todos os ranks de excitação dentro da base, manteve a precisão para esses estados agrupados.

Significância e Alegações
O artigo afirma que o FCIQMC estabelece um novo padrão para cálculos de estrutura nuclear ab initio:

1. Solucionador de Espaço Total Estocástico: É capaz de tratar sistemas além do alcance do NCSM convencional ao amostrar diretamente a função de onda exata no espaço de configuração completo.
2. Incertezas Quantificadas: Ao contrário dos métodos de expansão truncada onde a estimativa de incerteza é desafiadora, o FCIQMC fornece resultados com incertezas de muitos corpos quantificadas e sub-percentuais.
3. Referencial para a Teoria: Os resultados servem como benchmarks exatos para avaliar a precisão das interações quirais e a convergência de métodos de muitos corpos de escalonamento polinomial (MBPT, IMSRG, CC). O estudo demonstra explicitamente o papel crucial das correlações de ordem superior que frequentemente faltam em esquemas truncados.
4. Impacto Futuro: Os autores postulam que este método facilitará o desenvolvimento de métodos de muitos corpos de próxima geração com maior precisão e ajudará a resolver enigmas de longa data na teoria nuclear, como a evolução dos raios de carga na cadeia do cálcio e as correções de quebra de isospin em decaimentos Fermi superpermitidos, onde efeitos de muitos corpos de ordem superior são centrais.

O trabalho representa a primeira aplicação do FCIQMC a cálculos ab initio precisos de núcleos finitos em interações fortes, marcando um passo significativo para a compreensão das forças nucleares e simetrias fundamentais com precisão sem precedentes.