Alumínio comprimido expõe erro no modelo usado para ler matéria quente e densa

Nem sólido comum, nem plasma ideal

Comprima um metal até densidades acima da normal e aqueça-o a milhares de kelvin ao mesmo tempo. O que resulta não é um sólido, porque a rede cristalina se desfez, nem um plasma no sentido usual, porque os elétrons continuam apertados demais para se comportarem como partículas livres. Esse território intermediário recebeu nome próprio: matéria quente e densa, ou warm dense matter na literatura em inglês.

A dificuldade não é semântica, é aritmética. Duas grandezas medem o quanto um sistema desses escapa das aproximações confortáveis. O parâmetro de acoplamento de Coulomb, Γ, compara a energia potencial e a energia térmica por partícula. O parâmetro de degenerescência, Θ, compara a energia térmica e a energia de Fermi dos elétrons. Na matéria quente e densa os dois ficam da ordem de 1 — e é exatamente aí que nenhuma das duas tradições teóricas tem por onde começar, como resume o grupo de teoria de matéria densa do Laboratório Nacional de Los Alamos.

• Γ da ordem de 1: as partículas não estão nem livres, como em um plasma quente e rarefeito, nem presas em posições fixas, como em um cristal frio.
• Θ da ordem de 1: os elétrons não formam nem um gás clássico, sujeito à estatística de Boltzmann, nem um gás totalmente degenerado como o de um metal à temperatura ambiente.
• Sem parâmetro pequeno: as duas expansões perturbativas usuais — acoplamento fraco de um lado, degenerescência completa do outro — perdem o termo em que se apoiariam.

Bonitz e colaboradores descreveram o problema em 2020 como o aparecimento simultâneo de degenerescência quântica, correlações de Coulomb e efeitos térmicos, com a sobreposição das fases de plasma e condensada. Daí a importância de uma propriedade específica: a resposta da densidade eletrônica a uma perturbação externa. É ela que a teoria precisa acertar e é ela que o espalhamento Thomson de raios X mede de rotina, como registra a revisão de Dornheim e colegas publicada em 2023. Esse cruzamento é o que dá sentido a um teste como o feito no European XFEL.

Oito nanossegundos para criar o estado, 25 femtossegundos para fotografá-lo

A primeira pergunta diante de um resultado assim é onde ele foi produzido. Não foi em uma célula de bigornas de diamante, o aparelho que aperta uma amostra minúscula entre duas pontas de diamante e consegue mantê-la sob pressão por dias. Foram duas máquinas trabalhando em série: um laser de alta potência criou o estado, e um laser de elétrons livres de raios X o interrogou antes que ele se desfizesse.

O palco é o instrumento HED, de High Energy Density, do European XFEL, em Schenefeld, na Alemanha, operado junto com o consórcio de usuários HIBEF. Segundo a descrição técnica publicada por Zastrau e colegas em 2021, o instrumento recebe fótons de 5 a 25 keV do ondulador SASE2, em pulsos de 50 femtossegundos ou menos, e dispõe do DiPOLE 100-X, um laser bombeado por diodos capaz de entregar 100 joules em um pulso de 10 nanossegundos, a até 10 disparos por segundo.

1. Alvo: uma folha de alumínio de cerca de 50 micrômetros colada a um ablador de Kapton de cerca de 25 micrômetros.
2. Impulso: o DiPOLE, com a frequência dobrada para 515 nanômetros, entrega um pulso quadrado de 8 nanossegundos e cerca de 26 joules nominais, dos quais aproximadamente 19 são absorvidos. O Kapton vaporiza e o recuo empurra uma onda de choque para dentro do metal.
3. Cronometragem: um interferômetro VISAR registra a saída do choque pela face oposta em 6,5 ± 0,1 nanossegundo, valor estável ao longo de toda a série de disparos.
4. Sonda: em 6,3 nanossegundos, cerca de 200 picossegundos antes dessa saída, chega o pulso de raios X — 25 femtossegundos, 8.307 elétron-volts, em configuração autossemeada de banda estreita.
5. Coleta: um espectrômetro von Hámos com cristal HAPG orientado no plano (004) registra o espectro espalhado, com resolução efetiva em torno de 3,5 elétron-volts. Foram quatro ângulos de espalhamento e cerca de 40 disparos acumulados por geometria.
6. Calibração: a função de fonte e instrumento foi medida à parte, em alvos de alumínio frio, e descrita pela soma de dois perfis de Voigt.

