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 function 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. We assume a spin-polarization of electrons at high magnetic field (0, 0, B) oriented perpendicular to the quantum dot plane and use the Landau gauge. In the ED, described in detail in Ref.[4], the single electron wave functions used for construction of the Slater determinants are obtained via diagonalization of the single electron Hamiltonian in the multicenter basis [7,8,9] of M displaced lowest Landau level wave functionswhere α is treated as variational parameter. In the present UHF approach the one-electron orbitals (1) are optimized self-consistently. We study up to N = 4 electrons, use the material data of GaAs and a basis of 12 centers (x k , y k ) put on an ellipse with a size determined variationally. Basis (1) Classical system of three electrons in an elliptical confinement potential with hω x = 3 meV andhω y = 4 meV possesses two equivalent lowest-energy configurations [cf. inset of Fig. 1] and the quantum system undergoes parity transformations [4] with the magnetic field [cf. Fig. 1]. Superposition [4] of the two lowest-energy eigenstates2 yields a broken-symmetry (BS) charge density with a distinct electron separation. Fig. 1 shows that in contrast to the exact ground-state energy the UHF energy estimate is a smooth function of the magnetic field.The charge densities of considered states are shown in Figs. 2(a) and 2(b) for two magnetic field values corresponding to the even-odd energy crossing presented in Fig. 1. The phase φ in the BS state is chosen such that the electrons are localized at the classical Wigner molecule positions. Notice that in the UHF the separation of electrons is more pronounced than in the exact BS states. Fig. 2(c) shows the pair correlation function (PCF) [2] for the UHF and the exact BS state corresponding to the charge density of Fig. 2(a) with the position of one of the electrons fixed at two different locations: in the center and on the edge of the central charge puddle. In contrast to the exact BS state in the UHF wave function the two electrons are insensitive to the actual position of the third electron in its charge puddle. This is a consequence of the single-determinantal form of the UHF wave function, and can be easily explained for two electrons. In the spin-polarized two electron Wigner molecule the UHF spatial wave function is given by Ψ α (r 1 )Ψ β (r 2 )−Ψ β ...