We study the properties of Bose-Einstein Condensate (BEC) systems consisting of two scalars, focusing on both the case where the BEC is stellar scale as well as the case when it is galactic scale. After studying the stability of such systems and making contact with existing single scalar limits, we undertake a numerical study of the two interacting scalars using Einstein-Klein-Gordon (EKG) equations, including both non-gravitational self-interactions and interactions between the species. We show that the presence of extra scalars and possible interactions between them can leave unique imprints on the BEC system mass profile, especially when the system transitions from being dominated by one scalar to being dominated by the other. At stellar scales (nonlinear regime,) we observe that a repulsive interaction between the two scalars of the type $+phi_1^2 phi_2^2$ can stabilize the BEC system and support it up to high compactness, a phenomenon only known to exist in the $+phi^4$ system. We provide simple analytic understanding of this behavior and point out that it can lead to interesting gravitational wave signals at LIGO-Virgo. At galactic scales, on the other hand, we show that two-scalar BECs can address the scaling problem that arises when one uses ultralight dark matter mass profiles to fit observed galactic core mass profiles. In the end, we construct a particle model of two ultralight scalars with the repulsive $phi_1^2 phi_2^2$ interaction using collective symmetry breaking. We develop a fast numerical code that utilizes the relaxation method to solve the EKG system, which can be easily generalized to multiple scalars.
Download