Avaliação de diferentes potenciais interatômicos no cálculo do tensor de elasticidade do tungstato de zircônio
Abstract
O Tungstato de Zircônio (ZrW2O8) é um material que exibe Expansão Térmica Negativa (ETN), isotrópica em um amplo intervalo de temperatura (0,3 a 1050 K). Apesar de amplamente estudado, existem controvérsias acerca dos mecanismos microscópicos responsáveis por este comportamento. A fase cúbica deste composto, denominada a-ZrW2O8, já foi motivo de estudo através de simulações computacionais utilizando Potenciais Interatômicos (PI) e Dinâmica de Rede na Aproximação Quasi-Harmônica (DRQH). Nos dois PI distintos propostos na literatura conseguiu-se reproduzir a ETN da a-ZrW2O8, mas não a dependência com a temperatura do tensor de elasticidade. É partindo desta observação que este trabalho pretende avaliar o desempenho de PI existentes e de novos PI em simulações computacionais visando a descrição da dependência com a temperatura do tensor de elasticidade da a-ZrW2O8 entre 0 e 300 K. Utilizaram-se dados experimentais, tais como posições atômicas, parâmetros de rede e o tensor de elasticidade da a-ZrW2O8 em temperaturas entre 0 e 300 K e, em outra série de cálculos, a hipersuperfície de energia ab initio no limite atérmico para obter os parâmetros dos PI. Diferentes estratégias foram empregadas na busca pelos parâmetros dos PI incluindo minimização em linha, Newton-Raphson/BFGS e Algoritmo Genético (AG). Concluiu-se que não é possível descrever as propriedades estruturais e elásticas da a-ZrW2O8 em função da temperatura com PI simples e que esta incapacidade não está relacionada a qualquer limitação da DRQH ou dos parâmetros dos PI, mas à forma analítica dos PI empregados. Isto sugere que se deve ter cautela na interpretação de resultados obtidos com estes potencias já disponíveis na literatura. Como alternativas para a solução deste problema, pode-se considerar o uso de redes neurais para a representação da hipersuperfície de energia ab initio, o uso de PI mais sofisticados que levam em consideração a vizinhança atômica (bond order potentials) e, também, cálculo ab initio a T > 0, este último a um custo computacional muito mais elevado.