Table of contents

Warm inflation is, as of today, one of the best motivated mechanisms for explaining an early inflationary period. In this paper, we derive and analyze the current bounds on warm inflation with a monomial potential U ∝ ϕ^p, using the constraints from the PLANCK mission. In particular, we discuss the parameter space of the tensor-to-scalar ratio r and the potential coupling λ of the monomial warm inflation in terms of the number of e-folds. We obtain that the theoretical tensor-to-scalar ratio r ∼ 10⁻⁸ is much smaller than the current observational constrain r ≲ 0.12, despite a relatively large value of the field excursion Δ ϕ ∼ 0.1M_Pl. Warm inflation thus eludes the Lyth bound set on the tensor-to-scalar ratio by the field excursion.

We present a new N-body code, gevolution, for the evolution of large scale structure in the Universe. Our code is based on a weak field expansion of General Relativity and calculates all six metric degrees of freedom in Poisson gauge. N-body particles are evolved by solving the geodesic equation which we write in terms of a canonical momentum such that it remains valid also for relativistic particles. We validate the code by considering the Schwarzschild solution and, in the Newtonian limit, by comparing with the Newtonian N-body codes Gadget-2 and RAMSES. We then proceed with a simulation of large scale structure in a Universe with massive neutrinos where we study the gravitational slip induced by the neutrino shear stress. The code can be extended to include different kinds of dark energy or modified gravity models and going beyond the usually adopted quasi-static approximation. Our code is publicly available.

We present a new analytic approach to describe large scale structure formation in the mildly non-linear regime. The central object of the method is the time-dependent probability distribution function generating correlators of the cosmological observables at a given moment of time. Expanding the distribution function around the Gaussian weight we formulate a perturbative technique to calculate non-linear corrections to cosmological correlators, similar to the diagrammatic expansion in a three-dimensional Euclidean quantum field theory, with time playing the role of an external parameter. For the physically relevant case of cold dark matter in an Einstein-de Sitter universe, the time evolution of the distribution function can be found exactly and is encapsulated by a time-dependent coupling constant controlling the perturbative expansion. We show that all building blocks of the expansion are free from spurious infrared enhanced contributions that plague the standard cosmological perturbation theory. This paves the way towards the systematic resummation of infrared effects in large scale structure formation. We also argue that the approach proposed here provides a natural framework to account for the influence of short-scale dynamics on larger scales along the lines of effective field theory.

We study the decoupling process of neutrinos in the early universe in the presence of three-flavour oscillations. The evolution of the neutrino spectra is found by solving the corresponding momentum-dependent kinetic equations for the neutrino density matrix, including for the first time the proper collision integrals for both diagonal and off-diagonal elements. This improved calculation modifies the evolution of the off-diagonal elements of the neutrino density matrix and changes the deviation from equilibrium of the frozen neutrino spectra. However, it does not vary the contribution of neutrinos to the cosmological energy density in the form of radiation, usually expressed in terms of the effective number of neutrinos, N_eff. We find a value of N_eff = 3.045, in agreement with previous theoretical calculations and consistent with the latest analysis of Planck data. This result does not depend on the ordering of neutrino masses. We also consider the effect of non-standard neutrino-electron interactions (NSI), predicted in many theoretical models where neutrinos acquire mass. For two sets of NSI parameters allowed by present data, we find that N_eff can be reduced down to 3.040 or enhanced up to 3.059.

In this paper, we consider the Universe at the late stage of its evolution and deep inside the cell of uniformity. At these scales, we can consider the Universe to be filled with dust-like matter in the form of discretely distributed galaxies, a K-essence scalar field, playing the role of dark energy, and radiation as matter sources. We investigate such a Universe in the mechanical approach. This means that the peculiar velocities of the inhomogeneities (in the form of galaxies) as well as the fluctuations of the other perfect fluids are non-relativistic. Such fluids are designated as coupled because they are concentrated around the inhomogeneities. In the present paper, we investigate the conditions under which the K-essence scalar field with the most general form for its action can become coupled. We investigate at the background level three particular examples of the K-essence models: (i) the pure kinetic K-essence field, (ii) a K-essence with a constant speed of sound and (iii) the K-essence model with the Lagrangian bX+cX²−V(ϕ). We demonstrate that if the K-essence is coupled, all these K-essence models take the form of multicomponent perfect fluids where one of the component is the cosmological constant. Therefore, they can provide the late-time cosmic acceleration and be simultaneously compatible with the mechanical approach.

Recently, a family of hairy black holes in 4-dimensional Einstein gravity minimally coupled to a complex, massive scalar field was discovered [1]. Besides the mass M and spin angular momentum J, these objects are characterized by a Noether charge Q, measuring the amount of scalar hair, which is not associated to a Gauss law and cannot be measured at spatial infinity. Introducing a dimensionless scalar hair parameter q, ranging from 0 to 1, we recover (a subset of) Kerr black holes for q = 0 and a family of rotating boson stars for q = 1. In the present paper, we explore the possibility of measuring q for astrophysical black holes with current and future X-ray missions. We study the iron Kα line expected in the reflection spectrum of such hairy black holes and we simulate observations with Suzaku and eXTP. As a proof of concept, we point out, by analyzing a sample of hairy black holes, that current observations can already constrain the scalar hair parameter q, because black holes with q close to 1 would have iron lines definitively different from those we observe in the available data. We conclude that a detailed scanning of the full space of solutions, together with data from the future X-ray missions, like eXTP, will be able to put relevant constraints on the astrophysical realization of Kerr black holes with scalar hair.

As we are entering the era of precision cosmology, it is necessary to count on accurate cosmological predictions from any proposed model of dark matter. In this paper we present a novel approach to the cosmological evolution of scalar fields that eases their analytic and numerical analysis at the background and at the linear order of perturbations. The new method makes use of appropriate angular variables that simplify the writing of the equations of motion, and which also show that the usual field variables play a secondary role in the cosmological dynamics. We apply the method to a scalar field endowed with a quadratic potential and revisit its properties as dark matter. Some of the results known in the literature are recovered, and a better understanding of the physical properties of the model is provided. It is confirmed that there exists a Jeans wavenumber k_J, directly related to the suppression of linear perturbations at wavenumbers k > k_J, and which is verified to be k_J = a √mH. We also discuss some semi-analytical results that are well satisfied by the full numerical solutions obtained from an amended version of the CMB code CLASS. Finally we draw some of the implications that this new treatment of the equations of motion may have in the prediction of cosmological observables from scalar field dark matter models.

