Search

We searched for an isotropic stochastic gravitational wave background in the second data release of the International Pulsar Timing Array, a global collaboration synthesizing decadal-length pulsar-timing campaigns in North America, Europe, and Australia. In our reference search for a power law strain spectrum of the form $h_c = A(f/1\,\mathrm{yr}^{-1})^{\alpha}$, we found strong evidence for a spectrally-similar low-frequency stochastic process of amplitude $A = 3.8^{+6.3}_{-2.5}\times10^{-15}$ and spectral index $\alpha = -0.5 \pm 0.5$, where the uncertainties represent 95\% credible regions, using information from the auto- and cross-correlation terms between the pulsars in the array. For a spectral index of $\alpha = -2/3$, as expected from a population of inspiralling supermassive black hole binaries, the recovered amplitude is $A = 2.8^{+1.2}_{-0.8}\times10^{-15}$. Nonetheless, no significant evidence of the Hellings-Downs correlations that would indicate a gravitational-wave origin was found. We also analyzed the constituent data from the individual pulsar timing arrays in a consistent way, and clearly demonstrate that the combined international data set is more sensitive. Furthermore, we demonstrate that this combined data set produces comparable constraints to recent single-array data sets which have more data than the constituent parts of the combination. Future international data releases will deliver increased sensitivity to gravitational wave radiation, and significantly increase the detection probability., Comment: 17 pages, 12 figures, accepted in MNRAS

In Paper I of this series, we detected a significant value of the braking index ($n$) for 19 young, high-$\dot{E}$ radio pulsars using $\sim$ 10 years of timing observations from the 64-m Parkes radio telescope. Here we investigate this result in more detail using a Bayesian pulsar timing framework to model timing noise and to perform selection to distinguish between models containing exponential glitch recovery and braking index signatures. We show that consistent values of $n$ are maintained with the addition of substantial archival data, even in the presence of glitches. We provide strong arguments that our measurements are unlikely due to exponential recovery signals from unseen glitches even though glitches play a key role in the evolution of a pulsar's spin frequency. We conclude that, at least over decadal time scales, the value of $n$ can be significantly larger than the canonical 3 and discuss the implications for the evolution of pulsars., Comment: 15 pages, 10 figures and 5 tables

In this paper, we describe the International Pulsar Timing Array second data release, which includes recent pulsar timing data obtained by three regional consortia: the European Pulsar Timing Array, the North American Nanohertz Observatory for Gravitational Waves, and the Parkes Pulsar Timing Array. We analyse and where possible combine high-precision timing data for 65 millisecond pulsars which are regularly observed by these groups. A basic noise analysis, including the processes which are both correlated and uncorrelated in time, provides noise models and timing ephemerides for the pulsars. We find that the timing precisions of pulsars are generally improved compared to the previous data release, mainly due to the addition of new data in the combination. The main purpose of this work is to create the most up-to-date IPTA data release. These data are publicly available for searches for low-frequency gravitational waves and other pulsar science., Comment: Submitted to MNRAS and in review, 23 pages, 5 figures

The smooth spin-down of young pulsars is perturbed by two non-deterministic phenomenon, glitches and timing noise. Although the timing noise provides insights into nuclear and plasma physics at extreme densities, it acts as a barrier to high-precision pulsar timing experiments. An improved methodology based on Bayesian inference is developed to simultaneously model the stochastic and deterministic parameters for a sample of 85 high-$\dot{E}$ radio pulsars observed for $\sim$ 10 years with the 64-m Parkes radio telescope. Timing noise is known to be a red process and we develop a parametrization based on the red-noise amplitude ($A_{\rm red}$) and spectral index ($\beta$). We measure the median $A_{\rm red}$ to be $-10.4^{+1.8}_{-1.7}$ yr$^{3/2}$ and $\beta$ to be $-5.2^{+3.0}_{-3.8}$ and show that the strength of timing noise scales proportionally to $\nu^{1}|\dot{\nu}|^{-0.6\pm0.1}$, where $\nu$ is the spin frequency of the pulsar and $\dot{\nu}$ its spin-down rate. Finally, we measure significant braking indices for 19 pulsars, proper motions for two pulsars and discuss the presence of periodic modulation in the arrival times of five pulsars.

We search for continuous gravitational waves (CGWs) produced by individual super-massive black-hole binaries (SMBHBs) in circular orbits using high-cadence timing observations of PSR J1713$+$0747. We observe this millisecond pulsar using the telescopes in the European Pulsar Timing Array (EPTA) with an average cadence of approximately 1.6 days over the period between April 2011 and July 2015, including an approximately daily average between February 2013 and April 2014. The high-cadence observations are used to improve the pulsar timing sensitivity across the GW frequency range of $0.008-5$ $\mu$Hz. We use two algorithms in the analysis, including a spectral fitting method and a Bayesian approach. For an independent comparison, we also use a previously published Bayesian algorithm. We find that the Bayesian approaches provide optimal results and the timing observations of the pulsar place a 95 per cent upper limit on the sky-averaged strain amplitude of CGWs to be $\lesssim3.5\times10^{-13}$ at a reference frequency of 1 $\mu$Hz. We also find a 95 per cent upper limit on the sky-averaged strain amplitude of low-frequency CGWs to be $\lesssim1.4\times10^{-14}$ at a reference frequency of 20~nHz., Comment: 12 pages, 4 figures. Accepted for publication in MNRAS

Symmetries play an important role in modern theories of gravity. The strong equivalence principle (SEP) constitutes a collection of gravitational symmetries which are all implemented by general relativity. Alternative theories, however, are generally expected to violate some aspects of SEP. We test three aspects of SEP using observed change rates in the orbital period and eccentricity of binary pulsar J1713+0747: 1. the gravitational constant's constancy as part of locational invariance of gravitation; 2. the post-Newtonian parameter $\hat{\alpha}_3$ in gravitational Lorentz invariance; 3. the universality of free fall (UFF) for strongly self-gravitating bodies. Based on the pulsar timing result of the combined dataset from the North American Nanohertz Gravitational Observatory (NANOGrav) and the European Pulsar Timing Array (EPTA), we find $\dot{G}/G = (-0.1 \pm 0.9) \times 10^{-12}\,{\rm yr}^{-1}$, which is weaker than Solar system limits, but applies for strongly self-gravitating objects. Furthermore, we obtain the constraints $|\Delta|< 0.002$ for the UFF test and $-3\times10^{-20} < \hat{\alpha}_3 < 4\times10^{-20}$ at 95% confidence. These are the first direct UFF and $\hat{\alpha}_3$ tests based on pulsar binaries, and they overcome various limitations of previous tests., Comment: 11 pages, 6 figures, and 1 table; submitted to MNRAS

