Модели водородной плазмы теряют устойчивость при Γ ≥ 1

Исследователи Московского физико-технического института и Объединенного института высоких температур РАН выяснили, в каких условиях квазиклассические модели водородной плазмы перестают давать физически достоверные результаты. При сильном взаимодействии частиц такие расчеты неправильно воспроизводят квантовые эффекты, связанные с принципом Паули.

Из-за этого в моделируемой системе появляются несуществующие электрон-ионные кластеры. Одновременно плазма теряет устойчивость, хотя подобного поведения у реального вещества возникать не должно.

Водородную плазму моделируют при изучении широкого круга систем. К ним относятся теплое плотное вещество, а также разреженная и ультрахолодная плазма. Полностью квантовые методы позволяют подробно учитывать свойства электронов, но требуют значительных вычислительных ресурсов.

По этой причине ученые часто применяют классическую молекулярную динамику, дополненную квантовыми псевдопотенциалами. В таких расчетах протоны и электроны движутся по классическим уравнениям, а квантовые свойства приближенно задаются псевдопотенциалами, полученными при решении задач квантовой механики.

Такой подход корректирует взаимодействие частиц на малых расстояниях и заметно ускоряет вычисления. Однако использовать его можно не при любых параметрах плазмы: область применимости квазиклассического приближения ограничена.

Авторы проверили стандартную и улучшенную версии псевдопотенциала Кельбга. Расчеты провели для классической водородной плазмы в широком диапазоне температур и плотностей. В этой системе частицы с высокой точностью допустимо рассматривать как точечные.

Полученные данные сопоставили с эталонными квантовыми расчетами. Для этого применили метод Монте-Карло с интегралами по траекториям. При слабом и умеренном взаимодействии обе версии псевдопотенциала давали физически разумную картину, а улучшенный вариант хорошо совпадал с результатами квантового моделирования.

По мере усиления взаимодействия в расчетах начали проявляться артефакты. Их выраженность связана с безразмерным параметром неидеальности Γ. Чем выше его значение, тем сильнее взаимодействие частиц по сравнению с их тепловым движением.

При Γ выше примерно 0,5 в моделях появлялись нефизичные особенности. Когда параметр достигал значения 1 и превышал его, система становилась неустойчивой. Электроны с одинаковой проекцией спина начинали слишком близко сходиться и объединяться в кластеры, которые в реальной плазме возникнуть не могут.

«Нам было важно не просто очертить границу применимости известных псевдопотенциалов, а понять, что именно происходит при ее переходе. Мы показали, что проблема начинается там, где парное описание перестает воспроизводить обменные квантовые корреляции между электронами», — сообщил Александр Онегин, аспирант МФТИ и младший научный сотрудник ОИВТ РАН.

Источник неустойчивости связан с принципом Паули. Этот фундаментальный закон квантовой механики не позволяет электронам с одинаковой проекцией спина занимать одно и то же квантовое состояние. Поэтому такие частицы в определенном смысле избегают друг друга.

В квазиклассических расчетах это поведение пытаются передать с помощью парных псевдопотенциалов, зависящих от спина. Но при сильном взаимодействии одного парного описания уже недостаточно. Оно не может корректно воспроизвести многочастичные обменные корреляции между электронами.

Нарушение проявляется сразу в нескольких характеристиках плазмы. На радиальных функциях распределения возникают аномальные пики, которые указывают на чрезмерное сближение частиц. Кроме того, энергия системы без ограничений снижается по мере продолжения моделирования.

Такое уменьшение энергии не означает, что система пришла к равновесию. Оно свидетельствует о развитии нефизической неустойчивости и показывает, что расчет вышел за пределы корректного применения выбранной модели.

Результаты работы имеют значение для теории плазмы и для создания вычислительных методов. Исследование показывает, в каких условиях псевдопотенциальные квазиклассические модели еще дают надежное описание, а когда для расчетов требуются более строгие квантовые подходы.

Кроме того, ученые определили, каких физических эффектов не хватает моделям, основанным только на парном взаимодействии, при переходе к сильной неидеальности плазмы. Это позволяет точнее оценивать надежность подобных расчетов для разных систем.

«Квазиклассические модели по-прежнему востребованы, поскольку позволяют исследовать сложные плазменные системы без чрезмерных вычислительных затрат. Наш результат задает практическую границу применимости таких подходов: показывает, где они еще обеспечивают физически надежное описание, а где уже нет. При этом пока не существует основанного на первых принципах способа надежно расширить эту область применимости», — отметил Павел Левашов, доктор физико-математических наук и заведующий кафедрой физики высокотемпературных процессов МФТИ.

Работу опубликовали в журнале Physics of Plasmas. Исследование поддержал Российский научный фонд в рамках гранта №24-19-00746.