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.

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.

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.

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.

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.

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.

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.

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.

Supermassive Binary Black Holes (SMBBHs) with sub-pc separations form in the course of galaxy mergers, if both galaxies harbour massive black holes. Clear observational evidence for them however still eludes us. We propose a novel method of identifying these systems by means of reverberation mapping their circumbinary disk after a tidal disruption event has ionized it. The tidal disruption of a star at the secondary leads to strong asymmetries in the disk response. We model the shape of the velocity--delay maps for various toy disk models and more realistic gas distributions obtained by SPH simulations. The emissivity of the ionized disk is calculated with {\em Cloudy}. We find peculiar asymmetries in the maps for off center ionizing sources that may help us constrain geometrical parameters of a circumbinary disk such as semi-major axis and orbital phase of the secondary, as well as help strengthen the observational evidence for sub-parsec SMBBHs as such.

Dense stellar systems such as globular clusters and dense nuclear clusters are the breeding ground of sources of gravitational waves for the advanced detectors LIGO and Virgo. These systems deserve a close study to estimate rates and parameter distribution. This is not an easy task, since the evolution of a dense stellar cluster involves the integration of $N$ bodies with high resolution in time and space and including hard binaries and their encounters and, in the case of gravitational waves (GWs), one needs to take into account important relativistic corrections. In this work we present the first implementation of the effect of spin in mergers in a direct-summation code, NBODY6. We employ non-spinning post-Newtonian corrections to the Newtonian accelerations up to 3.5 post-Newtonian (PN) order as well as the spin-orbit coupling up to next-to-lowest order and the lowest order spin-spin coupling. We integrate spin precession and add a consistent treatment of mergers. We analyse the implementation by running a set of two-body experiments and then we run a set of 500 simulations of a relativistic stellar cluster. In spite of the large number of mergers in our tests, the application of the algorithm is robust. We find in particular the formation of a runaway star whose spin decays with the mass it wins, independently of the initial value of the spins of the stars. More remarkably, the subset of compact objects that do not undergo many mergers, and hence represent a more realistic system, has a correlation between the final absolute spin and the initial choice for the initial distribution, which could provide us with information about the evolution of spins in dense clusters once the first detections have started.

Hubble Space Telescope observations of globular clusters (GCs) in the Antennae galaxy show clusters of clusters, or regions in the galaxy that span hundreds of parsecs, where many of the GCs are doomed to collide, and eventually merge. Several such objects appear likely to present a significant range in ages, hence possibly metallicities, and their merger could plausibly lead to multi-metallic GCs. Here we explore this process with direct-summation N-body simulations with GPU hardware. Our results reveal that colliding GCs with different metallicities and ages can produce a GC with multiplicity and occupation fractions not unlike those observed in multi-metallic clusters. In our simulations, the merged clusters have a phase with a larger amount of flattening than average, as a consequence of rapid rotation- thus suggesting that relatively recent mergers may play a role in producing highly flattened, multi-metallic clusters. We additionally explore the role of the King parameter of the cluster in the occupation fractions with a set of 160 direct-summation simulations and find that for equal size clusters the King parameter of the progenitor clusters determines the occupation fractions in the merger product, while in unequal size mergers the size of the clusters dominates the distribution of stars in the new GC. In particular, we find that the observed distribution of populations in Omega Cen can be described to some extent with our dynamical models.

In our current interpretation of the hierarchical structure of the universe it is well established that galaxies collide and merge with each other during their lifetime. If massive black holes (MBHs) reside in galactic centres, we expect them to form binaries in galactic nuclei surrounded by a circumbinary disc. If cooling is efficient enough, the gas in the disc will clump and trigger stellar formation in situ. In this first paper we address the evolution of the binary under the influence of the newly formed stars, which form individually and also clustered. We use SPH techniques to evolve the gas in the circumbinary disc and to study the phase of star formation. When the amount of gas in the disc is negligible, we further evolve the system with a high-accurate direct-summation $N-$body code to follow the evolution of the stars, the innermost binary and tidal disruption events (TDEs). For this, we modify the direct N-body code to (i) include treatment of TDEs and to (ii) include "gas cloud particles" that mimic the gas, so that the stellar clusters do not disolve when we follow their infall on to the MBHs. We find that the amount of stars disrupted by either infalling stellar clusters or individual stars is as large as 10^{-4}/yr per binary, higher than expected for typical galaxies.

