Search

We investigate the possibility of reducing the number of degrees of freedom (d.o.f.) starting from generic metric theories of gravity by introducing multiple auxiliary constraints (ACs), under the restriction of retaining spatial covariance as a gauge symmetry. Arbitrary numbers of scalar-, vector- and tensor-type ACs are considered a priori, yet we find that no vector- and tensor-type constraints should be introduced, and that scalar-type ACs should be no more than four for the purpose of constructing minimally modified gravity (MMG) theories which propagate only two tensorial d.o.f., like general relativity (GR). Through a detailed Hamiltonian analysis, we exhaust all the possible classifications of ACs and find out the corresponding minimalizing and symmetrizing conditions for obtaining the MMG theories. In particular, no condition is required in the case of four ACs, hence in this case the theory can couple with matter consistently and naturally. To illustrate our formalism, we build a concrete model for this specific case by using the Cayley-Hamilton theorem and derive the dispersion relation of the gravitational waves, which is subject to constraints from the observations., 15 pages, 1 table, no figure

We investigate the properties of dark energy halos in models with a nonminimal coupling in the dark sector. We show, using a quasistatic approximation, that a coupling of the mass of dark matter particles to a standard quintessence scalar field $\phi$ generally leads to the formation of dark energy concentrations in and around compact dark matter objects. These are associated with regions where scalar field gradients are large and the dark energy equation of state parameter is close to $-1/3$. We find that the energy and radius of a dark energy halo are approximately given by $E_{\rm halo} \sim \boldsymbol{\beta}^2 \varphi \, m$ and $r_{\rm halo} \sim \sqrt{\boldsymbol{\beta} \,\varphi ({R}/{H})}$, where $\varphi=Gm/(R c^2)$, $m$ and $R$ are, respectively, the mass and radius of the associated dark matter object, $\boldsymbol{\beta} = -(8\pi G)^{-1/2} d \ln m/d \phi$ is the nonminimal coupling strength parameter, $H$ is the Hubble parameter, $G$ is the gravitational constant, and $c$ is the speed of light in vacuum. We further show that current observational limits on $\boldsymbol{\beta}$ over a wide redshift range lead to stringent constraints on $E_{\rm halo}/m$ and, therefore, on the impact of dark energy halos on the value of the dark energy equation of state parameter. We also briefly comment on potential backreaction effects that may be associated with the breakdown of the quasistatic approximation and determine the regions of parameter space where such a breakdown might be expected to occur., Comment: 6 pages [minor correction in Eq. (28)]

The detection of black-hole binaries at high redshifts is a cornerstone of the science case of third-generation gravitational-wave interferometers. The star-formation rate peaks at z~2 and decreases by orders of magnitude by z~10. Any confident detection of gravitational waves from such high redshifts would imply either the presence of stars formed from pristine material originating from cosmological nucleosynthesis (the so-called population III stars), or black holes that are the direct relics of quantum fluctuations in the early Universe (the so-called primordial black holes). Crucially, detecting sources at cosmological distances does not imply inferring that sources are located there, with the latter posing more stringent requirements. To this end, we present two figures of merit, which we refer to as "z-z plot" and "inference horizon", that quantify the largest redshift one can possibly claim a source to be beyond. We argue that such inference requirements, in addition to detection requirements, should be investigated when quantifying the scientific payoff of future gravitational-wave facilities., 6 pages, 4 figures. v2: matches version published on PRD Lett

Axions are considered as a candidate of dark matter. The axions form coherent clouds which delay and amplify gravitational waves (GWs) at a resonant frequency. That is, the coherent axion clouds produce secondary GWs following a primary wave. All GWs from compact binary mergers detected so far propagate in the Milky Way halo composed of dark matter so that such GWs are followed by the secondary GWs. Since the properties of the secondary GWs depend on the axion mass and coupling with the parity violating sector of gravity, we can search for a characteristic signal produced by axions. In our previous study, we have developed a search method optimized for the axion signals, and obtained about ten times stronger constraint on the coupling than the previous best one for the axion mass range, $[1.7 \times 10^{-13},\,8.5 \times 10^{-12}]\,\mathrm{eV}$. However, in our previous search, we assume that dark matter is homogeneously distributed in the Milky Way halo, which is unrealistic. In this paper, we extend the dark matter profile to a more realistic one called the core NFW profile, and find that the constraint in our previous study is robust to the dark matter profile., Comment: 12 pages, 8 figures. Minor English corrections

We study the observational implications of a class of inflationary models wherein the inflaton is coupled to the Einstein tensor through a generalized non-minimal derivative coupling (GNMDC). In particular, we explore whether these models can generate suitable features in the primordial spectrum of curvature perturbations as a possible explanation for the large-scale anomalies associated with the angular power spectrum of CMB temperature anisotropies. We derive model-independent constraints on the GNMDC function for such a scenario, considering both the scalar and tensor perturbations. We modify CosmoMC to accommodate our GNMDC framework and investigate different classes of inflationary models using a fully consistent numerical approach. We find that the hilltop-quartic model with a specific choice of the GNMDC function provides a considerable improvement over the best-fit reference $\Lambda$CDM model with a nearly scale-invariant power spectrum. While the large-scale structure observations should be able to provide independent constraints, future CMB experiments, such as CMB-S4 and CMB-Bharat, are expected to constrain further the parameter space of such beyond canonical single-field inflationary models., Comment: v1: 34 pages, 7 figures, 4 tables; v2: 26 pages, 8 figures, 4 tables, published in Phys. Rev. D

We consider a model of baryogenesis based on adding lepton number-violating quadratic mass terms to the inflaton potential of a non-minimally coupled inflation model. The $L$-violating mass terms generate a lepton asymmetry in a complex inflaton field via the mass term Affleck-Dine mechanism, which is transferred to the Standard Model (SM) sector when the inflaton decays to right-handed (RH) neutrinos. The model is minimal in that it requires only the SM sector, RH neutrinos, and a non-minimally coupled inflaton sector. We find that baryon isocurvature fluctuations can be observable in metric inflation but are negligible in Palatini inflation. The model is compatible with reheating temperatures that may be detectable in the observable primordial gravitational waves predicted by metric inflation., 18 pages, 6 figures. Published version