A assimetria entre as escalas de tempo é o que torna a medida possível. O pulso de raios X dura cerca de 250 mil vezes menos que o intervalo entre o início do choque e o instante em que ele é sondado; para efeito da medida, a matéria comprimida está parada. O preço é que cada disparo destrói o alvo — acumular 40 disparos em cada uma das quatro geometrias só é viável com um laser de bombeio capaz de repetir o ciclo várias vezes por segundo, e não uma vez por dia.

O estado alcançado não foi lido diretamente de um instrumento, e sim reconstruído. Simulações unidimensionais de hidrodinâmica radiativa, com o código HELIOS-CR, foram calibradas até reproduzir o tempo de saída visto no VISAR; a difração de raios X coletada no mesmo disparo verificou o fator de estrutura iônico contra simulações de dinâmica molecular por teoria do funcional da densidade. O resultado é um platô de cerca de 30 micrômetros de alumínio comprimido, com densidade entre 3,75 e 4,5 gramas por centímetro cúbico — o alumínio comum tem 2,70 —, a aproximadamente 50 gigapascais, equivalentes a 500 mil atmosferas segundo o HZDR, e temperatura eletrônica em torno de 0,6 elétron-volt, ou cerca de 7 mil kelvin.

O que um plásmon denuncia sobre os elétrons

O sinal que interessa chama-se plásmon. Os elétrons de um metal não respondem individualmente a uma perturbação: eles oscilam em conjunto, como um único fluido carregado, e essa oscilação coletiva tem frequência própria. Quando um fóton de raios X atravessa a amostra e é espalhado, pode ceder energia justamente para excitar essa oscilação. No espectro do feixe espalhado aparece então um pico deslocado, e a posição desse pico é a energia do plásmon.

A novidade do experimento está no adjetivo que abre o título do artigo: resolvido em momento. A energia do plásmon não é um número único — ela varia com o vetor de onda transferido no espalhamento, que por sua vez depende do ângulo em que o detector está posicionado. Medindo em quatro ângulos, a equipe cobriu vetores de onda de 0,99 a 2,57 por angstrom e obteve uma curva inteira, a relação de dispersão, em vez de um ponto. E não apenas a posição do pico: também a forma da linha, assimétrica e alargada pelo amortecimento de Landau, o processo pelo qual a oscilação coletiva transfere energia a elétrons individuais.

Uma curva é um teste bem mais duro que um ponto isolado: um modelo errado quase sempre pode ser ajustado para acertar um valor único, mas reproduzir posição e forma ao longo de toda a faixa de vetores de onda é outra exigência.

Nossas medidas são precisas o bastante para distinguir com clareza entre modelos concorrentes. Isso importa porque esses modelos são amplamente usados para diagnosticar estados extremos da matéria.

Os modelos erraram, e sempre para o mesmo lado

O procedimento padrão para ler um espectro de espalhamento Thomson de raios X passa pela decomposição de Chihara, que separa a contribuição dos elétrons ligados da dos elétrons livres, e por modelos de gás de elétrons uniforme para descrever a parte livre. É dessa cadeia que saem, rotineiramente, a densidade e a temperatura atribuídas a um experimento de matéria quente e densa. Foi essa cadeia que a medida colocou à prova.

A direção do erro é tão informativa quanto o tamanho. Os modelos de gás uniforme não erram para os dois lados de maneira aleatória: erram sempre por cima, o que caracteriza viés sistemático, não ruído estatístico.

