Abstract. -Numerically exact path-integral Monte Carlo data are presented for N ≤ 10 strongly interacting electrons confined in a 2D parabolic quantum dot, including a defect to break rotational symmetry. Low densities are studied, where an incipient Wigner molecule forms. A single impurity is found to cause drastic effects: (1) The standard shell-filling sequence with magic numbers N = 4, 6, 9, corresponding to peaks in the addition energy ∆(N ), is destroyed, with a new peak at N = 8, (2) spin gaps decrease, (3) for N = 8, sub-Hund's rule spin S = 0 is induced, and (4) spatial ordering of the electrons becomes rather sensitive to spin. We also comment on the recently observed bunching phenomenon.During the past few years, the electronic properties of 2D quantum dots have been the subject of intense theoretical studies [1,2]. A particularly interesting aspect of such artificial atoms is the wide experimental tunability of the electron number N and the Brueckner interaction strength parameter r s in high-quality semiconductor dots [3,4,6]. In particular, it has been established both theoretically [1,7] and experimentally [8,9] that a "Wigner molecule" of crystallized electrons forms already at surprisingly high densities corresponding to r s > ∼ 2 [10]. In this paper, we present numerically exact path-integral Monte Carlo (PIMC) [11] simulation results for N ≤ 10 electrons in the most challenging incipient Wigner molecule regime, choosing r s ≈ 4. Due to the breakdown of effective single-particle descriptions, this parameter regime is notoriously difficult to treat within approximation schemes such as Hartree-Fock theory [1], density functional theory [1,12,13], semiclassical [14] or classical [15] approaches. Especially concerning disordered dots, in the absence of benchmark calculations, their accuracy has so far been largely unclear.The incipient Wigner molecule regime also exhibits unexpected and interesting physics, in particular rather dramatic effects when including just one impurity in order to break rotational symmetry. Below we shall demonstrate this in three different ways. First, despite the presence of strong interactions, the clean system still exhibits shell structure, with exceptional stability of the dot for "magic numbers" N = 4, 6, and 9. The stability is reflected in peaks in the corresponding addition energy, ∆(N ) = E(N + 1) − 2E(N ) + E(N − 1),c EDP Sciences