Probing globular clusters parameters through gravitational wave lensing with stellar-mass black hole binaries

Resumo Técnico: Sondagem de Parâmetros de Aglomerados Globulares através de Lente Gravitacional de Ondas Gravitacionais com Binárias de Buracos Negros de Massa Estelar

Definição do Problema
Aglomerados globulares (AGs) são sistemas estelares densos que servem como laboratórios para o estudo da dinâmica estelar e da formação precoce de galáxias. Um parâmetro crítico para compreender sua dinâmica interna, formação e evolução é sua massa, tipicamente inferida a partir da dispersão de velocidade estelar central ($\sigma_v$). Os métodos atuais para determinar a massa de AGs dependem de observações eletromagnéticas (EM), como medições de dispersão de velocidade e perfis de brilho superficial. No entanto, esses métodos enfrentam limitações significativas, incluindo efeitos de projeção na linha de visão, aglomeração estelar em regiões centrais densas, ruído de disparo (shot noise) de estrelas gigantes luminosas e incertezas na identificação do verdadeiro centro do aglomerado. Esses desafios são exacerbados em aglomerados com colapso de núcleo ou sistemas que potencialmente abrigam buracos negros de massa intermediária (IMBHs). Além disso, embora os próximos levantamentos de campo amplo (ex: LSST, Euclid) expandirão os catálogos de AGs com posições e distâncias precisas, medir a cinemática interna como a dispersão de velocidade continua sendo difícil e tipicamente requer acompanhamento espectroscópico dedicado.

Os autores propõem uma sonda complementar: utilizar ondas gravitacionais (GW) de coalescências de binárias compactas (CBC) que são lentiadas por AGs. Diferente da lente forte por galáxias, que produz múltiplas imagens, a lenteção por AGs opera no regime de óptica de onda. Neste regime, a lente não cria múltiplos sinais, mas imprime uma amplificação dependente da frequência no sinal de GW. O artigo investiga se essas assinaturas de óptica de onda podem ser usadas para recuperar as propriedades intrínsecas do AG de lente, especificamente sua dispersão de velocidade central.

Metodologia
O estudo emprega um framework de simulação e estimação bayesiana de parâmetros para testar a recuperabilidade dos parâmetros do AG a partir de sinais de GW lentiados.

1. Modelagem da Lente: Os AGs são modelados como Esferas Isotérmicas Singulares (SIS). A distribuição de massa é aproximada por um perfil de densidade $\rho(r) = \sigma_v^2 / (2\pi G r^2)$, onde a dispersão de velocidade $\sigma_v$ é assumida como constante. A massa de lente efetiva ($M_{SIS}$) é derivada como uma função da dispersão de velocidade $\sigma_v$ e da distância da lente ($D_l$).
2. Simulação de Sinal: Os autores simulam sinais de GW semelhantes ao GW150914 (uma fusão de binária de buracos negros de massa estelar). Esses sinais são lentiados pelo fator de amplificação SIS derivado do integral de difração de Fresnel-Kirchhoff. A amplificação depende da frequência adimensional ($w$) e do parâmetro de impacto ($y$).
3. Estudo de Injeção: Os sinais lentiados são injetados em ruído gaussiano gerado a partir das densidades espectrais de potência (PSDs) de uma rede de três detectores (LIGO Hanford, LIGO Livingston e Virgo) do período O4.
4. Estimação de Parâmetros: Os sinais injetados são recuperados usando inferência bayesiana (via o pacote Bilby e dynesty de amostragem aninhada/nested sampling) com templates de forma de onda que incluem tanto parâmetros da fonte (15 parâmetros binários usando o aproximante IMRPhenomXPHM) quanto parâmetros da lente ($M_{SIS}$ da lente SIS e parâmetro de impacto $y$).
5. Identificação de Candidatos: A partir das distribuições posteriores recuperadas, uma região de localização no céu com 90% de credibilidade é construída. Esta região é cruzada com um catálogo de AGs conhecidos da Via Láctea (de [37]) para identificar o possível lente.
6. Recuperação da Dispersão de Velocidade: Uma vez identificado um candidato a AG, sua distância conhecida ($D_l$) do catálogo é usada para inverter a posterior de $M_{SIS}$ recuperada para estimar a dispersão de velocidade ($\sigma_v$) do aglomerado usando a relação $M_{SIS} \propto \sigma_v^4 D_l$.