A explicação apontada pelos autores não está em alguma correlação eletrônica que faltasse ao gás uniforme, e sim no que o gás uniforme simplesmente não tem: íons. O choque apaga a ordem cristalina de longo alcance e localiza os elétrons de valência em camadas de coordenação distorcidas em torno dos núcleos. A 7 mil kelvin — muito acima dos 660,3 graus Celsius em que o alumínio derrete à pressão ambiente — o metal já não é uma rede periódica, e sim um líquido denso e desordenado. É essa desordem, e não uma correção de segunda ordem no gás de elétrons, que fixa a dispersão e o amortecimento do plásmon.

Só quando levamos em conta a estrutura desordenada real é que teoria e experimento concordam.

O corolário prático é desconfortável para quem depende da rota rápida. Se o modelo de análise desloca o plásmon, a densidade e a temperatura extraídas dele saem enviesadas, e com elas as propriedades derivadas desse estado: opacidade, condutividade elétrica, transporte de energia. Os autores recomendam abordagens ab initio, como a TDDFT em temperatura finita, sempre que forem computacionalmente viáveis.

Da folha de alumínio ao interior de Júpiter

A ligação com os planetas não é direta, e convém dizer isso antes de qualquer outra coisa. Um ponto de alumínio a 50 gigapascais e 7 mil kelvin não é uma condição joviana: os experimentos de hidrogênio e hélio considerados representativos do interior de Júpiter operam entre 93 e 150 gigapascais e entre 4.700 e 10.200 kelvin. O que viaja da bancada para o planeta não é o estado, é o método — e, junto com ele, o viés do método.

A cadeia funciona assim: experimentos de compressão dinâmica produzem estados extremos; a leitura desses estados calibra equações de estado, que amarram pressão, densidade e temperatura de um material; as equações de estado alimentam modelos de estrutura interna; e esses modelos são invertidos contra observações — no caso de Júpiter, os momentos gravitacionais medidos pela sonda Juno. Um viés na primeira etapa não desaparece nas seguintes.

Júpiter e Saturno

O peso desse elo foi quantificado por Howard e colegas em 2023. Ao refazer a inversão dos dados da Juno com equações de estado de hidrogênio e hélio diferentes, e ao permitir modificações controladas nessas equações, o grupo obteve tanto soluções com núcleos diluídos extensos quanto — ao contrário de outros modelos recentes — soluções com núcleos diluídos pequenos, estendendo-se a 20% da massa do planeta. A conclusão dos autores é explícita: as equações de estado usadas nos modelos de Júpiter têm efeito crucial sobre a estrutura e a composição inferidas.

Não é um detalhe de segunda ordem. Um núcleo diluído que alcança 20% da massa e um núcleo diluído que se espalha pelo planeta implicam histórias de formação diferentes. E há um segundo problema no mesmo regime. Em 2021, Brygoo e colegas comprimiram com laser amostras de hidrogênio e hélio previamente apertadas em células de bigornas de diamante e encontraram uma região de imiscibilidade — hélio que se separa do hidrogênio — encerrando-se acima de 150 gigapascais a 10.200 kelvin, com uma mudança mais sutil acima de 93 gigapascais a 4.700 kelvin. É a base experimental da hipótese da chuva de hélio no interior de Júpiter e de Saturno.

Os sub-Netunos

Nos exoplanetas do tipo sub-Netuno — mais massivos que a Terra e menores que Netuno, uma faixa de tamanho que o Sistema Solar não possui — o problema fica ainda mais agudo, porque nem a divisão entre núcleo rochoso e envelope de hidrogênio resiste. Rogers, Young e Schlichting publicaram em novembro de 2025 o primeiro modelo evolutivo de interior de sub-Netuno que incorpora a miscibilidade entre magma de silicato e hidrogênio, miscíveis acima de cerca de 4.000 kelvin nas pressões relevantes. Segundo o trabalho, sub-Netunos jovens podem guardar várias dezenas de por cento de sua massa de hidrogênio dentro do próprio interior, hidrogênio que vai sendo exsolvido conforme o planeta esfria, com a maior parte já fora depois de alguns bilhões de anos.

