The dispersal of compact protoplanetary discs

Resumo Técnico: A Dispersão de Discos Protoplanetários Compactos

Declaração do Problema
Discos protoplanetários são os locais de formação planetária, contudo seus mecanismos de dispersão permanecem um tema de investigação ativa. Evidências observacionais, particularmente de estudos de evolução cor-cor em populações estelares jovens, indicam que a dispersão do disco procede preferencialmente de dentro para fora, criando uma cavidade central antes que o disco externo seja erodido. No entanto, modelos teóricos padrão frequentemente assumem discos com grandes raios externos (centenas de unidades astronômicas), muito além do raio gravitacional ($R_g$) onde a fotoevaporação é mais eficiente. Levantamentos recentes (e.g., ALMA) revelam que uma fração substancial de discos é "compacta", com raios externos comparáveis ou menores que $R_g$. Modelos existentes, que tipicamente assumem perfis de perda de massa não restritos estendendo-se a grandes raios, falham em reproduzir com precisão as trajetórias evolutivas e as escalas de tempo de dispersão desses sistemas compactos. Especificamente, prescrições padrão frequentemente preveem erosão de fora para dentro para discos compactos, contradizendo diagnósticos observacionais.

Metodologia
Os autores empregaram uma abordagem multietapa combinando simulações hidrodinâmicas de radiação de alta resolução com modelos de síntese populacional viscosa 1D:

1. Simulações Hidrodinâmicas de Radiação: Utilizando o código PLUTO acoplado ao framework de transferência radiativa MOCASSIN, os autores simularam a fotoevaporação impulsionada por raios X para uma estrela de $1 M_\odot$. Eles exploraram quatro raios de corte iniciais distintos do disco ($R_{cut}$): 10, 50, 100 e 200 ua. As simulações rastrearam a estrutura térmica, densidade e campos de velocidade até que um estado quase estacionário fosse alcançado. O raio gravitacional foi adotado como $R_g = 5$ ua.
2. Derivação de Relações de Escala: A partir das simulações hidrodinâmicas, os autores analisaram as distribuições cumulativas da taxa de perda de massa. Eles determinaram como a taxa total de perda de massa do vento depende do raio de corte externo do disco.
3. Síntese Populacional: Essas descobertas foram incorporadas a um código de evolução viscosa 1D (SPOCK). Os autores realizaram uma síntese populacional de 4.096 discos sintéticos usando uma estratégia de amostragem Quasi-Monte Carlo (sequências de Sobol) para explorar um espaço de parâmetros de 9 dimensões (massa estelar, luminosidade em raios X, massa inicial do disco, raio de corte, viscosidade, campo externo de FUV, etc.).
4. Prescrições Testadas: Os modelos compararam dois cenários:
   • Prescrição Padrão: Taxas de perda de massa por fotoevaporação interna aplicadas sem considerar o raio externo finito do disco (efetivamente assumindo que o disco se estende ao infinito ou além do pico de perda de massa).
   • Prescrição Dependente do Corte: A taxa de perda de massa é truncada e normalizada com base no raio de corte específico do disco, derivada das relações de escala hidrodinâmicas.
   • Fotoevaporação Externa: Alguns modelos incluíram fotoevaporação externa suave por FUV (usando a grade FRIED) para simular estrelas OB próximas.

Principais Contribuições e Resultados

• Escala Hidrodinâmica: As simulações revelaram que os perfis de perda de massa superficial são quase independentes do tamanho do disco; a morfologia do vento e as taxas locais de perda de massa nas regiões internas permanecem consistentes independentemente do corte externo. No entanto, a taxa total de perda de massa integrada é significativamente reduzida para discos compactos porque o vento é truncado no raio de corte. A taxa cumulativa de perda de massa escala com o raio, crescendo como $\dot{M}_w \propto R^3$ próximo ao raio gravitacional e achatando-se para $\dot{M}_w \propto R^{0.5}$ em raios maiores.
• Trajetórias Evolutivas:
  • Prescrição Padrão: Quando aplicada a discos compactos ($R_{cut} < 20-30$ ua), o modelo padrão prevê que a taxa de perda de massa do vento excede a taxa de espalhamento viscoso fora do raio de corte. Isso leva à dispersão de fora para dentro, onde as bordas externas são erodidas primeiro, contradizendo evidências observacionais.
  • Prescrição Dependente do Corte: Quando a taxa de perda de massa é escalada para o raio de corte, a taxa total do vento é reduzida suficientemente para permitir que o espalhamento viscoso domine inicialmente. Isso resulta em dispersão de dentro para fora, onde uma lacuna se abre próximo ao raio gravitacional e se expande para fora, consistente com diagnósticos observacionais.
• Resultados da Síntese Populacional:
  • A prescrição dependente do corte, combinada com uma taxa de fotoevaporação interna reduzida (escalada por um fator de 10 para alinhar com modelos químicos recentes), reproduz as frações de discos observadas e as taxas de acreção em função da idade do aglomerado sem ajuste de parâmetros.
  • A vida média do disco no modelo aprimorado é de $2,70 \pm 0,13$ Myr, o que se alinha bem com estimativas observacionais ($2,24 \pm 0,35$ Myr) e é significativamente mais longa do que as previsões de modelos sem a redução da taxa ($1,13 \pm 0,05$ Myr).
  • A inclusão de fotoevaporação externa suave não impede o espalhamento do disco quando a prescrição de raio de corte é utilizada, mas resulta em um aumento mais suave nos raios do disco ao longo do tempo e em uma fração maior de discos compactos persistindo até idades mais avançadas.

Significância
O artigo argumenta que a dispersão de discos protoplanetários compactos não pode ser capturada com precisão por modelos padrão de fotoevaporação não escalados. Os autores propõem uma relação de escala simples e motivada fisicamente que leva em conta o raio externo finito do disco. Ao incorporar essa dependência, os modelos reconciliam com sucesso as previsões teóricas com o esvaziamento observado de dentro para fora de populações estelares jovens e a prevalência de discos compactos. O estudo conclui que as prescrições de dispersão devem incluir explicitamente a dependência do raio de corte do disco para prever corretamente a evolução e a vida útil de discos compactos, fornecendo uma correção necessária aos modelos atuais de síntese populacional.