Full Uncertainty Quantification of Sign-Problem-Free Quantum Monte Carlo Methods and Nuclear Lattice Effective Field Theory Benchmarks

Resumen Técnico: Cuantificación Completa de la Incertidumbre de los Métodos de Monte Carlo de Cuantización de Signo Libre y Referencias de la Teoría de Campo Efectivo de Red Nuclear

Planteamiento del Problema
La teoría nuclear ab initio moderna ha logrado un éxito significativo en la reproducción de observables para núcleos ligeros y materia nuclear infinita combinando la teoría de campo efectivo (EFT) quiral con resolvedores de muchos cuerpos. Sin embargo, la extensión de estos métodos a núcleos de masa media y pesados (hasta $^{100}$Sn) enfrenta un cuello de botella crítico: la falta de una cuantificación rigurosa de la incertidumbre para los errores algorítmicos sistemáticos. Sin tal control, las discrepancias entre la teoría y el experimento suelen atribuirse a las fuerzas nucleares subyacentes, lo que conduce al "reajuste" de las constantes de baja energía (LEC) para compensar las deficiencias del resolvedor. Esto enmascara el verdadero origen de los errores teóricos y compromete el poder predictivo de los enfoques ab initio.

Específicamente, críticas recientes (Ref. [1]) han cuestionado la validez de las referencias de la Teoría de Campo Efectivo de Red Nuclear (NLEFT), citando problemas con los límites variacionales, las energías de correlación y los mecanismos de saturación nuclear. Estas críticas a menudo derivan de la comparación de resultados obtenidos de diferentes teorías de red reguladas (por ejemplo, mezclando formalismos de matriz de transferencia y de Hamiltoniano) o tratamientos inconsistentes de volumen finito, en lugar de fallos genuinos de los métodos de muchos cuerpos.

Metodología
Los autores emplean métodos de Monte Carlo de Cuántica (QMC) libres de problema de signo, específicamente el Monte Carlo de Campo Auxiliar (AFQMC), dentro del marco de NLEFT. El estudio se centra en la acción LAT-OPT1 recientemente desarrollada, un Hamiltoniano de red mínimo libre de problema de signo que incorpora una interacción de acoplamiento espín-órbita no perturbativa mientras preserva la simetría antiunitaria requerida para pesos de Monte Carlo positivos en núcleos pares-pares.

La metodología implica un análisis exhaustivo del presupuesto de error y un reexamen diagnóstico de las comparaciones de referencia:

1. Cuantificación de la Incertidumbre: Los autores cuantifican sistemáticamente las incertidumbres estadísticas, de tiempo imaginario, de volumen finito, de discretización temporal, inducidas por el campo auxiliar y de signo residual para núcleos que van desde $^4$He hasta $^{100}$Sn.
2. Tratamiento No Perturbativo del Acoplamiento Espín-Órbita: A diferencia de enfoques anteriores donde las fuerzas de espín-órbita se trataban como inserciones perturbativas, LAT-OPT1 integra la densidad de espín-órbita directamente en el Hamiltoniano de proyección en el tiempo euclidiano.
3. Verificaciones de Consistencia del Regulador: El artículo compara rigurosamente diferentes formalismos computacionales:
   • Matriz de Transferencia vs. Hamiltoniano: Analizando las diferencias entre los cálculos de la matriz de transferencia con espaciado de red temporal finito ($a_t$) y el límite del Hamiltoniano cuando $a_t \to 0$.
   • Condiciones de Contorno: Investigando el impacto de las Condiciones de Contorno Periódicas (PBC) frente a las Condiciones de Contorno de Borde Promedio (ATBC) en los efectos de capa de volumen finito en la materia nuclear.
   • Renormalización: Examinando el límite continuo y la renormalización de las interacciones de corto alcance.

Contribuciones Clave y Resultados

