Source author record

Stephen R. Taylor

Stephen R. Taylor appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

20works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

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.

preprint2022arXiv

A Parallelized Bayesian Approach To Accelerated Gravitational-Wave Background Characterization

The characterization of nanohertz-frequency gravitational waves (GWs) with pulsar-timing arrays requires a continual expansion of datasets and monitored pulsars. Whereas detection of the stochastic GW background is predicated on measuring a distinctive pattern of inter-pulsar correlations, characterizing the background's spectrum is driven by information encoded in the power spectra of the individual pulsars' time series. We propose a new technique for rapid Bayesian characterization of the stochastic GW background that is fully parallelized over pulsar datasets. This Factorized Likelihood (FL) technique empowers a modular approach to parameter estimation of the GW background, multi-stage model selection of a spectrally-common stochastic process and quadrupolar inter-pulsar correlations, and statistical cross-validation of measured signals between independent pulsar sub-arrays. We demonstrate the equivalence of this technique's efficacy with the full pulsar-timing array likelihood, yet at a fraction of the required time. Our technique is fast, easily implemented, and trivially allows for new data and pulsars to be combined with legacy datasets without re-analysis of the latter.

preprint2022arXiv

Gravitational Wave Statistics for Pulsar Timing Arrays: Examining Bias from Using a Finite Number of Pulsars

Recently, many different pulsar timing array (PTA) collaborations have reported strong evidence for a common stochastic process in their data sets. The reported amplitudes are in tension with previously computed upper limits. In this paper, we investigate how using a subset of a set of pulsars biases Bayesian upper limit recovery. We generate 500 simulated PTA data sets based on the NANOGrav 11-year data set with an injected stochastic gravitational wave background (GWB). We then compute upper limits by sampling individual pulsar likelihoods, and combine them through a factorized version of the PTA likelihood to obtain upper limits on the GWB amplitude using different numbers of pulsars. We find that it is possible to recover an upper limit (95\% credible interval) \textit{below} the injected value, and that it is significantly more likely for this to occur when using a subset of pulsars to compute the upper limit. When picking pulsars to induce the maximum possible bias, we find that the 95\% Bayesian upper limit recovered is below the injected value in 10.6\% (53 of 500) of realizations. Further, we find that if we choose a subset of pulsars in order to obtain a lower upper limit than when using the full set of pulsars, the distribution of upper limits obtained from these 500 realizations is shifted to lower amplitude values.

preprint2022arXiv

Quasars with Periodic Variability: Capabilities and Limitations of Bayesian Searches for Supermassive Black Hole Binaries in Time-Domain Surveys

Supermassive black hole binaries (SMBHBs) are an inevitable consequence of galaxy mergers. At subparsec separations, they are practically impossible to resolve, and the most promising technique is to search for quasars with periodic variability. However, searches for quasar periodicity in time-domain data are challenging due to the stochastic variability of quasars. In this paper, we used Bayesian methods to disentangle periodic SMBHB signals from intrinsic damped random walk (DRW) variability in active galactic nuclei light curves. We simulated a wide variety of realistic DRW and DRW+sine light curves. Their observed properties are modeled after the Catalina Real-time Transient Survey (CRTS) and expected properties of the upcoming Legacy Survey of Space and Time (LSST) from the Vera C. Rubin Observatory. Through a careful analysis of parameter estimation and Bayesian model selection, we investigated the range of parameter space for which binary systems can be detected. We also examined which DRW signals can mimic periodicity and be falsely classified as binary candidates. We found that periodic signals are more easily detectable if the period is short or the amplitude of the signal is large compared to the contribution of the DRW noise. We saw similar detection rates both in the CRTS and LSST-like simulations, while the false-detection rate depends on the quality of the data and is minimal in LSST. Our idealized simulations provide an excellent way to uncover the intrinsic limitations in quasar periodicity searches and set the stage for future searches for SMBHBs.

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.

preprint2021arXiv

Multi-messenger Time-domain Signatures Of Supermassive Black Hole Binaries

