It has been recently suggested \cite{Bezrukov:2017ike,Bezrukov:2018wvd} that a cosmic scalar field can completely change the keV-scale sterile neutrino production in the early Universe. Its effect may, for various parameter choices, either suppress sterile neutrino production and make moderate active-sterile mixing cosmologically acceptable, or increase the production and generate considerable dark matter component out of sterile neutrino with otherwise negligible mixing with SM. In this paper we provide analytic estimates complementing and providing details of the numerical calculations performed in \cite{Bezrukov:2018wvd} in the case of resonant amplification of the sterile neutrino production. We also discuss phenomenological and theoretical issues related to the successful implementation of this idea in fully realistic extensions of the Standard Model of particle physics.

We study the implications of the Planck temperature power spectrum at low multipoles, $\ell<1000$, and SPTPol data. We show that this combination predicts consistent lensing-induced smoothing of acoustic peaks within $\Lambda$CDM cosmology and yields the robust predictions of the cosmological parameters. Combining only the Planck large-scale temperature data and the SPTPol polarization and lensing measurements within $\Lambda$CDM model we found substantially lower values of linear matter density perturbation $\sigma_8$ which bring the late-time parameter $S_8=\sigma_8\sqrt{\Omega_m/0.3}=0.763\pm0.022$ into accordance with galaxy clustering and weak lensing measurements. It also raises up the Hubble constant $H_0=69.68\pm1.00{\rm \,\,km\,s^{-1}Mpc^{-1}}$ that reduces the Hubble tension to the $2.5\sigma$ level. We examine the residual tension in the Early Dark Energy (EDE) model which produces the brief energy injection prior to recombination. We implement both the background and perturbation evolutions of the scalar field which potential scales as $V(\phi)\propto \phi^{2n}$. Including cosmic shear measurements (KiDS, VIKING-450, DES) and local distance-ladder data (SH0ES) to the combined fit we found that EDE completely alleviates the Hubble tension while not degradating the fit to large-scale structure data. The EDE scenario significantly improves the goodness-of-fit by $2.9\sigma$ in comparison with the concordance $\Lambda$CDM model. The account for the intermediate-redshift data (the supernova dataset and baryon acoustic oscillation data) fits perfectly to our parameter predictions and indicates the preference of EDE over $\Lambda$CDM at $3\sigma$.

We discuss the gravitational creation of superheavy particles $\chi$ in an inflationary scenario with a quartic potential and a non-minimal coupling between the inflaton $\varphi$ and the Ricci curvature: $\xi \varphi^2 R/2$. We show that for large constants $\xi >> 1$, there can be abundant production of particles $\chi$ with masses largely exceeding the inflationary Hubble rate $H_{infl}$, up to $(a~few) \times \xi H_{infl}$, even if they are conformally coupled to gravity. We discuss two scenarios involving these gravitationally produced particles $\chi$. In the first scenario, the inflaton has only gravitational interactions with the matter sector and the particles $\chi$ reheat the Universe. In this picture, the inflaton decays only due to the cosmic expansion, and effectively contributes to dark radiation, which can be of the observable size. The existing limits on dark radiation lead to an upper bound on the reheating temperature. In the second scenario, the particles $\chi$ constitute Dark Matter, if substantially stable. In this case, their typical masses should be in the ballpark of the Grand Unification scale.

We calculate the signal rate of hypothetical heavy neutral leptons (HNL or sterile neutrinos) from kaon decays expected in the framework of the SHiP experiment. The kaons are produced in the hadronic shower initiated in the beam-dump mode by 400\,GeV protons from CERN SPS. For a sufficiently light HNL (when the decays are kinematically allowed) we find kaon decays to be a noticeably richer source of HNL as compared to $D$-meson decays adopted in previous studies of the HNL phenomenology at SHiP. In particular, SHiP is capable of fully exploring the central part of the kinematically allowed region of the HNL mass and mixing with active neutrinos down to the lower cosmological bound. The latter is associated with HNL decays in the early Universe to energetic products rescattering off and thus destroying light nuclei produced at the primordial nucleosynthesis. A consistency of the HNL model with smaller mixing would require either a non-standard cosmology or new ingredients in the HNL sector, closing the room for the minimal seesaw type I model with sterile neutrinos lighter than kaons.

We consider Dark Matter composed of an oscillating singlet scalar field. On top of the mass term, the scalar is equipped with a potential spontaneously breaking Z_2-symmetry. This potential dominates at early times and leads to the time-dependent expectation value of the scalar, which decreases in the expanding Universe. As it drops below some critical value, the symmetry gets restored, and the Dark Matter field starts to oscillate around zero. We arrange the spontaneous symmetry breaking through the interaction of the scalar with the Ricci curvature. In that way, superheavy Dark Matter can be produced at very early times. Depending on its mass, the production takes place at inflation (very large masses up to the Grand Unification scale), at preheating, or at radiation-dominated stage (masses 10^{6}-10^{7} Gev).

Baksan Experiment on Sterile Neutrino (BEST) is presently at the stage of production of the artificial neutrino source Cr-51, the gallium exposure will start in July and proceed for three months. While aiming specifically at investigating the Gallium neutrino anomaly (SAGE and GALLEX experiments), BEST can do more and it is tempting to estimate its ability in testing sterile neutrino explanation of antineutrino (reactor) anomalies. We observe a moderate sensitivity to the region in model parameter space (sterile neutrino mass and mixing with active electron neutrino) outlined by the old reactor antineutrino anomaly and the best fit of DANSS experiment, while the Neutrino-4 favorite region falls right in the BEST ballpark. In particular, by analyzing SAGE+GALLEX and Neutrino-4 $\chi^2$ distributions we find that Neutrino-4 results are fully consistent with the Gallium anomaly, the significance of the combined anomaly almost reaches 4$\sigma$ level. If the BEST confirms the Neutrino-4 results, the joint analysis will indicate more than 5$\sigma$ evidence for the sterile neutrino of eV-scale mass.

