Search

Accurate knowledge of the properties of hydrogen at high compression is crucial for astrophysics (e.g. planetary and stellar interiors, brown dwarfs, atmosphere of compact stars) and laboratory experiments, including inertial confinement fusion. There exists experimental data for the equation of state, conductivity, and Thomson scattering spectra. However, the analysis of the measurements at extreme pressures and temperatures typically involves additional model assumptions, which makes it difficult to assess the accuracy of the experimental data. rigorously. On the other hand, theory and modeling have produced extensive collections of data. They originate from a very large variety of models and simulations including path integral Monte Carlo (PIMC) simulations, density functional theory (DFT), chemical models, machine-learned models, and combinations thereof. At the same time, each of these methods has fundamental limitations (fermion sign problem in PIMC, approximate exchange-correlation functionals of DFT, inconsistent interaction energy contributions in chemical models, etc.), so for some parameter ranges accurate predictions are difficult. Recently, a number of breakthroughs in first principle PIMC and DFT simulations were achieved which are discussed in this review. Here we use these results to benchmark different simulation methods. We present an update of the hydrogen phase diagram at high pressures, the expected phase transitions, and thermodynamic properties including the equation of state and momentum distribution. Furthermore, we discuss available dynamic results for warm dense hydrogen, including the conductivity, dynamic structure factor, plasmon dispersion, imaginary-time structure, and density response functions. We conclude by outlining strategies to combine different simulations to achieve accurate theoretical predictions.