We investigate the dependency of Higgs inflation on the non-renormalisable matching between the low energy Standard Model limit and the inflationary regime at high energies. We show that for the top mass range m_t ≳ 171.8 GeV the scenario robustly predicts the spectral index n_s ≃ 0.97 and the tensor-to-scalar ratio r ≃ 0.003. The matching is however non-trivial, even the best-fit values m_h = 125.09 GeV and m_t = 173.21 GeV require a jump δλ ∼ 0.01 in the Higgs coupling below the inflationary scale. For m_t ≲ 171.8 GeV, the matching may generate a feature in the inflationary potential. In this case the predicted values of n_s and r vary but the model is still falsifiable. For example, a detection of negative running of spectral index at level α_s ≲ −0.01 would rule out Higgs inflation.

We present a new approximation to include fully general relativistic pressure and velocity in Newtonian hydrodynamics. The energy conservation, momentum conservation and two Poisson's equations are consistently derived from Einstein's gravity in the zero-shear gauge assuming weak gravity and action-at-a-distance limit. The equations show proper special relativity limit in the absence of gravity. Our approximation is complementary to the post-Newtonian approximation and the equations are valid in fully nonlinear situations.

We demonstrate how the image analysis technique of wavelet decomposition can be applied to the gamma-ray sky to separate emission on different angular scales. New structures on scales that differ from the scales of the conventional astrophysical foreground and background uncertainties can be robustly extracted, allowing a model-independent characterization with no presumption of exact signal morphology. As a test case, we generate mock gamma-ray data to demonstrate our ability to extract extended signals without assuming a fixed spatial template. For some point source luminosity functions, our technique also allows us to differentiate a diffuse signal in gamma-rays from dark matter annihilation and extended gamma-ray point source populations in a data-driven way.

We use the 2015 Planck likelihood in combination with the Bicep2/Keck likelihood (BKP and BK14) to constrain the chirality, χ, of primordial gravitational waves in a scale-invariant scenario. In this framework, the parameter χ enters theory always coupled to the tensor-to-scalar ratio, r, e.g. in combination of the form χ ⋅ r. Thus, the capability to detect χ critically depends on the value of r. We find that with present data sets χ is de facto unconstrained. We also provide forecasts for χ from future CMB experiments, including COrE+, exploring several fiducial values of r. We find that the current limit on r is tight enough to disfavor a neat detection of χ. For example, in the unlikely case in which r∼0.1(0.05), the maximal chirality case, i.e. χ = ±1, could be detected with a significance of ∼2.5(1.5)σ at best. We conclude that the two-point statistics at the basis of CMB likelihood functions is currently unable to constrain chirality and may only provide weak limits on χ in the most optimistic scenarios. Hence, it is crucial to investigate the use of other observables, e.g. provided by higher order statistics, to constrain these kinds of parity violating theories with the CMB.

We use classical lattice simulations in 3+1 dimensions to study the interplay between the resonant production of particles during preheating and the subsequent decay of these into a set of secondary species. We choose to work in a simplified version of Higgs inflation in which the Higgs field non-minimally coupled to gravity plays the role of the inflaton. Our numerical results extend the analytical estimates in the literature beyond the linear regime and shed some light on the limitations of the analytical techniques. The inclusion of fast and inefficient decays postpones the onset of parametric resonance by depleting the particles produced at the bottom of the potential. In spite of this delay, fermions are shown to play an important role on the destruction of the inflaton field. The limitations of our approach and its applications to a realistic Higgs inflation scenario are also discussed.

Motivated by the recent indications for a 750 GeV resonance in the di-photon final state at the LHC, in this work we analyse the compatibility of the excess with the broad photon excess detected at the Galactic Centre. Intriguingly, by analysing the parameter space of an effective models where a 750 GeV pseudoscalar particles mediates the interaction between the Standard Model and a scalar dark sector, we prove the compatibility of the two signals. We show, however, that the LHC mono-jet searches and the Fermi LAT measurements strongly limit the viable parameter space. We comment on the possible impact of cosmic antiproton flux measurement by the AMS-02 experiment.

We perform a detailed analysis of the synchrotron signals produced by dark matter annihilations and decays. We consider different set-ups for the propagation of electrons and positrons, the galactic magnetic field and dark matter properties. We then confront these signals with radio and microwave maps, including Planck measurements, from a frequency of 22 MHz up to 70 GHz. We derive two sets of constraints: conservative and progressive, the latter based on a modeling of the astrophysical emission. Radio and microwave constraints are complementary to those obtained with other indirect detection methods, especially for dark matter annihilating into leptonic channels.

We address the impact of consistent modifications of gravity on the largest observable scales, focusing on relativistic effects in galaxy number counts and the cross-correlation between the matter large scale structure (LSS) distribution and the cosmic microwave background (CMB). Our analysis applies to a very broad class of general scalar-tensor theories encoded in the Horndeski Lagrangian and is fully consistent on linear scales, retaining the full dynamics of the scalar field and not assuming quasi-static evolution. As particular examples we consider self-accelerating Covariant Galileons, Brans-Dicke theory and parameterizations based on the effective field theory of dark energy, using the hi class code to address the impact of these models on relativistic corrections to LSS observables. We find that especially effects which involve integrals along the line of sight (lensing convergence, time delay and the integrated Sachs-Wolfe effect—ISW) can be considerably modified, and even lead to O(1000%) deviations from General Relativity in the case of the ISW effect for Galileon models, for which standard probes such as the growth function only vary by O(10%). These effects become dominant when correlating galaxy number counts at different redshifts and can lead to ∼ 50% deviations in the total signal that might be observable by future LSS surveys. Because of their integrated nature, these deep-redshift cross-correlations are sensitive to modifications of gravity even when probing eras much before dark energy domination. We further isolate the ISW effect using the cross-correlation between LSS and CMB temperature anisotropies and use current data to further constrain Horndeski models. Forthcoming large-volume galaxy surveys using multiple-tracers will search for all these effects, opening a new window to probe gravity and cosmic acceleration at the largest scales available in our universe.

