A Parameterized YSO Accretion Disk Model with Increasing Accretion Rate: Predicted Outburst Lightcurves

Resumo Técnico: Um Modelo de Disco de Acreção de Objeto Estelar Jovem Parametrizado com Taxa de Acreção Crescente

Declaração do Problema
Objetos Estelares Jovens (YSOs), particularmente a subclasse FU Orionis (FU Ori), sofrem explosões de acreção dramáticas caracterizadas por aumentos de luminosidade de várias ordens de magnitude ao longo de meses a anos, seguidos por decaimentos ao longo de décadas. Embora a transição da acreção magnetosférica de estado baixo para a acreção de explosão de estado alto seja um fenômeno reconhecido, compreender a evolução específica dos componentes do sistema — fotosfera estelar, choques magnetosféricos, discos de gás viscosos e discos de poeira passivos — durante esses eventos permanece desafiador. Modelos existentes frequentemente dependem de simulações numéricas complexas e computacionalmente custosas (por exemplo, MHD) ou de suposições de estado estacionário que podem não capturar as mudanças morfológicas dependentes do tempo observadas em curvas de luz multibanda. Há uma necessidade de uma estrutura parametrizada que acople esses componentes físicos para simular a evolução temporal das distribuições de energia espectral (SEDs) e das curvas de luz em resposta a taxas de acreção variáveis.

Metodologia
Os autores apresentam um pipeline de modelagem parametrizado construído sobre o código YSOpy, projetado para simular fotometria e espectros de séries temporais para estrelas jovens em acreção. Diferentemente de simulações numéricas de transferência radiativa completa ou autoconsistentes, esta abordagem utiliza uma infraestrutura modular para calcular espectros para quatro componentes distintos do sistema separadamente antes de somá-los:

1. Fotosfera Estelar: Modelada como um corpo negro ou camada atmosférica.
2. Choques de Acreção Magnetosférica: Modelados usando uma aproximação de camada de hidrogênio onde a temperatura ($T_{hotspot}$), a densidade eletrônica ($n_e$) e a profundidade óptica ($\tau$) são derivadas de forma autoconsistente a partir da taxa de acreção ($\dot{M}$) e da velocidade de queda livre ($v_{ff}$).
3. Disco de Gás Aquecido Viscamente: Modelado usando um perfil de temperatura modificado de Shakura & Sunyaev (1973) para o disco interno, utilizando grades atmosféricas BT-Settl para temperaturas $>1400$ K e corpos negros para regiões mais frias.
4. Disco de Poeira Aquecido Passivamente: Modelado usando um perfil de temperatura de Chiang & Goldreich (1997) para poeira irradiada, com um limite interno determinado pelo raio de sublimação da poeira ($R_{sub}$).

O pipeline itera através de um perfil de taxa de acreção dependente do tempo ($\dot{M}(t)$). Em cada passo de tempo, o modelo calcula dinamicamente raios críticos, incluindo o raio de truncamento magnetosférico ($R_{trunc}$), o raio de corotação ($R_{co-rot}$) e o raio do disco interno de poeira ($R_{dust,inner}$). O raio do disco interno de gás ($R_{in}$) é determinado pelo equilíbrio entre pressões magnéticas e de acreção, transitando de $R_{co-rot}$ em baixos $\dot{M}$ para $R_{trunc}$ (ou $R_*$) em altos $\dot{M}$. O modelo leva em conta o "desligamento" da acreção magnetosférica poloidal quando o disco atinge a superfície estelar ($R_{in} \approx R_*$). Finalmente, os espectros de alta resolução resolvidos no tempo são convoluídos com funções de transmissão de filtros padrão (por exemplo, Gaia, 2MASS, WISE) para gerar curvas de luz e curvas de cor multibanda.

Principais Contribuições

• Estrutura Parametrizada Integrada: O artigo introduz um pipeline simplificado que acopla componentes estelares, de choque, de gás e de poeira em um único modelo que evolui no tempo, permitindo a simulação de ciclos completos de explosão sem o custo computacional de simulações MHD 3D.
• Física de Choque Autoconsistente: O modelo deriva temperaturas de manchas quentes e densidades eletrônicas diretamente da taxa de acreção e da velocidade de queda livre, garantindo conservação de energia ($L_{slab} = L_{blackbody}$) via um fator de escala $\chi$.
• Limites de Disco Dinâmicos: O modelo rastreia explicitamente a evolução do limite do disco interno de gás e da frente de sublimação da poeira como funções de $\dot{M}$, capturando a transição de um disco limitado pela corotação para um disco truncado magneticamente e, finalmente, para um disco tocando a superfície estelar.
• Ferramentas de Diagnóstico Multibanda: A infraestrutura gera curvas de luz e curvas de cor sintéticas em bandas ópticas e infravermelhas, fornecendo uma ferramenta de diagnóstico para interpretar morfologias de explosões observadas.