Um interior assim é, por definição, matéria quente e densa: uma mistura fortemente acoplada e parcialmente degenerada, cuja equação de estado é justamente o que se quer descobrir. Modelar planetas dessa classe depende de acertar exatamente o regime em que os modelos de gás uniforme acabaram de falhar. O mesmo vale para o outro grande consumidor de equações de estado nessa faixa: a fusão por confinamento inercial, em que uma cápsula de combustível é comprimida a estados quentes e densos antes de queimar. Em 5 de dezembro de 2022, o National Ignition Facility obteve pela primeira vez ganho maior que 1 — 2,05 megajoules de luz laser produziram 3,1 megajoules de fusão —, e cada passo desse programa depende de saber em que estado o combustível de fato está.

Duas rotas para o mesmo estado, e uma delas passa por Campinas

Há duas maneiras de levar matéria a esse regime, e elas têm defeitos complementares. Na compressão estática, com bigornas de diamante, a amostra pode permanecer sob pressão por dias, o que permite medidas longas e repetidas; a dificuldade é aquecê-la bastante sem perdê-la. Na compressão dinâmica, um laser gera uma onda de choque que alcança pressão e temperatura altas ao mesmo tempo; o preço é que o estado dura nanossegundos e cada disparo destrói o alvo.

As duas rotas não competem, se completam — e às vezes se somam no mesmo aparelho, como no experimento de hidrogênio e hélio descrito acima, que só alcançou seu canto do diagrama de fases porque combinou bigorna e laser.

O Brasil não aparece entre os autores do estudo do European XFEL — as 21 instituições envolvidas são europeias, norte-americanas e japonesas. Mas o país tem infraestrutura instalada na rota estática. A linha de luz EMA do Sirius, o síncrotron do Laboratório Nacional de Luz Síncrotron, em Campinas, combina difração, espectroscopia de absorção e espalhamento inelástico de raios X em amostras submetidas a extremos de temperatura, pressão e campo magnético, com fluxo de até 10¹⁴ fótons por segundo e foco de até 0,1 por 0,1 micrômetro. A própria descrição da linha coloca o objetivo nesses termos:

simular o interior de planetas e entender as propriedades da matéria que os compõe

São perguntas iguais com ferramentas diferentes. O que o experimento alemão acrescenta é uma advertência que vale para as duas rotas: em matéria quente e densa, o modelo usado para interpretar a medida entra no resultado e precisa ser testado junto com ela.

O que o teste não fecha

Um benchmark também precisa declarar o tamanho do que cobre, e este cobre pouco — de propósito, porque é a restrição que o torna utilizável como referência.

• Um material e um ponto. Alumínio, a cerca de 50 gigapascais e 0,6 elétron-volt. Nada garante que a mesma hierarquia entre os modelos se mantenha em hidrogênio, em berílio ou a temperaturas dez vezes maiores.
• O estado não é um ponto, é uma janela. A densidade ficou entre 3,75 e 4,5 gramas por centímetro cúbico, e a comparação com a teoria foi feita com média sobre essa faixa. Parte da incerteza citada vem dessa média, não do instrumento.
• A TDDFT também tem escolhas embutidas. Sua exatidão depende do funcional de troca e correlação adotado e das aproximações feitas no núcleo dinâmico. Concordar com o dado aqui não é o mesmo que estar certa em toda parte.
• O custo continua sendo o obstáculo. Cálculos ab initio não cabem na análise de rotina de cada disparo, e a recomendação dos autores é usá-los onde for computacionalmente viável — o que deixa em aberto o que fazer quando não for.
• O estado é tratado como congelado. A hipótese se sustenta porque o pulso de 25 femtossegundos é muito mais curto que a evolução do sistema, mas continua sendo uma hipótese, e limita o que se pode afirmar sobre a dinâmica mais lenta.

O que segue em aberto não é uma lei física ausente. É uma fronteira: até onde a aproximação barata permanece aceitável e a partir de que ponto ela começa a enganar. O experimento mostrou que, com desordem induzida por choque em alumínio, esse ponto já foi ultrapassado. Onde exatamente ele fica — em qual elemento, a qual temperatura, com quanto de ordem residual — é o que os próximos disparos, em outros materiais e outras condições, terão de mapear.