Аннотация:Оригинальная методика атомистического компьютерного моделирования твердых растворов замещения, разработанная в [1] была применена для изучения локальной структуры в изоморфной системе гранатов гроссуляр Са3Al2(SiO4)3 - уваровит Са3Cr2(SiO4)3. Параметры межатомных потенциалов были оптимизированы с использованием эксперимента данных о структурных, упругих и термодинамических характеристик этих минералов. Расчеты локальной структуры для семи различных составов твердого раствора проводились по программе GULP 3.4.9 на суперкомпьютере СКИФ «Чебышев» МГУ в сверхячейке 2×2×4 структурного типа граната, содержащей 2560 атомов. Внутри такой ячейки были сконструированы разупорядоченные конфигурации с соотношением Al:Cr = 1:7, 1:5, 1:3, 1:1, 3:1, 5:1 и 7:1, имитирующие свойства статистического распределения атомов по позициям в макроскопическом кристалле. Были построены гистограммы межатомных расстояний М-O и O-O в сверхячейке и оценены значения параметров релаксации октаэдров CrO6 и AlO6.
График изменений расстояний M-O (M=Al, Cr) в зависимости от состава твердого раствора Ca3(Al, Cr)2(SiO4)3 в сравнении с экспериментальными данными, представлен на рис. 1. Промежуточные точки определялись усреднением всех рассчитанных расстояний в сверхячейке после процедуры энергетической минимизации, а крайние точки – расчетом методом Мотта-Литтлтона в приближении изолированного дефекта. Как видно из рисунка, параметры релаксации CrO6 и AlO6 октаэдров, равны 0.86 и 0.83, соответственно, и практически не зависят от состава твердого раствора, что хорошо согласуется с имеющейся экспериментальной информацией, полученной из анализа спектроскопических данных