Table of contents

Chameleon scalar fields can screen their associated fifth forces from detection by changing their mass with the local density. These models are an archetypal example of a screening mechanism, and have become an important target for both cosmological surveys and terrestrial experiments. In particular there has been much recent interest in searching for chameleon fifth forces in the laboratory. It is known that the chameleon force is less screened around non-spherical sources, but only the field profiles around a few simple shapes are known analytically. In this work we introduce a numerical code that solves for the chameleon field around arbitrary shapes with azimuthal symmetry placed in a spherical vacuum chamber. We find that deviations from spherical symmetry can increase the chameleon acceleration experienced by a test particle, and that the least screened objects are those which minimize some internal dimension. For the shapes considered in this work, keeping the mass, density and background environment fixed, the accelerations due to the source varied by a factor of ∼ 3.

AMS-02 has reached the sensitivity to probe canonical thermal WIMPs by their annihilation into antiprotons. Due to the high precision of the data, uncertainties in the astrophysical background have become the most limiting factor for indirect dark matter detection. In this work we systematically quantify and—where possible—reduce uncertainties in the antiproton background. We constrain the propagation of charged cosmic rays through the combination of antiproton, B/C and positron data. Cross section uncertainties are determined from a wide collection of accelerator data and are—for the first time ever—fully taken into account. This allows us to robustly constrain even subdominant dark matter signals through their spectral properties. For a standard NFW dark matter profile we are able to exclude thermal WIMPs with masses up to 570 GeV which annihilate into bottom quarks. While we confirm a reported excess compatible with dark matter of mass around 80 GeV, its local (global) significance only reaches 2.2 σ (1.1 σ) in our analysis.

We consider the freeze-in production of 7 keV axino dark matter (DM) in the supersymmetric Dine-Fischler-Srednicki-Zhitnitsky (DFSZ) model in light of the 3.5 keV line excess. The warmness of such 7 keV DM produced from the thermal bath, in general, appears in tension with Ly-α forest data, although a direct comparison is not straightforward. This is because the Ly-α forest constraints are usually reported on the mass of the conventional warm dark matter (WDM), where large entropy production is implicitly assumed to occur in the thermal bath after WDM particles are decoupled. The phase space distribution of freeze-in axino DM varies depending on production processes and axino DM may alleviate the tension with the tight Ly-α forest constraints. By solving the Boltzmann equation, we first obtain the resultant phase space distribution of axinos produced by 2-body decay, 3-body decay, and 2-to-2 scattering respectively. The reduced collision term and resultant phase space distribution are useful for studying other freeze-in scenarios as well. We then calculate the resultant linear matter power spectra for such axino DM and directly compare them with the linear matter power spectra for the conventional WDM . In order to demonstrate realistic axino DM production, we consider benchmark points with Higgsino next-to-light supersymmetric particle (NLSP) and wino NLSP. In the case of Higgsino NLSP, the phase space distribution of axinos is colder than that in the conventional WDM case, so the most stringent Ly-α forest constraint can be evaded with mild entropy production from saxion decay inherent in the supersymmetric DFSZ axion model.

Gravitationally collapsed objects are known to be biased tracers of an underlying density contrast. Using symmetry arguments, generalised biasing schemes have recently been developed to relate the halo density contrast δ_h with the underlying density contrast δ, divergence of velocity θ and their higher-order derivatives. This is done by constructing invariants such as s, t, ψ,η. We show how the generating function formalism in Eulerian standard perturbation theory (SPT) can be used to show that many of the additional terms based on extended Galilean and Lifshitz symmetry actually do not make any contribution to the higher-order statistics of biased tracers. Other terms can also be drastically simplified allowing us to write the vertices associated with δ_h in terms of the vertices of δ and θ, the higher-order derivatives and the bias coefficients. We also compute the cumulant correlators (CCs) for two different tracer populations. These perturbative results are valid for tree-level contributions but at an arbitrary order. We also take into account the stochastic nature bias in our analysis. Extending previous results of a local polynomial model of bias, we express the one-point cumulants _N and their two-point counterparts, the CCs i.e. _pq, of biased tracers in terms of that of their underlying density contrast counterparts. As a by-product of our calculation we also discuss the results using approximations based on Lagrangian perturbation theory (LPT).

In a confined system of multiple Fermions, the particles are forced into high energy levels by the Pauli Exclusion Principle. We refer to this system as a Pauli tower. We pursue the investigation of a model for sub-Chandrasekhar supernovae Ia explosions (SNIa) in which the energy stored in the Pauli tower is released to trigger a nuclear deflagration. The simplest physical model for such a degeneracy breakdown and collapse of the Pauli tower is a phase transition to an exactly supersymmetric state in which the scalar partners of protons, neutrons, and leptons become degenerate with the familiar fermions of our world as in the supersymmetric standard model with susy breaking parameters relaxed to zero. We focus on the ability of the susy phase transition model to fit the total SNIa rate as well as the delay time distribution of SNIa after the birth of a progenitor white dwarf. We also study the ejected mass distribution and its correlation with delay time. Finally, we discuss the expected SNIa remnant in the form of a black hole of roughly Jupiter mass and the prospects for detecting such remnants.

The standard model of cosmology with nearly Gaussian, isotropic, scale invariant and adiabatic initial conditions describes the cosmological observations well. However, the study of any deviation from the mentioned conditions will open up a new horizon to the physics of early universe. In this work, we study the effect of the oscillatory and step-like features in potentials of inflationary models in late time large scale structure observations. Mainly we study the matter power spectrum, number density of the structures, dark matter halo bias and specifically CMB lensing. We show that the oscillatory models can introduce some degeneracy with late time effects on BAO scale. We also conclude that high frequency oscillatory models which are favored by Planck data do not have significant effect on the non linear structure formation. Finally we show that inflationary models with step functions which deviates from the standard model in small scales l ⩽ 1 Mpc can be constrained by future experiments via CMB lensing. We propose the idea that CMB lensing is a bias independent observation which can be used as a small scale physics probe due to distribution of the lenses in low redshifts. Meantime this model can alter the prediction of the cosmological model for the number density of small structures and can be used as a probable explanation for galactic scale crisis of ΛCDM.