Gravitons radiated from light, evaporating black holes contribute to the stochastic background of gravitational waves. The spectrum of such emission depends on both the mass and the spin of the black holes, as well as on the redshifting that occurs between the black hole formation and today. This, in turn, depends on the expansion history of the universe, which is largely unknown and unconstrained at times prior to the synthesis of light elements. Here, we explore the features of the stochastic background of gravitational waves from black hole evaporation under a broad range of possible early cosmological histories. The resulting gravitational wave signals typically peak at very high frequencies, and offer opportunities for proposed ultra-high frequency gravitational wave detectors. Lower-frequency peaks are also possible, albeit with a suppressed intensity that is likely well below the threshold of detectability. We find that the largest intensity peaks correspond to cosmologies dominated by fluids with equations of state that have large pressure-to-density ratios. Such scenarios can be constrained on the basis of violation of $\Delta N_{\rm eff}$ bounds., Comment: 13 pages, 7 figures

Bigravity is one of the natural extensions of general relativity and contains an additional massive spin-2 field which can be a good candidate for dark matter. To discuss the production of spin-2 dark matter, we study fixed point solutions of the background equations for axisymmetric Bianchi type-I Universes in two bigravity theories without Boulware-Deser ghost, i.e., Hassan-Rosen bigravity and Minimal Theory of Bigravity. We investigate the local and global stability of the fixed points and classify them. Based on the general analysis, we propose a new scenario where spin-2 dark matter is produced by the transition from an anisotropic fixed point solution to isotropic one. The produced spin-2 dark matter can account for all or a part of dark matter and can be directly detected by laser interferometers in the same way as gravitational waves., Comment: 15 pages, 3 figures

We explore the dynamics of FLRW cosmologies which consist of dark matter, radiation and dark energy with a quadratic equation of state. Standard cosmological singularities arise due to energy conditions which are violated by dark energy, therefore we focus our analysis on non-singular bouncing and cyclic cosmologies, in particular focusing on the possibility of closed models always having a bounce for any initial conditions. We analyse the range of dynamical behaviour admitted by the system, and find a class of closed models that admit a non-singular bounce, with early- and late-time accelerated expansion connected by a decelerating phase. In all cases, we find the bouncing models are only relevant when dark matter and radiation appear at a certain energy scale, and so require a period such as reheating. We then investigate imposing an upper bound on the dark matter and radiation, such that their energy densities cannot become infinite. We find that bounces are always the general closed model, and a class of models exist with early- and late-time acceleration, connected by a decelerating phase. We also consider parameter values for the dark energy component, such that the discrepancy between the observed value of $\Lambda$ and the theoretical estimates of the contributions to the effective cosmological constant expected from quantum field theory would be explained. However, we find that the class of models left does not allow for an early- and late-time accelerated expansion, connected by a decelerating period where large-scale structure could form. Nonetheless, our qualitative analysis serves as a basis for the construction of more realistic models with realistic quantitative behaviour., Comment: 26 pages, 21 figures. Updated version to reflect changes made for publication

In this paper, we have investigated the processes of evaporation and accretion of primordial black holes during the radiation-dominated era and the matter-dominated era. This subject is very important since usually these two processes are considered independent of each other. In other words, previous works consider them in such a way that they do not have a direct effect on each other, and as a result, their effects on the mass of primordial black holes are calculated separately. The calculations of this paper indicate that assuming these two processes independently of each other will lead to wrong results that only give correct answers within certain limits. In fact, in general, it is a mistake to consider the static state for the event horizon of primordial black holes and perform calculations related to their evaporation, while the radius of primordial black holes is constantly changing due to accretion. In addition, we have shown that considering the dynamic event horizon in some masses and in some times can lead to the shutdown of the Hawking evaporation process. This study is much more accurate and detailed than our previous study. These calculations show well the mass evolution of primordial black holes from the time of formation to the end of the matter-dominated era, taking into account both the main processes governing black holes, evaporation and accretion.

The linear matter power spectrum $P(k,z)$ connects theory with large scale structure observations in cosmology. Its scale dependence is entirely encoded in the matter transfer function $T(k)$, which can be computed numerically by Boltzmann solvers, and can also be computed semi-analytically by using fitting functions such as the well-known Bardeen-Bond-Kaiser-Szalay (BBKS) and Eisenstein-Hu (EH) formulae. However, both the BBKS and EH formulae have some significant drawbacks. On the one hand, although BBKS is a simple expression, it is only accurate up to $10\%$, which is well above the $1\%$ precision goal of forthcoming surveys. On the other hand, while EH is as accurate as required by upcoming experiments, it is a rather long and complicated expression. Here, we use the Genetic Algorithms (GAs), a particular machine learning technique, to derive simple and accurate fitting formulae for the transfer function $T(k)$. When the effects of massive neutrinos are also considered, our expression slightly improves over the EH formula, while being notably shorter in comparison., Comment: 12 pages, 7 figures, 2 tables. Changes match published version

We present first-order post-Newtonian (1PN) approximations of a general imperfect fluid and of an axion as a coherently oscillating massive scalar field, both in the cosmological context. For the axion, using the Klein transformation and Madelung transformation we derive the Schr\"odinger and Madelung hydrodynamic formulations, respectively, in exact covariant way and to 1PN order. Complete sets of equations for the 1PN formulations are derived without fixing the temporal gauge condition. We study the linear instability in cosmology and a static limit for both fluid and axion; these are presented independently of the gauge condition to 1PN order, thus are naturally gauge-invariant., Comment: 18 pages, no figure

