Wave dark matter ($\psi$DM) predicts a compact soliton core and a granular halo in every galaxy. This work presents the first simulation study of an elliptical galaxy by including both stars and $\psi$DM, focusing on the systematic changes of the central soliton and halo granules. With the addition of stars in the inner halo, we find the soliton core consistently becomes more prominent by absorbing mass from the host halo than that without stars, and the halo granules become "non-isothermal", "hotter" in the inner halo and "cooler" in the outer halo, as opposed to the isothermal halo in pure $\psi$DM cosmological simulations. Moreover, the composite (star+$\psi$DM) mass density is found to follow a $r^{-2}$ isothermal profile near the half-light radius in most cases. Most striking is the velocity dispersion of halo stars that increases rapidly toward the galactic center by a factor of at least 2 inside the half-light radius caused by the deepened soliton gravitational potential, a result that compares favorably with observations of elliptical galaxies and bulges in spiral galaxies. However in some rare situations we find a phase segregation turning a compact distribution of stars into two distinct populations with high and very low velocity dispersions; while the high-velocity component mostly resides in the halo, the very low-velocity component is bound to the interior of the soliton core, resembling stars in faint dwarf spheroidal galaxies.

We examine the nonlinear structure of gravitationally collapsed objects that form in our simulations of wavelike cold dark matter ($\psi$DM), described by the Schr\"{o}dinger-Poisson (SP) equation with a particle mass $\sim 10^{-22} {\rm eV}$. A distinct gravitationally self-bound solitonic core is found at the center of every halo, with a profile quite different from cores modeled in the warm or self-interacting dark matter scenarios. Furthermore, we show that each solitonic core is surrounded by an extended halo composed of large fluctuating dark matter granules which modulate the halo density on a scale comparable to the diameter of the solitonic core. The scaling symmetry of the SP equation and the uncertainty principle tightly relate the core mass to the halo specific energy, which, in the context of cosmological structure formation, leads to a simple scaling between core mass ($M_c$) and halo mass ($M_h$), $M_c \propto a^{-1/2} M_h^{1/3}$, where $a$ is the cosmic scale factor. We verify this scaling relation by (i) examining the internal structure of a statistical sample of virialized halos that form in our 3D cosmological simulations, and by (ii) merging multiple solitons to create individual virialized objects. Sufficient simulation resolution is achieved by adaptive mesh refinement and graphic processing units acceleration. From this scaling relation, present dwarf satellite galaxies are predicted to have kpc sized cores and a minimum mass of $\sim 10^8 {M_\odot}$, capable of solving the small-scale controversies in the cold dark matter model. Moreover, galaxies of $2\times10^{12} {M_\odot}$ at $z=8$ should have massive solitonic cores of $\sim 2\times10^9 {M_\odot}$ within $\sim 60 {\rm pc}$. Such cores can provide a favorable local environment for funneling the gas that leads to the prompt formation of early stellar spheroids and quasars.

An alternative bosonic dark matter model is examined in detail via high-resolution simulations. These bosons have particle mass of order $10^{-22}eV$ and are non-interacting. If they do exist and can account for structure formation, these bosons must be condensed into the Bose-Einstein state and described by a coherent wave function. This matter, also known as Fuzzy Dark Matter (Hu, Barkana & Gruzinov 2000),, is speculated to be able, first, to eliminate the sub-galactic halos to solve the problem of over-abundance of dwarf galaxies, and, second, to produce flat halo cores in galaxies suggested by some observations. Our simulation results show that although this extremely light bosonic dark matter indeed suppresses low-mass halos, it can, to the contrary of expectation, yield singular halo cores. The density profile of the singular halo is almost identical to the halo profile of Navarro, Frenk & White (1997). Such a profile seems to be universal, in that it can be produced via either accretion or merger.

Mounting evidence from $x$-ray observations reveals that bound objects should receive some form of energy in the past injected from non-gravitaional sources. We report that an instantaneous heating scheme, for which gases in dense regions were subjected to a temperature jump of 1keV at $z=2$ whereas those in rarified regions remained intact, can produce bound objects obeying the observed mass-temperature and luminosity-temperature relations. Such preheating lowers the peak Sunyaev-Zeldovich (SZ) power by a factor of 2 and exacerbates the need for the normalization of matter fluctuations $\sigma_8$ to assume an extreme high value $(\sim 1.1)$ for the SZ signals to account for the excess anisotropy on 5-arcminute scale detected by the Cosmic Background Imager in the cosmic microwave background radiation.