Can photoevaporation open gaps in protoplanetary discs?

Enunciado do Problema
A dispersão de discos protoplanetários é um fator crítico na determinação do tempo e das condições da formação planetária. Embora a fotoevaporação — a remoção de gás via ventos termicamente impulsionados por radiação estelar de alta energia — seja amplamente reconhecida como um mecanismo fundamental de dispersão, a base teórica de como ela abre e sustenta lacunas permanece incerta. Modelos padrão tipicamente dependem de cálculos de evolução viscosa unidimensionais (1D), onde a taxa de perda de massa por fotoevaporação, $\dot{\Sigma}(r)$, é tratada como um termo de sumidouro fixo derivado de simulações bidimensionais (2D) estáticas. Essa abordagem assume implicitamente que o perfil de perda de massa permanece constante mesmo conforme a estrutura do disco evolui. Os autores argumentam que isso negligencia um ciclo de feedback crucial: à medida que uma lacuna começa a se formar, a estrutura local do disco muda, potencialmente alterando as condições de lançamento do vento, reduzindo as taxas locais de perda de massa e permitindo que o influxo viscoso preencha a região depletada. Se a fotoevaporação for autolimitante desta maneira, ela pode ser menos eficiente em limpar o gás nos raios onde planetas gigantes se formam do que o previamente assumido, com implicações significativas para modelos de acreção de núcleo e migração.

Metodologia
Para investigar se a fotoevaporação sozinha pode abrir e sustentar lacunas quando o fluxo responde dinamicamente à estrutura do disco, os autores empregaram uma abordagem computacional de duas fases:

1. Fase I (Evolução 1D): Eles utilizaram o código de evolução viscosa spock para simular um disco passando por fotoevaporação por raios-X até o início da formação de uma lacuna. A simulação evoluiu um perfil de densidade de superfície de gás ($\Sigma$) ao longo de 2,1 Myr, incorporando um perfil de perda de massa de Picogna et al. (2019). Esta fase estabeleceu as condições iniciais para as simulações 2D, especificamente um estado onde uma depleção parcial na densidade de superfície havia se desenvolvido.
2. Fase II (Hidrodinâmica de Radiação 2D): O perfil de densidade de superfície da simulação 1D em 2,1 Myr foi mapeado em um sistema de coordenadas esféricas axissimétrico 2D ($r-\theta$) usando o código pluto. Estas simulações modelaram o fluxo de saída fotoevaporativo e sua interação com o disco.
   • Modelo p-g: Utilizou a densidade de superfície evoluída da Fase I, que incluía uma depleção de densidade preexistente.
   • Modelo p-ng: Utilizou um perfil de densidade de superfície inicial modificado onde a depleção foi suprimida (substituída por uma lei de potência) para avaliar o impacto da morfologia inicial na abertura de lacunas.
   • As simulações empregaram uma aproximação de aquecimento de linha de visada e assumiram equilíbrio radiativo, desacoplando a transferência radiativa da hidrodinâmica para reduzir o custo computacional.

Principais Resultados

• Formação de Lacuna Autolimitante: Nas simulações 2D, uma vez que uma depressão de densidade se formou, a taxa local de perda de massa diminuiu bruscamente dentro da região da lacuna. Isso ocorreu porque a região depletada permaneceu na sombra do disco interno, reduzindo a exposição à irradiação estelar. Consequentemente, a lacuna não se aprofundou em uma cavidade totalmente evacuada.
• Reenchimento Parcial: O influxo viscoso e o transporte de massa radial ao longo da superfície do disco (e através do vento) atuaram para preencher parcialmente a região depletada. O sistema evoluiu para uma zona de depleção parcial persistente em vez de um buraco interno limpo.
• Insensibilidade às Condições Iniciais: Comparar modelos com (p-g) e sem (p-ng) uma depleção inicial mostrou que o resultado evolutivo de longo prazo foi amplamente independente da morfologia inicial do disco. Em ambos os casos, o sistema se autorregulou para um estado de depleção parcial.
• Limitações de Prescrições 1D: Os autores testaram uma prescrição de primeira ordem para aproximar o comportamento 2D em modelos 1D, escalonando o perfil primordial de perda de massa com base na depleção local da densidade de superfície. Embora isso tenha melhorado a concordância em relação às prescrições estáticas, os modelos 1D ainda evoluíram muito rapidamente e produziram lacunas mais profundas do que as simulações 2D. Esta discrepância destaca a incapacidade dos modelos 1D de capturar o transporte de massa radial ao longo da superfície do disco e a complexa interação de sombreamento.
• Assinaturas Observacionais: Embora a lacuna de gás permaneça parcialmente preenchida, o máximo de pressão na borda externa da depressão ainda pode atuar como uma armadilha eficiente para grãos de poeira. Isso poderia produzir assinaturas observacionais semelhantes a discos de transição (cavidades internas na emissão de contínuo) mesmo sem a limpeza completa do gás.

Significância e Alegações
O artigo desafia o paradigma padrão de que a fotoevaporação esculpe eficientemente cavidades internas limpas e produz diretamente discos de transição. Em vez disso, os autores propõem que a fotoevoração pode levar a um estado quase estacionário caracterizado por uma zona de depleção parcial persistente. Isso sugere que discos exibindo depleção moderada sem limpeza interna completa podem representar sistemas influenciados pela fotoevoração, mas não totalmente dispersos.

Os autores fornecem uma prescrição modesta de primeira ordem para aproximar este comportamento autorregulado em modelos de evolução de disco 1D, que é adequada para uso em estudos de formação planetária e síntese de populações. No entanto, eles declaram explicitamente que esta prescrição permanece aproximada e não pode capturar totalmente a natureza multidimensional do fluxo, particularmente a redistribuição de massa radial. O artigo conclui que novas simulações multidimensionais e exploração do espaço de parâmetros são necessárias para derivar receitas robustas para modelos globais de discos e populações planetárias. O trabalho ressalta a necessidade de ir além dos tratamentos estáticos de perda de massa para compreender a interação entre ventos internos e caminhos de formação planetária.