We present an extremely deep CO(1-0) observation of a confirmed $z=1.62$ galaxy cluster. We detect two spectroscopically confirmed cluster members in CO(1-0) with $S/N>5$. Both galaxies have log(${\cal M_{\star}}$/\msol)$>11$ and are gas rich, with ${\cal M}_{\rm mol}$/(${\cal M_{\star}}+{\cal M}_{\rm mol}$)$\sim 0.17-0.45$. One of these galaxies lies on the star formation rate (SFR)-${\cal M_{\star}}$ sequence while the other lies an order of magnitude below. We compare the cluster galaxies to other SFR-selected galaxies with CO measurements and find that they have CO luminosities consistent with expectations given their infrared luminosities. We also find that they have comparable gas fractions and star formation efficiencies (SFE) to what is expected from published field galaxy scaling relations. The galaxies are compact in their stellar light distribution, at the extreme end for all high redshift star-forming galaxies. However, their SFE is consistent with other field galaxies at comparable compactness. This is similar to two other sources selected in a blind CO survey of the HDF-N. Despite living in a highly quenched proto-cluster core, the molecular gas properties of these two galaxies, one of which may be in the processes of quenching, appear entirely consistent with field scaling relations between the molecular gas content, stellar mass, star formation rate, and redshift. We speculate that these cluster galaxies cannot have any further substantive gas accretion if they are to become members of the dominant passive population in $z<1$ clusters., Comment: 19 pages, 11 figures, 3 tables. Accepted for publication in the Astrophysical Journal

We present a robust approach to incorporating models for the time-variable broadening of the pulse profile due to scattering in the ionized interstellar medium into profile-domain pulsar timing analysis. We use this approach to simultaneously estimate temporal variations in both the dispersion measure (DM) and scattering, together with a model for the pulse profile that includes smooth evolution as a function of frequency, and the pulsar's timing model. We show that fixing the scattering timescales when forming time-of-arrival estimates, as has been suggested in the context of traditional pulsar timing analysis, can significantly underestimate the uncertainties in both DM, and the arrival time of the pulse, leading to bias in the timing parameters. We apply our method using a new, publicly available, GPU accelerated code, both to simulations, and observations of the millisecond pulsar PSR J1643$-$1224. This pulsar is known to exhibit significant scattering variability compared to typical millisecond pulsars, and we find including low-frequency ($< 1$ GHz) data without a model for these scattering variations leads to significant periodic structure in the DM, and also biases the astrometric parameters at the $4\sigma$ level, for example, changing proper motion in right ascension by $0.50 \pm 0.12$. If low frequency observations are to be included when significant scattering variations are present, we conclude it is necessary to not just model those variations, but also to sample the parameters that describe the variations simultaneously with all other parameters in the model, a task for which profile domain pulsar timing is ideally suited., Comment: 12 pages, 6 figures, Accepted to MNRAS

We introduce a new method for performing robust Bayesian estimation of the three-dimensional spatial power spectrum at the Epoch of Reionization (EoR), from interferometric observations. The versatility of this technique allows us to present two approaches. First, when the observations span only a small number of independent spatial frequencies ($k$-modes) we sample directly from the spherical power spectrum coefficients that describe the EoR signal realisation. Second, when the number of $k$-modes to be included in the model becomes large, we sample from the joint probability density of the spherical power spectrum and the signal coefficients, using Hamiltonian Monte Carlo methods to explore this high dimensional ($\sim$ 20000) space efficiently. This approach has been successfully applied to simulated observations that include astrophysically realistic foregrounds in a companion publication (Sims et al. 2016). Here we focus on explaining the methodology in detail, and use simple foreground models to both demonstrate its efficacy, and highlight salient features. In particular, we show that including an arbitrary flat spectrum continuum foreground that is $10^8$ times greater in power than the EoR signal has no detectable impact on our parameter estimates of the EoR power spectrum recovered from the data., Comment: 16 pages, 6 figures. Submitted to MNRAS

We extend profile domain pulsar timing to incorporate wide-band effects such as frequency-dependent profile evolution and broadband shape variation in the pulse profile. We also incorporate models for temporal variations in both pulse width and in the separation in phase of the main pulse and interpulse. We perform the analysis with both nested sampling and Hamiltonian Monte Carlo methods. In the latter case we introduce a new parameterisation of the posterior that is extremely efficient in the low signal-to-noise regime and can be readily applied to a wide range of scientific problems. We apply this methodology to a series of simulations, and to between seven and nine yr of observations for PSRs J1713$+$0747, J1744$-$1134, and J1909$-$3744 with frequency coverage that spans 700-3600MHz. We use a smooth model for profile evolution across the full frequency range, and compare smooth and piecewise models for the temporal variations in DM. We find the profile domain framework consistently results in improved timing precision compared to the standard analysis paradigm by as much as 40% for timing parameters. Incorporating smoothness in the DM variations into the model further improves timing precision by as much as 30%. For PSR J1713+0747 we also detect pulse shape variation uncorrelated between epochs, which we attribute to variation intrinsic to the pulsar at a level consistent with previously published analyses. Not accounting for this shape variation biases the measured arrival times at the level of $\sim$30ns, the same order of magnitude as the expected shift due to gravitational-waves in the pulsar timing band., Comment: 22 pages, submitted to MNRAS

Pulsar timing has enabled some of the strongest tests of fundamental physics. Central to the technique is the assumption that the detected radio pulses can be used to accurately measure the rotation of the pulsar. Here we report on a broad-band variation in the pulse profile of the millisecond pulsar J1643-1224. A new component of emission suddenly appears in the pulse profile, decays over 4 months, and results in a permanently modified pulse shape. Profile variations such as these may be the origin of timing noise observed in other millisecond pulsars. The sensitivity of pulsar-timing observations to gravitational radiation can be increased by accounting for this variability., Comment: Accepted to ApJ Letters

