A simple model in which both electrons of the shallow-negative donor D ± ± are instantaneously localized approximately at the same distance with respect to the ion nucleus is considered. Within this model, the D ± ± Hamiltonian becomes separable, and the effect of the electron±electron interaction is essentially taken into account by the interaction of each of the electrons with a renormalized-nucleus effective charge which depends on the angle between the two electron vector positions. The trial function is taken as a product of two three-parameter one-particle wave functions to calculate the ground state energy as a function of the angle. Theoretical results for the D ± ± energy, including zero-point energy corrections within the harmonic approximation, are obtained for D ± ± states centered in GaAs±(Ga,Al)As quantum wells as a function of the well width, and for different values of the magnetic field, applied perpendicular to the interfaces of the heterostructure. Our results are in good agreement with experimental data and with theoretical calculations using variational or Monte Carlo techniques.Low-dimensional semiconductor microstructures with neutral and charged impurities and impurity complexes bring about new phenomena with a lot of challenges to both theoretical and experimental physics. Shallow-negative donors D ± ± (i.e. neutral shallow donors that bind an additional electron) have been experimentally observed [1,2] in GaAs±(Ga,Al)As quantum wells (QWs). It has been also demonstrated [3,4] by variational calculations that the two-dimensional D ± ± binding energy is about ten times larger than the three-dimensional one, and the impurity mean-value radius is reduced by half when its form is changed from spherical to planar. For this reason, optical excitations of the D ± ± complex may fall in the same spectral region [5] of the neutral D 0 state in a QW. The identification of the D ± ± features in a given experimental spectrum is, therefore, not straightforward and it suggests the need of a careful theoretical analysis. Any calculation of the D ± ± spectrum has to take into account the electronic correlation, which is responsible for the system stability. Variational [4 to 6] and Monte Carlo [7] methods have been used to find the D ± ± binding energy in a QW. Both types of calculations require a set of adjustable parameters and quite long computational times (for example, the Chandrasekhar-type [8] trial function with seven variational parameters [6] has been used). Due to the computational cost of such calculations, no detailed studies of the evolution of the D ± ± spectrum as a function of well width, barrier height, donor position, etc. were performed. On the other hand, in few-particle systems in