The capture of a compact object in a galactic nucleus by a massive black hole (MBH) is the best way to map space and time around it. Compact objects such as stellar black holes on a capture orbit with a very high eccentricity have been wrongly assumed to be lost for the system after an intense burst of radiation, which has been described as a "direct plunge". We prove that these very eccentric capture orbits spend actually a similar number of cycles in a LISA-like detector as those with lower eccentricities if the central MBH is spinning. Although the rates are higher for high-eccentricity EMRIs, the spin also enhances the rates of lower-eccentricity EMRIs. This last kind have received more attention because of the fact that high-eccentricity EMRIs were thought to be direct plunges and thus negligible. On the other hand, recent work on stellar dynamics has demonstrated that there seems to be a complot in phase space acting on these lower-eccentricity captures, since their rates decrease significantly by the presence of a blockade in the rate at which orbital angular momenta change takes place. This so-called "Schwarzschild barrier" is a result of the impact of relativistic precession on to the stellar potential torques, and thus it affects the enhancement on lower-eccentricity EMRIs that one would expect from resonant relaxation. We confirm and quantify the existence of this barrier using a statitical sample of 2,500 direct-summation N-body simulations using both a post-Newtonian but also, and for the first time, a geodesic approximation for the relativistic orbits. The existence of the barrier prevents "traditional EMRIs" from approaching the central MBH, but if the central MBH is spinning the rate will be anyway dominated by highly-eccentric extreme-mass ratio inspirals, which insolently ignore the presence of the barrier, because they are driven by two-body relaxation.

The capture of a compact object in a galactic nucleus by a massive black hole (MBH), an extreme-mass ratio inspiral (EMRI), is the best way to map space and time around it. Recent work on stellar dynamics has demonstrated that there seems to be a complot in phase space acting on low-eccentricity captures, since their rates decrease significantly by the presence of a blockade in the rate at which orbital angular momenta change takes place. This so-called "Schwarzschild barrier" is a result of the impact of relativistic precession on to the stellar potential torques, and thus it affects the enhancement on lower-eccentricity EMRIs that one would expect from resonant relaxation. We confirm and quantify the existence of this barrier using a statistical sample of 2,500 direct-summation N-body simulations using both a post-Newtonian and also for the first time in a direct-summation code a geodesic approximation for the relativistic orbits. The existence of the barrier prevents low-eccentricity EMRIs from approaching the central MBH, but high-eccentricity EMRIs, which have been wrongly classified as "direct plunges" until recently, ignore the presence of the barrier, because they are driven by two-body relaxation. Hence, since the rates are significantly affected in the case of low-eccentricity EMRIs, we predict that a LISA-like observatory such as eLISA will predominantly detect high-eccentricity EMRIs.

Measurements of gravitational waves from the inspiral of a stellar-mass compact object into a massive black hole are unique probes to test General Relativity (GR) and MBH properties, as well as the stellar distribution about these holes in galactic nuclei. Current data analysis techniques can provide us with parameter estimation with very narrow errors. However, an EMRI is not a two-body problem, since other stellar bodies orbiting nearby will influence the capture orbit. Any deviation from the isolated inspiral will induce a small, though observable deviation from the idealised waveform which could be misinterpreted as a failure of GR. Based on conservative analysis of mass segregation in a Milky Way like nucleus, we estimate that the possibility that another star has a semi-major axis comparable to that of the EMRI is non-negligible, although probably very small. This star introduces an observable perturbation in the orbit in the case in which we consider only loss of energy via gravitational radiation. When considering the two first-order non-dissipative post-Newtonian contributions (the periapsis shift of the orbit), the evolution of the orbital elements of the EMRI turns out to be chaotic in nature. The implications of this study are twofold. From the one side, the application to testing GR and measuring MBHs parameters with the detection of EMRIs in galactic nuclei with a millihertz mission will be even more challenging than believed. From the other side, this behaviour could in principle be used as a signature of mass segregation in galactic nuclei.