We study in detail the recently proposed mechanism of generating superheavy Dark Matter with the mass larger than the Hubble rate at the end of inflation. A real scalar field constituting Dark Matter linearly couples to the inflaton. As a result of this interaction, the scalar gets displaced from its zero expectation value. This offset feeds into the energy density of Dark Matter. This mechanism is universal and can be implemented in a generic inflationary scenario. Phenomenology of the model is comprised of Dark Matter decay into inflatons, which in turn decay into Standard Model species triggering cascades of high energy particles contributing to the cosmic ray flux. We evaluate the lifetime of Dark Matter and obtain limits on the inflationary scenarios, where this mechanism does not lead to the conflict with the Dark Matter stability considerations/studies of cosmic ray propagation.

We describe two new generation mechanisms for Dark Matter composed of sterile neutrinos with ${\cal O}(1)$ keV mass. The model contains a light scalar field which coherently oscillates in the early Universe and modulates the Majorana mass of the sterile neutrino. In a region of model parameter space, the oscillations between active and sterile neutrinos are resonantly enhanced. This mechanism allows us to produce sterile neutrino DM with small mixing angle with active neutrinos, thus evading the X-ray constraints. At the same time the spectrum of produced DM is much cooler, than in the case of ordinary oscillations in plasma, opening a window of lower mass DM, which is otherwise forbidden by structure formation considerations. In other regions of the model parameter space, where the resonance does not appear, another mechanism can operate: large field suppresses the active-sterile oscillations, but instead sterile neutrinos are produced by the oscillating scalar field when the effective fermion mass crosses zero. In this case DM component is cold, and even 1 keV neutrino is consistent with the cosmic structure formation.

We introduce a class of non-supersymmetric models explaining baryogenesis a la Affleck-Dine, which use a decay of two superheavy scalar fields with close masses. These scalars acquire non-zero expectation values during inflation through linear couplings to a function of an inflaton. After the inflaton decay, the model possesses approximate U(1)-invariance, explicitly broken by a small mass splitting. This splitting leads to the baryogenesis in the early Universe. Resulting baryon asymmetry is automatically small for the scalars with the masses about the Grand Unification scale and larger. It is fully determined by the inflaton dynamics and the Lagrangian parameters, i.e., is independent of initial pre-inflationary conditions for the scalars. As a consequence, baryon perturbations are purely adiabatic. We point out a possible origin of the mass splitting: masses of scalars degenerate at some large energy scale may acquire different loop corrections due to the interaction with the inflaton. Compared to electroweak baryogenesis and conventional Affleck-Dine scenarios, our mechanism generically leads to the proton decay suppressed by the powers of the superheavy scalar masses, which makes this scenario potentially testable.

We study a model of a keV-scale sterile neutrino with a relatively large mixing with the Standard Model sector. Usual considerations predict active generation of such particles in the early Universe, which leads to constraints from the total Dark Matter density and absence of X-ray signal from sterile neutrino decay. These bounds together may deem any attempt of creation of the keV scale sterile neutrino in the laboratory unfeasible. We argue that for models with a hidden sector coupled to the sterile neutrino these bounds can be evaded, opening new perspectives for the direct studies at neutrino experiments such as Troitsk $\nu$-mass and KATRIN. We estimate the generation of sterile neutrinos in scenarios with the hidden sector dynamics keeping the sterile neutrinos either massless or superheavy in the early Universe. In both cases the generation by oscillations from active neutrinos in plasma is suppressed.

We propose a universal mechanism of producing dark matter and baryon (lepton) charge at the stage of the quasi-de Sitter expansion of the Universe---inflation. The key ingredient of the mechanism is a linear coupling of the field, responsible for generation of dark matter or baryon (lepton) charge, to a function of the inflaton. During inflation this induces almost constant force dragging the corresponding field to the non-zero value. This force explicitly breaks quantum numbers associated with dark matter/baryon abundance at later stages. As a particular realization of the mechanism we introduce a super-heavy complex scalar field with the mass larger than the Hubble rate during the last e-folds of inflation. The global U(1)-symmetry is violated due to the linear coupling of the phase of the complex scalar to the inflaton. The symmetry breaking leads to the generation of a non-zero Noether charge. The latter is directly related to the dark matter abundance, or, alternatively, can be converted into baryon asymmetry, if the complex scalar carries the baryon charge.

It has been recently suggested~\cite{Berezhiani:2015yta} that a subdominant fraction of dark matter decaying after recombination may alleviate tension between high-redshift (CMB anisotropy) and low-redshift (Hubble constant, cluster counts) measurements. In this report, we continue our previous study~\cite{Chudaykin:2016yfk} of the decaying dark matter (DDM) model adding all available recent baryon acoustic oscillation (BAO) and redshift space distortions (RSD) measurements. We find, that the BAO/RSD measurements generically prefer the standard $\Lambda$CDM and combined with other cosmological measurements impose an upper limit on the DDM fraction at the level of $\sim$\,5\,\%, strengthening by a factor of 1.5 limits obtained in \cite{Chudaykin:2016yfk} mostly from CMB data. However, the numbers vary from one analysis to other based on the same Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 12 (DR12) galaxy sample. Overall, the model with a few percent DDM fraction provides a better fit to the combined cosmological data as compared to the $\Lambda$CDM: the cluster counting and direct measurements of the Hubble parameter are responsible for that. The improvement can be as large as 1.5\,$\sigma$ and grows to 3.3\,$\sigma$ when the CMB lensing power amplitude ${\rm A_L}$ is introduced as a free fitting parameter.

