Classical molecular dynamics simulations of hydrogen plasmas and development of an analytical statistical model for computational validity assessment

M. A. Gigosos,D. Gonzalez-Herrero,R. Florido,A. Calisti,S. Ferri,B. Talin

DOI: https://doi.org/10.1103/PhysRevE.98.033307

2024-02-05

Abstract:Classical molecular dynamics simulations of hydrogen plasmas have been performed with emphasis on the analysis of equilibration process. Theoretical basis of simulation model as well as numerically relevant aspects -- such as the proper choice and definition of simulation units -- are discussed in detail, thus proving a thorough implementation of the computer simulation technique. Because of lack of experimental data, molecular dynamics simulations are often considered as "idealized computational experiments" for benchmarking of theoretical models. However, these simulations are certainly challenging and consequently a validation procedure is also demanded. In this work we develop an analytical statistical equilibrium model for computational validity assessment of plasma particle dynamics simulations. Remarkable agreement between model and molecular dynamics results including a classical treatment of ionization-recombination mechanism is obtained for a wide range of plasma coupling parameter. Furthermore, the analytical model provides guidance to securely terminate simulation runs once the equilibrium stage has been reached, which in turn gives confidence on the statistics that potentially may be extracted from time-histories of simulated physical quantities.

Plasma Physics

What problem does this paper attempt to address?

The paper aims to address the following key issues:

1. **Analysis of Equilibrium Processes in Plasma Dynamics Simulations**: The paper primarily focuses on the classical molecular dynamics simulation of hydrogen plasma and discusses in detail the theoretical basis of the simulation model and its numerical aspects, such as the appropriate choice and definition of units. Through this work, the comprehensive implementation of computer simulation technology is demonstrated.

2. **Development of Methods for Verifying Computational Effectiveness**: Due to the lack of experimental data, molecular dynamics simulations are often regarded as idealized computational experiments for the benchmarking of theoretical models. However, these simulations are challenging, thus requiring a verification procedure. This paper develops an analytical statistical equilibrium model for evaluating the computational effectiveness of plasma particle dynamics simulations. The model shows significant consistency with molecular dynamics results over a wide range of plasma coupling parameters and can guide when to safely terminate the simulation run after reaching the equilibrium stage.

3. **Handling of Ionization-Recombination Mechanisms**: The model proposed in the paper explicitly includes the ionization-recombination process within the classical framework, making the recombined ion and neutral pairs part of the final ionization equilibrium and equilibrium state. This approach is suitable for studying strongly coupled plasmas and is not limited to fully ionized scenarios.

In summary, the focus of the paper is on developing a reliable method for plasma particle dynamics simulation and validating its effectiveness through an analytical model, providing new solutions particularly in handling ionization-recombination mechanisms.