Showing 1–36 of 36 results for author: Unal, C

The NANOGrav 15 yr Data Set: Looking for Signs of Discreteness in the Gravitational-wave Background

Authors: Gabriella Agazie, Paul T. Baker, Bence Bécsy, Laura Blecha, Adam Brazier, Paul R. Brook, Lucas Brown, Sarah Burke-Spolaor, J. Andrew Casey-Clyde, Maria Charisi, Shami Chatterjee, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Megan E. DeCesar, Paul B. Demorest, Heling Deng, Timothy Dolch, Elizabeth C. Ferrara, William Fiore, Emmanuel Fonseca, Gabriel E. Freedman, Nate Garver-Daniels , et al. (58 additional authors not shown)

Abstract: The cosmic merger history of supermassive black hole binaries (SMBHBs) is expected to produce a low-frequency gravitational wave background (GWB). Here we investigate how signs of the discrete nature of this GWB can manifest in pulsar timing arrays through excursions from, and breaks in, the expected $f_{\mathrm{GW}}^{-2/3}$ power-law of the GWB strain spectrum. To do this, we create a semi-analyt… ▽ More The cosmic merger history of supermassive black hole binaries (SMBHBs) is expected to produce a low-frequency gravitational wave background (GWB). Here we investigate how signs of the discrete nature of this GWB can manifest in pulsar timing arrays through excursions from, and breaks in, the expected $f_{\mathrm{GW}}^{-2/3}$ power-law of the GWB strain spectrum. To do this, we create a semi-analytic SMBHB population model, fit to NANOGrav's 15 yr GWB amplitude, and with 1,000 realizations we study the populations' characteristic strain and residual spectra. Comparing our models to the NANOGrav 15 yr spectrum, we find two interesting excursions from the power-law. The first, at $2 \; \mathrm{nHz}$, is below our GWB realizations with $p$-value significance $p = 0.05$ to $0.06$ ($\approx 1.8 σ- 1.9 σ$). The second, at $16 \; \mathrm{nHz}$, is above our GWB realizations with $p = 0.04$ to $0.15$ ($\approx 1.4 σ- 2.1 σ$). We explore the properties of a loud SMBHB which could cause such an excursion. Our simulations also show that the expected number of SMBHBs decreases by three orders of magnitude, from $\sim 10^6$ to $\sim 10^3$, between $2\; \mathrm{nHz}$ and $20 \; \mathrm{nHz}$. This causes a break in the strain spectrum as the stochasticity of the background breaks down at $26^{+28}_{-19} \; \mathrm{nHz}$, consistent with predictions pre-dating GWB measurements. The diminished GWB signal from SMBHBs at frequencies above the $26$~nHz break opens a window for PTAs to detect continuous GWs from individual SMBHBs or GWs from the early universe. △ Less

Submitted 10 April, 2024; originally announced April 2024.

Comments: 10 pages, 8 figures, 1 appendix, submitted to ApJ

Primordial black holes and their gravitational-wave signatures

Authors: Eleni Bagui, Sebastien Clesse, Valerio De Luca, Jose María Ezquiaga, Gabriele Franciolini, Juan García-Bellido, Cristian Joana, Rajeev Kumar Jain, Sachiko Kuroyanagi, Ilia Musco, Theodoros Papanikolaou, Alvise Raccanelli, Sébastien Renaux-Petel, Antonio Riotto, Ester Ruiz Morales, Marco Scalisi, Olga Sergijenko, Caner Unal, Vincent Vennin, David Wands

Abstract: In the recent years, primordial black holes (PBHs) have emerged as one of the most interesting and hotly debated topics in cosmology. Among other possibilities, PBHs could explain both some of the signals from binary black hole mergers observed in gravitational wave detectors and an important component of the dark matter in the Universe. Significant progress has been achieved both on the theory si… ▽ More In the recent years, primordial black holes (PBHs) have emerged as one of the most interesting and hotly debated topics in cosmology. Among other possibilities, PBHs could explain both some of the signals from binary black hole mergers observed in gravitational wave detectors and an important component of the dark matter in the Universe. Significant progress has been achieved both on the theory side and from the point of view of observations, including new models and more accurate calculations of PBH formation, evolution, clustering, merger rates, as well as new astrophysical and cosmological probes. In this work, we review, analyse and combine the latest developments in order to perform end-to-end calculations of the various gravitational wave signatures of PBHs. Different ways to distinguish PBHs from stellar black holes are emphasized. Finally, we discuss their detectability with LISA, the first planned gravitational-wave observatory in space. △ Less

Submitted 30 October, 2023; originally announced October 2023.

Comments: 161 pages, 47 figures, comments welcome

The NANOGrav 15-year data set: Search for Transverse Polarization Modes in the Gravitational-Wave Background

Authors: Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Jeremy Baier, Paul T. Baker, Bence Bécsy, Laura Blecha, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Rand Burnette, Robin Case, J. Andrew Casey-Clyde, Maria Charisi, Shami Chatterjee, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Kathryn Crowter, Megan E. DeCesar, Dallas DeGan, Paul B. Demorest , et al. (74 additional authors not shown)

Abstract: Recently we found compelling evidence for a gravitational wave background with Hellings and Downs (HD) correlations in our 15-year data set. These correlations describe gravitational waves as predicted by general relativity, which has two transverse polarization modes. However, more general metric theories of gravity can have additional polarization modes which produce different interpulsar correl… ▽ More Recently we found compelling evidence for a gravitational wave background with Hellings and Downs (HD) correlations in our 15-year data set. These correlations describe gravitational waves as predicted by general relativity, which has two transverse polarization modes. However, more general metric theories of gravity can have additional polarization modes which produce different interpulsar correlations. In this work we search the NANOGrav 15-year data set for evidence of a gravitational wave background with quadrupolar Hellings and Downs (HD) and Scalar Transverse (ST) correlations. We find that HD correlations are the best fit to the data, and no significant evidence in favor of ST correlations. While Bayes factors show strong evidence for a correlated signal, the data does not strongly prefer either correlation signature, with Bayes factors $\sim 2$ when comparing HD to ST correlations, and $\sim 1$ for HD plus ST correlations to HD correlations alone. However, when modeled alongside HD correlations, the amplitude and spectral index posteriors for ST correlations are uninformative, with the HD process accounting for the vast majority of the total signal. Using the optimal statistic, a frequentist technique that focuses on the pulsar-pair cross-correlations, we find median signal-to-noise-ratios of 5.0 for HD and 4.6 for ST correlations when fit for separately, and median signal-to-noise-ratios of 3.5 for HD and 3.0 for ST correlations when fit for simultaneously. While the signal-to-noise-ratios for each of the correlations are comparable, the estimated amplitude and spectral index for HD are a significantly better fit to the total signal, in agreement with our Bayesian analysis. △ Less

Submitted 18 October, 2023; originally announced October 2023.

Comments: 11 pages, 5 figures

The NANOGrav 12.5-year data set: A computationally efficient eccentric binary search pipeline and constraints on an eccentric supermassive binary candidate in 3C 66B

Authors: Gabriella Agazie, Zaven Arzoumanian, Paul T. Baker, Bence Bécsy, Laura Blecha, Harsha Blumer, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, J. Andrew Casey-Clyde, Maria Charisi, Shami Chatterjee, Belinda D. Cheeseboro, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Megan E. DeCesar, Paul B. Demorest, Lankeswar Dey, Timothy Dolch, Justin A. Ellis, Robert D. Ferdman, Elizabeth C. Ferrara , et al. (63 additional authors not shown)

Abstract: The radio galaxy 3C 66B has been hypothesized to host a supermassive black hole binary (SMBHB) at its center based on electromagnetic observations. Its apparent 1.05-year period and low redshift ($\sim0.02$) make it an interesting testbed to search for low-frequency gravitational waves (GWs) using Pulsar Timing Array (PTA) experiments. This source has been subjected to multiple searches for contin… ▽ More The radio galaxy 3C 66B has been hypothesized to host a supermassive black hole binary (SMBHB) at its center based on electromagnetic observations. Its apparent 1.05-year period and low redshift ($\sim0.02$) make it an interesting testbed to search for low-frequency gravitational waves (GWs) using Pulsar Timing Array (PTA) experiments. This source has been subjected to multiple searches for continuous GWs from a circular SMBHB, resulting in progressively more stringent constraints on its GW amplitude and chirp mass. In this paper, we develop a pipeline for performing Bayesian targeted searches for eccentric SMBHBs in PTA data sets, and test its efficacy by applying it on simulated data sets with varying injected signal strengths. We also search for a realistic eccentric SMBHB source in 3C 66B using the NANOGrav 12.5-year data set employing PTA signal models containing Earth term-only as well as Earth+Pulsar term contributions using this pipeline. Due to limitations in our PTA signal model, we get meaningful results only when the initial eccentricity $e_0<0.5$ and the symmetric mass ratio $η>0.1$. We find no evidence for an eccentric SMBHB signal in our data, and therefore place 95% upper limits on the PTA signal amplitude of $88.1\pm3.7$ ns for the Earth term-only and $81.74\pm0.86$ ns for the Earth+Pulsar term searches for $e_0<0.5$ and $η>0.1$. Similar 95% upper limits on the chirp mass are $(1.98 \pm 0.05) \times 10^9\,M_{\odot}$ and $(1.81 \pm 0.01) \times 10^9\,M_{\odot}$. These upper limits, while less stringent than those calculated from a circular binary search in the NANOGrav 12.5-year data set, are consistent with the SMBHB model of 3C 66B developed from electromagnetic observations. △ Less

