No Arabic abstract
Recently, we have shown how current cosmological N-body codes already follow the fine grained phase-space information of the dark matter fluid. Using a tetrahedral tesselation of the three-dimensional manifold that describes perfectly cold fluids in six-dimensional phase space, the phase-space distribution function can be followed throughout the simulation. This allows one to project the distribution function into configuration space to obtain highly accurate densities, velocities, and velocity dispersions. Here, we exploit this technique to show first steps on how to devise an improved particle-mesh technique. At its heart, the new method thus relies on a piecewise linear approximation of the phase space distribution function rather than the usual particle discretisation. We use pseudo-particles that approximate the masses of the tetrahedral cells up to quadrupolar order as the locations for cloud-in-cell (CIC) deposit instead of the particle locations themselves as in standard CIC deposit. We demonstrate that this modification already gives much improved stability and more accurate dynamics of the collisionless dark matter fluid at high force and low mass resolution. We demonstrate the validity and advantages of this method with various test problems as well as hot/warm-dark matter simulations which have been known to exhibit artificial fragmentation. This completely unphysical behaviour is much reduced in the new approach. The current limitations of our approach are discussed in detail and future improvements are outlined.
If dark matter is mainly composed of axions, the density distribution can be nonuniformly distributed, being clumpy instead. By solving the Einstein-Klein-Gordon system of a scalar field with the potential energy density of an axionlike particle, we obtain the maximum mass of the self-gravitating system made of axions, called axion stars. The collision of axion stars with neutron stars may release the energy of axions due to the conversion of axions into photons in the presence of the neutron stars magnetic field. We estimate the energy release and show that it should be much less than previous estimates.Future data from femtolensing should strongly constrain this scenario.
Using $N$-body simulations ($Nsim 10^6 - 10^7$), we examine how a non-axisymmetric dark halo affects the dynamical evolution of the structure in collisionless (stellar) discs. We demonstrate how the model parameters such as mass of the halo, initial conditions in the disc and the halo axes ratio affect morphology and kinematics of the stellar discs. We show that a non-axisymmetric halo can generate a large-scale spiral density pattern in the embedded stellar disc. The pattern is observed in the disc for many periods of its revolution, even if the disc is gravitationally over-stable. The growth of the spiral arms is not accompanied by significant dynamical heating of the disc, irrelevant to its initial parameters. We also investigate transformation of the dark halos shape driven by the long-lived spiral pattern in the disc . We show that the analysis of the velocity field in the stellar disc and in the spiral pattern gives us a possibility to figure out the spatial orientation of the triaxial-shaped dark halo and to measure the triaxiality.
We consider a cosmological scenario where the dark sector is described by two perfect fluids that interact through a velocity-dependent coupling. This coupling gives rise to an interaction in the dark sector driven by the relative velocity of the components, thus making the background evolution oblivious to the interaction and only the perturbed Euler equations are affected at first order. We obtain the equations governing this system with the Schutz-Sorkin Lagrangian formulation for perfect fluids and derive the corresponding stability conditions to avoid ghosts and Laplacian instabilities. As a particular example, we study a model where dark energy behaves as a radiation fluid at high redshift while it effectively becomes a cosmological constant in the late Universe. Within this scenario, we show that the interaction of both dark components leads to a suppression of the dark matter clustering at late times. We also argue the possibility that this suppression of clustering together with the additional dark radiation at early times can simultaneously alleviate the $sigma_8$ and $H_0$ tensions.
Directional detection can provide unambiguous observation of Dark Matter interactions even in presence of insidious backgrounds. The DM-TPC collaboration is developing a detector with the goal of measuring the direction and sense of nuclear recoils produced in Dark Matter interactions. The detector consists of a Time Projection Chamber with optical readout filled with CF$_4$ gas at low pressure. A collision between a WIMP and a gas molecule results in a nuclear recoil of 1-2 mm. The measurement of the energy loss along the recoil allows us to determine the sense and the direction of the recoil. Results from a prototype detector operated in a low-energy neutron beam clearly demonstrate the suitability of this approach to measure directionality. A full-scale module with an active volume of about one cubic meter is now being designed. This detector, which will be operated underground in 2009, will allow us to set competitive limits on spin-dependent Dark Matter interactions using a directional detector.
The cosmic neutrino background is an important component of the Universe that is difficult to include in cosmological simulations due to the extremely large velocity dispersion of neutrino particles. We develop a new approach to simulate cosmic neutrinos that decomposes the Fermi-Dirac phase space into shells of constant speed and then evolves those shells using hydrodynamic equations. These collisionless hydrodynamic equations are chosen to match linear theory, free particle evolution and allow for superposition. We implement this method into the information-optimized cosmological $N$-body code CUBE and demonstrate that neutrino perturbations can be accurately resolved to at least $ksim1 h/$Mpc. This technique allows for neutrino memory requirements to be decreased by up to $sim 10^3$ compared to traditional $N$-body methods.