O resultado saiu em 16 de setembro de 2026 no The Astronomical Journal, assinado por doze autores de instituições dos Estados Unidos, da Alemanha, do Chile, da França e do Reino Unido, e carrega no título a palavra que resume a ambição do trabalho: um benchmark, isto é, um padrão de referência, um conjunto de medidas independentes contra o qual qualquer modelo de atmosfera estelar pode ser testado no infravermelho próximo. A escolha do termo não é modesta. Modelos de atmosfera são a engrenagem escondida de boa parte da astrofísica moderna, e o que a equipe do CHARA encontrou é que essa engrenagem gira com uma folga maior do que se supunha.

Para entender por que o limbo de uma estrela é mais escuro que o centro, basta pensar no que o olho faz ao mirar uma piscina. Olhando de cima, a vista atravessa a água até o fundo; olhando de raspão, a luz percorre um caminho longo dentro do líquido e a visão se perde nas camadas superficiais. Com a fotosfera acontece algo parecido. No centro do disco, a linha de visada entra perpendicularmente no gás e alcança camadas profundas e quentes antes de ser absorvida. Na borda, a mesma linha entra em ângulo rasante e só chega a camadas rasas, mais frias. Como gás mais frio brilha menos, a borda parece apagada. O perfil do esmaecimento, do centro ao limbo, é portanto um registro de como a temperatura cai com a altura e de quão opaco o gás é em cada comprimento de onda.

Esse registro vale ouro, porque quase tudo o que se mede sobre uma estrela passa por ele. Quando um telescópio determina o diâmetro angular de um astro, precisa corrigir o efeito de um disco que não é uniformemente brilhante. Quando um planeta transita diante de sua estrela e bloqueia uma fração da luz, a profundidade do trânsito depende de quão escura é a região do disco que ele cobre em cada instante. Binárias eclipsantes, lentes gravitacionais e a espectroscopia de atmosferas de exoplanetas com o Telescópio Espacial James Webb enfrentam a mesma dependência. Na prática, quase todos esses trabalhos usam coeficientes de escurecimento de limbo calculados por modelos, e não medidos, porque medir exige enxergar o disco da estrela, e o disco de quase qualquer estrela é pequeno demais para isso.

O problema é que a precisão das observações espaciais já ultrapassou a precisão desses modelos. Dados do Webb, como os de Zafar Rustamkulov e colaboradores em 2023, alcançam um nível em que a incerteza sistemática na lei de escurecimento adotada deixa de ser desprezível, e trabalhos de Louis-Philippe Coulombe e de Kuldeep Verma, ambos em 2024, mostraram o quanto a escolha dessa lei afeta o que se extrai de um espectro de trânsito. Ao mesmo tempo, as grades de modelos amplamente usadas discordam entre si, e medidas empíricas vêm apontando desencontros com as previsões desde o levantamento de Antonio Claret e Steven Bloemen, em 2011, até as simulações tridimensionais de Zazralt Magic, em 2015, e as medições interferométricas de Pierre Kervella, em 2017.

A maioria das medidas empíricas vem de trânsitos, binárias eclipsantes e eventos de microlente. Todas sondam o disco estelar por via indireta, ocultando uma parte dele ou ampliando regiões diferentes da superfície em momentos diferentes. A interferometria óptica de longa linha de base faz algo mais direto. Em vez de formar uma imagem, ela combina a luz captada por telescópios separados por dezenas ou centenas de metros e mede o contraste das franjas de interferência, a chamada visibilidade. A visibilidade é a transformada de Fourier da distribuição de brilho no céu, e seu quadrado, o V², a grandeza que os instrumentos registram, amostra o perfil radial de brilho da estrela. Um disco de bordas escuras e um disco uniforme produzem curvas de visibilidade diferentes.