The effect of Quantum Gravity (QG) may bring a tiny light speed variation as v(E)=c(1−E/E_LV), where E is the photon energy and E_LV is a Lorentz violation scale. A remarkable regularity was suggested in previous studies to look for the light speed variation from high energy photon events of Gamma Ray Bursts (GRBs). We provide a general analysis on the data of 25 bright GRBs observed by the Fermi Gamma-ray Space Telescope (FGST). Such method allows a completed scan over all possibilities in a more clean and impartial way without any bias compared to previous intuitive analysis. The results show that with the increase in the intrinsic energies of photons, such regularity truly emerges and gradually becomes significant. For photons with intrinsic energies higher than 40 GeV, the regularity exists at a significance of 3–5 σ with E_LV=3.6× 10¹⁷ GeV determined by the GRB data.

Anisotropic inflation is an interesting model with an \U(1) gauge field and it predicts the statistical anisotropy of the curvature perturbation characterized by a parameter g_*. However, we find that the background gauge field does not follow the classical attractor solution due to the stochastic effect. We develop the stochastic formalism of a vector field and solve Langevin and Fokker-Planck equations. It is shown that this model is excluded by the CMB constraint g_*≤ 10⁻² with a high probability about 99.999%.

We consider the cosmological consequences of a special scalar-tensor-vector theory of gravity, known as MOG (for MOdified Gravity), proposed to address the dark matter problem. This theory introduces two scalar fields G(x) and μ(x), and one vector field ϕ_α(x), in addition to the metric tensor. We set the corresponding self-interaction potentials to zero, as in the standard form of MOG. Then using the phase space analysis in the flat Friedmann-Robertson-Walker background, we show that the theory possesses a viable sequence of cosmological epochs with acceptable time dependency for the cosmic scale factor. We also investigate MOG's potential as a dark energy model and show that extra fields in MOG cannot provide a late time accelerated expansion. Furthermore, using a dynamical system approach to solve the non-linear field equations numerically, we calculate the angular size of the sound horizon, i.e. θ_s, in MOG. We find that 8× 10⁻³rad<θ_s<8.2× 10⁻³ rad  which is way outside the current observational bounds. Finally, we generalize MOG to a modified form called mMOG, and we find that mMOG passes the sound-horizon constraint. However, mMOG also cannot be considered as a dark energy model unless one adds a cosmological constant, and more importantly, the matter dominated era is still slightly different from the standard case.

We argue that the acoustic damping of the matter power spectrum is not a generic feature of the kinetic decoupling of dark matter, but even the enhancement (overshooting) can be realized depending on the nature of the kinetic decoupling when compared to that in the standard cold dark matter model. We consider a model that exhibits a sudden kinetic decoupling and investigate cosmological perturbations in the standard cosmological background numerically in the model. We also give an analytic discussion in a simplified setup. Our results indicate that the nature of the kinetic decoupling could have a great impact on small scale density perturbations.

Studying the physics of compact objects in modified theories of gravity is important for understanding how future observations can test alternatives to General Relativity. We consider a subset of vector-tensor Galileon theories of gravity characterized by new symmetries, which can prevent the propagation of the vector longitudinal polarization, even in absence of Abelian gauge invariance. We investigate new spherically symmetric and slowly rotating solutions for these systems, including an arbitrary matter Lagrangian. We show that, under certain conditions, there always exist stealth configurations whose geometry coincides with solutions of Einstein gravity coupled with the additional matter. Such solutions have a non-trivial profile for the vector field, characterized by independent integration constants, which extends to asymptotic infinity. We interpret our findings in terms of the symmetries and features of the original vector-tensor action, and on the number of degrees of freedom that it propagates. These results are important to eventually describe gravitationally bound configurations in modified theories of gravity, such as black holes and neutron stars, including realistic matter fields forming or surrounding the object.

In this paper, we study the effects of general relativistic corrections on the observed galaxy power spectrum in thawing class of cubic Galileon model with linear potential that preserves the shift symmetry. In this scenario, the observed galaxy power spectrum differs from the standard matter power spectrum mainly due to redshift space distortion (RSD) factor and relativistic effects. The RSD term enhances the matter power spectrum both at larger and smaller scales whereas the relativistic terms further enhance the matter power spectrum only at larger scales. In comparison with ΛCDM, the observed galaxy power spectrum is always suppressed at large scales in this scenario although this suppression is always small compared to the canonical quintessence scenario. We also study the small scale non linear power spectrum for Galileon model using the publicly available HMcode. Our study shows that at the nonlinear regime, the small scale power is enhanced in Galileon model compared to the ΛCDM . But this enhancement in Galileon model is lower than the quintessence case.

The recent detection of the gravitational wave signal GW170817 together with an electromagnetic counterpart GRB 170817A from the merger of two neutron stars puts a stringent bound on the tensor propagation speed. This constraint can be automatically satisfied in the framework of massive gravity. In this work, we consider a general SO(3)-invariant massive gravity with five propagating degrees of freedom and derive the conditions for the absence of ghosts and Laplacian instabilities in the presence of a matter perfect fluid on the flat Friedmann-Lemaître-Robertson-Walker (FLRW) cosmological background. The graviton potential containing the dependence of three-dimensional metrics and a fiducial metric coupled to a temporal scalar field gives rise to a scenario of the late-time cosmic acceleration in which the dark energy equation of state w_DE is equivalent to −1 or varies in time. We find that the deviation from the value w_DE=−1 provides important contributions to the quantities associated with the stability conditions of tensor, vector, and scalar perturbations. In concrete models, we study the dynamics of dark energy arising from the graviton potential and show that there exist viable parameter spaces in which neither ghosts nor Laplacian instabilities are present for both w_DE>−1 and w_DE<−1. We also generally obtain the effective gravitational coupling G_eff with non-relativistic matter as well as the gravitational slip parameter η_s associated with the observations of large-scale structures and weak lensing. We show that, apart from a specific case, the two quantities G_eff and η_s are similar to those in general relativity for scalar perturbations deep inside the sound horizon.

We analyze Bose-Einstein condensates on three types of spherically symmetric and static charged black-hole spacetimes: the Reissner-Nordström spacetime, Hoffmann's Born-Infeld black-hole spacetime, and the regular Ayón-Beato-García spacetime. The Bose-Einstein condensate is modeled in terms of a massive scalar field that satisfies a Klein-Gordon equation with a self-interaction term. The scalar field is assumed to be uncharged and not self-gravitating. If the mass parameter of the scalar field is chosen sufficiently small, there are quasi-bound states of the scalar field that may be interpreted as dark matter clouds. We estimate the size and the total energy of such clouds around charged supermassive black holes and we investigate if their observable features can be used for discriminating between the different types of charged black holes.