Três AGs específicos do catálogo da Via Láctea foram selecionados para o estudo baseando-se em suas distintas dispersões de velocidade centrais: Liller-1 ($\sigma_v \approx 19,6$ km s$^{-1}$), NGC 6093 ($\sigma_v \approx 11,1$ km s$^{-1}$) e NGC 2419 ($\sigma_v \approx 5,1$ km s$^{-1}$). Para cada aglomerado, 20 injeções foram realizadas com variando os parâmetros de impacto ($y \in [0,001, 5]$).

Resultos Principais

• Recuperação da Massa de Lente e Dispersão de Velocidade: Para alinhamentos fonte-lente favoráveis (pequenos parâmetros de impacto), os valores injetados de dispersão de velocidade foram recuperados com sucesso dentro dos intervalos de credibilidade de 68%. Para Liller-1, o valor mediano recuperado foi de 19,7 km s$^{-1}$, correspondendo de perto ao valor injetado de 19,6 km s$^{-1}$.
• Impacto do Parâmetro de Impacto: A precisão da recuperação é fortemente dependente do parâmetro de impacto $y$. À medida que $y$ aumenta, a amplificação diminui, enfraquecendo a impressão da lenteção no sinal. Isso leva a uma amostragem inferior da massa de lenteção e ao aumento da incerteza na dispersão de velocidade inferida. Para $y$ maiores, os valores medianos recuperados podem se deslocar do valor real devido a realizações de ruído.
• Localização no Céu e Cruzamento de Catálogo: O aglomerado injetado foi corretamente identificado dentro da região de localização no céu de 90% de credibilidade em todos os casos testados. O método distinguiu com sucesso o aglomerado de lenteção de outros no catálogo, desde que a fonte estivesse posicionada atrás do aglomerado alvo específico.
• Degenerescência de Parâmetros: O estudo observou uma anti-correlação entre a massa de lenteção recuperada ($M_{SIS}$) e o parâmetro de impacto ($y$) para $y \sim 1$, consistente com achados anteriores para lentes de massa pontual.

Significância e Alegações
O artigo demonstra que GWs lentiadas podem servir como uma sonda complementar para a dinâmica de AGs, permitindo especificamente a estimativa da dispersão de velocidade do aglomerado e, por extensão, da massa. Isso oferece um novo método para superar as limitações inerentes à determinação de massa baseada em EM, particularmente em regiões densas ou de núcleo colapsado.

Os autores observam que a probabilidade de um AG da Via Láctea atuar como uma lente para um evento de GW é baixa ($\tau \sim 10^{-8}$), destacando o potencial de detectores de próxima geração para descobrir tais eventos. A metodologia é apresentada como uma ferramenta para:

1. Identificar confiavelmente eventos de GW como sinais lentiados se um AG da Via Láctea for a lente.
2. Restringir a dispersão de velocidade e a massa do aglomerado de lenteção.
3. Potencialmente indicar a presença de IMBHs nos núcleos dos aglomerados.

O estudo é de escopo modesto, utilizando o modelo SIS e focando em AGs da Via Láctea onde as distâncias são bem conhecidas. Os autores reconhecem que, para AGs extragalácticos (onde os catálogos podem ser incompletos), existiria uma degenerescência entre a distância da lente e a dispersão de velocidade, exigindo uma estimação de parâmetros mais extensa. Eles também observam que, embora o trabalho atual utilize o modelo SIS, a metodologia pode ser estendida para modelos mais realistas (ex: esferas isotérmicas com núcleo, modelos de King) em trabalhos futuros. O artigo conclui que uma rede global de pelo menos três detectores de GW com sensibilidade comparável é essencial para aumentar a precisão da localização no céu necessária para que esta técnica seja eficaz.