Submitted 15 January, 2024; v1 submitted 29 September, 2023; originally announced September 2023.

Comments: 27 Pages, 10 Figures, 1 Table, Accepted for publication in ApJ

How to Detect an Astrophysical Nanohertz Gravitational-Wave Background

Authors: Bence Bécsy, Neil J. Cornish, Patrick M. Meyers, Luke Zoltan Kelley, Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Paul T. Baker, Laura Blecha, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, J. Andrew Casey-Clyde, Maria Charisi, Shami Chatterjee, Katerina Chatziioannou, Tyler Cohen, James M. Cordes, Fronefield Crawford, H. Thankful Cromartie, Kathryn Crowter, Megan E. DeCesar, Paul B. Demorest, Timothy Dolch , et al. (71 additional authors not shown)

Abstract: Analysis of pulsar timing data have provided evidence for a stochastic gravitational wave background in the nHz frequency band. The most plausible source of such a background is the superposition of signals from millions of supermassive black hole binaries. The standard statistical techniques used to search for such a background and assess its significance make several simplifying assumptions, nam… ▽ More Analysis of pulsar timing data have provided evidence for a stochastic gravitational wave background in the nHz frequency band. The most plausible source of such a background is the superposition of signals from millions of supermassive black hole binaries. The standard statistical techniques used to search for such a background and assess its significance make several simplifying assumptions, namely: i) Gaussianity; ii) isotropy; and most often iii) a power-law spectrum. However, a stochastic background from a finite collection of binaries does not exactly satisfy any of these assumptions. To understand the effect of these assumptions, we test standard analysis techniques on a large collection of realistic simulated datasets. The dataset length, observing schedule, and noise levels were chosen to emulate the NANOGrav 15-year dataset. Simulated signals from millions of binaries drawn from models based on the Illustris cosmological hydrodynamical simulation were added to the data. We find that the standard statistical methods perform remarkably well on these simulated datasets, despite their fundamental assumptions not being strictly met. They are able to achieve a confident detection of the background. However, even for a fixed set of astrophysical parameters, different realizations of the universe result in a large variance in the significance and recovered parameters of the background. We also find that the presence of loud individual binaries can bias the spectral recovery of the background if we do not account for them. △ Less

Submitted 1 December, 2023; v1 submitted 8 September, 2023; originally announced September 2023.

Comments: 14 pages, 8 figures, version matching published paper

Journal ref: ApJ 959 9 (2023)

Comparing recent PTA results on the nanohertz stochastic gravitational wave background

Authors: The International Pulsar Timing Array Collaboration, G. Agazie, J. Antoniadis, A. Anumarlapudi, A. M. Archibald, P. Arumugam, S. Arumugam, Z. Arzoumanian, J. Askew, S. Babak, M. Bagchi, M. Bailes, A. -S. Bak Nielsen, P. T. Baker, C. G. Bassa, A. Bathula, B. Bécsy, A. Berthereau, N. D. R. Bhat, L. Blecha, M. Bonetti, E. Bortolas, A. Brazier, P. R. Brook, M. Burgay , et al. (220 additional authors not shown)

Abstract: The Australian, Chinese, European, Indian, and North American pulsar timing array (PTA) collaborations recently reported, at varying levels, evidence for the presence of a nanohertz gravitational wave background (GWB). Given that each PTA made different choices in modeling their data, we perform a comparison of the GWB and individual pulsar noise parameters across the results reported from the PTA… ▽ More The Australian, Chinese, European, Indian, and North American pulsar timing array (PTA) collaborations recently reported, at varying levels, evidence for the presence of a nanohertz gravitational wave background (GWB). Given that each PTA made different choices in modeling their data, we perform a comparison of the GWB and individual pulsar noise parameters across the results reported from the PTAs that constitute the International Pulsar Timing Array (IPTA). We show that despite making different modeling choices, there is no significant difference in the GWB parameters that are measured by the different PTAs, agreeing within $1σ$. The pulsar noise parameters are also consistent between different PTAs for the majority of the pulsars included in these analyses. We bridge the differences in modeling choices by adopting a standardized noise model for all pulsars and PTAs, finding that under this model there is a reduction in the tension in the pulsar noise parameters. As part of this reanalysis, we "extended" each PTA's data set by adding extra pulsars that were not timed by that PTA. Under these extensions, we find better constraints on the GWB amplitude and a higher signal-to-noise ratio for the Hellings and Downs correlations. These extensions serve as a prelude to the benefits offered by a full combination of data across all pulsars in the IPTA, i.e., the IPTA's Data Release 3, which will involve not just adding in additional pulsars, but also including data from all three PTAs where any given pulsar is timed by more than as single PTA. △ Less

Submitted 1 September, 2023; originally announced September 2023.

Comments: 21 pages, 9 figures, submitted to ApJ

The NANOGrav 12.5-year Data Set: Search for Gravitational Wave Memory

Authors: Gabriella Agazie, Zaven Arzoumanian, Paul T. Baker, Bence Bécsy, Laura Blecha, Harsha Blumer, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Rand Burnette, Robin Case, J. Andrew Casey-Clyde, Maria Charisi, Shami Chatterjee, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Megan E. DeCesar, Dallas DeGan, Paul B. Demorest, Timothy Dolch, Brendan Drachler, Justin A. Ellis , et al. (65 additional authors not shown)

Abstract: We present the results of a Bayesian search for gravitational wave (GW) memory in the NANOGrav 12.5-yr data set. We find no convincing evidence for any gravitational wave memory signals in this data set (Bayes factor = 2.8). As such, we go on to place upper limits on the strain amplitude of GW memory events as a function of sky location and event epoch. These upper limits are computed using a sign… ▽ More We present the results of a Bayesian search for gravitational wave (GW) memory in the NANOGrav 12.5-yr data set. We find no convincing evidence for any gravitational wave memory signals in this data set (Bayes factor = 2.8). As such, we go on to place upper limits on the strain amplitude of GW memory events as a function of sky location and event epoch. These upper limits are computed using a signal model that assumes the existence of a common, spatially uncorrelated red noise in addition to a GW memory signal. The median strain upper limit as a function of sky position is approximately $3.3 \times 10^{-14}$. We also find that there are some differences in the upper limits as a function of sky position centered around PSR J0613$-$0200. This suggests that this pulsar has some excess noise which can be confounded with GW memory. Finally, the upper limits as a function of burst epoch continue to improve at later epochs. This improvement is attributable to the continued growth of the pulsar timing array. △ Less

Submitted 25 July, 2023; originally announced July 2023.

Comments: 29 pages, 5 figures

Axion-Gauge Dynamics During Inflation as the Origin of Pulsar Timing Array Signals and Primordial Black Holes

Authors: Caner Unal, Alexandros Papageorgiou, Ippei Obata

Abstract: We demonstrate that the recently announced signal for a stochastic gravitational wave background (SGWB) from pulsar timing array (PTA) observations, if attributed to new physics, is compatible with primordial GW production due to axion-gauge dynamics during inflation. More specifically we find that axion-$U(1)$ models may lead to sufficient particle production to explain the signal while simultane… ▽ More We demonstrate that the recently announced signal for a stochastic gravitational wave background (SGWB) from pulsar timing array (PTA) observations, if attributed to new physics, is compatible with primordial GW production due to axion-gauge dynamics during inflation. More specifically we find that axion-$U(1)$ models may lead to sufficient particle production to explain the signal while simultaneously source some fraction of sub-solar mass primordial black holes (PBHs) as a signature. Moreover there is a parity violation in GW sector, hence the model suggests chiral GW search as a concrete target for future. We further analyze the axion-$SU(2)$ coupling signatures and find that in the low/mild backreaction regime, it is incapable of producing PTA evidence and the tensor-to-scalar ratio is low at the peak, hence it overproduces scalar perturbations and PBHs. △ Less

Submitted 20 July, 2023; v1 submitted 5 July, 2023; originally announced July 2023.

Comments: 5 pages, 2 figures, total 7 pages, comments welcome

The NANOGrav 15-year Gravitational-Wave Background Analysis Pipeline

Authors: Aaron D. Johnson, Patrick M. Meyers, Paul T. Baker, Neil J. Cornish, Jeffrey S. Hazboun, Tyson B. Littenberg, Joseph D. Romano, Stephen R. Taylor, Michele Vallisneri, Sarah J. Vigeland, Ken D. Olum, Xavier Siemens, Justin A. Ellis, Rutger van Haasteren, Sophie Hourihane, Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Laura Blecha, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Bence Bécsy, J. Andrew Casey-Clyde , et al. (71 additional authors not shown)

