Pressure-Temperature Phase Diagram and $\lambda$-Transition in Liquid Sulfur

Utilizando simulações de dinâmica molecular com potenciais de aprendizado de máquina, o estudo investiga microscopicamente a transição $\lambda$ do enxofre, reconstruindo seu diagrama de fases e demonstrando que a polimerização é iniciada por anéis não-$S_8$ pouco antes do ponto de fusão.

O essencial

O Mistério do Enxofre: De Anéis de Diamante a Correntes de Ferro

Imagine que o enxofre é como um conjunto de peças de LEGO. Na maior parte do tempo, essas peças estão montadas em pequenos anéis perfeitos (como se fossem pequenas correntes de bicicleta fechadas). Mas, dependendo de quanto você "esquenta" ou "aperta" esse material, essas peças decidem mudar de forma de um jeito muito dramático.

Este estudo científico usou supercomputadores para observar esse "drama" em nível microscópico.

1. O "Baile de Máscaras" (A Transição $\lambda$)

Imagine uma festa de gala onde todos os convidados estão usando máscaras de anéis perfeitos. De repente, a temperatura sobe e a música fica muito agitada. Alguns convidados começam a tirar as máscaras, a perder o equilíbrio e, em vez de ficarem em seus pequenos grupos circulares, eles começam a dar as mãos e formar filas gigantescas que atravessam o salão.

Na ciência, chamamos isso de transição $\lambda$. O enxofre deixa de ser um líquido de "anéis" (pequenas moléculas isoladas) e se transforma em um líquido de "correntes" (polímeros longos). É como se o líquido passasse de uma água leve para algo mais viscoso, como um mel ou um xarope, porque essas correntes longas ficam se emaranhando umas nas outras.

2. O "Culpado" da Confusão (Os Anéis Imperfeitos)

Os cientistas queriam saber: como a primeira corrente começa? É como perguntar como uma briga começa em uma sala silenciosa.

Eles descobriram que o culpado não são os anéis perfeitos de 8 átomos (os $S_8$). O problema são os "anéis rebeldes" (com 5, 6 ou 7 átomos). Esses anéis são instáveis, como uma peça de LEGO mal encaixada. Eles se quebram facilmente e servem como o "primeiro dominó" a cair, desencadeando a formação das grandes correntes.

3. O "Esmagamento" (Pressão vs. Temperatura)

Agora, imagine que você está tentando fazer essa festa acontecer, mas alguém começa a apertar as paredes do salão (isso é a pressão).

• Com pouca pressão: O enxofre primeiro derrete (os anéis soltam-se do cristal) e só depois, com mais calor, eles viram correntes.
• Com muita pressão: As paredes apertam tanto que não há espaço para o "baile" acontecer de forma organizada. O enxofre já começa a formar as correntes dentro do próprio cristal, antes mesmo de derreter! É como se as pessoas começassem a dar as mãos ainda sentadas em suas cadeiras, antes mesmo de levantarem para dançar.

4. Como eles fizeram isso? (O "Simulador de Realidade")

Os cientistas não usaram tubos de ensaio comuns, mas sim um "Simulador de Realidade Virtual" ultra avançado (chamado de Machine Learning Interatomic Potential).

Em vez de tentar calcular cada átomo do zero (o que levaria mil anos), eles "treinaram" o computador com inteligência artificial para que ele aprendesse as regras do jogo: "se o átomo chegar perto assim, ele gruda; se esquentar assim, ele solta". Isso permitiu que eles vissem o enxofre se transformando em tempo real, como se estivessem assistindo a um filme em câmera lenta de um mundo invisível.

Resumo da Ópera:

O estudo criou um mapa (diagrama de fase) que mostra exatamente como o enxofre se comporta. Ele revela que o enxofre é um camaleão: ele pode ser um cristal organizado, um líquido de anéis leves ou um "emaranhado" de correntes longas, dependendo de quão quente ou quão apertado o ambiente estiver.

Resumo técnico

Resumo Técnico: Diagrama de Fases Pressão-Temperatura e Transição $\lambda$ no Enxofre Líquido

Título Original: Pressure–Temperature Phase Diagram and $\lambda$-Transition in Liquid Sulfur
Autores: Sonia Salomoni, Frédéric Datchi, A. Marco Saitta e Arthur France-Lanord.

---

1. O Problema (Contexto e Motivação)

