Source author record

A. Sesana

A. Sesana 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

30works
8topics
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

30 published item(s)

preprint2022arXiv

The International Pulsar Timing Array second data release: Search for an isotropic Gravitational Wave Background

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})^α$, 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 $α= -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 $α= -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.

preprint2021arXiv

Multi-messenger-Athena Synergy White Paper

In this paper we explore the scientific synergies between Athena and some of the key multi-messenger facilities that should be operative concurrently with Athena. These facilities include LIGO A+, Advanced Virgo+ and future detectors for ground-based observation of gravitational waves (GW), LISA for space-based observations of GW, IceCube and KM3NeT for neutrino observations, and CTA for very high energy observations. These science themes encompass pressing issues in astrophysics, cosmology and fundamental physics such as: the central engine and jet physics in compact binary mergers, accretion processes and jet physics in Super-Massive Binary Black Holes (SMBBHs) and in compact stellar binaries, the equation of state of neutron stars, cosmic accelerators and the origin of Cosmic Rays (CRs), the origin of intermediate and high-Z elements in the Universe, the Cosmic distance scale and tests of General Relativity and the Standard Model. Observational strategies for implementing the identified science topics are also discussed. A significant part of the sources targeted by multi-messenger facilities is of transient nature. We have thus also discussed the synergy of \textsl{Athena} with wide-field high-energy facilities, taking THESEUS as a case study for transient discovery. This discussion covers all the Athena science goals that rely on follow-up observations of high-energy transients identified by external observatories, and includes also topics that are not based on multi-messenger observations, such as the search for missing baryons or the observation of early star populations and metal enrichment at the cosmic dawn with Gamma-Ray Bursts (GRBs).

preprint2019arXiv

A pulsar-based timescale from the International Pulsar Timing Array

We have constructed a new timescale, TT(IPTA16), based on observations of radio pulsars presented in the first data release from the International Pulsar Timing Array (IPTA). We used two analysis techniques with independent estimates of the noise models for the pulsar observations and different algorithms for obtaining the pulsar timescale. The two analyses agree within the estimated uncertainties and both agree with TT(BIPM17), a post-corrected timescale produced by the Bureau International des Poids et Mesures (BIPM). We show that both methods could detect significant errors in TT(BIPM17) if they were present. We estimate the stability of the atomic clocks from which TT(BIPM17) is derived using observations of four rubidium fountain clocks at the US Naval Observatory. Comparing the power spectrum of TT(IPTA16) with that of these fountain clocks suggests that pulsar-based timescales are unlikely to contribute to the stability of the best timescales over the next decade, but they will remain a valuable independent check on atomic timescales. We also find that the stability of the pulsar-based timescale is likely to be limited by our knowledge of solar-system dynamics, and that errors in TT(BIPM17) will not be a limiting factor for the primary goal of the IPTA, which is to search for the signatures of nano-Hertz gravitational waves.

preprint2016arXiv

From Spin Noise to Systematics: Stochastic Processes in the First International Pulsar Timing Array Data Release

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.

preprint2016arXiv

High-precision timing of 42 millisecond pulsars with the European Pulsar Timing Array

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.

preprint2016arXiv

Selection bias in dynamically-measured super-massive black hole samples: consequences for pulsar timing arrays

Supermassive black hole -- host galaxy relations are key to the computation of the expected gravitational wave background (GWB) in the pulsar timing array (PTA) frequency band. It has been recently pointed out that standard relations adopted in GWB computations are in fact biased-high. We show that when this selection bias is taken into account, the expected GWB in the PTA band is a factor of about three smaller than previously estimated. Compared to other scaling relations recently published in the literature, the median amplitude of the signal at $f=1$yr$^{-1}$ drops from $1.3\times10^{-15}$ to $4\times10^{-16}$. Although this solves any potential tension between theoretical predictions and recent PTA limits without invoking other dynamical effects (such as stalling, eccentricity or strong coupling with the galactic environment), it also makes the GWB detection more challenging.