We construct foreground simulations comprising spatially correlated extragalactic and diffuse Galactic emission components and calculate the `intrinsic' (instrument-free) two-dimensional spatial power spectrum and the cylindrically and spherically averaged three-dimensional k-space power spectra of the Epoch of Reionization (EoR) and our foreground simulations using a Bayesian power spectral estimation framework. This leads us to identify a model dependent region of optimal signal estimation for our foreground and EoR models, within which the spatial power in the EoR signal relative to foregrounds is maximised. We identify a target field dependent region, in k-space, of intrinsic foreground power spectral contamination at low k_perp and k_parallel and a transition to a relatively foreground-free intrinsic EoR window in the complement to this region. The contaminated region of k-space demonstrates that simultaneous estimation of the EoR and foregrounds is important for obtaining statistically robust estimates of the EoR power spectrum; biased results will be obtained from methodologies that ignore their covariance. Using simulated observations with frequency dependent uv-coverage and primary beam, with the former derived for HERA in 37-antenna and 331-antenna configuration, we recover instrumental power spectra consistent with their intrinsic counterparts. We discuss the implications of these results for optimal strategies for unbiased estimation of the EoR power spectrum., Comment: 25 pages, 16 figures, accepted for publication in MNRAS

We present two methods for determining the significance of a stochastic gravitational-wave (GW) background affecting a pulsar-timing array, where detection is based on evidence for quadrupolar spatial correlations between pulsars. Rather than constructing noise simulations, we eliminate the GWB spatial correlations in the true datasets to assess detection significance with all real data features intact. In our first method, we perform random phase shifts in the signal-model basis functions. This phase shifting eliminates signal phase coherence between pulsars, while keeping the statistical properties of the pulsar timing residuals intact. We then explore a method to null correlations between pulsars by using a "scrambled" overlap-reduction function in the signal model for the array. This scrambled function is orthogonal to what we expect of a real GW background signal. We demonstrate the efficacy of these methods using Bayesian model selection on a set of simulated datasets that contain a stochastic GW signal, timing noise, undiagnosed glitches, and uncertainties in the Solar system ephemeris. Finally, we introduce an overarching formalism under which these two techniques are naturally linked. These methods are immediately applicable to all current pulsar-timing array datasets, and should become standard tools for future analyses., Comment: 14 pages, 8 figures, 1 table. Updated to match published version

We present evidence for a small glitch in the spin evolution of the millisecond pulsar J0613$-$0200, using the EPTA Data Release 1.0, combined with Jodrell Bank analogue filterbank TOAs recorded with the Lovell telescope and Effelsberg Pulsar Observing System TOAs. A spin frequency step of 0.82(3) nHz and frequency derivative step of ${-1.6(39) \times 10^{-19}\,\text{Hz} \ \text{s}^{-1}}$ are measured at the epoch of MJD 50888(30). After PSR B1821$-$24A, this is only the second glitch ever observed in a millisecond pulsar, with a fractional size in frequency of ${\Delta \nu/\nu=2.5(1) \times 10^{-12}}$, which is several times smaller than the previous smallest glitch. PSR J0613$-$0200 is used in gravitational wave searches with pulsar timing arrays, and is to date only the second such pulsar to have experienced a glitch in a combined 886 pulsar-years of observations. We find that accurately modelling the glitch does not impact the timing precision for pulsar timing array applications. We estimate that for the current set of millisecond pulsars included in the International Pulsar Timing Array, there is a probability of $\sim 50$% that another glitch will be observed in a timing array pulsar within 10 years., Comment: accepted for publication in MNRAS

We present proper motions for 15 pulsars which are observed regularly by the Nanshan 25-m radio telescope. Two methods, the frequentist method (Coles et al.2011) and the Bayesian (Lentati et al. 2014) method, are used and the results are compared. We demonstrate that the two methods can be applied to young pulsar data sets that exhibit large amounts of timing noise with steep spectral exponents and give consistent results. The measured positions also agree with very-long-baseline interferometric positions. Proper motions for four pulsars are obtained for the first time, and improved values are obtained for five pulsars., Comment: 8 pages, 2 figures and 7 tables. Accepted for publication in MNRAS

Pulsars show two classes of rotational irregularities that can be used to understand neutron-star interiors and magnetospheres: glitches and timing noise. Here we present an analysis of the Vela pulsar spanning nearly 21 yr of observation and including 8 glitches. We identify the relative pulse number of all of the observations between glitches, with the only pulse-number ambiguities existing over glitch events. We use the phase coherence of the timing solution to simultaneously model the timing noise and glitches in a Bayesian framework, allowing us to select preferred models for both. We find the glitches can be described using only permanent and transient changes in spin frequency, i.e., no step changes in frequency derivative. For all of the glitches, we only need two exponentially decaying changes in spin frequency to model the transient components. In contrast to previous studies, we find that the dominant transient components decay on a common $\approx$ 1300 d time scale, and that a larger fraction ( $\gtrsim 25\%$) of glitch amplitudes are associated with these transient components. We also detect shorter-duration transient components of $\approx$ 25 d, as previously observed, but are limited in sensitivity to events with shorter durations by the cadence of our observations. The timing noise is well described by a steep power-law process that is independent of the glitches and subdominant to the glitch recovery. The braking index is constrained to be $<$ 8 with 95% confidence. This methodology can be used to robustly measure the properties of glitches and timing noise in other pulsars., Comment: Accepted by MNRAS. 10 pages, 1 Figure, 2 tables

We report on 22 yrs of radio timing observations of the millisecond pulsar J1024$-$0719 by the telescopes participating in the European Pulsar Timing Array (EPTA). These observations reveal a significant second derivative of the pulsar spin frequency and confirm the discrepancy between the parallax and Shklovskii distances that has been reported earlier. We also present optical astrometry, photometry and spectroscopy of 2MASS J10243869$-$0719190. We find that it is a low-metallicity main-sequence star (K7V spectral type, $\mathrm{[M/H]}=-1.0$, $T_\mathrm{eff}=4050\pm50$ K) and that its position, proper motion and distance are consistent with those of PSR J1024$-$0719. We conclude that PSR J1024$-$0719 and 2MASS J10243869$-$0719190 form a common proper motion pair and are gravitationally bound. The gravitational interaction between the main-sequence star and the pulsar accounts for the spin frequency derivatives, which in turn resolves the distance discrepancy. Our observations suggest that the pulsar and main-sequence star are in an extremely wide ($P_\mathrm{b}>200$ yr) orbit. Combining the radial velocity of the companion and proper motion of the pulsar, we find that the binary system has a high spatial velocity of $384\pm45$ km s$^{-1}$ with respect to the local standard of rest and has a Galactic orbit consistent with halo objects. Since the observed main-sequence companion star cannot have recycled the pulsar to millisecond spin periods, an exotic formation scenario is required. We demonstrate that this extremely wide-orbit binary could have evolved from a triple system that underwent an asymmetric supernova explosion, though find that significant fine-tuning during the explosion is required. Finally, we discuss the implications of the long period orbit on the timing stability of PSR J1024$-$0719 in light of its inclusion in pulsar timing arrays., Comment: 16 pages, 12 figures, 2 tables, submitted for publication in MNRAS

We report on the high-precision timing of 42 radio millisecond pulsars (MSPs) observed by the European Pulsar Timing Array (EPTA). This EPTA Data Release 1.0 extends up to mid-2014 and baselines range from 7-18 years. It forms the basis for the stochastic gravitational-wave background, anisotropic background, and continuous-wave limits recently presented by the EPTA elsewhere. The Bayesian timing analysis performed with TempoNest yields the detection of several new parameters: seven parallaxes, nine proper motions and, in the case of six binary pulsars, an apparent change of the semi-major axis. We find the NE2001 Galactic electron density model to be a better match to our parallax distances (after correction from the Lutz-Kelker bias) than the M2 and M3 models by Schnitzeler (2012). However, we measure an average uncertainty of 80\% (fractional) for NE2001, three times larger than what is typically assumed in the literature. We revisit the transverse velocity distribution for a set of 19 isolated and 57 binary MSPs and find no statistical difference between these two populations. We detect Shapiro delay in the timing residuals of PSRs J1600$-$3053 and J1918$-$0642, implying pulsar and companion masses $m_p=1.22_{-0.35}^{+0.5} \text{M}_{\odot}$, $m_c = 0.21_{-0.04}^{+0.06} \text{M}_{\odot }$ and $m_p=1.25_{-0.4}^{+0.6} \text{M}_{\odot}$, $m_c = 0.23_{-0.05}^{+0.07} \text{M}_{\odot }$, respectively. Finally, we use the measurement of the orbital period derivative to set a stringent constraint on the distance to PSRs J1012$+$5307 and J1909$-$3744, and set limits on the longitude of ascending node through the search of the annual-orbital parallax for PSRs J1600$-$3053 and J1909$-$3744., Comment: 42 pages, 11 figures. Accepted for publication in MNRAS

We analyse the stochastic properties of the 49 pulsars that comprise the first International Pulsar Timing Array (IPTA) data release. We use Bayesian methodology, performing model selection to determine the optimal description of the stochastic signals present in each pulsar. In addition to spin-noise and dispersion-measure (DM) variations, these models can include timing noise unique to a single observing system, or frequency band. We show the improved radio-frequency coverage and presence of overlapping data from different observing systems in the IPTA data set enables us to separate both system and band-dependent effects with much greater efficacy than in the individual PTA data sets. For example, we show that PSR J1643$-$1224 has, in addition to DM variations, significant band-dependent noise that is coherent between PTAs which we interpret as coming from time-variable scattering or refraction in the ionised interstellar medium. Failing to model these different contributions appropriately can dramatically alter the astrophysical interpretation of the stochastic signals observed in the residuals. In some cases, the spectral exponent of the spin noise signal can vary from 1.6 to 4 depending upon the model, which has direct implications for the long-term sensitivity of the pulsar to a stochastic gravitational-wave (GW) background. By using a more appropriate model, however, we can greatly improve a pulsar's sensitivity to GWs. For example, including system and band-dependent signals in the PSR J0437$-$4715 data set improves the upper limit on a fiducial GW background by $\sim 60\%$ compared to a model that includes DM variations and spin-noise only., Comment: 29 pages. 16 figures. Accepted for publication in MNRAS

The highly stable spin of neutron stars can be exploited for a variety of (astro-)physical investigations. In particular arrays of pulsars with rotational periods of the order of milliseconds can be used to detect correlated signals such as those caused by gravitational waves. Three such "Pulsar Timing Arrays" (PTAs) have been set up around the world over the past decades and collectively form the "International" PTA (IPTA). In this paper, we describe the first joint analysis of the data from the three regional PTAs, i.e. of the first IPTA data set. We describe the available PTA data, the approach presently followed for its combination and suggest improvements for future PTA research. Particular attention is paid to subtle details (such as underestimation of measurement uncertainty and long-period noise) that have often been ignored but which become important in this unprecedentedly large and inhomogeneous data set. We identify and describe in detail several factors that complicate IPTA research and provide recommendations for future pulsar timing efforts. The first IPTA data release presented here (and available online) is used to demonstrate the IPTA's potential of improving upon gravitational-wave limits placed by individual PTAs by a factor of ~2 and provides a 2-sigma limit on the dimensionless amplitude of a stochastic GWB of 1.7x10^{-15} at a frequency of 1 yr^{-1}. This is 1.7 times less constraining than the limit placed by (Shannon et al. 2015), due mostly to the more recent, high-quality data they used., Comment: 25 pages, 6 tables, 5 figures. Accepted for publication in MNRAS

The sensitivity of Pulsar Timing Arrays to gravitational waves depends on the noise present in the individual pulsar timing data. Noise may be either intrinsic or extrinsic to the pulsar. Intrinsic sources of noise will include rotational instabilities, for example. Extrinsic sources of noise include contributions from physical processes which are not sufficiently well modelled, for example, dispersion and scattering effects, analysis errors and instrumental instabilities. We present the results from a noise analysis for 42 millisecond pulsars (MSPs) observed with the European Pulsar Timing Array. For characterising the low-frequency, stochastic and achromatic noise component, or "timing noise", we employ two methods, based on Bayesian and frequentist statistics. For 25 MSPs, we achieve statistically significant measurements of their timing noise parameters and find that the two methods give consistent results. For the remaining 17 MSPs, we place upper limits on the timing noise amplitude at the 95% confidence level. We additionally place an upper limit on the contribution to the pulsar noise budget from errors in the reference terrestrial time standards (below 1%), and we find evidence for a noise component which is present only in the data of one of the four used telescopes. Finally, we estimate that the timing noise of individual pulsars reduces the sensitivity of this data set to an isotropic, stochastic GW background by a factor of >9.1 and by a factor of >2.3 for continuous GWs from resolvable, inspiralling supermassive black-hole binaries with circular orbits., Comment: Accepted for publication by the Monthly Notices of the Royal Astronomical Society

Gravitational waves are expected to be radiated by supermassive black hole binaries formed during galaxy mergers. A stochastic superposition of gravitational waves from all such binary systems will modulate the arrival times of pulses from radio pulsars. Using observations of millisecond pulsars obtained with the Parkes radio telescope, we constrain the characteristic amplitude of this background, $A_{\rm c,yr}$, to be < $1.0\times10^{-15}$ with 95% confidence. This limit excludes predicted ranges for $A_{\rm c,yr}$ from current models with 91-99.7% probability. We conclude that binary evolution is either stalled or dramatically accelerated by galactic-center environments, and that higher-cadence and shorter-wavelength observations would result in an increased sensitivity to gravitational waves., Comment: Published in Science 25 September 2015

We extend the recently introduced Bayesian framework `Generative Pulsar Timing Analysis' to incorporate both pulse jitter (high frequency variation in the arrival time of the pulse) and epoch to epoch stochasticity in the shape of the pulse profile. This framework allows for a full timing analysis to be performed on the folded profile data, rather than the site arrival times as is typical in most timing studies. We apply this extended framework both to simulations, and to an 11 yr, 10 cm data set for PSR J1909$-$3744. Using simulations, we show that temporal profile variation can induce timing noise in the residuals that when performing a standard timing analysis is highly covariant with the signal expected from a gravitational wave (GW) background. When working in the profile domain, these variations are de-correlated from the expected GW signal, resulting in significant improvement in the obtained upper limits. Using the PSR J1909$-$3744 data set from the Parkes Pulsar Timing Array project, we find significant evidence for systematic high-frequency profile variation resulting from non-Gaussian noise in the oldest observing system, but no evidence for either detectable pulse jitter, or low-frequency profile shape variation. Using our profile domain framework we therefore obtain upper limits on a red noise process with a spectral index of $\gamma = 13/3$ of $1\times10^{-15}$, consistent with previously published limits., Comment: 16 pages, 13 figures, 2 tables

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}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., Comment: 16 pages, 11 figures, accepted for publication in MNRAS