Abstract: This paper presents rigorous tests of pulsar timing array methods and software, examining their consistency across a wide range of injected parameters and signal strength. We discuss updates to the 15-year isotropic gravitational-wave background analyses and their corresponding code representations. Descriptions of the internal structure of the flagship algorithms \texttt{Enterprise} and \texttt{P… ▽ More This paper presents rigorous tests of pulsar timing array methods and software, examining their consistency across a wide range of injected parameters and signal strength. We discuss updates to the 15-year isotropic gravitational-wave background analyses and their corresponding code representations. Descriptions of the internal structure of the flagship algorithms \texttt{Enterprise} and \texttt{PTMCMCSampler} are given to facilitate understanding of the PTA likelihood structure, how models are built, and what methods are currently used in sampling the high-dimensional PTA parameter space. We introduce a novel version of the PTA likelihood that uses a two-step marginalization procedure that performs much faster when the white noise parameters remain fixed. We perform stringent tests of consistency and correctness of the Bayesian and frequentist analysis software. For the Bayesian analysis, we test prior recovery, injection recovery, and Bayes factors. For the frequentist analysis, we test that the cross-correlation-based optimal statistic, when modified to account for a non-negligible gravitational-wave background, accurately recovers the amplitude of the background. We also summarize recent advances and tests performed on the optimal statistic in the literature from both GWB detection and parameter estimation perspectives. The tests presented here validate current and future analyses of PTA data. △ Less

Submitted 7 July, 2023; v1 submitted 28 June, 2023; originally announced June 2023.

Comments: 30 pages, 10 figures, 1 table; Companion paper to "The NANOGrav 15-year Data Set: Evidence for a Gravitational-Wave Background"; For questions or comments, please email comments@nanograv.org

The NANOGrav 15-year Data Set: Bayesian Limits on Gravitational Waves from Individual Supermassive Black Hole Binaries

Authors: Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Paul T. Baker, Bence Bécsy, Laura Blecha, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Robin Case, J. Andrew Casey-Clyde, Maria Charisi, Shami Chatterjee, Tyler Cohen, James M. Cordes, Neil Cornish, Fronefield Crawford, H. Thankful Cromartie, Kathryn Crowter, Megan DeCesar, Paul B. Demorest, Matthew C. Digman, Timothy Dolch, Brendan Drachler , et al. (74 additional authors not shown)

Abstract: Evidence for a low-frequency stochastic gravitational wave background has recently been reported based on analyses of pulsar timing array data. The most likely source of such a background is a population of supermassive black hole binaries, the loudest of which may be individually detected in these datasets. Here we present the search for individual supermassive black hole binaries in the NANOGrav… ▽ More Evidence for a low-frequency stochastic gravitational wave background has recently been reported based on analyses of pulsar timing array data. The most likely source of such a background is a population of supermassive black hole binaries, the loudest of which may be individually detected in these datasets. Here we present the search for individual supermassive black hole binaries in the NANOGrav 15-year dataset. We introduce several new techniques, which enhance the efficiency and modeling accuracy of the analysis. The search uncovered weak evidence for two candidate signals, one with a gravitational-wave frequency of $\sim$4 nHz, and another at $\sim$170 nHz. The significance of the low-frequency candidate was greatly diminished when Hellings-Downs correlations were included in the background model. The high-frequency candidate was discounted due to the lack of a plausible host galaxy, the unlikely astrophysical prior odds of finding such a source, and since most of its support comes from a single pulsar with a commensurate binary period. Finding no compelling evidence for signals from individual binary systems, we place upper limits on the strain amplitude of gravitational waves emitted by such systems. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 23 pages, 13 figures, 2 tables. Accepted for publication in Astrophysical Journal Letters as part of Focus on NANOGrav's 15-year Data Set and the Gravitational Wave Background. For questions or comments, please email comments@nanograv.org

The NANOGrav 15-year Data Set: Search for Anisotropy in the Gravitational-Wave Background

Authors: Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Paul T. Baker, Bence Bécsy, Laura Blecha, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, J. Andrew Casey-Clyde, Maria Charisi, Shami Chatterjee, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Kathryn Crowter, Megan E. DeCesar, Paul B. Demorest, Timothy Dolch, Brendan Drachler, Elizabeth C. Ferrara, William Fiore , et al. (68 additional authors not shown)

Abstract: The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) has reported evidence for the presence of an isotropic nanohertz gravitational wave background (GWB) in its 15 yr dataset. However, if the GWB is produced by a population of inspiraling supermassive black hole binary (SMBHB) systems, then the background is predicted to be anisotropic, depending on the distribution of these… ▽ More The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) has reported evidence for the presence of an isotropic nanohertz gravitational wave background (GWB) in its 15 yr dataset. However, if the GWB is produced by a population of inspiraling supermassive black hole binary (SMBHB) systems, then the background is predicted to be anisotropic, depending on the distribution of these systems in the local Universe and the statistical properties of the SMBHB population. In this work, we search for anisotropy in the GWB using multiple methods and bases to describe the distribution of the GWB power on the sky. We do not find significant evidence of anisotropy, and place a Bayesian $95\%$ upper limit on the level of broadband anisotropy such that $(C_{l>0} / C_{l=0}) < 20\%$. We also derive conservative estimates on the anisotropy expected from a random distribution of SMBHB systems using astrophysical simulations conditioned on the isotropic GWB inferred in the 15-yr dataset, and show that this dataset has sufficient sensitivity to probe a large fraction of the predicted level of anisotropy. We end by highlighting the opportunities and challenges in searching for anisotropy in pulsar timing array data. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 19 pages, 11 figures; submitted to Astrophysical Journal Letters as part of Focus on NANOGrav's 15-year Data Set and the Gravitational Wave Background. For questions or comments, please email comments@nanograv.org

The NANOGrav 15-year Data Set: Constraints on Supermassive Black Hole Binaries from the Gravitational Wave Background

Authors: Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Paul T. Baker, Bence Bécsy, Laura Blecha, Alexander Bonilla, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Rand Burnette, Robin Case, J. Andrew Casey-Clyde, Maria Charisi, Shami Chatterjee, Katerina Chatziioannou, Belinda D. Cheeseboro, Siyuan Chen, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Kathryn Crowter, Curt J. Cutler , et al. (89 additional authors not shown)

Abstract: The NANOGrav 15-year data set shows evidence for the presence of a low-frequency gravitational-wave background (GWB). While many physical processes can source such low-frequency gravitational waves, here we analyze the signal as coming from a population of supermassive black hole (SMBH) binaries distributed throughout the Universe. We show that astrophysically motivated models of SMBH binary popul… ▽ More The NANOGrav 15-year data set shows evidence for the presence of a low-frequency gravitational-wave background (GWB). While many physical processes can source such low-frequency gravitational waves, here we analyze the signal as coming from a population of supermassive black hole (SMBH) binaries distributed throughout the Universe. We show that astrophysically motivated models of SMBH binary populations are able to reproduce both the amplitude and shape of the observed low-frequency gravitational-wave spectrum. While multiple model variations are able to reproduce the GWB spectrum at our current measurement precision, our results highlight the importance of accurately modeling binary evolution for producing realistic GWB spectra. Additionally, while reasonable parameters are able to reproduce the 15-year observations, the implied GWB amplitude necessitates either a large number of parameters to be at the edges of expected values, or a small number of parameters to be notably different from standard expectations. While we are not yet able to definitively establish the origin of the inferred GWB signal, the consistency of the signal with astrophysical expectations offers a tantalizing prospect for confirming that SMBH binaries are able to form, reach sub-parsec separations, and eventually coalesce. As the significance grows over time, higher-order features of the GWB spectrum will definitively determine the nature of the GWB and allow for novel constraints on SMBH populations. △ Less

Submitted 18 July, 2023; v1 submitted 28 June, 2023; originally announced June 2023.

Comments: Accepted by Astrophysical Journal Letters as part of Focus on NANOGrav's 15-year Data Set and the Gravitational Wave Background. For questions or comments, please email comments@nanograv.org. Edited to fix two equation typos (Eq.13 & 21), and minor text typos

The NANOGrav 15-year Data Set: Search for Signals from New Physics

Authors: Adeela Afzal, Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Paul T. Baker, Bence Bécsy, Jose Juan Blanco-Pillado, Laura Blecha, Kimberly K. Boddy, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Rand Burnette, Robin Case, Maria Charisi, Shami Chatterjee, Katerina Chatziioannou, Belinda D. Cheeseboro, Siyuan Chen, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie , et al. (98 additional authors not shown)

