Сравнение потенциальных функций для молекулярно-динамического моделирования сорбции метана в силикалите

С целью оптимального выбора, в работе исследованы два вида силового поля для оценки возможности их применения к моделированию процессов поглощения и переноса в системе силикалит (синтетический цеолит) – метан; для этого использован метод молекулярной динамики, реализованный в пакете LAMMPS. Приоритетными критериями выбора были сохранение целостности структур при длительных расчетах, стабильность системы, возможность моделирования отдельного кластера цеолита большого размера без использования периодических граничных условий для размножения решетки в пространстве. Для потенциала, удовлетворившего этим требованиям, проведен расчет структурных, термодинамических и транспортных характеристик системы. Рассчитаны радиальная и угловая функции распределения частиц для различных пар и троек атомов для чистого цеолита и для двухкомпонетной системы силикалит – метан. Достигнуто хорошее согласие всех полученных характеристик с литературными данными численных и натурных экспериментов.