The Starobinsky model of modified gravity suggested to explain dark energy may be also considered in the astrophysical context. Recently it has been pointed out that in contracting regions curvature oscillations around the GR value may lead to the production of high energy particles which contribute to the cosmic ray flux. We revisit these calculations in the Einstein frame and show that the continuous approximation for the matter density used in the original calculations is not valid. We show that this problem is generic in $F(R)$-gravity models introduced to describe the dark energy. We go beyond the approximation and find the rate of particle production to be negligible.

Small evaporating black holes were proposed to be dangerous inducing fast decay of the electroweak false vacuum. We observe that the flat-spectrum matter perturbations growing at the post-inflationary matter dominated stage can produce such black holes in a tiny amount which may nevertheless be sufficient to destroy the vacuum in the visible part of the Universe via the induced process. If the decay probability in the vicinity of Planck-mass black holes was of order one as suggested in literature, the absence of such objects in the early Universe would put severe constraints on inflation and subsequent stages thus excluding many well-motivated models (e.g. the R^2-inflation) and supporting the need of new physics in the Higgs sector. We give a qualitative argument, however, that exponential suppression of the probability should persist in the limit of small black hole masses. This suppression relaxes our cosmological constraints, and, if sufficiently strong, may cancel them.

String scenarios generically predict that we live in a so called axiverse: the Universe with about a hundred of light axion species which are decoupled from the Standard Model particles. However, the axions can couple to the inflaton which leads to their production after inflation. Then, these axions remain in the expanding Universe contributing to the dark radiation component, which is severely bounded from present cosmological data. We place a general constraint on the axion production rate and apply it to several variants of reasonable inflaton-to-axion couplings. The limit merely constrains the number of ultralight axions and the relative strength of inflaton-to-axion coupling. It is valid in both large and small field inflationary models irrespectively of the axion energy scales and masses. Thus, the limit is complementary to those associated with the Universe overclosure and axion isocurvature fluctuations. In particular, a hundred of axions is forbidden if inflaton universally couples to all the fields at reheating. In the case of gravitational sector being responsible for the reheating of the Universe (which is a natural option in all inflationary models with modified gravity), the axion production can be efficient. We find that in the Starobinsky $R^2$-inflation even a single axion (e.g. the standard QCD-axion) is in tension with the Planck data, making the model inconsistent with the axiverse. The general conclusion is that an inflation with inefficient reheating mechanism and low reheating temperature may be in tension with the presence of light scalars.

We present a comprehensive review of keV-scale sterile neutrino Dark Matter, collecting views and insights from all disciplines involved - cosmology, astrophysics, nuclear, and particle physics - in each case viewed from both theoretical and experimental/observational perspectives. After reviewing the role of active neutrinos in particle physics, astrophysics, and cosmology, we focus on sterile neutrinos in the context of the Dark Matter puzzle. Here, we first review the physics motivation for sterile neutrino Dark Matter, based on challenges and tensions in purely cold Dark Matter scenarios. We then round out the discussion by critically summarizing all known constraints on sterile neutrino Dark Matter arising from astrophysical observations, laboratory experiments, and theoretical considerations. In this context, we provide a balanced discourse on the possibly positive signal from X-ray observations. Another focus of the paper concerns the construction of particle physics models, aiming to explain how sterile neutrinos of keV-scale masses could arise in concrete settings beyond the Standard Model of elementary particle physics. The paper ends with an extensive review of current and future astrophysical and laboratory searches, highlighting new ideas and their experimental challenges, as well as future perspectives for the discovery of sterile neutrinos.

It has been recently suggested~\cite{Berezhiani:2015yta} that emerging tension between cosmological parameter values derived in high-redshift (CMB anisotropy) and low-redshift (cluster counts, Hubble constant) measurements can be reconciled in a model which contains subdominant fraction of dark matter decaying after recombination. We check the model against the CMB Planck data. We find that lensing of the CMB anisotropies by the large-scale structure gives strong extra constraints on this model, limiting the fraction as $F<8\%$ at 2\,$\sigma$ confidence level. However, investigating the combined data set of CMB and conflicting low-$z$ measurements, we obtain that the model with $F\approx2\!-\!5$\% exhibits better fit (by 1.5-3\,$\sigma$ depending on the lensing priors) compared to that of the concordance $\Lambda$CDM cosmological model.

In light of the Higgs boson discovery we reconsider generation of the baryon asymmetry in the non-minimal split Supersymmetry model with an additional singlet superfield in the Higgs sector. We find that successful baryogenesis during the first order electroweak phase transition is possible within phenomenologically viable part of the model parameter space. We discuss several phenomenological consequences of this scenario, namely, predictions for the electric dipole moments of electron and neutron and collider signatures of light charginos and neutralinos.

A classical evolution in chaotic inflationary models starts at high energy densities with semi-classical initial conditions presumably consistent with universal quantum nature of all the fundamental forces. That is each quantum contributes the same amount to the energy density. We point out the upper limit on this amount inherent in this approach, so that all the quanta are inside the weak-coupling domain. We discuss this issue in realistic models with modified gravity, $R^2$- and Higgs-inflations, emphasizing the specific change of the initial conditions with metric frame, while all the quanta still contribute equal parts. The analysis can be performed straightforwardly in any model with modified gravity ($F(R)$-gravity, scalars with non-minimal couplings to gravity, etc).

Numerous anomalous results in neutrino oscillation experiments can be attributed to interference of ~1 eV sterile neutrino. The specially designed to fully explore the Gallium anomaly Baksan Experiment on Sterile Transitions (BEST) starts next year. We investigate the sensitivity of BEST in searches for sterile neutrino mixed with electron neutrino. Then, performing the combined analysis of all the Gallium experiments (SAGE, GALLEX, BEST) we find the regions in model parameter space (sterile neutrino mass and mixing angle), which will be excluded if BEST agrees with no sterile neutrino hypothesis. For the opposite case, if BEST observes the signal as it follows from the sterile neutrino explanation of the Gallium (SAGE and GALLEX) anomaly, we show how BEST will improve upon the present estimates of the model parameters.