Black hole superradiance provides a window into the dynamics of light scalar fields and their interactions close to a rotating black hole. Due to the rotation of the black hole, the amplitude of the scalar field becomes magnified, leading to a "black hole bomb" effect. Recent work has demonstrated that rotating black holes in dynamical Chern-Simons gravity possess unique structures, the "Chern-Simons caps," which may influence the behavior of matter near the black hole. Motivated by the presence of these caps, we study superradiance in dynamical Chern-Simons gravity in the context of a slowly rotating black hole. We find that additional modes are excited and contribute to the superradiance beyond what is expected for a Kerr black hole. Studying the superradiant spectrum of perturbations, we find that the Chern-Simons contributions give rise to small corrections to the angular dependence of the resulting scalar cloud. Finally, we comment on potential observable consequences and future avenues for investigation., Comment: Error in Eq. 51-53 and 69-71 corrected, updated discussion

We discuss, in the context of alternative theories of gravity with nonminimal coupling between matter and curvature, if inflationary solutions driven by a single scalar field can be reconciled with the swampland conjectures about the emergence of de Sitter solutions in string theory. We find that the slow-roll conditions are incompatible with the swampland conjectures for a fairly generic inflationary solution in such alternative theories of gravity., 12 pages

Future ground-based gravitational wave observatories will be ideal probes of the environments surrounding black holes with masses $1 - 10\,\mathrm{M_\odot}$. Binary black hole mergers with mass ratios of order $q=m_2/m_1\lesssim10^{-3}$ can remain in the frequency band of such detectors for months or years, enabling precision searches for modifications of their gravitational waveforms with respect to vacuum inspirals. As a concrete example of an environmental effect, we consider here a population of binary primordial black holes which are expected to be embedded in dense cold dark matter spikes. We provide a viable formation scenario for these systems compatible with all observational constraints, and predict upper and lower limits on the merger rates of small mass ratio pairs. Given a detected signal of one such system by either Einstein Telescope or Cosmic Explorer, we show that the properties of the binary and of the dark matter spike can be measured to excellent precision with one week's worth of data, if the effect of the dark matter spike on the waveform is taken into account. However, we show that there is a risk of biased parameter inference or missing the events entirely if the effect of the predicted dark matter overdensity around these objects is not properly accounted for., Comment: 14 pages plus appendices, 15 figures

We study the properties of the Newtonian gravitational potential in a spherical Universe for different topologies. For this, we use the non-Euclidean Newtonian theory developed in Vigneron [2022, Class. & Quantum Gravity, 39, 155006] describing Newtonian gravitation in a spherical or hyperbolic Universe. The potential is calculated for a point mass in all the globally homogeneous regular spherical topologies, i.e. whose fundamental domain is unique and is a platonic solid. We provide the exact solution and the Taylor expansion series of the potential at a test position near the point mass. We show that the odd terms of the expansion can be interpreted as coming from the presence of a non-zero spatial scalar curvature, while the even terms relate to the closed nature of the topological space. A consequence is that, compared to the point mass solution in a 3-torus, widely used in Newtonian cosmological simulations, the spherical cases all feature an additional attractive first order term dependent solely on the spatial curvature. The choice of topology only affects the potential at second order and higher. For typical estimates of cosmological scales (curvature and topology), the strongest topological effect occurs in the case of the Poincaré dodecahedral space, but in general the effect of curvature dominates over topology. We also provide the set of equations that can be used to perform $N$-body simulations of structure formation in spherical topologies., Accepted in Physical Review D, 12 pages, 4 tables, calculation source code at https://codeberg.org/boud/topoaccel

In the former part, we study the gravitational collapse of pressureless dust and find special solutions, where, in both the physical and fiducial sectors, the exterior and interior spacetime geometries are given by the Schwarzschild spacetimes and the Friedmann-Lema\^itre-Robertson-Walker universes dominated by pressureless dust, respectively, with specific time slicings. In the case where the Lagrange multipliers are trivial and have no jump across the matter interfaces in both the physical and fiducial sectors, the junction conditions across them remain the same as those in general relativity (GR). For simplicity, we foliate the interior geometry by homogeneous and isotropic spacetimes. For a spatially flat interior universe, we foliate the exterior geometry by a time-independent flat space, while for a spatiallycurved interior universe, we foliate the exterior geometry by a time-independent space with deficit solid angle. Despite the rather restrictive choice of foliations, we find interesting classes of exact solutions that represent gravitational collapse in MTBG. In the spatially flat case, under a certain tuning of the initial condition, we find exact solutions of matter collapse in which the two sectors evolve independently. In the spatially closed case, once the matter energy densities and the Schwarzschild radii are tuned between the two sectors, we find exact solutions that correspond to the Oppenheimer-Snyder model in GR. In the latter part, we study odd-parity perturbations of the Schwarzschild$-$de Sitter solutions written in the spatially flat coordinates. For the higher-multipole modes $\ell\geq2$, we find that, in general, the system reduces to that of four physical modes, where two of them are dynamical and the remaining two are shadowy, i.e., satisfying only elliptic equations., Comment: 28 pages, published version

We put constraints on the peak absolute magnitude, $M_B$ of type Ia supernova using the Pantheon sample for type Ia supernova observations and the cosmic chronometers data for the Hubble parameter by a model independent and non-parametric approach. Our analysis is based on the Gaussian process regression. We find percent level bounds on the peak absolute magnitude given as $M_B=-19.384\pm0.052$. For completeness and to check the consistency of the results, we also include the Baryon acoustic oscillation data and the prior of the comoving sound horizon from Planck 2018 cosmic microwave background observations. The inclusion of these two data gives tighter constraints on $M_B$ at the sub-percent level. We obtain constraints on $M_B$ from the combination of pantheon compilation of type Ia supernova observations and baryon acoustic oscillation observations given as $M_B=-19.396\pm0.016$. When adding the cosmic chronometer observations with these observations, we find $M_B=-19.395\pm0.015$. The mean values of peak absolute magnitude from all these data are consistent with each other and the values are approximately equal to $-19.4$., Accepted for publication in PRD; detailed calculation and discussion of methodology added; appendices added; 13+10 pages; 10+2 figures; 1+6 tables; comments are welcome

