No Arabic abstract
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.
We review the physics and phenomenology of wave dark matter: a bosonic dark matter candidate lighter than about 30 eV. Such particles have a de Broglie wavelength exceeding the average inter-particle separation in a galaxy like the Milky Way, and are well described as classical waves. We outline the particle physics motivations for them, including the QCD axion and ultra-light axion-like-particles such as fuzzy dark matter. The wave nature of the dark matter implies a rich phenomenology: (1) Wave interference leads to order unity density fluctuations on de Broglie scale. A manifestation is vortices where the density vanishes and around which the velocity circulates. There is one vortex ring per de Broglie volume on average. (2) For sufficiently low masses, soliton condensation occurs at centers of halos. The soliton oscillates and random walks, another manifestation of wave interference. The halo/subhalo abundance is suppressed at small masses, but the precise prediction from numerical wave simulations remains to be determined. (3) For ultra-light ~$10^{-22}$ eV dark matter, the wave interference substructures can be probed by tidal streams/gravitational lensing. The signal can be distinguished from that due to subhalos by the dependence on stream orbital radius/image separation. (4) Axion detection experiments are sensitive to interference substructures for moderately light masses. The stochastic nature of the waves affects the interpretation of experiments and motivates the measurement of correlation functions. Current constraints and open questions, covering detection experiments and cosmological/galactic/black-hole observations, are discussed.
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.
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.
It is widely accepted that dark matter contributes about a quarter of the critical mass-energy density in our Universe. The nature of dark matter is currently unknown, with the mass of possible constituents spanning nearly one hundred orders of magnitude. The ultralight scalar field dark matter, consisting of extremely light bosons with $m sim 10^{-22}$ eV and often called fuzzy dark matter, provides intriguing solutions to some challenges at sub-Galactic scales for the standard cold dark matter model. As shown by Khmelnitsky and Rubakov, such a scalar field in the Galaxy would produce an oscillating gravitational potential with nanohertz frequencies, resulting in periodic variations in the times of arrival of radio pulses from pulsars. The Parkes Pulsar Timing Array (PPTA) has been monitoring 20 millisecond pulsars at two to three weeks intervals for more than a decade. In addition to the detection of nanohertz gravitational waves, PPTA offers the opportunity for direct searches for fuzzy dark matter in an astrophysically feasible range of masses. We analyze the latest PPTA data set which includes timing observations for 26 pulsars made between 2004 and 2016. We perform a search in this data set for evidence of ultralight dark matter in the Galaxy using Bayesian and Frequentist methods. No statistically significant detection has been made. We therefore place upper limits on the local dark matter density. Our limits, improving on previous searches by a factor of two to five, constrain the dark matter density of ultralight bosons with $m leq 10^{-23}$ eV to be below $6,text{GeV},text{cm}^{-3}$ with 95% confidence in the Earth neighborhood. Finally, we discuss the prospect of probing the astrophysically favored mass range $m gtrsim 10^{-22}$ eV with next-generation pulsar timing facilities.
The black hole merging rates inferred after the gravitational-wave detection by Advanced LIGO/VIRGO and the relatively high mass of the progenitors are consistent with models of dark matter made of massive primordial black holes (PBH). PBH binaries emit gravitational waves in a broad range of frequencies that will be probed by future space interferometers (LISA) and pulsar timing arrays (PTA). The amplitude of the stochastic gravitational-wave background expected for PBH dark matter is calculated taking into account various effects such as initial eccentricity of binaries, PBH velocities, mass distribution and clustering. It allows a detection by the LISA space interferometer, and possibly by the PTA of the SKA radio-telescope. Interestingly, one can distinguish this background from the one of non-primordial massive binaries through a specific frequency dependence, resulting from the maximal impact parameter of binaries formed by PBH capture, depending on the PBH velocity distribution and their clustering properties. Moreover, we find that the gravitational wave spectrum is boosted by the width of PBH mass distribution, compared with that of the monochromatic spectrum. The current PTA constraints already rule out broad-mass PBH models covering more than three decades of masses, but evading the microlensing and CMB constraints due to clustering.