A diferença, porém, não aparece em qualquer lugar da curva. Para um disco uniforme, a visibilidade cai da unidade até zero num ponto que depende do diâmetro angular, o primeiro nulo, e volta a crescer num segundo lóbulo bem mais fraco, depois num terceiro. O diâmetro é fixado pela posição dos nulos. A forma do esmaecimento, por sua vez, está escondida na curvatura e na altura dos lóbulos seguintes. Quem observa uma estrela apenas no primeiro lóbulo mede o tamanho dela e mais nada; quem alcança o segundo e o terceiro lóbulo começa a ler o perfil da borda. Isso exige linhas de base longas, estrelas grandes no céu e uma calibração precisa ao nível de uma parte em mil.

É exatamente essa combinação que o CHARA oferece. O conjunto, operado pela Universidade Estadual da Geórgia no Observatório de Mount Wilson, na Califórnia, reúne seis telescópios de 1 metro dispostos em Y e separados por até 331 metros, o que dá uma resolução angular da ordem de meio milissegundo de arco na banda H e de 0,7 milissegundo de arco na banda K. Para dimensionar o que isso representa, um milissegundo de arco é o tamanho aparente de uma moeda de um real vista a cerca de 5.500 quilômetros de distância. Em operação científica desde meados dos anos 2000, foi esse conjunto que produziu, em 2007, a primeira imagem resolvida da superfície de uma estrela da sequência principal além do Sol, Altair, publicada por John Monnier, um dos autores do novo trabalho.

Dois instrumentos trabalham juntos no laboratório de combinação de feixes do conjunto. O MIRC-X combina a luz dos seis telescópios na banda H, entre 1,5 e 1,75 micrômetro, e o MYSTIC faz o mesmo na banda K, entre 2,0 e 2,37 micrômetros. Os dois enxergam a mesma estrela no mesmo instante, pelas mesmas óticas e com os mesmos calibradores, o que elimina boa parte das diferenças sistemáticas que atormentam comparações entre observações feitas em noites ou instrumentos distintos. Cada exposição rende 15 medidas de V², uma por par de telescópios, e 10 fases de fechamento, ou closure phases, combinações de fases ao longo de triângulos de telescópios que sobrevivem à turbulência atmosférica e são sensíveis a assimetrias no objeto observado. O MYSTIC funcionou como rastreador de franjas, estabilizando a interferência enquanto o MIRC-X gravava os dados científicos.

A lista de alvos começou grande. A equipe partiu do Catálogo de Diâmetros Estelares do JMMC, o centro francês Jean-Marie Mariotti, e selecionou estrelas com diâmetros entre 2 e 6 milissegundos de arco, brilho inferior à quarta magnitude em H ou K, cores compatíveis com estrelas normais, declinação acima de 15 graus sul e nenhuma suspeita de binaridade ou pulsação. Sobraram 744 candidatas observáveis de Mount Wilson. Dessas, 50 foram observadas, e 19 caíram na peneira: cinco mostraram sinais de estrutura no brilho superficial, duas se revelaram binárias e doze não tinham cobertura ou qualidade de dados suficientes para prender o perfil de escurecimento. A amostra final tem 31 estrelas.

São estrelas brilhantes e próximas, dos tipos espectrais F ao M, e quase todas evoluídas: subgigantes, gigantes, gigantes luminosas e supergigantes. As temperaturas efetivas vão de 3.690 a 6.570 kelvin, a gravidade superficial cobre quase quatro ordens de grandeza, e os raios vão de 2,05 vezes o raio do Sol, no caso de Procyon, a pouco mais de 300 raios solares, no caso de V424 Lacertae, uma supergigante K5 que, posta no lugar do Sol, engoliria a órbita da Terra. Entre os nomes, há velhos conhecidos dos atlas celestes: Mirfak, a estrela mais brilhante de Perseu; Schedar, em Cassiopeia; Mebsuta, em Gêmeos; Ain, nas Híades; Rasalas, na cabeça do Leão; Errai, em Cefeu; e Muphrid, vizinha de Arcturus no Boieiro.

