Ionisation Calculations using Classical Molecular Dynamics
arXiv:2409.01078 · doi:10.1103/PhysRevE.111.015204
Abstract
By performing an ensemble of molecular dynamics simulations, the model-dependent ionisation state is computed for strongly interacting systems self-consistently. This is accomplished through a free energy minimisation framework based on the technique of thermodynamic integration. To illustrate the method, two simple models applicable to partially ionised hydrogen plasma are presented in which pair potentials are employed between ions and neutral particles. Within the models, electrons are either bound in the hydrogen ground state or distributed in a uniform charge-neutralising background. Particular attention is given to the transition between atomic gas and ionised plasma, where the effect of neutral interactions is explored beyond commonly used models in the chemical picture. Furthermore, pressure ionisation is observed when short range repulsion effects are included between neutrals. The developed technique is general, and we discuss the applicability to a variety of molecular dynamics models for partially ionised warm dense matter.
14 pages, 8 figures
References in corpus (21)
- The Uniform Electron Gas at Warm Dense Matter Conditions
- Ab initio simulation of warm dense matter
- Temperature-dependent quantum pair potentials and their application to dense partially ionized hydrogen plasmas
- First principles simulations of dense hydrogen
- Carbon ionization at Gbar pressures: an ab initio perspective on astrophysical high-density plasmas
- Structural characteristics of strongly coupled ions in a dense quantum plasma
- Free-energy model for fluid helium at high density
- Dynamical structure factor of strongly coupled ions in a dense quantum plasma
- The equation of state of partially ionized hydrogen and deuterium plasma revisited
- Screening of a test charge in a free-electron gas at warm dense matter and dense non-ideal plasma conditions
- Direct Evaluation of the Phase Diagrams of Dense Multicomponent Plasmas by Integration of the Clapeyron Equations
- Ionization and transport in partially ionized multicomponent plasmas: Application to atmospheres of hot Jupiters
- An accurate equation of state for the one component plasma in the low coupling regime
- Improved calculations of mean ionization states with an average-atom model
- Classical molecular dynamics simulations of hydrogen plasmas and development of an analytical statistical model for computational validity assessment
- One--Component Plasma of a Million Particles via angular--averaged Ewald potential: A Monte Carlo study
- Investigating Mechanisms of State Localization in Highly-Ionized Dense Plasmas
- Thermodynamics of ionization and dissociation in hydrogen plasmas including fluctuations and magnetic fields
- Ion Modes in Dense Ionized Plasmas through Non-Adiabatic Molecular Dynamics
- Development of a new quantum trajectory molecular dynamics framework
- Modelling of warm dense hydrogen via explicit real time electron dynamics: Dynamic structure factors