Supermassive black hole binaries (SMBHBs) are a natural outcome of galaxy mergers and should form frequently in galactic nuclei. Sub-parsec binaries can be identified from their bright electromagnetic emission, e.g., Active Galactic Nuclei (AGN) with Doppler shifted broad emission lines or AGN with periodic variability, as well as from the emission of strong gravitational radiation. The most massive binaries (with total mass >10^8 M_sol) emit in the nanohertz band and are targeted by Pulsar Timing Arrays (PTAs). Here we examine the synergy between electromagnetic and gravitational wave signatures of SMBHBs. We connect both signals to the orbital dynamics of the binary and examine the common link between them, laying the foundation for joint multi-messenger observations. We find that periodic variability arising from relativistic Doppler boost is the most promising electromagnetic signature to connect with GWs. We delineate the parameter space (binary total mass/chirp mass versus binary period/GW frequency) for which joint observations are feasible. Currently multi-messenger detections are possible only for the most massive and nearby galaxies, limited by the sensitivity of PTAs. However, we demonstrate that as PTAs collect more data in the upcoming years, the overlapping parameter space is expected to expand significantly.

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.

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.

preprint2020arXiv

Pulsar timing array signals induced by black hole binaries in relativistic eccentric orbits

Individual supermassive black hole binaries in non-circular orbits are possible nanohertz gravitational wave sources for the rapidly maturing Pulsar Timing Array experiments. We develop an accurate and efficient approach to compute Pulsar Timing Array signals due to gravitational waves from inspiraling supermassive black hole binaries in relativistic eccentric orbits. Our approach employs a Keplerian-type parametric solution to model third post-Newtonian accurate precessing eccentric orbits while a novel semi-analytic prescription is provided to model the effects of quadrupolar order gravitational wave emission. These inputs lead to a semi-analytic prescription to model such signals, induced by non-spinning black hole binaries inspiralling along arbitrary eccentricity orbits. Additionally, we provide a fully analytic prescription to model Pulsar Timing Array signals from black hole binaries inspiraling along moderately eccentric orbits, influenced by Boetzelet al.[Phys. Rev. D 96,044011(2017)]. These approaches are being incorporated into Enterprise and TEMPO2 for searching the presence of such binaries in Pulsar Timing Array datasets.

preprint2015arXiv

Detection of eccentric supermassive black hole binaries with pulsar timing arrays: Signal-to-noise ratio calculations

We present a detailed analysis of the expected signal-to-noise ratios of supermassive black hole binaries on eccentric orbits observed by pulsar timing arrays. We derive several analytical relations that extend the results of Peters and Mathews [Phys. Rev. D 131, 435 (1963)] to quantify the impact of eccentricity in the detection of single resolvable binaries in the pulsar timing array band. We present ready-to-use expressions to compute the increase/loss in signal-to-noise ratio of eccentric single resolvable sources whose dominant harmonic is located in the low/high frequency sensitivity regime of pulsar timing arrays. Building upon the work of Phinney (arXiv:astro-ph/0108028) and Enoki and Nagashima [Prog. Theor. Phys. 117, 241 (2007)], we present an analytical framework that enables the construction of rapid spectra for a stochastic gravitational-wave background generated by a cosmological population of eccentric sources. We confirm previous findings which indicate that, relative to a population of quasicircular binaries, the strain of a stochastic, isotropic gravitational-wave background generated by a cosmological population of eccentric binaries will be suppressed in the frequency band of pulsar timing arrays. We quantify this effect in terms of signal-to-noise ratios in a pulsar timing array.

preprint2015arXiv

Estimating the sensitivity of pulsar timing arrays

The sensitivity curve of a canonical pulsar timing array is calculated for two types of source: a monochromatic wave and a stochastic background. These calculations are performed in both a Bayesian and frequentist framework, using both analytical and numerical methods. These calculations are used to clarify the interpretation of the sensitivity curves and to illustrate the sometimes overlooked fact that the sensitivity curve depends not only on the properties of the pulse time-of-arrival data set but also on the properties of the source being observed. The Bayesian and frequentist frameworks were found to give consistent results and the analytic and numerical calculations were also found to be in good agreement.

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.

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.

