Impact of neutron captures on the r-process: I. Propagation of statistical model uncertainties

El universo es una vasta fábrica química, pero durante décadas, los científicos han luchado por explicar cómo se forjaron aproximadamente la mitad de todos los elementos más pesados que el hierro. Estos elementos pesados, que incluyen desde el oro hasta el uranio, se crean en los eventos cósmicos más violentos imaginables, como la colisión de dos estrellas de neutrones. En estos entornos extremos, los núcleos atómicos son bombardeados por un flujo de neutrones, construyendo rápidamente átomos cada vez más pesados en un proceso conocido como el proceso de captura rápida de neutrones. Aunque los astrónomos han confirmado que estas colisiones ocurren y producen los tipos correctos de elementos, la física precisa que rige cómo cambian los átomos durante la explosión sigue siendo en gran medida un misterio. Debido a que no podemos recrear estas condiciones en un laboratorio, los investigadores deben confiar en complejos modelos computacionales para simular el viaje de la materia desde un átomo simple hasta un elemento pesado. La precisión de estas simulaciones depende enteramente de los datos introducidos en ellas, específicamente de cómo entendemos los estados de energía interna de los átomos inestables y ricos en neutrones que existen solo por un instante fugaz.

Un equipo de investigadores liderado por Yu An en la Universidad Técnica de Darmstadt ha realizado una inmersión profunda en una de las piezas más críticas de este rompecabezas: los modelos estadísticos utilizados para predecir cómo estos átomos inestables capturan neutrones. En su estudio, no intentaron inventar nueva física, sino que probaron cuánto cambian los resultados finales de una explosión cósmica cuando los supuestos subyacentes sobre la estructura atómica se alteran ligeramente. Se centraron en dos ingredientes teóricos específicos: la densidad de niveles nucleares, que describe de cuántas formas diferentes puede un núcleo atómico vibrar o almacenar energía, y la función de fuerza de rayos gamma, que mide la probabilidad de que el núcleo emita energía lumínica mientras se estabiliza. Al ejecutar miles de simulaciones con diferentes combinaciones de estos modelos teóricos, el equipo mapeó exactamente dónde nuestra falta de conocimiento crea las brechas más grandes en nuestra comprensión de la historia química del universo.

Los investigadores comenzaron simulando las secuelas de la fusión de dos estrellas de neutrones, un escenario donde dos estrellas muertas chocan y eyectan una nube de material superdenso. Ejecutaron esta simulación utilizando sesenta combinaciones diferentes de los modelos teóricos disponibles, creando una amplia gama de posibles resultados para la mezcla final de elementos. Los resultados revelaron una jerarquía clara de incertidumbre. Los modelos que describen cómo los núcleos almacenan energía, conocidos como densidades de niveles nucleares, tuvieron un impacto mucho mayor en la composición química final que los modelos que describen cómo emiten luz. Mientras que los picos principales de elementos pesados, como los correspondientes al oro y al platino, permanecieron relativamente estables independientemente del modelo utilizado, las regiones entre estos picos mostraron diferencias dramáticas. En algunas simulaciones, la abundancia de elementos en un rango de masa específico alrededor de 172 a 178 unidades atómicas se triplicaba o desaparecía por completo dependiendo de qué modelo de energía se eligiera.

Para entender por qué sucede esto, el equipo observó de cerca la cronología de los eventos dentro de la nube de escombros en expansión. Tras la explosión inicial, la temperatura desciende rápidamente. Al principio, los átomos están tan calientes que los neutrones son capturados y luego expulsados constantemente en una danza equilibrada. Eventualmente, la temperatura cae lo suficiente como para que este equilibrio se rompa y los átomos comiencen a establecerse en sus formas finales. Los investigadores descubrieron que la elección del modelo de energía se vuelve crítica en la fase "fría" de este proceso, que ocurre cuando la temperatura desciende a unos 0.45 mil millones de grados. En esta ventana, si el modelo teórico predice que los átomos capturan neutrones lentamente, el material permanece en una región específica de la tabla periódica, creando un subpico distintivo de elementos. Si el modelo predice una captura más rápida, ese material pasa de largo para convertirse en elementos aún más pesados, dejando un vacío tras de sí. Esta sensibilidad no es solo una peculiaridad de una sola simulación; el equipo encontró que incluso al promediar los resultados de más de 11,000 escenarios de explosión diferentes para imitar una mezcla galáctica realista, estos huecos y picos específicos persistían.

El estudio también exploró una vasta cuadrícula de condiciones hipotéticas, variando la densidad inicial y la relación de protones a electrones en el material para ver cómo reaccionan diferentes entornos a estas incertidumbres. Descubrieron que el impacto de los modelos de energía no es uniforme en todos los eventos cósmicos. En entornos donde el suministro de neutrones es escaso, la elección del modelo importa muy poco porque los átomos simplemente se quedan sin combustible antes de que puedan ser significativamente alterados. Sin embargo, en entornos con un suministro rico de neutrones y condiciones termodinámicas específicas, la diferencia entre modelos puede determinar si un elemento particular se produce en abundancia o apenas existe. Los investigadores identificaron que las incertidumbres más significativas surgen en escenarios "fríos", donde los átomos continúan capturando neutrones mucho después de que el equilibrio inicial se haya roto. En estas condiciones, pequeños cambios en la velocidad predicha de captura de neutrones pueden desplazar todo el flujo de materia, decidiendo si el proceso se detiene en el segundo pico de elementos pesados o si avanza hacia el tercero.

En última instancia, este trabajo sirve como un mapa crucial para futuros experimentos. Indica a los físicos nucleares exactamente qué regiones de la tabla atómica son más sensibles a nuestras brechas teóricas actuales. El estudio confirma que, si bien tenemos una imagen general de cómo se fabrican los elementos pesados, los detalles finos del paisaje químico dependen en gran medida de qué tan bien entendamos los estados de energía interna de los núcleos inestables. Los investigadores concluyen que, para refinar nuestra comprensión de la evolución química del universo, necesitamos mejores datos experimentales y descripciones teóricas más precisas de las densidades de niveles nucleares, particularmente para los isótopos específicos que se encuentran entre los picos principales de los elementos pesados. Sin estas mejoras, nuestras simulaciones del cosmos continuarán cargando con una sombra de incertidumbre, incapaces de explicar plenamente la abundancia precisa de los elementos pesados que conforman nuestro mundo.