This white paper summarizes the activities of the Brazilian community concerning dark matter physics and highlights the importance of financial support to Brazilian groups that are deeply involved in experimental endeavours. The flagships of the Brazilian dark matter program are the Cherenkov Telescope Array, DARKSIDE, SBN and LHC experiments, but we emphasize that smaller experiments such as DAMIC and CONNIE constitute important probes to dark sectors as well and should receive special attention. Small experimental projects showing the potential to probe new regions of parameter space of dark matter models are encouraged. On the theoretical and phenomenological side, some groups are devoted to astrophysical aspects such as the dark matter density profile while others explore the signature of dark matter models at colliders, direct and indirect detection experiments. In summary, the Brazilian dark matter community that was born not long ago has grown tremendously in the past years and now plays an important role in the hunt for a dark matter particle.

The recent detection of the cosmic microwave background polarimeter experiment BICEP2 of tensor fluctuations in the B-mode power spectrum basically excludes all plausible axion models where its decay constant is above $10^{13}$ GeV. Moreover, there are strong theoretical, astrophysical, and cosmological motivations for models involving, in addition to the axion, also axion-like particles (ALPs), with decay constants in the intermediate scale range, between $10^9$ GeV and $10^{13}$ GeV. Here, we present a general analysis of models with an axion and further ALPs and derive bounds on the relative size of the axion and ALP photon (and electron) coupling. We discuss what we can learn from measurements of the axion and ALP photon couplings about the fundamental parameters of the underlying ultraviolet completion of the theory. For the latter we consider extensions of the Standard Model in which the axion and the ALP(s) appear as pseudo Nambu-Goldstone bosons from the breaking of global chiral $U(1)$ (Peccei-Quinn (PQ)) symmetries, occuring accidentally as low energy remnants from exact discrete symmetries. In such models, the axion and the further ALP are protected from disastrous explicit symmetry breaking effects due to Planck-scale suppressed operators. The scenarios considered exploit heavy right handed neutrinos getting their mass via PQ symmetry breaking and thus explain the small mass of the active neutrinos via a seesaw relation between the electroweak and an intermediate PQ symmetry breaking scale. We show some models that can accommodate simultaneously an axion dark matter candidate, an ALP explaining the anomalous transparency of the universe for $\gamma$-rays, and an ALP explaining the recently reported 3.55 keV gamma line from galaxies and clusters of galaxies, if the respective decay constants are of intermediate scale.