For the Warm Dark Matter (WDM) candidates the momentum distribution of particles becomes important, since it can be probed with observations of Lyman-$\alpha$ forest structures and confronted with coarse grained phase space density in galaxy clusters. We recall the calculation bt Kaplinghat (2005) of the spectrum in case of dark matter non-thermal production in decays of heavy particles emphasizing on the inherent applicability conditions, which are rather restrictive and sometimes ignored in literature. Cold part of the spectrum requires special care when WDM is considered.

We address the cosmological role of an additional ${\cal O}(1)$ eV sterile neutrino in modified gravity models. We confront the present cosmological data with predictions of the FLRW cosmological model based on a variant of $f(R)$ modified gravity proposed by one of the authors previously. This viable cosmological model which deviation from general relativity with a cosmological constant $\Lambda$ decreases as $R^{-2n}$ for large, but not too large values of the Ricci scalar $R$ provides an alternative explanation of present dark energy and the accelerated expansion of the Universe. Various up-to-date cosmological data sets exploited include Planck CMB anisotropy, CMB lensing potential, BAO, cluster mass function and Hubble constant measurements. We find that the CMB+BAO constraints strongly the sum of neutrino masses from above. This excludes values $\lambda\sim 1$ for which distinctive cosmological features of the model are mostly pronounced as compared to the $\Lambda$CDM model, since then free streaming damping of perturbations due to neutrino rest masses is not sufficient to compensate their extra growth occurring in $f(R)$ gravity. Thus, we obtain $\lambda>8.2$ ($2\sigma$) with cluster systematics and $\lambda>9.4$ ($2\sigma$) without that. In the latter case we find for the sterile neutrino mass $0.47\,\,\rm{eV}$$\,<\,$$m_{\nu,\,\rm{sterile}}$$\,<\,$$1\,\,\rm{eV}$ ($2\sigma$) assuming the active neutrinos are massless, not significantly larger than in the standard $\Lambda$CDM with the same data set: $0.45\,\,\rm{eV}$$\,<\,$$m_{\nu,\,\rm{sterile}}$$\,<\,$$0.92\,\,\rm{eV}$ ($2\sigma$). However, a possible discovery of a sterile neutrino with the mass $m_{\nu,\,\rm{sterile}} \approx 1.5\,$eV motivated by various anomalies in neutrino oscillation experiments would favor cosmology based on $f(R)$ gravity rather than the $\Lambda$CDM model.

We address the question of unlikeness of $R^2$- and Higgs inflations exhibiting exponentially flat potentials and hence apparently violating the inherent in a chaotic inflation initial condition when kinetic, gradient and potential terms are all of order one in Planck units. Placing the initial conditions in the Jourdan frame we find both models not worse than any other models with unbounded from above potentials: the terms in the Einstein frame are all of the same order, though appropriately smaller.

We study the mechanism of generation of 7 keV sterile neutrino Dark Matter (DM) in the model with light inflaton $\chi$, which serves as a messenger of scale invariance breaking. In this model the inflaton, in addition to providing reheating to the Standard Model (SM) particles, decays directly into sterile neutrinos. The latter are responsible for the active neutrino oscillations via seesaw type I like formula. While the two sterile neutrinos may also produce the lepton asymmetry in the primordial plasma and hence explain the baryon asymmetry of the Universe, the third one being the lightest may be of 7 keV and serve as DM. For this mechanism to work, the mass of the inflaton is bound to be light (0.1-1 GeV) and uniquely determines its properties, which allows to test the model. For particle physics experiments these are: inflaton lifetime ($10^{-5}-10^{-12}$ s), partial decay width of B-meson to kaon and inflaton ($10^{-6}-10^{-4}$) and inflaton branching ratios into light SM particles like it would be for the SM Higgs boson of the same mass. For cosmological experiments these are: spectral index of scalar perturbations ($n_s\simeq0.957-0.967$), and amount of tensor perturbations produced at inflation (tensor-to-scalar ratio $r\simeq0.15-0.005$).

We study the creation of high energy SM particles in the Starobinsky model of dark energy (a variant of F(R)-gravity) inside the regions contracting due to the Jeans instability. In this modification of gravity the additional degree of freedom -- scalaron -- behaves as a particle with mass depending on matter density. So when the mass changes light scalarons could be created at a non-adiabatic stage. Later scalaron mass grows and could reach large values, even the value of 10^{13} GeV, favored by early-time inflation. Heavy scalarons decay contributing to the cosmic ray flux. We analytically calculated the number density of created particles for the case of exponential (Jeans) contraction and found it negligibly small provided the phenomenologically viable and cosmologically interesting range of model parameters. We expect similar results for a generic model of F(R)-gravity mimicking cosmological constant.

Recent cosmological data favour R^2-inflation and some amount of non-standard dark radiation in the Universe. We show that a framework of high energy scale invariance can explain these data. The spontaneous breaking of this symmetry provides gravity with the Planck mass and particle physics with the electroweak scale. We found that the corresponding massless Nambu--Goldstone bosons -- dilatons -- are produced at reheating by the inflaton decay right at the amount needed to explain primordial abundances of light chemical elements and anisotropy of the cosmic microwave background. Then we extended the discussion on the interplay with Higgs-inflation and on general class of inflationary models where dilatons are allowed and may form the dark radiation. As a result we put a lower limit on the reheating temperature in a general scale invariant model of inflation.