We provide a fast algorithm to diagnose any directional dependence in the cosmological parameters by calculating maps of local cosmological parameter estimates and their joint errors. The technique implements a fast quadratic estimator technique based on Wiener filtering and convolution of the sky with a patch shape. It uses only three map-resolution spherical harmonic transforms per parameter and applies to any data set with full sky or a partial sky coverage. We apply this method to Planck SMICA-2015 and obtain fluctuation map for six cosmological parameters. Our estimate shows that the Planck data is consistent with a single global value of the cosmological parameters and is not influenced by any severe local contaminations. This method is applicable also to other angular or 3D data sets of future missions to scrutinize any local variation in the cosmological parameters.

We study the power spectrum of quasi-single field inflation where strong coupling is considered. The contribution from the massive propagator can be divided into local and non-local contributions. The local one is the leading contribution and is power-law suppressed as a function of mass, while the non-local contribution is exponentially suppressed in the large mass limit. For the local contribution, it is possible to use the effective field theory approach to study the power spectrum in the strongly coupled region of the parameter space. For the non-local contribution, we developed a partial effective field theory method to simplify the calculation: when there are multiple massive propagators, one can fully compute it after integrating out all but one massive propagator by effective field theory. The result retains the "standard clock" signal, which is interesting for probing the expansion history of the primordial universe and the physics of a "cosmological collider". The error involved compared to the full calculation is power law suppressed by the effective mass of the heavy field.

We consider cosmological dynamics in the theory of gravity with the scalar field possessing the nonminimal kinetic coupling to curvature given as η G^μνϕ_,μϕ_,ν, where η is an arbitrary coupling parameter, and the scalar potential V(ϕ) which assumed to be as general as possible. With an appropriate dimensionless parametrization we represent the field equations as an autonomous dynamical system which contains ultimately only one arbitrary function χ (x)= 8 π | η | V(x/√8 π) with x=√8 πϕ. Then, assuming the rather general properties of χ(x), we analyze stationary points and their stability, as well as all possible asymptotical regimes of the dynamical system. It has been shown that for a broad class of χ(x) there exist attractors representing three accelerated regimes of the Universe evolution, including de Sitter expansion (or late-time inflation), the Little Rip scenario, and the Big Rip scenario. As the specific examples, we consider a power-law potential V(ϕ)=M⁴(ϕ/ϕ₀)^σ, Higgs-like potential V(ϕ)=λ/4(ϕ²−ϕ₀²)², and exponential potential V(ϕ)=M⁴ e^−ϕ/ϕ₀.

Cosmological inflation provides the initial conditions for the structure formation. However, the origin of large-scale magnetic fields can not be addressed in this framework. The key issue for this long-standing problem is the conformal invariance of the electromagnetic (EM) field in 4-D. While many approaches have been proposed in the literature for breaking conformal invariance of the EM action, here, we provide a completely new way of looking at the modifications to the EM action and generation of primordial magnetic fields during inflation. We explicitly construct a higher derivative EM action that breaks conformal invariance by demanding three conditions—theory be described by vector potential A^μ and its derivatives, Gauge invariance be satisfied, and equations of motion be linear in second derivatives of vector potential. The unique feature of our model is that appreciable magnetic fields are generated at small wavelengths while tiny magnetic fields are generated at large wavelengths that are consistent with current observations.

We consider scalar field models of dark energy interacting with dark matter through a coupling proportional to the contraction of the four-derivative of the scalar field with the four-velocity of the dark matter fluid. The coupling is realized at the Lagrangian level employing the formalism of Scalar-Fluid theories, which use a consistent Lagrangian approach for relativistic fluid to describe dark matter. This framework produces fully covariant field equations, from which we can derive unequivocal cosmological equations at both background and linear perturbations levels. The background evolution is analyzed in detail applying dynamical systems techniques, which allow us to find the complete asymptotic behavior of the universe given any set of model parameters and initial conditions. Furthermore we study linear cosmological perturbations investigating the growth of cosmic structures within the quasi-static approximation. We find that these interacting dark energy models give rise to interesting phenomenological dynamics, including late-time transitions from dark matter to dark energy domination, matter and accelerated scaling solutions and dynamical crossing of the phantom barrier. Moreover we obtain possible deviations from standard ΛCDM behavior at the linear perturbations level, which have an impact on the dynamics of structure formation and might provide characteristic observational signatures.

Recently there has been much interest in the spatial distribution of light scalar dark matter, especially axions, throughout the universe. When the local gravitational interactions between the scalar modes are sufficiently rapid, it can cause the field to re-organize into a BEC of gravitationally bound clumps. While these clumps are stable when only gravitation is included, the picture is complicated by the presence of the axion's attractive self-interactions, which can potentially cause the clumps to collapse. Here we perform a detailed stability analysis to determine under what conditions the clumps are stable. In this paper we focus on spherical configurations, leaving aspherical configurations for future work. We identify branches of clump solutions of the axion-gravity-self-interacting system and study their stability properties. We find that clumps that are (spatially) large are stable, while clumps that are (spatially) small are unstable and may collapse. Furthermore, there is a maximum number of particles that can be in a clump. We map out the full space of solutions, which includes quasi-stable axitons, and clarify how a recent claim in the literature of a new ultra-dense branch of stable solutions rests on an invalid use of the non-relativistic approximation. We also consider repulsive self-interactions that may arise from a generic scalar dark matter candidate, finding a single stable branch that extends to arbitrary particle number.

We study inflation in models with many interacting fields subject to randomly generated scalar potentials. We use methods from non-equilibrium random matrix theory to construct the potentials and an adaption of the `transport method' to evolve the two-point correlators during inflation. This construction allows, for the first time, for an explicit study of models with up to 100 interacting fields supporting a period of `approximately saddle-point' inflation. We determine the statistical predictions for observables by generating over 30,000 models with 2–100 fields supporting at least 60 efolds of inflation. These studies lead us to seven lessons: i) Manyfield inflation is not single-field inflation, ii) The larger the number of fields, the simpler and sharper the predictions, iii) Planck compatibility is not rare, but future experiments may rule out this class of models, iv) The smoother the potentials, the sharper the predictions, v) Hyperparameters can transition from stiff to sloppy, vi) Despite tachyons, isocurvature can decay, vii) Eigenvalue repulsion drives the predictions. We conclude that many of the `generic predictions' of single-field inflation can be emergent features of complex inflation models.

