The origin of complex behavior of liquid carbon: an insight from computer simulation

Título: A Origem do Comportamento Complexo do Carbono Líquido: Uma Perspectiva a partir de Simulação Computacional

Autores: Yu. D. Fomin
Instituição: Instituto Vereshchagin de Física de Alta Pressão, Academia Russa de Ciências.

---

1. O Problema

A fusão do grafite permanece um dos problemas mais controversos e de longa data na física de altas temperaturas. Embora experimentos tenham sido realizados desde o início do século XX, existe uma dispersão significativa e insatisfatória nos valores reportados para o ponto de fusão do grafite, variando de 4000 K a 6700 K.

As principais dificuldades identificadas no artigo incluem:

• Incerteza nos Pontos Críticos e Triplos: As estimativas existentes para o ponto crítico (transição líquido-gás) e o ponto triplo (sólido-líquido-gás) do carbono baseiam-se em suposições arbitrárias (ex: salto de volume de 20% na fusão) e dados experimentais inconsistentes.
• Falhas nas Suposições Termodinâmicas: Muitos estudos assumem que a capacidade calorífica e o coeficiente de expansão térmica do carbono líquido são constantes ou semelhantes aos do grafite sólido.
• Ignorância da Fase Gasosa: A maioria das análises experimentais não leva em conta a influência da fase gasosa (sublimação) e das flutuações críticas nas medições de calorimetria e densidade.
• Contradições Recentes: Medidas recentes sugerem capacidades caloríficas altas e comportamentos que indicam que o "carbono líquido" pode estar, na verdade, em um estado supercrítico, onde as propriedades termodinâmicas variam drasticamente.

2. Metodologia

O autor utilizou simulações de Dinâmica Molecular (DM) para investigar as propriedades termodinâmicas do carbono líquido e localizar seu ponto crítico.

• Potencial de Interação: Foi utilizado o potencial de aprendizado de máquina GAP-20 (Gaussian Approximation Potential), ajustado para reproduzir fases condensadas do carbono.
• Configuração do Sistema: Simulação de um sistema de 4096 átomos de carbono em uma caixa cúbica com condições de contorno periódicas.
• Conjunto Termodinâmico: Ensemble Canônico (NVT).
• Parâmetros de Simulação:
  • Passo de tempo: 2,0 fs.
  • Período de equilíbrio: 10 ps.
  • Período de produção: 20 ps.
• Varredura de Estados: Simulações realizadas em uma vasta gama de densidades (0,05 a 2,0 g/ml) e temperaturas (4000 K a 7000 K, com foco em 4000-5000 K).
• Análise de Dados: Os dados de pressão e energia interna foram ajustados a uma equação polinomial para calcular derivadas termodinâmicas, incluindo capacidade calorífica isobárica ($c_P$), coeficiente de expansão térmica ($\alpha_P$) e compressibilidade isotérmica ($\beta_T$). O ponto crítico foi identificado como o ponto de inflexão nas isotermas.

3. Contribuições Principais e Resultados

A. Localização do Ponto Crítico

A simulação com o potencial GAP-20 resultou em um ponto crítico significativamente mais baixo do que as estimativas teóricas anteriores (como a de Leider et al., que previa $T_c \approx 6810$ K):

• Temperatura Crítica ($T_c$): 4235 K
• Densidade Crítica ($\rho_c$): 0,1887 g/ml
• Pressão Crítica ($P_c$): 106,61 bar

Este valor de $T_c$ (4235 K) é notavelmente próximo às estimativas de fusão do grafite, sugerindo que a fusão ocorre em condições onde o sistema está muito próximo ou já no regime supercrítico.

B. Linhas de Widom e Flutuações

O estudo mapeou as Linhas de Widom (linhas de máximo das funções de resposta termodinâmica no regime supercrítico):

• Observou-se que a capacidade calorífica isobárica ($c_P$) e o coeficiente de expansão térmica ($\alpha_P$) atingem máximos acentuados próximos ao ponto crítico.
• A $c_P$ pode atingir valores de até $25 k_B$ por partícula próximo ao crítico, caindo rapidamente à medida que se afasta.
• A compressibilidade isotérmica ($\beta_T$) também exibe um máximo, mas desaparece mais rapidamente com o aumento da temperatura.

C. Cristalização Espontânea

Durante as simulações, observou-se cristalização espontânea do carbono líquido em grafite a 4100 K e 4200 K, indicando que a linha de fusão está muito próxima do ponto crítico no plano densidade-temperatura.

D. Explicação para a Dispersão Experimental

O trabalho propõe que a grande variação nos resultados experimentais do ponto de fusão do grafite deve-se ao fato de que, durante a fusão, o sistema entra em um estado supercrítico ou próximo a ele.

• Flutuações Termodinâmicas: Perto do ponto crítico, pequenas incertezas experimentais em pressão ou temperatura levam a grandes erros na determinação de densidade e capacidade calorífica.
• Falha de Modelos Simplificados: As suposições comuns de capacidade calorífica constante ou expansão térmica linear falham completamente nesta região, levando a interpretações errôneas dos dados de calorimetria e pirometria.

4. Significado e Limitações

Significado:
O artigo oferece uma explicação física plausível para a inconsistência histórica nos dados de fusão do grafite. Ele sugere que a proximidade entre a linha de fusão e o ponto crítico líquido-gás gera flutuações termodinâmicas extremas que confundem as medições experimentais. O trabalho destaca a necessidade de considerar fenômenos críticos e a presença de fase gasosa (supercrítica) ao analisar dados de altas temperaturas.

Limitações e Cuidados:

• Validade do Potencial GAP-20: O potencial foi ajustado para fases condensadas (sólido/líquido denso) e pode não ser preciso o suficiente para descrever a fase gasosa diluída ou o ponto crítico real.
• Deslocamento de Temperatura: O autor admite que o $T_c$ calculado (4235 K) pode ser subestimado em relação ao valor real do carbono. É provável que o ponto crítico real esteja em temperaturas mais altas (talvez próximo aos 6000 K), mas a qualidade qualitativa do comportamento (presença de máximos termodinâmicos e flutuações) permanece válida.
• Tamanho do Sistema e Tempo: Devido ao custo computacional do GAP-20, o tamanho do sistema (4096 átomos) e o tempo de simulação (20 ps) podem ser insuficientes para capturar totalmente as flutuações críticas de longo alcance, embora sejam suficientes para identificar tendências qualitativas.

Conclusão

O estudo conclui que a complexidade do comportamento do carbono líquido e a dificuldade em determinar seu ponto de fusão exato estão intrinsecamente ligadas à proximidade da fusão com o ponto crítico líquido-gás. Novos modelos de potencial, ajustados simultaneamente para fases gasosas, líquidas e sólidas, são necessários para simulações mais precisas, mas a ideia central de que flutuações críticas distorcem as medições experimentais é um avanço significativo na compreensão do problema.