preprint2016arXiv

The International Pulsar Timing Array: First Data Release

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.

preprint2016arXiv

The noise properties of 42 millisecond pulsars from the European Pulsar Timing Array and their impact on gravitational wave searches

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.

preprint2015arXiv

Infalling clouds onto super-massive black hole binaries - I. Formation of discs, accretion and gas dynamics

There is compelling evidence that most -if not all- galaxies harbour a super-massive black hole (SMBH) at their nucleus, hence binaries of these massive objects are an inevitable product of the hierarchical evolution of structures in the universe, and represent an important but thus-far elusive phase of galaxy evolution. Gas accretion via a circumbinary disc is thought to be important for the dynamical evolution of SMBH binaries, as well as in producing luminous emission that can be used to infer their properties. One plausible source of the gaseous fuel is clumps of gas formed due to turbulence and gravitational instabilities in the interstellar medium, that later fall toward and interact with the binary. In this context, we model numerically the evolution of turbulent clouds in near-radial infall onto equal-mass SMBH binaries, using a modified version of the SPH code GADGET-3. We present a total of 12 simulations that explore different possible pericentre distances and relative inclinations, and show that the formation of circumbinary discs and discs around each SMBH ('mini-discs') depend on those parameters. We also study the dynamics of the formed discs, and the variability of the feeding rate onto the SMBHs in the different configurations.

preprint2015arXiv

Limits on anisotropy in the nanohertz stochastic gravitational-wave background

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.

preprint2014arXiv

Gravitational wave astronomy with the SKA

On a time scale of years to decades, gravitational wave (GW) astronomy will become a reality. Low frequency (nanoHz) GWs are detectable through long-term timing observations of the most stable pulsars. Radio observatories worldwide are currently carrying out observing programmes to detect GWs, with data sets being shared through the International Pulsar Timing Array project. One of the most likely sources of low frequency GWs are supermassive black hole binaries (SMBHBs), detectable as a background due to a large number of binaries, or as continuous or burst emission from individual sources. No GW signal has yet been detected, but stringent constraints are already being placed on galaxy evolution models. The SKA will bring this research to fruition. In this chapter, we describe how timing observations using SKA1 will contribute to detecting GWs, or can confirm a detection if a first signal already has been identified when SKA1 commences observations. We describe how SKA observations will identify the source(s) of a GW signal, search for anisotropies in the background, improve models of galaxy evolution, test theories of gravity, and characterise the early inspiral phase of a SMBHB system. We describe the impact of the large number of millisecond pulsars to be discovered by the SKA; and the observing cadence, observation durations, and instrumentation required to reach the necessary sensitivity. We describe the noise processes that will influence the achievable precision with the SKA. We assume a long-term timing programme using the SKA1-MID array and consider the implications of modifications to the current design. We describe the possible benefits from observations using SKA1-LOW. Finally, we describe GW detection prospects with SKA1 and SKA2, and end with a description of the expectations of GW astronomy.

preprint2014arXiv

NANOGrav Limits on Gravitational Waves from Individual Supermassive Black Hole Binaries in Circular Orbits

The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) project currently observes 43 pulsars using the Green Bank and Arecibo radio telescopes. In this work we use a subset of 17 pulsars timed for a span of roughly five years (2005--2010). We analyze these data using standard pulsar timing models, with the addition of time-variable dispersion measure and frequency-variable pulse shape terms. Within the timing data, we perform a search for continuous gravitational waves from individual supermassive black hole binaries in circular orbits using robust frequentist and Bayesian techniques. We find that there is no evidence for the presence of a detectable continuous gravitational wave; however, we can use these data to place the most constraining upper limits to date on the strength of such gravitational waves. Using the full 17 pulsar dataset we place a 95% upper limit on the sky-averaged strain amplitude of $h_0\lesssim 3.8\times 10^{-14}$ at a frequency of 10 nHz. Furthermore, we place 95% \emph{all sky} lower limits on the luminosity distance to such gravitational wave sources finding that the $d_L \gtrsim 425$ Mpc for sources at a frequency of 10 nHz and chirp mass $10^{10}{\rm M}_{\odot}$. We find that for gravitational wave sources near our best timed pulsars in the sky, the sensitivity of the pulsar timing array is increased by a factor of $\sim$4 over the sky-averaged sensitivity. Finally we place limits on the coalescence rate of the most massive supermassive black hole binaries.

