The formation and evolution of superdense clumps (or subhalos) is studied. Such clumps of dark matter (DM) can be produced by many mechanisms, most notably by spiky features in the spectrum of inflationary perturbations and by cosmological phase transitions. Being produced very early during the radiation dominated epoch, superdense clumps evolve as isolated objects. They do not belong to hierarchical structures for a long time after production, and therefore they are not destroyed by tidal interactions during the formation of larger structures. For DM particles with masses close to the electroweak (EW) mass scale, superdense clumps evolve towards a power-law density profile $\rho(r) \propto r^{-1.8}$ with a central core. Superdense clumps cannot be composed of standard neutralinos, since their annihilations would overproduce the diffuse gamma radiation. If the clumps are constituted of superheavy DM particles and develop a sufficiently large central density, the evolution of their central part can lead to a 'gravithermal catastrophe.' In such a case, the initial density profile turns into an isothermal profile with $\rho \propto r^{-2}$ and a new, much smaller core in the center. Superdense clumps can be bserved by gamma radiation from DM annihilations and by gravitational wave detectors, while the production of primordial black holes and cascade nucleosynthesis constrain this scenario.

Superheavy dark matter (SHDM) exchanges energy with its environment much slower than particles with masses close to the electroweak (EW) scale and has therefore different small-scale clustering properties. Using the neutralino as candidate for the SHDM, we find that free-streaming allows the formation of DM clumps of all masses down to $\sim 260 m_\chi$ in the case of bino. If small-scale clumps evolve from a non-standard, spiky spectrum of perturbations, DM clumps may form during the radiation dominated era. These clumps are not destroyed by tidal interactions and can be extremely dense. In the case of a bino, a "gravithermal catastrophe" can develop in the central part of the most dense clumps, increasing further the central density and thus the annihilation signal. In the case of a higgsino, the annihilation signal is enhanced by the Sommerfeld effect. As a result annihilations of superheavy neutralinos in dense clumps may lead to observable fluxes of annihilation products in the form of ultrahigh energy particles, for both cases, higgsinos and binos, as lightest supersymmetric particles.

Interpretations of indirect searches for dark matter (DM) require theoretical predictions for the annihilation or decay rates of DM into stable particles of the standard model. These predictions include usually only final states accessible as lowest order tree-level processes, with electromagnetic bremsstrahlung and the loop-suppressed two gamma-ray line as exceptions. We show that this restriction may lead to severely biased results for DM tailored to produce only leptons in final states and with mass in the TeV range. For such models, unavoidable electroweak bremsstrahlung of Z and W-bosons has a significant influence both on the branching ratio and the spectral shape of the final state particles. We work out the consequences for two situations: Firstly, the idealized case where DM annihilates at tree level with 100% branching ratio into neutrinos. For a given cross section, this leads eventually to "minimal yields" of photons, electrons, positrons and antiprotons. Secondly, the case where the only allowed two-body final states are electrons. The latter case is typical of models aimed at fitting cosmic ray e^- and e^+ data. We find that the multimessenger signatures of such models can be significantly modified with respect to results presented in the literature.