Search

Warm dense matter (WDM) -- an exotic state of highly compressed matter -- has attracted high interest in recent years in astrophysics and for dense laboratory systems. At the same time, this state is extremely difficult to treat theoretically. This is due to the simultaneous appearance of quantum degeneracy, Coulomb correlations and thermal effects, as well as the overlap of plasma and condensed phases. Recent breakthroughs are due to the successful application of density functional theory (DFT) methods which, however, often lack the necessary accuracy and predictive capability for WDM applications. The situation has changed with the availability of the first \textit{ab initio} data for the exchange-correlation free energy of the warm dense uniform electron gas (UEG) that were obtained by quantum Monte Carlo (QMC) simulations, for recent reviews, see Dornheim \textit{et al.}, Phys. Plasmas \textbf{24}, 056303 (2017) and Phys. Rep. \textbf{744}, 1-86 (2018). In the present article we review recent further progress in QMC simulations of the warm dense UEG: namely, \textit{ab initio} results for the static local field correction $G(q)$ and for the dynamic structure factor $S(q,\omega)$. These data are of key relevance for the comparison with x-ray scattering experiments at free electron laser facilities and for the improvement of theoretical models. In the second part of this paper we discuss simulations of WDM out of equilibrium. The theoretical approaches include Born-Oppenheimer molecular dynamics, quantum kinetic theory, time-dependent DFT and hydrodynamics. Here we analyze strengths and limitations of these methods and argue that progress in WDM simulations will require a suitable combination of all methods. A particular role might be played by quantum hydrodynamics, and we concentrate on problems, recent progress, and possible improvements of this method.

The dynamical structure factor (DSF) of strongly coupled ions in dense plasmas with partially and strongly degenerate electrons is investigated. The main focus is on the impact of electronic correlations (non-ideality) on the ionic DSF. The latter is computed by carrying out molecular dynamics (MD) simulations with a screened ion-ion interaction potential. The electronic screening is taken into account by invoking the Singwi-Tosi-Land-Sj\"olander approximation, and compared to the MD simulation data obtained considering the electronic screening in the random phase approximation and using the Yukawa potential. %This allows us to gain insight into the impact of the electronic non-ideality on the ionic DSF. We find that electronic correlations lead to lower values of the ion-acoustic mode frequencies and to an extension of the applicability limit with respect to the wave-number of a hydrodynamic description. Moreover, we show that even in the limit of weak electronic coupling, electronic correlations have a non-negligible impact on the ionic longitudinal sound speed. Additionally, the applicability of the Yukawa potential with an adjustable screening parameter is discussed, which will be of interest, e.g., for the interpretation of experimental results for the ionic DSF of dense plasmas.

Solids facing a plasma are a common situation in many astrophysical systems and laboratory setups. Moreover, many plasma technology applications rely on the control of the plasma-surface interaction. However, presently often a fundamental understanding of them is missing, so most technological applications are being developed via trial and error. In the majority of plasma simulations surface processes are either neglected or treated via phenomenological parameters such as sticking coefficients, sputter rates or secondary electron emission coefficients. However, those parameters are known only in some cases and with very limited accuracy. Similarly, while surface physics simulations have often studied the impact of single ions or neutrals, so far, the influence of a plasma medium and correlations between successive impacts have not been taken into account. Such an approach cannot have predictive power. In this paper we discuss in some detail the physical processes a the plasma-solid interface which brings us to the necessity of coupled plasma-solid simulations. We briefly summarize relevant theoretical methods from solid state and surface physics that are suitable to contribute to such an approach and identify four methods. The first are mesoscopic simulations such as kinetic Monte Carlo (KMC) and molecular dynamics (MD) that are able to treat complex processes on large scales but neglect electronic effects. The second are quantum kinetic methods based on the quantum Boltzmann equation that give access to a more accurate treatment of surface processes using simplifying models for the solid. The third approach are ab initio simulations of surface process that are based on density functional theory (DFT) and time-dependent DFT. The fourths are nonequilibrium Green functions that able to treat correlation effects in the material and at the interface.

