Abstract: Wigner molecules formed at high magnetic fields in circular and elliptic quantum dots are studied by exact diagonalization (ED) and unrestricted Hartree-Fock (UHF) methods with multicenter basis of displaced lowest Landau level wave functions. The broken symmetry states with semi-classical charge density constructed from superpositions of the ED solutions are compared to the UHF results. UHF overlooks the dependence of the few-electron wave functions on the actual relative positions of electrons localized in different charge puddles and partially compensates for this neglect by an exaggerated separation of charge islands which are more strongly localized than in the exact broken-symmetry states. (C) 2004 Elsevier B.V. All rights reserved.