Физики из МФТИ выяснили, где квазиклассические модели плазмы теряют устойчивость

Исследователи из МФТИ и Объединенного института высоких температур РАН определили условия, при которых квазиклассические модели водородной плазмы начинают давать нефизичные результаты. Ученые показали, что при достаточно сильном взаимодействии частиц такие модели перестают корректно воспроизводить квантовые эффекты, связанные с принципом Паули. В расчетах начинают возникать несуществующие электрон-ионные кластеры, а сама моделируемая система становится неустойчивой. Работа опубликована в журнале Physics of Plasmas. Исследование поддержано Российским научным фондом (грант № 24-19-00746).

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

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

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

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

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

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

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

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

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

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

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

Источник: Физики из МФТИ выяснили, когда квазиклассические модели плазмы перестают работать — За науку