Procyon merece um parêntese. A estrela mais brilhante do Cão Menor é uma subgigante F5 a apenas 11,5 anos-luz, com uma anã branca como companheira, e é a única da amostra com gravidade superficial de estrela da sequência principal. Com 5,417 milissegundos de arco de diâmetro, ela é também um dos alvos mais bem resolvidos: os dados do CHARA cruzam pelo menos três nulos de visibilidade na banda H, o que permite distinguir as diferentes leis de escurecimento diretamente no domínio da medida. Ela ancora a ponta quente da amostra, e seu perfil serve de contraste para as gigantes frias que dominam a lista.

A preferência por gigantes tem razão técnica. Para medir o escurecimento, e não apenas o tamanho, os dados precisam ir além do primeiro nulo e entrar no segundo lóbulo da curva de visibilidade. Na amostra final, todas as 31 estrelas têm cobertura até o segundo nulo na banda H, e 24 delas também na banda K, onde o comprimento de onda maior reduz a resolução. Simulações feitas pela própria equipe, com dados sintéticos gerados na mesma configuração de seis telescópios, mostraram que a recuperação do coeficiente de escurecimento só se torna robusta a partir de cerca de 1,8 milissegundo de arco em H e 2,3 em K; o menor alvo da amostra mede 2,2 milissegundos de arco, e o ajuste conjunto das duas bandas segura a solução mesmo nos casos limítrofes.

As observações foram feitas entre 2025 e 2026, depois de uma correção importante no próprio CHARA. No início de 2025, a equipe do observatório eliminou um problema de tremulação temporal nas linhas de atraso, os trilhos ópticos que compensam a diferença de caminho da luz entre os telescópios, reduzindo os erros residuais de rastreamento para menos de 20 nanômetros. Esse número, que corresponde a perdas instrumentais de visibilidade abaixo de 1% nas duas bandas, é o que permite afirmar que a curvatura medida nos lóbulos vem da estrela e não da máquina. O protocolo de observação alternou estrelas científicas e calibradoras numa cadência fixa, com varreduras coerentes de cinco minutos amostradas a cada 4 milissegundos.

A calibração recebeu cuidado proporcional ao que estava em jogo. Os 42 calibradores, estrelas de disco minúsculo que servem de régua para a resposta do instrumento, foram triados com os indicadores astrométricos do Gaia EDR3 para excluir binárias e tiveram seus diâmetros refinados calibrando-os uns contra os outros nas mesmas noites. A escala de comprimento de onda, que fixa onde caem os nulos e, portanto, o diâmetro, foi monitorada com um módulo de etalon chamado ARMADA, que injeta nos feixes um sinal parecido com o de uma binária artificial. Os fatores de correção resultantes diferem da unidade por frações de 1% e carregam incertezas de poucas partes em dez mil, o tipo de minúcia que separa uma medida de 0,3% de uma de 3%.

Com os dados em mãos, a equipe recorreu ao PMOIRED, pacote Python de modelagem paramétrica de dados interferométricos desenvolvido por Antoine Mérand, do Observatório Europeu do Sul e coautor do estudo. Para cada estrela, as visibilidades das bandas H e K foram ajustadas em conjunto com um único diâmetro angular e um único fator de normalização, mas com coeficientes de escurecimento independentes em cada banda. A escolha é consequente: o diâmetro físico da estrela é o mesmo nas duas cores, enquanto o perfil de brilho não precisa ser. Diâmetros de disco uniforme em cada banda foram ajustados à parte, como referência, e as incertezas vieram de reamostragens do tipo bootstrap, em que as varreduras de dez minutos de cada alvo são sorteadas com reposição e reajustadas milhares de vezes.

