El estudio resuelve la discrepancia relacionada con el exceso de galaxias a z≥10 detectadas por JWST en comparación con el modelo estándar ΛCDM. Se considera un escenario con un pico (bump) en el espectro de potencia de las perturbaciones a una escala de ~10¹⁰ M⊙, que sin embargo desplaza la época de reionización a un pasado demasiado temprano. Se demuestra que añadiendo supresión en escalas menores (mediante reducción de la amplitud de las perturbaciones primordiales o introduciendo materia oscura tibia) se reduce el número de halos de baja masa y, en consecuencia, la cantidad de radiación UV ionizante, resolviendo la contradicción. Utilizando la teoría extendida de Press-Schechter (EPS) se obtuvo la distribución de las tasas de fusión de halos en presencia del pico y se verificó con simulaciones directas de N cuerpos. A partir de esta distribución y la suposición de que las fusiones desencadenan estallidos de formación estelar, se calculó la función de luminosidad UV de las galaxias tempranas, consistente con los datos observacionales. En este modelo, la fracción de galaxias que se forman mediante un mecanismo monolítico (colapso a gran escala único) aumenta considerablemente, y las primeras estrellas aparecen directamente en halos con masas ≥10⁸ M⊙, denominándose "estrellas de población IV" (con composición química primordial).
El modelo cosmológico estándar del Big Bang con materia oscura fría (ΛCDM) describe con éxito la estructura a gran escala, pero enfrenta dificultades a escalas galácticas. Tras el lanzamiento del James Webb, se descubrió que a z>10 hay aproximadamente entre 10 y 100 veces más galaxias masivas de lo predicho. Esto obliga a reconsiderar ya sea las propiedades de la materia oscura, o la forma del espectro de perturbaciones primordiales proveniente de la época inflacionaria. El cumplimiento simultáneo de tres condiciones —el exceso de galaxias, la concordancia con la época de reionización según los datos del fondo cósmico de microondas y el ajuste de la función de luminosidad UV— se convierte en una prueba crucial para las nuevas teorías.
Los autores parametrizaron el espectro de potencia de las perturbaciones de densidad como el espectro ΛCDM estándar multiplicado por un «pico» gaussiano con amplitud A=20 y escala característica k0=4.69 Mpc⁻¹, y añadieron una supresión exponencial a escalas menores que Mfs=10⁹ M⊙ (imitando el efecto de libre vuelo de las partículas de materia oscura caliente). Dentro del formalismo extendido de Press-Schechter (EPS), se calcularon las funciones de masa de halos y las probabilidades condicionales de fusión. Para comprobar las fórmulas analíticas de la tasa de acreción, se realizaron simulaciones de N-cuerpos con GADGET-2 para espectros con el pico. La función de luminosidad UV se obtuvo vinculando la tasa de formación estelar, inducida por la colisión de discos de gas durante la fusión, con la acreción de halos pequeños sobre los grandes, y teniendo en cuenta la baja probabilidad (~3%) de que los discos colisionen directamente en un tiempo dinámico. Se utilizaron datos fotométricos del JWST y confirmaciones espectroscópicas para la comparación.
El modelo con el «pico» y el corte produce, a z~9, un aumento en el número de halos masivos (M~10¹⁰ M⊙) de aproximadamente un orden de magnitud, y a z~12, de dos órdenes en comparación con ΛCDM, lo que resuelve directamente el problema del exceso de galaxias. Además, la supresión de las perturbaciones a pequeña escala evita la formación demasiado temprana de estrellas y la ionización: como se muestra en la figura 3, la curva de ionización x(z) prácticamente coincide con la estándar, y la ionización del 50% se alcanza a z~10, valor cercano a 7.68±0.79 según los datos de Planck del fondo cósmico de microondas. La fórmula analítica para la distribución de la tasa de crecimiento de masa de los halos, derivada del EPS con una corrección por tiempo dinámico, concuerda bien con las simulaciones (fig. 4). A partir de esta distribución, las funciones de luminosidad UV calculadas (fig. 6) para z=7, 12 y 17 reproducen los puntos observacionales con los parámetros de eficiencia de formación estelar indicados en la tabla. El modelo predice que la fracción de galaxias que se forman de manera monolítica es notablemente mayor, y que las primeras estrellas (población IV) nacen directamente en halos masivos con masas del orden de 10⁸ M⊙ y superiores, sin pasar por una etapa de enriquecimiento previo en metales.
El trabajo demuestra que la introducción de dos características en el espectro de perturbaciones iniciales —un refuerzo a la escala de las protogalaxias y un corte a escalas menores— puede resolver de un plumazo varias tensiones observacionales del modelo ΛCDM: el exceso de galaxias masivas tempranas, el problema del exceso de satélites y el desplazamiento de la época de reionización. Esto refuerza la motivación para buscar desviaciones del espectro de ley de potencias simple que predicen los modelos de inflación más sencillos, y abre la vía a sondear la física inflacionaria y las propiedades de la materia oscura mediante observaciones del Universo lejano. El concepto de «población IV» subraya que la evolución química podría haber comenzado directamente en sistemas masivos, lo cual es importante para interpretar los datos de los primeros agujeros negros supermasivos.
El desarrollo natural de la idea será la construcción de modelos inflacionarios autoconsistentes que generen un espectro con pico y corte, así como simulaciones hidrodinámicas detalladas de la formación de galaxias en estas cosmologías. Resulta de especial interés el estudio de la formación conjunta de las primeras estrellas y las semillas de los agujeros negros supermasivos en nubes que colapsan monolíticamente. Las próximas campañas de observación del JWST y futuros instrumentos (por ejemplo, el Telescopio Espacial Roman) permitirán precisar la función de luminosidad a z>12 y verificar las predicciones de un menor número de galaxias enanas satélite, lo que constituirá una prueba independiente del modelo con corte en el espectro.
El escenario descrito influye directamente en la interpretación de los datos de JWST y Planck, impone nuevas restricciones a los modelos de materia oscura (caliente vs fría) e inflación, así como a la teoría de la evolución química del Universo.
Se requieren simulaciones especializadas de N-cuerpos de gran volumen con un espectro que contenga simultáneamente el pico y el corte, para refinar las tasas de fusión y la distribución espacial de los halos. Paralelamente, es necesario investigar la influencia de la retroalimentación de las estrellas de población IV y la acreción sobre los agujeros negros en las galaxias masivas tempranas.
El modelo propuesto aborda directamente problemas no resueltos de la física: la naturaleza de la materia oscura (partículas con masas del orden de keV o un potencial inflacionario no estándar), el mecanismo de generación de las fluctuaciones semilla durante la inflación, así como la explicación del rápido crecimiento de los agujeros negros supermasivos a z>6. Además, ofrece una visión unificada de la formación «monolítica» y jerárquica de las galaxias.
🎯 El término «estrellas de población IV» fue introducido por los autores para designar una hipotética primera generación de astros nacidos en halos masivos con composición química primordial —en contraste con las clásicas estrellas de población III, que se forman en minihalos con masas ~10⁵–10⁶ M⊙. Si existieron tales estrellas, pudieron ser brillantes pero de corta vida, y sus explosiones enriquecieron rápidamente el medio intergaláctico con metales.