We study a new procedure to measure the sound horizon scale via Baryonic Acoustic Oscillations (BAO). Instead of fitting the measured power spectrum (PS) to a theoretical model containing the cosmological informations and all the nonlinear effects, we define a procedure to project out (or to "extract") the oscillating component from a given nonlinear PS. We show that the BAO scale extracted in this way is extremely robust and, moreover, can be reproduced by simple theoretical models at any redshift. By using N-body simulations, we discuss the effect of the nonlinear evolution of the matter field, of redshift space distortions and of scale-dependent halo bias, showing that all these effects can be reproduced with sub-percent accuracy. We give a one-parameter theoretical model based on a simple (IR) modification of 1-loop perturbation theory, which reproduces the BAO scale from measurements of halo clustering in redshift space at better than 0.1% level and does not need any external UV input, such as coefficients measured from N-body simulations.

We develop a method of reconstructing the lensing field from lensed CMB temperature and polarization maps in real space as an alternative to the harmonic space estimators currently in use by extending an existing real space lensing estimator for temperature to polarization. Real space estimators have the advantage of being local in nature and they are thus equipped to deal with the nonuniform sky coverage, especially galactic cuts and point source excisions, found in experimental data. We characterize some of the properties and limitations of these estimators and test them on simulated maps with Planck, AdvACT and CMB-S4 noise. We show that the reconstructions for large-scale lensing fields are accurate, and that the polarization reconstructions improve on those from CMB temperature maps for future experiments as expected. High-fidelity lensing maps can be reconstructed with futuristic experiments like CMB-S4.

Exact solutions describing rotating black holes can offer important tests for alternative theories of gravity, motivated by the dark energy and dark matter problems. We present an analytic rotating black hole solution for a class of vector-tensor theories of modified gravity, valid for arbitrary values of the rotation parameter. The new configuration is characterised by parametrically large deviations from the Kerr-Newman geometry, controlled by non-minimal couplings between vectors and gravity. It has an oblate horizon in Boyer-Lindquist coordinates, and it can rotate more rapidly and have a larger ergosphere than black holes in General Relativity (GR) with the same asymptotic properties. We analytically investigate the features of the innermost stable circular orbits for massive objects on the equatorial plane, and show that stable orbits lie further away from the black hole horizon with respect to rotating black holes in GR. We also comment on possible applications of our findings for the extraction of rotational energy from the black hole.

It is shown that the vacuum-like energy of the Higgs potential at non-zero temperatures leads, in the course of the cosmological expansion, to a small but non-negligible rise of the entropy density in the comoving volume. This increase is calculated in the frameworks of the minimal standard model. The result can have a noticeable effect on the outcome of baryo-through-leptogenesis.

We investigate inflation models in Jordan frame supergravity, in which an inflaton non-minimally couples to the scalar curvature. By imposing the condition that an inflaton would have the canonical kinetic term in the Jordan frame, we construct inflation models with asymptotically flat potential through pole inflation technique and discuss their relation to the models based on Einstein frame supergravity. We also show that the model proposed by Ferrara et al. has special position and the relation between the Kähler potential and the frame function is uniquely determined by requiring that scalars take the canonical kinetic terms in the Jordan frame and that a frame function consists only of a holomorphic term (and its anti-holomorphic counterpart) for symmetry breaking terms. Our case corresponds to relaxing the latter condition.

Matter bounces refer to scenarios wherein the universe contracts at early times as in a matter dominated epoch until the scale factor reaches a minimum, after which it starts expanding. While such scenarios are known to lead to scale invariant spectra of primordial perturbations after the bounce, the challenge has been to construct completely symmetric bounces that lead to a tensor-to-scalar ratio which is small enough to be consistent with the recent cosmological data. In this work, we construct a model involving two scalar fields (a canonical field and a non-canonical ghost field) to drive the symmetric matter bounce and study the evolution of the scalar perturbations in the model. We find that the model can be completely described in terms of a single parameter, viz. the ratio of the scale associated with the bounce to the value of the scale factor at the bounce. We evolve the scalar perturbations numerically across the bounce and evaluate the scalar power spectra after the bounce. We show that, while the scalar and tensor perturbation spectra are scale invariant over scales of cosmological interest, the tensor-to-scalar ratio proves to be much smaller than the current upper bound from the observations of the cosmic microwave background anisotropies by the Planck mission. We also support our numerical analysis with analytical arguments.

We consider minimal U(1) extensions of the Standard Model in which one of the right-handed neutrinos is charged under the new gauge symmetry and plays the role of dark matter. In particular, we perform a detailed phenomenological study for the case of a U(1)_(B−L)3 flavoured B−L symmetry. If perturbativity is required up to high-scales, we find an upper bound on the dark matter mass of m_χ≲2 TeV, significantly stronger than that obtained in simplified models. Furthermore, if the U(1)_(B−L)3 breaking scalar has significant mixing with the SM Higgs, there are already strong constraints from direct detection. On the other hand, there remains significant viable parameter space in the case of small mixing, which may be probed in the future via LHC Z^' searches and indirect detection. We also comment on more general anomaly-free symmetries consistent with a TeV-scale RH neutrino dark matter candidate, and show that if two heavy RH neutrinos for leptogenesis are also required, one is naturally led to a single-parameter class of U(1) symmetries.

We work out the one-loop and order κ² m_ϕ² UV divergent contributions, coming from Unimodular Gravity and General Relativity, to the S matrix element of the scattering process ϕ + ϕ→ ϕ + ϕ in a λ ϕ⁴ theory with mass m_ϕ. We show that both Unimodular Gravity and General Relativity give rise to the same UV divergent contributions in Dimensional Regularization. This seems to be at odds with the known result that in a multiplicative MS dimensional regularization scheme the General Relativity corrections, in the de Donder gauge, to the beta function, β_λ, of the λ coupling do not vanish, whereas the Unimodular Gravity corrections, in a certain gauge, do vanish. Actually, by comparing the UV divergent contributions calculated in this paper with those which give rise to the non-vanishing gravitational corrections to β_λ, one readily concludes that the UV divergent contributions that yield the just mentioned non-vanishing gravitational corrections to β_λ do not contribute to the UV divergent behaviour of the S matrix element of ϕ + ϕ→ ϕ + ϕ. This shows that any physical consequence—such as the existence of asymptotic freedom due to gravitational interactions—drawn from the value of β_λ is not physically meaningful.