The paucity of observed supermassive black hole binaries (SMBHBs) may imply that the gravitational wave background (GWB) from this population is anisotropic, rendering existing analyses sub-optimal. We present the first constraints on the angular distribution of a nanohertz stochastic GWB from circular, inspiral-driven SMBHBs using the $2015$ European Pulsar Timing Array data [Desvignes et al. (in prep.)]. Our analysis of the GWB in the $\sim 2 - 90$ nHz band shows consistency with isotropy, with the strain amplitude in $l>0$ spherical harmonic multipoles $\lesssim 40\%$ of the monopole value. We expect that these more general techniques will become standard tools to probe the angular distribution of source populations., Comment: 6 pages, 2 figures, 1 table. Accepted for publication in Physical Review Letters

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 $\Omega_\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\mu/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\mu/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 $\Omega^\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., Comment: 24 pages, 5 tables, 17 figures

A new Bayesian method for the analysis of folded pulsar timing data is presented that allows for the simultaneous evaluation of evolution in the pulse profile in either frequency or time, along with the timing model and additional stochastic processes such as red spin noise, or dispersion measure variations. We model the pulse profiles using `shapelets' - a complete ortho-normal set of basis functions that allow us to recreate any physical profile shape. Any evolution in the profiles can then be described as either an arbitrary number of independent profiles, or using some functional form. We perform simulations to compare this approach with established methods for pulsar timing analysis, and to demonstrate model selection between different evolutionary scenarios using the Bayesian evidence. %s The simplicity of our method allows for many possible extensions, such as including models for correlated noise in the pulse profile, or broadening of the pulse profiles due to scattering. As such, while it is a marked departure from standard pulsar timing analysis methods, it has clear applications for both new and current datasets, such as those from the European Pulsar Timing Array (EPTA) and International Pulsar Timing Array (IPTA)., Comment: 11 pages, 6 Figures. arXiv admin note: text overlap with arXiv:1405.2460

The broad spectral bandwidth at mm and cm-wavelengths provided by the recent upgrades to the Karl G. Jansky Very Large Array (VLA) has made it possible to conduct unbiased searches for molecular CO line emission at redshifts, z > 1.31. We present the discovery of a gas-rich, star-forming galaxy at z = 2.48, through the detection of CO(1-0) line emission in the COLDz survey, through a sensitive, Ka-band (31 to 39 GHz) VLA survey of a 6.5 square arcminute region of the COSMOS field. We argue that the broad line (FWHM ~570 +/- 80 km/s) is most likely to be CO(1-0) at z=2.48, as the integrated emission is spatially coincident with an infrared-detected galaxy with a photometric redshift estimate of z = 3.2 +/- 0.4. The CO(1-0) line luminosity is L'_CO = (2.2 +/- 0.3) x 10^{10} K km/s pc^2, suggesting a cold molecular gas mass of M_gas ~ (2 - 8)x10^{10}M_solar depending on the assumed value of the molecular gas mass to CO luminosity ratio alpha_CO. The estimated infrared luminosity from the (rest-frame) far-infrared spectral energy distribution (SED) is L_IR = 2.5x10^{12} L_solar and the star-formation rate is ~250 M_solar/yr, with the SED shape indicating substantial dust obscuration of the stellar light. The infrared to CO line luminosity ratio is ~114+/-19 L_solar/(K km/s pc^2), similar to galaxies with similar SFRs selected at UV/optical to radio wavelengths. This discovery confirms the potential for molecular emission line surveys as a route to study populations of gas-rich galaxies in the future.

A new Bayesian method for performing an image domain search for line-emitting galaxies is presented. The method uses both spatial and spectral information to robustly determine the source properties, employing either simple Gaussian, or other physically motivated models whilst using the evidence to determine the probability that the source is real. In this paper, we describe the method, and its application to both a simulated data set, and a blind survey for cold molecular gas using observations of the Hubble Deep Field North taken with the Plateau de Bure Interferometer. We make a total of 6 robust detections in the survey, 5 of which have counterparts in other observing bands. We identify the most secure detections found in a previous investigation, while finding one new probable line source with an optical ID not seen in the previous analysis. This study acts as a pilot application of Bayesian statistics to future searches to be carried out both for low-$J$ CO transitions of high redshift galaxies using the JVLA, and at millimeter wavelengths with ALMA, enabling the inference of robust scientific conclusions about the history of the molecular gas properties of star-forming galaxies in the Universe through cosmic time., Comment: 12 pages, 4 Figures, 5 Tables

We introduce a method for performing a robust Bayesian analysis of non-Gaussianity present in pulsar timing data, simultaneously with the pulsar timing model, and additional stochastic parameters such as those describing red spin noise and dispersion measure variations. The parameters used to define the presence of non-Gaussianity are zero for Gaussian processes, giving a simple method of defining the strength of non-Gaussian behaviour. We use simulations to show that assuming Gaussian statistics when the noise in the data is drawn from a non-Gaussian distribution can significantly increase the uncertainties associated with the pulsar timing model parameters. We then apply the method to the publicly available 15 year Parkes Pulsar Timing Array data release 1 dataset for the binary pulsar J0437$-$4715. In this analysis we present a significant detection of non-Gaussianity in the uncorrelated non-thermal noise, but we find that it does not yet impact the timing model or stochastic parameter estimates significantly compared to analysis performed assuming Gaussian statistics. The methods presented are, however, shown to be of immediate practical use for current European Pulsar Timing Array (EPTA) and International Pulsar Timing Array (IPTA) datasets., Comment: 17 pages, 9 Figures

We present detections of the CO(J= 1-0) emission line in a sample of four massive star-forming galaxies at z~1.5-2.2 obtained with the Karl G. Jansky Very Large Array (VLA). Combining these observations with previous CO(2-1) and CO(3-2) detections of these galaxies, we study the excitation properties of the molecular gas in our sample sources. We find an average line brightness temperature ratios of R_{21}=0.70+\-0.16 and R_{31}=0.50+\-0.29, based on measurements for three and two galaxies, respectively. These results provide additional support to previous indications of sub-thermal gas excitation for the CO(3-2) line with a typically assumed line ratio R_{31}~0.5. For one of our targets, BzK-21000, we present spatially resolved CO line maps. At the resolution of 0.18 arcsec (1.5 kpc), most of the emission is resolved out except for some clumpy structure. From this, we attempt to identify molecular gas clumps in the data cube, finding 4 possible candidates. We estimate that <40 % of the molecular gas is confined to giant clumps (~1.5 kpc in size), and thus most of the gas could be distributed in small fainter clouds or in fairly diffuse extended regions of lower brightness temperatures than our sensitivity limit., Comment: Accepted 26 April 2014; Title modified to match published version

In order to reach the sensitivity required to detect gravitational waves, pulsar timing array experiments need to mitigate as much noise as possible in timing data. A dominant amount of noise is likely due to variations in the dispersion measure. To correct for such variations, we develop a statistical method inspired by the maximum likelihood estimator and optimal filtering. Our method consists of two major steps. First, the spectral index and amplitude of dispersion measure variations are measured via a time-domain spectral analysis. Second, the linear optimal filter is constructed based on the model parameters found in the first step, and is used to extract the dispersion measure variation waveforms. Compared to current existing methods, this method has better time resolution for the study of short timescale dispersion variations, and generally produces smaller errors in waveform estimations. This method can process irregularly sampled data without any interpolation because of its time-domain nature. Furthermore, it offers the possibility to interpolate or extrapolate the waveform estimation to regions where no data is available. Examples using simulated data sets are included for demonstration., Comment: 15 pages, 15 figures, submitted 15th Sept. 2013, accepted 2nd April 2014 by MNRAS. MNRAS, 2014

We present Karl G. Jansky Very Large Array (VLA) observations of 44 GHz continuum and CO J=2-1 line emission in BR1202-0725 at z=4.7 (a starburst galaxy and quasar pair) and BRI1335-0417 at z=4.4 (also hosting a quasar). With the full 8 GHz bandwidth capabilities of the upgraded VLA, we study the (rest-frame) 250 GHz thermal dust continuum emission for the first time along with the cold molecular gas traced by the Low-J CO line emission. The measured CO J=2-1 line luminosities of BR1202-0725 are L'(CO) = (8.7+/-0.8)x10^10 K km/s pc^2 and L'(CO) = (6.0+/-0.5)x10^10 K km/s pc^2 for the submm galaxy (SMG) and quasar, which are equal to previous measurements of the CO J=5-4 line luminosities implying thermalized line emission and we estimate a combined cold molecular gas mass of ~9x10^10 Msun. In BRI1335-0417 we measure L'(CO) = (7.3+/-0.6)x10^10 K km/s pc^2. We detect continuum emission in the SMG BR1202-0725 North (S(44GHz) = 51+/-6 microJy), while the quasar is detected with S(44GHz) = 24+/-6 microJy and in BRI1335-0417 we measure S(44GHz) = 40+/-7 microJy. Combining our continuum observations with previous data at (rest-frame) far-infrared and cm-wavelengths, we fit three component models in order to estimate the star-formation rates. This spectral energy distribution fitting suggests that the dominant contribution to the observed 44~GHz continuum is thermal dust emission, while either thermal free-free or synchrotron emission contributes less than 30%., Comment: 16 pages , 4 figures. Accepted for publication in ApJ

We present direct constraints on the CO luminosity function at high redshift and the resulting cosmic evolution of the molecular gas density, $\rho_{\rm H2}$(z), based on a blind molecular line scan in the Hubble Deep Field North (HDF-N) using the IRAM Plateau de Bure Interferometer. Our line scan of the entire 3mm window (79-115 GHz) covers a cosmic volume of ~7000 Mpc$^3$, and redshift ranges z<0.45, 1.012. We use the rich multiwavelength and spectroscopic database of the HDF-N to derive some of the best constraints on CO luminosities in high redshift galaxies to date. We combine the blind CO detections in our molecular line scan (presented in a companion paper) with stacked CO limits from galaxies with available spectroscopic redshifts (slit or mask spectroscopy from Keck and grism spectroscopy from HST) to give first blind constraints on high-z CO luminosity functions and the cosmic evolution of the H2 mass density $\rho_{\rm H2}$(z) out to redshifts z~3. A comparison to empirical predictions of $\rho_{\rm H2}$(z) shows that the securely detected sources in our molecular line scan already provide significant contributions to the predicted $\rho_{\rm H2}$(z) in the redshift bins ~1.5 and ~2.7. Accounting for galaxies with CO luminosities that are not probed by our observations results in cosmic molecular gas densities $\rho_{\rm H2}$(z) that are higher than current predictions. We note however that the current uncertainties (in particular the luminosity limits, number of detections, as well as cosmic volume probed) are significant, a situation that is about to change with the emerging ALMA observatory., Comment: Accepted for publication in ApJ

We present a molecular line scan in the Hubble Deep Field North (HDF-N) that covers the entire 3mm window (79-115 GHz) using the IRAM Plateau de Bure Interferometer. Our CO redshift coverage spans z<0.45, 12. We reach a CO detection limit that is deep enough to detect essentially all z>1 CO lines reported in the literature so far. We have developed and applied different line searching algorithms, resulting in the discovery of 17 line candidates. We estimate that the rate of false positive line detections is ~2/17. We identify optical/NIR counterparts from the deep ancillary database of the HDF-N for seven of these candidates and investigate their available SEDs. Two secure CO detections in our scan are identified with star-forming galaxies at z=1.784 and at z=2.047. These galaxies have colors consistent with the `BzK' color selection and they show relatively bright CO emission compared with galaxies of similar dust continuum luminosity. We also detect two spectral lines in the submillimeter galaxy HDF850.1 at z=5.183. We consider an additional 9 line candidates as high quality. Our observations also provide a deep 3mm continuum map (1-sigma noise level = 8.6 $\mu$Jy/beam). Via a stacking approach, we find that optical/MIR bright galaxies contribute only to <50% of the SFR density at 1