• Control Cuantitativo de Errores Sistemáticos:
  El estudio demuestra que la incertidumbre computacional acumulada de muchos cuerpos en las energías del estado fundamental para núcleos doblemente mágicos hasta $^{100}$Sn está muy por debajo del nivel porcentual.

  • Tiempo Imaginario y Volumen Finito: Las extrapolaciones al tiempo de proyección infinito ($\tau \to \infty$) y al volumen infinito ($L \to \infty$) son estables. Se demuestra que un tamaño de caja de $L=11$ (con $a=1.32$ fm) es suficiente para la convergencia de los núcleos estudiados.
  • Discretización Temporal: El uso de un formalismo de matriz de transferencia con $a_t$ finito induce un Hamiltoniano efectivo ($H_{eff}$) que difiere del Hamiltoniano objetivo ($H$). Los autores muestran que, aunque $H_{eff}$ depende de $a_t$, el error de energía escala como $O(a_t^2)$, permitiendo un control preciso.
  • Términos Inducidos por el Campo Auxiliar: El completado de cuadrados requerido para mantener el muestreo libre de problema de signo para el término de espín-órbita genera términos inducidos de dos cuerpos ($\rho_s^2$) y tres cuerpos. Estos se cuantifican y se encuentran como correcciones pequeñas (pocos MeV), confirmando la validez de la construcción.
  • Problema de Signo Residual: Para núcleos de $A$ impar y de par impar, persiste un leve problema de signo residual. Sin embargo, la fase promedio permanece cercana a la unidad (0.94–0.98), y las incertidumbres estadísticas siguen siendo comparables a las de los sistemas pares-pares.

• Rol No Perturbativo del Acoplamiento Espín-Órbita:
  Un hallazgo crítico es que la interacción de espín-órbita debe incluirse de forma no perturbativa en la proyección del tiempo euclidiano. Tratarla solo como una inserción perturbativa sobre una función de onda con simetría SU(4) omite contribuciones de estructura de capa de varios a diez MeV en núcleos de masa media. La inclusión de esta en la propagación genera dinámicamente la estructura de capa correcta.

• Resolución de Discrepancias de Referencia:
  El artículo aborda las críticas recientes (Ref. [1]) identificando errores conceptuales y técnicos en la comparación de diferentes teorías de red:

  • Comparaciones Inequivalentes: Comparar los límites de Hartree-Fock (HF) de un Hamiltoniano con $a_t \to 0$ con los resultados de AFQMC de una matriz de transferencia con $a_t$ finito es inválido. La matriz de transferencia con $a_t$ finito define una teoría regulada distinta con interacciones inducidas. Cuando se utiliza la misma teoría regulada de manera consistente, el principio variacional se cumple (las energías HF limitan las energías proyectadas).
  • Energías de Correlación: Las pequeñas energías de correlación en materia de neutrones para la interacción de "Elementos Esenciales" (EE) no implican funciones de onda triviales. Un análisis detallado de las densidades irreducibles revela una superfluidez multimodal significativa (correlaciones de onda s, onda p y de cuarteto) a pesar de las pequeñas diferencias de energía.
  • Saturación Nuclear: El aplanamiento de la energía de interacción a altas densidades no es un artefacto de empaquetamiento de la red. Derivaciones analíticas en espacio continuo con reguladores invariantes de Galileo muestran que el aplanamiento de la energía de interacción es una consecuencia natural de la regulación en el espacio de momentos. Las propiedades de saturación de las interacciones de NLEFT siguen la escala del regulador inducido por el esparcimiento, no la escala de ocupación de los sitios de la red.
  • Límite Continuo: La afirmación de que las interacciones de red conducen al colapso del sistema en el límite continuo es refutada. El argumento de escala utilizado en la crítica se basa en una renormalización incorrecta de nivel Born; una renormalización adecuada hacia observables de dispersión físicos produce un comportamiento de escala diferente ($C_c \propto a$) que evita el colapso.

Significado y Reivindicaciones
El artículo afirma establecer una base sólida para el análisis completo de errores en las simulaciones de NLEFT con fuerzas quirales completas. Al cuantificar rigurosamente las incertidumbres y demostrar que los métodos libres de problema de signo pueden alcanzar una precisión sub-porcentual para núcleos pesados, este trabajo valida a NLEFT como una referencia fiable para los métodos de muchos cuerpos nucleares.

Los autores sostienen que las aparentes discrepancias en la literatura reciente se resuelven cuando las comparaciones se realizan entre teorías reguladas equivalentes con parámetros, volúmenes y condiciones de contorno consistentes. La acción LAT-OPT1 se presenta no meramente como una herramienta fenomenológica, sino como un modelo de referencia no perturbativo y exactamente resoluble que captura la dinámica esencial de la estructura de capa. Este marco permite la separación de las deficiencias genuinas de la fuerza nuclear de los errores del resolvedor de muchos cuerpos, proporcionando una vía para estudios predictivos de la evolución de las capas, la espectroscopía y la materia densa en núcleos de masa media y pesados.