.jpg)
Физики из Объединенного института высоких температур РАН и МФТИ научились проводить более устойчивые и физически корректные компьютерные эксперименты со слабо вырожденной водородной плазмой в области низких температур.
В недрах звезд, в установках управляемого термоядерного синтеза или в экспериментах с лазерным сжатием водород существует в форме плазмы. Чтобы предсказать поведение такой субстанции, необходимо знать ее фундаментальные свойства, включая уравнение состояния, степень ионизации, механизм образования атомов и молекул. Теоретически рассчитать эти параметры для неидеальной плазмы чрезвычайно сложно из-за необходимости учитывать квантовые эффекты и коллективные взаимодействия между огромным числом частиц.
Золотым стандартом считается квантовый метод Монте-Карло с интегралами по траекториям, который дает очень точные результаты. Однако он предполагает колоссальные вычислительные затраты при попытке описать системы, состоящие из большого количества частиц.
Для изучения плазмы в режиме, где квантовое вырождение электронов мало, ученые часто используют более быстрые квазиклассические методы. В этом случае частицы движутся по классическим траекториям, но их взаимодействие описывается не обычным кулоновским потенциалом, а учитывает квантовый принцип неопределенности и обменное взаимодействие.
Ученые поставили перед собой задачу модифицировать метод молекулярной динамики так, чтобы сохранить вычислительную эффективность. В новом подходе принят во внимание характерный размер электронов при их взаимодействии. В результате эффективное отталкивание между электронами с одинаковой проекцией начинает проявляться на больших расстояниях между частицами, чем в исходной модели. Такая модификация предотвращает образование в моделировании искусственных связанных состояний, нарушающих принцип Паули.
Ученые провели масштабные серии расчетов методом молекулярной динамики с модифицированным псевдопотенциалом для невырожденной водородной плазмы с малым параметром вырождения. Как и ожидалось, при температурах выше 50 тысяч градусов расхождения между предложенным и более точным квантовым методом Монте-Карло оказались минимальными, составив менее 1% для энергии и давления. Однако по мере снижения температуры расхождения нарастали, и в самой холодной рассмотренной точке энергия из модели молекулярной динамики оказалась занижена на десятки процентов.
Анализ выявил важное методическое ограничение используемой модели. Модификация силового взаимодействия позволила устранить нефизичное сближение электронов с одинаковой проекцией спина, однако дальнейшее сравнение с квантово статистическими расчетами показало, что для более точного воспроизведения свойств системы требуется уточнение формы эффективного парного потенциала и его градиента.
Предложенная модификация сил стабилизирует моделирование и позволяет исследовать структуру водородной плазмы в ранее труднодоступном низкотемпературном режиме.
Метод может использоваться для масштабных параметрических расчетов уравнения состояния водородной и дейтериевой плазмы, важных для задач астрофизики и инерциального термоядерного синтеза, тогда как дальнейшие исследования направлены на проверку гипотезы о расхождении градиента псевдопотенциала с точным взаимодействием и на обобщение предложенной модификации сил на более сложные многокомпонентные плазменные системы, включая симуляции с переменным числом частиц для более надежного определения термодинамического предела.