The neutrino mass hierarchy is one of the neutrino fundamental properties yet to be determined. We introduce a method to determine neutrino mass hierarchy by comparing the interaction rate of neutral current (NC) interactions, ν() + p → ν() + p, and inverse beta decays (IBD), _e + p → n + e⁺, of supernova neutrinos in scintillation detectors. Neutrino flavor conversions inside the supernova are sensitive to neutrino mass hierarchy. Due to Mikheyev-Smirnov-Wolfenstein effects, the full swapping of _e flux with the _x (x = μ, τ) one occurs in the inverted hierarchy, while such a swapping does not occur in the normal hierarchy. As a result, more high energy IBD events occur in the detector for the inverted hierarchy than the high energy IBD events in the normal hierarchy. By comparing IBD interaction rate with the mass hierarchy independent NC interaction rate, one can determine the neutrino mass hierarchy.

We obtain a closed system of equations for scalar perturbations in a multi-component braneworld. Our braneworld possesses a phantom-like equation of state at late times, w_eff < −1, but no big-rip future singularity. In addition to matter and radiation, the braneworld possesses a new effective degree of freedom—the `Weyl fluid' or `dark radiation'. Setting initial conditions on super-Hubble spatial scales at the epoch of radiation domination, we evolve perturbations of radiation, pressureless matter and the Weyl fluid until the present epoch. We observe a gradual decrease in the amplitude of the Weyl-fluid perturbations after Hubble-radius crossing, which results in a negligible effect of the Weyl fluid on the evolution of matter perturbations on spatial scales relevant for structure formation. Consequently, the quasi-static approximation of Koyama and Maartens provides a good fit to the exact results during the matter-dominated epoch. We find that the late-time growth of density perturbations on the brane proceeds at a faster rate than in ΛCDM. Additionally, the gravitational potentials Φ and Ψ evolve differently on the brane than in ΛCDM, for which Φ = Ψ. On the brane, by contrast, the ratio Φ/Ψ exceeds unity during the late matter-dominated epoch (z ≲ 50). These features emerge as smoking gun tests of phantom brane cosmology and allow predictions of this scenario to be tested against observations of galaxy clustering and large-scale structure.

We derive the contribution to relativistic galaxy number count fluctuations from vector and tensor perturbations within linear perturbation theory. Our result is consistent with the the relativistic corrections to number counts due to scalar perturbation, where the Bardeen potentials are replaced with line-of-sight projection of vector and tensor quantities. Since vector and tensor perturbations do not lead to density fluctuations the standard density term in the number counts is absent. We apply our results to vector perturbations which are induced from scalar perturbations at second order and give numerical estimates of their contributions to the power spectrum of relativistic galaxy number counts.

We employ an asteroseismic model of L19-2, a relatively massive (M_⋆ ∼ 0.75 M_⊙) and hot (T_eff ∼ 12 100 K) pulsating DA (H-rich atmosphere) white dwarf star (DAV or ZZ Ceti variable), and use the observed values of the temporal rates of period change of its dominant pulsation modes (Π ∼ 113 s and Π ∼ 192 s), to derive a new constraint on the mass of the axion, the hypothetical non-barionic particle considered as a possible component of the dark matter of the Universe. If the asteroseismic model employed is an accurate representation of L19-2, then our results indicate hints of extra cooling in this star, compatible with emission of axions of mass m_a cos²β ≲ 25 meV or an axion-electron coupling constant of g_ae ≲ 7 × 10⁻¹³.

The infrared dynamics of a light, minimally coupled scalar field in de Sitter spacetime with Ricci curvature R = 12H², averaged over horizon sized regions of physical volume V_H = (4π/3)(1/H)³, can be interpreted as Brownian motion in a medium with de Sitter temperature T_DS = ℏH/2π. We demonstrate this by directly deriving the effective action of scalar field fluctuations with wavelengths larger than the de Sitter curvature radius and generalizing Starobinsky's seminal results on stochastic inflation. The effective action describes stochastic dynamics and the fluctuating force drives the field to an equilibrium characterized by a thermal Gibbs distribution at temperature T_DS which corresponds to a de Sitter invariant state. Hence, approach towards this state can be interpreted as thermalization. We show that the stochastic kinetic energy of the coarse-grained description corresponds to the norm of ∂_μϕ and takes a well defined value per horizon volume ½⟨(∇ϕ)²⟩ = − ½T_DS/V_H. This approach allows for the non-perturbative computation of the de Sitter invariant stress energy tensor ⟨T_μν⟩ for an arbitrary scalar potential.

We present, for the first time in the literature, a full reconstruction of the total (linear and non-linear) ISW/Rees-Sciama effect in the presence of massive neutrinos, together with its cross-correlations with CMB-lensing and weak-lensing signals. The present analyses make use of all-sky maps extracted via ray-tracing across the gravitational potential distribution provided by the ``Dark Energy and Massive Neutrino Universe'' (DEMNUni) project, a set of large-volume, high-resolution cosmological N-body simulations, where neutrinos are treated as separate collisionless particles. We correctly recover, at 1–2% accuracy, the linear predictions from CAMB. Concerning the CMB-lensing and weak-lensing signals, we also recover, with similar accuracy, the signal predicted by Boltzmann codes, once non-linear neutrino corrections to HALOFIT are accounted for. Interestingly, in the ISW/Rees-Sciama signal, and its cross correlation with lensing, we find an excess of power with respect to the massless case, due to free streaming neutrinos, roughly at the transition scale between the linear and non-linear regimes. The excess is ∼ 5 – 10% at l ∼ 100 for the ISW/Rees-Sciama auto power spectrum, depending on the total neutrino mass M_ν, and becomes a factor of ∼ 4 for M_ν = 0.3 eV, at l ∼ 600, for the ISW/Rees-Sciama cross power with CMB-lensing. This effect should be taken into account for the correct estimation of the CMB temperature bispectrum in the presence of massive neutrinos.

