Accretion of Primordial Black Holes in Stellar Interiors

Resumo Técnico: Acreção de Buracos Negros Primordiais em Interiores Estelares

Enunciado do Problema
Buracos Negros Primordiais (PBHs) permanecem como um candidato viável à matéria escura, particularmente na faixa de massa sub-lunar mal restringida ($10^{-16} M_\odot \lesssim M_{BH} \lesssim 10^{-10} M_\odot$). As estrelas fornecem um laboratório único para sondar PBHs neste regime, uma vez que um PBH capturado e incorporado ao núcleo estelar ocorreria acreção do gás circundante, potencialmente alterando a evolução da estrela. Estudos anteriores (ex: Bellinger et al. 2023) modelaram este cenário utilizando códigos de evolução estelar 1D, mas basearam-se em suposições derivadas da teoria de disco de acreção fino. Especificamente, assumiram uma eficiência radiativa fixa $\eta \approx 0,08$ e aplicaram o limite de luminosidade de Eddington para todas as massas. No entanto, para PBHs nesta faixa de massa, o raio de Bondi ($r_B$) é ordens de magnitude menor que o raio de circularização estelar, impedindo a formação de um disco. Em vez disso, a acreção é esférica, e o gás é comprimido e aquecido a temperaturas relativísticas ($T \sim 10^{10}–10^{12}$ K) próximo ao raio de Schwarzschild ($r_S$). Neste regime, a microfísica (bremsstrahlung relativística, processos de par, espalhamento Klein–Nishina) e a profundidade óptica do fluxo diferem fundamentalmente do paradigma do disco fino, necessitando de um cálculo autoconsistente da taxa de acreção ($\dot{M}$), eficiência radiativa ($\eta$) e estrutura do fluxo.

Metodologia
Os autores resolvem as equações de Euler dependentes do tempo com um termo de fonte de resfriamento implícito para modelar o fluxo de acreção sobre um PBH incorporado em um núcleo estelar do tipo solar.

• Microfísica: O termo de resfriamento $\epsilon(\rho, T)$ inclui bremsstrahlung relativística (elétron-íon e elétron-elétron) e aniquilação elétron-pósitron. A opacidade contabiliza o espalhamento Thomson, absorção free-free e supressão Klein–Nishina em altas energias fotônicas.
• Abordagem Numérica: Os autores utilizam o radbondi, um pacote Python de código aberto, para evoluir o fluxo de um perfil de Bondi adiabático inicial para um estado estacionário autoconsistente. Esta abordagem dependente do tempo é necessária porque, no regime dominado pelo resfriamento, o ponto sônico torna-se um foco (autovalores complexos) em vez de um ponto de sela, tornando os métodos padrão de shooting de ODE em estado estacionário inválidos.
• Tratamento de Regimes:
  • Opticamente Fino ($M_{BH} \lesssim 5 \times 10^{-13} M_\odot$): O solver de Euler computa diretamente a eficiência radiativa física e a taxa de acreção.
  • Opticamente Espesso ($M_{BH} \gtrsim 5 \times 10^{-13} M_\odot$): A esfera de Bondi torna-se opticamente espessa, entrando em um regime de aprisionamento de fótons. Os autores adotam a prescrição analítica de Begelman (1979), onde fótons produzidos localmente no interior do raio de aprisionamento ($r_{trap}$) são advectados de volta ao buraco negro, enquanto a luminosidade que escapa é autorregulada em $L_{esc} \lesssim 0,6 L_{Edd}$.
• Feedback: O estudo modela autoconsistentemente o feedback radiativo, onde a luminosidade depositada fora da esfera de Bondi aquece o plasma ambiente. Este aquecimento é moderado pelo transporte convectivo (Teoria do Comprimento de Mistura), evitando os aumentos extremos de temperatura previstos pela difusão radiativa pura.

Resultos Principais
O estudo identifica três regimes distintos de acreção esférica e deriva os seguintes resultados:

