Ionization calculations using classical molecular dynamics

By performing an ensemble of molecular dynamics simulations, the model-dependent ionization state is computed for strongly interacting systems self-consistently. This is accomplished through a free energy minimization framework based on the technique of thermodynamic integration. To illustrate the method, two simple models applicable to partially ionized 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-neutralizing background. Particular attention is given to the transition between atomic gas and ionized plasma, where the effect of neutral interactions is explored beyond commonly used models in the chemical picture. Furthermore, pressure ionization 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 ionized warm dense matter.