Abstract: The 15-year pulsar timing data set collected by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) shows positive evidence for the presence of a low-frequency gravitational-wave (GW) background. In this paper, we investigate potential cosmological interpretations of this signal, specifically cosmic inflation, scalar-induced GWs, first-order phase transitions, cosmic string… ▽ More The 15-year pulsar timing data set collected by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) shows positive evidence for the presence of a low-frequency gravitational-wave (GW) background. In this paper, we investigate potential cosmological interpretations of this signal, specifically cosmic inflation, scalar-induced GWs, first-order phase transitions, cosmic strings, and domain walls. We find that, with the exception of stable cosmic strings of field theory origin, all these models can reproduce the observed signal. When compared to the standard interpretation in terms of inspiraling supermassive black hole binaries (SMBHBs), many cosmological models seem to provide a better fit resulting in Bayes factors in the range from 10 to 100. However, these results strongly depend on modeling assumptions about the cosmic SMBHB population and, at this stage, should not be regarded as evidence for new physics. Furthermore, we identify excluded parameter regions where the predicted GW signal from cosmological sources significantly exceeds the NANOGrav signal. These parameter constraints are independent of the origin of the NANOGrav signal and illustrate how pulsar timing data provide a new way to constrain the parameter space of these models. Finally, we search for deterministic signals produced by models of ultralight dark matter (ULDM) and dark matter substructures in the Milky Way. We find no evidence for either of these signals and thus report updated constraints on these models. In the case of ULDM, these constraints outperform torsion balance and atomic clock constraints for ULDM coupled to electrons, muons, or gluons. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 74 pages, 31 figures, 4 tables; published in Astrophysical Journal Letters as part of Focus on NANOGrav's 15-year Data Set and the Gravitational Wave Background. For questions or comments, please email comments@nanograv.org

The NANOGrav 15-Year Data Set: Detector Characterization and Noise Budget

Authors: Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Paul T. Baker, Bence Bécsy, Laura Blecha, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Maria Charisi, Shami Chatterjee, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Kathryn Crowter, Megan E. Decesar, Paul B. Demorest, Timothy Dolch, Brendan Drachler, Elizabeth C. Ferrara, William Fiore, Emmanuel Fonseca , et al. (66 additional authors not shown)

Abstract: Pulsar timing arrays (PTAs) are galactic-scale gravitational wave detectors. Each individual arm, composed of a millisecond pulsar, a radio telescope, and a kiloparsecs-long path, differs in its properties but, in aggregate, can be used to extract low-frequency gravitational wave (GW) signals. We present a noise and sensitivity analysis to accompany the NANOGrav 15-year data release and associated… ▽ More Pulsar timing arrays (PTAs) are galactic-scale gravitational wave detectors. Each individual arm, composed of a millisecond pulsar, a radio telescope, and a kiloparsecs-long path, differs in its properties but, in aggregate, can be used to extract low-frequency gravitational wave (GW) signals. We present a noise and sensitivity analysis to accompany the NANOGrav 15-year data release and associated papers, along with an in-depth introduction to PTA noise models. As a first step in our analysis, we characterize each individual pulsar data set with three types of white noise parameters and two red noise parameters. These parameters, along with the timing model and, particularly, a piecewise-constant model for the time-variable dispersion measure, determine the sensitivity curve over the low-frequency GW band we are searching. We tabulate information for all of the pulsars in this data release and present some representative sensitivity curves. We then combine the individual pulsar sensitivities using a signal-to-noise-ratio statistic to calculate the global sensitivity of the PTA to a stochastic background of GWs, obtaining a minimum noise characteristic strain of $7\times 10^{-15}$ at 5 nHz. A power law-integrated analysis shows rough agreement with the amplitudes recovered in NANOGrav's 15-year GW background analysis. While our phenomenological noise model does not model all known physical effects explicitly, it provides an accurate characterization of the noise in the data while preserving sensitivity to multiple classes of GW signals. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 67 pages, 73 figures, 3 tables; published in Astrophysical Journal Letters as part of Focus on NANOGrav's 15-year Data Set and the Gravitational Wave Background. For questions or comments, please email comments@nanograv.org

The NANOGrav 15-year Data Set: Observations and Timing of 68 Millisecond Pulsars

Authors: Gabriella Agazie, Md Faisal Alam, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Paul T. Baker, Laura Blecha, Victoria Bonidie, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Bence Bécsy, Christopher Chapman, Maria Charisi, Shami Chatterjee, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Kathryn Crowter, Megan E. DeCesar, Paul B. Demorest, Timothy Dolch, Brendan Drachler , et al. (75 additional authors not shown)

Abstract: We present observations and timing analyses of 68 millisecond pulsars (MSPs) comprising the 15-year data set of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). NANOGrav is a pulsar timing array (PTA) experiment that is sensitive to low-frequency gravitational waves. This is NANOGrav's fifth public data release, including both "narrowband" and "wideband" time-of-arrival… ▽ More We present observations and timing analyses of 68 millisecond pulsars (MSPs) comprising the 15-year data set of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). NANOGrav is a pulsar timing array (PTA) experiment that is sensitive to low-frequency gravitational waves. This is NANOGrav's fifth public data release, including both "narrowband" and "wideband" time-of-arrival (TOA) measurements and corresponding pulsar timing models. We have added 21 MSPs and extended our timing baselines by three years, now spanning nearly 16 years for some of our sources. The data were collected using the Arecibo Observatory, the Green Bank Telescope, and the Very Large Array between frequencies of 327 MHz and 3 GHz, with most sources observed approximately monthly. A number of notable methodological and procedural changes were made compared to our previous data sets. These improve the overall quality of the TOA data set and are part of the transition to new pulsar timing and PTA analysis software packages. For the first time, our data products are accompanied by a full suite of software to reproduce data reduction, analysis, and results. Our timing models include a variety of newly detected astrometric and binary pulsar parameters, including several significant improvements to pulsar mass constraints. We find that the time series of 23 pulsars contain detectable levels of red noise, 10 of which are new measurements. In this data set, we find evidence for a stochastic gravitational-wave background. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 90 pages, 74 figures, 6 tables; published in Astrophysical Journal Letters as part of Focus on NANOGrav's 15-year Data Set and the Gravitational Wave Background. For questions or comments, please email comments@nanograv.org

The NANOGrav 15-year Data Set: Evidence for a Gravitational-Wave Background

Authors: Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Paul T. Baker, Bence Becsy, Laura Blecha, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Rand Burnette, Robin Case, Maria Charisi, Shami Chatterjee, Katerina Chatziioannou, Belinda D. Cheeseboro, Siyuan Chen, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Kathryn Crowter, Curt J. Cutler, Megan E. DeCesar , et al. (89 additional authors not shown)

Abstract: We report multiple lines of evidence for a stochastic signal that is correlated among 67 pulsars from the 15-year pulsar-timing data set collected by the North American Nanohertz Observatory for Gravitational Waves. The correlations follow the Hellings-Downs pattern expected for a stochastic gravitational-wave background. The presence of such a gravitational-wave background with a power-law-spectr… ▽ More We report multiple lines of evidence for a stochastic signal that is correlated among 67 pulsars from the 15-year pulsar-timing data set collected by the North American Nanohertz Observatory for Gravitational Waves. The correlations follow the Hellings-Downs pattern expected for a stochastic gravitational-wave background. The presence of such a gravitational-wave background with a power-law-spectrum is favored over a model with only independent pulsar noises with a Bayes factor in excess of $10^{14}$, and this same model is favored over an uncorrelated common power-law-spectrum model with Bayes factors of 200-1000, depending on spectral modeling choices. We have built a statistical background distribution for these latter Bayes factors using a method that removes inter-pulsar correlations from our data set, finding $p = 10^{-3}$ (approx. $3σ$) for the observed Bayes factors in the null no-correlation scenario. A frequentist test statistic built directly as a weighted sum of inter-pulsar correlations yields $p = 5 \times 10^{-5} - 1.9 \times 10^{-4}$ (approx. $3.5 - 4σ$). Assuming a fiducial $f^{-2/3}$ characteristic-strain spectrum, as appropriate for an ensemble of binary supermassive black-hole inspirals, the strain amplitude is $2.4^{+0.7}_{-0.6} \times 10^{-15}$ (median + 90% credible interval) at a reference frequency of 1/(1 yr). The inferred gravitational-wave background amplitude and spectrum are consistent with astrophysical expectations for a signal from a population of supermassive black-hole binaries, although more exotic cosmological and astrophysical sources cannot be excluded. The observation of Hellings-Downs correlations points to the gravitational-wave origin of this signal. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 30 pages, 18 figures. Published in Astrophysical Journal Letters as part of Focus on NANOGrav's 15-year Data Set and the Gravitational Wave Background. For questions or comments, please email comments@nanograv.org

Superradiance Properties of Light Black Holes and $10^{-12}$-$10^{21}$ eV Bosons

Authors: Caner Unal