Quatro leis analíticas descreveram o perfil de brilho, cada uma com sua linhagem. A lei linear remonta a Karl Schwarzschild, em 1906; a quadrática, a Zdeněk Kopal, em 1950; a lei de potência, a Daniel Hestroffer, em 1997; e a chamada power-2, uma lei de potência com dois parâmetros, a Pierre Maxted, em 2018. A equipe elegeu a lei de potência de um parâmetro como referência para as comparações. Nela, o brilho em cada ponto do disco é o cosseno do ângulo entre a linha de visada e a normal à superfície, elevado a um expoente alfa. Alfa funciona como um botão de intensidade do efeito: zero corresponde a um disco uniforme, e valores maiores escurecem a borda com mais vigor.

Do outro lado da balança ficaram cinco grades de modelos de atmosfera, escolhidas para cobrir filosofias diferentes de construção. A grade Kurucz ATLAS, com 741 modelos, representa a abordagem clássica unidimensional e plano-paralela, em que a atmosfera é uma pilha de camadas horizontais em equilíbrio. As grades MPS1 e MPS2, com 34.160 modelos cada, são versões modernas da mesma arquitetura, calculadas por Nadiia Kostogryz e colaboradores a partir de 2022, com misturas de abundâncias e tratamentos de convecção distintos. A grade Stagger, com 99 modelos, vem de simulações tridimensionais de radiação e hidrodinâmica, nas quais a convecção aparece como células reais de gás subindo e descendo. E a grade SATLAS, de Hilding Neilson e John Lester, é esférica, pensada para gigantes cuja atmosfera estendida não cabe na aproximação de camadas planas.

Comparar medida e modelo exige mais do que abrir uma tabela de coeficientes. Em vez de usar valores pré-calculados, a equipe pegou os perfis de intensidade de cada grade, integrou-os sobre as bandas passantes exatas do MIRC-X e do MYSTIC e ajustou a eles as mesmas leis analíticas usadas nos dados. Para propagar a incerteza dos parâmetros estelares, cada estrela foi sorteada 10 mil vezes com desvios de 100 kelvin na temperatura e de 0,2 dex na gravidade e na metalicidade, e o coeficiente adotado foi a mediana dessas realizações. Os perfis foram ajustados até o limite nativo de cada grade, bem perto do limbo geométrico, porque é lá que os dados do CHARA são sensíveis; truncar o ajuste antes da borda, como se faz em algumas tabelas para estabilizar os coeficientes, deixaria de ser uma comparação de igual para igual.

Havia ainda uma armadilha sutil. Um coeficiente obtido ajustando uma lei diretamente ao perfil de intensidade de um modelo não é a mesma coisa que o coeficiente que o CHARA mediria se aquela estrela tivesse exatamente aquele perfil, porque o ajuste no domínio das visibilidades pesa as regiões do disco de outro jeito e só enxerga o perfil nas frequências espaciais que foram amostradas. A equipe contornou o problema com um procedimento batizado de SVAM, sigla em inglês para modelo de atmosfera em visibilidades sintéticas: converteu cada perfil de modelo em visibilidades sintéticas, amostrou-as exatamente nos mesmos pontos do plano de frequências, nos mesmos comprimentos de onda e com as mesmas incertezas dos dados reais, e refez o ajuste. É como julgar duas fotografias pela mesma lente, com a mesma abertura, em vez de comparar uma foto com uma pintura.

Os diâmetros vieram primeiro, e vieram com precisão formal mediana entre 0,2% e 0,3%. Procyon ficou em 5,417 ± 0,005 milissegundos de arco; Schedar, em 5,664 ± 0,010; ζ Cephei, em 5,239 ± 0,013. Convertidos em raios físicos com as distâncias do Gaia, os números dão 2,05 raios solares para Procyon, 43 para Schedar e 172 para ζ Cephei, com a incerteza dominada pela paralaxe, não pela interferometria. Como esperado, os diâmetros com escurecimento são maiores que os de disco uniforme, porque um disco de bordas apagadas precisa ser um pouco maior para produzir a mesma curva de visibilidade; a correção mediana foi de 3,4% em H e 3,2% em K, e diminui para estrelas mais quentes.

