Об одном алгоритме расчета движений молекул двухатомных газов

Рассматривается проблема моделирования свойств двухатомных газов методамимолекулярной динамики. Подобные исследования являются традиционными дляфизики вещества. В настоящее время наблюдается повышенный интерес к даннойпроблеме в связи с развитием нанотехнологий и их внедрением в различные отрасли промышленности. В предлагаемой работе рассмотрен вопрос об уточнении исходной классической модели динамики Ньютона. В частности, обсуждается методика учета дополнительных степеней свободы, характеризующих вращательныедвижения двухатомных молекул. Это связано с необходимостью корректного вычисления теплоемкости. Для решения данной задачи предложено добавить в модель молекулярной динамики уравнения для моментов импульса и вращательныхскоростей молекул. Для такой расширенной постановки разработан специальныйчисленный алгоритм, обобщающий схему Верле. На основе предложенного алгоритма разработана расчетная программа. С ее помощью вычислена кривая теплоемкости для азота в диапазоне температур 100-400 К при давлении 1 атм. Полученные расчетные данные согласуются с известными данными из справочников.