At large scales and for sufficiently early times, dark matter is described as a pressureless perfect fluid---dust---non-interacting with Standard Model fields. These features are captured by a simple model with two scalars: a Lagrange multiplier and another playing the role of the velocity potential. That model arises naturally in some gravitational frameworks, e.g., the mimetic dark matter scenario. We consider an extension of the model by means of higher derivative terms, such that the dust solutions are preserved at the background level, but there is a non-zero sound speed at the linear level. We associate this {\it Modified Dust} with dark matter, and study the linear evolution of cosmological perturbations in that picture. The most prominent effect is the suppression of their power spectrum for sufficiently large cosmological momenta. This can be relevant in view of the problems that cold dark matter faces at sub-galactic scales, e.g., the missing satellites problem. At even shorter scales, however, perturbations of Modified Dust are enhanced compared to the predictions of more common particle dark matter scenarios. This is a peculiarity of their evolution in radiation dominated background. We also briefly discuss clustering of Modified Dust. We write the system of equations in the Newtonian limit, and sketch the possible mechanism which could prevent the appearance of caustic singularities. The same mechanism may be relevant in light of the core-cusp problem.

We simulate the adiabatic contraction of a dark matter (DM) distribution during the process of the star formation, paying particular attention to the phase space distribution of the DM particles after the contraction. Assuming the initial uniform density and Maxwellian distribution of DM velocities, we find that the number $n(r)$ of DM particles within the radius $r$ scales like $n(r) \propto r^{1.5}$, leading to the DM density profile $\rho\propto r^{-1.5}$, in agreement with the Liouville theorem and previous numerical studies. At the same time, the number of DM particles $\nu(r)$ with periastra smaller than $r$ is parametrically larger, $\nu(r) \propto r$, implying that many particles contributing at any given moment into the density $\rho(r)$ at small $r$ have very elongated orbits and spend most of their time at distances larger than $r$. This has implications for the capture of DM by stars in the process of their formation. As a concrete example we consider the case of primordial black holes (PBH). We show that accounting for very eccentric orbits boosts the amount of captured PBH by a factor of up to $2\times 10^3$ depending on the PBH mass, improving correspondingly the previously derived constraints on the PBH abundance.

In arXiv:1401.3025 a new mechanism of energy loss by primordial black holes passing through neutron stars has been proposed, which is more efficient, by many orders of magnitude, than the existing ones. In this comment we point out a problem in the calculations of arXiv:1401.3025 that may invalidate the result.

We investigate constraints on primordial black holes (PBHs) as dark matter candidates that arise from their capture by neutron stars (NSs). If a PBH is captured by a NS, the star is accreted onto the PBH and gets destroyed in a very short time. Thus, mere observations of NSs put limits on the abundance of PBHs. High DM densities and low velocities are required to constrain the fraction of PBHs in DM. Such conditions may be realized in the cores of globular clusters if the latter are of a primordial origin. Assuming that cores of globular clusters possess the DM densities exceeding several hundred GeV/cm$^3$ would imply that PBHs are excluded as comprising all of the dark matter in the mass range $3\times 10^{18} \text{g} \lesssim m_\text{BH}\lesssim 10^{24} \text{g}$. At the DM density of $2\times 10^3$ GeV/cm$^3$ that has been found in simulations in the corresponding models, less than 5% of the DM may consist of PBH for these PBH masses.

By considering adiabatic contraction of the dark matter (DM) during star formation, we estimate the amount of DM trapped in stars at their birth. If the DM consists partly of primordial black holes (PBHs), they will be trapped together with the rest of the DM and will be finally inherited by a star compact remnant --- a white dwarf (WD) or a neutron star (NS), which they will destroy in a short time. Observations of WDs and NSs thus impose constraints on the abundance of PBH. We show that the best constraints come from WDs and NSs in globular clusters which exclude the DM consisting entirely of PBH in the mass range $10^{16}{\rm g} - 3\times 10^{22}{\rm g}$, with the strongest constraint on the fraction $\Omega_{\rm PBH} /\Omega_{\rm DM}\lesssim 10^{-2}$ being in the range of PBH masses $10^{17}{\rm g} - 10^{18}$ g.