O coração do trabalho, porém, está nos coeficientes. Em todas as 31 estrelas, o alfa medido na banda H é maior que o alfa medido na banda K. A direção é a esperada: no infravermelho, a opacidade e a emissão do gás variam menos com a temperatura, de modo que a diferença entre centro e borda se suaviza à medida que o comprimento de onda cresce. O que não era esperado foi o tamanho do degrau. Na amostra do CHARA, o coeficiente cai, em mediana, 38,15% de H para K, com dois terços das estrelas entre 33% e 46%. As grades de modelos, submetidas ao mesmo procedimento, preveem quedas entre 17% e 22%. A grade MPS2, usada como referência por sua cobertura densa e regular, prevê 19,78%, pouco mais da metade do valor medido.

Olhando banda por banda, o desencontro se concentra em H. Ali, os coeficientes medidos são sistematicamente maiores que os de todas as cinco grades: 21% acima da MPS2 e da MPS1, 27% acima da Kurucz, 29% acima da Stagger e 13% acima da SATLAS, em valores medianos. Em K, as diferenças encolhem e mudam de sinal conforme a grade: o CHARA fica 5% abaixo da MPS2, 5% abaixo da Kurucz, 6% abaixo da MPS1, 1% acima da Stagger e 12% abaixo da SATLAS. Em números absolutos, a mediana do alfa medido em H supera a previsão da MPS2 em 0,0465, e em K fica 0,0107 abaixo. Nenhuma das grades acerta ao mesmo tempo o nível de H, o nível de K e o salto entre os dois.

A estrela υ Persei, uma gigante K3 de 4.374 kelvin e 3,615 milissegundos de arco, ilustra bem o que os números dizem. Observada em várias épocas e resolvida com folga além do primeiro nulo, ela permite sobrepor aos dados a curva do ajuste empírico e as curvas previstas pelos perfis fixos das grades Stagger e SATLAS, com o diâmetro deixado livre em cada caso. No segundo e no terceiro lóbulos, onde a curvatura responde ao limbo, a lei de potência ajustada acompanha os pontos com qui-quadrado reduzido de 1,65, enquanto o perfil da SATLAS fica em 2,93 e o da Stagger em 3,18. Na banda K, as três curvas quase se confundem. É a assinatura do problema inteiro condensada numa única estrela: discrepância em H, concordância em K.

Antes de culpar os modelos, a equipe tentou culpar a si mesma. A primeira suspeita natural é que companheiras não detectadas ou manchas na superfície estivessem distorcendo a curvatura das visibilidades. As fases de fechamento resolvem essa dúvida: numa fonte simétrica, elas são zero, e qualquer assimetria as desloca. Para cada estrela, o PMOIRED varreu um campo de 100 a 270 milissegundos de arco à procura de companheiras e testou modelos com uma mancha na fotosfera. Nada apareceu acima de 3 sigma. O limite de detecção, de cerca de 5 magnitudes, implica que qualquer companheira na região teria de ser mais fraca que 1% do brilho da primária, e manchas isoladas acima desse mesmo 1% teriam sido vistas.

A segunda suspeita recai sobre a própria máquina de ajuste. Os coeficientes se mostraram reprodutíveis quando os dados foram reagrupados em médias de 30 segundos em vez de 10 minutos. Conjuntos sintéticos gerados com a cobertura real do CHARA e perfis de entrada tirados da MPS2 devolveram os coeficientes de entrada dentro das incertezas declaradas. Como metade dos calibradores não está rigorosamente irresolvida nas linhas de base mais longas da banda H, a equipe repetiu a análise com os dados de uma noite específica, 4 de abril de 2026, em que todos os calibradores mediam menos de 0,4 milissegundo de arco, e obteve o mesmo comportamento. Restringir os dados de H ao mesmo intervalo de frequências espaciais coberto em K tampouco mudou as tendências.