O enxofre elemental possui um diagrama de fases extremamente rico devido à sua diversidade molecular, variando de anéis cíclicos (como o $S_8$) a cadeias poliméricas lineares. Um fenômeno central é a transição $\lambda$, uma transição induzida por temperatura no estado líquido onde os anéis de $S_8$ se abrem e se polimerizam em cadeias longas.

Apesar de ser bem documentada experimentalmente (mudanças bruscas de viscosidade, cor e capacidade térmica), a transição $\lambda$ apresenta desafios teóricos e computacionais significativos:

• Escala de Tempo e Tamanho: Simulações de Dinâmica Molecular (MD) baseadas em Teoria do Funcional da Densidade (DFT) são computacionalmente proibitivas para capturar a polimerização, que exige sistemas grandes e tempos de equilíbrio longos.
• Complexidade de Ligação: Potenciais empíricos tradicionais têm dificuldade em descrever com precisão a quebra e formação de ligações químicas.
• Lacunas de Conhecimento: Não se conhecia, em nível atomístico, a sequência completa de transformações desde o cristal $\alpha$-$S_8$ até o líquido polimérico, nem o mecanismo exato de como a polimerização ocorre quando a linha de fusão e a de polimerização se encontram em altas pressões.

2. Metodologia

Os autores utilizaram uma abordagem de ponta combinando aprendizado de máquina com física de primeira ordem:

• Potencial Interatômico Aprendido por Máquina (MLIP): Desenvolveram um modelo Allegro (um modelo de passagem de mensagens equivariante local) treinado em um conjunto de dados de DFT (nível GGA-PBE).
• Aprendizado Ativo (Active Learning): Utilizaram a estratégia Query by Committee (QBC) para selecionar as configurações de treinamento mais informativas, garantindo que o potencial cobrisse tanto fases moleculares quanto poliméricas.
• Simulações de Dinâmica Molecular (MD): Realizaram simulações de grande escala utilizando o software LAMMPS, empregando "rampas de temperatura" (T-ramps) para observar as transições de fase de forma espontânea sob diferentes pressões e taxas de aquecimento.
• Análise Estrutural: Utilizaram descritores de estrutura local (como o vetor SOAP - Smooth Overlap of Atomic Position) e análise de coordenação atômica para rastrear a evolução das espécies moleculares.

3. Principais Contribuições e Resultados

O estudo fornece uma visão microscópica inédita da transição $\lambda$ e do diagrama de fases do enxofre:

• Mecanismo da Transição $\lambda$ (Baixa Pressão): Identificaram que a polimerização não começa apenas pela abertura direta de anéis $S_8$, mas é desencadeada pela formação de anéis não-$S_8$ (como $S_7$). Esses anéis menores atuam como "centros reativos" ou sementes que facilitam a abertura de cadeias e a subsequente catenação.
• Dependência da Taxa de Aquecimento: As simulações reproduziram com sucesso a dependência da temperatura de transição ($T_\lambda$) em relação à taxa de aquecimento, um comportamento crítico observado experimentalmente. O modelo também reproduziu o pico característico de capacidade térmica ($C_s$) na forma de $\lambda$.
• Regime de Alta Pressão (Fusão e Polimerização Coincidentes): Em pressões intermediárias (acima de ~0,1 GPa), o estudo revelou que a polimerização começa dentro da fase cristalina, ligeiramente antes da fusão. Observaram-se anéis não-$S_8$ e cadeias poliméricas que ainda mantêm a conformação da rede cristalina (polimerização topoquímica). Quando a temperatura aumenta o suficiente, a fusão ocorre a partir dessas porções poliméricas já existentes.
• Reconstrução do Diagrama de Fases: O trabalho reconstrói o diagrama P-T, mostrando que a temperatura de polimerização diminui moderadamente com a pressão até encontrar a linha de fusão em um ponto crítico, após o qual o cristal derrete diretamente em um líquido polimérico.

4. Significância

Este trabalho é altamente significativo por:

1. Resolver um problema clássico: Fornecer a primeira descrição atomística completa da sequência de transições de fase do enxofre.
2. Validar novas tecnologias: Demonstrar o poder de potenciais interatômicos aprendidos por máquina (MLIPs) para simular processos químicos complexos (quebra/formação de ligações) em escalas de tempo e tamanho antes impossíveis para a DFT pura.
3. Avanço Teórico: Desafiar modelos fenomenológicos antigos (como o de Tobolsky-Eisenberg) ao mostrar o papel crucial das espécies moleculares intermediárias (anéis não-$S_8$) no mecanismo de polimerização.