A self-consistent Higgs-portal framework for dark matter--admixed neutron stars: Collider-motivated benchmarks meet multimessenger constraints

Resumo Técnico: Um Modelo de Portal de Higgs Autoconsistente para Estrelas de Nêutrons com Admissão de Matéria Escura

Definição do Problema
A modelagem teórica de estrelas de nêutrons com admissão de Matéria Escura (DANSs) dentro de um framework de fluido único tem dependido historicamente de prescrições fenomenológicas que carecem de autoconsistência astrofísica. Especificamente, abordagens convencionais tratam o momento de Fermi local da DM ($k^\chi_F$) ou a densidade local de DM ($n_\chi$) como parâmetros externamente fixos, desacoplando a distribuição de DM do ambiente bariônico estelar. Esse congelamento ad hoc da densidade de DM negligencia os processos físicos de captura, termalização e assentamento gravitacional dentro da estrela. Além disso, embora o modelo de portal de Higgs seja uma realização de partícula-física amplamente utilizada para interações de DM, sua implementação em estudos de estrelas de nêutrons (NS) frequentemente falha em acoplar dinamicamente os setores de matéria escura (DM) e bariônica (BM), levando a uma descrição incompleta da equação de estado (EOS) e da resposta termodinâmica.

Metodologia
Para abordar essas limitações, os autores desenvolvem um framework de campo médio relativístico (RMF) de fluido único e autoconsistente que estende o modelo padrão de portal de Higgs. A metodologia envolve as seguintes etapas principais:

1. Extensão do Modelo: Os autores introduzem um mediador vetorial neutro massivo, $X_\mu$ (identificado como um bóson dark $Z'$), associado a uma simetria de calibre $U(1)_X$ oculta. Este mediador gera um campo médio vetorial repulsivo dependente da densidade que acopla dinamicamente os setores BM e DM ao modificar seus potenciais químicos. Esta abordagem é análoga ao papel do méson $\omega$ na teoria RMF padrão, mas aplicada ao setor de DM.
2. Acoplamento Autoconsistente: Em vez de prescrever um $k^\chi_F$ fixo, o framework impõe uma fração global de partículas, $F_\chi = N_\chi/N_B$, como o único parâmetro astrofísico livre. A densidade local de DM é derivada de forma autoconsistente como $n_\chi = F_\chi n_B$, garantindo que a distribuição de DM evolua naturalmente com o perfil de densidade bariônica.
3. Benchmarks de Física de Partículas: Os parâmetros do modelo são restringidos por benchmarks motivados por colisores (ATLAS) e limites de detecção direta. A massa da DM é definida como $m_\chi = 200$ GeV (paradigma WIMP), a massa do mediador para $m_X = 1.5$ TeV, e os acoplamentos são fixados em $g_{qX} = 0.25$ e $g_{\chi X} = 1$. O acoplamento Higgs-DM é definido como $g_h = 0.07$.
4. Derivação da Equação de Estado (EOS): Os autores derivam expressões analíticas exatas para a densidade de energia total e pressão, incorporando contribuições do portal escalar de Higgs e da interação vetorial $Z'$. Crucialmente, eles derivam expressões analíticas rigorosas para a velocidade de som adiabática ($c_s^2$), sua derivada de densidade e o índice adiabático ($\Gamma$). Estas derivações contabilizam explicitamente a dependência da densidade de $n_\chi$ e $k^\chi_F$, distinguindo o framework de modelos com $k^\chi_F$ fixo.
5. Implementação Numérica: O framework é aplicado utilizando três EOSs RMF distintas para a matéria bariônica: NL3$\omega\rho$, FSU2R e DD2. As equações de Tolman–Oppenheimer–Volkoff (TOV) são resolvidas para determinar a estrutura macroscópica (massa, raio, energia de ligação) das estrelas para variados valores de $F_\chi$ (0 a 0,2%).

Principais Contribuições

• Formalismo Autoconsistente: O artigo estabelece um framework onde a densidade local de DM não é uma entrada externa, mas um resultado autoconsistente da fração global de partículas e da densidade bariônica local ($n_\chi = F_\chi n_B$). Isso elimina a necessidade de prescrições arbitrárias do momento de Fermi da DM.
• Benchmarks Termodinâmicos Analíticos: Os autores fornecem expressões analíticas exatas para $c_s^2$, $dc_s^2/dn_B$ e $\Gamma$ no contexto de sistemas acoplados BM-DM. Estas expressões identificam as contribuições microscópicas específicas (térmicas, escalares, vetoriais e de interação) para a resposta termodinâmica, oferecendo um benchmark rigoroso para avaliar a causalidade e a estabilidade.
• Decomposição de Interação: O estudo decompõe quantitativamente a energia de interação BM-DM em componentes escalares (Higgs) e vetoriais ($Z'$). Demonstra que, para os benchmarks adotados e motivados por colisores, o canal vetorial repulsivo domina amplamente o canal escalar (por um fator de $>20$), fornecendo uma explicação microscópica para os efeitos macroscópicos observados.
• Conexão Micro-Macro: O trabalho une restrições de física de partículas (limites de colisor para a massa e acoplamentos de $Z'$) com observáveis astrofísicos (relações massa-raio, deformabilidade de maré), criando um framework unificado para restrições multimensageiras.

Resultos

• Suavização da EOS: O aumento da fração de DM ($F_\chi$) suaviza sistematicamente a EOS nuclear. Isso ocorre porque as partículas pesadas de DM ($m_\chi \gg m_N$) contribuem significativamente para a densidade de energia, mas fornecem suporte de pressão negligenciável. Consequentemente, a massa máxima da NS ($M_{\text{max}}$) diminui aproximadamente 25–29% conforme $F_\chi$ aumenta de 0 para 0,2%, sendo a redução mais pronunciada para EOSs mais suaves (FSU2R).
• Compactação Estelar: A presença de DM leva a configurações estelares mais compactas. Para uma estrela canônica de 1,4 $M_\odot$, o raio diminui de 2 a 3 km, e a densidade bariônica central aumenta em até ~20%.
• Resposta Termodinâmica: A velocidade de som adiabática ($c_s^2$) é sistematicamente suprimida, particularmente em altas densidades, devido à contribuição dinâmica do setor de DM para as derivadas termodinâmicas. A derivada de densidade de $c_s^2$ também mostra uma diminuição monotônica com o aumento de $F_\chi$, indicando uma resposta termodinâmica mais gradual.
• Anomalias de Energia de Ligação: O estudo relata que, para altos valores de $F_\chi$, a energia de ligação gravitacional convencional ($E_b = M_B - M$) pode tornar-se negativa. Os autores esclarecem que isso não implica instabilidade hidrostática (que é determinada por $dM/d\varepsilon_c > 0$), mas reflete o desacoplamento do orçamento de massa bariônica ($M_B$) da massa de repouso total ($M_0$), que inclui a substancial massa de repouso dos WIMPs pesados.
• Propriedades Microscópicas: A DM permanece profundamente não-relativística em todo o interior estelar ($k^\chi_F \lesssim 100$ MeV, $\mu_\chi \approx m_\chi$). O deslocamento do potencial químico é impulsionado primariamente pela interação vetorial dependente da densidade, enquanto a contribuição escalar permanece subdominante.

Significância
O artigo afirma que impor a autoconsistência no tratamento do setor escuro é essencial para obter EOSs de fluido único fisicamente significativas. Ao derivar expressões analíticas exatas para quantidades termodinâmicas, o framework fornece um benchmark geral para avaliar a estabilidade e a causalidade da matéria com admissão de DM. O trabalho estabelece um elo direto e quantitativo entre modelos WIMP motivados por colisores e astrofísica multimensageira. Demonstra que estrelas compactas podem servir como laboratórios complementares aos experimentos terrestres, onde restrições observacionais de ondas gravitacionais (ex: GW170817) e tempo de raios-X (ex: NICER) podem ser traduzidas em limites sobre a quantidade de DM acumulada dentro de estrelas de nêutrons para um cenário fixo de física de partículas. Os resultados sugerem que pulsares massivos fornecem um limite superior rigoroso sobre o conteúdo de DM, e que as assinaturas específicas de admissão de WIMPs pesados (como a suavização da EOS e perfis modificados de $c_s^2$) podem ser distinguidas das incertezas da EOS nuclear.