A violation of the null energy condition (NEC) during inflation in a single-field inflation model will naturally enhance the amplitude of the parity violation effect (defined by $\Delta\chi$) of inflationary primordial gravitational waves (GWs), provided the inflaton is non-minimally coupled to a gravitational Chern-Simons term. After going through the NEC-violating phase, the Universe enters subsequent slow-roll inflation with a higher energy scale (i.e., a greater Hubble parameter $H$), which results in an enhanced nearly scale-invariant power spectrum (i.e., $P_{\rm T}$) of inflationary primordial GWs in the high-frequency band, while $P_{\rm T}$ remains consistent with observations in the frequency band of the cosmic microwave background. Therefore, the violation of NEC during inflation will amplify the observability (i.e., $P_{\rm T}\cdot\Delta\chi$) of the parity violation effect on small scales. Intriguingly, our model has particular oscillatory features on $\Delta\chi$ that may not be mimicked by others., Comment: 24 pages, 16 figures; new references added, updated to match the published version

Single-photon atom gradiometry is a powerful experimental technique that can be employed to search for the oscillation of atomic transition energies induced by ultralight scalar dark matter (ULDM). In the sub-Hz regime the background is expected to be dominated by gravity gradient noise (GGN), which arises as a result of mass fluctuations around the experiment. In this work we model the GGN as surface Rayleigh waves, and we construct a likelihood-based analysis that consistently folds GGN into the sensitivity estimates of vertical atom gradiometers in the frequency window between 1 mHz and 1 Hz. We show that in certain geological settings GGN can be significantly mitigated when operating a multigradiometer configuration, which consists of three or more atom interferometers in the same baseline. Multigradiometer experiments, such as future versions of AION and MAGIS-100, have the potential to probe regions of scalar ULDM parameter space in the sub-Hz regime that have not been excluded by existing experiments., Comment: 28 pages plus appendices, 11 figures, comments welcome. v3: matches published version

Galilean Genesis is generically plagued with a strong coupling problem, but this can be avoided depending on the hierarchy between a classical energy scale of genesis and a strong coupling scale. In this paper, we investigate whether or not the models of Galilean Genesis without the strong coupling problem can explain the statistical properties of the observed CMB fluctuations based on two unified frameworks of Galilean Genesis. By focusing on the class in which the propagation speeds of the scalar and tensor perturbations are constant, we show that the models avoiding strong coupling and allowing a slightly red-tilted scalar power spectrum suffer from an overproduction of a scalar non-Gaussianity., 16 pages, 1 figure

In this paper, we consider the hypothesis that fractions of dark matter could be constituted by primordial black holes (PBHs). To test this possibility, we work out the observational properties of a static black hole embedded in the dark matter envelope made of a PBH source. The corresponding modifications of geometry due to such a physical system are investigated, with a particular focus on the accretion disk luminosity in spiral galaxies. The impact of the PBH presence is analyzed through modification of the disk luminosity and kinematic quantities. Thus, we discuss possible constraints on the PBH abundance in view of the most recent theoretical bounds. The results of our study indicate that suitable PBH masses are $M_\text{PBH}\in[10^6,10^{12}]M_\odot$ for PBH fractions $f_\text{PBH}\in[10^{-3},1]$. In particular, a comparison with the predictions of the exponential sphere density profile for dark matter suggests that the best-matching configuration is achieved for $f_\text{PBH}=1$ and $M_\text{PBH}=10^6 M_\odot$. Consequences with respect to the current knowledge on primordial black hole physics are discussed., Comment: 10 pages, 6 figures, Accepted for publication in Phys. Rev. D

Our present understanding of the reheating phase is incomplete due to a lack of observations. Apart from its cosmological implications, the reheating should play a vital role in particle physics and inflation model building. Conventionally reheating dynamics are modeled by invoking arbitrary coupling among the inflaton and daughter fields. Such an approach lacks robust cosmological predictions due to its arbitrary couplings and is difficult to verify through observation. In this paper, we propose a minimal reheating scenario where the inflaton is coupled with all the daughter fields only gravitationally. Besides being successful in reheating the Universe, the scenario offers a strong cosmological prediction of the primordial gravitational wave spectrum and discards a large number of possible models of dark matter and inflation that are otherwise consistent with Planck., 7 pages, 6 figures; v3 matches the published version in PRD

A tantalizing hint of isotropic cosmic birefringence has been found in the $E B$ cross-power spectrum of the cosmic microwave background (CMB) polarization data with a statistical significance of $3\sigma$. A pseudoscalar field coupled to the CMB photons via the Chern-Simons term can explain this observation. The same field may also be responsible for early dark energy (EDE), which alleviates the so-called Hubble tension. Since the EDE field evolves significantly during the recombination epoch, the conventional formula that relates $E B$ to the difference between the $E$- and $B$-mode auto-power spectra is no longer valid. Solving the Boltzmann equation for polarized photons and the dynamics of the EDE field consistently, we find that currently favored parameter space of the EDE model yields a variety of shapes of the $EB$ spectrum, which can be tested by CMB experiments., Comment: 7 pages, 4 figures, 1 table. (v2) Expanded discussion on forecasts for SO and CMB-S4. Accepted for publication in Phys. Rev. D

Gravitational waves offer a new window to probe the nature of gravity, including answering if the mediating particle, graviton, has a non-zero mass or not. Pulsar timing arrays measure stochastic gravitational wave background (SGWB) at $\sim1-100$~nanohertz. Recently, the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) collaboration reported an uncorrelated common-spectrum process in their 12.5-year data set with no substantial evidence that the process comes from the SGWB predicted by general relativity. In this work, we explore the possibility of an SGWB from massive gravity in the data set and find that a massless graviton is preferred because of the relatively larger Bayes factor. Without statistically significant evidence for dispersion-related correlations predicted by massive gravity, we place upper limits on the amplitude of the SGWB for graviton mass smaller than $10^{-23}$~eV as $A_{\rm{MG}}, 7 pages, 2 figures, 1 table; accepted for publication in Phys. Rev. D