Here we present a Bayesian method of including discrete measurements of dispersion measure due to the interstellar medium in the direction of a pulsar as prior information in the analysis of that pulsar. We use a simple simulation to show the efficacy of this method, where the inclusion of the additional measurements results in both a significant increase in the precision with which the timing model parameters can be obtained, and an improved upper limit on the amplitude of any red noise in the dataset. We show that this method can be applied where no multi-frequency data exists across much of the dataset, and where there is no simultaneous multi-frequency data for any given observing epoch. Including such information in the analysis of upcoming International Pulsar Timing Array (IPTA) and European Pulsar Timing Array (EPTA) data releases could therefore prove invaluable in obtaining the most constraining limits on gravitational wave signals within those datasets., Comment: 7 pages, 1 Table, 3 Figures. arXiv admin note: substantial text overlap with arXiv:1310.2120

A new Bayesian software package for the analysis of pulsar timing data is presented in the form of TempoNest which allows for the robust determination of the non-linear pulsar timing solution simultaneously with a range of additional stochastic parameters. This includes both red spin noise and dispersion measure variations using either power law descriptions of the noise, or through a model-independent method that parameterises the power at individual frequencies in the signal. We use TempoNest to show that at noise levels representative of current datasets in the European Pulsar Timing Array (EPTA) and International Pulsar Timing Array (IPTA) the linear timing model can underestimate the uncertainties of the timing solution by up to an order of magnitude. We also show how to perform Bayesian model selection between different sets of timing model and stochastic parameters, for example, by demonstrating that in the pulsar B1937+21 both the dispersion measure variations and spin noise in the data are optimally modelled by simple power laws. Finally we show that not including the stochastic parameters simultaneously with the timing model can lead to unpredictable variation in the estimated uncertainties, compromising the robustness of the scientific results extracted from such analysis., Comment: 21 pages, 12 Figures

