Path-integral Monte Carlo simulations for interacting few-electron quantum dots with spin-orbit coupling

We develop path-integral Monte Carlo simulations for a parabolic two-dimensional (2D) quantum dot containing N interacting electrons in the presence of Dresselhaus and/or Rashba spin-orbit couplings. Our method solves in a natural way the spin contamination problem and allows for numerically exact finite-temperature results at weak spin-orbit coupling. For N < 10 electrons, we present data for the addition energy, the particle density, and the total spin S in the Wigner molecule regime of strong Coulomb interactions. We identify magic numbers at N=3 and N=7 via a peak in the addition energy. These magic numbers differ both from weak interaction and classical predictions, and are stable with respect to (weak) spin-orbit couplings.