We point out that recently discovered 3.5 keV line in X-ray spectra from various galaxy clusters and the Andromeda galaxy can be naturally explained by physics of a sector responsible for spontaneous supersymmetry breaking in models with R-parity. In this scenario the source of this line could be decay of sgoldstino - scalar superpartner of massive gravitino. At the same time the dominant dark matter component is stable gravitino whose mass is predicted to be about 0.25 GeV.

We introduce conformal coupling of the Standard Model Higgs field to gravity and discuss the subsequent modification of $R^2$-inflation. The main observation is a lower temperature of reheating which happens mostly through scalaron decays into gluons due to the conformal (trace) anomaly. This modifies all predictions of the original $R^2$-inflation. To the next-to-leading order in slow roll parameters we calculate amplitudes and indices of scalar and tensor perturbations produced at inflation. The results are compared to the next-to-leading order predictions of $R^2$-inflation with minimally coupled Higgs field and of Higgs-inflation. We discuss additional features in gravity wave signal that may help to distinguish the proposed variant of $R^2$-inflation. Remarkably, the features are expected in the region available for study at future experiments like BBO and DECIGO. Finally, we check that (meta)stability of electroweak vacuum in the cosmological model is consistent with recent results of searches for the Higgs boson at LHC.

We search for ultra-high energy photons by analyzing geometrical properties of shower fronts of events registered by the Telescope Array surface detector. By making use of an event-by-event statistical method, we derive upper limits on the absolute flux of primary photons with energies above 10^19, 10^19.5 and 10^20 eV based on the first three years of data taken.

We show that the cusp in the dark matter (DM) distribution required to explain the recently found excess in the gamma-ray spectrum at energies ~130 GeV in terms of the DM annihilations cannot survive the tidal forces if it is offset by ~1.5^\circ from the Galactic centre as suggested by observations.