We propose the disappearance of "the hyperon puzzle" in neutron star (NS) by invoking two new-physics prescriptions: modified gravity theory and braneworld scenario. By assuming that NS lives on a 3-brane within a 5d empty AdS bulk, gravitationally governed by Eddington-inspired Born-Infeld (EiBI) theory, the field equations can be effectively cast into the usual Einstein's with "apparent" anisotropic energy-momentum tensor. Solving the corresponding brane-TOV equations numerically, we study its mass-radius relation. It is known that the appearance of finite brane tension λ reduces the compactness of the star. The compatibility of the braneworld results with observational constraints of NS mass and radius can be restored in our model by varying the EiBI's coupling constant, κ. We found that within the astrophysically-accepted range of parameters (0<κ<6×10⁶m² and λ≫1 MeV⁴) the NS can have mass ∼2.1 M_⊙ and radius ∼10 km.

Starting from the static, spherically symmetric black hole solutions in massless Einstein-Maxwell-dilaton (EMD) theories, we build a "skeleton" action, that is, we phenomenologically replace black holes by an appropriate effective point particle action, which is well suited to the formal treatment of the many-body problem in EMD theories. We find that, depending crucially on the value of their scalar cosmological environment, black holes can undergo steep "scalarization" transitions, inducing large deviations to the general relativistic two-body dynamics, as shown, for example, when computing the first post-Keplerian Lagrangian of EMD theories.

There are indications that γ-ray dark objects such as supernovae (SNe) with choked jets, and the cores of active galactic nuclei may contribute to the diffuse flux of astrophysical neutrinos measured by the IceCube observatory. In particular, stripped-envelope SNe have received much attention since they are capable of producing relativistic jets and could explain the diversity in observations of collapsar explosions (e.g., gamma-ray bursts (GRBs), low-luminosity GRBs, and Type Ibc SNe). We use an unbinned maximum likelihood method to search for spatial and temporal coincidences between Type Ibc core-collapse SNe, which may harbor a choked jet, and muon neutrinos from a sample of IceCube up-going track-like events measured from May 2011–May 2012. In this stacking analysis, we find no significant deviation from a background-only hypothesis using one year of data, and are able to place upper limits on the total amount of isotropic equivalent energy that choked jet core-collapse SNe deposit in cosmic rays _cr and the fraction of core-collapse SNe which have a jet pointed towards Earth f_jet. This analysis can be extended with yet to be made public IceCube data, and the increased amount of optically detected core-collapse SNe discovered by wide field-of-view surveys such as the Palomar Transient Factory and All-Sky Automated Survey for Supernovae. The choked jet SNe/high-energy cosmic neutrino connection can be more tightly constrained in the near future.

We study constant-roll inflation using the β-function formalism. We show that the constant rate of the inflaton roll is translated into a first order differential equation for the β-function which can be solved easily. The solutions to this equation correspond to the usual constant-roll models. We then construct, by perturbing these exact solutions, more general classes of models that satisfy the constant-roll equation asymptotically. In the case of an asymptotic power law solution, these corrections naturally provide an end to the inflationary phase. Interestingly, while from a theoretical point of view (in particular in terms of the holographic interpretation) these models are intrinsically different from standard slow-roll inflation, they may have phenomenological predictions in good agreement with present cosmological data.

We extend the transport framework for numerically evaluating the power spectrum and bispectrum in multi-field inflation to the case of a curved field-space metric. This method naturally accounts for all sub- and super-horizon tree level effects, including those induced by the curvature of the field-space. We present an open source implementation of our equations in an extension of the publicly available PyTransport code. Finally we illustrate how our technique is applied to examples of inflationary models with a non-trivial field-space metric.

We analyze the ability of galaxy and CMB lensing surveys to constrain massive neutrinos and new models of dark radiation. We present a Fisher forecast analysis for neutrino mass constraints with the LSST galaxy survey and the CMB S4 survey. A joint analysis of the three galaxy and shear 2-point functions from LSST, along with key systematics parameters and Planck priors, can constrain the neutrino masses to ∑ m_ν = 0.041 eV at 1-σ level, comparable to constraints expected from Stage 4 CMB lensing. If low redshift information from upcoming spectroscopic surveys like DESI is included, the constraint becomes ∑ m_ν = 0.032 eV . These constraints are derived having marginalized over the number of relativistic species (N_eff), which is somewhat degenerate with the neutrino mass. We also explore the gain by combining LSST and CMB S4, that is, using the five relevant auto- and cross-correlations of the two datasets. We conclude that advances in modeling the nonlinear regime and the measurements of other parameters are required to ensure a neutrino mass detection. Using the same datasets, we explore the ability of LSST-era surveys to test "nonstandard" models with dark radiation. We find that if evidence for dark radiation is found from N_eff measurements, the mass of the dark radiation candidate can be measured at a 1-σ level of 0.162 eV for fermionic dark radiation, and 0.137 eV for bosonic dark radiation, for Δ N_eff = 0.15. We also find that the NNaturalness model of Arkani-Hamed et al. [1], with extra light degrees of freedom, has a sub-percent effect on the power spectrum: even more ambitious surveys than the ones considered here will be needed to test such models.

In this paper we present a calculation of the expected flux of the mono-energetic 14.4 keV solar axions emitted by the M1 type nuclear transition of ⁵⁷Fe in the Sun. These axions can be detected, e.g., by inverse coherent Bragg-Primakoff conversion in single-crystal TeO₂ bolometers. The ingredients of this calculation are i) the axion nucleon coupling, estimated in several popular axion models and ii)the nuclear spin matrix elements involving realistic shell model calculations with both proton and neutron excitations. For the benefit of the experiments we have also calculated the branching ratio involving axion and photon emission. We find the solar axion flux on Earth to be Φ_a = 0.703×10⁹cm^-2s^-1 (10⁷ GeV/f_a)² and the branching ratio of axion to photon for the same model to be: w_a/w_γ = 0.229×10^-15 ≈ 2×10^-16.

