Path integral Monte Carlo simulation results are presented which allow to c
haracterize the equilibrium properties of a few-electron system in a 2D har
monic trap in rigorous manner. In particular, Wigner crystallization is obs
erved and investigated in a broad range of temperature and confinement stre
ngth. An analytical expression for the phase boundary of the Wigner crystal
is obtained.