Drifters on the edge of town: $λ$ Boötis stars in clusters

Resumo Técnico: Errantes na periferia da cidade: estrelas λ Boötis em aglomerados

Definição do Problema
As estrelas λ Boötis (λ Boo) são estrelas do tipo A quimicamente peculiares, caracterizadas por abundâncias de elementos leves (C, N, O, S) semelhantes às do Sol, mas com deficiências significativas de elementos do pico do ferro. A principal explicação teórica postula que estas estrelas accrecionam gás prístino e pobre em metais do Meio Interestelar (ISM). No entanto, este cenário enfrenta um desafio significativo: descobertas recentes de estrelas λ Boo dentro ou nas periferias de aglomerados estelares contradizem a expectativa de que a intensa radiação ionizante de estrelas intermediárias e massivas (> 5 M⊙) em tais ambientes impediria a acreção de gás prístino do ISM ou removeria os discos circunestelares necessários para o processo. Além disso, existe uma tensão entre as taxas de acreção necessárias para produzir abundâncias de λ Boo (10⁻¹⁴ a 10⁻⁹ M⊙ yr⁻¹) e as taxas muito menores inferidas de observações de anãs brancas (< 10⁻¹⁷ M⊙ yr⁻¹). O problema específico abordado é se as histórias dinâmicas de estrelas A em aglomerados em evolução permitem que elas escapem do núcleo irradiante do aglomerado, acrecionem gás prístino e, posteriormente, retornem, reconciliando assim a presença de estrelas λ Boo em ambientes de aglomerados com a hipótese de acreção do ISM.

Metodologia
Os autores utilizam simulações de N-corpos puramente gravitacionais para modelar a formação e a evolução dinâmica de regiões de formação estelar ao longo de 100 Myr.

• Configuração da Simulação: Quatro conjuntos de simulações foram conduzidos com números totais de estrelas ($N_*$) de 150, 500, 1000 e 3000. As condições iniciais incluíram subestrutura espacial e cinemática gerada pelo método fractal de caixa (dimensão fractal $D=1.6$) e uma razão virial de $\alpha_{vir} = 0.3$, levando ao colapso inicial e à formação de aglomerados centralmente concentrados. As massas estelares foram amostradas de uma Função de Massa Inicial (IMF) de Maschberger (2013) com um limite superior de 50 M⊙. Binários primordiais foram excluídos.
• Análise de Pós-Processamento: As simulações foram analisadas post-hoc para determinar:
  1. Potencial de Acreção: As estrelas A foram rastreadas para identificar quando se moviam além do raio de Jacobi (tidal) ($r_J$), definido como o limite onde poderiam encontrar gás prístino do ISM. Quando fora de $r_J$, as taxas de acreção de Bondi-Hoyle-Lyttleton (BHL) foram calculadas com base na velocidade da estrela em relação a uma nuvem de ISM estacionária ($\rho = 10$ cm⁻³).
  2. Ambiente de Radiação: Os fluxos de UV distante (FUV) e UV extremo (EUV) incidentes em cada estrela A foram calculados com base na distância para estrelas massivas (> 5 M⊙ e > 20 M⊙, respectivamente).
  3. Sobrevivência do Disco: Utilizando a grade FRIED de modelos de fotoevaporação e prescrições analíticas, os autores modelaram a perda de massa dos discos circunestelares das estrelas A. Um disco foi considerado "destruído" se a massa total perdida excedesse a massa inicial do disco (assumida como 10% da massa estelar). Uma estrela era considerada candidata à formação de λ Boo apenas se acrescionasse material do ISM enquanto retivesse seu disco.
  4. Restrições Temporais: Os cálculos foram repetidos com uma restrição que impedia a acreção após 40 Myr para contabilizar o potencial enriquecimento do ISM por supernovas do Tipo Ia.

Principais Resultados

• Históricos Dinâmicos: Em todas as realizações das simulações, a maioria (90–98%) das estrelas A viajou além do raio de Jacobi durante a evolução de 100 Myr. Isso ocorreu tanto através da ejeção do núcleo do aglomerado quanto através de "deriva", onde as estrelas se moveram para fora do raio tidal e posteriormente retornaram às regiões internas.
• Taxas de Acreção: As estrelas A que se moveram além do raio de Jacobi atingiram taxas de acreção BHL variando de $10^{-12}$ a $10^{-5}$ M⊙ yr⁻¹. Estrelas ejetadas a altas velocidades apresentaram taxas menores ($10^{-12}–10^{-10}$ M⊙ yr⁻¹), enquanto aquelas em órbitas centradas no aglomerado que derivaram para fora e voltaram apresentaram taxas mais altas ($10^{-9}–10^{-5}$ M⊙ yr⁻¹), sobrepondo-se à faixa necessária para explicar as abundâncias de λ Boo.
• Sobrevivência do Disco vs. Radiação: A sobrevivência dos discos circunestelares dependeu fortemente do número total de estrelas (e, portanto, do número de fontes ionizantes) no aglomerado.
  • Em regiões de baixa massa ($N_* = 150$), ~55% das estrelas A que se aventuraram além do raio de Jacobi retiveram seus discos por > 10 Myr.
  • Em regiões de alta massa ($N_* = 3000$), 0% das estrelas A retiveram seus discos devido aos intensos campos de FUV/EUV, o que impede efetivamente a formação de λ Boo nestes ambientes.
  • Em regiões intermediárias ($N_* = 500, 1000$), uma fração das estrelas A (8,8% e 4,8%, respectivamente) reteve discos enquanto realizava a acreção.
• Sistemas Binários: O estudo descobriu que sistemas binários largos (> 100 au) contendo estrelas A podem se formar via captura durante a dissolução de regiões de formação estelar. Crucialmente, uma estrela poderia acrescionar material prístino do ISM, desenvolver abundâncias de λ Boo e, então, capturar um companheiro (ou vice-versa). Este mecanismo permite a existência de sistemas binários onde os componentes possuem abundâncias químicas significativamente diferentes. A distribuição de separação simulada de tais binários intercepta observados binários largos contendo estrelas λ Boo (ex: HD 87304 + CD-33 6615B).

Significância e Alegações
O artigo alega demonstrar que a presença de estrelas λ Boo em ou perto de aglomerados estelares não invalida inerentemente o cenário de acreção do ISM. Os autores argumentam que:

1. A Fuga Dinâmica é Suficiente: As estrelas A não precisam ser permanentemente ejetadas para acrescionar gás prístino; elas podem derivar além do raio tidal, acrescionar e retornar ao centro do aglomerado.
2. Limiares de Radiação: A viabilidade do cenário de acreção está fortemente correlacionada com a massa do aglomerado. Aglomerados de baixa massa (semelhantes a Theia 139 e HSC 1640, onde estrelas λ Boo foram encontradas recentemente) possuem campos de radiação suficientemente fracos para permitir a sobrevivência do disco e a acreção, enquanto aglomerados massivos não permitem.
3. Diagnóstico Binário: A observação de sistemas binários com abundâncias díspares entre os componentes não invalida a hipótese de acreção. Tais sistemas podem se formar dinamicamente após o evento de acreção ter ocorrido.

Os autores mantêm-se modestos quanto à validade absoluta das taxas de acreção, observando que a tensão entre as taxas requeridas e os limites de anãs brancas permanece não resolvida. Eles concluem que, embora seus modelos dinâmicos apoiem a possibilidade de formação de estrelas λ Boo em aglomerados, pesquisas adicionais são necessárias para validar as taxas de acreção específicas e a interação física entre discos em acreção e o campo de radiação.