BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Chiara M. F. Mingarelli

Chiara M. F. Mingarelli appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
gr-qcastro-ph.GAastro-ph.HEastro-ph.COastro-ph.IMhep-phphysics.data-an

Catalog footprint

What is connected

20works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

20 published item(s)

preprint2026arXiv

The NANOGrav 15 yr Data Set: Piecewise Power-Law Reconstruction of the Gravitational-Wave Background

The NANOGrav 15-year (NG15) data set provides evidence for a gravitational-wave background (GWB) signal at nanohertz frequencies, which is expected to originate either from a cosmic population of inspiraling supermassive black-hole binaries or new particle physics in the early Universe. A firm identification of the source of the NG15 signal requires an accurate reconstruction of its frequency spectrum. In this paper, we provide such a spectral characterization of the NG15 signal based on a piecewise power-law (PPL) ansatz that strikes a balance between existing alternatives in the literature. Our PPL reconstruction is more flexible than the standard constant-power-law model, which describes the GWB spectrum in terms of only two parameters: an amplitude A and a spectral index gamma. Concurrently, it better approximates physically realistic GWB spectra -- especially those of cosmological origin -- than the free spectral model, since the latter allows for arbitrary variations in the GWB amplitude from one frequency bin to the next. Our PPL reconstruction of the NG15 signal relies on individual PPL models with a fixed number of internal nodes (i.e., constant power law, broken power law, doubly broken power law, etc.) that are ultimately combined in a Bayesian model average. The data products resulting from our analysis provide the basis for fast refits of spectral GWB models.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

UGC 4211: A Confirmed Dual Active Galactic Nucleus in the Local Universe at 230 pc Nuclear Separation

We present multi-wavelength high-spatial resolution (~0.1'', 70 pc) observations of UGC 4211 at z=0.03474, a late-stage major galaxy merger at the closest nuclear separation yet found in near-IR imaging (0.32'', ~230 pc projected separation). Using Hubble Space Telescope/STIS, VLT/MUSE+AO, Keck/OSIRIS+AO spectroscopy, and ALMA observations, we show that the spatial distribution, optical and NIR emission lines, and millimeter continuum emission are all consistent with both nuclei being powered by accreting supermassive black holes (SMBHs). Our data, combined with common black hole mass prescriptions, suggests that both SMBHs have similar masses, log MBH~8.1 (south) and log MBH~8.3 (north), respectively. The projected separation of 230 pc (~6X the black hole sphere of influence) represents the closest-separation dual AGN studied to date with multi-wavelength resolved spectroscopy and shows the potential of nuclear (<50 pc) continuum observations with ALMA to discover hidden growing SMBH pairs. While the exact occurrence rate of close-separation dual AGN is not yet known, it may be surprisingly high, given that UGC 4211 was found within a small, volume-limited sample of nearby hard X-ray detected AGN. Observations of dual SMBH binaries in the sub-kpc regime at the final stages of dynamical friction provide important constraints for future gravitational wave observatories.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Searching For Gravitational Waves From Cosmological Phase Transitions With The NANOGrav 12.5-year dataset

We search for a first-order phase transition gravitational wave signal in 45 pulsars from the NANOGrav 12.5 year dataset. We find that the data can be modeled in terms of a strong first order phase transition taking place at temperatures below the electroweak scale. However, we do not observe any strong preference for a phase-transition interpretation of the signal over the standard astrophysical interpretation in terms of supermassive black holes mergers; but we expect to gain additional discriminating power with future datasets, improving the signal to noise ratio and extending the sensitivity window to lower frequencies. An interesting open question is how well gravitational wave observatories could separate such signals.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A quasar-based supermassive black hole binary population model: implications for the gravitational-wave background

