On the strength of a double Nobel prize winning experiment (Super)Kamiokande and an extremely successful long baseline neutrino programme, the third generation Water Cherenkov detector, Hyper-Kamiokande, is being developed by an international collaboration as a leading worldwide experiment based in Japan. The Hyper-Kamiokande detector will be hosted in the Tochibora mine, about 295 km away from the J-PARC proton accelerator research complex in Tokai, Japan. The currently existing accelerator will be steadily upgraded to reach a MW beam by the start of the experiment. A suite of near detectors will be vital to constrain the beam for neutrino oscillation measurements. A new cavern will be excavated at the Tochibora mine to host the detector. The experiment will be the largest underground water Cherenkov detector in the world and will be instrumented with new technology photosensors, faster and with higher quantum efficiency than the ones in Super-Kamiokande. The science that will be developed will be able to shape the future theoretical framework and generations of experiments. Hyper-Kamiokande will be able to measure with the highest precision the leptonic CP violation that could explain the baryon asymmetry in the Universe. The experiment also has a demonstrated excellent capability to search for proton decay, providing a significant improvement in discovery sensitivity over current searches for the proton lifetime. The atmospheric neutrinos will allow to determine the neutrino mass ordering and, together with the beam, able to precisely test the three-flavour neutrino oscillation paradigm and search for new phenomena. A strong astrophysical programme will be carried out at the experiment that will detect supernova neutrinos and will measure precisely solar neutrino oscillation.

In this work we combine information from relic abundance, direct detection, cosmic microwave background, positron fraction, gamma rays, and colliders to explore the existing constraints on couplings between Dark Matter and Standard Model constituents when no underlying model or correlation is assumed. For definiteness, we include independent vector-like effective interactions for each Standard Model fermion. Our results show that low Dark Matter masses below 20 GeV are disfavoured at the $3 \sigma$ level with respect to higher masses, due to the tension between the relic abundance requirement and upper constraints on the Dark Matter couplings. Furthermore, large couplings are typically only allowed in combinations which avoid effective couplings to the nuclei used in direct detection experiments.

We carry out a time-series analysis of the combined data from three experiments measuring the cosmic muon flux at the Gran Sasso laboratory, at a depth of 3800 m.w.e. These data, taken by the MACRO, LVD and Borexino experiments, span a period of over 20 years, and correspond to muons with a threshold energy, at sea level, of around 1.3 TeV. We compare the best-fit period and phase of the full muon data set with the combined DAMA/NaI and DAMA/LIBRA data, which spans the same time period, as a test of the hypothesis that the cosmic ray muon flux is responsible for the annual modulation detected by DAMA. We find in the muon data a large-amplitude fluctuation with a period of around one year, and a phase that is incompatible with that of the DAMA modulation at 5.2 sigmas. Aside from this annual variation, the muon data also contains a further significant modulation with a period between 10 and 11 years and a power well above the 99.9% C.L threshold for noise, whose phase corresponds well with the solar cycle: a surprising observation for such high energy muons. We see no corresponding long-period oscillation in the stratospheric temperature data.

The popular freeze-out paradigm for Dark Matter (DM) production, relies on DM-baryon couplings of the order of the weak interactions. However, different search strategies for DM have failed to provide a conclusive evidence of such (non-gravitational) interactions, while greatly reducing the parameter space of many representative models. This motivates the study of alternative mechanisms for DM genesis. In the freeze-in framework, the DM is slowly populated from the thermal bath while never reaching equilibrium. In this work, we analyse in detail the possibility of producing a frozen-in DM via a mediator particle which acts as a portal. We give analytical estimates of different freeze-in regimes and support them with full numerical analyses, taking into account the proper distribution functions of bath particles. Finally, we constrain the parameter space of generic models by requiring agreement with DM relic abundance observations.

Asymmetric Dark Matter (ADM) models invoke a particle-antiparticle asymmetry, similar to the one observed in the Baryon sector, to account for the Dark Matter (DM) abundance. Both asymmetries are usually generated by the same mechanism and generally related, thus predicting DM masses around 5 GeV in order to obtain the correct density. The main challenge for successful models is to ensure efficient annihilation of the thermally produced symmetric component of such a light DM candidate without violating constraints from collider or direct searches. A common way to overcome this involves a light mediator, into which DM can efficiently annihilate and which subsequently decays into Standard Model particles. Here we explore the scenario where the light mediator decays instead into lighter degrees of freedom in the dark sector that act as radiation in the early Universe. While this assumption makes indirect DM searches challenging, it leads to signals of extra radiation at BBN and CMB. Under certain conditions, precise measurements of the number of relativistic species, such as those expected from the Planck satellite, can provide information on the structure of the dark sector. We also discuss the constraints of the interactions between DM and Dark Radiation from their imprint in the matter power spectrum.

This white paper addresses the hypothesis of light sterile neutrinos based on recent anomalies observed in neutrino experiments and the latest astrophysical data.

We discuss aidnogenesis, the generation of a dark matter asymmetry via new sphaleron processes associated to an extra non-abelian gauge symmetry common to both the visible and the dark sectors. Such a theory can naturally produce an abundance of asymmetric dark matter which is of the same size as the lepton and baryon asymmetries, as suggested by the similar sizes of the observed baryonic and dark matter energy content, and provide a definite prediction for the mass of the dark matter particle. We discuss in detail a minimal realization in which the Standard Model is only extended by dark matter fermions which form "dark baryons" through an SU(3) interaction, and a (broken) horizontal symmetry that induces the new sphalerons. The dark matter mass is predicted to be approximately 6 GeV, close to the region favored by DAMA and CoGeNT. Furthermore, a remnant of the horizontal symmetry should be broken at a lower scale and can also explain the Tevatron dimuon anomaly.

Geometrical tests such as the combination of the Hubble parameter H(z) and the angular diameter distance d_A(z) can, in principle, break the degeneracy between the dark energy equation of state parameter w(z), and the spatial curvature Omega_k in a direct, model-independent way. In practice, constraints on these quantities achievable from realistic experiments, such as those to be provided by Baryon Acoustic Oscillation (BAO) galaxy surveys in combination with CMB data, can resolve the cosmic confusion between the dark energy equation of state parameter and curvature only statistically and within a parameterized model for w(z). Combining measurements of both H(z) and d_A(z) up to sufficiently high redshifts around z = 2 and employing a parameterization of the redshift evolution of the dark energy equation of state are the keys to resolve the w(z)-Omega_k degeneracy.

We generalize to non-flat geometries the formalism of Simon et al. (2005) to reconstruct the dark energy potential. This formalism makes use of quantities similar to the Horizon-flow parameters in inflation, can, in principle, be made non-parametric and is general enough to be applied outside the simple, single scalar field quintessence. Since presently available and forthcoming data do not allow a non-parametric and exact reconstruction of the potential, we consider a general parametric description in term of Chebyshev polynomials. We then consider present and future measurements of H(z), Baryon Acoustic Oscillations surveys and Supernovae type 1A surveys, and investigate their constraints on the dark energy potential. We find that, relaxing the flatness assumption increases the errors on the reconstructed dark energy evolution but does not open up significant degeneracies, provided that a modest prior on geometry is imposed. Direct measurements of H(z), such as those provided by BAO surveys, are crucially important to constrain the evolution of the dark energy potential and the dark energy equation of state, especially for non-trivial deviations from the standard LambdaCDM model.