Abstract. -We study electron molecules in realistic vertically coupled quantum dots in a strong magnetic field. Computing the energy spectrum, pair correlation functions, and dynamical form factor as a function of inter-dot coupling via diagonalization of the many-body Hamiltonian, we identify structural transitions between different phases, some of which do not have a classical counterpart. The calculated Raman cross section shows how such phases can be experimentally singled out.Electron systems form a Wigner crystal at sufficiently low density or high magnetic field B [1]. Theoretical [2,3] and experimental [4,5] studies suggest that lowering dimensionality favors localization: in this perspective interacting electrons confined in a quantum dot (QD) [6], sometimes called Wigner molecules [7], are interesting in their own right [8,9], due to the interplay between the electron-electron repulsion and the confining potential. This leads to a complex zero temperature phase diagram [10], as compared to the infinite layer case, as well as to complex melting mechanisms [9]. The formation of coupled QDs (artificial molecules) introduces qualitatively new physics [11]. New energy scales appear -inter-dot tunneling, inter-vs intra-dot Coulomb correlation -, whose balance controls the phase diagram [12]. Significantly, these parameters can be tuned by inter-dot distance d and/or electron density, so that the nature of these few-particle systems and their phases can be explored experimentally.In this Letter we discuss quantum mechanical calculations of N electrons in a coupled QD structure in the strong field regime, where localization is ensured in the parent isolated QDs. Monitoring the spatial correlation functions, we identify different ground states depending on the inter-dot coupling. More interestingly, we find that some of the phases do not have a classical counterpart [13], and are ascribed to a two-three-two dimensional (2D-3D-2D) transition of the electronic system. Such phases were not identified in previous studies of the coupled layer system due to the neglect of tunneling and/or finite width of the layers [14-17], i.e., of the 3D nature of the system, which turns out to be essential for their formation. From the analysis of the dynamical form factor, we associate to each phase peculiar collective charge