preprint2013arXiv

Gravitational wave emission from binary supermassive black holes

Massive black hole binaries (MBHBs) are unavoidable outcomes of the hierarchical structure formation process, and according to the theory of general relativity are expected to be the loudest gravitational wave (GW) sources in the Universe. In this article I provide a broad overview of MBHBs as GW sources. After reviewing the basics of GW emission from binary systems and of MBHB formation, evolution and dynamics, I describe in some details the connection between binary properties and the emitted gravitational waveform. Direct GW observations will provide an unprecedented wealth of information about the physical nature and the astrophysical properties of these extreme objects, allowing to reconstruct their cosmic history, dynamics and coupling with their dense stellar and gas environment. In this context I describe ongoing and future efforts to make a direct detection with space based interferometry and pulsar timing arrays, highlighting the invaluable scientific payouts of such enterprises.

preprint2013arXiv

Insights on the astrophysics of supermassive black hole binaries from pulsar timing observations

Pulsar timing arrays (PTAs) are designed to detect the predicted gravitational wave (GW) background produced by a cosmological population of supermassive black hole (SMBH) binaries. In this contribution I review the physics of such GW background, highlighting its dependence on the overall binary population, the relation between SMBHs and their hosts, and their coupling with the stellar and gaseous environment. The latter is particularly relevant when it drives the binaries to extreme eccentricities (e>0.9), which might be the case for stellar-driven systems. This causes a substantial suppression of the low frequency signal, potentially posing a serious threat to the effectiveness of PTA observations. A future PTA detection will allow to directly observe for the first time subparsec SMBH binaries on their way to the GW driven coalescence, providing important answers of the outstanding questions related to the physics underlying the formation and evolution of these spectacular sources.

preprint2013arXiv

Pathway to the Square Kilometre Array - The German White Paper -

The Square Kilometre Array (SKA) is the most ambitious radio telescope ever planned. With a collecting area of about a square kilometre, the SKA will be far superior in sensitivity and observing speed to all current radio facilities. The scientific capability promised by the SKA and its technological challenges provide an ideal base for interdisciplinary research, technology transfer, and collaboration between universities, research centres and industry. The SKA in the radio regime and the European Extreme Large Telescope (E-ELT) in the optical band are on the roadmap of the European Strategy Forum for Research Infrastructures (ESFRI) and have been recognised as the essential facilities for European research in astronomy. This "White Paper" outlines the German science and R&D interests in the SKA project and will provide the basis for future funding applications to secure German involvement in the Square Kilometre Array.

preprint2013arXiv

Pulsar Timing Array Analysis for Black Hole Backgrounds

An astrophysical population of supermassive black hole binaries is thought to be the strongest source of gravitational waves in the frequency range covered by Pulsar Timing Arrays (PTAs). A potential cause for concern is that the standard cross-correlation method used in PTA data analysis assumes that the signals are isotropically distributed and Gaussian random, while the signals from a black hole population are likely to be anisotropic and deterministic. Here we argue that while the conventional analysis is not optimal, it is not hopeless either, as the standard Hellings-Downs correlation curve turns out to hold for point sources, and the small effective number of signal samples blurs the distinction between Gaussian and deterministic signals. Possible improvements to the standard cross-correlation analysis that account for the anisotropy of the signal are discussed.

preprint2013arXiv