We present ALMA observations of the [C II] 158 micron fine structure line and dust continuum emission from the host galaxies of five redshift 6 quasars. We also report complementary observations of 250 GHz dust continuum and CO (6-5) line emission from the z=6.00 quasar SDSS J231038.88+185519.7. The ALMA observations were carried out in the extended array at 0.7" resolution. We have detected the line and dust continuum in all five objects. The derived [C II] line luminosities are 1.6x10^{9} to 8.8x10^{9} Lsun and the [C II]-to-FIR luminosity ratios are 3.0-5.6x10^{-4}, which is comparable to the values found in other high-redshift quasar-starburst systems and local ultra-luminous infrared galaxies. The sources are marginally resolved and the intrinsic source sizes (major axis FWHM) are constrained to be 0.3" to 0.6" (i.e., 1.7 to 3.5 kpc) for the [C II] line emission and 0.2" to 0.4" (i.e., 1.2 to 2.3 kpc) for the continuum. These measurements indicate that there is vigorous star formation over the central few kpc in the quasar host galaxies. The ALMA observations also constrain the dynamical properties of the atomic gas in the starburst nuclei. The intensity-weighted velocity maps of three sources show clear velocity gradients. Such velocity gradients are consistent with a rotating, gravitationally bound gas component, although they are not uniquely interpreted as such. Under the simplifying assumption of rotation, the implied dynamical masses within the [C II]-emitting regions are of order 10^{10} to 10^{11} Msun. Given these estimates, the mass ratios between the SMBHs and the spheroidal bulge are an order of magnitude higher than the mean value found in local spheroidal galaxies, which is in agreement with results from previous CO observations of high redshift quasars., Comment: 25 pages, 5 figures, accepted for publication in ApJ