preprint2015arXiv

Mapping gravitational-wave backgrounds in modified theories of gravity using pulsar timing arrays

We extend our previous work on applying CMB techniques to the mapping of gravitational-wave backgrounds to backgrounds which have non-GR polarisations. Our analysis and results are presented in the context of pulsar-timing array observations, but the overarching methods are general, and can be easily applied to LIGO or eLISA observations using appropriately modified response functions. Analytic expressions for the pulsar-timing response to gravitational waves with non-GR polarisation are given for each mode of a spin-weighted spherical-harmonic decomposition of the background, which permit the signal to be mapped across the sky to any desired resolution. We also derive the pulsar-timing overlap reduction functions for the various non-GR polarisations, finding analytic forms for anisotropic backgrounds with scalar-transverse ("breathing") and vector-longitudinal polarisations, and a semi-analytic form for scalar-longitudinal backgrounds. Our results indicate that pulsar-timing observations will be completely insensitive to scalar-transverse mode anisotropies in the polarisation amplitude beyond dipole, and anisotropies in the power beyond quadrupole. Analogously to our previous findings that pulsar-timing observations lack sensitivity to tensor-curl modes for a transverse-traceless tensor background, we also find insensitivity to vector-curl modes for a vector-longitudinal background.

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.

preprint2013arXiv

Searching For Anisotropic Gravitational-wave Backgrounds Using Pulsar Timing Arrays

We present the results of simulated injections testing the first Bayesian search-pipeline capable of investigating the angular-structure of a gravitational-wave (GW) background influencing pulsar signals. A stochastic background of GWs from the incoherent superposition of many inspiraling supermassive black hole binaries at nHz frequencies is likely to be the dominant GW signal detectable by pulsar timing arrays (PTAs). Even though one might expect a background composed of a high-redshift cosmological population of sources to be fairly isotropic, deviations from isotropy may be indicative of local GW hotspots or some form of continuous anisotropy in the angular-distribution of GW-power. A GWB induces time-of-arrival deviations in pulsar signals which are correlated between separated pulsars. In an isotropic background this cross-correlation follows a distinctive relationship, known as the Hellings and Downs curve, that depends only on the angular separation of the pulsars. If the background is anisotropic, the cross-correlation is different, but predictable, and also depends on the absolute position of the pulsars. By simulating datasets containing GWBs with various anisotropic configurations, we have explored the prospects for constraining anisotropy using near future data. We find that at moderate to high signal to noise ratio the assumption of isotropy is no longer an appropriate description of the simulated background. Furthermore, we can recover the nature of the injected anisotropy in a Bayesian parameter-estimation search, and propose a prior on the anisotropy search-space motivated by the physicality of the implied distribution of sources.

preprint2012arXiv

Cosmology with the lights off: Standard sirens in the Einstein Telescope era

We explore the prospects for constraining cosmology using gravitational-wave (GW) observations of neutron-star binaries by the proposed Einstein Telescope (ET), exploiting the narrowness of the neutron-star mass function. Double neutron-star (DNS) binaries are expected to be one of the first sources detected after "first-light" of Advanced LIGO and are expected to be detected at a rate of a few tens per year in the advanced era. However the proposed ET could catalog tens of thousands per year. Combining the measured source redshift distributions with GW-network distance determinations will permit not only the precision measurement of background cosmological parameters, but will provide an insight into the astrophysical properties of these DNS systems. Of particular interest will be to probe the distribution of delay times between DNS-binary creation and subsequent merger, as well as the evolution of the star-formation rate density within ET's detection horizon. Keeping H_0, Ω_{m,0} and Ω_{Λ,0} fixed and investigating the precision with which the dark-energy equation-of-state parameters could be recovered, we found that with 10^5 detected DNS binaries we could constrain these parameters to an accuracy similar to forecasted constraints from future CMB+BAO+SNIa measurements. Furthermore, modeling the merger delay-time distribution as a power-law, and the star-formation rate (SFR) density as a parametrized version of the Porciani and Madau SF2 model, we find that the associated astrophysical parameters are constrained to within ~ 10%. All parameter precisions scaled as 1/sqrt(N), where N is the number of cataloged detections. We also investigated how precisions varied with the intrinsic underlying properties of the Universe and with the distance reach of the network (which may be affected by the low-frequency cutoff of the detector).