In the context of classical mechanics, we study the conditions under which higher-order derivative theories can evade the so-called Ostrogradsky instability. More precisely, we consider general Lagrangians with second order time derivatives, of the form L(̈ϕ^a, ^a, ϕ^a; qⁱ, qⁱ) with a = 1,⋯,n and i = 1,⋯,m. For n = 1, assuming that the qⁱ's form a nondegenerate subsystem, we confirm that the degeneracy of the kinetic matrix eliminates the Ostrogradsky instability. The degeneracy implies, in the Hamiltonian formulation of the theory, the existence of a primary constraint, which generates a secondary constraint, thus eliminating the Ostrogradsky ghost. For n > 1, we show that, in addition to the degeneracy of the kinetic matrix, one needs to impose extra conditions to ensure the presence of a sufficient number of secondary constraints that can eliminate all the Ostrogradsky ghosts. When these conditions that ensure the disappearance of the Ostrogradsky instability are satisfied, we show that the Euler-Lagrange equations, which involve a priori higher order derivatives, can be reduced to a second order system.

In this paper, we consider rolling tachyon, with steep run-away type of potentials non-minimally coupled to massive neutrino matter. The coupling dynamically builds up at late times as neutrino matter turns non-relativistic. In case of scaling and string inspired potentials, we have shown that non-minimal coupling leads to minimum in the field potential. Given a suitable choice of model parameters, it is shown to give rise to late-time acceleration with the desired equation of state.

Inverse Compton (IC) scattering of the anisotropic CMB fluctuations off cosmic electron plasmas generates a polarization of the associated Sunyaev-Zel'dovich (SZ) effect. The polarized SZ effect has important applications in cosmology and in astrophysics of galaxy clusters. However, this signal has been studied so far mostly in the non-relativistic regime which is valid only in the very low electron temperature limit for a thermal electron population and, as such, has limited astrophysical applications. Partial attempts to extend this calculation to the IC scattering of a thermal electron plasma in the relativistic regime have been done but these cannot be applied to a more general or mildly relativistic electron distribution. In this paper we derive a general form of the SZ effect polarization that is valid in the full relativistic approach for both thermal and non-thermal electron plasmas, as well as for a generic combination of various electron population which can be co-spatially distributed in the environments of galaxy clusters or radiogalaxy lobes. We derive the spectral shape of the Stokes parameters induced by the IC scattering of every CMB multipole for both thermal and non-thermal electron populations, focussing in particular on the CMB quadrupole and octupole that provide the largest detectable signals in cosmic structures (like galaxy clusters). We found that the CMB quadrupole induced Stoke parameter Q is always positive with a maximum amplitude at a frequency ≈ 216 GHz which increases non-linearly with increasing cluster temperature. On the contrary, the CMB octupole induced Q spectrum shows a cross-over frequency which depends on the cluster electron temperature in a linear way, while it shows a non-linear dependence on the minimum momentum p₁ of a non-thermal power-law spectrum as well as a linear dependence on the power-law spectral index of the non-thermal electron population. We discuss some of the possibilities to disentangle the quadrupole-induced Q spectrum from the octupole-induced one which will allow to measure these important cosmological quantities through the SZ effect polarization at different cluster locations in the universe. We finally apply our model to the Bullet cluster and derive the visibility windows of the total, quandrupole-induced and octupole-induced Stoke parameter Q in the frequency ranges accessible to SKA, ALMA, MILLIMETRON and CORE++ experiments.

We study the evolution of cosmological perturbations in an anti-de-Sitter (AdS) bulk through a cosmological singularity by mapping the dynamics onto the boundary conformal fields theory by means of the AdS/CFT correspondence. We consider a deformed AdS space-time obtained by considering a time-dependent dilaton which induces a curvature singularity in the bulk at a time which we call t = 0, and which asymptotically approaches AdS both for large positive and negative times. The boundary field theory becomes free when the bulk curvature goes to infinity. Hence, the evolution of the fluctuations is under better controle on the boundary than in the bulk. To avoid unbounded particle production across the bounce it is necessary to smooth out the curvature singularity at very high curvatures. We show how the bulk cosmological perturbations can be mapped onto boundary gauge field fluctuations. We evolve the latter and compare the spectrum of fluctuations on the infrared scales relevant for cosmological observations before and after the bounce point. We find that the index of the power spectrum of fluctuations is the same before and after the bounce.

The detection of the magnetic type B-mode polarization is the main goal of future cosmic microwave background (CMB) experiments. In the standard model, the B-mode map is a strong non-gaussian field due to the CMB lensing component. Besides the two-point correlation function, the other statistics are also very important to dig the information of the polarization map. In this paper, we employ the Minkowski functionals to study the morphological properties of the lensed B-mode maps. We find that the deviations from Gaussianity are very significant for both full and partial-sky surveys. As an application of the analysis, we investigate the morphological imprints of the foreground residuals in the B-mode map. We find that even for very tiny foreground residuals, the effects on the map can be detected by the Minkowski functional analysis. Therefore, it provides a complementary way to investigate the foreground contaminations in the CMB studies.

We use time-sliced perturbation theory (TSPT) to give an accurate description of the infrared non-linear effects affecting the baryonic acoustic oscillations (BAO) present in the distribution of matter at very large scales. In TSPT this can be done via a systematic resummation that has a simple diagrammatic representation and does not involve uncontrollable approximations. We discuss the power counting rules and derive explicit expressions for the resummed matter power spectrum up to next-to leading order and the bispectrum at the leading order. The two-point correlation function agrees well with N-body data at BAO scales. The systematic approach also allows to reliably assess the shift of the baryon acoustic peak due to non-linear effects.