1. Três Regimes de Acreção:

   • Bondi Quente ($M_{BH} \lesssim 10^{-14} M_\odot$): O resfriamento é dinamicamente desprezível. O fluxo é quase adiabático, atingindo temperaturas $T \sim 10^{11}$ K próximo a $r_S$. O fluxo interno torna-se colisional antes de atingir temperaturas de produção de pares, suprimindo o equilíbrio de pares. O bremsstrahlung domina o resfriamento.
   • Resfriamento por Bremsstrahlung ($10^{-14} \lesssim M_{BH} \lesssim 5 \times 10^{-13} M_\odot$): O resfriamento torna-se eficiente, drenando a energia térmica antes que o gás atinja temperaturas de produção de pares. O fluxo aproxima-se de um estado isotérmico ($T \approx T_\infty$), e a taxa de acreção é aumentada por um fator de $\sim 2–7$ em relação à taxa de Bondi adiabática.
   • Aprisionamento de Fótons ($M_{BH} \gtrsim 5 \times 10^{-13} M_\odot$): A esfera de Bondi é opticamente espessa. Embora a luminosidade emitida localmente aumente rapidamente com a massa, a luminosidade que escapa é limitada pelo limite de Begelman ($L_{esc} \lesssim 0,6 L_{Edd}$). Crucialmente, a taxa de acreção permanece próxima ao valor de Bondi ($\dot{M} \propto M_{BH}^2$) porque não há um limite de Eddington esférico para o influxo; a pressão de radiação apenas modula a gravidade efetiva.

2. Eficiência Radiativa: A eficiência radiativa computada autoconsistentemente é $\eta \approx 10^{-2}$ em toda a faixa opticamente fina, um valor um fator de ordem inferior à canônica eficiência de disco fino ($\eta \approx 0,08$) assumida na literatura anterior. A eficiência atinge o pico próximo a $10^{-15} M_\odot$ mas cai nos limites de resfriamento e isotérmico.

3. Trajetórias de Crescimento: Devido ao fato de o limite de Eddington esférico não se aplicar, o crescimento do PBH permanece super-exponencial ($\dot{M} \propto M_{BH}^2$) durante todo o regime de aprisionamento de fótons. O gargalo de crescimento é determinado inteiramente pela lenta fase de Bondi Quente nas massas mais baixas.

4. Massa Inicial Crítica: A massa inicial crítica necessária para que um PBH consuma uma estrela de massa solar dentro de um tempo de Hubble é $M_{0,crit} \sim 10^{-16} M_\odot$. Isto é cinco ordens de magnitude menor do que as estimativas anteriores ($\sim 10^{-11} M_\odot$) que assumiam maior eficiência e um gargalo de Eddington.

5. Campos Magnéticos: Os autores descobrem que o acúmulo de fluxo magnético é desprezível. O campo magnético ambiente no núcleo solar é fraco demais para deter dinamicamente o fluxo, e a reconexão magnética limita qualquer acúmulo de fluxo a níveis muito abaixo da importância dinâmica ($P_{mag}/P_{ram} \ll 1$).

Significância e Alegações
O artigo afirma fornecer o primeiro cálculo autoconsistente de eficiência radiativa e taxas de acreção para a acreção esférica de PBHs através da faixa de massa sub-lunar. Ao substituir as suposições de disco fino pelo tratamento microfísico de fluxo esférico, os autores demonstram que:

• A eficiência radiativa é significativamente menor do que o assumido anteriormente.
• Não existe uma fase limitada por Eddington para a acreção esférica; o fluxo permanece super-Eddington em termos de fluxo de massa sem ser limitado pela pressão de radiação.
• Consequentemente, as restrições sobre a matéria escura de PBH a partir de observações estelares são muito mais rigorosas do que se pensava anteriormente, pois PBHs com massas iniciais tão baixas quanto $10^{-16} M_\odot$ podem crescer para consumir suas estrelas hospedeiras dentro de um tempo de Hubble.

Os autores observam que seu trabalho foca na microfísica da acreção e não inclui a reação reversa na estrutura estelar global (o que exigiria o acoplamento a um código de evolução estelar como o MESA), mas eles fornecem as relações necessárias de $\eta(M_{BH})$ e $\dot{M}(M_{BH})$ para permitir tais estudos futuros. Eles também deixam a discussão sobre a formação de discos (ocorrendo em massas mais altas ou em estrelas com rotação mais rápida) para um artigo complementar.