Structures of quantum 2D electron-hole plasmas

We investigate structures of 2D quantum electron-hole (e-h) plasmas by the direct path integral Monte Carlo method (PIMC) in a wide range of temperature, density and hole-to-electron mass ratio. Our simulation includes a region of appearance and decay of the bound states (excitons and biexcitons), the Mott transition from the neutral e-h plasma to metallic-like clusters, formation from clusters the hexatic-like liquid and formation of the crystal-like lattice.

💡 Research Summary

This paper presents a comprehensive first‑principles investigation of two‑dimensional (2D) quantum electron‑hole (e‑h) plasmas using the direct Path‑Integral Monte‑Carlo (PIMC) method. By varying three key parameters—temperature (T), the Brueckner parameter (r_s, which measures the ratio of the mean inter‑particle distance to the 3D exciton Bohr radius), and the hole‑to‑electron mass ratio (M)—the authors map out the structural evolution of the plasma from dilute, bound‑state dominated regimes to dense, metallic‑like phases.

The theoretical framework starts from the canonical partition function for a neutral two‑component plasma with equal numbers of electrons and holes (N_e = N_h). The density matrix ρ(q, r, σ; β) is expressed as a path integral that fully incorporates Coulomb interactions, quantum exchange (via permutation operators), and spin statistics for both species. This formulation is valid in 1D, 2D, and 3D, but the simulations focus on the 2D case relevant to semiconductor quantum wells and atomically thin materials.

In the simulations, the authors consider a wide range of r_s values (from 6 down to 0.25) and a fixed low reduced temperature T/Ry ≈ 0.007 (Ry being the 3D e‑h binding energy). The mass ratio M is varied up to 800, effectively rendering holes as heavy, quasi‑classical particles while electrons remain highly quantum‑delocalized. The PIMC approach represents each particle by a “polymer” of beads; the spatial spread of these beads reflects the particle’s quantum delocalization.

Key findings can be grouped into four regimes:

1. Bound‑state formation at low density and low temperature – For r_s ≈ 6, the bead clouds of electrons are several times larger than those of holes, indicating that electrons are loosely bound around relatively localized holes. Pair‑distribution functions g_ee, g_hh, and g_eh display pronounced peaks at distances of roughly half to one Bohr radius, confirming the presence of excitons, biexcitons, and larger multi‑particle clusters.

2. Temperature‑induced ionization (Mott transition) – Raising the temperature at fixed density reduces the height of the bound‑state peaks and increases the fraction of free carriers. This reflects a classic Mott transition: thermal de‑localization of electron beads makes the electron‑hole overlap insufficient for binding, leading to a plasma of essentially free carriers.

3. Density‑driven structural transitions of the heavy component – As the density increases (r_s → 2, 0.5, 0.25), the electron beads become so delocalized that they span the entire simulation cell, while the heavy holes begin to order. At r_s ≈ 0.5, holes arrange in a six‑fold (hexatic) liquid with local triangular coordination but without long‑range translational order. At the highest density examined (r_s ≈ 0.25), holes crystallize into a Wigner‑type triangular lattice, while the electron subsystem remains a degenerate, highly delocalized gas. The transition from a liquid‑like hexatic phase to a true crystal is driven by the competition between kinetic energy (still significant for electrons) and the dominant Coulomb repulsion among the massive holes.

4. Effect of extreme mass anisotropy – For very large M (e.g., 800), the hole de‑Broglie wavelength is much smaller than the inter‑hole spacing even at high densities, allowing holes to behave essentially as classical point charges. Consequently, the hole lattice observed at r_s ≈ 0.25 closely resembles the classic 2D Wigner crystal. However, finite‑size effects and periodic boundary conditions introduce filament‑like bead structures at the cluster edges, reminiscent of defect lines or grain boundaries in real crystals.

Overall, the study demonstrates that PIMC can capture subtle quantum‑statistical effects (exchange, spin, and delocalization) while providing clear visualizations of real‑space structures. The authors identify precise density–temperature windows where bound complexes exist, delineate the Mott density for ionization, and show that, in two‑component plasmas with large mass asymmetry, the heavy component can undergo a sequence of phase transitions: from bound‑state clusters to a hexatic liquid and finally to a Wigner crystal embedded in a degenerate electron sea.

These results have direct relevance for experimental systems such as coupled quantum wells, transition‑metal dichalcogenide heterostructures, and other 2D semiconductor platforms where electron‑hole interactions are strong and mass ratios can be engineered (e.g., via strain or band‑structure design). Understanding the interplay of binding, ionization, and crystallization is essential for designing devices that exploit excitonic condensation, charge‑density waves, or exotic quantum phases in low‑dimensional materials.