We study the 1-loop renormalization group equation running in the simplest singlet Majoron model constructed by us earlier to accommodate the dark radiation and dark matter content in the universe. A comprehensive numerical study was performed to explore the whole model parameter space. A smaller effective number of neutrinos △ N_eff∼ 0.05, or a Majoron decoupling temperature higher than the charm quark mass, is preferred. We found that a heavy scalar dark matter, ρ, of mass 1.5–4 TeV is required by the stability of the scalar potential and an operational type-I see-saw mechanism for neutrino masses. A neutral scalar, S, of mass in the 10–100 GeV range and its mixing with the standard model Higgs as large as 0.1 is also predicted. The dominant decay modes are S into b and/or ωω. A sensitive search will come from rare Z decays via the chain Z → S+ f, where f is a Standard Model fermion, followed by S into a pair of Majoron and/or b-quarks. The interesting consequences of dark matter bound state due to the sizable Sρ ρ-coupling are discussed as well. In particular, shower-like events with an apparent neutrino energy at M_ρ could contribute to the observed effective neutrino flux in underground neutrino detectors such as IceCube.

The interdependence of luminosity distance, D_L and angular diameter distance, D_A given by the distance duality relation (DDR) is very significant in observational cosmology. It is very closely tied with the temperature-redshift relation of Cosmic Microwave Background (CMB) radiation. Any deviation from η(z)≡ D_L/D_A (1+z)² =1 indicates a possible emergence of new physics. Our aim in this work is to check the consistency of these relations using a non-parametric regression method namely, LOESS with SIMEX. This technique avoids dependency on the cosmological model and works with a minimal set of assumptions. Further, to analyze the efficiency of the methodology, we simulate a dataset of 020 points of η (z) data based on a phenomenological model η(z)= (1+z). The error on the simulated data points is obtained by using the temperature of CMB radiation at various redshifts. For testing the distance duality relation, we use the JLA SNe Ia data for luminosity distances, while the angular diameter distances are obtained from radio galaxies datasets. Since the DDR is linked with CMB temperature-redshift relation, therefore we also use the CMB temperature data to reconstruct η (z). It is important to note that with CMB data, we are able to study the evolution of DDR upto a very high redshift  z = 2.418. In this analysis, we find no evidence of deviation from η=1 within a 1σ region in the entire redshift range used in this analysis (0 < z ⩽ 2.418).

Dwarf satellite galaxies are ideal laboratories for identifying particle Dark Matter signals. When setting limits on particle Dark Matter properties from null searches, it becomes however crucial the level at which the Dark Matter density profile within these systems is constrained by observations. In the limit in which the spherical Jeans equation is assumed to be valid for a given tracer stellar population, we study the solution of this equation having the Dark Matter mass profile as an output rather than as a trial parametric input. Within our new formulation, we address to what level dwarf spheroidal galaxies feature a reliable mass estimator. We assess then possible extrapolation of the density profiles in the inner regions and—keeping explicit the dependence on the orbital anisotropy profile of the tracer population—we derive general trends on the line-of-sight integral of the density profile squared, a quantity commonly dubbed J-factor and crucial to estimate fluxes from prompt Dark Matter pair annihilations.

Taking Ursa Minor as a study case among Milky Way satellites, we perform Bayesian inference using the available kinematical data for this galaxy. Contrary to all previous studies, we avoid marginalization over quantities poorly constrained by observations or by theoretical arguments. We find minimal J-factors to be about 2 to 4 times smaller than commonly quoted estimates, approximately relaxing by the same amount the limit on Dark Matter pair annihilation cross section from gamma-ray surveys of Ursa Minor. At the same time, if one goes back to a fixed trial parametric form for the density, e.g. using a NFW or Burkert profile, we show that the minimal J can hardly be reduced by more than a factor of 1.5.

We discuss dynamics of a model of an energy transfer between dark energy (DE) and dark matter (DM) . The energy transfer is determined by a non-conservation law resulting from a diffusion of dark matter in an environment of dark energy. The relativistic invariance defines the diffusion in a unique way. The system can contain baryonic matter and radiation which do not interact with the dark sector. We treat the Friedman equation and the conservation laws as a closed dynamical system. The dynamics of the model is examined using the dynamical systems methods for demonstration how solutions depend on initial conditions. We also fit the model parameters using astronomical observation: SNIa, H(z), BAO and Alcock-Paczynski test. We show that the model with diffuse DM-DE is consistent with the data.

Recently, a braking index for the pulsar PSR J1640-4631 has been measured. With a braking index of n = 3.15 ± 0.03, this pulsar has the highest braking index ever measured. As it is well known, a pure magnetic dipole brake yields n = 3, whereas a pure gravitational wave (GW) brake yields n = 5. Therefore, each of these mechanisms alone can not account for the braking index found for PSR J1640-4631. Here we consider in detail that such a braking index could be accounted for if the spindown model combines magnetic dipole and GW brakes. Then, we briefly discuss the detectability of this pulsar by aLIGO and the planned Einstein Telescope. In particular, we show that the amplitude of the GW that comes from our model is around a factor four lower than the amplitude modeled exclusively by GW energy loss. Another interesting outcome of our modeling is that it is possible to obtain the ellipticity from the braking index and other pulsar parameters.

We present a novel method to detect the effects of dynamical friction in observed galaxy clusters. Following accretion into clusters, massive satellite galaxies will backsplash to systematically smaller radii than less massive satellites, an effect that may be detected by stacking the number density profiles of galaxies around clusters. We show that this effect may be understood using a simple toy model which reproduces the trends with halo properties observed in simulations. We search for this effect using SDSS redMaPPer clusters with richness 010 < λ < 2, and find that bright (M_i < −21.5) satellites have smaller splashback radii than fainter (0M_i > −2) satellites at 99% confidence.