Resultados
Os autores aplicaram o modelo a uma taxa de acreção linearmente crescente e a quatro perfis de acreção específicos imitando objetos FU Ori observados: HBC 722 (linear de três partes), V890 Aur (tangente hiperbólica), V960 Mon (exponencial de dois lados) e Gaia 17bpi (logístico $\times$ Gaussiano). As principais descobertas incluem:

• Evolução da Dominância de Componentes:

  • Baixo $\dot{M}$: O sistema é dominado pela fotosfera estelar e pelo disco de poeira passivo; choques magnetosféricos são negligenciáveis ou mais frios que a fotosfera.
  • $\dot{M}$ Crescente: À medida que $\dot{M}$ aumenta, as manchas quentes magnetosféricas brilham rapidamente, dominando primeiro as bandas UV e ópticas azuis (por exemplo, banda u).
  • Pico/Alto $\dot{M}$: O aquecimento viscoso do disco interno de gás torna-se a fonte dominante de fluxo, particularmente no infravermelho próximo (bandas Gaia, 2MASS). Se $\dot{M}$ exceder $\sim 10^{-5} - 10^{-4} M_\odot$ yr$^{-1}$, o raio de truncamento do disco move-se para dentro até a superfície estelar, desligando a acreção poloidal e a emissão de manchas quentes.
  • Comportamento no Infravermelho Médio: Bandas no infravermelho médio (WISE) respondem principalmente ao aquecimento e expansão da frente de sublimação da poeira, frequentemente brilhando antes das bandas ópticas à medida que o raio do disco interno de poeira se move para fora.

• Morfologia da Curva de Luz:

  • Curvas de luz ópticas vermelhas e no infravermelho próximo geralmente acompanham o perfil de taxa de acreção de entrada.
  • Curvas de luz no infravermelho médio são mais sensíveis à geometria e ao aquecimento do disco interno de poeira, às vezes exibindo um brilho "pré-cursor" ou temporização distinta em relação às bandas ópticas.
  • A evolução da cor varia: sistemas frequentemente tornam-se "mais azuis quando mais brilhantes" durante a subida devido ao aquecimento de manchas quentes e do disco interno, mas podem avermelhar durante o decaimento à medida que o disco viscoso esfria e a fotosfera reaparece.

• Sensibilidade aos Parâmetros:

  • Força do Campo Magnético: Campos mais fortes atrasam a aproximação do disco à superfície estelar, prolongando a fase de acreção magnetosférica e criando curvas de luz com múltiplos picos no UV/óptico.
  • Período de Rotação: Variações no período de corotação afetam principalmente o raio do disco interno no estado baixo, mas têm impacto negligenciável nas curvas de luz durante a fase de explosão de alta acreção.

Significado e Alegações
Os autores posicionam este trabalho como uma ferramenta ilustrativa para interpretar os mecanismos físicos por trás das explosões de FU Ori, em vez de um modelo preditivo para objetos individuais específicos. Eles afirmam que sua abordagem parametrizada demonstra com sucesso como diferentes componentes do sistema dominam o fluxo em vários comprimentos de onda ao longo de um ciclo de explosão.

O artigo destaca que o temporamento do brilho em diferentes comprimentos de onda oferece poder de diagnóstico: o brilho UV/óptico azul indica o início do aquecimento de choque magnetosférico, enquanto o brilho no infravermelho médio sinaliza a expansão da frente de sublimação da poeira. Os autores sugerem que futuras instalações como LSST e UVEX poderiam utilizar essas assinaturas previstas para identificar explosões em estágios iniciais. Além disso, o modelo fornece uma estrutura para inferir perfis de acreção subjacentes e taxas de acreção de estado baixo a partir de curvas de luz e curvas de cor observadas, preenchendo a lacuna entre descrições fenomenológicas simples e simulações físicas complexas. Os autores reconhecem que, embora seu modelo simplifique a reestruturação do campo magnético e as camadas limite de acreção equatorial, ele captura a física essencial da transição de estados de acreção semelhantes a T-Tauri para estados semelhantes a FU Ori.