Abstract: We derive the minimum spin value for the light black holes, $10^{-33}-1 \; M_\odot$ ($1-10^{33} \; {\rm gr}$), to experience superradiance via scalar, vector and tensor perturbations, corresponding to boson mass range $10^{-12}-10^{21}$ eV. We find that superradiance instability can happen even for very low spin values, ${\tilde a} \sim 10^{-5}-10^{-2}$. Since light black holes (BHs) are very unst… ▽ More We derive the minimum spin value for the light black holes, $10^{-33}-1 \; M_\odot$ ($1-10^{33} \; {\rm gr}$), to experience superradiance via scalar, vector and tensor perturbations, corresponding to boson mass range $10^{-12}-10^{21}$ eV. We find that superradiance instability can happen even for very low spin values, ${\tilde a} \sim 10^{-5}-10^{-2}$. Since light black holes (BHs) are very unstable to these perturbations and sensitive probes of bosonic particles, a single moderately spinning BH can probe/cover 2-5 orders of magnitude scalar (axion), vector (dark photon and/or photon with effective mass) and spin-2 mass. If superradiance exists, this drives the spin of the BH to almost zero immediately, independent of the BH formation mechanism. In the case that superradiance is not observed due to self-interactions, we find limits on the axion decay constant and energy density. We finally discuss briefly the implications on the Standard Model bosons and Higgs self-interaction. △ Less

Submitted 14 March, 2023; v1 submitted 19 January, 2023; originally announced January 2023.

Comments: 5 pages, 5 figures. Evaporating BHs included, BH mass range extended : 1gr to e33 gr=Msun, boson mass range extended : e-12 to e21 eV, refs added

Probing ultralight scalar, vector and tensor dark matter with pulsar timing arrays

Authors: Caner Unal, Federico R. Urban, Ely D. Kovetz

Abstract: Pulsar timing arrays (PTAs) are sensitive to oscillations in the gravitational potential along the line-of-sight due to ultralight particle pressure. We calculate the probing power of PTAs for ultralight bosons across all frequencies, from those larger than the inverse observation time to those smaller than the inverse distance to the pulsar. We show that since the signal amplitude grows comparabl… ▽ More Pulsar timing arrays (PTAs) are sensitive to oscillations in the gravitational potential along the line-of-sight due to ultralight particle pressure. We calculate the probing power of PTAs for ultralight bosons across all frequencies, from those larger than the inverse observation time to those smaller than the inverse distance to the pulsar. We show that since the signal amplitude grows comparably to the degradation in PTA sensitivity at frequencies smaller than inverse observation time, the discovery potential can be extended towards lower masses by over three decades, maintaining high precision. We demonstrate that, in the mass range $10^{-26} -10^{-23}$ eV, existing 15-year PTA data can robustly detect or rule out an ultralight component down to $O(1 - 10)\%$ of the total dark matter. Non-detection, together with other bounds in different mass ranges, will imply that ultralight scalar/axion can comprise at most $1-10\%$ of dark matter in the $10^{-30}\!-\!10^{-17}$ eV range. With 30 years of observation, current PTAs can extend the reach down to $0.1-1 \%$, while next-generation PTAs such as SKA can attain the $0.01-0.1\%$ precision. We generalize the analysis and derive predictions for ultralight spin-1 vector (i.e. dark photon) and spin-2 tensor dark components. △ Less

Submitted 20 June, 2024; v1 submitted 6 September, 2022; originally announced September 2022.

Comments: 6 pages, 3 figures, accepted by PLB

Journal ref: Physics Letters B, Volume 855, id.138830 (2024)

Unveiling cosmological information on small scales with line intensity mapping

Authors: Sarah Libanore, Caner Unal, Debanjan Sarkar, Ely D. Kovetz

Abstract: One of the toughest challenges in modern cosmology is to probe the small scales $k \sim 0.5\,$Mpc$^{-1}$ in the matter power spectrum and clustering. We show that such small scales will be accessible via upcoming line-intensity mapping surveys, with carbon monoxide (CO) emission from star-forming galaxies at high redshifts as an example. While these galaxies cannot be individually detected and the… ▽ More One of the toughest challenges in modern cosmology is to probe the small scales $k \sim 0.5\,$Mpc$^{-1}$ in the matter power spectrum and clustering. We show that such small scales will be accessible via upcoming line-intensity mapping surveys, with carbon monoxide (CO) emission from star-forming galaxies at high redshifts as an example. While these galaxies cannot be individually detected and the two-point correlations of the intensity fluctuation maps is not accessible at these scales, the voxel intensity distribution (VID) of the highly non-Gaussian intensity maps is sensitive to the integrated emission from faint sources. As we describe, the main limitations are due to uncertainties related with the halo mass function and the CO emission model. We show that via the VID, high-redshift next-generation experiments can probe deviations from $Λ$CDM of order unity, while stage-3 experiments will be able to probe deviations as small as $\lesssim10\%$ at least up to $k\sim 10\,$Mpc$^{-1}$. △ Less

Submitted 20 December, 2022; v1 submitted 2 August, 2022; originally announced August 2022.

Comments: 10 pages, 5 figures. Version accepted for publication

Journal ref: Phys. Rev. D 106, 123512 (2022)

Cosmology with the Laser Interferometer Space Antenna

Authors: Pierre Auclair, David Bacon, Tessa Baker, Tiago Barreiro, Nicola Bartolo, Enis Belgacem, Nicola Bellomo, Ido Ben-Dayan, Daniele Bertacca, Marc Besancon, Jose J. Blanco-Pillado, Diego Blas, Guillaume Boileau, Gianluca Calcagni, Robert Caldwell, Chiara Caprini, Carmelita Carbone, Chia-Feng Chang, Hsin-Yu Chen, Nelson Christensen, Sebastien Clesse, Denis Comelli, Giuseppe Congedo, Carlo Contaldi, Marco Crisostomi , et al. (155 additional authors not shown)

Abstract: The Laser Interferometer Space Antenna (LISA) has two scientific objectives of cosmological focus: to probe the expansion rate of the universe, and to understand stochastic gravitational-wave backgrounds and their implications for early universe and particle physics, from the MeV to the Planck scale. However, the range of potential cosmological applications of gravitational wave observations exten… ▽ More The Laser Interferometer Space Antenna (LISA) has two scientific objectives of cosmological focus: to probe the expansion rate of the universe, and to understand stochastic gravitational-wave backgrounds and their implications for early universe and particle physics, from the MeV to the Planck scale. However, the range of potential cosmological applications of gravitational wave observations extends well beyond these two objectives. This publication presents a summary of the state of the art in LISA cosmology, theory and methods, and identifies new opportunities to use gravitational wave observations by LISA to probe the universe. △ Less

Submitted 11 April, 2022; originally announced April 2022.

Report number: LISA CosWG-22-03

Astrophysics with the Laser Interferometer Space Antenna

Authors: Pau Amaro Seoane, Jeff Andrews, Manuel Arca Sedda, Abbas Askar, Quentin Baghi, Razvan Balasov, Imre Bartos, Simone S. Bavera, Jillian Bellovary, Christopher P. L. Berry, Emanuele Berti, Stefano Bianchi, Laura Blecha, Stephane Blondin, Tamara Bogdanović, Samuel Boissier, Matteo Bonetti, Silvia Bonoli, Elisa Bortolas, Katelyn Breivik, Pedro R. Capelo, Laurentiu Caramete, Federico Cattorini, Maria Charisi, Sylvain Chaty , et al. (134 additional authors not shown)

Abstract: The Laser Interferometer Space Antenna (LISA) will be a transformative experiment for gravitational wave astronomy, and, as such, it will offer unique opportunities to address many key astrophysical questions in a completely novel way. The synergy with ground-based and space-born instruments in the electromagnetic domain, by enabling multi-messenger observations, will add further to the discovery… ▽ More The Laser Interferometer Space Antenna (LISA) will be a transformative experiment for gravitational wave astronomy, and, as such, it will offer unique opportunities to address many key astrophysical questions in a completely novel way. The synergy with ground-based and space-born instruments in the electromagnetic domain, by enabling multi-messenger observations, will add further to the discovery potential of LISA. The next decade is crucial to prepare the astrophysical community for LISA's first observations. This review outlines the extensive landscape of astrophysical theory, numerical simulations, and astronomical observations that are instrumental for modeling and interpreting the upcoming LISA datastream. To this aim, the current knowledge in three main source classes for LISA is reviewed; ultracompact stellar-mass binaries, massive black hole binaries, and extreme or intermediate mass ratio inspirals. The relevant astrophysical processes and the established modeling techniques are summarized. Likewise, open issues and gaps in our understanding of these sources are highlighted, along with an indication of how LISA could help making progress in the different areas. New research avenues that LISA itself, or its joint exploitation with upcoming studies in the electromagnetic domain, will enable, are also illustrated. Improvements in modeling and analysis approaches, such as the combination of numerical simulations and modern data science techniques, are discussed. This review is intended to be a starting point for using LISA as a new discovery tool for understanding our Universe. △ Less

Submitted 25 May, 2023; v1 submitted 11 March, 2022; originally announced March 2022.

Journal ref: Living Reviews in Relativity, Volume 26, Article number: 2 (2023)

Properties of Ultralight Bosons from Spins of Heavy Quasars via Superradiance

Authors: Caner Unal, Fabio Pacucci, Abraham Loeb