We study cosmological perturbations in mimetic gravity in the presence of classified higher derivative terms which can make the mimetic perturbations stable. We show that the quadratic higher derivative terms which are independent of curvature and the cubic higher derivative terms which come from curvature corrections are sufficient to remove instabilities in mimetic perturbations. The classified higher derivative terms have the same dimensions but they contribute differently in the background and perturbed equations. Therefore, we can control both the background and the perturbation equations allowing us to construct the higher derivative extension of mimetic dark matter and the mimetic nonsingular bouncing scenarios. The latter can be thought as a new higher derivative effective action for the loop quantum cosmology scenario in which the equations of motion coincide with those suggested by loop quantum cosmology. We investigate a possible connection between the mimetic cosmology and the Randall-Sundrum cosmology.

The spectral index of scalar perturbations is an important observable that allows us to learn about inflationary physics. In particular, a detection of a significant deviation from a constant spectral index could enable us to rule out the simplest class of inflation models. We investigate whether future observations could rule out canonical single-field slow-roll inflation given the parameters allowed by current observational constraints. We find that future measurements of a constant running (or running of the running) of the spectral index over currently available scales are unlikely to achieve this. However, there remains a large region of parameter space (especially when considering the running of the running) for falsifying the assumed class of slow-roll models if future observations accurately constrain a much wider range of scales.

Dark energy equation of state can be effectively described by that of a barotropic fluid. The barotropic fluid model describes the background evolution and the functional form of the equation of state parameter is well constrained by the observations. Equally viable explanations of dark energy are via scalar field models, both canonical and non-canonical; these scalar field models being low energy descriptions of an underlying high energy theory. In this paper, we attempt to reconcile the two approaches to dark energy by way of reconstructing the evolution of the scalar field potential. For this analysis, we consider canonical quintessence scalar field and the phantom field for this reconstruction. We attempt to understand the analytical or semi-analytical forms of scalar field potentials corresponding to typical well behaved parameterisations of dark energy using the constraints from recent observations.

We have developed a phenomenological effective quantum-field theoretical model describing the "hadron gas" of the lightest pseudoscalar mesons, scalar σ-meson and σ-vacuum, i.e. the expectation value of the σ-field, at finite temperatures. The corresponding thermodynamic approach was formulated in terms of the generating functional derived from the effective Lagrangian providing the basic thermodynamic information about the "meson plasma + QCD condensate" system. This formalism enables us to study the QCD transition from the hadron phase with direct implications for cosmological evolution. Using the hypothesis about a positively-definite QCD vacuum contribution stochastically produced in early universe, we show that the universe could undergo a series of oscillations during the QCD epoch before resuming unbounded expansion.

The nonlinear interaction between one graviton and two scalars is enhanced in specific inflationary models, potentially leading to distinguishable signatures in the bispectrum of the cosmic microwave background (CMB) anisotropies. We develop the tools to examine such bispectrum signatures, and show a first application using WMAP temperature data. We consider several ℓ-ranges, estimating the g_tss amplitude parameter, by means of the so-called separable modal methodology. We do not find any evidence of a tensor-scalar-scalar signal at any scale. Our tightest bound on the size of the tensor-scalar-scalar correlator is derived from our measurement including all the multipoles in the range  2 ⩽ ℓ ⩽ 500 and it reads g_tss = −48 ± 28 (68%CL). This is the first direct observational constraint on the primordial tensor-scalar-scalar correlation, and it will be cross-checked and improved by applying the same pipeline to high-resolution temperature and polarization data from Planck and forthcoming CMB experiments.

We study a fermionic dark matter model in which the interaction of the dark and visible sectors is mediated by Higgs portal type couplings. Specifically, we consider the mixing of a dark sector scalar with the scalars of a Two Higgs Doublet Model extension of the Standard Model. Given that scalar exchange will result in a spin-independent dark matter-nucleon scattering cross section, such a model is potentially subject to stringent direct detection constraints. Moreover, the addition of new charged scalars introduce non-trivial flavour constraints. Nonetheless, this model allows more freedom than a standard Higgs portal scenario involving a single Higgs doublet, and much of the interesting parameter space is not well approximated by a Simplified Model with a single scalar mediator. We perform a detailed parameter scan to determine the mass and coupling parameters which satisfy direct detection, flavour, precision electroweak, stability, and perturbativity constraints, while still producing the correct relic density through thermal freezeout.

We study ghost-free multimetric theories for (N+1) tensor fields with a coupling to matter and maximal global symmetry group S_N×(Z₂)^N. Their mass spectra contain a massless mode, the graviton, and N massive spin-2 modes. One of the massive modes is distinct by being the heaviest, the remaining (N−1) massive modes are simply identical copies of each other. All relevant physics can therefore be understood from the case N=2. Focussing on this case, we compute the full perturbative action up to cubic order and derive several features that hold to all orders in perturbation theory. The lighter massive mode does not couple to matter and neither of the massive modes decay into massless gravitons. We propose the lighter massive particle as a candidate for dark matter and investigate its phenomenology in the parameter region where the matter coupling is dominated by the massless graviton. The relic density of massive spin-2 can originate from a freeze-in mechanism or from gravitational particle production, giving rise to two different dark matter scenarios. The allowed parameter regions are very different from those in scenarios with only one massive spin-2 field and more accessible to experiments.

In the concordance model of the Universe, the matter distribution—as observed in galaxy number counts or the intensity of line emission (such as the 21cm line of neutral hydrogen) —should have a kinematic dipole due to the Sun's motion relative to the CMB rest-frame. This dipole should be aligned with the kinematic dipole in the CMB temperature. Accurate measurement of the direction of the matter dipole will become possible with future galaxy surveys, and this will be a critical test of the foundations of the concordance model. The amplitude of the matter dipole is also a potential cosmological probe. We derive formulas for the amplitude of the kinematic dipole in galaxy redshift and intensity mapping surveys, taking into account the Doppler, aberration and other relativistic effects. The amplitude of the matter dipole can be significantly larger than that of the CMB dipole. Its redshift dependence encodes information on the evolution of the Universe and on the tracers, and we discuss possible ways to determine the amplitude.

