Dynamic Slowdown and Spatial Correlations in Viscous Silica Melt: Perspectives from Dynamic Disorder

Resumo Técnico: Desaceleração Dinâmica e Correlações Espaciais em Fundido Viscoso de Sílica

Declaração do Problema
A origem microscópica da dramática desaceleração dinâmica em líquidos formadores de vidro permanece um desafio central na ciência da matéria condensada. Embora medidas estruturais (por exemplo, funções de distribuição radial) mostrem mudança mínima ao resfriar, os tempos de relaxação e as viscosidades aumentam em ordens de grandeza. Esse desacoplamento é particularmente complexo em formadores de vidro fortes como a sílica amorfa (SiO₂), que exibe comportamento quase Arrhenius, mas apresenta dinâmicas de relaxação distintas e dependentes da espécie entre átomos de silício e oxigênio. Estudos anteriores identificaram eventos de "salto" como rearranjos estruturais elementares, mas as restrições estruturais específicas e a natureza da desordem dinâmica que governam esses saltos na sílica — particularmente como diferem de outros líquidos tetraédricos como a água — permanecem pouco compreendidas.

Metodologia
Os autores empregaram simulações extensivas de dinâmica molecular (DM) utilizando dois potenciais de interação distintos: o modelo Beest–Kramer–Santen (BKS) (incorporando eletrostática de longo alcance) e o modelo Coslovich–Pastore (CP) (um potencial de curto alcance simplificado). As simulações foram conduzidas em uma ampla faixa de temperatura (2500 K a 6100 K para BKS; unidades reduzidas de 0,29 a 0,60 para CP) usando o pacote LAMMPS.

Estruturas analíticas-chave incluíram:

1. Identificação de Saltos: Uma função de salto $h(t)$ foi definida para distinguir o movimento vibracional localizado de saltos intermitentes de ruptura de gaiola. Os limiares ($h^*$) foram determinados a partir de distribuições cumulativas para separar estados de gaiola e de salto.
2. Análise de Desordem Dinâmica: Os autores analisaram distribuições de tempo de residência e probabilidades de sobrevivência ($C_S(t)$) para quantificar desvios da estatística de Poisson. Eles empregaram um quadro de taxa dependente do tempo para distinguir entre limites de flutuação rápida (semelhante a Poisson) e flutuação lenta (desordem dinâmica).
3. Correladores Estruturais: A divergência de Kullback-Leibler (KL) foi usada para rastrear a evolução de variáveis estruturais locais (distâncias interatômicas e distribuições angulares) à medida que as partículas se aproximavam do limiar de salto.
4. Identificação de Variáveis Lentas: As probabilidades de sobrevivência foram recalculadas no "limite de flutuação lenta" tratando distâncias específicas de vizinhos (por exemplo, $r_{SiO4}$, $r_{SiSi4}$, $r_{OSi2}$) como variáveis lentas para determinar quais flutuações estruturais modulam mais fortemente as taxas de salto.
5. Correlações Ponto-para-Conjunto (PTS): Comprimentos de correlação estática PTS ($\xi_{PTS}$) foram computados usando simulações de cavidade cúbica com fronteiras congeladas para quantificar a extensão espacial de rearranjos cooperativos. Funções de sobreposição resolvidas por espécie foram analisadas para distinguir a cooperatividade do silício e do oxigênio.
6. Simulações Confinadas: Para verificar a causalidade, vizinhos específicos identificados como variáveis lentas foram seletivamente congelados em simulações de cavidade para observar a supressão da relaxação local.

Principais Resultados

• Dinâmica Não-Poissoniana e Desordem Dinâmica: À medida que a temperatura diminui, a dinâmica de salto desvia-se da estatística de Poisson, evidenciada pelo aumento dos parâmetros de aleatoriedade ($R$) e probabilidades de sobrevivência esticadas. Isso indica o surgimento de desordem dinâmica onde variáveis estruturais que evoluem lentamente modulam as taxas de salto.
• Variáveis Lentas Dependentes da Espécie:
  • Silício: Os saltos de silício são governados principalmente pelo deslocamento do quarto vizinho de oxigênio ($r_{SiO4}$). Nas temperaturas mais baixas, o quarto vizinho de silício ($r_{SiSi4}$) também se torna uma variável lenta relevante, indicando um rearranjo cooperativo de dimensão superior.
  • Oxigênio: Os saltos de oxigênio são dominados pelo deslocamento do segundo vizinho de silício ($r_{OSi2}$). Diferentemente do silício, a dinâmica do oxigênio não requer variáveis lentas adicionais em baixas temperaturas; a influência de vizinhos mais distantes (por exemplo, sexto vizinho de oxigênio) é negligenciável para a modulação da taxa.
• Correlações Espaciais Assimétricas: A análise ponto-para-conjunto revela uma assimetria marcante nas escalas de comprimento cooperativo. O $\xi_{PTS}$ para o silício aumenta moderadamente com o resfriamento, estendendo-se além da quarta casca de oxigênio. Em contraste, o $\xi_{PTS}$ para o oxigênio satura em um valor constante (limitado à segunda casca de silício) e não mostra crescimento com o resfriamento adicional.
• Comparação com a Água: Diferentemente da água super-resfriada, onde a desordem dinâmica é impulsionada pelo deslocamento lento de átomos de oxigênio de quarto (e posteriormente quinto) vizinhos dentro de uma rede de ligações de hidrogênio, a sílica exibe mecanismos específicos da espécie. A rede covalente rígida Si-O impõe caminhos de relaxação distintos: o silício requer rearranjos coletivos da rede, enquanto o movimento do oxigênio é mais localizado.

Significância e Afirmações
O artigo afirma fornecer uma estrutura microscópica que liga desordem dinâmica, restrições dependentes da espécie e correlações cooperativas em formadores de vidro fortes. Os autores afirmam que:

1. A desaceleração dinâmica na sílica decorre de restrições estruturais específicas e dependentes da espécie, e não de um efeito de rede uniforme.
2. A relaxação na sílica não pode ser descrita por uma única escala de tempo característica ou por uma única coordenada de reação; em vez disso, envolve um espaço de dimensão superior de variáveis lentas que cresce com o resfriamento para o silício, mas permanece localizado para o oxigênio.
3. Os comportamentos distintos do silício e do oxigênio explicam o caráter geral "forte" (Arrhenius) da sílica, onde o comportamento localizado e do tipo "forte" do oxigênio domina a dinâmica total, mascarando a tendência ligeiramente mais do tipo "frágil" do silício.
4. O estudo estabelece uma conexão direta entre a origem microscópica da desordem dinâmica (flutuações de taxa impulsionadas por reorganizações específicas de vizinhos) e a ordem amorfa estática (medida por $\xi_{PTS}$), sugerindo que estas são manifestações complementares do processo de relaxação.

Os autores concluem que essas descobertas oferecem uma visão mais profunda sobre a desaceleração de redes formadoras de vidro fortes e destacam a necessidade de tratamentos teóricos resolvidos por espécie para líquidos de rede, distinguindo-os de modelos como a água, onde uma única classe de variáveis lentas domina.