Numerical evidence for `multi-scalar stars'
  We present a class of general relativistic soliton-like solutions composed of multiple minimally coupled, massive, real scalar fields which interact only through the gravitational field. We describe a two-parameter family of solutions we call ``phase-shifted boson stars'' (parameterized by central density rho_0 and phase delta), which are obtained by solving the ordinary differential equations associated with boson stars and then altering the phase between the real and imaginary parts of the field. These solutions are similar to boson stars as well as the oscillating soliton stars found by Seidel and Suen [E. Seidel and W.M. Suen, Phys. Rev. Lett. 66, 1659 (1991)]; in particular, long-time numerical evolutions suggest that phase-shifted boson stars are stable. Our results indicate that scalar soliton-like solutions are perhaps more generic than has been previously thought. 
Scott H. Hawley, Matthew W. Choptuik