Получение изотермических характеристик, параметров уравнения состояния для PETN методами реакционной молекулярной динамики и темодинамики

Методом молекулярной динамики (МД) в программном пакете LAMMPS (Large-scale Atomic/Molecular Massively Parallel Simulator) с использованием реакционного силового поля ReaxFF-lg были рассчитаны изотермы нереагирующего вещества PETN в диапазоне давлений до 30 ГПа. Были получены значения коэффициента модуля всестороннего сжатия $K_0=\mathrm{9,6}$ ГПа и производной модуля сжатия $K_0$ по давлению $K_0^{\text{'}}=\mathrm{8,0}$, которые можно использовать как параметры термического уравнения Берча–Мурнагана 3-го порядка. Были подобраны коэффициенты уравнения состояния (УрС) в форме Ми–Грюнайзена, часто используемого для моделирования теплофизических свойств вещества, в том числе при статическом и ударно-волновом сжатии. Для нахождения коэффициентов УрС применялся метод построения изохорно-изотермического потенциала твердых веществ в форме квазигармонического приближения Эйнштейна. Верификация полученных результатов показала хорошее согласие с экспериментальными данными в широком диапазоне изменения давления и температуры, в том числе вдоль ударной адиабаты.