The nanohertz gravitational wave background (GWB) is believed to be dominated by GW emission from supermassive black hole binaries (SMBHBs). Observations of several dual active galactic nuclei (AGN) strongly suggest a link between AGN and SMBHBs, given that these dual AGN systems will eventually form bound binary pairs. Here we develop an exploratory SMBHB population model based on empirically constrained quasar populations, allowing us to decompose the GWB amplitude into an underlying distribution of SMBH masses, SMBHB number density, and volume enclosing the GWB. Our approach also allows us to self-consistently predict the number of local SMBHB systems from the GWB amplitude. Interestingly, we find the local number density of SMBHBs implied by the common-process signal in the NANOGrav 12.5-yr dataset to be roughly five times larger than previously predicted by other models. We also find that at most $\sim 25 \%$ of SMBHBs can be associated with quasars. Furthermore, our quasar-based approach predicts $\gtrsim 95\%$ of the GWB signal comes from $z \lesssim 2.5$, and that SMBHBs contributing to the GWB have masses $\gtrsim 10^8 M_\odot$. We also explore how different empirical galaxy-black hole scaling relations affect the local number density of GW sources, and find that relations predicting more massive black holes decrease the local number density of SMBHBs. Overall, our results point to the important role that a measurement of the GWB will play in directly constraining the cosmic population of SMBHBs, as well as their connections to quasars and galaxy mergers.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Insights into searches for anisotropies in the nanohertz gravitational-wave background

Within the next several years pulsar timing arrays (PTAs) are positioned to detect the stochastic gravitational-wave background (GWB) likely produced by the collection of inspiralling super-massive black holes binaries, and potentially constrain some exotic physics. So far most of the pulsar timing data analysis has focused on the monopole of the GWB, assuming it is perfectly isotropic. The natural next step is to search for anisotropies in the GWB. In this paper, we use the recently developed PTA Fisher matrix to gain insights into optimal search strategies for GWB anisotropies. For concreteness, we apply our results to EPTA data, using realistic noise characteristics of its pulsars. We project the detectability of a GWB whose angular dependence is assumed to be a linear combination of predetermined maps, such as spherical harmonics or coarse pixels. We find that the GWB monopole is always statistically correlated with these maps, implying a loss of sensitivity to the monopole when searching simultaneously for anisotropies. We then derive the angular distributions of the GWB intensity to which a PTA is most sensitive, and illustrate how one may use these "principal maps" to approximately reconstruct the angular dependence of the GWB. Since the principal maps are neither perfectly anisotropic nor uncorrelated with the monopole, we also develop a frequentist criterion to specifically search for anisotropies in the GWB without any prior knowledge about their angular distribution. Lastly, we show how to recover existing EPTA results with our Fisher formalism, and clarify their meaning. The tools presented here will be valuable in guiding and optimizing the computationally demanding analyses of pulsar timing data.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The NANOGrav 12.5-year Data Set: Search For An Isotropic Stochastic Gravitational-Wave Background

We search for an isotropic stochastic gravitational-wave background (GWB) in the $12.5$-year pulsar timing data set collected by the North American Nanohertz Observatory for Gravitational Waves. Our analysis finds strong evidence of a stochastic process, modeled as a power-law, with common amplitude and spectral slope across pulsars. The Bayesian posterior of the amplitude for an $f^{-2/3}$ power-law spectrum, expressed as the characteristic GW strain, has median $1.92 \times 10^{-15}$ and $5\%$--$95\%$ quantiles of $1.37$--$2.67 \times 10^{-15}$ at a reference frequency of $f_\mathrm{yr} = 1 ~\mathrm{yr}^{-1}$. The Bayes factor in favor of the common-spectrum process versus independent red-noise processes in each pulsar exceeds $10,000$. However, we find no statistically significant evidence that this process has quadrupolar spatial correlations, which we would consider necessary to claim a GWB detection consistent with general relativity. We find that the process has neither monopolar nor dipolar correlations, which may arise from, for example, reference clock or solar system ephemeris systematics, respectively. The amplitude posterior has significant support above previously reported upper limits; we explain this in terms of the Bayesian priors assumed for intrinsic pulsar red noise. We examine potential implications for the supermassive black hole binary population under the hypothesis that the signal is indeed astrophysical in nature.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Accelerated orbital decay of supermassive black hole binaries in merging nuclear star clusters