Abstract: The mass and the spin of accreting and jetted black holes, at the center of Active Galactic Nuclei (AGNs), can be probed by analyzing their electromagnetic spectra. For this purpose, we use the Spin-Modified Fundamental Plane of black hole activity, which non-linearly connects the following four variables (in the source frame): radio luminosity, X-ray or optical luminosity (via the [OIII] emission… ▽ More The mass and the spin of accreting and jetted black holes, at the center of Active Galactic Nuclei (AGNs), can be probed by analyzing their electromagnetic spectra. For this purpose, we use the Spin-Modified Fundamental Plane of black hole activity, which non-linearly connects the following four variables (in the source frame): radio luminosity, X-ray or optical luminosity (via the [OIII] emission line), black hole mass and spin. Taking into account the uncertainties in luminosity measurements, conversion factors, relativistic beaming and physical properties of the AGN system, we derive lower bounds on the spins of a group of heavy, jetted AGNs. Using these results, we study the direct implications on the mass spectrum of the ultra-light particles of scalar (axion-like), vector (dark photon) and tensor types (additional spin-2 particles). We close unexplored gap in the parameter space $10^{-20}-10^{-19}$eV. We obtain upper bounds on the axion decay constant (equivalently lower bounds on the self-interaction strength) considering self-interactions could prevent the axion particles entering the instability, and be the reason for non-observation of superradiance. Assuming axion is described by mass and decay constant, we obtain upper limits on what fraction of dark matter can be formed by ultra-light particles and find that single spieces axion-like light particle can constitute at most $10\%$ of the dark matter in the mass range: $ 10^{-21} < μ\, (\mathrm{eV}) < 10^{-17}$. △ Less

Submitted 19 April, 2021; v1 submitted 23 December, 2020; originally announced December 2020.

Comments: 15 pages, 7 figures, version accepted for publication to JCAP

Journal ref: JCAP 05, 007 (2021)

Challenges and Opportunities of Gravitational Wave Searches at MHz to GHz Frequencies

Authors: N. Aggarwal, O. D. Aguiar, A. Bauswein, G. Cella, S. Clesse, A. M. Cruise, V. Domcke, D. G. Figueroa, A. Geraci, M. Goryachev, H. Grote, M. Hindmarsh, F. Muia, N. Mukund, D. Ottaway, M. Peloso, F. Quevedo, A. Ricciardone, J. Steinlechner, S. Steinlechner, S. Sun, M. E. Tobar, F. Torrenti, C. Unal, G. White

Abstract: The first direct measurement of gravitational waves by the LIGO and Virgo collaborations has opened up new avenues to explore our Universe. This white paper outlines the challenges and gains expected in gravitational wave searches at frequencies above the LIGO/Virgo band, with a particular focus on Ultra High-Frequency Gravitational Waves (UHF-GWs), covering the MHz to GHz range. The absence of kn… ▽ More The first direct measurement of gravitational waves by the LIGO and Virgo collaborations has opened up new avenues to explore our Universe. This white paper outlines the challenges and gains expected in gravitational wave searches at frequencies above the LIGO/Virgo band, with a particular focus on Ultra High-Frequency Gravitational Waves (UHF-GWs), covering the MHz to GHz range. The absence of known astrophysical sources in this frequency range provides a unique opportunity to discover physics beyond the Standard Model operating both in the early and late Universe, and we highlight some of the most promising gravitational sources. We review several detector concepts which have been proposed to take up this challenge, and compare their expected sensitivity with the signal strength predicted in various models. This report is the summary of the workshop "Challenges and opportunities of high-frequency gravitational wave detection" held at ICTP Trieste, Italy in October 2019, that set up the stage for the recently launched Ultra-High-Frequency Gravitational Wave (UHF-GW) initiative. △ Less

Submitted 13 December, 2021; v1 submitted 24 November, 2020; originally announced November 2020.

Comments: Published in Living Reviews in Relativity

Report number: CERN-TH-2020-185, HIP-2020-28/TH, DESY 20-195

Journal ref: Living Reviews in Relativity volume 24, Article number: 4 (2021)

Multi-messenger Probes of Inflationary Fluctuations and Primordial Black Holes

Authors: Caner Unal, Ely D. Kovetz, Subodh P. Patil

Abstract: Next generation cosmic microwave background spectral distortion and pulsar timing array experiments have the potential to probe primordial fluctuations at small scales with remarkable sensitivity. We demonstrate the potential of these probes to either detect signatures of primordial black holes (PBHs) sourced from primordial overdensities within the standard thermal history of the universe over a… ▽ More Next generation cosmic microwave background spectral distortion and pulsar timing array experiments have the potential to probe primordial fluctuations at small scales with remarkable sensitivity. We demonstrate the potential of these probes to either detect signatures of primordial black holes (PBHs) sourced from primordial overdensities within the standard thermal history of the universe over a 13-decade mass range ${\cal O}(0.1-10^{12})M_\odot$, or constrain their existence to a negligible abundance. Our conclusions are based only on global cosmological signals, and are robust under changes in i) the statistical properties of the primordial density fluctuations (whether Gaussian or non-Gaussian), ii) the merger and accretion history of the PBHs and assumptions about associated astrophysical processes, and iii) clustering statistics. Any positive detection of enhanced primordial fluctuations at small scales would have far-reaching implications from the content of dark matter to origin of BHs in the centers of galaxies, and to the field content of the inflation. On the other hand, their non-detection would also have important corollaries. For example, non-detection up to forecast sensitivities would tell us that PBHs larger than a fraction of a solar mass can constitute no more than a negligible fraction of dark matter. Moreover, non-detection will also rule out the scenario that PBHs generated by primordial overdensities could be the progenitors of super-massive black holes (SMBHs), of topical interest as there are only a few widely accepted proposals for the formation of SMBHs, an even more pressing question after the detection of active galactic nuclei over a billion solar masses at redshifts $z \geq 7$. Finally, non-detection sets the strongest bounds on the amplitude of small scale inflationary fluctuations for over 6 decades. △ Less

Submitted 1 March, 2021; v1 submitted 25 August, 2020; originally announced August 2020.

Comments: 9 pages, 2 figures, accepted version

Journal ref: Phys. Rev. D 103, no.6, 063519 (2021)

On Spin Dependence of the Fundamental Plane of Black Hole Activity

Authors: Caner Unal, Abraham Loeb

Abstract: The Fundamental Plane (FP) of Black Hole (BH) Activity in galactic nuclei relates X-ray and radio luminosities to BH mass and accretion rate. However, there is a large scatter exhibited by the data, which motivated us for a new variable. We add BH spin as a new variable and estimate the spin dependence of the jet power and disk luminosity in terms of radio and X-ray luminosities. We assume the Bla… ▽ More The Fundamental Plane (FP) of Black Hole (BH) Activity in galactic nuclei relates X-ray and radio luminosities to BH mass and accretion rate. However, there is a large scatter exhibited by the data, which motivated us for a new variable. We add BH spin as a new variable and estimate the spin dependence of the jet power and disk luminosity in terms of radio and X-ray luminosities. We assume the Blandford-Znajek process as the main source of the outflow, and find that the jet power depends on BH spin stronger than quadratically at moderate and large spin values. We perform a statistical analysis for 10 AGNs which have sub-Eddington accretion rates and whose spin values are measured independently via the reflection or continuum-fitting methods, and find that the spin-dependent relation describes the data significantly better. This analysis, if supported with more data, could imply not only the spin dependence of the FP relation, but also the Blandford-Znajek process in AGN jets. △ Less

Submitted 21 April, 2020; v1 submitted 26 February, 2020; originally announced February 2020.

Comments: Accepted for publication in MNRAS

Nonlinear perturbations from axion-gauge fields dynamics during inflation

Authors: Alexandros Papageorgiou, Marco Peloso, Caner Unal

Abstract: We study a variant of the Chromo-Natural Inflation (CNI) mechanism in which the inflaton interacts only gravitationally with the CNI fields. Integrating out all the non-dynamical scalar fields of the model results in a coupling between the perturbations of the inflaton and of the CNI pseudo-scalar which is significantly greater than the one obtained in the absence of the gauge CNI dynamics. We com… ▽ More We study a variant of the Chromo-Natural Inflation (CNI) mechanism in which the inflaton interacts only gravitationally with the CNI fields. Integrating out all the non-dynamical scalar fields of the model results in a coupling between the perturbations of the inflaton and of the CNI pseudo-scalar which is significantly greater than the one obtained in the absence of the gauge CNI dynamics. We compute how this greater coupling impacts the power spectrum of the inflaton perturbations that are sourced nonlinearly by the unstable (tensor) gauge CNI modes, and we require that the amplitude of these modes is well below that of the linear perturbations. Combining this result with various constraints, including backreaction effects, the requirement of having observable and dominant sourced gravitational waves (GW), and the current upper bound on the tensor-to-scalar ratio, significantly constrains the range of parameter space where this model can produce an interesting GW signal. △ Less

Submitted 3 September, 2019; v1 submitted 2 April, 2019; originally announced April 2019.

Comments: 42 pages, 8 figures

The Ineludible non-Gaussianity of the Primordial Black Hole Abundance

Authors: V. De Luca, G. Franciolini, A. Kehagias, M. Peloso, A. Riotto, C. Ünal