Accelerating black holes, connected to cosmic strings could evolve to supermassive black holes. However, if they are going to take part in structure formation and resides at the center of galaxies, their acceleration should be small. This slow acceleration does not change the shadow or image position in gravitational lensing effect significantly. However we show that the time delay associated to these images change significantly. This is in contrast with when the theory governing the strong gravitational field around the black hole is different from general relativity, where not only the differential time delays but the angular position of images would be different. We conclude that, if the observed angular position of images are compatible with the prediction of general relativity, a possible deviation in the differential time delays between the observed values and those predicted by general relativity, could be due to the acceleration of the black hole., v3: 14 pages, 12 figures, extended version of the upcoming PRL paper

We resolve the potential-restriction problem in K/G inflation by introducing nonminimal coupling. In this context, Higgs field successfully drives inflation satisfying CMB observations while enhancing curvature perturbations at small scales, which in turn accounts for primordial black holes (PBHs) and scalar induced gravitational waves (SIGWs). We then uncover the effect of the non-canonical kinetic coupling function in more detail and study its the observational constraint. Besides, we also give the gauge invariant expression for the integral kernel of SIGWs, which is related to terms propagating with the speed of light. Finally, the non-Gaussian effect on PBH abundance and SIGWs is studied. We find that non-Gaussianity makes PBHs form more easily, but its effect on the energy density of SIGWs is negligible.

The final mass distribution of primordial black holes is sensitive to the equation of state of the Universe at the scales accessible by the power spectrum. Motivated by the presence of phase transitions in several beyond the Standard Model theories, some of which are strongly coupled, we analyze the production of primordial black holes during such phase transitions, which we model using the gauge/gravity duality. We focus in the (often regarded as physically uninteresting) case for which the phase transition is just a smooth crossover. We find an enhancement of primordial black hole production in the range $M_{\rm{PBH}}\in[10^{-16},10^{-6}]M_{\odot}$., Comment: 5 pages + references, 4 figures. v3: minor changes, updated to match the published version in PRD

We implement adaptive mesh refinement (AMR) simulations of global topological strings using the public code, GRChombo. We perform a quantitative investigation of massive radiation from single sinusoidally displaced string configurations, studying a range of string widths defined by the coupling parameter $\lambda$ over two orders of magnitude, effectively varying the mass of radiated particles $m_H \sim \sqrt{\lambda}$. We perform an in-depth investigation into the effects of AMR on massive radiation emission, including radiation trapping and the refinement required to resolve high frequency modes. We use quantitative diagnostic tools to determine the eigenmode decomposition, showing a complex superposition of high frequency propagating modes with different phase and group velocities. We conclude that massive radiation is generally strongly suppressed relative to the preferred massless channel, with suppression increasing at lower amplitudes and higher $\lambda$. Only in extreme nonlinear regimes (e.g.\ with relative amplitude $\varepsilon \sim 1.5$ and $\lambda < 1$) do we observe massive and massless radiation to be emitted at comparable magnitude. We find that massive radiation is emitted in distinct high harmonics of the fundamental frequency of the string, and we demonstrate that, for the sinusoidal configurations studied, massive radiation is exponentially suppressed with $\sqrt{\lambda}$ (i.e. the particle mass). Finally, we place these results in the context of axions and gravitational waves produced by cosmological cosmic string networks, and note that AMR provides a significant opportunity to explore higher $\lambda$ (thin string) regimes whilst using fewer computational resources., Comment: 26 pages, 21 figures. Version published Phys. Rev. D

We derive the weak field limit of scalar-Gauss-Bonnet theory and place novel bounds on the parameter space using terrestrial and space-based experiments. In order to analyze the theory in the context of a wide range of experiments, we compute the deviations from Einstein gravity around source masses with planar, cylindrical, and spherical symmetry. We find a correction to the Newtonian potential around spherical and cylindrical sources that can be larger than PPN corrections sufficiently close to the source. We use this to improve on laboratory constraints on the scalar-Gauss-Bonnet coupling parameter $\Lambda$ by two orders of magnitude. Present laboratory and Solar System bounds reported here are superseded by tests deriving from black holes., Comment: 18 pages