Magnetized plasmas are well known to exhibit a rich spectrum of collective modes. Here, we focus on the density modes in dense or cold plasmas, where strong coupling effects alter the mode spectrum known from traditional weakly coupled plasmas. In particular, we study the dynamic structure factor (DSF) of the magnetized one-component plasma with molecular dynamics simulations. Extending our previous results [H.~K\"ahlert and M.~Bonitz, Phys. Rev. Research \textbf{2022}, 4, 013197], it is shown that Bernstein modes can be observed in the weakly magnetized regime, where they are found below the upper hybrid frequency, provided the coupling strength is sufficiently low. We investigate the DSF for a variety of different wave numbers and plasma parameters and show that even small magnetization can give rise to a strong zero-frequency mode perpendicular to the magnetic field and change the dispersion as well as the damping of the upper hybrid mode.

The study of matter under extreme densities and temperatures as they occur e.g. in astrophysical objects and nuclear fusion applications has emerged as one of the most active frontiers in physics, material science, and related disciplines. In this context, a key quantity is given by the dynamic structure factor $S(\mathbf{q},\omega)$, which is probed in scattering experiments -- the most widely used method of diagnostics at these extreme conditions. In addition to its crucial importance for the study of warm dense matter, the modelling of such dynamic properties of correlated quantum many-body systems constitutes one of the most fundamental theoretical challenges of our time. Here we report a hitherto unexplained \emph{roton feature} in $S(\mathbf{q},\omega)$ of the warm dense electron gas, and introduce a microscopic explanation in terms of a new \emph{electronic pair alignment} model. This new paradigm will be highly important for the understanding of warm dense matter, and has a direct impact on the interpretation of scattering experiments. Moreover, we expect our results to give unprecedented insights into the dynamics of a number of correlated quantum many-body systems such as ultracold helium, dipolar supersolids, and bilayer heterostructures.

Plasmas in strong magnetic fields have been mainly studied in two distinct limiting cases--that of weak and strong nonideality with very different physical properties. While the former is well described by the familiar theory of Braginskii, the latter regime is closer to the behavior of a Coulomb liquid. Here we study in detail the transition between both regimes. We focus on the evolution of the dynamic structure factor of the magnetized one-component plasma from weak to strong coupling, which is studied with first-principle molecular dynamics simulations. The simulations show the vanishing of Bernstein modes and the emergence of higher harmonics of the upper hybrid mode across the magnetic field, a redistribution of spectral power between the two main collective modes under oblique angles, and a suppression of plasmon damping along the magnetic field. Comparison with results from various models, including the random phase approximation, a Mermin-type dielectric function, and the Quasi-Localized Charge Approximation show that none of the theories is capable of reproducing the crossover that occurs when the coupling parameter is on the order of unity. The findings are relevant to the scattering spectra, stopping power, and transport coefficients of correlated magnetized plasmas.

The self-diffusion phenomenon in a two-dimensional dusty plasma at extremely strong (effective) magnetic fields is studied experimentally and by means of molecular dynamics simulations. In the experiment the high magnetic field is introduced by rotating the particle cloud and observing the particle trajectories in a co-rotating frame, which allows reaching effective magnetic fields up to 3000 Tesla. The experimental results confirm the predictions of the simulations: (i) super-diffusive behavior is found at intermediate time-scales and (ii) the dependence of the self-diffusion coefficient on the magnetic field is well reproduced., Comment: accepted by Physical Review E

This work examines the formation of wake fields caused by ions streaming around a charged dust particle, using three-dimensional particle-in-cell (PIC) simulations with charge-neutral collisions included. The influence of an external driving electric field, which leads to a non-Maxwellian distribution of ions, is investigated in detail. The wake features formed for non-Maxwellian ions exhibit significant deviations from those observed within the model of a shifted Maxwellian distribution. The dependence of the peak amplitude and position of the wake potential upon the degree of collisionality is analyzed for a wide range of streaming velocities (Mach numbers). In contrast to a shifted Maxwellian distribution of ions, the drift-driven non-Maxwellian distribution exhibits an increase of the wake amplitude of the first attractive peak with increase in collisionality for high streaming velocities. At very low Mach numbers, collision-induced amplification is observed for Maxwellian as well as non-Maxwellian distributions.

Comparative studies of capacitively coupled radio-frequency discharges in helium and argon at pressures between 10 and 80 Pa are presented applying two different fluid modelling approaches as well as two independently developed particle-in-cell/Monte Carlo collision (PIC/MCC) codes. The focus is on the analysis of the range of applicability of a recently proposed fluid model including an improved drift-diffusion approximation for the electron component as well as its comparison with fluid modelling results using the classical drift-diffusion approximation and benchmark results obtained by PIC/MCC simulations. Main features of this time- and space-dependent fluid model are given. It is found that the novel approach shows generally quite good agreement with the macroscopic properties derived by the kinetic simulations and is largely able to characterize qualitatively and quantitatively the discharge behaviour even at conditions when the classical fluid modelling approach fails. Furthermore, the excellent agreement between the two PIC/MCC simulation codes using the velocity Verlet method for the integration of the equations of motion verifies their accuracy and applicability.

The electrostatic potential of a moving dust grain in a complex plasma with magnetized ions is computed using linear response theory, thereby extending our previous work for unmagnetized plasmas [P. Ludwig et al., New J. Phys. 14, 053016 (2012)]. In addition to the magnetic field, our approach accounts for a finite ion temperature as well as ion-neutral collisions. Our recently introduced code \texttt{Kielstream} is used for an efficient calculation of the dust potential. Increasing the magnetization of the ions, we find that the shape of the potential crucially depends on the Mach number $M$. In the regime of subsonic ion flow ($M<1$), a strong magnetization gives rise to a potential distribution that is qualitatively different from the unmagnetized limit, while for $M>1$ the magnetic field effectively suppresses the plasma wakefield.

Molecular Dynamics (MD) simulations of a strongly coupled binary ionic mixture have revealed the appearance of sharp minima in the species resolved dynamical density fluctuation spectra. This phenomenon is reminiscent of the well-known Fano anti-resonance, occurring in various physical processes. We give a theoretical analysis using the Quasi Localized Charge Approximation, and demonstrate that the observed phenomenon in the equilibrium spectrum is a novel manifestation of the Fano mechanism, that occurs at characteristic frequencies of the system different from the conventional classical Fano frequencies.

A rotating dusty plasma apparatus was constructed to provide the possibility of experimental emulation of extremely high magnetic fields by means of the Coriolis force, observable in a co-rotating measurement frame. We present collective excitation spectra for different rotation rates with a magnetic induction equivalent of up to 3200 Tesla. We identify the onset of magnetoplasmon-equivalent mode dispersion in the rotating macroscopic two-dimensional single-layer dusty plasma. The experimental results are supported by molecular dynamics simulations of 2D magnetized Yukawa systems., Comment: 5 pages, 6 figures, accepted by PRL

Dust Dynamics Simulations utilizing a dynamical screening approach are performed to study the effect of ion-streaming on the self-organized structures in a three-dimensional spherically confined complex (dusty) plasma. Varying the Mach number M - the ratio of ion drift velocity to the sound velocity, the simulations reproduce the experimentally observed cluster configurations in the two limiting cases: at M=0 strongly correlated crystalline structures consisting of nested spherical shells (Yukawa balls) and, for M\geq1, flow-aligned dust chains, respectively. In addition, our simulations reveal a discontinuous transition between these two limits. It is found that already a moderate ion drift velocity (M\approx0.1) destabilizes the highly ordered Yukawa balls and initiates an abrupt melting transition. The critical value of M is found to be independent of the cluster size.

The theoretical description of complex (dusty) plasmas requires multiscale concepts that adequately incorporate the correlated interplay of streaming electrons and ions, neutrals, and dust grains. Knowing the effective dust-dust interaction, the multiscale problem can be effectively reduced to a one-component plasma model of the dust subsystem. The goal of the present publication is a systematic evaluation of the electrostatic potential distribution around a dust grain in the presence of a streaming plasma environment by means of two complementary approaches: (i) a high precision computation of the dynamically screened Coulomb potential from the dynamic dielectric function, and (ii) full 3D particle-in-cell simulations, which self-consistently include dynamical grain charging and non-linear effects. The applicability of these two approaches is addressed.

A theoretical description of the radial density profile for charged particles with Yukawa interaction in a harmonic trap is described. At strong Coulomb coupling shell structure is observed in both computer simulations and experiments. Correlations responsible for such shell structure are described here using a recently developed model based in density functional theory. A wide range of particle number, Coulomb coupling, and screening lengths is considered within the fluid phase. A hypernetted chain approximation shows the formation of shell structure, but fails to give quantitative agreement with Monte Carlo simulation results at strong coupling. Significantly better agreement is obtained within the hypernetted chain structure using a renormalized coupling constant, representing bridge function corrections., Comment: 9 pages, 10 figures

A system of N classical Coulomb charges trapped in a harmonic potential displays shell structure and orientational ordering. The local density profile is well understood from theory, simulation, and experiment. Here, pair correlations are considered for this highly inhomogeneous system for both the fluid and ordered states. In the former, it is noted that there is a close relationship to pair correlations in the uniform OCP. For the ordered state, it is shown that the disordered "tiling" is closely related to the ground state Thomson sites for a single sphere., Comment: 8 pages, 3 figues. Results presented at "Strongly Coupled Coulomb Systems", Budapest, Hungary, July 2011. Conference proceedings to be published in "Contribution to Plasma Physics"