We present further analysis of the [CII] 158$\mu$m fine structure line and thermal dust continuum emission from the archetype extreme starburst/AGN group of galaxies in the early Universe, BRI 1202-0725 at $z=4.7$, using the Atacama Large Millimeter Array. The group is long noted for having a closely separated (26kpc in projection) FIR-hyperluminous quasar host galaxy and an optically obscured submm galaxy (SMG). A short ALMA test observation reveals a rich laboratory for the study of the myriad processes involved in clustered massive galaxy formation in the early Universe. Strong [CII] emission from the SMG and the quasar have been reported earlier by Wagg et al. (2012) based on these observations. In this letter, we examine in more detail the imaging results from the ALMA observations, including velocity channel images, position-velocity plots, and line moment images. We present detections of [CII] emission from two Ly$\alpha$-selected galaxies in the group, demonstrating the relative ease with which ALMA can detect the [CII] emission from lower star formation rate galaxies at high redshift. Imaging of the [CII] emission shows a clear velocity gradient across the SMG, possibly indicating rotation or a more complex dynamical system on a scale $\sim 10$kpc. There is evidence in the quasar spectrum and images for a possible outflow toward the southwest, as well as more extended emission (a 'bridge'), between the quasar and the SMG, although the latter could simply be emission from Ly$\alpha$-1 blending with that of the quasar at the limited spatial resolution of the current observations. These results provide an unprecedented view of a major merger of gas rich galaxies driving extreme starbursts and AGN accretion during the formation of massive galaxies and supermassive black holes within 1.3 Gyr of the Big Bang., Comment: for the astrophysical journal. 5 figs