The coalescence of supermassive black holes (SMBHs) should generate the strongest sources of gravitational waves (GWs) in the Universe. However, the dynamics of their coalescence is the subject of much debate. In this study, we use a suite of $N$-body simulations to follow the merger of two nuclear star clusters (NSCs), each hosting a SMBH in their centre. We find that the presence of distinct star clusters around each SMBH has important consequences for the dynamical evolution of the SMBH binary: (i) The separation between the SMBHs decreases by a few orders of magnitude in the first few Myrs by the combined effects of dynamical friction and a drag force caused by tidally stripped stars. In fact, this is a significant speedup for equal mass ratio binaries, and becomes extreme for unequal mass ratios, e.g. 1:10 or 1:100, which traditional dynamical friction alone would not permit to bind. (ii) The subsequent binary hardening is driven by the gravitational slingshots between the SMBH binary and stars, and also depends on the mass ratio between the SMBHs. Thus, with this additional drag force, we find that all SMBHs in our suite coalesce within a Hubble time. Given that about 50% of Milky Way sized galaxies host NSCs, our results are encouraging for upcoming GW observations with the Laser Interferometer Space Antenna -- LISA -- which will detect SMBH coalescence in the $10^4-10^7$$M_{\rm \odot}$ mass range.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Constraining Galactic Structure with the LISA White Dwarf Foreground

White dwarfs comprise 95% of all stellar remnants, and are thus an excellent tracer of old stellar populations in the Milky Way. Current and planned telescopes are not able to directly probe the white dwarf population in its entirety due to its inherently low luminosity. However, the Galactic population of double white dwarf binaries gives rise to a millihertz gravitational-wave foreground detectable by the Laser Interferometer Space Antenna (LISA). Here we show how characterizing this foreground's angular power spectrum will enable us to probe the Galactic structure in a novel way and measure the vertical scale height of the Galaxy's oldest stellar populations. We do this using a binary population synthesis study that incorporates different Galactic spatial distributions for the double white dwarf population. We find that the level of anisotropy in the white dwarf foreground's angular power spectrum is strongly dependent on the vertical scale height of the population. Finally, we show that LISA can probe the vertical scale height of the Galactic double white dwarf population with an accuracy of 50 pc-200pc, depending on angular resolution limits, using the angular power spectrum of the white dwarf foreground.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Fisher formalism for anisotropic gravitational-wave background searches with pulsar timing arrays

Pulsar timing arrays (PTAs) are currently the only experiments directly sensitive to gravitational waves with decade-long periods. Within the next five to ten years, PTAs are expected to detect the stochastic gravitational-wave background (SGWB) collectively sourced by inspiralling supermassive black hole binaries. It is expected that this background is mostly isotropic, and current searches focus on the monopole part of the SGWB. Looking ahead, anisotropies in the SGWB may provide a trove of additional information both on known and unknown astrophysical and cosmological sources. In this paper, we build a simple yet realistic Fisher formalism for anisotropic SGWB searches with PTAs. Our formalism is able to accommodate realistic properties of PTAs, and allows simple and accurate forecasts. We illustrate our approach with an idealized PTA consisting of identical, isotropically distributed pulsars. In a companion paper, we apply our formalism to current PTAs and show that it can be a powerful tool to guide and optimize real data analysis.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Multi-Messenger Gravitational Wave Searches with Pulsar Timing Arrays: Application to 3C66B Using the NANOGrav 11-year Data Set

