An empirical potential with fluctuating charges is proposed for modelling (MgO)n clusters in both the molecular (small n) and bulk (n → ∞) regimes. Vectorial polarization forces are explicitely taken into account in the self-consistent determination of the charges. Our model predicts cuboid cluster structures, in agreement with previous experimental and theoretical results. The effective charge transferred between magnesium and oxygen smoothly increases from 1 to 2, with an estimated crossover size above 300 MgO molecules.PACS numbers: 36.40. Cg,36.40.Ei Crystalline magnesium oxyde is a purely ionic compound in which the Mg Z+ and O Z− ions carry a charge Z around ±2.1 In gas phase, the oxyde anion O 2− is unstable and spontaneously decays into O − +e due to the strong electron-electron repulsion. As a result, the effective atomic charge in the MgO molecule is much smaller than 2: independent ab initio calculations by Ziemann and Castleman 2 and by Recio et al.3 found Z ∼ 0.8. In the intermediate sizes regime, (MgO) n clusters are thus expected to show intriguing properties due to a partially covalent character of the chemical bonding. Beyond condensed matter or molecular physics, these clusters received some special attention in the astrophysics community, where they have been involved in the nucleation process of dust in circumstellar shells around M-stars.Despite the vast amount of experimental 2,4,5 and theoretical 2,3,6,7,8,9,10,11,12 investigations on neutral or charged clusters, the way and the rate at which chemical bonding evolves from ionic/covalent at small sizes toward purely ionic in the bulk remains essentially unexplored. Because the electrostatic field created by the ions does not vanish in finite systems, the highly polarizable oxyde anion has a rather floppy and deformable outer electron cloud, which could be responsible for a partial screening of the repulsion between cations. However, the situation is complicated by the possible coordination-dependence of the charge transferred locally.Theoretical studies of (MgO) n clusters can be effectively separated into two groups. Ab initio or densityfunctional theory (DFT) based calculations have been performed on specific geometries, in a rather limited size range. 2,3,6,7,8,9 These works predict that small clusters exhibit cuboid-like shapes similar to NaCl rocksalt clusters. The apparent charge transferred, as estimated from Mulliken populations, is indeed size-and coordination-dependent, and lies between 1 and 1.5 for 2 ≤ n ≤ 13.3,9 More empirical methods have also been used to predict optimal structures. None of these empirical potentials account for the different charges transferred in MgO clusters. Only in Ref. 10 the authors explicitely employed a size-dependent value of the charge Z, using an arbitrary law Z(n) = (2ζn + 1)/(ζn + 1). ζ was taken such that the crossover n * between ionic/covalent and purely ionic, for which Z equals 1.5, occurs approximately at n * = 20. However, the low energy structures found by Köhler and coworkers signifi...