Monte Carlo simulations of dense quantum plasmas

Fermionic path-integral Monte Carlo results for the uniform electron gas at finite temperature

The uniform electron gas (UEG) at finite temperature has recently attracted substantial interest due to the experimental progress in the field of warm dense matter. To explain the experimental data, accurate theoretical models for high-density plasmas are needed that depend crucially on the quality of the thermodynamic properties of the quantum degenerate nonideal electrons and of the treatment of their interaction with the positive background. Recent fixed-node path-integral Monte Carlo (RPIMC) data are believed to be the most accurate for the UEG at finite temperature, but they become questionable at high degeneracy when the Brueckner parameter rs=a/aB--the ratio of the mean interparticle distance to the Bohr radius--approaches 1. The validity range of these simulations and their predictive capabilities for the UEG are presently unknown. This is due to the unknown quality of the used fixed nodes and of the finite-size scaling from N=33 simulated particles (per spin projection) to the macroscopic limit. To analyze these questions, we present alternative direct fermionic path integral Monte Carlo (DPIMC) simulations that are independent from RPIMC. Our simulations take into account quantum effects not only in the electron system but also in their interaction with the uniform positive background. Also, we use substantially larger particle numbers (up to three times more) and perform an extrapolation to the macroscopic limit. We observe very good agreement with RPIMC, for the polarized electron gas, up to moderate densities around rs=4, and larger deviations for the unpolarized case, for low temperatures. For higher densities (high electron degeneracy), rs≲1.5, both RPIMC and DPIMC are problematic due to the increased fermion sign problem.

Correlation effects in partially ionized mass asymmetric electron-hole plasmas

The effects of strong Coulomb correlations in dense three-dimensional electron-hole plasmas are studied by means of unbiased direct path integral Monte Carlo simulations. The formation and dissociation of bound states, such as excitons and biexcitons, is analyzed and the density-temperature region of their appearance is identified. At high density, the Mott transition to the fully ionized metallic state (electron-hole liquid) is detected. Particular attention is paid to the influence of the hole to electron mass ratio M on the properties of the plasma. Above a critical value of about M=80 formation of a hole Coulomb crystal was recently verified [Bonitz, Phys. Rev. Lett. 95, 235006 (2005)] which is supported by additional results. Results are related to the excitonic phase diagram of intermediate valent Tm[Se,Te], where large values of M have been observed experimentally.