Physics of Eclipsing Binaries. VI. Hot, compact stars

Resumo Técnico: Física de Binárias em Eclipse. VI. Estrelas Quentes e Compactas

Definição do Problema
A modelagem precisa de curvas de luz de binárias em eclipse requer a atribuição de intensidades emergentes apropriadas aos elementos superficiais dos componentes binários. Embora códigos como o PHOEBE utilizem atmosferas de modelos para derivar essas intensidades com base na temperatura efetiva ($T_{\text{eff}}$), gravidade superficial ($\log g$) e composição química, as implementações existentes foram historicamente limitadas a tipos estelares específicos. Versões anteriores dependiam dos modelos de Castelli & Kurucz (CK2004) e PHOENIX, que são mais adequados para estrelas da sequência principal e gigantes de baixa temperatura. Consequentemente, para estrelas quentes e compactas (como anãs brancas quentes e estrelas subanãs O), os modelos frequentemente recorrem a aproximações de corpo negro. Esta aproximação é demonstradamente deficiente para esses objetos, levando a imprecisões significativas nos parâmetros físicos derivados. Além disso, versões anteriores do código falhariam inteiramente se elementos superficiais específicos (devido à irradiação, distorção de maré ou rotação rápida) caíssem fora dos limites de parâmetros das grades de modelos atmosféricos disponíveis.

Metodologia
Este artigo detalha a implementação de novas capacidades no lançamento do PHOEBE 2.5, focando na integração de atmosferas de modelos não-LTE para estrelas quentes e compactas e no desenvolvimento de um tratamento robusto para parâmetros fora dos limites.

1. Integração de Novas Atmosferas:

   • Grades TMAP: Os autores incorporaram atmosferas de modelos de plano paralelo do Tübingen Model Atmosphere Package (TMAP). Estas cobrem quatro grades específicas: sdO, DA, DAO e DO, abrangendo $T_{\text{eff}}$ de 20 kK a 200 kK e $\log g$ de 4,0 a 9,5. Um novo parâmetro, loghefrac (fração numérica logarítmica de Hélio para Hidrogênio), foi introduzido para definir abundâncias para estas grades, distinguindo-as das definições baseadas em metalicidade dos modelos anteriores.
   • Grades Montréal/Tremblay: O lançamento também inclui modelos LTE Montréal/Tremblay para anãs brancas DA, cobrindo $T_{\text{eff}}$ de 3,75 kK a 60 kK. Os autores excluiram explicitamente os modelos não-LTE TLUSTY de Claret et al. (2020a) devido a erros de cálculo identificados nas intensidades fornecidas naquele trabalho.

2. Mistura (Blending) e Extrapolação:
   Para resolver o problema de elementos superficiais que caem fora dos limites da grade, o PHOEBE 2.5 implementa dois novos modos operacionais:

   • Extrapolação: Os usuários podem selecionar 'mais próximo' (usando o nó definido mais próximo) ou 'linear' (extrapolação linear a partir do hipercubo totalmente definido mais próximo) para estimar valores para parâmetros fora dos limites.
   • Mistura (Blending): Este modo cria uma transição suave entre os valores de atmosfera de modelo extrapolados e os valores de corpo negro. A razão de mistura depende da distância da grade de atmosfera de modelo, evitando que o código falhe quando apenas alguns elementos superficiais estão fora dos limites.

3. Funcionalidades Adicionais:

   • Boosting Doppler: A funcionalidade de boosting relativístico (beaming) foi reativada. Diferente de versões anteriores, os usuários devem agora fornecer fatores de boosting calculados externamente (por exemplo, da literatura), em vez de depender de interpolação interna.
   • Recursos Definidos pelo Usuário: Uma nova arquitetura permite que os usuários injetem código personalizado na lógica do PHOEBE. Isso possibilita modificações na malha, valores locais (velocidades radiais, temperaturas) ou no modelo sintético final, facilitando recursos como manchas migratórias, rotação diferencial e terceira luz variável.

Resultados Principais e Validação

• Consistência Atmosférica: Os autores verificaram a consistência entre as grades TMAP e Tremblay. Comparações de intensidades normais e coeficientes de escurecimento de limbo na banda SDSS $g'$ mostram concordância dentro de um desvio padrão da unidade (razão $1,01 \pm 0,10$). A transição entre as quatro grades TMAP (sdO, DA, DAO, DO) mostrou-se contínua em relação aos coeficientes de escurecimento de limbo e intensidade normal como função da abundância.
• Impacto da Aproximação de Corpo Negro: Um sistema binário sintético (uma pré-anã branca de 127 kK com uma companheira K5V) foi modelado. O estudo descobriu que o uso de uma aproximação de corpo negro para o componente quente resulta em uma subestimação da temperatura efetiva de aproximadamente 25% (exigindo uma redução para ~94 kK para coincidir com a curva de luz TMAP). Além disso, reproduzir a amplitude do efeito de irradiação com um corpo negro exigiu um albedo bolométrico não físico de 2,0 para a companheira, comparado a um valor físico de 0,6 no modelo TMAP.
• Coeficientes de Escurecimento de Limbo: Os autores geraram e forneceram tabelas de coeficientes de escurecimento de limbo derivados dos novos modelos TMAP e Tremblay para conjuntos de filtros padrão (SDSS e Johnson), adequados para uso em outros códigos de modelagem binária.

Significância
O artigo afirma que a inclusão das atmosferas TMAP e Tremblay, juntamente com as capacidades de mistura/extrapolação, melhora significativamente a fidelidade do PHOEBE para estrelas quentes e compactas. A principal significância reside na eliminação da dependência da aproximação de corpo negro para esses objetos, o que os autores demonstram poder levar a erros substanciais nos parâmetros estelares derivados (especificamente temperatura efetiva e luminosidade). Ao permitir que o código lide graciosamente com elementos superficiais que caem fora dos limites da grade, o novo lançamento permite a modelagem de sistemas complexos envolvendo forte irradiação, distorção de maré e rotação rápida sem recorrer a estados de erro ou aproximações imprecisas. Adicionalmente, a reativação do boosting Doppler e a introdução de recursos definidos pelo usuário expandem a utilidade do código para a análise de binárias compactas e fenômenos estelares dinâmicos.