A Maximum-Entropy Method for Zero-Skewness Valence GPDs Constrained by Nucleon Electromagnetic Form Factors

Resumen Técnico: Un Método de Máxima Entropía para GPDs de Valencia de Sesgo Cero Restringidos por Factores de Forma Electromagnéticos del Nucleón

Planteamiento del Problema
La estructura espacial y de momento del nucleón se describe mediante dos conjuntos complementarios pero incompletos de observables: los factores de forma (FF) electromagnéticos elásticos y las funciones de distribución de partones (PDF). Los FF codifican la distribución espacial transversal del nucleón, pero integran sobre todas las fracciones de momento de los partones ($x$), mientras que las PDF describen la distribución de momento longitudinal pero carecen de resolución espacial transversal. Las Distribuciones de Partones Generalizadas (GPDs) unifican estas descripciones, aunque la reconstrucción de la estructura transversal dependiente de $x$ de los quarks de valencia a partir únicamente de los FF elásticos constituye un problema inverso severamente subdeterminado. Específicamente, un único FF restringe solo los momentos integrados en $x$ de las GPDs, dejando arbitrario el perfil transversal local dependiente de $x$. El desafío consiste en construir un marco de GPDs de valencia de sesgo cero, físicamente motivado, estable y de menor sesgo, que satisfaga cuatro restricciones de FF nucleares ($F_1^p, F_1^n, F_2^p, F_2^n$) y normalizaciones de límite frontal sin introducir oscilaciones no físicas o estructuras espurias.

Metodología
Los autores formulan un marco de Método de Máxima Entropía (MEM) de perfil reducido para resolver este problema inverso. El enfoque se basa en los siguientes componentes clave:

1. Ansatz y Restricciones: Las GPDs de valencia de sesgo cero, $H_v^q(x,t)$ y $E_v^q(x,t)$, se modelan utilizando un ansatz de perfil en $t$ exponencial:
   $$H_v^q(x, t) = q_v(x) \exp[t f_H^q(x)], \quad E_v^q(x, t) = e_v^q(x) \exp[t f_E^q(x)]$$
   Aquí, $q_v(x)$ y $e_v^q(x)$ son las PDF de valencia del límite frontal y las distribuciones de tipo Pauli del límite frontal, respectivamente. Las funciones desconocidas $f_H^q(x)$ y $f_E^q(x)$ representan las funciones de perfil transversal positivas. Las restricciones incluyen:

   • Momentos Elásticos: Los primeros momentos de Mellin de las GPDs deben reproducir los factores de forma Dirac y Pauli separados por sabor derivados de los FF experimentales del protón y el neutrón.
   • Normalizaciones del Límite Frontal: En $t=0$, las GPDs deben reducirse a las PDF de entrada ($q_v$) y satisfacer las normalizaciones del momento magnético anómalo de sabor ($\kappa^q$).

2. Variedad de Perfil Reducido: Para regularizar el problema inverso subdeterminado, los autores evitan una inversión directa basada en rejilla de las funciones de perfil. En su lugar, restringen los perfiles a una variedad positiva de baja dimensión utilizando una forma funcional específica:
   $$f(x) = 0.05 + (1-x)^2 \exp(c_0 + c_1 x + c_2 x^2)$$
   Esto reduce los grados de libertad de un gran conjunto de puntos a 12 parámetros de forma (tres coeficientes para cada uno de los cuatro perfiles: $f_H^u, f_H^d, f_E^u, f_E^d$). Esta reducción suprime los modos locales que los momentos elásticos no pueden restringir.

3. Regularización por Entropía: Se maximiza el funcional de MEM $Q = \alpha S - \frac{1}{2}\chi^2$. El término de verosimilitud $\chi^2$ mide la desviación entre los FF reconstruidos y los datos de entrada. El término de entropía $S$ (Shannon-Jaynes) penaliza las desviaciones de un modelo por defecto suave y positivo, asegurando que la solución sea el perfil de "menor sesgo" compatible con las restricciones. El hiperparámetro $\alpha$ equilibra la fidelidad de los datos frente a la suavidad previa.

4. Implementación Numérica: El cálculo utiliza parametrizaciones de la función beta analítica para las entradas frontales y una parametrización suave basada en dipolos para los FF elásticos como prueba de concepto. La optimización se realiza mediante el algoritmo L-BFGS-B en una rejilla discreta, con las entradas frontales normalizadas en la misma rejilla para asegurar la satisfacción exacta de las reglas de suma con precisión de máquina en $t=0$.

Resultos Clave

• Estabilidad de la Reconstrucción: La solución de perfil reducido es altamente estable bajo variaciones del peso de la entropía $\alpha$ (probado en el rango de $10^{-4}$ a $10^{-2}$). La desviación máxima puntual en las funciones de perfil permanece por debajo del 0.01%, lo que indica que la solución está impulsada por las restricciones de los factores de forma más que por el modelo previo de entropía.
• Reproducción de Factores de Forma: Las GPDs reconstruidas reproducen los factores de forma de Sachs de entrada ($G_E^p, G_M^p, G_E^n, G_M^n$) al nivel del percentual en el rango $0 \le Q^2 \le 5 \text{ GeV}^2$. Las normalizaciones frontales para los números de quarks de valencia y los momentos magnéticos anómalos se satisfacen exactamente en la rejilla discreta.
• Interpretación Física: Los perfiles resultantes $f^q(x)$ son positivos y disminuyen conforme $x \to 1$. Esto conduce a distribuciones de parámetro de impacto donde los partones de gran $x$ están más localizados cerca del centro de movimiento de la parte transversal (encogimiento transversal), lo cual es consistente con las expectativas físicas.
• Verificaciones Cualitativas: Se construyó un proxy de la asimetría de espín de haz ($A_{LU}$) en la electroproducción exclusiva de fotones utilizando los perfiles reconstruidos. El proxy reproduce con éxito el patrón azimutal $\sin \phi$ dominante observado en los datos de CLAS12, proporcionando una verificación de consistencia cualitativa más allá de las restricciones elásticas.

Significancia y Reivindicaciones
El artículo afirma proporcionar una línea base de sesgo cero controlada para conectar las restricciones de los factores de forma elásticos con los perfiles transversales dependientes de $x$. Los autores enfatizan que la novedad no reside en el ansatz de perfil en $t$ exponencial en sí mismo (que es estándar), sino en la formulación del problema de los cuatro factores de forma como un único proceso de MEM acoplado. Este proceso integra las restricciones de momentos elásticos, las normalizaciones de PDF frontales y las restricciones de momento anómalo dentro de un regularizador basado en la entropía que preserva la positividad.

El trabajo se presenta como una prueba de concepto que demuestra que se puede extraer un perfil transversal de baja dimensión, físicamente motivado, a partir de datos elásticos e información frontal sin sobreajuste. Los autores declaran explícitamente que los resultados son condicionales a la variedad de perfil reducido elegida y a las entradas analíticas específicas utilizadas. No pretenden realizar una extracción fenomenológica global, sino más bien un punto de partida estable para futuros análisis que puedan incorporar ajustes empíricos de factores de forma con información de covarianza, conjuntos de PDF modernos, entradas de QCD en la red y observables exclusivos duros. El método está diseñado para ser agnóstico a la entrada, permitiendo que la maquinaria de reconstrucción se aplique a diferentes parametrizaciones sin modificaciones estructurales.