Exact Neutron-Proton Wavefunctions Using the Phase Function Method

Resumen Técnico: Funciones de Onda Exactas Neutrón-Protón mediante el Método de la Función de Fase

Planteamiento del Problema
El objetivo principal de la física nuclear teórica en problemas de dispersión es, a menudo, derivar la función de onda, a partir de la cual se pueden extraer otros observables como los desplazamientos de fase, las secciones eficaces y los parámetros de baja energía. Sin embargo, los experimentales no miden directamente las funciones de onda; miden secciones eficaces diferenciales y totales, que están relacionadas con los desplazos de fase. Aunque existen diversas técnicas de aproximación (por ejemplo, la aproximación de Born, la aproximación de Brysk) para resolver la ecuación de Schrödinger para los desplazos de fase, este trabajo se centra en obtener funciones de onda radiales, funciones de amplitud y desplazos de fase dependientes de la distancia exactos para la dispersión neutrón-protón ($n-p$). El estudio aborda específicamente los canales no acoplados S, P y D, con el fin de proporcionar una descripción transparente y precisa de estos estados utilizando el Método de la Función de Fase (PFM, por sus siglas en inglés) combinado con potenciales inversos derivados de la función de Morse.

Metodología
El estudio emplea el Método de la Función de Fase (PFM), una técnica que transforma la ecuación diferencial de segundo orden de Schrödinger en una ecuación diferencial de primer orden de tipo Riccati para el desplazamiento de fase radial, $\delta_\ell(k, r)$.

1. Modelo de Potencial: La interacción se modela mediante el potencial de Morse, $V(r) = V_0(e^{-2(r-r_m)/a_m} - 2e^{-(r-r_m)/a_m})$. Este potencial fue elegido por su solvabilidad exacta, su invarianza de forma y la existencia de una expresión analítica exacta para el estado $^1S_0$.
2. Optimización de Parámetros: Los parámetros del potencial de Morse ($V_0$, $r_m$, $a_m$) no fueron ajustados arbitrariamente. En su lugar, se optimizaron utilizando el análisis exhaustivo de ondas parciales de GRANADA, que comprende 6.713 puntos de datos experimentales de desplazos de fase $n-p$ desde 1950 hasta 2013. La optimización minimizó el error cuadrático medio (MSE) entre el modelo y los datos experimentales.
3. Solución Numérica:
   • La ecuación de Riccati para el desplazamiento de fase $\delta_\ell(k, r)$ se resuelve numéricamente mediante el método de Runge-Kutta de 5º orden (RK-5) con la condición inicial $\delta_\ell(0) = 0$.
   • Se utilizan relaciones de recurrencia específicas para adaptar las funciones de Riccati-Bessel y Riccati-Neumann para diferentes estados de momento angular ($\ell = 0, 1, 2$), lo que produce formas específicas de la ecuación PFM para las ondas S, P y D.
   • Una vez determinada la función de fase, la función de amplitud $A_\ell(r)$ y la función de onda radial exacta $u_\ell(r)$ se calculan utilizando ecuaciones diferenciales acopladas de primer orden derivadas dentro del marco del PFM.
4. Alcance: Los cálculos se realizaron para energías de laboratorio $E_{lab} = [1, 10, 50, 100, 150, 250, 350]$ MeV, con una dependencia radial detallada presentada hasta los 5 fm. El estudio se restringe a canales no acoplados ($^1S_0, ^1P_1, ^3P_0, ^3P_1, ^1D_2, ^3D_2$), excluyendo explícitamente los efectos de canales acoplados (como la mezcla $^3S_1-^3D_1$ debida a fuerzas tensoriales).

Contribuciones Clave

• Funciones de Onda Exactas: El artículo proporciona funciones de onda radiales $u(r)$, funciones de amplitud $A(r)$ y funciones de fase $\delta(r)$ explícitas y dependientes de la distancia para varios canales de dispersión $n-p$ no acoplados, en lugar de solo desplazos de fase asintóticos.
• Aplicación del Potencial Inverso: Demuestra la utilidad de utilizar potenciales de Morse inversos, optimizados frente a un conjunto masivo de datos experimentales (GRANADA), como una referencia de orden cero para cálculos de dispersión de alta precisión.
• Validación contra Nijmegen-II: Las funciones de onda calculadas se comparan directamente con los resultados de la interacción de alta precisión Nijmegen-II, sirviendo como un punto de referencia para la exactitud del enfoque PFM.

Resultados
Los resultados se presentan para energías de laboratorio de 10, 50, 150 y 250 MeV a través de seis canales no acoplados:

• Estado $^1S_0$: El desplazamiento de fase radial $\delta_0(r)$ muestra un comportamiento dependiente de la energía, acumulándose gradualmente a bajas energías pero exhibiendo una reducción en la acumulación y una reversión parcial a radios pequeños conforme aumenta la energía, lo que refleja el núcleo repulsivo de corto alcance. Las funciones de onda muestran un excelente acuerdo con Nijmegen-II, particularmente a energías más altas, con un emparejamiento asintótico suave.
• Estado $^1P_1$: Se encuentra que el potencial inverso es puramente repulsivo. En consecuencia, el desplazamiento de fase permanece negativo en todo el rango de interacción. Las funciones de onda exhiben una amplitud suprimida a distancias cortas, característica de las interacciones repulsivas de onda P, y se alinean bien con los resultados de Nijmegen-II.
• Estados $^3P_0$ y $^3P_1$: Estos canales muestran comportamientos distintos. El estado $^3P_0$ muestra una transición de desplazos de fase negativos a positivos, indicando una interacción entre la repulsión de corto alcance y la atracción de rango intermedio. El estado $^3P_1$ está dominado por la repulsión, manteniendo un desplazamiento de fase negativo. Ambos muestran un buen acuerdo con Nijmegen-II, con una precisión que mejora a energías más altas.
• Estados $^1D_2$ y $^3D_2$: Debido a la barrera centrífuga, las ondas D son menos sensibles al núcleo de corto alcance. Ambos estados son predominantemente atractivos, lo que resulta en desplazos de fase positivos. Las funciones de onda muestran un muy buen acuerdo con Nijmegen-II, con discrepancias menores a distancias cortas atribuidas a la forma simplificada del potencial de Morse en lugar del método PFM en sí.

Significancia y Reivindicaciones
El artículo afirma que el Método de la Función de Fase, cuando se combina con potenciales de Morse inversos optimizados frente a datos experimentales, proporciona un marco fiable, preciso y físicamente transparente para describir los estados de dispersión neutrón-protón no acoplados.

• Exactitud: Las funciones de onda obtenidas muestran un "excelente acuerdo" con los resultados de alta precisión de Nijmegen-II, validando el enfoque PFM para generar funciones de onda de dispersión realistas.
• Transparencia: El método ofrece una interpretación física clara de cómo las regiones atractivas y repulsivas del potencial contribuyen al desplazamiento de fase total (acumulación positiva vs. negativa).
• Robustez: La ausencia de inestabilidades numéricas u oscilaciones no físicas en las funciones de amplitud confirma la estabilidad de la formulación PFM de primer orden.
• Limitaciones y Trabajo Futuro: Los autores señalan modestamente que el presente trabajo se restringe a canales no acoplados. Reconocen que los efectos de canales acoplados (por ejemplo, las interacciones tensoriales en el canal $^3S_1-^3D_1$) no están incluidos y requieren una formulación de matriz más compleja. La extensión a canales acoplados y otros sistemas de dos cuerpos (por ejemplo, $n\alpha, p\alpha$) se identifica como una continuación natural de este estudio.