Disk-Regulated Mass Transfer Between Rotating Non-Degenerate Stars: Insights from Be and sdOB Binaries

Resumo Técnico: Transferência de Massa Regulada por Disco entre Estrelas Não Degeneradas em Rotação

Definição do Problema
A transferência de massa entre estrelas não degeneradas é um processo crítico, porém ainda pouco compreendido, na evolução binária. Uma incerteza central é a eficiência de acreção — a fração da massa transferida que é retida pelo acretor. As prescrições padrão de síntese de populações binárias (BPS) e simulações detalhadas de evolução estelar (por exemplo, usando o MESA) tipicamente empregam uma prescrição de acreção "rotacionalmente limitada". Neste arcabouço, a acreção é restringida assim que o acretor se aproxima de sua taxa de rotação crítica de ruptura, o que frequentemente leva a eficiências de acreção muito baixas.

No entanto, restrições observacionais de sistemas binários de estrelas Be, particularmente sistemas Be + subdwarf O/B-type (sdOB), sugerem eficiências de acreção significativamente mais altas (frequentemente $\gtrsim 30\%$ a $>50\%$). Os modelos rotacionalmente limitados existentes subestimam sistematicamente as massas das estrelas Be nesses sistemas e falham em reproduzir o espaço de parâmetros observado de binárias Be+sdOB. Além disso, estudos anteriores indicam que nem as prescrições rotacionalmente limitadas nem as termicamente limitadas, isoladamente, podem satisfazer os requisitos de eficiência derivados das populações observadas de Be+sdOB.

Metodologia
Os autores propõem uma prescrição de acreção de massa de "acoplamento disco-estrela" motivada fisicamente para abordar essas discrepâncias. O mecanismo físico central postula que um disco de acreção ou de decréscimo regula o momento angular transportado para o acretor, permitindo que a acreção de massa continue mesmo quando a estrela rotaciona próxima de seu limite crítico.

1. Regulação do Momento Angular: Baseando-se nas soluções de disco fino em estado estacionário de Popham & Narayan (1991), os autores modelam o momento angular específico acretado pela estrela ($j_{acc}$) como uma função da taxa de rotação estelar ($\Omega/\Omega_{crit}$).

   • Em taxas de rotação baixas, $j_{acc}$ é aproximadamente o valor Kepleriano na superfície estelar.
   • À medida que a estrela se aproxima de uma taxa de rotação limiar (adotada como $\Omega/\Omega_{crit} \approx 0,9$), a tensão de cisalhamento do disco permite que o momento angular seja transportado para fora da estrela para o disco. Consequentemente, $j_{acc}$ diminui bruscamente, podendo tornar-se negativo, permitindo a continuação da acreção de massa sem que a estrela exceda sua rotação crítica.
   • Os autores implementam uma forma funcional suavizada para $j_{acc}$ para garantir a estabilidade numérica em simulações 1D, evitando descontinuidades que surgiriam de um limiar estrito.

2. Leis de Conservação e Cálculo de Eficiência: O modelo assume que o excesso de momento angular extraído da estrela é transportado através do disco e removido via perda de massa na borda externa do disco ($R_{out}$). Ao aplicar a conservação de massa e de momento angular, a eficiência de acreção ($\eta = \dot{M}_{acc}/\dot{M}_{tr}$) é calculada de forma autoconsistente. A eficiência depende da razão entre o momento angular específico na borda do disco e na borda de acreção.

3. Estrutura de Simulação: Os autores implementaram esta prescrição dentro do framework POSYDON, utilizando o MESA (Modules for Experiments in Stellar Astrophysics) para evoluir modelos binários.

   • Eles desativaram o esquema de vento implícito usado para prevenir a rotação supercrítica, dependendo, em vez disso, de ventos rotacionalmente aumentados explícitos e da taxa de acreção reduzida no regime supercrítico para restaurar a rotação subcrítica.
   • Uma grade de modelos binários foi simulada com metalicidade solar, variando as massas primárias ($2-9 M_\odot$), razões de massa ($0,4-0,95$) e períodos orbitais.
   • O estudo também investigou o impacto dos tratamentos de overshooting do núcleo (comparando as configurações padrão do POSYDON com uma prescrição de overshooting reduzido) sobre as massas finais dos componentes.