Gravitational wave (GW) source counts have been recently shown to be able to test how gravitational radiation propagates with the distance from the source. Here, we extend this formalism to cosmological scales, i.e. the high redshift regime, and we discuss the complications of applying this methodology to high redshift sources. We also allow for models with compactified extra dimensions like in the Kaluza-Klein model. Furthermore, we also consider the case of intermediate redshifts, i.e. 0 < z ≲ 1, where we show it is possible to find an analytical approximation for the source counts dN/d(S/N). This can be done in terms of cosmological parameters, such as the matter density Ω_m,0 of the cosmological constant model or the cosmographic parameters for a general dark energy model. Our analysis is as general as possible, but it depends on two important factors: a source model for the black hole binary mergers and the GW source to galaxy bias. This methodology also allows us to obtain the higher order corrections of the source counts in terms of the signal-to-noise S/N. We then forecast the sensitivity of future observations in constraining GW physics but also the underlying cosmology by simulating sources distributed over a finite range of signal-to-noise with a number of sources ranging from 10 to 500 sources as expected from future detectors. We find that with 500 events it will be possible to provide constraints on the matter density parameter at present Ω_m,0 on the order of a few percent and with the precision growing fast with the number of events. In the case of extra dimensions we find that depending on the degeneracies of the model, with 500 events it may be possible to provide stringent limits on the existence of the extra dimensions if the aforementioned degeneracies can be broken.

This is the third paper in a series establishing a quantitative relation between inflationary scalar field potential landscapes and the relic perturbations left by the collision between bubbles produced during eternal inflation. We introduce a new method for computing cosmological observables from numerical relativity simulations of bubble collisions in one space and one time dimension. This method tiles comoving hypersurfaces with locally-perturbed Friedmann-Robertson-Walker coordinate patches. The method extends previous work, which was limited to the spacetime region just inside the future light cone of the collision, and allows us to explore the full bubble-collision spacetime. We validate our new methods against previous work, and present a full set of predictions for the comoving curvature perturbation and local negative spatial curvature produced by identical and non-identical bubble collisions, in single scalar field models of eternal inflation. In both collision types, there is a non-zero contribution to the spatial curvature and cosmic microwave background quadrupole. Some collisions between non-identical bubbles excite wall modes, giving extra structure to the predicted temperature anisotropies. We comment on the implications of our results for future observational searches. For non-identical bubble collisions, we also find that the surfaces of constant field can readjust in the presence of a collision to produce spatially infinite sections that become nearly homogeneous deep into the region affected by the collision. Contrary to previous assumptions, this is true even in the bubble into which the domain wall is accelerating.

The Beyond Horndeski class of alternative gravity theories allow for Self-accelerating de-Sitter cosmologies with no need for a cosmological constant. This makes them viable alternatives to ΛCDM and so testing their small-scale predictions against General Relativity is of paramount importance. These theories generically predict deviations in both the Newtonian force law and the gravitational lensing of light inside extended objects. Therefore, by simultaneously fitting the X-ray and lensing profiles of galaxy clusters new constraints can be obtained. In this work, we apply this methodology to the stacked profiles of 58 high-redshift (0.1 < z < 1.2) clusters using X-ray surface brightness profiles from the XMM Cluster Survey and weak lensing profiles from CFHTLenS. By performing a multi-parameter Markov chain Monte Carlo analysis, we are able to place new constraints on the parameters governing deviations from Newton's law Υ₁ = −0.11^+0.93_−0.67 and light bending Υ₂ = −0.22^+1.22_−1.19. Both constraints are consistent with General Relativity, for which Υ₁ = Υ₂ = 0. We present here the first observational constraints on Υ₂, as well as the first extragalactic measurement of both parameters.

We present a generalization of the effective field theory (EFT) formalism for dark energy and modified gravity models to include operators with higher order spatial derivatives. This allows the extension of the EFT framework to a wider class of gravity theories such as Hořava gravity. We present the corresponding extended action, both in the EFT and the Arnowitt-Deser-Misner (ADM) formalism, and proceed to work out a convenient mapping between the two, providing a self contained and general procedure to translate a given model of gravity into the EFT language at the basis of the Einstein-Boltzmann solver EFTCAMB. Putting this mapping at work, we illustrate, for several interesting models of dark energy and modified gravity, how to express them in the ADM notation and then map them into the EFT formalism. We also provide for the first time, the full mapping of GLPV models into the EFT framework. We next perform a thorough analysis of the physical stability of the generalized EFT action, in absence of matter components. We work out viability conditions that correspond to the absence of ghosts and modes that propagate with a negative speed of sound in the scalar and tensor sector, as well as the absence of tachyonic modes in the scalar sector. Finally, we extend and generalize the phenomenological basis in terms of α-functions introduced to parametrize Horndeski models, to cover all theories with higher order spatial derivatives included in our extended action. We elaborate on the impact of the additional functions on physical quantities, such as the kinetic term and the speeds of propagation for scalar and tensor modes.

We derive the exact third-order analytic solution of the matter density fluctuation in the proper-time hypersurface in a ΛCDM universe, accounting for the explicit time-dependence and clarifying the relation to the initial condition. Furthermore, we compare our analytic solution to the previous calculation in the comoving gauge, and to the standard Newtonian perturbation theory by providing Fourier kernels for the relativistic effects. Our results provide an essential ingredient for a complete description of galaxy bias in the relativistic context.

We perform a Hamiltonian analysis of a large class of scalar-tensor Lagrangians which depend quadratically on the second derivatives of a scalar field. By resorting to a convenient choice of dynamical variables, we show that the Hamiltonian can be written in a very simple form, where the Hamiltonian and the momentum constraints are easily identified. In the case of degenerate Lagrangians, which include the Horndeski and beyond Horndeski quartic Lagrangians, our analysis confirms that the dimension of the physical phase space is reduced by the primary and secondary constraints due to the degeneracy, thus leading to the elimination of the dangerous Ostrogradsky ghost. We also present the Hamiltonian formulation for nondegenerate theories and find that they contain four degrees of freedom, including a ghost, as expected. We finally discuss the status of the unitary gauge from the Hamiltonian perspective.

