Formation and X-ray emission from hot bubbles in planetary nebulae - III. The impact of [Wolf-Rayet]-type winds

Resumo Técnico: Formação e Emissão de Raios X de Bolhas Quentes em Nebulosas Planetárias - III. O Impacto de Ventos do tipo [Wolf-Rayet]

Definição do Problema
As nebulosas planetárias (PNe) exibem bolhas quentes, emissoras de raios X, formadas pela interação de ventos rápidos de pós-AGB (pós-Ramo da Ascensão no Gigante Assintótico) com ejectas de AGB mais lentas e densas. Enquanto modelos padrão assumem estrelas centrais ricas em hidrogênio (H-ricas), uma fração significativa de PNe que emitem raios X hospeda estrelas centrais do tipo [Wolf-Rayet] ([WR]), deficientes em hidrogênio. Estudos anteriores sugeriram que a alta metalicidade (rica em C) dos ventos [WR] aumenta o resfriamento radiativo, podendo atrasar a formação da bolha e alterar as estruturas térmicas. No entanto, os modelos existentes frequentemente basearam-se em simulações 1D com aumentos artificiais de vento ou assumiram a condução térmica como o principal mecanismo para resfriar o plasma para as temperaturas observadas de raios X suaves ($T_X \approx 1-3 \times 10^6$ K). O impacto específico dos parâmetros de vento [WR] autossistentes, derivados de prescrições empíricas recentes de perda de massa, sobre a evolução hidrodinâmica e a emissão sintética de raios X das bolhas quentes ainda carece de uma exploração completa em um contexto multidimensional.

Metodologia
Os autores utilizam simulações de radiação-hidrodinâmica usando o código PLUTO para modelar o meio circunestelar de PNe. A metodologia envolve uma abordagem de múltiplos estágios:

1. Evolução Estelar: Modelos de evolução para progenitores de $1, 2$e$3 M_\odot$ são gerados usando o MESA. Duas trajetórias evolutivas são computadas para cada massa: uma trajetória de pós-AGB padrão rica em H (usando a prescrição de perda de massa de Miller Bertolami 2016) e uma trajetória [WR] utilizando uma prescrição recente de taxa de perda de massa para estrelas deficientes em H, derivada de Toalá, Todt & Sander (2024).
2. Configuração Hidrodinâmica:
   • Simulações 1D: Realizadas em grades radiais para estudar a formação e evolução global da bolha, incorporando curvas de resfriamento específicas para abundâncias ricas em H e deficientes em H (ricas em C) geradas via Cloudy.
   • Simulações 2D: Realizadas em grades esféricas para capturar instabilidades hidrodinâmicas (camada fina de Vishniac e Rayleigh-Taylor) na interface de interação vento-vento. Estas simulações rastreiam a mistura do material AGB e do vento pós-AGB usando um parâmetro escalar advectado.
3. Síntese de Raios X: Espectros sintéticos de raios X (0,3–2,0 keV) são calculados via pós-processamento usando o banco de dados atômico CHIANTI. A Medida de Emissão Diferencial (DEM) é computada separadamente para as fases de gás ricas em H e pobres em H baseada no escalar advectado, permitindo o cálculo das temperaturas de plasma ($T_X$) e luminosidades ($L_X$) ponderadas pela emissividade, que contabilizam a composição mista da bolha.

Principais Contribuições

• Evolução [WR] Autossistente: Ao contrário de trabalhos anteriores que aumentavam artificialmente os parâmetros de vento para uma única massa estelar, este estudo utiliza trajetórias MESA com taxas de perda de massa [WR] autossistentes através de uma gama de massas de progenitores ($1-3 M_\odot$).
• Mistura Multidimensional: O estudo evita explicitamente a condução térmica, demonstrando que as instabilidades hidrodinâmicas em simulações 2D geram naturalmente camadas de mistura turbulenta que regulam as temperaturas das bolhas.
• Análise de Raios X de Componente Dupla: Os autores introduzem um método rigoroso para separar e recombinar as contribuições de raios X das fases de gás ricas em H e pobres em H, reconhecendo que as bolhas [WR] contêm uma mistura de materiais quimicamente distintos.

Resultados

• Atraso na Formação da Bolha: O resfriamento radiativo aumentado nos ventos [WR] ricos em C atrasa significativamente a formação de bolhas quentes em comparação com as contrapartes ricas em H. Para um modelo de $1 M_\odot$, a formação da bolha é atrasada em aproximadamente 6.000 anos após o início do pós-AGB, enquanto modelos de massa superior mostram atrasos menos dramáticos. Este atraso é atribuído à necessidade de que a velocidade do vento exceda um limiar crítico ($v_{crit}$) para superar o rápido resfriamento radiativo.
• Luminosidade de Raios X Aumentada: Modelos que incorporam ventos [WR] exibem luminosidades de raios X ($L_X$) significativamente maiores do que os modelos ricos em H. Isso é impulsionado pela maior luminosidade mecânica do vento ($L_{wind}$) resultante das elevadas taxas de perda de massa das estrelas [WR].
• Convergência da Temperatura do Plasma: Apesar da maior potência do vento e das diferentes composições químicas, a temperatura do plasma ponderada pela emissividade ($T_X$) do gás emissor de raios X converge para valores de $1-3 \times 10^6$ K para ambos os modelos [WR] e ricos em H, uma vez que a mistura é estabelecida.
• Papel da Mistura vs. Condução Térmica: As simulações 2D mostram que a mistura turbulenta transporta eficientemente o material nebular mais frio para dentro da bolha quente, regulando a temperatura para os níveis observados sem a necessidade de condução térmica. Comparações com modelos 1D revelam que as simulações 1D superestimam sistematicamente $T_X$ (por fatores de 1,2–1,6, até ~5 em modelos massivos) devido à ausência deste mecanismo de mistura.
• Brilho de Superfície: O brilho de superfície de raios X atinge o pico precocemente e declina rapidamente (em cerca de 4 ordens de magnitude) conforme a bolha se expande, sugerindo que a emissão difusa detectável de raios X é restrita a uma janela evolutiva específica, onde a bolha se expandiu o suficiente, mas antes que ela resfrie e se dilua demais.

Significância e Alegações
O artigo afirma reforçar a hipótese de que a mistura hidrodinâmica é um mecanismo chave para gerar a emissão de raios X suave em PNe, mesmo para sistemas que hospedam estrelas centrais [WR]. Os autores afirmam que seus resultados:

1. Explicam a ausência de emissão de raios X em PNe [WR] jovens e compactas como uma consequência do atraso na formação da bolha devido ao forte resfriamento radiativo.
2. Demonstram que as temperaturas de raios X suaves observadas ($T_X \approx 1-3 \times 10^6$ K) são um resultado natural da mistura multidimensional, independentemente de a estrela central seguir um caminho evolutivo rico em H ou [WR].
3. Fornecem uma alternativa robusta aos modelos de condução térmica, sugerindo que a emissão de raios X suaves pode surgir mesmo na presença de campos magnéticos que poderiam suprimir a condução térmica de elétrons.
4. Validam o uso de prescrições empíricas de perda de massa [WR] na modelagem hidrodinâmica, mostrando que reproduzem as propriedades do vento estelar observadas ($L, T_{eff}, \dot{M}, v_\infty$) de CSPNe emissoras de raios X.

Os autores concluem que, embora sua aproximação de vento seja prática, trabalhos futuros deveriam idealmente modelar explicitamente o cenário "born-again" para capturar totalmente a transição para ventos deficientes em H. No entanto, a abordagem atual fornece um arcabouço fisicamente motivado para compreender as propriedades de raios X de PNe [WR].