Restava a possibilidade de uma diferença instrumental entre MIRC-X e MYSTIC. Os dois combinadores recebem a luz dos mesmos feixes, passam pelo mesmo trem óptico, observam os mesmos calibradores na mesma cadência e são reduzidos no mesmo arcabouço, e nenhuma noite mostrou um deslocamento diferencial capaz de reproduzir o padrão de toda a amostra. Para sete estrelas observadas em mais de uma época, Mirfak, υ Persei, α Vulpeculae, λ Herculis, 109 Herculis, δ Boötis e β Draconis, os diâmetros se repetiram ao nível de poucos décimos de por cento e os coeficientes dentro das barras de erro. A dispersão de estrela para estrela, bem maior que isso, é portanto real, e não ruído de medida.

Há um ponto em que a equipe pede cautela na leitura. Os coeficientes são, nas palavras do artigo, resumos empíricos compactos da curvatura das visibilidades, e não descrições físicas únicas da atmosfera. Trocar a lei analítica desloca o valor absoluto de alfa, e mesmo o procedimento SVAM altera os coeficientes de modelo num nível perceptível, da ordem de um centésimo. Só que esses deslocamentos são pequenos diante das diferenças entre CHARA e modelos, e, decisivamente, a mesma tendência qualitativa reaparece quando se usam as outras leis. A discrepância não nasce da parametrização; ela sobrevive a ela.

E o diâmetro? Aqui vem a parte tranquilizadora do estudo. Quando o perfil empírico foi substituído pelo perfil fixo da MPS2, deixando livre apenas a escala angular, o diâmetro mudou em mediana 0,08%, com dispersão de 0,61%. Trocar a grade de referência também pouco altera: SATLAS, Stagger, Kurucz e MPS1 produzem diâmetros que diferem da MPS2 por 0,365%, 0,153%, 0,003% e 0,004%, respectivamente. As leis analíticas alternativas se comportam de modo parecido, exceto a linear, que subestima a curvatura e encolhe o diâmetro em 1,1%. A explicação é geométrica: o tamanho está preso nos nulos, que todas as grades reproduzem bem, enquanto o escurecimento vive nos lóbulos altos, onde elas tropeçam.

A conclusão prática tem duas faces. Para quem usa a interferometria para medir raios estelares, e deles deriva temperaturas efetivas e luminosidades ao combinar o diâmetro com a paralaxe e o fluxo bolométrico, os diâmetros do CHARA são robustos, e o estudo entrega uma escala homogênea de 31 gigantes e supergigantes brilhantes, com incertezas formais abaixo de 0,3%. Para quem usa modelos de atmosfera para prever o escurecimento de limbo de uma estrela evoluída no infravermelho próximo, o recado é desconfortável: a forma do perfil, em particular sua dependência com o comprimento de onda, não está sendo reproduzida, e o erro na banda H chega a um quinto do valor do coeficiente.

De onde vem o desencontro? A equipe examina os suspeitos sem apontar um culpado. As maiores discrepâncias em H aparecem de preferência entre as estrelas mais frias e mais evoluídas, abaixo de cerca de 5.000 kelvin, num eco do que Jayshil Patel e Néstor Espinoza encontraram em 2022 nos trânsitos do TESS, onde os coeficientes quadráticos empíricos de estrelas frias se afastam sistematicamente das previsões PHOENIX e ATLAS. Mas a temperatura sozinha não explica tudo: Mirfak, uma supergigante F5 de 6.193 kelvin, exibe um dos desvios mais marcantes da amostra. Classe de luminosidade, gravidade superficial, extensão da atmosfera e estruturas na superfície entram na lista de fatores prováveis.