When galaxies merge, the supermassive black holes in their centers may form binaries and, during the process of merger, emit low-frequency gravitational radiation in the process. In this paper we consider the galaxy 3C66B, which was used as the target of the first multi-messenger search for gravitational waves. Due to the observed periodicities present in the photometric and astrometric data of the source of the source, it has been theorized to contain a supermassive black hole binary. Its apparent 1.05-year orbital period would place the gravitational wave emission directly in the pulsar timing band. Since the first pulsar timing array study of 3C66B, revised models of the source have been published, and timing array sensitivities and techniques have improved dramatically. With these advances, we further constrain the chirp mass of the potential supermassive black hole binary in 3C66B to less than $(1.65\pm0.02) \times 10^9~{M_\odot}$ using data from the NANOGrav 11-year data set. This upper limit provides a factor of 1.6 improvement over previous limits, and a factor of 4.3 over the first search done. Nevertheless, the most recent orbital model for the source is still consistent with our limit from pulsar timing array data. In addition, we are able to quantify the improvement made by the inclusion of source properties gleaned from electromagnetic data to `blind' pulsar timing array searches. With these methods, it is apparent that it is not necessary to obtain exact a priori knowledge of the period of a binary to gain meaningful astrophysical inferences.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Gravitational-wave cosmology across 29 decades in frequency

Quantum fluctuations of the gravitational field in the early Universe, amplified by inflation, produce a primordial gravitational-wave background across a broad frequency band. We derive constraints on the spectrum of this gravitational radiation, and hence on theories of the early Universe, by combining experiments that cover 29 orders of magnitude in frequency. These include Planck observations of cosmic microwave background temperature and polarization power spectra and lensing, together with baryon acoustic oscillations and big bang nucleosynthesis measurements, as well as new pulsar timing array and ground-based interferometer limits. While individual experiments constrain the gravitational-wave energy density in specific frequency bands, the combination of experiments allows us to constrain cosmological parameters, including the inflationary spectral index, $n_t$, and the tensor-to-scalar ratio, $r$. Results from individual experiments include the most stringent nanohertz limit of the primordial background to date from the Parkes Pulsar Timing Array, $Ω_{\rm gw}(f)<2.3\times10^{-10}$. Observations of the cosmic microwave background alone limit the gravitational-wave spectral index at 95\% confidence to $n_t\lesssim5$ for a tensor-to-scalar ratio of $r = 0.11$. However, the combination of all the above experiments limits $n_t<0.36$. Future Advanced LIGO observations are expected to further constrain $n_t<0.34$ by 2020. When cosmic microwave background experiments detect a non-zero $r$, our results will imply even more stringent constraints on $n_t$ and hence theories of the early Universe.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

European Pulsar Timing Array Limits On An Isotropic Stochastic Gravitational-Wave Background

We present new limits on an isotropic stochastic gravitational-wave background (GWB) using a six pulsar dataset spanning 18 yr of observations from the 2015 European Pulsar Timing Array data release. Performing a Bayesian analysis, we fit simultaneously for the intrinsic noise parameters for each pulsar, along with common correlated signals including clock, and Solar System ephemeris errors, obtaining a robust 95$\%$ upper limit on the dimensionless strain amplitude $A$ of the background of $A<3.0\times 10^{-15}$ at a reference frequency of $1\mathrm{yr^{-1}}$ and a spectral index of $13/3$, corresponding to a background from inspiralling super-massive black hole binaries, constraining the GW energy density to $Ω_\mathrm{gw}(f)h^2 < 1.1\times10^{-9}$ at 2.8 nHz. We also present limits on the correlated power spectrum at a series of discrete frequencies, and show that our sensitivity to a fiducial isotropic GWB is highest at a frequency of $\sim 5\times10^{-9}$~Hz. Finally we discuss the implications of our analysis for the astrophysics of supermassive black hole binaries, and present 95$\%$ upper limits on the string tension, $Gμ/c^2$, characterising a background produced by a cosmic string network for a set of possible scenarios, and for a stochastic relic GWB. For a Nambu-Goto field theory cosmic string network, we set a limit $Gμ/c^2<1.3\times10^{-7}$, identical to that set by the {\it Planck} Collaboration, when combining {\it Planck} and high-$\ell$ Cosmic Microwave Background data from other experiments. For a stochastic relic background we set a limit of $Ω^\mathrm{relic}_\mathrm{gw}(f)h^2<1.2 \times10^{-9}$, a factor of 9 improvement over the most stringent limits previously set by a pulsar timing array.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

European Pulsar Timing Array Limits on Continuous Gravitational Waves from Individual Supermassive Black Hole Binaries

We have searched for continuous gravitational wave (CGW) signals produced by individually resolvable, circular supermassive black hole binaries (SMBHBs) in the latest EPTA dataset, which consists of ultra-precise timing data on 41 millisecond pulsars. We develop frequentist and Bayesian detection algorithms to search both for monochromatic and frequency-evolving systems. None of the adopted algorithms show evidence for the presence of such a CGW signal, indicating that the data are best described by pulsar and radiometer noise only. Depending on the adopted detection algorithm, the 95\% upper limit on the sky-averaged strain amplitude lies in the range $6\times 10^{-15}<A<1.5\times10^{-14}$ at $5{\rm nHz}<f<7{\rm nHz}$. This limit varies by a factor of five, depending on the assumed source position, and the most constraining limit is achieved towards the positions of the most sensitive pulsars in the timing array. The most robust upper limit -- obtained via a full Bayesian analysis searching simultaneously over the signal and pulsar noise on the subset of ours six best pulsars -- is $A\approx10^{-14}$. These limits, the most stringent to date at $f<10{\rm nHz}$, exclude the presence of sub-centiparsec binaries with chirp mass $\cal{M}_c>10^9$M$_\odot$ out to a distance of about 25Mpc, and with $\cal{M}_c>10^{10}$M$_\odot$ out to a distance of about 1Gpc ($z\approx0.2$). We show that state-of-the-art SMBHB population models predict $<1\%$ probability of detecting a CGW with the current EPTA dataset, consistent with the reported non-detection. We stress, however, that PTA limits on individual CGW have improved by almost an order of magnitude in the last five years. The continuing advances in pulsar timing data acquisition and analysis techniques will allow for strong astrophysical constraints on the population of nearby SMBHBs in the coming years.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Fast Radio Bursts and Radio Transients from Black Hole Batteries

Most black holes (BHs) will absorb a neutron star (NS) companion fully intact, without tidal disruption, suggesting the pair will remain dark to telescopes. Even without tidal disruption, electromagnetic luminosity is generated from the battery phase of the binary when the BH interacts with the NS magnetic field. Originally the luminosity was expected in high-energy X-rays or gamma-rays, however we conjecture that some of the battery power is emitted in the radio bandwidth. While the luminosity and timescale are suggestive of fast radio bursts (FRBs; millisecond-scale radio transients) NS--BH coalescence rates are too low to make these a primary FRB source. Instead, we propose the transients form a FRB sub-population, distinguishable by a double peak with a precursor. The rapid ramp-up in luminosity manifests as a precursor to the burst which is $20\%-80\%$ as luminous, given 0.5~ms timing resolution. The main burst is from the peak luminosity before merger. The post-merger burst follows from the NS magnetic field migration to the BH, causing a shock. NS--BH pairs are especially desirable for ground-based gravitational wave (GW) observatories since the pair might not otherwise be detected, with electromagnetic counterparts greatly augmenting the scientific leverage beyond the GW signal. Valuably, the electromagnetic signal can break degeneracies in the parameters encoded in the GW as well as probe the NS magnetic field strength, yielding insights into open problems in NS magnetic field decay.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Phase-coherent mapping of gravitational-wave backgrounds using ground-based laser interferometers

We extend the formalisms developed in Gair et al. and Cornish and van Haasteren to create maps of gravitational-wave backgrounds using a network of ground-based laser interferometers. We show that in contrast to pulsar timing arrays, which are insensitive to half of the gravitational-wave sky (the curl modes), a network of ground-based interferometers is sensitive to both the gradient and curl components of the background. The spatial separation of a network of interferometers, or of a single interferometer at different times during its rotational and orbital motion around the Sun, allows for recovery of both components. We derive expressions for the response functions of a laser interferometer in the small-antenna limit, and use these expressions to calculate the overlap reduction function for a pair of interferometers. We also construct maximum-likelihood estimates of the + and x-polarization modes of the gravitational-wave sky in terms of the response matrix for a network of ground-based interferometers, evaluated at discrete times during Earth's rotational and orbital motion around the Sun. We demonstrate the feasibility of this approach for some simple simulated backgrounds (a single point source and spatially-extended distributions having only grad or curl components), calculating maximum-likelihood sky maps and uncertainty maps based on the (pseudo)inverse of the response matrix. The distinction between this approach and standard methods for mapping gravitational-wave power is also discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Mapping gravitational-wave backgrounds using methods from CMB analysis: Application to pulsar timing arrays

We describe an alternative approach to the analysis of gravitational-wave backgrounds, based on the formalism used to characterise the polarisation of the cosmic microwave background. In contrast to standard analyses, this approach makes no assumptions about the nature of the background and so has the potential to reveal much more about the physical processes that generated it. An arbitrary background can be decomposed into modes whose angular dependence on the sky is given by gradients and curls of spherical harmonics. We derive the pulsar timing overlap reduction functions for the individual modes, which are given by simple combinations of spherical harmonics evaluated at the pulsar locations. We show how these can be used to recover the components of an arbitrary background, giving explicit results for both isotropic and anisotropic uncorrelated backgrounds. We also find that the response of a pulsar timing array to curl modes is identically zero, so half of the gravitational-wave sky will never be observed using pulsar timing, no matter how many pulsars are included in the array. An isotropic, unpolarised and uncorrelated background can be accurately represented using only three modes, and so a search of this type will be only slightly more complicated than the standard cross-correlation search using the Hellings and Downs overlap reduction function. However, by measuring the components of individual modes of the background and checking for consistency with isotropy, this approach has the potential to reveal much more information. Each individual mode on its own describes a background that is correlated between different points on the sky. A measurement of the components that indicates the presence of correlations in the background on large angular scales would suggest startling new physics.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The effect of small inter-pulsar distance variations in stochastic gravitational wave background searches with Pulsar Timing Arrays

One of the primary objectives for Pulsar Timing Arrays (PTAs) is to detect a stochastic background generated by the incoherent superposition of gravitational waves (GWs), in particular from the cosmic population of supermassive black hole binaries. Current stochastic background searches assume that pulsars in a PTA are separated from each other and the Earth by many GW wavelengths. As more millisecond pulsars are discovered and added to PTAs, some may be separated by only a few radiation wavelengths or less, resulting in correlated GW phase changes between close pulsars in the array. Here we investigate how PTA overlap reduction functions (ORFs), up to quadrupole order, are affected by these additional correlated phase changes, and how they are in turn affected by relaxing the assumption that all pulsars are equidistant from the solar system barycenter. We find that in the low frequency GW background limit of $f\sim10^{-9}$~Hz, and for pulsars at varying distances from the Earth, that these additional correlations only affect the ORFs by a few percent for pulsar pairs at large angular separations, as expected. However when nearby (order 100 pc) pulsars are separated by less than a few degrees, the correlated phase changes can introduce variations of a few tens of percent in the magnitude of the isotropic ORF, and much larger fractional differences in the anisotropic ORFs-- up to 188 in the $m=0$, $l=2$ ORF for equidistant pulsars separated by 3 degrees. In fact, the magnitude of most of the anisotropic ORFs is largest at small, but non-zero, pulsar separations. Finally, we write down a small angle approximation for the correlated phase changes which can easily be implemented in search pipelines, and for completeness, examine the behavior of the ORFs for pulsars which lie at a radiation wavelength from the Earth.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Analysis of the first IPTA Mock Data Challenge by the EPTA timing data analysis working group

This is a summary of the methods we used to analyse the first IPTA Mock Data Challenge (MDC), and the obtained results. We have used a Bayesian analysis in the time domain, accelerated using the recently developed ABC-method which consists of a form of lossy linear data compression. The TOAs were first processed with Tempo2, where the design matrix was extracted for use in a subsequent Bayesian analysis. We used different noise models to analyse the datasets: no red noise, red noise the same for all pulsars, and individual red noise per pulsar. We sampled from the likelihood with four different samplers: "emcee", "t-walk", "Metropolis-Hastings", and "pyMultiNest". All but emcee agreed on the final result, with emcee failing due to artefacts of the high-dimensionality of the problem. An interesting issue we ran into was that the prior of all the 36 (red) noise amplitudes strongly affects the results. A flat prior in the noise amplitude biases the inferred GWB amplitude, whereas a flat prior in log-amplitude seems to work well. This issue is only apparent when using a noise model with individually modelled red noise for all pulsars. Our results for the blind challenges are in good agreement with the injected values. For the GWB amplitudes we found h_c = 1.03 +/- 0.11 [10^{-14}], h_c = 5.70 +/- 0.35 [10^{-14}], and h_c = 6.91 +/- 1.72 [10^{-15}], and for the GWB spectral index we found gamma = 4.28 +/- 0.20, gamma = 4.35 +/- 0.09, and gamma = 3.75 +/- 0.40. We note that for closed challenge 3 there was quite some covariance between the signal and the red noise: if we constrain the GWB spectral index to the usual choice of gamma = 13/3, we obtain the estimates: h_c = 10.0 +/- 0.64 [10^{-15}], h_c = 56.3 +/- 2.42 [10^{-15}], and h_c = 4.83 +/- 0.50 [10^{-15}], with one-sided 2 sigma upper-limits of: h_c <= 10.98 [10^{-15}], h_c <= 60.29 [10^{-15}], and h_c <= 5.65 [10^{-15}].

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Characterising gravitational wave stochastic background anisotropy with Pulsar Timing Arrays

Detecting a stochastic gravitational wave background, particularly radiation from individually unresolvable super-massive black hole binary systems, is one of the primary targets for Pulsar Timing Arrays. Increasingly more stringent upper limits are being set on these signals under the assumption that the background radiation is isotropic. However, some level of anisotropy may be present and the characterisation of the power at different angular scales carries important information. We show that the standard analysis for isotropic backgrounds can be generalised in a conceptually straightforward way to the case of generic anisotropic background radiation by decomposing the angular distribution of the gravitational wave power on the sky into multipole moments. We introduce the concept of generalised overlap reduction functions which characterise the effect of the anisotropy multipoles on the correlation of the timing residuals from the pulsars timed by a Pulsar Timing Array. In a search for a signal characterised by a generic anisotropy, the generalised overlap reduction functions play the role of the so-called Hellings and Downs curve used for isotropic radiation. We compute the generalised overlap reduction functions for a generic level of anisotropy and Pulsar Timing Array configuration. We also provide an order of magnitude estimate of the level of anisotropy that can be expected in the background generated by super-massive black hole binary systems.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

A generalized Heckmann-Schucking cosmological solution in the presence of a negative cosmological constant

An exact solution of the Einstein equations for a Bianchi-I universe in the presence of dust, stiff matter and a negative cosmological constant, generalising the well-known Heckmann-Schucking solution is presented. This solution describes a universe existing during a finite period of cosmic time, where the beginning and the end of its evolution are characterized by the presence of Kasner type cosmological singularities.

0 citations0 reviewsNo code linkedOpen access