We investigate the three-point correlation between the Lyman-α forest and the CMB weak lensing (δ_F δ_F κ) expressed as the cross-correlation between the CMB weak lensing field and local variations in the forest power spectrum. In addition to the standard gravitational bispectrum term, we note the existence of a non-standard systematic term coming from mis-estimation of the mean flux over the finite length of Lyman-α skewers. We numerically calculate the angular cross-power spectrum and discuss its features. We integrate it into zero-lag correlation function and compare our predictions with recent results by Doux et al. We find that our predictions are statistically consistent with the measurement, and including the systematic term improves the agreement with the measurement. We comment on the implication of the response of the Lyman-α forest power spectrum to the long-wavelength density perturbations.

We consider a modified gravity model with a massive graviton, but which nevertheless only propagates two gravitational degrees of freedom and which is free of ghosts. We show that non-singular bouncing cosmological background solutions can be generated. In addition, the mass term for the graviton prevents anisotropies from blowing up in the contracting phase and also suppresses the spectrum of gravitational waves compared to that of the scalar cosmological perturbations. This addresses two of the main problems of the matter bounce scenario.

The presence of additional particles during inflation leads to non-Gaussianity in late-time correlators of primordial curvature perturbations. The shape and amplitude of this signal depend on the mass and spin of the extra particles. Constraints on this distinct form of primordial non-Gaussianity, therefore, provide a wealth of information on the particle content during inflation. We investigate the potential of upcoming galaxy surveys in constraining such a signature through its impact on the observed galaxy power spectrum. Primordial non-Gaussianity of various shapes induces a scale-dependent bias on tracers of large-scale structure, such as galaxies. Using this signature we obtain constraints on massive particles during inflation, which can have non-zero spins. In particular, we show that the prospects for constraining particles with spins 0 and 1 are promising, while constraining particles with spin 2 from power spectrum alone seems challenging. We show that the multi-tracer technique can significantly improve the constraints from the power spectrum by at least an order of magnitude. Furthermore, we analyze the effect of non-linearities due to gravitational evolution on the forecasted constraints on the masses of the extra particles and the amplitudes of the imprinted non-Gaussian signal. We find that gravitational evolution affects the constraints by less than a factor of 2.

We perform a quantitative comparison between N-body simulations and the Schrö\-dinger-Poisson system in 1+1 dimensions. In particular, we study halo formation with different initial conditions. We observe the convergence of various observables in the Planck constant ℏ and also test virialization. We discuss the generation of higher order cumulants of the particle distribution function which demonstrates that the Schrödinger-Poisson equations should not be perceived as a generalization of the dust model with quantum pressure but rather as one way of sampling the phase space of the Vlasov-Poisson system—just as N-body simulations. Finally, we quantitatively recover the scaling behavior of the halo density profile from N-body simulations.

We consider a recently proposed model in which dark matter interacts with a thermal background of dark radiation. Dark radiation consists of relativistic degrees of freedom which allow larger values of the expansion rate of the universe today to be consistent with CMB data (H₀-problem). Scattering between dark matter and radiation suppresses the matter power spectrum at small scales and can explain the apparent discrepancies between ΛCDM predictions of the matter power spectrum and direct measurements of Large Scale Structure LSS (σ₈-problem). We go beyond previous work in two ways: 1. we enlarge the parameter space of our previous model and allow for an arbitrary fraction of the dark matter to be interacting and 2. we update the data sets used in our fits, most importantly we include LSS data with full k-dependence to explore the sensitivity of current data to the shape of the matter power spectrum. We find that LSS data prefer models with overall suppressed matter clustering due to dark matter - dark radiation interactions over ΛCDM at 3–4 σ. However recent weak lensing measurements of the power spectrum are not yet precise enough to clearly distinguish two limits of the model with different predicted shapes for the linear matter power spectrum. In two appendices we give a derivation of the coupled dark matter and dark radiation perturbation equations from the Boltzmann equation in order to clarify a confusion in the recent literature, and we derive analytic approximations to the solutions of the perturbation equations in the two physically interesting limits of all dark matter weakly interacting or a small fraction of dark matter strongly interacting.

CMB temperature fluctuation observations provide a precise measurement of the primordial power spectrum on large scales, corresponding to wavenumbers 10⁻³ Mpc⁻¹ ≲ k ≲ 0.1  Mpc⁻¹, [1-7, 11]. Luminous red galaxies and galaxy clusters probe the matter power spectrum on overlapping scales (0.02 Mpc⁻¹ ≲ k ≲ 0.7 Mpc⁻¹; [10, 12-20]), while the Lyman-alpha forest reaches slightly smaller scales (0.3 Mpc⁻¹ ≲ k ≲ 3 Mpc⁻¹; [22]). These observations indicate that the primordial power spectrum is nearly scale-invariant with an amplitude close to 2 × 10⁻⁹, [5, 23-28]. These observations strongly support Inflation and motivate us to obtain observations and constraints reaching to smaller scales on the primordial curvature power spectrum and by implication on Inflation. We are able to obtain limits to much higher values of k ≲ 10⁵ Mpc⁻¹ and with less sensitivity even higher k ≲ 10¹⁹− 10²³ Mpc⁻¹ using limits from CMB spectral distortions and other limits on ultracompact minihalo objects (UCMHs) and Primordial Black Holes (PBHs). PBHs are one of the known candidates for the Dark Matter (DM). Due to their very early formation, they could give us valuable information about the primordial curvature perturbations. These are complementary to other cosmological bounds on the amplitude of the primordial fluctuations. In this paper, we revisit and collect all the published constraints on both PBHs and UCMHs. We show that unless one uses the CMB spectral distortion, PBHs give us a very relaxed bounds on the primordial curvature perturbations. UCMHs, on the other hand, are very informative over a reasonable k range (3 ≲ k ≲ 10⁶ Mpc⁻¹) and lead to significant upper-bounds on the curvature spectrum. We review the conditions under which the tighter constraints on the UCMHs could imply extremely strong bounds on the fraction of DM that could be PBHs in reasonable models. Failure to satisfy these conditions would lead to over production of the UCMHs which is inconsistent with the observations. Therefore, we can almost rule out PBH within their overlap scales with the UCMHs. We compare the UCMH bounds coming from those experiments which are sensitive to the nature of the DM, such as γ-rays, Neutrinos and Reionization, with those which are insensitive to the type of the DM, e.g. the pulsar-timing as well as CMB spectral distortion. We explicitly show that they lead to comparable results which are independent of the type of DM. These bounds however do depend on the required initial density perturbation, i.e. δ_min. It could be either a constant or a scale-dependent function. As we will show, the constraints differ by three orders of magnitude depend on our choice of required initial perturbations.