We use sensitive observations of three high redshift sources; [CII] fine structure and CO(2-1) rotational transitions for the z=6.4 Quasar host galaxy (QSO) J1148+5251, and [CII] and CO(5-4) transitions from the QSO BR1202-0725 and its sub-millimeter companion (SMG) galaxy at z=4.7. We use these observations to place constraints on the quantity Dz = z(CO) - z(CII) for each source where z(CO) and z(CII) are the observed redshifts of the CO rotational transition and [CII] fine structure transition respectively, using a combination of approaches; 1) By modelling the emission line profiles using `shapelets' to compare both the emission redshifts and the line profiles themselves, in order to make inferences about the intrinsic velocity differences between the molecular and atomic gas, and 2) By performing a marginalisation over all model parameters in order to calculate a non-parametric estimate of Dz. We derive 99% confidence intervals for the marginalised posterior of Dz of (-1.9 pm 1.3) x10^-3, (-3 pm 8) x10^-4 and (-2 pm 4) x10^-3 for J1148+5251, and the BR1202-0725 QSO and SMG respectively. We show the [CII] and CO(2-1) line profiles for J1148+5251 are consistent with each other within the limits of the data, whilst the [CII] and CO(5-4) line profiles from the BR1202-0725 QSO and SMG respectively have 65 and >99.9% probabilities of being inconsistent, with the CO(5-4) lines ~ 30% wider than the [CII] lines. Therefore whilst the observed values of Dz can correspond to variations in the quantity Delta F/F with cosmic time, where F=alpha^2/mu, with alpha the fine structure constant, and mu the proton-to-electron mass ratio, of both (-3.3 pm 2.3) x10^-4 for a look back time of 12.9 Gyr and of (-5 pm 15) x10^-5 for a look back time of 12.4 Gyr we propose that they are the result of the two species of gas being spatially separated as indicated by the inconsistencies in their line profiles., Comment: 11 pages, 10 figures

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., Comment: 20 pages, 13 figures, 6 tables. Updated with additional studies to reconstruct the Hellings and Downs curve, and to recover the GW-background spectrum. Submitted to PRD

A new model independent method is presented for the analysis of pulsar timing data and the estimation of the spectral properties of an isotropic gravitational wave background (GWB). We show that by rephrasing the likelihood we are able to eliminate the most costly aspects of computation normally associated with this type of data analysis. When applied to the International Pulsar Timing Array Mock Data Challenge data sets this results in speedups of approximately 2 to 3 orders of magnitude compared to established methods. We present three applications of the new likelihood. In the low signal to noise regime we sample directly from the power spectrum coefficients of the GWB signal realization. In the high signal to noise regime, where the data can support a large number of coefficients, we sample from the joint probability density of the power spectrum coefficients for the individual pulsars and the GWB signal realization. Critically in both these cases we need make no assumptions about the form of the power spectrum of the GWB, or the individual pulsars. Finally we present a method for characterizing the spatial correlation between pulsars on the sky, making no assumptions about the form of that correlation, and therefore providing the only truly general Bayesian method of confirming a GWB detection from pulsar timing data., Comment: 9 pages, 4 figures

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., Comment: 6 pages, 1 figure, 1 table

We present Karl G. Jansky Very Large Array (VLA) observations of the CO(2-1) emission in the z=4.05 submillimeter galaxy (SMG) GN20. These high-resolution data allow us to image the molecular gas at 1.3 kpc resolution just 1.6 Gyr after the Big Bang. The data reveal a clumpy, extended gas reservoir, 14 +/- 4 kpc in diameter, in unprecedented detail. A dynamical analysis shows that the data are consistent with a rotating disk of total dynamical mass 5.4 +/- 2.4 X 10^11 M_sun. We use this dynamical mass estimate to constrain the CO-to-H_2 mass conversion factor (alpha_CO), finding alpha_CO=1.1 +/- 0.6 M_sun (K km s^-1 pc^2)^-1. We identify five distinct molecular gas clumps in the disk of GN20 with masses a few percent of the total gas mass, brightness temperatures of 16-31K, and surface densities of >3,200-4,500 X (alpha_CO/0.8) M_sun pc^-2. Virial mass estimates indicate they could be self-gravitating, and we constrain their CO-to-H_2 mass conversion factor to be <0.2-0.7 M_sun (K km s^-1 pc^2)^-1. A multiwavelength comparison demonstrates that the molecular gas is concentrated in a region of the galaxy that is heavily obscured in the rest-frame UV/optical. We investigate the spatially-resolved gas excitation and find that the CO(6-5)/CO(2-1) ratio is constant with radius, consistent with star formation occuring over a large portion of the disk. We discuss the implications of our results in the context of different fueling scenarios for SMGs., Comment: 15 pages, 9 figures, accepted for publication in ApJ

We report results from a deep Jansky Very Large Array (JVLA) search for CO 1-0 line emission from galaxies in a candidate galaxy cluster at z~1.55 in the COSMOS field. We target 4 galaxies with optical spectroscopic redshifts in the range z=1.47-1.59. Two of these 4 galaxies, ID51613 and ID51813, are nominally detected in CO line emission at the 3-4 sigma level. We find CO luminosities of 2.4x10^10 K km/s pc^2 and 1.3x10^10 K km/s pc^2, respectively. Taking advantage from the clustering and 2-GHz bandwidth of the JVLA, we perform a search for emission lines in the proximity of optical sources within the field of view of our observations. We limit our search to galaxies with K<23.5 (AB) and z_phot=1.2-1.8. We find 2 bright optical galaxies to be associated with significant emission line peaks (>4 sigma) in the data cube, which we identify with the CO line emission. To test the reliability of the line peaks found, we performed a parallel search for line peaks using a Bayesian inference method. Monte Carlo simulations show that such associations are statistically significant, with probabilities of chance association of 3.5% and 10.7% for ID 51207 and ID 51380, respectively. Modeling of their optical/IR SEDs indicates that the CO detected galaxies and candidates have stellar masses and SFRs in the range (0.3-1.1)x10^11 M_sun and 60-160 M_sun/yr, with SFEs comparable to that found in other star-forming galaxies at similar redshifts. By comparing the space density of CO emitters derived from our observations with the space density derived from previous CO detections at z~1.5, and with semi-analytic predictions for the CO luminosity function, we suggest that the latter tend to underestimate the number of CO galaxies detected at high-redshift. Finally, we argue about the benefits of future blind CO searches in clustered fields with upcoming submm/radio facilities., Comment: Accepted for publication in MNRAS. Abstract has been slightly shortened compared to original pdf version

The Hubble Deep Field (HDF) is a region in the sky that provides one of the deepest multi-wavelength views of the distant universe and has led to the detection of thousands of galaxies seen throughout cosmic time. An early map of the HDF at a wavelength of 850 microns that is sensitive to dust emission powered by star formation revealed the brightest source in the field, dubbed HDF850.1. For more than a decade, this source remained elusive and, despite significant efforts, no counterpart at shorter wavelengths, and thus no redshift, size or mass, could be identified. Here we report, using a millimeter wave molecular line scan, an unambiguous redshift determination for HDF850.1 of z=5.183. This places HDF850.1 in a galaxy overdensity at z~5.2 in the HDF, corresponding to a cosmic age of only 1.1 Gyr after the Big Bang. This redshift is significantly higher than earlier estimates and higher than most of the >100 sub-millimeter bright galaxies identified to date. The source has a star formation rate of 850 M_sun/yr and is spatially resolved on scales of 5 kpc, with an implied dynamical mass of ~1.3x10^11 M_sun, a significant fraction of which is present in the form of molecular gas. Despite our accurate redshift and position, a counterpart arising from starlight remains elusive., Comment: Nature, in press