We estimate the thermal masses and damping rates of active (m < eV) and sterile (M ∼ GeV) neutrinos with thermal momenta k∼ 3T at temperatures below the electroweak crossover (5 GeV < T < 160 GeV) . These quantities fix the equilibration or ``washout'' rates of Standard Model lepton number densities. Sterile neutrinos interact via direct scatterings mediated by Yukawa couplings, and via their overlap with active neutrinos. Including all leading-order reactions we find that the washout rate generally exceeds the Hubble rate for 5 GeV < T < 30 GeV . Therefore it is challenging to generate a large lepton asymmetry facilitating dark matter computations operating at T < 5 GeV, whereas the generation of a baryon asymmetry at T > 130 GeV remains an option. Our differential rates are tabulated in a form suitable for studies of specific scenarios with given neutrino Yukawa matrices.

We investigate the inflaton decay and reheating period after the end of inflation in the non-minimal derivative coupling (NDC) model with chaotic potential. In general, this model is known to provide an enhanced slow-roll inflation caused by gravitationally enhanced friction. We find violent oscillations of Hubble parameter which induces oscillations of the sound speed squared, implying the Lagrangian instability of curvature perturbation ζ under the comoving gauge φ = 0. Also, it is shown that the curvature perturbation blows up at  = 0, leading to the breakdown of the comoving gauge at  = 0. Therefore, we use the Newtonian gauge to perform the perturbation analysis where the Newtonian potential is employed as a physical variable. The curvature perturbation is not considered as a physical variable which describes a relevant perturbation during reheating.

Small-scale structure is studied in the context of dissipative dark matter, arising for instance in models with a hidden unbroken Abelian sector, so that dark matter couples to a massless dark photon. The dark sector interacts with ordinary matter via gravity and photon-dark photon kinetic mixing. Mirror dark matter is a theoretically constrained special case where all parameters are fixed except for the kinetic mixing strength, . In these models, the dark matter halo around spiral and irregular galaxies takes the form of a dissipative plasma which evolves in response to various heating and cooling processes. It has been argued previously that such dynamics can account for the inferred cored density profiles of galaxies and other related structural features. Here we focus on the apparent deficit of nearby small galaxies (``missing satellite problem"), which these dissipative models have the potential to address through small-scale power suppression by acoustic and diffusion damping. Using a variant of the extended Press-Schechter formalism, we evaluate the halo mass function for the special case of mirror dark matter. Considering a simplified model where M_baryons ∝ M_halo, we relate the halo mass function to more directly observable quantities, and find that for  ≈ 2 × 10⁻¹⁰ such a simplified description is compatible with the measured galaxy luminosity and velocity functions. On scales M_halo ≲ 10⁸ M_⊙, diffusion damping exponentially suppresses the halo mass function, suggesting a nonprimordial origin for dwarf spheroidal satellite galaxies, which we speculate were formed via a top-down fragmentation process as the result of nonlinear dissipative collapse of larger density perturbations. This could explain the planar orientation of satellite galaxies around Andromeda and the Milky Way.

The covariant and gauge invariant calculation of the current expectation value in the homogeneous electric field in 1+3 dimensional de Sitter spacetime is shown. The result accords with previous work obtained by using adiabatic subtraction scheme. We therefore conclude the counterintuitive behaviors of the current in the infrared (IR) regime such as IR hyperconductivity and negative current are not artifacts of the renormalization scheme, but are real IR effects of the spacetime.

There is ample evidence from rotation curves that dark matter halos around disk galaxies have nontrivial dynamics. Of particular significance are: a) the cored dark matter profile of disk galaxies, b) correlations of the shape of rotation curves with baryonic properties, and c) Tully-Fisher relations. Dark matter halos around disk galaxies may have nontrivial dynamics if dark matter is strongly self interacting and dissipative. Multicomponent hidden sector dark matter featuring a massless `dark photon' (from an unbroken dark U(1) gauge interaction) which kinetically mixes with the ordinary photon provides a concrete example of such dark matter. The kinetic mixing interaction facilitates halo heating by enabling ordinary supernovae to be a source of these `dark photons'. Dark matter halos can expand and contract in response to the heating and cooling processes, but for a sufficiently isolated halo could have evolved to a steady state or `equilibrium' configuration where heating and cooling rates locally balance. This dynamics allows the dark matter density profile to be related to the distribution of ordinary supernovae in the disk of a given galaxy. In a previous paper a simple and predictive formula was derived encoding this relation. Here we improve on previous work by modelling the supernovae distribution via the measured UV and Hα fluxes, and compare the resulting dark matter halo profiles with the rotation curve data for each dwarf galaxy in the LITTLE THINGS sample. The dissipative dark matter concept is further developed and some conclusions drawn.

We explore Schwinger effect of spin 1/2 charged particles with static electric field in 1+3 dimensional de Sitter spacetime. We analytically calculate the vacuum expectation value of the spinor current which is induced by the produced particles in the electric field. The renormalization is performed with the adiabatic subtraction scheme. We find that the current becomes negative, namely it flows in the direction opposite to the electric field, if the electric field is weaker than a certain threshold value depending on the fermion mass, which is also known to happen in the case of scalar charged particles in 1+3 de Sitter spacetime. Contrary to the scalar case, however, the IR hyperconductivity is absent in the spinor case.

We consider the steady-state regime describing the density profile of a dark matter halo, if dark matter is treated as a Bose-Einstein condensate. We first solve the fluid equation for ``canonical'' cold dark matter, obtaining a class of density profiles which includes the Navarro-Frenk-White profile, and which diverge at the halo core. We then solve numerically the equation obtained when an additional ``quantum pressure'' term is included in the computation of the density profile. The solution to this latter case is finite at the halo core, possibly avoiding the ``cuspy halo problem'' present in some cold dark matter theories. Within the model proposed, we predict the mass of the cold dark matter particle to be of the order of M_χ c² ≈ 10⁻²⁴ eV, which is of the same order of magnitude as that predicted in ultra-light scalar cold dark matter models. Finally, we derive the differential equation describing perturbations in the density and the pressure of the dark matter fluid.

In this article we study Kaluza-Klein (KK) dimensional reduction of massive Abelian gauge theories with charged matter fields on a circle. Since local gauge transformations change position dependence of the charged fields, the decomposition of the charged matter fields into KK modes is gauge dependent. While whole KK mass spectrum is independent of the gauge choice, the mode number depends on the gauge. The masses of the KK modes also depend on the field value of the zero-mode of the extra dimensional component of the gauge field. In particular, one of the KK modes in the KK tower of each massless 5D charged field becomes massless at particular values of the extra-dimensional component of the gauge field. When the extra-dimensional component of the gauge field is identified with the inflaton, this structure leads to recursive cosmological particle productions.

