Исследователи из МФТИ и Объединенного института высоких температур РАН определили, при каких условиях квазиклассические модели водородной плазмы начинают давать нефизичные результаты. Ученые показали, что при достаточно сильном взаимодействии частиц такие модели перестают правильно учитывать квантовые эффекты, связанные с принципом Паули. В результате в расчетах возникают нефизичные электрон-ионные кластеры, а моделируемая система теряет устойчивость. Работа опубликована в журнале Physics of Plasmas. Исследование поддержано Российским научным фондом (грант №24-19-00746).
Поведение водородной плазмы моделируют при изучении самых разных систем — от теплого плотного вещества до разреженной и ультрахолодной плазмы. Полностью квантовые расчеты позволяют детально учитывать свойства электронов, однако требуют очень больших вычислительных ресурсов. Поэтому исследователи часто используют классическую молекулярную динамику с квантовыми псевдопотенциалами. В таких расчетах протоны и электроны движутся согласно классическим уравнениям, а квантовые эффекты приближенно учитываются с помощью полученных из решения уравнений квантовой механики псевдопотенциалов. Это позволяет скорректировать взаимодействие частиц на малых расстояниях и значительно ускорить моделирование. Однако область применимости такого приближения ограничена.
Ученые исследовали стандартный и улучшенный варианты псевдопотенциала Кельбга в широком диапазоне температур и плотностей для классической водородной плазмы, в которой частицы с высокой точностью можно считать точечными. Результаты сравнивались с эталонными квантовыми расчетами методом Монте-Карло с интегралами по траекториям. При слабом и умеренном взаимодействии частиц оба псевдопотенциала давали физически разумные результаты, причем улучшенный вариант хорошо согласовывался с квантовыми расчетами. Однако по мере усиления взаимодействия в системе начинали возникать заметные артефакты.
Силу взаимодействия характеризует безразмерный параметр неидеальности Γ: чем он больше, тем сильнее взаимодействие частиц по сравнению с их тепловым движением. Исследование показало, что при Γ выше примерно 0,5 в расчетах появляются нефизичные особенности, а при Γ ≥ 1 модель становится неустойчивой. Электроны с одинаковой проекцией спина начинают чрезмерно сближаться и образовывать кластеры, которых в реальной плазме возникать не должно.
«Нам было важно не просто очертить границу применимости известных псевдопотенциалов, а понять, что именно происходит при ее переходе. Мы показали, что проблема начинается там, где парное описание перестает воспроизводить обменные квантовые корреляции между электронами»,— говорит Александр Онегин, аспирант МФТИ и младший научный сотрудник ОИВТ РАН.
Причина неустойчивости связана с принципом Паули. Этот фундаментальный закон квантовой механики запрещает электронам с одинаковой проекцией спина занимать одно и то же квантовое состояние. Благодаря этому такие электроны в определенном смысле избегают друг друга. В квазиклассических моделях такое поведение пытаются воспроизвести с помощью спин-зависимых парных псевдопотенциалов.
Однако при сильном взаимодействии частиц одного парного описания оказывается недостаточно: оно уже не может правильно передать многочастичные обменные корреляции. Нарушение проявляется сразу в нескольких характеристиках моделируемой системы. В радиальных функциях распределения возникают аномальные пики, указывающие на чрезмерное сближение частиц. Одновременно энергия системы неограниченно уменьшается с течением времени моделирования. Такое поведение связано не с установлением равновесия, а с развитием нефизической неустойчивости.
Полученные результаты важны как для теории плазмы, так и для развития вычислительных методов. Они показывают, в какой области квазиклассические модели с псевдопотенциалами дают надежное описание, а при каких условиях необходимо использовать более строгие квантовые подходы. Кроме того, работа указывает, каких именно физических эффектов не хватает парным моделям при переходе к сильному взаимодействию частиц.
«Квазиклассические модели по-прежнему востребованы, поскольку позволяют исследовать сложные плазменные системы без чрезмерных вычислительных затрат. Наш результат задает практическую границу применимости таких подходов: показывает, где они еще обеспечивают физически надежное описание, а где уже нет. При этом пока не существует основанного на первых принципах способа надежно расширить эту область применимости»,— отмечает Павел Левашов, доктор физико-математических наук, заведующий кафедрой физики высокотемпературных процессов МФТИ.
Исследование поддержано Российским научным фондом (грант №24-19-00746).
Научная статья: Onegin A. S., Demyanov G. S., Levashov P. R. The role of Pauli principle in simulations of classical plasma // Physics of Plasmas, 2026; 33, 042706. DOI: https://doi.org/10.1063/5.0315414
