A numerical study of the Schrodinger-Newton equations

In this paper we present a numerical study of the Schrodinger-Newton (SN) equations. We begin by considering the linear stability of the spherically symmetric stationary solutions found numerically by Moroz et al (1998 Class. Quantum Grav. 15 2733-42) and Bernstein et al (1998 Mod. Phys. Lett. A13 2327-36). The ground state, characterized as the state of lowest energy, turns out to he linearly stable, with only imaginary eigenvalues. The (n + 1)th state is linearly unstable having n quadruples of complex eigenvalues (as well as imaginary eigenvalues), Where a quadruple consists of {lambda,(λ) over bar, -lambda, -(λ) over bar} for complex lambda. Then, we consider the time-dependent SN equations in three-dimensions with three kinds of symmetry: spherically symmetric, axially symmetric and translationally symmetric. We find that the solutions show a balance between the dispersive tendencies of the Schrodinger equation and the gravitional attraction from the Poisson equation. Only the ground state is stable, and lumps of probability attract each other gravitationally before dispersing.