A grade esférica SATLAS é o caso mais instrutivo. Por abandonar a hipótese de camadas planas, ela é a única a tratar a extensão geométrica da atmosfera de uma gigante, e de fato produz o menor desvio absoluto em H, de 13%. Só que ela passa do ponto em K, onde superestima o escurecimento em 12%, e com isso o contraste entre as bandas continua longe do medido. Geometria esférica, portanto, mexe na escala absoluta dos coeficientes, mas não resolve a dependência com o comprimento de onda. Algo mais, na física das camadas que emitem em H, precisa estar diferente do que as grades assumem.

Os candidatos são conhecidos de quem modela atmosferas frias. A convecção aparece na superfície de uma gigante como grânulos gigantescos, que modulam a temperatura e alteram o brilho de cada ponto do disco de um modo que nenhum modelo unidimensional captura por construção. Opacidade molecular, sobretudo de água e monóxido de carbono, muda a altura das camadas de onde a luz escapa em cada banda. Campos magnéticos e estruturas de superfície de baixo contraste podem achatar ou acentuar o perfil. E a própria quebra das hipóteses unidimensional e plano-paralela nas gravidades baixíssimas de supergigantes como V424 Lacertae, com log g de 0,32, coloca a maioria das grades fora de sua zona de validade: a Kurucz só cobre gravidades de estrelas anãs, e as MPS param em log g igual a 3.

Esse último ponto exige honestidade sobre o alcance da comparação. Para as estrelas de gravidade mais baixa, as previsões das grades Kurucz e MPS são extrapolações de borda, não modelos casados com a estrela. A Stagger alcança gravidades menores, mas é esparsa e irregular em temperatura, gravidade e metalicidade. A SATLAS assume composição solar e usa a massa, não a metalicidade, como parâmetro. A equipe prefere, por isso, ler as diferenças entre grades como sensibilidade dos modelos, e não como um ranking de qualidade. Os modelos PHOENIX ficaram de fora porque os ajustes das leis analíticas até o limbo não se mostraram estáveis, e grades como MARCS e NewEra são apontadas como extensões naturais para gigantes frias.

O CHARA não está sozinho em seu diagnóstico. Em 2017, Pierre Kervella e colaboradores, incluindo Alexandre Gallenne, coautor do novo artigo, usaram o instrumento PIONIER no interferômetro do VLT para medir o escurecimento de Alfa Centauri A e B na banda H e encontraram perfis mais suaves que os previstos por modelos unidimensionais e tridimensionais. Em 2026, Chowhan e colaboradores, com o instrumento PAVO no próprio CHARA, viram que a grade Stagger superestima o escurecimento da gigante κ Cygni na banda R. E as análises de trânsitos do Kepler, do TESS e do Webb, de Maxted em 2023 a Kostogryz em 2024, vêm acumulando desvios semelhantes. O que o novo trabalho acrescenta é um padrão de referência para estrelas evoluídas, no infravermelho próximo, em duas bandas simultâneas e para uma amostra inteira.

A relação com os trânsitos é de complementaridade, não de concorrência. Os estudos de trânsito restringem, por natureza, estrelas compactas da sequência principal com planetas de período curto; gigantes e subgigantes têm raios grandes, probabilidade de trânsito baixa e trânsitos rasos, o que torna o escurecimento quase impossível de extrair por ocultação. Binárias eclipsantes com componentes evoluídas exigem modelar distorção de maré, escurecimento gravitacional e reflexão. Até onde a equipe sabe, não existe nenhuma restrição publicada de escurecimento de limbo por trânsito, no TESS ou no Webb, para qualquer das 31 estrelas. Para esse regime de atmosfera, a interferometria é, por enquanto, a única régua direta disponível.

Há também limites que a própria equipe faz questão de listar. Não existe, para a amostra completa, monitoramento fotométrico ou espectroscópico simultâneo que permita quantificar atividade superficial e variabilidade temporal em escalas relevantes para convecção, manchas ou pulsação. As fases de fechamento rejeitam assimetrias fortes e companheiras, mas não descartam estruturas de baixo contraste ou quase simétricas em torno do eixo, que poderiam alterar o perfil inferido e contribuir para a dispersão de estrela para estrela. O próximo passo sugerido é vasculhar as curvas de luz do TESS dessas estrelas brilhantes, seguindo o método de Rudrasingam e colaboradores publicado em 2026, à procura de assinaturas de convecção, manchas ou pulsação que possam ser confrontadas com a interferometria.