preprint2012arXiv

Hubble without the Hubble: cosmology using advanced gravitational-wave detectors alone

We investigate a novel approach to measuring the Hubble constant using gravitational-wave (GW) signals from compact binaries by exploiting the narrowness of the distribution of masses of the underlying neutron-star population. Gravitational-wave observations with a network of detectors will permit a direct, independent measurement of the distance to the source systems. If the redshift of the source is known, these inspiraling double-neutron-star binary systems can be used as standard sirens to extract cosmological information. Unfortunately, the redshift and the system chirp mass are degenerate in GW observations. Thus, most previous work has assumed that the source redshift is obtained from electromagnetic counterparts. In this paper, we explore what we can learn about the background cosmology and the mass distribution of neutron stars from the set of neutron-star (NS) mergers detected by such a network. We use a Bayesian formalism to analyze catalogs of NS-NS inspiral detections. We find that it is possible to constrain the Hubble constant, H_0, and the parameters of the NS mass function using gravitational-wave data alone, without relying on electromagnetic counterparts. Under reasonable assumptions, we will be able to determine H_0 to +/- 10% using ~100 observations, provided the Gaussian half-width of the underlying double NS mass distribution is less than 0.04 solar masses. The expected precision depends linearly on the intrinsic width of the NS mass function, but has only a weak dependence on H_0 near the default parameter values. Finally, we consider what happens if, for some fraction of our data catalog, we have an electromagnetically measured redshift. The detection, and cataloging, of these compact-object mergers will allow precision astronomy, and provide a determination of H_0 which is independent of the local distance scale.

preprint2012arXiv

Using Swarm Intelligence To Accelerate Pulsar Timing Analysis

We provide brief notes on a particle swarm-optimisation approach to constraining the properties of a stochastic gravitational-wave background in the first International Pulsar Timing Array data-challenge. The technique employs many computational-agents which explore parameter space, remembering their most optimal positions and also sharing this information with all other agents. It is this sharing of information which accelerates the convergence of all agents to the global best-fit location in a very short number of iterations. Error estimates can also be provided by fitting a multivariate Gaussian to the recorded fitness of all visited points.

preprint2012arXiv

Weighing The Evidence For A Gravitational-Wave Background In The First International Pulsar Timing Array Data Challenge

We describe an analysis of the First International Pulsar Timing Array Data Challenge. We employ a robust, unbiased Bayesian framework developed by van Haasteren to study the three Open and Closed datasets, testing various models for each dataset and using MultiNest to recover the evidence for the purposes of Bayesian model-selection. The parameter constraints of the favoured model are confirmed using an adaptive MCMC technique. Our results for Closed1 favoured a gravitational-wave background with strain amplitude at f=1 yr-1, A, of (1.1 +/- 0.1) x 10^{-14}, power spectral-index gamma=4.30 +/- 0.15 and no evidence for red-timing noise or single-sources. The evidence for Closed2 favours a gravitational-wave background with A=(6.1 +/- 0.3) x 10^{-14}, gamma=4.34 +/- 0.09 with no red-timing noise or single-sources. Finally, the evidence for Closed3 favours the presence of red-timing noise and a gravitational-wave background, with no single-sources. The properties of the background were A=(5 +/- 1) x 10^{-15} and gamma=4.23 +/- 0.35, while the properties of the red-noise were N_{red}=(12 +/- 4) ns and gamma_{red}=1.5 +/- 0.3. In all cases the redness of the recovered background is consistent with a source-population of inspiraling supermassive black-hole binaries. We also investigate the effect that down-sampling of the datasets has on parameter constraints and run-time. Finally we provide a proof-of-principle study of the ability of the Bayesian framework used in this paper to reconstruct the angular correlation of gravitational-wave background induced timing-residuals, comparing this to the Hellings and Downs curve.