No Arabic abstract
Wave dark matter ($psi$DM), which satisfies the Schrodinger-Poisson equation, has recently attracted substantial attention as a possible dark matter candidate. Numerical simulations have in the past provided a powerful tool to explore this new territory of possibility. Despite their successes to reveal several key features of $psi$DM, further progress in simulations is limited, in that cosmological simulations so far can only address formation of halos below $sim 2times 10^{11} M_odot$ and substantially more massive halos have become computationally very challenging to obtain. For this reason, the present work adopts a different approach in assessing massive halos by constructing wave-halo solutions directly from the wave distribution function. This approach bears certain similarity with the analytical construction of particle-halo (cold dark matter model). Instead of many collisionless particles, one deals with one single wave that has many non-interacting eigenstates. The key ingredient in the wave-halo construction is the distribution function of the wave power, and we use several halos produced by structure formation simulations as templates to determine the wave distribution function. Among different models, we find the fermionic King model presents the best fits and we use it for our wave-halo construction. We have devised an iteration method for constructing the nonlinear halo, and demonstrate its stability by three-dimensional simulations. A Milky-Way-sized halo has also been constructed, and the inner halo is found flatter than the NFW profile. These wave-halos have small-scale interferences both in space and time producing time-dependent granules. While the spatial scale of granules varies little, the correlation time is found to increase with radius by one order of magnitude across the halo.
We present a wave generalization of the classic Schwarzschild method for constructing self-consistent halos -- such a halo consists of a suitable superposition of waves instead of particle orbits, chosen to yield a desired mean density profile. As an illustration, the method is applied to spherically symmetric halos. We derive an analytic relation between the particle distribution function and the wave superposition amplitudes, and show how it simplifies in the high energy (WKB) limit. We verify the stability of such constructed halos by numerically evolving the Schrodinger-Poisson system. The algorithm provides an efficient and accurate way to simulate the time-dependent halo substructures from wave interference. We use this method to construct halos with a variety of density profiles, all of which have a core from the ground-state wave function, though the core-halo relation need not be the standard one.
This papers explores the self similar solutions of the Vlasov-Poisson system and their relation to the gravitational collapse of dynamically cold systems. Analytic solutions are derived for power law potential in one dimension, and extensions of these solutions in three dimensions are proposed. Next the self similarity of the collapse of cold dynamical systems is investigated numerically. The fold system in phase space is consistent with analytic self similar solutions, the solutions present all the proper self-similar scalings. An additional point is the appearance of an $x^{-(1/2)}$ law at the center of the system for initial conditions with power law index larger than $-(1/2)$. It is found that the first appearance of the $x^{-(1/2)}$ law corresponds to the formation of a singularity very close to the center. Finally the general properties of self similar multi dimensional solutions near equilibrium are investigated. Smooth and continuous self similar solutions have power law behavior at equilibrium. However cold initial conditions result in discontinuous phase space solutions, and the smoothed phase space density looses its auto similar properties. This problem is easily solved by observing that the probability distribution of the phase space density $P$ is identical except for scaling parameters to the probability distribution of the smoothed phase space density $P_S$. As a consequence $P_S$ inherit the self similar properties of $P$. This particular property is at the origin of the universal power law observed in numerical simulation for ${rho}/{sigma^3}$. The self similar properties of $P_S$ implies that other quantities should have also an universal power law behavior with predictable exponents. This hypothesis is tested using a numerical model of the phase space density of cold dark matter halos, an excellent agreement is obtained.
Galactic disks in triaxial dark matter halos become deformed by the elliptical potential in the plane of the disk in such a way as to counteract the halo ellipticity. We develop a technique to calculate the equilibrium configuration of such a disk in the combined disk-halo potential, which is based on the method of Jog (2000) but accounts for the radial variation in both the halo potential and the disk ellipticity. This crucial ingredient results in qualitatively different behavior of the disk: the disk circularizes the potential at small radii, even for a reasonably low disk mass. This effect has important implications for proposals to reconcile cuspy halo density profiles with low surface brightness galaxy rotation curves using halo triaxiality. The disk ellipticities in our models are consistent with observational estimates based on two-dimensional velocity fields and isophotal axis ratios.
Wave Dark Matter (WaveDM) has recently gained attention as a viable candidate to account for the dark matter content of the Universe. In this paper we explore the extent to which dark matter halos in this model, and under what conditions, are able to reproduce strong lensing systems. First, we analytically explore the lensing properties of the model -- finding that a pure WaveDM density profile, a soliton profile, produces a weaker lensing effect than other similar cored profiles. Then we analyze models with a soliton embedded in an NFW profile, as has been found in numerical simulations of structure formation. We use a benchmark model with a boson mass of $m_a=10^{-22}{rm eV}$, for which we see that there is a bi-modality in the contribution of the external NFW part of the profile, and actually some of the free parameters associated with it are not well constrained. We find that for configurations with boson masses $10^{-23}$ -- $10^{-22}{rm eV}$, a range of masses preferred by dwarf galaxy kinematics, the soliton profile alone can fit the data but its size is incompatible with the luminous extent of the lens galaxies. Likewise, boson masses of the order of $10^{-21}{rm eV}$, which would be consistent with Lyman-$alpha$ constraints and consist of more compact soliton configurations, necessarily require the NFW part in order to reproduce the observed Einstein radii. We then conclude that lens systems impose a conservative lower bound $m_a > 10^{-24}$ and that the NFW envelope around the soliton must be present to satisfy the observational requirements.
We have constructed realistic, self-consistent models of triaxial elliptical galaxies embedded in triaxial dark matter halos. Self-consistent solutions by means of the standard orbital superposition technique introduced by Schwarzschild were found in each of the three cases studied. Chaotic orbits were found to be important in all of the models, and their presence was shown to imply a possible slow evolution of the shapes of the halos. The equilibrium velocity distribution is reproduced by a Lorentzian function better than by a Gaussian. Our results demonstrate for the first time that triaxial dark matter halos can co-exist with triaxial galaxies.