The Gravitational Universe

The last century has seen enormous progress in our understanding of the Universe. We know the life cycles of stars, the structure of galaxies, the remnants of the big bang, and have a general understanding of how the Universe evolved. We have come remarkably far using electromagnetic radiation as our tool for observing the Universe. However, gravity is the engine behind many of the processes in the Universe, and much of its action is dark. Opening a gravitational window on the Universe will let us go further than any alternative. Gravity has its own messenger: Gravitational waves, ripples in the fabric of spacetime. They travel essentially undisturbed and let us peer deep into the formation of the first seed black holes, exploring redshifts as large as z ~ 20, prior to the epoch of cosmic re-ionisation. Exquisite and unprecedented measurements of black hole masses and spins will make it possible to trace the history of black holes across all stages of galaxy evolution, and at the same time constrain any deviation from the Kerr metric of General Relativity. eLISA will be the first ever mission to study the entire Universe with gravitational waves. eLISA is an all-sky monitor and will offer a wide view of a dynamic cosmos using gravitational waves as new and unique messengers to unveil The Gravitational Universe. It provides the closest ever view of the early processes at TeV energies, has guaranteed sources in the form of verification binaries in the Milky Way, and can probe the entire Universe, from its smallest scales around singularities and black holes, all the way to cosmological dimensions.

preprint2012arXiv

Detecting massive black hole binaries and unveiling their cosmic history with gravitational wave observations

Space based gravitational wave astronomy will open a completely new window on the Universe and massive black holes binaries are expected to be among the primary actors on this upcoming stage. The New Gravitational-wave Observatory (NGO) is a space interferometer proposal derived from the former Laser Interferometer Space Antenna (LISA) concept. We describe here its capabilities of observing massive black hole binaries throughout the Universe, measuring their relevant parameters (masses, spins, distance to the observer) to high precision. The statistical properties of the population of detected systems can be used to constrain the massive black hole cosmic history, providing deep insights into the faint, high redshift Universe.

preprint2012arXiv

Systematic investigation of the expected gravitational wave signal from supermassive black hole binaries in the pulsar timing band

In this letter we carry out the first systematic investigation of the expected gravitational wave (GW) background generated by supermassive black hole (SMBH) binaries in the nHz frequency band accessible to pulsar timing arrays (PTAs). We take from the literature several estimates of the redshift dependent galaxy mass function and of the fraction of close galaxy pairs to derive a wide range of galaxy merger rates. We then exploit empirical black hole-host relations to populate merging galaxies with SMBHs. The result of our procedure is a collection of a large number of phenomenological SMBH binary merger rates consistent with current observational constraints on the galaxy assembly at z<1.5. For each merger rate we compute the associated GW signal, eventually producing a large set of estimates of the nHz GW background that we use to infer confidence intervals of its expected amplitude. When considering the most recent SMBH-host relations, accounting for ultra-massive black holes in brightest cluster galaxies, we find that the nominal $1σ$ interval of the expected GW signal is only a factor of 3-to-10 below current PTA limits, implying a non negligible chance of detection in the next few years.

preprint2011arXiv

A practical guide to the massive black hole cosmic history

I review our current understanding of massive black hole (MBH) formation and evolution along the cosmic history. After a brief introductory overview of the relevance of MBHs in the hierarchical structure formation paradigm, I discuss the main viable channels for seed BH formation at high redshift and for their subsequent mass growth and spin evolution. The emerging hierarchical picture, where MBHs grow through merger triggered accretion episodes, acquiring their mass while shining as quasars, is overall robust, but too simplistic to explain the diversity observed in MBH phenomenology. I briefly discuss which future observations will help to shed light on the MBH cosmic history in the near future, paying particular attention to the upcoming gravitational wave window.

preprint2011arXiv

Blindly detecting orbital modulations of jets from merging supermassive black holes