Abstract: We study the formation of primordial black holes when they are generated by the collapse of large overdensities in the early universe. Since the density contrast is related to the comoving curvature perturbation by a nonlinear relation, the overdensity statistics is unavoidably non-Gaussian. We show that the abundance of primordial black holes at formation may not be captured by a perturbative app… ▽ More We study the formation of primordial black holes when they are generated by the collapse of large overdensities in the early universe. Since the density contrast is related to the comoving curvature perturbation by a nonlinear relation, the overdensity statistics is unavoidably non-Gaussian. We show that the abundance of primordial black holes at formation may not be captured by a perturbative approach which retains the first few cumulants of the non-Gaussian probability distribution. We provide two techniques to calculate the non-Gaussian abundance of primordial black holes at formation, one based on peak theory and the other on threshold statistics. Our results show that the unavoidable non-Gaussian nature of the inhomogeneities in the energy density makes it harder to generate PBHs. We provide simple (semi-)analytical expressions to calculate the non-Gaussian abundances of the primordial black holes and show that for both narrow and broad power spectra the gaussian case from threshold statistics is reproduced by increasing the amplitude of the power spectrum by a factor ${\cal O}(2÷3)$. △ Less

Submitted 30 September, 2019; v1 submitted 1 April, 2019; originally announced April 2019.

Comments: 26 pages, 10 figures, matching published version

Electromagnetic probes of primordial black holes as dark matter

Authors: Y. Ali-Haimoud, S. Clesse, J. Garcia-Bellido, A. Kashlinsky, L. Wyrzykowski, A. Achucarro, L. Amendola, J. Annis, A. Arbey, R. G. Arendt, F. Atrio-Barandela, N. Bellomo, K. Belotsky, J-L. Bernal, S. Bird, V. Bozza, C. Byrnes, S. Calchi Novati, F. Calore, B. J. Carr, J. Chluba, I. Cholis, A. Cieplak, P. Cole, I. Dalianis , et al. (69 additional authors not shown)

Abstract: The LIGO discoveries have rekindled suggestions that primordial black holes (BHs) may constitute part to all of the dark matter (DM) in the Universe. Such suggestions came from 1) the observed merger rate of the BHs, 2) their unusual masses, 3) their low/zero spins, and 4) also from the independently uncovered cosmic infrared background (CIB) fluctuations signal of high amplitude and coherence wit… ▽ More The LIGO discoveries have rekindled suggestions that primordial black holes (BHs) may constitute part to all of the dark matter (DM) in the Universe. Such suggestions came from 1) the observed merger rate of the BHs, 2) their unusual masses, 3) their low/zero spins, and 4) also from the independently uncovered cosmic infrared background (CIB) fluctuations signal of high amplitude and coherence with unresolved cosmic X-ray background (CXB). Here we summarize the prospects to resolve this important issue with electromagnetic observations using the instruments and tools expected in the 2020's. These prospects appear promising to make significant, and potentially critical, advances. We demonstrate that in the next decade, new space- and ground-borne electromagnetic instruments, combined with concurrent theoretical efforts, should shed critical light on the long-considered link between primordial BHs and DM. Specifically the new data and methodologies under this program will involve: I) Probing with high precision the spatial spectrum of source-subtracted CIB with Euclid and WFIRST, and its coherence with unresolved cosmic X-ray background using eROSITA and Athena, II) Advanced searches for microlensing of Galactic stars by the intervening Galactic Halo BHs with OGLE, Gaia, LSST and WFIRST, III) Supernovae (SNe) lensing in the upcoming surveys with WFIRST, LSST and also potentially with Euclid and JWST, IV) Advanced theoretical work to understand the details of PBH accretion and evolution and their influence on cosmic microwave background (CMB) anisotropies in light of the next generation CMB experiments, V) Better new samples and theoretical understanding involving stability and properties of ultra faint dwarf galaxies, pulsar timing, and cosmological quasar lensing. △ Less

Submitted 12 March, 2019; v1 submitted 11 March, 2019; originally announced March 2019.

Comments: Science whitepaper submitted to the Astro2020 Decadal Survey

Imprints of Primordial Non-Gaussianity on Gravitational Wave Spectrum

Authors: Caner Unal

Abstract: Although Cosmic Microwave Background and Large Scale Structure probe the largest scales of our universe with ever increasing precision, our knowledge about the smaller scales is still very limited other than the bounds on Primordial Black Holes. We show that the statistical properties of the small scale quantum fluctuations can be probed via the stochastic gravitational wave background, which is i… ▽ More Although Cosmic Microwave Background and Large Scale Structure probe the largest scales of our universe with ever increasing precision, our knowledge about the smaller scales is still very limited other than the bounds on Primordial Black Holes. We show that the statistical properties of the small scale quantum fluctuations can be probed via the stochastic gravitational wave background, which is induced as the scalar modes re-enter the horizon. We found that even if scalar curvature fluctuations have a subdominant non-Gaussian component, these non-Gaussian perturbations can source a dominant portion of the induced GWs. Moreover, the GWs sourced by non-Gaussian scalar fluctuations peaks at a higher frequency and this can result in distinctive observational signatures. We found that the sensitive next-generation-interferometers, which will/could reach $Ω_{GW}h^2 \sim 10^{-15}$ (such as PTA-SKA, LISA, DECIGO, BBO, CE, ET), can probe $f_{NL} \sim 0.5$ which is even better than the predictions of the next generation CMB experiments. If the induced GW background is detected, but not the signatures arising from the non-Gaussian component, $ζ= ζ_G + f_{\rm NL} \, ζ_G^{2}$, this translates into bounds on $f_{\rm NL}$ depending on the amplitude and the width of the GW signal. If the induced GW background is not detected at all, this translates into bounds on scalar fluctuations. The results are independent from the fact that whether PBH are DM or completely negligible part of the current energy density. △ Less

Submitted 7 January, 2020; v1 submitted 22 November, 2018; originally announced November 2018.

Comments: 7 pages, 10 figures, added figures of merit to stress that interferometers can probe nonGaussianity sensitively (fnl~ 0.5)

Journal ref: Phys. Rev. D 99, 041301 (2019)

Nonlinear perturbations from the coupling of the inflaton to a non-Abelian gauge field, with a focus on Chromo-Natural Inflation

Authors: Alexandros Papageorgiou, Marco Peloso, Caner Unal

Abstract: Several models of inflation employing a triplet of SU(2) vectors with spatially orthogonal vacuum expectation values (VEVs) have been recently proposed. One (tensor) combination $t$ of the vector modes is amplified in some momentum range during inflation. Due to the vector VEVs, this combination mixes with gravitational waves (GW) at the linear level, resulting in a GW amplification that has been… ▽ More Several models of inflation employing a triplet of SU(2) vectors with spatially orthogonal vacuum expectation values (VEVs) have been recently proposed. One (tensor) combination $t$ of the vector modes is amplified in some momentum range during inflation. Due to the vector VEVs, this combination mixes with gravitational waves (GW) at the linear level, resulting in a GW amplification that has been well studied in the literature. Scalar perturbations in this class of models have been so far studied only at the linear level. We perform a first step toward the nonlinear computation using as an example the original model of Chromo-Natural Inflation. We compute the contribution to the scalar power spectrum arising from the coupling of the combination $t$ to the inflaton. This contribution is mostly controlled by a single parameter of the model (namely, the ratio between the mass of the fluctuations of the vector field and the Hubble rate), and, for a wide range of this parameter, it can significantly affect the phenomenology obtained from the linear theory. This nonlinear contribution is significantly blue, improving the comparison between the two-point function and the Cosmic Microwave Background (CMB) data. This growth can be also relevant for smaller scale phenomenology, such as large scale structure, CMB distortions, and primordial black holes. △ Less

Submitted 13 September, 2018; v1 submitted 21 June, 2018; originally announced June 2018.

Comments: 30 pages, 8 figures, correction to loop integral due to extra term, main results unchanged, match with the JCAP version

Report number: UMN-TH 3723/18

Gravitational Wave signatures of inflationary models from Primordial Black Hole Dark Matter

Authors: Juan Garcia-Bellido, Marco Peloso, Caner Unal

Abstract: Primordial Black Holes (PBH) could be the cold dark matter of the universe. They could have arisen from large (order one) curvature fluctuations produced during inflation that reentered the horizon in the radiation era. At reentry, these fluctuations source gravitational waves (GW) via second order anisotropic stresses. These GW, together with those (possibly) sourced during inflation by the same… ▽ More Primordial Black Holes (PBH) could be the cold dark matter of the universe. They could have arisen from large (order one) curvature fluctuations produced during inflation that reentered the horizon in the radiation era. At reentry, these fluctuations source gravitational waves (GW) via second order anisotropic stresses. These GW, together with those (possibly) sourced during inflation by the same mechanism responsible for the large curvature fluctuations, constitute a primordial stochastic GW background (SGWB) that unavoidably accompanies the PBH formation. We study how the amplitude and the range of frequencies of this signal depend on the statistics (Gaussian versus $χ^2$) of the primordial curvature fluctuations, and on the evolution of the PBH mass function due to accretion and merging. We then compare this signal with the sensitivity of present and future detectors, at PTA and LISA scales. We find that this SGWB will help to probe, or strongly constrain, the early universe mechanism of PBH production. The comparison between the peak mass of the PBH distribution and the peak frequency of this SGWB will provide important information on the merging and accretion evolution of the PBH mass distribution from their formation to the present era. Different assumptions on the statistics and on the PBH evolution also result in different amounts of CMB $μ$-distortions. Therefore the above results can be complemented by the detection (or the absence) of $μ$-distortions with an experiment such as PIXIE. △ Less