Pulsar timing offers an independent avenue to test general relativity and alternative gravity theories. This requires an understanding of how metric polarizations beyond the familiar transverse tensor ones imprint as a stochastic gravitational wave background and correlate the arrival time of radio pulses from a pair of millisecond pulsars. In this work, we focus on an isotropic stochastic gravitational wave background and present a straightforward, self-contained formalism for obtaining the power spectrum and the overlap reduction function, the relevant physical observable in a pulsar timing array, for generic gravitational degrees of freedom featuring both transverse and longitudinal modes off the light cone. We additionally highlight our consideration of finite pulsar distances, which we find significant in two ways: first, making all the modes well defined, and second, keeping the small scale power that is contained by pulsars of subdegree separations in the sky. We discuss this for tensor, vector, and scalar polarizations, for each one focusing on the angular power spectrum and the overlap reduction function for an isotropic stochastic gravitational wave background. Our results pave the road for an efficient numerical method for examining the gravitational wave induced spatial correlations across millisecond pulsars in a pulsar timing array., Comment: 28 pages, 5 figures, discussion improved, to appear in PRD, see implementation `PTAfast' https://github.com/reggiebernardo/PTAfast

Since general relativity is the unique theory of massless spin 2 particles at large distances, the most reasonable way to have significant modifications is to introduce one or more light scalars that mediate a new long-range force. Most existing studies of such scalars invoke models that exhibit some kind of "screening" at short distances to hide the force from solar system tests. However, as is well known, such modifications also exhibit superluminality, which can be interpreted as a form of acausality. In this work we explore explicitly subluminal and causal scalar field models. In particular, we study a conformally coupled scalar $\phi$, with a small coupling to matter to obey solar system bounds, and a non-canonical kinetic term $K(X)$ ($X=(\partial\phi)^2/2$) that obeys all subluminality constraints and is hyperbolic. We consider $K(X)$ that is canonical for small $X$, but beyond some nonlinear scale enters a new scaling regime of power $p$, with $1/2, Comment: 44 pages, 6 figures. V3: Updated toward version published in PRD

The interaction of binary black hole mergers with their environments can be studied using numerical relativity simulations. These start only a short finite time before merger, at which point appropriate initial conditions must be imposed. A key task is therefore to identify the configuration that is appropriate for the binary and its environment at this stage of the evolution. In this work we study the behaviour of wave dark matter around equal mass black hole binaries, finding that there is a preferred, quasi-stationary profile that persists and grows over multiple orbits, in contrast to heavier mass dark matter where any overdensity tends to be dispersed by the binary motion. Whilst different initial configurations converge to the preferred quasi-stationary one after several orbits, unwanted transient oscillations are generated in the process, which may impact on the signal in short simulation runs. We also point out that naively superimposing the matter onto a circular binary results in artificially eccentric orbits due to the matter backreaction, which is an effect of the initial conditions and not a signature of dark matter. We discuss the further work required so that comparison of waveforms obtained with environments to vacuum cases can be done in a meaningful way., 9 pages, 5 figures, 1 appendix, 2 movies: https://youtube.com/playlist?list=PLSkfizpQDrcb8Zd7C9yEfnMpnKQk7I9ts Comments welcome!

We derive the fixed-Λ and unimodular propagators using the path integral formalism as applied to the Einstein-Cartan action. The simplicity of the action (which is linear in the lapse function) allows for an exact integration starting from the lapse function and the enforcement of the Hamiltonian constraint, leading to a product of Chern-Simons states if the connection is fixed at the endpoints. No saddle point approximation is needed. Should the metric be fixed at the endpoints, then, depending on the contour chosen for the connection, Hartle-Hawking or Vilenkin propagators are obtained. Thus, in this approach one trades a choice of contour in the lapse function for one in the connection, where appropriate. The unimodular propagators are also trivial to obtain via the path integral, and the previously derived expressions are recovered.

Could new physics break the mirror symmetry of the Universe? Utilizing recent measurements of the parity-odd four-point correlation function of BOSS galaxies, we probe the physics of inflation by placing constraints on the amplitude of a number of parity-violating models. Within canonical models of (single-field, slow-roll) inflation, no parity-asymmetry can occur; however, it has recently been shown that breaking of the standard assumptions can lead to parity violation within the Effective Field Theory of Inflation (EFTI). In particular, we consider the Ghost Condensate and Cosmological Collider scenarios -- the former for the leading and subleading operators in the EFTI and the latter for different values of mass and speed of an exchanged spin-$1$ particle -- for a total of $18$ models. Each instance yields a definite prediction for the inflationary trispectrum, which we convert to a late-time galaxy correlator prediction (through a highly non-trivial calculation) and constrain using the observed data. We find no evidence for inflationary parity-violation (with each of the $18$ models having significances below $2\sigma$), and place the first constraints on the relevant coupling strengths, at a level comparable with the theoretical perturbativity bounds. This is also the first time Cosmological Collider signatures have directly been searched for in observational data. We further show that possible secondary parity-violating signatures in galaxy clustering can be systematically described within the Effective Field Theory of Large-Scale Structure. We argue that these late-time contributions are subdominant compared to the primordial parity-odd signal for a vast region of parameter space. In summary, the results of this paper disfavor the notion that the recent hints of parity-violation observed in the distribution of galaxies are due to new physics., Comment: 20+9 spooky pages, 4 figures, accepted by Phys. Rev. D. Code available at https://github.com/oliverphilcox/Parity-Odd-4PCF

We systematically investigate the microlensing effect of charged spherically symmetric wormhole, where the light source is remote from the throat. Remarkably, there will be at most three images by considering the charge part. We study all situations including three images, two images, and one image, respectively. The numerical result shows that the range of total magnification is from $10^5$ to $10^{-2}$ depending on various metrics. In the case of three images, there will be two maximal values of magnification (a peak, and a gentle peak) when the contribution via mass is much less than that of charge. However, we cannot distinguish the case that forms three images or only one image as the total magnification is of order $10^5$. Finally, our theoretical investigation could shed new light on exploring the wormhole with the microlensing effect., 10 pages, 9 figures

We formulate the transition from decelerated to accelerated expansion as a bounce in connection space and study its quantum cosmology, knowing that reflections are notorious for bringing quantum effects to the fore. We use a formalism for obtaining a time variable via the demotion of the constants of Nature to integration constants, and focus on a toy Universe containing only radiation and a cosmological constant $\Lambda$ for its simplicity. We find that, beside the usual factor ordering ambiguities, there is an ambiguity in the order of the quantum equation, leading to two distinct theories: one second, the other first order. In both cases two time variables may be defined, conjugate to $\Lambda$ and to the radiation constant of motion. We make little headway with the second-order theory, but are able to produce solutions to the first-order theory. They exhibit the well-known "ringing" whereby incident and reflected waves interfere, leading to oscillations in the probability distribution even for well-peaked wave packets. We also examine in detail the probability measure within the semiclassical approximation. Close to the bounce, the probability distribution becomes double-peaked, with one peak following a trajectory close to the classical limit but with a Hubble parameter slightly shifted downwards, and the other with a value of $b$ stuck at its minimum. An examination of the effects still closer to the bounce, and within a more realistic model involving matter and $\Lambda$, is left to future work., Comment: 15 pages, 13 figures; v3: very minor changes in presentation (mostly typographical), close to version published in Phys.Rev.D

There are many single field inflationary models that are consistent with the recent Planck 2018 measurements of the spectral index $n_s$ and tensor-to-scalar ratio $r$. Despite good agreement with observational data some of these models suffer from having unregularized potentials which would produce a collapsing universe shortly after the end of inflation. In this paper we show that how one chooses to correct the behaviour potential towards the end of inflation can have a significant effect on the inflationary predictions of the model, specifically in the case of quartic hilltop and radiatively corrected Higgs inflation., Comment: 25 pages, 25 figures

The measurement of the inflationary stochastic gravitational-wave background (SGWB) is one of the main goals of future GW experiments. In direct GW experiments, an obstacle to achieving it is the isolation of the inflationary SGWB from the other types of SGWB. In this paper, as a distinguishable signature of the inflationary SGWB, we argue the detectability of its universal property: antipodal correlations, i.e., correlations of GWs from the opposite directions, as a consequence of the horizon re-entry. A phase-coherent method has been known to be of no use for detecting the angular correlations in SGWB due to a problematic phase factor that erases the signal. We thus investigate whether we can construct a phase-incoherent estimator of the antipodal correlations in the intensity map. We found that the conclusion depends on whether the inflationary GWs have statistical isotropy or not. In the standard inflationary models with statistical homogeneity and isotropy, there is no estimator that is sensitive to the antipodal correlations but does not suffer from the problematic phase factor. On the other hand, it is possible to find a non-vanishing estimator of the antipodal correlations for inflationary models with statistical anisotropy. SGWB from anisotropic inflation is distinguishable from the other components., Comment: 8+4 pages, 2 figures; published version, added discussions on the detectability for anisotropic inflation

We consider ${\cal{R}}^2$-inflation in Palatini gravity, in the presence of scalar fields coupled to gravity. These theories, in the Einstein frame, and for one scalar field $h$, share common features with $K$ - inflation models. We apply this formalism for the study of single-field inflationary models, whose potentials are monomials, $ V \sim h^{n} $, with $ n $ a positive even integer. We also study the Higgs model non-minimally coupled to gravity. With ${\cal{R}}^2$-terms coupled to gravity as $\sim \alpha {\cal{R}}^2 $, with $\alpha$ constant, the instantaneous reheating temperature $T_{ins}$, is bounded by $ T_{ins} \leq { 0.290 \, m_{Planck}} / {\, \alpha^{1/4}} $, with the upper bound being saturated for large $\alpha$. For such large $\alpha$ need go beyond slow-roll to calculate reliably the cosmological parameters, among these the end of inflation through which $T_{ins}$ is determined. In fact, as inflaton rolls towards the end of inflation point, the quartic in the velocity terms, unavoidable in Palatini gravity, play a significant role and can not be ignored. The values of $\alpha$, and other parameters, are constrained by cosmological data, setting bounds on the inflationary scale $M_{s} \sim 1/\sqrt{\alpha}$ and the reheating temperature of the Universe., Comment: 35 pages, 18 figures, PDFLaTeX, text and references added, typos corrected

This paper contains a detailed study of the properties of a simple model attempting to explain dark energy as originated from quantum fluctuations of a light spectator scalar field in inflation. In [1] we recently outlined how Starobinsky's stochastic formalism can be used to study the spatial correlations imprinted on dark energy by its quantum origin in this model and we studied their possible role in relieving the Hubble tension. Here we provide a more comprehensive derivation of the results in [1] and we refine some of our estimates, comparing to the approximate results obtained previously. Among the main results, we analyze the non-coincident correlators predicted by a full field theoretical treatment and their relation with those computed within the stochastic formalism. We find that in the region where stochastic theory predicts significant sub-Hubble correlators it is in disagreement with field theoretical predictions. However, agreement can be restored by introducing a reduced speed of sound for the scalar field. We also discuss an alternative approach to the problem of studying correlators within the stochastic formalism based directly on the evolution of probability distributions. We find that the two approaches give the same answer for 2-point functions of the field, but not for 4-point functions relevant to density correlators and we discuss the behaviour of the two methods with respect to Wick's theorem., v2: version submitted to Phys. Rev. D

We investigate the impact of radiation from primordial black holes (PBHs), in the mass range of $10^{15} \lesssim M_{\rm PBH} \lesssim 10^{17}~\rm g$ and $10^{2} \lesssim M_{\rm PBH} \lesssim 10^{4}~M_{\odot}$, on the 21-cm angular-power spectrum in the dark ages. PBHs in the former mass range effect the 21-cm angular-power spectrum through the evaporation known as Hawking radiation, while the radiation from the accretion process in the latter mass range. In the dark ages, radiation from PBHs can increase the ionization fraction and temperature of the intergalactic medium, change the global 21-cm differential brightness temperature and then effect the 21-cm angular-power spectrum. Taking into account the effects of PBHs, we find that in the dark ages, $30 \lesssim z \lesssim 100$, the amplitude of the 21-cm angular-power spectrum is decreased depending on the mass and mass fraction of PBHs. We also investigate the potential constraints on the mass fraction of PBHs in the form of dark matter for the future radio telescope in lunar obit or on the farside surface of the Moon., 8 pages, 6 figures. Title changed slightly. Discussions and references added. Published in Physical Review D

Hydrodynamic and magnetohydrodynamic turbulence in the early Universe can drive gravitational waves (GWs) and imprint their spectrum onto that of GWs, which might still be observable today. We study the production of the GW background from freely decaying magnetohydrodynamic turbulence from helical and nonhelical initial magnetic fields. To understand the produced GW spectra, we develop a simple model on the basis of the evolution of the magnetic stress tensor. We find that the GW spectra obtained in this model reproduce those obtained in numerical simulations if we consider the detailed time evolution of the low-frequency tail of the stress spectrum from numerical simulations. We also show that the shapes of the produced GW frequency spectra are different for helical and nonhelical cases for the same initial magnetic energy spectra. Such differences can help distinguish helical and nonhelical initial magnetic fields from a polarized background of GWs -- especially when the expected circular polarization cannot be detected directly., 12 pages, 12 figures, 3 tables, Accepted for publication in Phys. Rev. D

Recently, the cross-correlation between $E$- and $B$-mode polarization of the cosmic microwave background (CMB), which is well explained by cosmic birefringence with rotation angle $\beta\approx 0.3$ deg, has been found in CMB polarization data. We carefully investigate the possibility of explaining the observed $EB$ correlation by the primordial chiral gravitational waves (CGWs), which can be generated in the parity-violating theories in the primordial Universe. We found that the CGWs scenario does not work due to the overproduction of the $BB$ auto-correlation which far exceeds the observed one by SPTPol and POLARBEAR., Comment: 5 pages, 2 figures

We have recently proposed a simple relativistic theory which reduces to modified Newtonian dynamics for the weak-field quasistatic situations applied to galaxies, and to cosmological behavior as in the $\Lambda$CDM model, yielding a realistic cosmology in line with observations. A key requirement of any such model is that Minkowski space is stable against linear perturbations. We expand the theory action to second order in perturbations on a Minkowski background and show that it leads to healthy dispersion relations involving propagating massive modes in the vector and the scalar sector. We use Hamiltonian methods to eliminate constraints present, demonstrate that the massive modes have Hamiltonian bounded from below and show that a nonpropagating mode with a linear time dependence may have unbounded Hamiltonian for wave numbers $k< \mu$ and bounded otherwise. The scale $\mu$ is estimated to be $\lesssim \mathrm{Mpc}^{-1}$ so that the low momenta instability may only play a role on cosmological scales., Comment: Change of paper title; included an example of a free function; a new appendix included to indicate the effect of the MOND term on results; convention for the definition of canonical momenta in Fourier space has been revised without change to the final results

Some cosmological models with non-negligible dark energy fractions in particular windows of the pre-recombination epoch are capable of alleviating the Hubble tension quite efficiently, while keeping the good description of the data that are used to build the cosmic inverse distance ladder. There has been an intensive discussion in the community on whether these models enhance the power of matter fluctuations, leading {\it de facto} to a worsening of the tension with the large-scale structure measurements. We address this pivotal question in the context of several early dark energy (EDE) models, considering also in some cases a coupling between dark energy and dark matter, and the effect of massive neutrinos. We fit them using the Planck 2018 likelihoods, the supernovae of Type Ia from the Pantheon compilation and data on baryon acoustic oscillations. We find that ultra-light axion-like (ULA) EDE can actually alleviate the $H_0$ tension without increasing the values of $\sigma_{12}$ with respect to those found in the $\Lambda$CDM, whereas EDE with an exponential potential does not have any impact on the tensions. A coupling in the dark sector tends to enhance the clustering of matter, and the data limit a lot the influence of massive neutrinos, since the upper bounds on the sum of their masses are too close to those obtained in the standard model. We find that in the best case, namely ULA, the Hubble tension is reduced to $\sim 2\sigma$., Comment: 27 pages, 16 figures, 8 tables. Version accepted for publication in Phys. Rev. D

We revisit D-term inflation with the bounds of the cosmic string tension from gravitational wave observations and consider the possible deviation of the spectral index compared with the $\Lambda$CDM model in light of prerecombination resolutions of the Hubble tension. D-term inflation requires very small coupling constants under these constraints. We show that natural coupling constants ($g=\lambda=0.1$) can be achieved in the case of D-term inflation on the brane., Comment: 11 pages, to be published in Physical Review D

Composed of ultralight bosons, fuzzy dark matter provides an intriguing solution to challenges that the standard cold dark matter model encounters on sub-galactic scales. The ultralight dark matter with mass $m\sim10^{-23} \rm{eV}$ will induce a periodic oscillation in gravitational potentials with a frequency in the nanohertz band, leading to observable effects in the arrival times of radio pulses from pulsars. Unlike scalar dark matter, pulsar timing signals induced by the vector dark matter are dependent on the oscillation direction of the vector fields. In this work, we search for ultralight vector dark matter in the mass range of $[2\times 10^{-24}, 2\times 10^{-22}]{\rm{eV}}$ through its gravitational effect in the Parkes Pulsar Timing Array (PPTA) second data release. Since no statistically significant detection is made, we place $95\%$ upper limits on the local dark matter density as $\rho_{\rm{\tiny{VF}}} \lesssim 5{\rm{GeV/cm^{3}}}$ for $m\lesssim 10^{-23}{\rm{eV}}$. As no preferred direction is found for the vector dark matter, these constraints are comparable to those given by the scalar dark matter search with an earlier 12-year data set of PPTA., Comment: 9 pages, 3 figures, 2 tables; accepted for publication as a Letter in Phys. Rev. D

In this work, we introduce a class of extended Minimal Theories of Massive Gravity (eMTMG), without requiring a priori that the theory should admit the same homogeneous and isotropic cosmological solutions as the de Rham-Gabadadze-Tolley massive gravity. The theory is constructed as to have only two degrees of freedom in the gravity sector. In order to perform this step we first introduce a precursor theory endowed with a general graviton mass term, to which, at the level of the Hamiltonian, we add two extra constraints as to remove the unwanted degrees of freedom, which otherwise would typically lead to ghosts and/or instabilities. On analyzing the number of independent constraints and the properties of tensor mode perturbations, we see that the gravitational waves are the only propagating gravitational degrees of freedom which do acquire a non-trivial mass, as expected. In order to understand how the effective gravitational force works for this theory we then investigate cosmological scalar perturbations in the presence of a pressureless fluid. We then restrict the whole class of models by imposing the following conditions at all times: 1) it is possible to define an effective gravitational constant, $G_{{\rm eff}}$; 2) the value $G_{\text{eff}}/G_{N}$ is always finite but not always equal to unity (as to allow some non-trivial modifications of gravity, besides the massive tensorial modes); and 3) the square of mass of the graviton is always positive. These constraints automatically make also the ISW-effect contributions finite at all times. Finally we focus on a simple subclass of such theories, and show they already can give a rich and interesting phenomenology., Comment: 20 pages+appendices, 1 figure

This work investigates the formation of primordial black holes within a radiation fluid with an anisotropic pressure. We focus our attention on the initial conditions describing cosmological perturbations in the super horizon regime, using a covariant form of the equation of state in terms of pressure and energy density gradients. The effect of the anisotropy is to modify the initial shape of the cosmological perturbations with respect to the isotropic case. Using the dependence of the threshold $\delta_\mathrm{c}$ for primordial black holes with respect to the shape of cosmological perturbations, we estimate here how the threshold is varying with respect to the amplitude of the anisotropy. If this variation is large enough it could lead to a significant variation of the abundance of PBHs., Comment: Published in Physical Review D, Version 4: 16 pages, 7 figures, minor changes in Sec. III