In the last few years before merger, supermassive black hole binaries will rapidly inspiral and precess in a magnetic field imposed by a surrounding circumbinary disk. Multiple simulations suggest this relative motion will convert some of the local energy to a Poynting-dominated outflow, with a luminosity 10^{43} erg/s * (B/10^4 G)^2(M/10^8 Msun)^2 (v/0.4 c)^2, some of which may emerge as synchrotron emission at frequencies near 1 GHz where current and planned wide-field radio surveys will operate. On top of a secular increase in power on the gravitational wave inspiral timescale, orbital motion will produce significant, detectable modulations, both on orbital periods and (if black hole spins are not aligned with the binary's total angular momenta) spin-orbit precession timescales. Because the gravitational wave merger time increases rapidly with separation, we find vast numbers of these transients are ubiquitously predicted, unless explicitly ruled out (by low efficiency $ε$) or obscured (by accretion geometry f_{geo}). If the fraction of Poynting flux converted to radio emission times the fraction of lines of sight accessible $f_{geo}$ is sufficiently large (f_{geo} ε> 2\times 10^{-4} for a 1 year orbital period), at least one event is accessible to future blind surveys at a nominal 10^4 {deg}^2 with 0.5 mJy sensitivity. Our procedure generalizes to other flux-limited surveys designed to investigate EM signatures associated with many modulations produced by merging SMBH binaries.

preprint2011arXiv

Cosmography with strong lensing of LISA gravitational wave sources

LISA might detect gravitational waves from mergers of massive black hole binaries strongly lensed by intervening galaxies (Sereno et al. 2010). The detection of multiple gravitational lensing events would provide a new tool for cosmography. Constraints on cosmological parameters could be placed by exploiting either lensing statistics of strongly lensed sources or time delay measurements of lensed gravitational wave signals. These lensing methods do not need the measurement of the redshifts of the sources and the identification of their electromagnetic counterparts. They would extend cosmological probes to redshift z <= 10 and are then complementary to other lower or higher redshift tests, such as type Ia supernovae or cosmic microwave background. The accuracy of lensing tests strongly depends on the formation history of the merging binaries, and the related number of total detectable multiple images. Lensing amplification might also help to find the host galaxies. Any measurement of the source redshifts would allow to exploit the distance-redshift test in combination with lensing methods. Time-delay analyses might measure the Hubble parameter H_0 with accuracy of >= 10 km s^{-1}Mpc^{-1}. With prior knowledge of H_0, lensing statistics and time delays might constrain the dark matter density (delta Omega_M >= 0.08, due to parameter degeneracy). Inclusion of our methods with other available orthogonal techniques might significantly reduce the uncertainty contours for Omega_M and the dark energy equation of state.

preprint2011arXiv

Massive black hole binaries: dynamical evolution and observational signatures

The study of the dynamical evolution of massive black hole pairs in mergers is crucial in the context of a hierarchical galaxy formation scenario. The timescales for the formation and the coalescence of black hole binaries are still poorly constrained, resulting in large uncertainties in the expected rate of massive black hole binaries detectable in the electromagnetic and gravitational wave spectra. Here we review the current theoretical understanding of the black hole pairing in galaxy mergers, with a particular attention to recent developments and open issues. We conclude with a review of the expected observational signatures of massive binaries, and of the candidates discussed in literature to date.

preprint2011arXiv

Multimessenger astronomy with pulsar timing and X-ray observations of massive black hole binaries

We demonstrate that very massive (>10^8\msun), cosmologically nearby (z<1) black hole binaries (MBHBs), which are primary targets for ongoing and upcoming pulsar timing arrays (PTAs), are particularly appealing multimessenger carriers. According to current models for massive black hole formation and evolution, the planned Square Kilometer Array (SKA) will collect gravitational wave signals from thousands of such massive systems, being able to individually resolve and locate in the sky several of them (maybe up to a hundred). By employing a standard model for the evolution of MBHBs in circumbinary discs, with the aid of dedicated numerical simulations, we characterize the gas-binary interplay, identifying possible electromagnetic signatures of the PTA sources. We concentrate our investigation on two particularly promising scenarios in the high energy domain, namely, the detection of X-ray periodic variability and of double broad Kαiron lines. Up to several hundreds of periodic X-ray sources with a flux >10^-13 erg s^-1 cm^-2 will be in the reach of upcoming X-ray observatories. Double relativistic Kαlines may be observable in a handful of low redshift (z<0.3) sources by proposed deep X-ray probes, such as Athena. (Abridged)

preprint2010arXiv

Gravitational waves and pulsar timing: stochastic background, individual sources and parameter estimation

Massive black holes are key ingredients of the assembly and evolution of cosmic structures. Pulsar Timing Arrays (PTAs) currently provide the only means to observe gravitational radiation from massive black hole binary systems with masses >10^7 solar masses. The whole cosmic population produces a signal consisting of two components: (i) a stochastic background resulting from the incoherent superposition of radiation from the all the sources, and (ii) a handful of individually resolvable signals that raise above the background level and are produced by sources sufficiently close and/or massive. Considering a wide range of massive black hole binary assembly scenarios, we investigate both the level and shape of the background and the statistics of resolvable sources. We predict a characteristic background amplitude in the interval h_c(f = 10^-8 Hz)~5*10^-16 - 5*10^-15, within the detection range of the complete Parkes PTA. We also quantify the capability of PTAs of measuring the parameters of individual sources, focusing on monochromatic signals produced by binaries in circular orbit. We investigate how the results depend on the number and distribution of pulsars in the array, by computing the variance-covariance matrix of the parameter measurements. For plausible Square Kilometre Array (SKA) observations (100 pulsars uniformly distributed in the sky), and assuming a coherent signal-to-noise ratio of 10, the sky position of massive black hole binaries can be located within a ~40deg^2 error box, opening promising prospects for detecting a putative electromagnetic counterpart to the gravitational wave emission. The planned SKA, can plausibly observe these unique systems, although the number of detections is likely to be small. (Abridged)

preprint2010arXiv

Self consistent model for the evolution of eccentric massive black hole binaries in stellar environments: implications for gravitational wave observations

We construct evolutionary tracks for massive black hole binaries (MBHBs) embedded in a surrounding distribution of stars. The dynamics of the binary is evolved by taking into account the erosion of the central stellar cusp bound to the massive black holes, the scattering of unbound stars feeding the binary loss cone, and the emission of gravitational waves (GWs). Stellar dynamics is treated in a hybrid fashion by coupling the results of numerical 3-body scattering experiments of bound and unbound stars to an analytical framework for the evolution of the stellar density distribution and for the efficiency of the binary loss cone refilling. Our main focus is on the behaviour of the binary eccentricity, in the attempt of addressing its importance in the merger process and its possible impact for GW detection with the planned Laser Interferometer Space Antenna ({\it LISA}), and ongoing and forthcoming pulsar timing array (PTA) campaigns. We produce a family of evolutionary tracks extensively sampling the relevant parameters of the system which are the binary mass, mass ratio and initial eccentricity, the slope of the stellar density distribution, its normalization and the efficiency of loss cone refilling. We find that, in general, stellar dynamics causes a dramatic increase of the MBHB eccentricity, especially for initially already mildly eccentric and/or unequal mass binaries. When applied to standard MBHB population models, our results predict eccentricities in the ranges $10^{-3}-0.2$ and $0.03-0.3$ for sources detectable by {\it LISA} and PTA respectively. Such figures may have a significant impact on the signal modelling, on source detection, and on the development of parameter estimation algorithms.

preprint2010arXiv

Strong lensing of gravitational waves as seen by LISA

We discuss strong gravitational lensing of gravitational waves from merging of massive black hole binaries in the context of the LISA mission. Detection of multiple events would provide invaluable information on competing theories of gravity, evolution and formation of structures and, with complementary observations, constraints on H_0 and other cosmological parameters. Most of the optical depth for lensing is provided by intervening massive galactic halos, for which wave optics effects are negligible. Probabilities to observe multiple events are sizable for a broad range of formation histories. For the most optimistic models, up to 4 multiple events with a signal to noise ratio >= 8 are expected in a 5-year mission. Chances are significant even for conservative models with either light (<= 60%) or heavy (<= 40%) seeds. Due to lensing amplification, some intrinsically too faint signals are brought over threshold (<= 2 per year).

preprint2009arXiv

Enhanced tidal disruption rates from massive black hole binaries

"Hard" massive black hole (MBH) binaries embedded in steep stellar cusps can shrink via three-body slingshot interactions. We show that this process will inevitably be accompanied by a burst of stellar tidal disruptions, at a rate that can be several orders of magnitude larger than that appropriate for a single MBH. Our numerical scattering experiments reveal that: 1) a significant fraction of stars initially bound to the primary hole are scattered into its tidal disruption loss cone by gravitational interactions with the secondary hole, an enhancement effect that is more pronounced for very unequal-mass binaries; 2) about 25% (40%) of all strongly interacting stars are tidally disrupted by a MBH binary of mass ratio q=1/81 (q=1/243) and eccentricity 0.1; and 3) two mechanisms dominate the fueling of the tidal disruption loss cone, a Kozai non-resonant interaction that causes the secular evolution of the stellar angular momentum in the field of the binary, and the effect of close encounters with the secondary hole that change the stellar orbital parameters in a chaotic way. For a hard MBH binary of 10^7 solar masses and mass ratio 0.01, embedded in an isothermal stellar cusp of velocity dispersion sigma*=100 km/s, the tidal disruption rate can be as large as 1/yr. This is 4 orders of magnitude higher than estimated for a single MBH fed by two-body relaxation. When applied to the case of a putative intermediate-mass black hole inspiraling onto Sgr A*, our results predict tidal disruption rates ~0.05-0.1/yr.