O caminho adiante passa por mais comprimentos de onda. Um teste interno já foi feito: a equipe dividiu cada banda em duas metades, uma mais azul e outra mais vermelha, e refez os ajustes. Em H, a metade de comprimento de onda menor produz coeficientes maiores que a metade de comprimento de onda maior, exatamente o sentido esperado se o escurecimento cresce para o azul; em K, a diferença some dentro do ruído. Observações simultâneas em J, H e K com os mesmos dois instrumentos testariam se o desencontro com os modelos cresce rumo ao infravermelho mais curto, e o instrumento de luz visível SPICA, do próprio CHARA, onde o escurecimento é intrinsecamente mais forte, daria um teste ainda mais severo.

Também há a questão das estrelas mais frias. Gigantes vermelhas de tipo M tardio, abaixo de 3.700 kelvin, vivem no regime em que extensão esférica, opacidade molecular e superfície convectiva tornam o perfil de brilho mais sensível ao modelo adotado, e as grades esféricas recentes, como a NewEra de Peter Hauschildt e colaboradores, publicada em 2025, indicam que as diferenças em relação a modelos plano-paralelos são maiores justamente para gigantes. A amostra atual para em 3.690 kelvin, com α Vulpeculae. Estender a comparação para baixo, com mais épocas e cobertura espectral do R ao K, dirá se a tensão aumenta com temperaturas menores, gravidades menores e comprimentos de onda mais curtos, como as tendências atuais sugerem.

Tudo o que sustenta o trabalho é público. Os arquivos calibrados de visibilidade estão no banco de dados OIDB do JMMC; a tabela legível por máquina com os ajustes de todas as leis, os coeficientes de todas as grades e os diâmetros de perfil fixo acompanha a versão on-line do artigo; e o código, as figuras completas, 42 pares de estrela e época em 84 gráficos por banda, e os resíduos de fase de fechamento estão depositados no Zenodo. A ideia explícita dos autores é que outros grupos possam testar suas grades de atmosfera contra esses 31 perfis sem repetir as observações, e que a própria comunidade de exoplanetas use os números como teste de estresse para as prescrições que adota.

Há uma elegância silenciosa na forma como esse resultado foi obtido. Ninguém fotografou a borda de Mirfak ou de Schedar. O que o CHARA registrou foram padrões de interferência, oscilações de contraste em franjas de luz combinadas sobre centenas de metros de Mount Wilson, e foi na curvatura discreta do segundo e do terceiro lóbulos dessas oscilações que o perfil da borda se revelou. É uma medida feita quase inteiramente no espaço das frequências, traduzida depois para o espaço das imagens. As figuras de disco que acompanham o artigo, com υ Persei brilhando mais no centro e esmaecendo na borda, são reconstruções de modelo; o dado em si é mais abstrato, e por isso mesmo mais difícil de enganar.

Volte à folha de papel com o Sol projetado. A borda escura que você vê ali é a mesma física que agora foi medida em 31 estrelas a dezenas e centenas de anos-luz, com um resultado que não cabe nos modelos: entre a banda H e a banda K, o limbo dessas estrelas muda quase duas vezes mais do que a teoria prevê. Os diâmetros resistem, as escalas de raio se mantêm, e nada disso derruba a astrofísica das atmosferas estelares. Mas algo nas camadas que emitem a 1,6 micrômetro, nas gigantes frias em especial, está acontecendo de um modo que ainda não sabemos descrever. O limbo de uma estrela é, por definição, o lugar onde ela termina. Para os modelos, é onde a verificação começa.