Аннотация. В статье обсуждается влияние замены ионов кальция Ca²⁺ в определенных позициях гидроксиапатита на ионы железа Fe²⁺ и Fe³⁺. Показано, как методом моделирования молекулярной динамики с применением косоугольной ячейки в периодических граничных условиях найти плотность гидроксиапатита. Полученные значения плотности хорошо согласуются с известными данными и результатами расчетом методами теории функционала плотности. Также описан процесс выпаривания ионов, плавления и сублимации разных образцов гидроксиапатита. Результаты молекулярно-динамического моделирования показывают заметные различия в температурах структурных перестроек для гидроксиапатита с замещениями ионами Fe²⁺ и Fe³⁺. Так, образец гидроксиапатита с ионами Fe²⁺ испытывает значительные структурные перестройки уже в диапазоне температур 727–827 °С, что согласуется с экспериментальными данными по распаду катион-замещенного гидроксиапатита с ионами железа при температурах около 800 °С.

 

Ключевые слова: метод молекулярной динамики, алгоритм поиска среднеквадратичных отклонений атомов, гидроксиапатит