We update the allowed parameter space of the simple chaotic inflationary model with quartic potential and light inflaton [Bezrukov,Gorbunov'2009] taking into account recent results from cosmology (CMB observations from SPT, ACT and WMAP) and from particle physics (LHC hints of the SM Higgs boson). The non-minimal (yet small) coupling to gravity of the inflaton becomes essential to fit the observational data. The inflaton has mass above 300 MeV and can be searched for at B-factories in B-meson two-body decays to kaon and inflaton. The inflaton lifetime depends on the model parameters, resulting in various inflaton signatures: either a missing energy, or a displaced vertex from the B-meson decay position, or a resonance in the Dalitz plot of a three particle decay. We also discuss the implementation of the inflaton model to the nuMSM, where the inflaton can be responsible for production of the dark matter sterile neutrino in the early Universe.

Standard model of particle physics fails to explain neutrino oscillations, dark matter and baryon asymmetry of the Universe. All these problems can be solved with three sterile neutrinos added to the model. Quite remarkably, if sterile neutrino masses are well below the electroweak scale, this modification --- Neutrino Minimal Standard Model (nuMSM) --- can be tested experimentally. We discuss a new experiment on search for decays of GeV-scale sterile neutrinos, which in the framework of the nuMSM are responsible for the matter-antimatter asymmetry generation and for the active neutrino masses. If lighter than 2 GeV, these particles can be produced in decays of charm mesons generated by high energy protons in a target of a proton beam-dump experiment, and subsequently decay into light SM particles. In order to fully explore this sector of the nuMSM, the new experiment requires data obtained with at least 10^{20} incident protons on target (achievable at CERN SPS in the future), and a big volume detector constructed from a large amount of identical single modules, with a total sterile neutrino decay length of few kilometers. The preliminary feasibility study for the proposed experiment shows that it has sensitivity which may either lead to the discovery of new particles below the Fermi scale --- right-handed partners of neutrinos --- or rule out see-saw sterile neutrinos with masses below 2 GeV.

Gravity takes care of both inflation and subsequent reheating in Starobinsky's R^2-model. The latter is due to inflaton gravitation decays dominated by scalar particle production. It is tempting to suggest that dark matter particles are also produced in this process. Since free scalars being too hot cannot serve as viable dark matter (Phys.Lett.B700:157-162,2011), we further study the issue considering two options: scalars with non-minimal coupling to gravity and superheavy scalars generated at inflationary stage. We found that the first option allows for viable warm or cold dark matter if scalar mass exceeds 0.8 MeV. The second option implies supercold dark matter with particle mass 10^16 GeV, which production is saturated at the end of inflation when inflaton-dependent scalar mass rapidly changes and violates adiabaticity. Similar result holds for superheavy fermion dark matter.

We present three features which can be used to distinguish the R^2-inflation Higgs-inflation from with ongoing, upcoming and planned experiments, assuming no new physics (apart form sterile neutrinos) up to inflationary scale. (i) Slightly different tilt of the scalar perturbation spectrum n_s and ratio r of scalar-to-tensor perturbation amplitudes. (ii) Gravity waves produced within R^2-model by collapsing, merging and evaporating scalaron clumps formed in the post-inflationary Universe. (iii) Different ranges of the possible Standard Model Higgs boson masses, where the electroweak vacuum remains stable while the Universe evolves after inflation. Specifically, in the R^2-model Higgs boson can be as light as 116 GeV. These effects mainly rely on the lower reheating temperature in the R^2-inflation.

The Telescope Array (TA) experiment, located in the western desert of Utah,USA, is designed for observation of extensive air showers from extremely high energy cosmic rays. The experiment has a surface detector array surrounded by three fluorescence detectors to enable simultaneous detection of shower particles at ground level and fluorescence photons along the shower track. The TA surface detectors and fluorescence detectors started full hybrid observation in March, 2008. In this article we describe the design and technical features of the TA surface detector.

Superluminal neutrinos are expected to lose energy due to bremstrauhlung. It is dominated by e+e--pair production if kinematically allowed. The same signature was used in searches for 3-body decays of hypothetical heavy sterile neutrinos. From the absence of these processes in CERN PS191 and CHARM experiments we set upper limits on the neutrino velocity in the energy range from 0.2 GeV to 280 GeV. Our limits are well below the neutrino velocity favored by the recent OPERA results. For energy-independent neutrino velocity the limits obtained in this paper are stronger than those coming from ICARUS experiment and observations of Supernova SN1987a.

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

The Telescope Array's Middle Drum fluorescence detector was instrumented with telescopes refurbished from the High Resolution Fly's Eye's HiRes-1 site. The data observed by Middle Drum in monocular mode was analyzed via the HiRes-1 profile-constrained geometry reconstruction technique and utilized the same calibration techniques enabling a direct comparison of the energy spectra and energy scales between the two experiments. The spectrum measured using the Middle Drum telescopes is based on a three-year exposure collected between December 16, 2007 and December 16, 2010. The calculated difference between the spectrum of the Middle Drum observations and the published spectrum obtained by the data collected by the HiRes-1 site allows the HiRes-1 energy scale to be transferred to Middle Drum. The HiRes energy scale is applied to the entire Telescope Array by making a comparison between Middle Drum monocular events and hybrid events that triggered both Middle Drum and the Telescope Array's scintillator Ground Array.

Since 2007, the Telescope Array (TA) experiment, based in Utah, USA, has been observing ultra high energy cosmic rays to understand their origins. The experiment involves a surface detector (SD) array and three fluorescence detector (FD) stations. FD stations, installed surrounding the SD array, measure the air fluorescence light emitted from extensive air showers (EASs) for precise determination of their energies and species. The detectors employed at one of the three FD stations were relocated from the High Resolution Fly's Eye experiment. At the other two stations, newly designed detectors were constructed for the TA experiment. An FD consists of a primary mirror and a camera equipped with photomultiplier tubes. To obtain the EAS parameters with high accuracies, understanding the FD optical characteristics is important. In this paper, we report the characteristics and installation of new FDs and the performances of the FD components. The results of the monitored mirror reflectance during the observation time are also described in this report.

The analysis of theories with non-minimal coupling of Higgs field to gravity revealed that they enter into strong coupling regime above certain Higgs-dependent cutoff, which may be considerably below the Planck scale. Assuming that the effective theory, complementing the Standard Model or its minimal extension--the nuMSM--contains a set of higher dimensional operators suppressed by the Higgs-dependent cutoff, we analyse the reheating of the Universe after the Higgs inflation. We show that extra terms do not spoil the Higgs inflation, but can lead to baryogenesis and to warm sterile neutrino dark matter production at the reheating stage of the Universe expansion. They can also result in neutrino mass generation and proton decay.

In R^2-inflation scalaron slow roll is responsible for the inflationary stage, while its oscillations reheat the Universe. We find that the same scalaron decays induced by gravity can also provide the dark matter production and leptogenesis. With R^2-term and three Majorana fermions added to the Standard Model, we arrive at the phenomenologically complete theory capable of simultaneously explaining neutrino oscillations, inflation, reheating, dark matter and baryon asymmetry of the Universe. Besides the seesaw mechanism in neutrino sector, we use only gravity, which solves all the problems by exploiting scalaron.

Comparing the signals measured by the surface and underground scintillator detectors of the Yakutsk Extensive Air Shower Array, we place upper limits on the integral flux and the fraction of primary cosmic-ray photons with energies E > 10^18 eV, E > 2*10^18 eV and E > 4*10^18 eV. The large collected statistics of the showers measured by large-area muon detectors provides a sensitivity to photon fractions < 10^-2, thus achieving precision previously unreachable at ultra-high energies.

We study the phenomenology of a realistic version of the chaotic inflationary model, which can be fully and directly explored in particle physics experiments. The inflaton mixes with the Standard Model Higgs boson via the scalar potential, and no additional scales above the electroweak scale are present in the model. The inflaton-to-Higgs coupling is responsible for both reheating in the Early Universe and the inflaton production in particle collisions. We find the allowed range of the light inflaton mass, 270 MeV<~m_chi<~1.8 GeV, and discuss the ways to find the inflaton. The most promising are two-body kaon and B-meson decays with branching ratios of orders 10^{-9} and 10^{-6}, respectively. The inflaton is unstable with the lifetime 10^{-9}--10^{-10} s. The inflaton decays can be searched for in a beam-target experiment, where, depending on the inflaton mass, from several billions to several tenths of millions inflatons can be produced per year with modern high-intensity beams.

We make use of the phase space density approach to discuss gravitino as a warm dark matter candidate. Barring fine tuning between the reheat temperature in the Universe and superparticle masses, we find that \emph{warm} gravitinos have both appropriate total mass density, Omega_{\tilde G} = \Omega_{DM} \simeq 0.2, and suitable primordial phase space density at low momenta provided that their mass is in the range 1 keV \lesssim \mg \lesssim 15 keV, the reheat temperature in the Universe is low, T_R \lesssim 10 TeV, and masses of some of the superparticles are sufficiently small, M \lesssim 350 GeV. The latter property implies that the gravitino warm dark matter scenario will be either ruled out or supported by the LHC experiments.

We apply phase-space density considerations to obtain lower bounds on the mass of sterile neutrino as dark matter candidate. The bounds are different for non-resonant production, resonant production in the presence of lepton asymmetry and production in decays of heavier particles. In the former case our bound is comparable to, but independent of the Lyman-alpha bound, and together with X-ray upper limit it disfavors non-resonantly produced sterile neutrino dark matter. An interesting feature of the latter case is that warm dark matter may be composed of heavy particles.

We propose a cosmological model which exhibits the phenomenon of self-acceleration: the Universe is attracted to the phase of accelerated expansion at late times even in the absence of the cosmological constant. The self-acceleration is inevitable in the sense that it cannot be neutralized by any negative explicit cosmological constant. The model is formulated in the framework of brane-world theories with a warped extra dimension. The key ingredient of the model is the brane-bulk energy transfer which is carried by bulk vector fields with a sigma-model-like boundary condition on the brane. We explicitly find the 5-dimensional metric corresponding to the late-time de Sitter expansion on the brane; this metric describes an AdS_5 black hole with growing mass. The present value of the Hubble parameter implies the scale of new physics of order 1 TeV, where the proposed model has to be replaced by putative UV-completion. The mechanism leading to the self-acceleration has AdS/CFT interpretation as occurring due to specific dynamics of conformal matter interacting with external "electric" fields. The Universe expansion history predicted by the model is distinguishable from the standard LambdaCDM cosmology.

A natural interpretation of the correlation between nearby Active Galactic Nuclei (AGN) and the highest-energy cosmic rays observed recently by the Pierre Auger Collaboration is that the sources of the cosmic rays are either AGN or other objects with a similar spatial distribution (the ``AGN hypothesis''). We question this interpretation. We calculate the expected distribution of the arrival directions of cosmic rays under the AGN hypothesis and argue that it is not supported by the data, one of manifestations of the discrepancy being the deficit of events from the direction of the Virgo supercluster. We briefly discuss possible alternative explanations including the origin of a significant part of the observed events from Cen A.

We analyse a sample of 33 extensive air showers (EAS) with estimated primary energies above 2\cdot 10^{19} eV and high-quality muon data recorded by the Yakutsk EAS array. We compare, event-by-event, the observed muon density to that expected from CORSIKA simulations for primary protons and iron, using SIBYLL and EPOS hadronic interaction models. The study suggests the presence of two distinct hadronic components, ``light'' and ``heavy''. Simulations with EPOS are in a good agreement with the expected composition in which the light component corresponds to protons and the heavy component to iron-like nuclei. With SYBILL, simulated muon densities for iron primaries are a factor of \sim 1.5 less than those observed for the heavy component, for the same electromagnetic signal. Assuming two-component proton-iron composition and the EPOS model, the fraction of protons with energies E>10^{19} eV is 0.52^{+0.19}_{-0.20} at 95% confidence level.

We argue that the data published by the Pierre Auger Collaboration (arXiv:0711.2256) disfavor at 99% confidence level their hypothesis that most of the highest-energy cosmic rays are protons from nearby astrophysical sources, either Active Galactic Nuclei or other objects with a similar spatial distribution.

The most common way to simplify extensive Monte-Carlo simulations of air showers is to use the thinning approximation. We study its effect on the physical parameters reconstructed from simulated showers. To this end, we have created a library of showers simulated without thinning with energies from 10^17 eV to 10^18 eV, various zenith angles and primaries. This library is publicly available. Physically interesting applications of the showers simulated without thinning are discussed. Observables reconstructed from these showers are compared to those obtained with the thinning approximation. The amount of artificial fluctuations introduced by thinning is estimated. A simple method, multisampling, is suggested which results in a controllable suppression of artificial fluctuations and at the same time requires less demanding computational resources as compared to the usual thinning.

We suggest a method which improves the precision of studies of the primary composition of ultra-high-energy cosmic rays. Two principal ingredients of the method are (1) comparison of the observed and simulated parameters for individual showers, without averaging over arrival directions and (2) event-by-event selection of simulated showers by the physical observables and not by the reconstructed primary parameters. A detailed description of the algorithm is presented and illustrated by several examples.

By making use of the data on the total signal and on the muon component of the air showers detected by the Yakutsk array, we analyze, in the frameworks of the recently suggested event-by-event approach, how large the fraction of primary gamma-rays at ultra-high energies can be. We derive upper limits on the photon fraction in the integral flux of primary cosmic rays. At the 95% confidence level (CL), these limits are 22% for primary energies E_0>4\cdot 10^{19}eV and 12% for E_0>2\cdot 10^{19}eV. Despite the presence of muonless events, the data are consistent with complete absence of photons at least at 95% CL. The sensitivity of the results to systematic uncertainties, in particular to those of the energy determination for non-photon primaries, is discussed.

We present the interpretation of the muon and scintillation signals of ultra-high-energy air showers observed by AGASA and Yakutsk extensive air shower array experiments. We consider case-by-case ten highest energy events with known muon content and conclude that at the 95% confidence level (C.L.) none of them was induced by a primary photon. Taking into account statistical fluctuations and differences in the energy estimation of proton and photon primaries, we derive an upper limit of 36% at 95% C.L. on the fraction of primary photons in the cosmic-ray flux above 10^20 eV. This result disfavors the Z-burst and superheavy dark-matter solutions to the GZK-cutoff problem.

It is shown by an explicit calculation that the excitations about the self-accelerating cosmological solution of the Dvali--Gabadaze--Porrati model contain a ghost mode. This raises serious doubts about viability of this solution. Our analysis reveals the similarity between the quadratic theory for the perturbations around the self-accelerating Universe and an Abelian gauge model with two Stueckelberg fields.

The existing correlation between BL Lacertae objects (BL Lacs) and cosmic-ray events observed by HiRes experiment provide sufficient information to formulate quantitatively the hypothesis about the flux of neutral cosmic-ray particles originated from BL Lacs. We determine the potential of future cosmic ray experiments to test this hypothesis by predicting the number of coincidences between arrival directions of cosmic rays and positions of BL Lacs on the celestial sphere, which should be observed in the future datasets. We find that the early Pierre Auger data will not have enough events to address this question. On the contrary, the final Pierre Auger data and the early Telescope Array data will be sufficient to fully test this hypothesis. If confirmed, it would imply the existence of highest-energy neutral particles coming from cosmological distances.

We found significant correlations between the arrival directions of the highest energy photons (E>10 GeV) observed by EGRET and positions of the BL Lac type objects (BL Lacs). The observed correlations imply that not less than three per cent of extragalactic photons at these energies originate from BL Lacs. Some of the correlating BL Lacs have no counterparts in the EGRET source catalog, i.e. do not coincide with strong emitters of gamma-rays at lower energy. The study of correlating BL Lacs suggests that they may form a subset which is statistically different from the total BL Lac catalog; we argue that they are prominent candidates for TeV gamma-ray sources. Our results demonstrate that the analysis of positional correlations is a powerful approach indispensable in cases when low statistics limits or even prohibits the standard case-by-case identification.

We consider various classes of persistent extragalactic astrophysical sources which have been suggested in literature as possible emitters of ultra-high-energy cosmic rays (UHECR). We compare the strength of the claimed correlations by a uniform procedure for all classes of sources by making use of the AGASA, Yakutsk and HiRes stereo data. BL Lac type objects correlate with the cosmic rays detected by all three independent experiments and are more probably, compared to other astrophysical sources, related to the UHECR origin. With the account of the Galactic magnetic field (not possible for the HiRes data at the moment), apart of BL Lac type objects, unidentified gamma-ray sources may be correlated with AGASA and Yakutsk cosmic rays.

Previously suggested correlations of BL Lac type objects with the arrival directions of the ultra-high-energy cosmic ray primaries are tested by making use of the HiRes stereoscopic data. The results of the study support the conclusion that BL Lacs may be the cosmic ray sources and suggest the presence of a small (a few percent) fraction of neutral primaries at E>10^{19} eV.

We study the large-scale distribution of the arrival directions of the highest energy cosmic rays observed by various experiments. Despite clearly insufficient statistics, we find a deficit of cosmic rays at energies higher than 10^{20} eV from a large part of the sky around the celestial North Pole. We speculate on possible explanations of this feature.

We calculate the cosmic microwave background (CMB) anisotropy spectrum in models with millicharged particles of electric charge q\sim 10^{-6}-10^{-1} in units of electron charge. We find that a large region of the parameter space for the millicharged particles exists where their effect on the CMB spectrum is similar to the effect of baryons. Using WMAP data on the CMB anisotropy and assuming Big Bang nucleosynthesis value for the baryon abundance we find that only a small fraction of cold dark matter, Omega_{mcp}h_0^2 < 0.007 (at 95% CL), may consists of millicharged particles with the parameters (charge and mass) from this region. This bound significantly narrows the allowed range of the parameters of millicharged particles. In models without paraphoton millicharged particles are now excluded as a dark matter candidate. We also speculate that recent observation of 511 keV gamma-rays from the Galactic bulge may be an indication that a (small) fraction of CDM is comprised of the millicharged particles.

Secondary photons from decays of metastable neutralinos can contribute to the ultra-high energy cosmic ray flux. The neutralino production rate is too low in acceleration mechanisms to affect the cosmic ray spectrum without emitting enormous energy in photons and neutrinos. However, in top-down models with sources not concentrated in galactic halos, neutralino decays change the spectrum significantly. We estimate the parameters of a model in which photons from neutralino decays are responsible for cosmic ray events with energies above 10^20 eV, and figure out distinctive experimental signatures for this model.

The Z-burst mechanism invoked to explain ultra-high energy cosmic rays is severely constrained by measurements of the cosmic gamma-ray background by EGRET. We discuss the case of optically thick sources and show that jets and hot spots of active galaxies cannot provide the optical depth required to suppress the photon flux. Other extragalactic accelerators (AGN cores and sites of gamma ray bursts), if they are optically thick, could be tested by future measurements of the secondary neutrino flux.

A set of potentially gamma-ray--loud BL Lac objects is selected by intersecting the EGRET and BL Lac catalogs. Of the resulting 14 objects, eight are found to correlate with arrival directions of ultra--high-energy cosmic rays (UHECRs), with significance of the order of 5 sigma. This suggests that gamma-ray emission can be used as a distinctive feature of those BL Lac objects that are capable of producing UHECR.

If Ultra High Energy Cosmic Rays (UHECRs) with E>4 10^{19} eV originate from BL Lacertae at cosmological distances as suggested by recent studies, the absence of the GZK cutoff can not be reconciled with Standard-Model particle properties. Axions would escape the GZK cutoff, but even the coherent conversion and back-conversion between photons and axions in large-scale magnetic fields is not enough to produce the required flux. However, one may construct models of other novel (pseudo)scalar neutral particles with properties that would allow for sufficient rates of particle production in the source and shower production in the atmosphere to explain the observations. As an explicit example for such particles we consider SUSY models with light sgoldstinos.

We consider a possibility to explain the observed suppression of the second acoustic peak in the anisotropy spectrum of the Cosmic Microwave Background (CMB) by interaction between a fraction of non-baryonic Cold Dark Matter (CDM) and normal baryonic matter. This scenario does not require any modifications in the standard Big Bang Nucleosynthesis (BBN). We estimate the required values of the cross-section-to-mass ratio for elastic scattering of CDM particles off baryons. In case of velocity-independent elastic scattering (in the velocity interval 10^{-5}--10^{-3}) we find that such particles do not contradict observational limits if they are heavier than \sim 10^5 GeV or lighter than \sim 0.5 GeV. Another candidate, which may appear in the models with infinite extra dimensions, is a quasistable charged particle decaying through tunneling into extra dimensions. Finally a millicharged particle with the electric charge ranging from \sim 10^{-4} to \sim 10^{-1} and with mass M\sim 0.1 GeV--1 TeV also may be responsible for the suppression of the second acoustic peak. As a byproduct we point out that CMB measurements set new limits on the allowed parameter space for the millicharged particles.