Detection of a dark matter subhalo in the strongly lensed system PJ011646

Resumo Técnico: Detecção de um Subhalo de Matéria Escura no Sistema Fortemente Lenteado PJ011646

Declaração do Problema
A natureza da matéria escura permanece um problema fundamental não resolvido na cosmologia. Enquanto o modelo de Matéria Escura Fria (CDM) prevê uma formação hierárquica de halos auto-ligados que se estendem a massas muito baixas, os modelos de Matéria Escura Quente (WDM) preveem uma supressão de halos de baixa massa abaixo de uma escala de corte característica. Distinguir entre esses modelos requer medições estatísticas da população de subhalos, particularmente no regime de massa onde as previsões divergem ($M_{200} \sim 10^7 - 10^9 M_\odot$). Um desafio primário na detecção dessas subestruturas via lente gravitacional forte é a "degenerescência" entre sinais genuínos de subhalos e complexidade angular não modelada na distribuição de massa suave da galáxia lente. Modelos simplificados de lei de potência elíptica (PL) frequentemente falham em capturar desvios em galáxias de tipo precoce (ETGs), levando a detecções espúrias de subhalos ou ao mascaramento de reais. Além disso, detecções anteriores foram limitadas pela resolução angular de telescópios ópticos (por exemplo, HST) e pela contaminação pela luz da galáxia lente.

Metodologia
Os autores apresentam uma análise de lenteamento forte do sistema PJ011646 (lente $z=0.555$, fonte $z=2.125$) utilizando observações de continuum de poeira de alta resolução ($\sim 0.1$ arco-segundo) do Array Large Millimeter/submillimeter de Atacama (ALMA) na Banda 7. A análise utiliza o conjunto de software de código aberto PyAutoLens, PyAutoFit e PyAutoGalaxy.

A metodologia prossegue em três etapas principais:

1. Modelagem Macro: A distribuição de massa da lente é modelada usando um perfil de lei de potência elíptica (EPL) com cisalhamento externo. Para abordar as limitações de modelos elípticos simples, os autores incorporam multipolos de Fourier de ordem superior (ordens $m=3$ e $m=4$, e adicionalmente $m=1$) para capturar desvios angulares na densidade de massa. A convergência desses multipolos é parametrizada para evitar degenerescências associadas a fases periódicas.
2. Reconstrução da Fonte: O brilho superficial da fonte é reconstruído usando um método de inversão bayesiana pixelada (Suyu et al. 2006) com pixelização em malha de Voronoi e interpolação de vizinho natural. Isso é realizado diretamente no plano $uv$ para aproveitar os dados interferométricos completos.
3. Busca por Subhalos: Uma busca baseada em grade é realizada no plano da imagem para detectar perturbações. Um perfil esférico de Navarro-Frenk-White (NFW) é adicionado como perturbador em cada célula da grade, e a diferença de evidência bayesiana logarítmica ($\Delta \ln Z$) é calculada em relação ao melhor modelo macro ajustado. Um limiar de $\Delta \ln Z > 10$ é adotado para identificar candidatos robustos. Os autores também testam um perfil NFW generalizado (gNFW) com uma inclinação interna livre ($\gamma$) para verificar desvios do perfil CDM padrão.

Principais Resultados

• Melhoria do Modelo Macro: A inclusão de multipolos de terceira e quarta ordem ($m=3, 4$) melhora significativamente o ajuste aos dados, resultando em um aumento de evidência bayesiana de $\Delta \ln Z = 52.1$ em relação a um modelo PL puro. As amplitudes desses multipolos são $\simeq 1.5\%$ da convergência, consistente com observações de ETGs próximas. Um modelo incluindo um multipolo de primeira ordem ($m=1$) fornece uma melhoria adicional marginal ($\Delta \ln Z = 6.5$), sugerindo um companheiro físico possível, mas não impulsionando o sinal primário.
• Detecção de Subhalo: Usando o melhor modelo macro PL+MP ajustado, a busca por grade identifica duas regiões no plano da imagem com $\Delta \ln Z > 10$. No entanto, uma análise de acompanhamento revela que essas duas regiões mapeiam para a mesma localização no plano da fonte, indicando uma degenerescência. Um modelo de dois subhalos não melhora a evidência sobre um modelo de único subhalo, e a posterior para o segundo candidato produz concentrações não físicas ($c_{200} > 1000$). Consequentemente, os autores concluem que apenas um subhalo é fisicamente real, localizado à direita do arco norte.
• Propriedades do Subhalo: O subhalo detectado tem uma massa $M_{200} = 2.78^{+0.43}_{-0.66} \times 10^{10} M_\odot$ e uma concentração $c_{200} = 30^{+5}_{-7}$. Esta concentração é consistente com expectativas para um subhalo NFW despojado por maré dentro de um halo pai. A significância da detecção é de $\sim 5.8\sigma$ em relação ao modelo macro PL+MP.
• Restrições de Massa Robustas: Devido à degenerescência entre o raio de escala NFW e a normalização, a massa total $M_{200}$ não é estritamente restringida. Em vez disso, os autores definem um "raio robusto" ($R_{sub} \approx 2$ kpc) onde a incerteza fracionária na massa contida projetada é minimizada. Dentro deste raio, a massa contida é $M_{sub} = 3.57^{+0.16}_{-0.14} \times 10^9 M_\odot$.
• Limites de Sensibilidade: Em células da grade onde nenhum subhalo é detectado ($\Delta \ln Z < 10$), os autores derivam limites superiores para massas de subhalos detectáveis. Nas regiões mais sensíveis (ao longo dos arcos lenteados brilhantes), a massa mínima detectável com significância de $3\sigma$ é $M_{200} \approx 8 \times 10^8 M_\odot$, assumindo a relação massa-concentração de Ludlow et al. (2016).

Significância e Alegações
O artigo afirma que este trabalho demonstra a capacidade da imagem de continuum do ALMA em resolução sub-arco-segundo para sondar subestrutura de matéria escura em um regime de massa ($M_{200} \sim 10^9 M_\odot$) onde os modelos CDM e WDM preveem abundâncias diferentes. Ao modelar explicitamente a complexidade angular da galáxia lente via multipolos, os autores mitigam incertezas sistemáticas que historicamente afetaram buscas por subhalos. A detecção de um subhalo com propriedades consistentes com previsões do CDM, juntamente com a derivação de limites superiores rigorosos em regiões de não detecção, fornece um teste robusto da formação de estrutura em pequena escala. Os autores notam que, embora a resolução angular seja comparável à do HST, a alta magnificação do sistema e a ausência de contaminação por luz de primeiro plano nos dados do ALMA permitem a sondagem de massas mais baixas do que anteriormente alcançável com dados ópticos. O estudo posiciona-se como um passo em direção a buscas sistemáticas e estatisticamente significativas por subestrutura de matéria escura usando a amostra crescente de galáxias formadoras de estrelas empoeiradas (DSFGs) fortemente lenteadas observadas com o ALMA.