Проблема установления взаимосвязи между химическим составом, атомной структурой и физическими свойствами кристаллов до сих пор является открытой. Квантомеханический подход к решению данной проблемы сопряжен с большими вычислительными затратами. В связи с этим, большую пользу в исследовании проблемы «структура-свойство» несет выявление общих структурных характеристик, влияющих на физические свойства кристаллов. Одной из таких общих характеристик является симметрия.
Симметрия кристалла определяет набор физических свойств, которыми может обладать кристалл. В силу того, что кристаллическая решетка периодична в пространстве, характерными операциями для пространственной симметрии являются инверсия, повороты (на 180°, 90°, 60° и 120°), отражение в плоскости, инверсионные повороты и параллельный перенос, а также различные комбинации данных преобразований. Инвариантность атомной структуры кристалла относительно одной из 230-ти пространственных групп симметрии называется федоровской симметрией.
Однако для некоторых кристаллов характерно, что большая часть их атомной структуры инвариантна относительно надгруппы группы симметрии кристалла как целого. Такое свойство кристаллов называют псевдосимметрией, а кристаллы называют псевдосимметричными.
Псевдосимметрия влияет на многие физические свойства кристаллов: генерация второй оптической гармоники, пьезоэффект, пироэффект, сегнетоэлектрический эфффект. Псевдосимметричное расположения атомов может приводить к закономерному погасанию отдельных групп рефлексов на дифракционных картинах от псевдосимметричных кристаллов. Наличие псевдосимметрии в кристалле является признаком возможного фазового перехода второго рода и т. д.
Рис. 1. Генерация второй оптической гармоники
Рис. 2. Атомный кластер
Наибольший интерес в настоящее время представляют атомные кластеры разупорядоченных кристаллов. Разупорядоченные кристаллы характеризуются наличием одной или нескольких смешанных кристаллографических позиций или позиций с дефицитной заселенностью – вакансиями. Разупорядоченность влияет на многие физические свойства кристаллов, в частности на оптические.
Для исследования разупорядоченных атомных кластеров реализован набор алгоритмов, позволяющих вычислять степень инвариантности функции электронной плотности атомного кластера, варьируя размер кластера по всем направлениям, заселенность разупорядоченных позиций и количество повторных вычислений. При моделировании кластера разупорядоченная позиция заселяется методом Монте-Карло. Каждый из алгоритмов представляет собой набор независимых экспериментов, отличающихся друг от друга варьируемыми параметрами. Распараллеливание данных алгоритмов реализовано для вычислительных систем с разделенной памятью, используя технологию MPI, где каждый из параллельных процессов выполняет свой набор экспериментов.
На основе предложенных параллельных алгоритмов реализован программный комплекс, позволяющий исследовать псевдосимметрию атомных структур кристаллов и атомных кластеров. Метод вычисления степени инвариантности функции электронной плотности имеет параллельные реализации для центрального процессора и графического ускорителя. Для исследования разупорядоченных атомных кластеров реализован ряд параллельных алгоритмов количественной оценки псевдосимметрии кристалла, варьируя размер атомного кластера и заселенность разупорядоченных позиций.
На рис. 3 представлен график зависимости времени вычисления функционала (2) от количества атомов и объема элементарной ячейки. Серия расчетов производилась на модельных кристаллах, в которых количество атомов и объем были изменены вручную.
Рис. 3. Время работы различных реализаций алгоритма в зависимости от количества атомов и объема ячейки.
Рис. 4. Ускорение алгоритма исследования псевдосимметрии атомных кластеров.