The structural properties of strongly coupled ions in dense plasmas with moderately to strongly degenerate electrons are investigated in the framework of the one-component plasma model of ions interacting through a screened pair interaction potential. Special focus is put on the description of the electronic screening in the Singwi-Tosi-Land-Sjoelander (STLS) approximation. Different crosschecks and analyses using ion potentials obtained from ground-state quantum Monte Carlo data, the random phase approximation (RPA), and existing analytical models are presented for the computation of the structural properties, such as the pair distribution and the static structure factor, of strongly coupled ions. The results are highly sensitive to the features of the screened pair interaction potential. This effect is particularly visible in the static structure factor. The applicability range of the screened potential computed from STLS is identified in terms of density and temperature of the electrons. It is demonstrated that at r_s > 1, where rs is the ratio of the mean inter-electronic distance to the Bohr radius, electronic correlations beyond RPA have a non-negligible effect on the structural properties. Additionally, the applicability of the hypernetted chain approximation for the calculation of the structural properties using the screened pair interaction potential is analyzed employing the effective coupling parameter approach.

Strongly coupled plasmas in which the interaction energy exceeds the kinetic energy play an important role in many astrophysical and laboratory systems including compact stars, laser plasmas and dusty plasmas. They exhibit many unusual collective properties, such as liquid or crystalline behaviour, peculiar oscillation spectra and transport properties. Recently, strongly coupled plasmas were studied in the presence of a strong magnetic field by computer simulations, and strong modifications of their transport properties and oscillation spectra were observed. While strong magnetization is common in stellar systems it is practically impossible to achieve in complex plasmas due to the large mass of the dust particles. Here we discuss a recently demonstrated approach to achieve very strong "magnetization" by a rotation of the neutral gas, and we present new results for macroscopic two-dimensional systems.

We propose and demonstrate a concept that mimics the magnetization of the heavy dust particles in a complex plasma while leaving the properties of the light species practically unaffected. It makes use of the frictional coupling between a complex plasma and the neutral gas, which allows to transfer angular momentum from a rotating gas column to a well-controlled rotation of the dust cloud. This induces a Coriolis force that acts exactly as the Lorentz force in a magnetic field. Experimental normal mode measurements for a small dust cluster with four particles show excellent agreement with theoretical predictions for a magnetized plasma.

A first-principle study of the collective oscillation spectrum of a strongly correlated one-component plasma in a strong magnetic field is presented. The spectrum consists of six fundamental modes which are found to be in good agreement with results from the Quasi-Localized Charge approximation. At high frequencies additional modes are observed that include Bernstein-type oscillations and their transverse counterparts which are of importance for the high-frequency optical and transport properties of these plasmas.

The excitation spectra of two-dimensional strongly coupled (liquid-like) Yukawa systems under the influence of a magnetic field perpendicular to the plane are found to sustain additional high-frequency modes at multiples of the magnetoplasmon. These modes are reminiscent of the well-known Bernstein modes but show a number of important differences due to the strong coupling of the particles.