Submitted 8 July, 2017; originally announced July 2017.

Comments: 32 pages, 12 figures

Report number: IFT--UAM/CSIC-17-056, UMN-TH 3630/17

Gravitational waves at interferometer scales and primordial black holes in axion inflation

Authors: Juan Garcia-Bellido, Marco Peloso, Caner Unal

Abstract: We study the prospects of detection at terrestrial and space interferometers, as well as at pulsar timing array experiments, of a stochastic gravitational wave background which can be produced in models of axion inflation. This potential signal, and the development of these experiments, open a new window on inflation on scales much smaller than those currently probed with Cosmic Microwave Backgrou… ▽ More We study the prospects of detection at terrestrial and space interferometers, as well as at pulsar timing array experiments, of a stochastic gravitational wave background which can be produced in models of axion inflation. This potential signal, and the development of these experiments, open a new window on inflation on scales much smaller than those currently probed with Cosmic Microwave Background and Large Scale Structure measurements. The sourced signal generated in axion inflation is an ideal candidate for such searches, since it naturally grows at small scales, and it has specific properties (chirality and non-gaussianity) that can distinguish it from an astrophysical background. We study under which conditions such a signal can be produced at an observable level, without the simultaneous overproduction of scalar perturbations in excess of what is allowed by the primordial black hole limits. We also explore the possibility that scalar perturbations generated in a modified version of this model may provide a distribution of primordial black holes compatible with the current bounds, that can act as a seeds of the present black holes in the universe. △ Less

Submitted 20 December, 2016; v1 submitted 12 October, 2016; originally announced October 2016.

Comments: 29 pages, 10 figures. Version accepted for publication in JCAP. Added PTA bounds and few minor corrections

Report number: IFT-UAM/CSIC-16-100, UMN-TH-3607/16

Rolling axions during inflation: perturbativity and signatures

Authors: Marco Peloso, Lorenzo Sorbo, Caner Unal

Abstract: The motion of a pseudo-scalar field $X$ during inflation naturally induces a significant amplification of the gauge fields to which it is coupled. The amplified gauge fields can source characteristic scalar and tensor primordial perturbations. Several phenomenological implications have been discussed in the cases in which (i) $X$ is the inflation, and (ii) $X$ is a field different from the inflati… ▽ More The motion of a pseudo-scalar field $X$ during inflation naturally induces a significant amplification of the gauge fields to which it is coupled. The amplified gauge fields can source characteristic scalar and tensor primordial perturbations. Several phenomenological implications have been discussed in the cases in which (i) $X$ is the inflation, and (ii) $X$ is a field different from the inflation, that experiences a temporary speed up during inflation. In this second case, visible sourced gravitational waves (GW) can be produced at the CMB scales without affecting the scalar perturbations, even if the scale of inflation is several orders of magnitude below what is required to produce a visible vacuum GW signal. Perturbativity considerations can be used to limit the regime in which these results are under perturbative control. We revised limits recently claimed for the case (i), and we extend these considerations to the case (ii). We show that, in both cases, these limits are satisfied by the applications that generate signals at CMB scales. Applications that generate gravitational waves and primordial black holes at much smaller scales are at the limit of the validity of this perturbativity analysis, so we expect those results to be valid up to possibly order one corrections. △ Less

Submitted 1 June, 2016; originally announced June 2016.

Comments: 37 pages, 9 figures

Report number: ACFI-T16-16, UMN--TH--3531/16

Scale-dependent gravitational waves from a rolling axion

Authors: Ryo Namba, Marco Peloso, Maresuke Shiraishi, Lorenzo Sorbo, Caner Unal

Abstract: We consider a model in which a pseudo-scalar field $σ$ rolls for some e-folds during inflation, sourcing one helicity of a gauge field. These fields are only gravitationally coupled to the inflaton, and therefore produce scalar and tensor primordial perturbations only through gravitational interactions. These sourced signals are localized on modes that exit the horizon while the roll of $σ$ is sig… ▽ More We consider a model in which a pseudo-scalar field $σ$ rolls for some e-folds during inflation, sourcing one helicity of a gauge field. These fields are only gravitationally coupled to the inflaton, and therefore produce scalar and tensor primordial perturbations only through gravitational interactions. These sourced signals are localized on modes that exit the horizon while the roll of $σ$ is significant. We focus our study on cases in which the model can simultaneously produce (i) a large gravitational wave signal, resulting in observable B-modes of the CMB polarizations, and (ii) sufficiently small scalar perturbations, so to be in agreement with the current limits from temperature anisotropies. Different choice of parameters can instead lead to a localized and visible departure from gaussianity in the scalar sector, either at CMB or LSS scales. △ Less

Submitted 24 September, 2015; originally announced September 2015.

Comments: 41+1 pages, 29 figures

Trajectories with suppressed tensor-to-scalar ratio in Aligned Natural Inflation

Authors: Marco Peloso, Caner Unal

Abstract: In Aligned Natural Inflation, an alignment between different potential terms produces an inflaton excursion greater than the axion scales in the potential. We show that, starting from a general potential of two axions with two aligned potential terms, the effective theory for the resulting light direction is characterized by four parameters: an effective potential scale, an effective axion constan… ▽ More In Aligned Natural Inflation, an alignment between different potential terms produces an inflaton excursion greater than the axion scales in the potential. We show that, starting from a general potential of two axions with two aligned potential terms, the effective theory for the resulting light direction is characterized by four parameters: an effective potential scale, an effective axion constant, and two extra parameters (related to ratios of the axion scales and the potential scales in the $2-$field theory). For all choices of these extra parameters, the model can support inflation along valleys (in the $2-$field space) that end in minima of the potential. This leads to a phenomenology similar to that of single field Natural Inflation. For a significant range of the extra two parameters, the model possesses also higher altitude inflationary trajectories passing through saddle points of the $2-$field potential, and disconnected from any minimum. These plateaus end when the heavier direction becomes unstable, and therefore all of inflation takes place close to the saddle point, where - due to the higher altitude - the potential is flatter (smaller $ε$ parameter). As a consequence, a tensor-to-scalar ratio $r = {\rm O } \left( 10^{-4} - 10^{-2} \right)$ can be easily achieved in the allowed $n_s$ region, well within the latest $1 σ$ CMB contours. △ Less

Submitted 10 April, 2015; originally announced April 2015.

Report number: UMN--TH--3429/15

The expected anisotropy in solid inflation

Authors: Nicola Bartolo, Marco Peloso, Angelo Ricciardone, Caner Unal

Abstract: Solid inflation is an effective field theory of inflation in which isotropy and homogeneity are accomplished via a specific combination of anisotropic sources (three scalar fields that individually break isotropy). This results in specific observational signatures that are not found in standard models of inflation: a non-trivial angular dependence for the squeezed bispectrum, and a possibly long p… ▽ More Solid inflation is an effective field theory of inflation in which isotropy and homogeneity are accomplished via a specific combination of anisotropic sources (three scalar fields that individually break isotropy). This results in specific observational signatures that are not found in standard models of inflation: a non-trivial angular dependence for the squeezed bispectrum, and a possibly long period of anisotropic inflation (to drive inflation, the "solid" must be very insensitive to any deformation, and thus background anisotropies are very slowly erased). In this paper we compute the expected level of statistical anisotropy in the power spectrum of the curvature perturbations of this model. To do so, we account for the classical background values of the three scalar fields that are generated on large (superhorizon) scales during inflation via a random walk sum, as the perturbation modes leave the horizon. Such an anisotropy is unavoidably generated, even starting from perfectly isotropic classical initial conditions. The expected level of anisotropy is related to the duration of inflation and to the amplitude of the squeezed bispectrum. If this amplitude is close to its current observational limit (so that one of the most interesting predictions of the model can be observed in the near future), we find that a level of statistical anisotropy $\gtrsim 3\%$ in the power spectrum is to be expected, if inflation lasted $\gtrsim 20-30$ e-folds more than the final $50-60$ efolds required to generare the CMB modes. We also comment and point out various similarities between solid inflation and models of inflation where a suitable coupling of the inflaton to a vector kinetic term $F^{2}$ gives frozen and scale invariant vector perturbations on superhorizon scales. △ Less

Submitted 30 July, 2014; originally announced July 2014.

Comments: 12 pages, 2 figures

Report number: UMN--TH--3342/14

Journal ref: JCAP 1411 (2014) 11, 009