Laboratory astroparticle physics

Abstract:

We present a novel scheme for modelling quantum plasmas in the warm dense matter (WDM) regime via a hybrid smoothed particle hydrodynamic - molecular dynamic treatment, here referred to as ‘Bohm SPH’. This treatment is founded upon Bohm’s interpretation of quantum mechanics for partially degenerate fluids, does not apply the Born-Oppenheimer approximation, and is computationally tractable, capable of modelling dynamics over ionic timescales at electronic time resolution. Bohm SPH is also capable of modelling non-Gaussian electron wavefunctions. We present an overview of our methodology, validation tests of the single particle case including the hydrogen 1s wavefunction, and comparisons to simulations of a warm dense hydrogen system performed with wave packet molecular dynamics.

Abstract:

Modeling many-body quantum systems with strong interactions is one of the core challenges of modern physics. A range of methods has been developed to approach this task, each with its own idiosyncrasies, approximations, and realm of applicability. Perhaps the most successful and ubiquitous of these approaches is density functional theory (DFT). Its Kohn-Sham formulation has been the basis for many fundamental physical insights, and it has been successfully applied to fields as diverse as quantum chemistry, condensed matter and dense plasmas. Despite the progress made by DFT and related schemes, however, there remain many problems that are intractable for existing methods. In particular, many approaches face a huge computational barrier when modeling large numbers of coupled electrons and ions at finite temperature. Here, we address this shortfall with a new approach to modeling many-body quantum systems. Based on the Bohmian trajectories formalism, our new method treats the full particle dynamics with a considerable increase in computational speed. As a result, we are able to perform large-scale simulations of coupled electron-ion systems without employing the adiabatic Born-Oppenheimer approximation.

Abstract:

Our understanding of physical systems generally depends on our ability to match complex computational modelling with measured experimental outcomes. However, simulations with large parameter spaces suffer from inverse problem instabilities, where similar simulated outputs can map back to very different sets of input parameters. While of fundamental importance, such instabilities are seldom resolved due to the intractably large number of simulations required to comprehensively explore parameter space. Here we show how Bayesian machine learning can be used to address inverse problem instabilities, and apply it to two popular experimental diagnostics in plasma physics. We find that the extraction of information from measurements simply on the basis of agreement with simulations is unreliable, and leads to a significant underestimation of uncertainties. We describe how to statistically quantify the effect of unstable inverse models, and describe an approach to experimental design that mitigates its impact.

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.