Pontus Svensson

Abstract:

The capability for molecular dynamics simulations to treat relativistic dynamics is extended by the inclusion of relativistic kinetic energy. In particular, relativistic corrections to the diffusion coefficient are considered for an argon gas modeled with a Lennard-Jones interaction. Forces are transmitted instantaneously without being retarded, an approximation that is allowed due to the short-range nature of the Lennard-Jones interaction. At a mass density of 1.4g/cm3, significant deviations from classical results are observed at temperatures above kBT≈0.05mc2, corresponding to an average thermal velocity of 32% of the speed of light. For temperatures approaching kBT≈mc2, the semirelativistic simulations agree with analytical results for hard spheres, which is seen to be a good approximation as far as diffusion effects are concerned.

Abstract:

Data generated for the figures in 'Development of a new quantum trajectory molecular dynamics framework' at https://dx.doi.org/10.1098/rsta.2022.0325 (and at https://doi.org/10.48550/arXiv.2211.08560) and statically compiled version of the code.

Abstract:

High energy density science is central for astrophysical and human-made fusion applications but is characterised by non-ideal plasma behaviour due to strong particle interactions, quantum effects, and relativistic corrections. In this thesis, two molecular dynamics (MD) formulations are presented along with their implementation, which address quantum and relativistic effects, respectively. First, an extension to wave packet molecular dynamics using anisotropic Gaussian states is presented, which is designed to model electron dynamics over ionic time scales in warm dense matter. Long-range interactions are treated with a generalised Ewald summation, and exchange effects are treated within a pairwise approximation. The MD formulation has been used to investigate electron dynamic structure factors (DSFs) and x-ray Thomson scattering, where electron and ion time scale features are extracted from a single computation. A semi-classical form for the DSF, that corrects for known quantum constraints, is provided. This method has been tested against explicit computations of the density response function in MD. The DSF is further discussed within a two-fluid model, parameterised by the equation of state and transport properties. By comparison with MD results - facilitated by Bayesian inference - the electron transport properties for a test system of warm dense hydrogen are extracted.

Second, relativistic corrections are investigated both due to kinematics and interactions. The velocity-dependent inertia of relativistic particles is seen to reduce diffusive transport for one-component plasmas, in line with analytical results. However, long-range electromagnetic interactions are modified due to the finite speed of light. This is accounted for in the MD model by time-evolving the long-range fields while the highly fluctuating short-range fields are approximated in a field-less description using either the electrostatic or Darwin approximation.