We consider a dark matter halo (DMH) of a spherical galaxy as a Bose-Einstein condensate of the ultra-light axions interacting with the baryonic matter. In the mean-field limit, we have derived the integro-differential equation of the Hartree-Fock type for the spherically symmetrical wave function of the DMH component. This equation includes two independent dimensionless parameters: (i) \b{eta}- the ratio of baryon and axion total mases and (ii) {\xi}- the ratio of characteristic baryon and axion spatial parameters. We extended our dissipation algorithm for studying numerically the ground state of the axion halo in the gravitational field produced by the baryonic component. We calculated the characteristic size, Xc, of DMH as a function of \b{eta} and {\xi} and obtained an analytical approximation for Xc.

We consider a two-component dark matter halo (DMH) of a galaxy containing ultra-light axions (ULA) of different mass. The DMH is described as a Bose-Einstein condensate (BEC) in its ground state. In the mean-field (MF) limit we have derived the integro-differential equations for the spherically symmetrical wave functions of the two DMH components. We studied numerically the radial distribution of the mass density of ULA and constructed the parameters which could be used to distinguish between the two- and one-component DMH. We also discuss an interesting connection between the BEC MF ground state of a one-component DMH and Hawking temperature and entropy of a corresponding Black Hole, and Unruh temperature.