Las estrellas de neutrones con mezcla de materia oscura se estudian mediante un formalismo de dos fluidos usando ecuaciones de estado agnósticas para la materia hadrónica y la materia oscura. La materia oscura se modela como un gas de Fermi definido solo por su ecuación de estado a baja densidad y su masa. La materia hadrónica se ancla a bajas densidades mediante teoría de campo efectivo quiral y a altas densidades por QCD perturbativa; una parametrización de la velocidad del sonido cubre densidades intermedias para hadrones y altas densidades para materia oscura, asegurando consistencia termodinámica sin sesgos. La materia oscura ligera forma halos extendidos que aumentan la deformabilidad de marea, mientras que la materia oscura pesada forma núcleos compactos que la disminuyen. Las restricciones observacionales se desplazan entonces desde la deformabilidad de marea de ondas gravitacionales para materia oscura ligera a datos de masa-radio de NICER para materia oscura pesada. A 1σ, la fracción de materia oscura está limitada a f_DM ≲ 0.11 para materia oscura ligera. Este marco agnóstico produce límites conservadores y ampliamente aplicables, y las estrellas de neutrones con masas similares pero deformabilidades de marea muy diferentes podrían ser una firma irrefutable de materia oscura.
Las estrellas de neutrones —restos de supernovas— son de los objetos más densos del Universo. Su densidad central supera varias veces la nuclear, y su gravedad solo es superada por los agujeros negros. Por eso, las estrellas de neutrones son un laboratorio ideal para probar la física fundamental: desde el comportamiento de la materia a densidades supernucleares hasta las propiedades de la materia oscura. La materia oscura, predicha por Vera Rubin y originada tras el Big Bang, constituye un cuarto de la energía del Universo, pero aún no se ha detectado directamente. Las estrellas de neutrones pueden capturar partículas de materia oscura, lo que altera sus parámetros observables —masa, radio y deformabilidad de marea, medida mediante ondas gravitacionales. Datos actuales de púlsares, descubiertos por Jocelyn Bell Burnell y estudiados con espectroscopía de rayos X (misión NICER), y de detectores de ondas gravitacionales (LIGO, Rainer Weiss y otros) permiten imponer límites a la fracción de materia oscura en estos objetos, sin atarse a modelos concretos más allá del Modelo Estándar.
Para evitar sesgos, los científicos aplicaron un enfoque agnóstico a las ecuaciones de estado de la materia nuclear y la oscura. La parte nuclear se ancla en cálculos de teoría de campo efectiva quiral a bajas densidades y en cromodinámica cuántica perturbativa a altas, mientras que la región intermedia se parametriza con perfiles aleatorios de la velocidad del sonido (c²_s), limitados solo por la causalidad (la velocidad del sonido no puede superar la velocidad de la luz) y la estabilidad termodinámica. Para la materia oscura se usó una parametrización similar, pero sin anclaje a altas densidades: por debajo de 0.1 de la densidad nuclear se describió como un gas de Fermi degenerado con masa de partícula m_D, y por encima, mediante interpolación de la velocidad del sonido. La estrella se trató como un sistema de dos fluidos: los dos componentes interactúan solo por gravedad, y para cada uno se resuelve su propia ecuación de equilibrio hidrostático (la ecuación de Tolman-Oppenheimer-Volkoff modificada). Se generaron conjuntos de ~10⁵ secuencias masa-radio para distintos m_D (0.2–1.1 GeV) y fracciones de materia oscura f_DM (0.01–0.15), que luego se compararon con datos observacionales.
El análisis mostró que la influencia de la materia oscura depende radicalmente de la masa de sus partículas. La materia oscura ligera (m_D ~ 0.2–0.3 GeV) forma un halo difuso alrededor de la estrella, aumentando su radio efectivo. Esto incrementa la deformabilidad de marea Λ (proporcional a R⁵), y estas configuraciones están fuertemente restringidas por los datos del evento de ondas gravitacionales GW170817. Por el contrario, la materia oscura pesada (m_D > 0.7 GeV) se concentra en un núcleo compacto, reduciendo el tamaño de la estrella y disminuyendo Λ. Aquí, las principales restricciones provienen de las mediciones de masa y radio de NICER. Para masas de partículas intermedias (0.5–0.7 GeV) son posibles ambas opciones, así como estructuras de transición 'núcleo-halo'. El análisis estadístico con todos los datos disponibles (NICER, GW170817 y el requisito de masa máxima ≥ 2.01 M⊙) permitió establecer que, para materia oscura ligera, la fracción permitida no supera el 11%. Para la pesada, las restricciones son más fuertes a altas fracciones, pero en general el espacio de parámetros sigue siendo mucho más amplio. Una predicción intrigante es que el hallazgo de dos estrellas de neutrones con masas similares pero deformabilidades de marea muy distintas podría ser una 'prueba irrefutable' —un indicio de la presencia de materia oscura en una de ellas.
Los resultados obtenidos son independientes del modelo y aplicables a una amplia clase de teorías de materia oscura. Demuestran que las observaciones multimensajero actuales pueden limitar significativamente el contenido de materia oscura en estrellas de neutrones incluso sin conocer su naturaleza microscópica. El tipo de restricción dominante —deformabilidad de marea o masa-radio— actúa como herramienta diagnóstica: el predominio de las restricciones de ondas gravitacionales apunta a estructuras de halo, mientras que el de rayos X indica configuraciones de núcleo de materia oscura.
En el futuro, nuevas mediciones de NICER para otros púlsares, así como datos de detectores de ondas gravitacionales de próxima generación (incluidos los espaciales), permitirán refinar estos límites. El desarrollo ulterior del método agnóstico podría incluir materia oscura bosónica y la rotación estelar, acercando los modelos a la realidad. También es prometedor estudiar la dependencia de la densidad de anclaje elegida para el sector oscuro.
Los resultados influirán en la física de la materia oscura, la astrofísica de estrellas de neutrones y la astronomía de ondas gravitacionales, proporcionando una nueva herramienta para buscar candidatos ligeros a partículas de materia oscura.
Los próximos pasos incluyen el análisis de datos con mayor precisión de futuras misiones como STROBE-X o Athena, y la repetición del estudio para materia oscura bosónica.
La investigación está directamente relacionada con el problema no resuelto de la naturaleza de la materia oscura y la cuestión de la composición de las estrellas de neutrones a densidades extremas, donde podría surgir materia exótica.
🎯 Una cucharadita de materia de una estrella de neutrones pesaría alrededor de mil millones de toneladas, aproximadamente lo mismo que el monte Everest.