The results of a recent fluid theory for the multipole modes of a Yukawa plasma in a spherical confinement [H. K\"{a}hlert and M. Bonitz, Phys. Rev. E \textbf{82}, 036407 (2010)] are compared with molecular dynamics simulations and the exact $N$-particle eigenmodes in the crystalline phase. Simulations confirm the existence of high order modes found in cold fluid theory. We investigate the influence of screening, coupling and friction on the mode spectra in detail. Good agreement between theory and simulation is found for weak to moderate screening and low order modes. The relations between the breathing mode in the fluid theory, simulation and the crystal eigenmode are investigated in further detail.

The normal modes of a three-dimensional Yukawa plasma in an isotropic, harmonic confinement are investigated by solving the linearized cold fluid equations. The eigenmodes are found analytically and expressed in terms of hypergeometric functions. It is found that the mode frequencies solely depend on the dimensionless plasma parameter $\xi=\kappa R$, where $R$ is the plasma radius and $\kappa$ the inverse screening length. The eigenfrequencies increase monotonically with $\xi$ and saturate in the limit $\xi\to\infty$. Compared with the results in the Coulomb limit~[D. H. E. Dubin, Phys. Rev. Lett. \textbf{66}, 2076 (1991)], we find a new class of modes characterized by the number $n$ which determines the number of radial nodes in the perturbed potential. These modes originate from the degenerate bulk modes of the Coulomb system. Analytical formulas for the eigenfrequencies are derived for limiting cases.

The existence of plasma oscillations at multiples of the magnetoplasmon frequency in a strongly coupled two-dimensional magnetized Yukawa plasma is reported, based on extensive molecular dynamics simulations. These modes are the analogues of Bernstein modes which are renormalized by strong interparticle correlations. Their properties are theoretically explained by a dielectric function incorporating the combined effect of a magnetic field, strong correlations and finite temperature.

A theoretical description of shell structure for charged particles in a harmonic trap is explored at strong coupling conditions of $\Gamma$ = 50 and 100. The theory is based on an extension of the hypernetted chain approximation to confined systems plus a phenomenological representation of associated bridge functions. Predictions are compared to corresponding Monte Carlo simulations and quantitative agreement for the radial density profile is obtained., Comment: 9 pages, 5 figures. Presented at the 13th International Conference on the Physics of Non-Ideal Plasmas (PNP 13) held in Chernogolovka, Russia (September 13-18, 2009). Proceedings to be published in "Contributions to Plasma Physics" (Dec. 2009-Jan. 2010)

The correlation buildup and the formation dynamics of the shell structure in a spherically confined one-component plasma are studied. Using Langevin dynamics simulations the relaxation processes and characteristic time scales and their dependence on the pair interaction and dissipation in the plasma are investigated. While in systems with Coulomb interaction (e.g. trapped ions) in a harmonic confinement shell formation starts at the plasma edge and proceeds inward, this trend is significantly weakened for dusty plasmas with Yukawa interaction. With a suitable change of the confinement conditions the crystallization scenario can be externally controlled.

A theoretical description for the radial density profile of a finite number of identical charged particles confined in a harmonic trap is developed for application over a wide range of Coulomb coupling (or, equivalently, temperatures) and particle numbers. A simple mean field approximation neglecting correlations yields a density profile which is monotonically decreasing with radius for all temperatures, in contrast to molecular dynamics simulations and experiments showing shell structure at lower temperatures. A more complete theoretical description including charge correlations is developed here by an extension of the hypernetted chain approximation, developed for bulk fluids, to the confined charges. The results reproduce all of the qualitative features observed in molecular dynamics simulations and experiments. These predictions are then tested quantitatively by comparison with new benchmark Monte Carlo simulations. Quantitative accuracy of the theory is obtained for the selected conditions by correcting the hypernetted chain approximation with a representation for the associated bridge functions., Comment: 10 figures, submitted to Physical Review E

A combined theoretical and experimental analysis of the normal modes of three-dimensional spherially confined Yukawa clusters is presented. Particular attention is paid to the breathing mode and the existence of multiple monopole oscillations in Yukawa systems. Finally, the influence of dissipation on the mode spectrum is investigated.

Small three-dimensional strongly coupled clusters of charged particles in a spherical confinement potential arrange themselves in nested concentric shells. If the particles are immersed into a background plasma the interaction is screened. The cluster shell configuration is known to be sensitive to the screening strength. With increased screening an increased population of the inner shell(s) is observed. Here, we present a detailed analysis of the ground state shell configurations and configuration changes in a wide range of screening parameters for clusters with particle numbers $N$ in the range of 11 to 60. We report three types of anomalous behaviors which are observed upon increase of screening, at fixed $N$ or for an increase of $N$ at fixed screening. The results are obtained by means of extensive first principle molecular dynamics simulations.

Recently the occurrence probabilities of ground- and metastable states of three-dimensional Yukawa clusters with 27 and 31 particles have been analyzed in dusty plasma experiments [Block et al., Physics of Plasmas 15, 040701 (2008)]. There it was found that, in many cases, the ground state appeared substantially less frequently than excited states. Here we analyze this question theoretically by means of molecular dynamics (MD) and Monte Carlo simulations and an analytical method based on the canonical partition function. We confirm that metastable states can occur with a significantly higher probability than the ground state. The results strongly depend on the screening parameter of the Yukawa interaction and the damping coefficient used in the MD simulations. The analytical method allows one to gain insight into the mechanisms being responsible for the occurrence probabilities of metastable states in strongly correlated finite systems.