Strong field gravitational lensings are dramatically disparate from those in the weak field by representing relativistic images due to light winds one to infinity loops around a lens before escaping. We study such a lensing caused by a charged Galileon black hole, which is expected to have possibility to evade no-hair theorem. We calculate the angular separations and time delays between different relativistic images of the charged Galileon black hole. All these observables can potentially be used to discriminate a charged Galileon black hole from others. We estimate the magnitudes of these observables for the closest supermassive black hole Sgr A*. The strong field lensing observables of the charged Galileon black hole can be close to those of a tidal Reissner-Nordström black hole or those of a Reissner-Nordström black hole. It will be helpful to distinguish these black holes if we can separate the outermost relativistic images and determine their angular separation, brightness difference and time delay, although it requires techniques beyond the current limit.

We study inflation in the Brans-Dicke gravity as a special model of the scalar-tensor gravity. We obtain the inflationary observables containing the scalar spectral index, the tensor-to-scalar ratio, the running of the scalar spectral index and the equilateral non-Gaussianity parameter in terms of the general form of the potential in the Jordan frame. Then, we compare the results for various inflationary potentials in light of the Planck 2015 data. Our study shows that in the Brans-Dicke gravity, the power-law, inverse power-law and exponential potentials are ruled out by the Planck 2015 data. But, the hilltop, Higgs, Coleman-Weinberg and natural potentials can be compatible with Planck 2015 TT,TE,EE+lowP data at 95% CL. Moreover, the D-brane, SB SUSY and displaced quadratic potentials can be in well agreement with the observational data since their results can lie inside the 68% CL region of Planck 2015 TT,TE,EE+lowP data.

We investigate a simple setup in which an excess in the di-photon invariant mass distribution around 750 GeV, as seen by the ATLAS and CMS collaborations, is originated through a pair of collimated photon pairs. In this framework a scalar state s decays into two light pseudo-Goldstone bosons a, each of which subsequently decays into a pair of collimated photons which are misidentified as a single photon. In a minimal context of spontaneous symmetry breaking, we show that coupling a complex scalar field Φ=(s+ia)/√2 to a fermionic dark matter candidate χ, also responsible for generating its mass, allows for the correct relic density in a large region of the parameter space, while not being excluded by the direct or indirect detection experiments. Moreover, the correct relic abundance can naturally co-exist with a relatively large width for the resonant field s.

We extend the DFSZ invisible axion model to simultaneously explain small Dirac neutrino masses and cosmic matter-antimatter asymmetry. After the Peccei-Quinn and electroweak symmetry breaking, the effective Yukawa couplings of the Dirac neutrinos to the standard model Higgs scalar can be highly suppressed by the ratio of the vacuum expectation value of an iso-triplet Higgs scalar over the masses of some heavy gauge-singlet fermions, iso-doublet Higgs scalars or iso-triplet fermions. The iso-triplet fields can carry a zero or nonzero hypercharge. Through the decays of the heavy gauge-singlet fermions, iso-doublet scalars or iso-triplet fermions, we can obtain a lepton asymmetry in the left-handed leptons and an opposite lepton asymmetry in the right-handed neutrinos. Since the right-handed neutrinos do not participate in the sphaleron processes, the left-handed lepton asymmetry can be partially converted to a baryon asymmetry.

We construct a class of generalized non-local gravity (GNLG) model which is the modified theory of general relativity (GR) obtained by adding a term m²ⁿ⁻² R□⁻ⁿR to the Einstein-Hilbert action. Concretely, we not only study the gravitational equation for the GNLG model by introducing auxiliary scalar fields, but also analyse the classical stability and examine the cosmological consequences of the model for different exponent n. We find that the half of the scalar fields are always ghost-like and the exponent n must be taken even number for a stable GNLG model. Meanwhile, the model spontaneously generates three dominant phases of the evolution of the universe, and the equation of state parameters turn out to be phantom-like. Furthermore, we clarify in another way that exponent n should be even numbers by the spherically symmetric static solutions in Newtonian gauge. It is worth stressing that the results given by us can include ones in refs. [28, 34] as the special case of n=2.

Generalized fluid equations, using sound speed c_eff² and viscosity c_vis² as effective parameters, provide a convenient phenomenological formalism for testing the relic neutrino "null hypothesis," i.e. that that neutrinos are relativistic and free-streaming prior to recombination. In this work, we relax the relativistic assumption and ask "to what extent can the generalized fluid equations accommodate finite neutrino mass?" We consider both the mass of active neutrinos, which are largely still relativistic at recombination m² / T² ∼ 0.2, and the effect of a semi-relativistic sterile component. While there is no one-to-one mapping between mass/mixing parameters and c_eff² and c_vis², we demonstrate that the existence of a neutrino mass could induce a bias to measurements of c_eff² and c_vis² at the level of 0.01 m² / T² ∼ 10⁻³.

Observations of the kinetic Sunyaev-Zel'dovich (kSZ) effect measure the density-weighted velocity field, a potentially powerful cosmological probe. This paper presents an analytical method to predict the power spectrum and two-point correlation function of the density-weighted velocity in redshift space, the direct observables in kSZ surveys. We show a simple relation between the density power spectrum and the density-weighted velocity power spectrum that holds for both dark matter and halos. Using this relation, we can then extend familiar perturbation expansion techniques to the kSZ power spectrum. One of the most important features of density-weighted velocity statistics in redshift space is the change in sign of the cross-correlation between the density and density-weighted velocity at mildly small scales due to nonlinear redshift space distortions. Our model can explain this characteristic feature without any free parameters. As a result, our results can precisely predict the non-linear behavior of the density-weighted velocity field in redshift space up to ∼ 30 h⁻¹ Mpc for dark matter particles at the redshifts of z=0.0, 0.5, and 1.0.

Please see the PDF file for details.