Principais Contribuições e Resultados

• Estimativas Analíticas: Os autores derivaram expressões analíticas para a eficiência de acreção. Eles mostram que, no regime de acreção por disco, a eficiência depende da rotação estelar, da razão de massa e da razão entre o raio estelar e o raio do lóbulo de Roche. No regime de impacto direto, ela depende da rotação e da razão de massa. As eficiências de acreção ponderadas pela massa global são encontradas próximas aos valores esperados perto da taxa de rotação limiar, onde o momento angular específico acretado declina bruscamente.
• Comparação de Evolução Binária:
  • Modelo Rotacionalmente Limitado: Neste modelo padrão, os acretores atingem a rotação crítica após acrecer apenas uma pequena fração da massa transferida ($\sim 0,2 M_\odot$ em um caso de teste), levando a eficiências baixas ($<10\%$ para a maioria dos sistemas). Isso resulta em estrelas Be sistematicamente menos massivas do que o observado.
  • Modelo de Acoplamento Disco-Estrela: Este modelo permite que os acretores acretem significativamente mais massa (por exemplo, $\sim 1,6 M_\odot$ na fase rápida do caso de teste) com eficiências atingindo $\sim 40-60\%$. O disco extrai o momento angular, permitindo que a estrela se expanda e acreça enquanto permanece próxima do limiar de rotação crítica.
• Reprodução de Observações:
  • Quando aplicado a sistemas Be+sdOB, o modelo rotacionalmente limitado falha em produzir estrelas Be com massas consistentes com as observações, particularmente para sistemas de órbita larga (Caso B).
    • O modelo de acoplamento disco-estrela produz estrelas Be com massas que cobrem uma parte significativa do intervalo observado.
  • Os autores descobriram que as massas finais dos componentes são sensíveis ao overshooting assumido. Um modelo combinando a prescrição de acoplamento disco-estrela com overshooting reduzido produz massas de componentes para sistemas Be+sdOB que estão em maior concordância com as observações do que os modelos padrão de alto overshooting.
• Períodos Orbitais: Embora o modelo de acoplamento disco-estrela aumente a retenção de massa, os períodos orbitais resultantes são comparáveis ou ligeiramente mais longos do que os do modelo rotacionalmente limitado. Isso ocorre porque a transferência de massa mais conservativa reduz a massa total perdida do sistema, levando a razões de massa mais extremas que aceleram a expansão orbital, o que compensa parcialmente a perda de momento angular pela borda do disco.

Significância e Alegações
O artigo afirma que a discrepância entre as propriedades observadas de binárias Be+sdOB e os resultados atuais de simulação pode ser resolvida adotando um modelo de transferência de massa regulado por disco. Os autores argumentam que:

1. As prescrições de acreção rotacionalmente limitadas subestimam sistematicamente as massas das estrelas Be em sistemas Be+sdOB.
2. Um mecanismo de acoplamento disco-estrela, que permite a regulação do transporte de momento angular, possibilita eficiências de acreção mais altas, consistentes com as restrições observacionais (especificamente o requisito de eficiências $>50\%$ em muitos sistemas).
3. A combinação desta nova prescrição de acreção com o overshooting do núcleo reduzido fornece uma explicação mais fisicamente consistente para as massas e períodos orbitais observados de binárias Be+sdOB.

Os autores reconhecem que a eficiência real depende de fatores incertos, como a estrutura do disco, campos magnéticos e as condições específicas de formação de disco e perda de massa (por exemplo, via pontos L2 ou hot spots). No entanto, concluem que o modelo de acoplamento disco-estrela oferece um mecanismo físico necessário para explicar a existência de estrelas Be massivas em binárias de órbita larga, onde o frenagem (spin-down) por maré é ineficiente.