Configuration Path Integral Monte Carlo Approach to the Static Density Response of the Warm Dense Electron Gas

Precise knowledge of the static density response function (SDRF) of the uniform electron gas (UEG) serves as key input for numerous applications, most importantly for density functional theory beyond generalized gradient approximations. Here we extend the configuration path integral Monte Carlo (CPIMC) formalism that was previously applied to the spatially uniform electron gas to the case of an inhomogeneous electron gas by adding a spatially periodic external potential. This procedure has recently been successfully used in permutation blocking path integral Monte Carlo simulations (PB-PIMC) of the warm dense electron gas [Dornheim \textit{et al.}, Phys. Rev. E in press, arXiv:1706.00315], but this method is restricted to low and moderate densities. Implementing this procedure into CPIMC allows us to obtain exact finite temperature results for the SDRF of the electron gas at \textit{high to moderate densities} closing the gap left open by the PB-PIMC data. In this paper we demonstrate how the CPIMC formalism can be efficiently extended to the spatially inhomogeneous electron gas and present the first data points. Finally, we discuss finite size errors involved in the quantum Monte Carlo results for the SDRF in detail and present a solution how to remove them that is based on a generalization of ground state techniques

Ion-Streaming Induced Order Transition in 3D Dust Clusters

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

Quantum Breathing Mode of Interacting Particles in a One-dimensional Harmonic Trap

Extending our previous work, we explore the breathing mode---the [uniform] radial expansion and contraction of a spatially confined system. We study the breathing mode across the transition from the ideal quantum to the classical regime and confirm that it is not independent of the pair interaction strength (coupling parameter). We present the results of time-dependent Hartree-Fock simulations for 2 to 20 fermions with Coulomb interaction and show how the quantum breathing mode depends on the particle number. We validate the accuracy of our results, comparing them to exact Configuration Interaction results for up to 8 particles