Cosmological constraints on the multiscalar field dark matter model

The main aim of this paper is to provide cosmological constraints on the multiscalar field dark matter model, in which we assume the dark matter is made up of different ultralight scalar fields. As a first approximation, we consider they are real and do not interact with each other. We study the equations for both the background and perturbations for N-fields and present the evolution of the density parameters, the mass power spectrum, and the cosmic microwave background spectrum. In particular, we focus on two scalar fields with several combinations for the potentials V () =1/2m_^2^2, V () =m_^2f^21+cos (/f), and V () =m_^2f^2cosh (/f) -1, however, the work, along with the code, could be easily extended to more fields. We use the data from baryon acoustic oscillation, big bang nucleosynthesis, Lyman- forest, and supernovae to find constraints on the sampling parameters for the cases of a single field and double field, along with the Bayesian evidence. We found that some combinations of the potentials get penalized through the evidence, however, for others there is a preference as good as for the cold dark matter.