preprint2009arXiv

The international pulsar timing array project: using pulsars as a gravitational wave detector

The International Pulsar Timing Array project combines observations of pulsars from both Northern and Southern hemisphere observatories with the main aim of detecting ultra-low frequency (~10^-9 to 10^-8 Hz) gravitational waves. Here we introduce the project, review the methods used to search for gravitational waves emitted from coalescing supermassive binary black-hole systems in the centres of merging galaxies and discuss the status of the project.

preprint2007arXiv

Interaction of massive black hole binaries with their stellar environment: II. Loss-cone depletion and binary orbital decay

We study the long-term evolution of massive black hole binaries (MBHBs) at the centers of galaxies using detailed scattering experiments to solve the full three-body problem. Ambient stars drawn from a isotropic Maxwellian distribution unbound to the binary are ejected by the gravitational slingshot. We construct a minimal, hybrid model for the depletion of the loss cone and the orbital decay of the binary, and show that secondary slingshots - stars returning on small impact parameter orbits to have a second super-elastic scattering with the MBHB - may considerably help the shrinking of the pair in the case of large binary mass ratios. In the absence of loss-cone refilling by two-body relaxation or other processes, the mass ejected before the stalling of a MBHB is half the binary reduced mass. About 50% of the ejected stars are expelled ejected in a "burst" lasting ~1E4 yrs M_6^1/4, where M_6 is the binary mass in units of 1E6 Msun. The loss cone is completely emptied in a few bulge crossing timescales, 1E7 yrs M_6^1/4. Even in the absence of two-body relaxation or gas dynamical processes, unequal mass and/or eccentric binaries with M_6 >0.1 can shrink to the gravitational wave emission regime in less than a Hubble time, and are therefore "safe" targets for the planned Laser Interferometer Space Antenna (LISA).