In this paper we study the vacuum quantum fluctuations of the stationary modes of an uncharged scalar field with mass m around a Schwarzschild black hole with mass M, at zero and non-zero temperatures. The procedure consists of calculating the energy eigenvalues starting from the exact solutions found for the dynamics of the scalar field, considering a frequency cutoff in which the particle is not absorbed by the black hole. From this result, we obtain the exterior contributions for the vacuum energy associated to the stationary states of the scalar field, by considering the half-summing of the levels of energy and taking into account the respective degeneracies, in order to better capture the nontrivial topology of the black hole spacetime. Then we use the Riemann's zeta function to regularize the vacuum energy thus found. Such a regularized quantity is the Casimir energy, whose analytic computation we show to yield a convergent series. The Casimir energy obtained does not take into account any boundaries artificially imposed on the system, just the nontrivial spacetime topology associated to the source and its singularity. We suggest that this latter manifests itself through the vacuum tension calculated on the event horizon. We also investigate the problem by considering the thermal corrections via Helmholtz free energy calculation, computing the Casimir internal energy, the corresponding tension on the event horizon, the Casimir entropy, and the thermal capacity of the regularized quantum vacuum, analyzing their behavior at low and high temperatures, pointing out the thermodynamic instability of the system in the considered regime, i.e. mM≪ 1.

Within the framework of braneworlds, several interesting physical effects can be described in a wide range of energy scales, starting from high-energy physics to cosmology and low-energy physics. An usual way to generate a thick braneworld model relies in coupling a bulk scalar field to higher dimensional warped gravity. Quite recently, a novel braneworld was generated with the aid of a tachyonic bulk scalar field, having several remarkable properties. It comprises a regular and stable solution that contains a relevant 3-brane with de Sitter induced metric, arising as an exact solution to the 5D field equations, describing the inflationary eras of our Universe. Besides, it is asymptotically flat, despite of the presence of a negative 5D cosmological constant, which is an interesting feature that contrasts with most of the known, asymptotically either dS or AdS models. Moreover, it encompasses a graviton spectrum with a single massless bound state, accounting for 4D gravity localized on the brane, separated from the continuum of Kaluza-Klein massive graviton modes by a mass gap that makes the 5D corrections to Newton's law to decay exponentially. Finally, gauge, scalar and fermion fields are also shown to be localized on this braneworld. In this work, we show that this tachyonic braneworld allows for a nontrivial solution with a vanishing 5D cosmological constant that preserves all the above mentioned remarkable properties with a less amount of parameters, constituting an important contribution to the construction of a realistic cosmological braneworld model.

The model in which Primordial Black Holes (PBHs) constitute a non-negligible fraction of the dark matter has (re)gained popularity after the first detections of binary black hole mergers. Most of the observational constraints to date have been derived assuming a single mass for all the PBHs, although some more recent works tried to generalize constraints to the case of extended mass functions. Here we derive a general methodology to obtain constraints for any PBH Extended Mass Distribution (EMD) and any observables in the desired mass range. Starting from those obtained for a monochromatic distribution, we convert them into constraints for EMDs by using an equivalent, effective mass M_eq that depends on the specific observable. We highlight how limits of validity of the PBH modelling affect the EMD parameter space. Finally, we present converted constraints on the total abundance of PBH from microlensing, stellar distribution in ultra-faint dwarf galaxies and CMB accretion for Lognormal and Power Law mass distributions, finding that EMD constraints are generally stronger than monochromatic ones.

We derive and numerically implement an algorithm for estimating the 3D power spectrum of the Lyman-α (Lyα) forest flux fluctuations. The algorithm exploits the unique geometry of Lyα forest data to efficiently measure the cross-spectrum between lines of sight as a function of parallel wavenumber, transverse separation and redshift. We start by approximating the global covariance matrix as block-diagonal, where only pixels from the same spectrum are correlated. We then compute the eigenvectors of the derivative of the signal covariance with respect to cross-spectrum parameters, and project the inverse-covariance-weighted spectra onto them. This acts much like a radial Fourier transform over redshift windows. The resulting cross-spectrum inference is then converted into our final product, an approximation of the likelihood for the 3D power spectrum expressed as second order Taylor expansion around a fiducial model. We demonstrate the accuracy and scalability of the algorithm and comment on possible extensions. Our algorithm will allow efficient analysis of the upcoming Dark Energy Spectroscopic Instrument dataset.

We study the problem of initial conditions for slow-roll inflation along a plateau-like scalar potential within the framework of fluctuation-dissipation dynamics. We consider, in particular, that inflation was preceded by a radiation-dominated epoch where the inflaton is coupled to light degrees of freedom and may reach a near-equilibrium state. We show that the homogeneous field component can be sufficiently localized at the origin to trigger a period of slow-roll if the interactions between the inflaton and the thermal degrees of freedom are sufficiently strong and argue that this does not necessarily spoil the flatness of the potential at the quantum level. We further conclude that the inflaton can still be held at the origin after its potential begins to dominate the energy balance, leading to a period of thermal inflation. This then suppresses the effects of nonlinear interactions between the homogeneous and inhomogeneous field modes that could prevent the former from entering a slow-roll regime. Finally, we discuss the possibility of an early period of chaotic inflation, at large field values, followed by a first stage of reheating and subsequently by a second inflationary epoch along the plateau about the origin. This scenario could prevent an early overclosure of the Universe, at the same time yielding a low tensor-to-scalar ratio in agreement with observations.

The Daya Bay Experiment consists of eight identically designed detectors located in three underground experimental halls named as EH1, EH2, EH3, with 250, 265 and 860 meters of water equivalent vertical overburden, respectively. Cosmic muon events have been recorded over a two-year period. The underground muon rate is observed to be positively correlated with the effective atmospheric temperature and to follow a seasonal modulation pattern. The correlation coefficient α, describing how a variation in the muon rate relates to a variation in the effective atmospheric temperature, is found to be α_EH1 = 0.362±0.031, α_EH2 = 0.433±0.038 and α_EH3 = 0.641±0.057 for each experimental hall.