In Einstein's general relativity, gravity is mediated by a massless metric field. The extension of general relativity to consistently include a mass for the graviton has profound implications for gravitation and cosmology. Salient features of various massive gravity theories can be captured by Galileon models, the simplest of which is the cubic Galileon. The presence of the Galileon field leads to additional gravitational radiation in binary pulsars where the Vainshtein mechanism is less suppressed than its fifth-force counterpart, which deserves a detailed confrontation with observations. We prudently choose fourteen well-timed binary pulsars, and from their intrinsic orbital decay rates we put a new bound on the graviton mass, $m_g \lesssim 2 \times 10^{-28}\,{\rm eV}/c^2$ at the 95% confidence level, assuming a flat prior on $\ln m_g$. It is equivalent to a bound on the graviton Compton wavelength $\lambda_g \gtrsim 7 \times 10^{21}\,{\rm m}$. Furthermore, we extensively simulate times of arrival for pulsars in orbit around stellar-mass black holes and the supermassive black hole at the Galactic center, and investigate their prospects in probing the cubic Galileon theory in the near future.

General relativity is a fully conservative theory, but there exist other possible metric theories of gravity. We consider non-conservative ones with a parameterized post-Newtonian (PPN) parameter, $\zeta_2$. A non-zero $\zeta_2$ induces a self-acceleration for the center of mass of an eccentric binary pulsar system, which contributes to the second time derivative of the pulsar spin frequency, $\ddot{\nu}$. In our work, using the method in Will (1992), we provide an improved analysis with four well-timed, carefully-chosen binary pulsars. In addition, we extend Will's method and derive $\zeta_2$'s effect on the third time derivative of the spin frequency, $\dddot{\nu}$. For PSR B1913+16, the constraint from $\dddot{\nu}$ is even tighter than that from $\ddot{\nu}$. We combine multiple pulsars with Bayesian inference, and obtain an upper limit, $\left|\zeta_{2}\right|<1.3\times10^{-5}$ at 95% confidence level, assuming a flat prior in $\log_{10} \left| \zeta_{2}\right|$. It improves the existing bound by a factor of three. Moreover, we propose an analytical timing formalism for $\zeta_2$. Our simulated times of arrival with simplified assumptions show binary pulsars' capability in limiting $\zeta_{2}$, and useful clues are extracted for real data analysis in future. In particular, we discover that for PSRs B1913+16 and J0737$-$3039A, $\dddot{\nu}$ can yield more constraining limits than $\ddot{\nu}$.

We describe system verification tests and early science results from the pulsar processor (PTUSE) developed for the newly-commissioned 64-dish SARAO MeerKAT radio telescope in South Africa. MeerKAT is a high-gain (~2.8 K/Jy) low-system temperature (~18 K at 20cm) radio array that currently operates from 580-1670 MHz and can produce tied-array beams suitable for pulsar observations. This paper presents results from the MeerTime Large Survey Project and commissioning tests with PTUSE. Highlights include observations of the double pulsar J0737-3039A, pulse profiles from 34 millisecond pulsars from a single 2.5h observation of the Globular cluster Terzan 5, the rotation measure of Ter5O, a 420-sigma giant pulse from the Large Magellanic Cloud pulsar PSR J0540-6919, and nulling identified in the slow pulsar PSR J0633-2015. One of the key design specifications for MeerKAT was absolute timing errors of less than 5 ns using their novel precise time system. Our timing of two bright millisecond pulsars confirm that MeerKAT delivers exceptional timing. PSR J2241-5236 exhibits a jitter limit of <4 ns per hour whilst timing of PSR J1909-3744 over almost 11 months yields an rms residual of 66 ns with only 4 min integrations. Our results confirm that the MeerKAT is an exceptional pulsar telescope. The array can be split into four separate sub-arrays to time over 1000 pulsars per day and the future deployment of S-band (1750-3500 MHz) receivers will further enhance its capabilities.

The gravitational strong equivalence principle (SEP) is a cornerstone of the general theory of relativity (GR). The extreme difference in binding energy between neutron stars and white dwarfs allows for precision tests of the SEP via the technique of pulsar timing. To date, the best limit on the validity of SEP under strong-field conditions was obtained with a unique pulsar in a triple stellar system, PSR J0337+1715. We report here on an improvement of this test using an independent data set acquired over 6 years with the Nan\c cay radio telescope (NRT). The improvements arise from a uniformly sampled data set, a theoretical analysis, and a treatment that fixes some short-comings in the previously published results, leading to better precision and reliability of the test. In contrast to the previously published test, we use a different long-term timing data set, developed a new timing model and an independent numerical integration of the motion of the system, and determined the masses and orbital parameters with a different methodology that treats the parameter $\Delta$, describing a possible strong-field SEP violation, identically to all other parameters. We obtain a violation parameter $\Delta = (+0.5 \pm 1.8) \times 10^{-6}$ at 95\% confidence level, which is compatible with and improves upon the previous study by 30\%. This result is statistics-limited and avoids limitation by systematics as previously encountered. We find evidence for red noise in the pulsar spin frequency, which is responsible for up to 10\% of the reported uncertainty. We use the improved limit on SEP violation to place constraints on a class of well-studied scalar-tensor theories, in particular we find $\omega_{\rm BD} > 140\,000$ for the Brans-Dicke parameter. The conservative limits presented here fully take into account current uncertainties in the equation for state of neutron-star matter.

Realistic synthetic observations of theoretical source models are essential for our understanding of real observational data. In using synthetic data, one can verify the extent to which source parameters can be recovered and evaluate how various data corruption effects can be calibrated. These studies are important when proposing observations of new sources, in the characterization of the capabilities of new or upgraded instruments, and when verifying model-based theoretical predictions in a comparison with observational data. We present the SYnthetic Measurement creator for long Baseline Arrays (SYMBA), a novel synthetic data generation pipeline for Very Long Baseline Interferometry (VLBI) observations. SYMBA takes into account several realistic atmospheric, instrumental, and calibration effects. We used SYMBA to create synthetic observations for the Event Horizon Telescope (EHT), a mm VLBI array, which has recently captured the first image of a black hole shadow. After testing SYMBA with simple source and corruption models, we study the importance of including all corruption and calibration effects. Based on two example general relativistic magnetohydrodynamics (GRMHD) model images of M87, we performed case studies to assess the attainable image quality with the current and future EHT array for different weather conditions. The results show that the effects of atmospheric and instrumental corruptions on the measured visibilities are significant. Despite these effects, we demonstrate how the overall structure of the input models can be recovered robustly after performing calibration steps. With the planned addition of new stations to the EHT array, images could be reconstructed with higher angular resolution and dynamic range. In our case study, these improvements allowed for a distinction between a thermal and a non-thermal GRMHD model based on salient features in reconstructed images.

Radio pulsars in short-period eccentric binary orbits can be used to study both gravitational dynamics and binary evolution. The binary system containing PSR J1141$-$6545 includes a massive white dwarf (WD) companion that formed before the gravitationally bound young radio pulsar. We observe a temporal evolution of the orbital inclination of this pulsar that we infer is caused by a combination of a Newtonian quadrupole moment and Lense-Thirring precession of the orbit resulting from rapid rotation of the WD. Lense-Thirring precession, an effect of relativistic frame-dragging, is a prediction of general relativity. This detection is consistent with the evolutionary scenario in which the WD accreted matter from the pulsar progenitor, spinning up the WD to a period $< 200$ seconds.

The grand challenges of contemporary fundamental physics---dark matter, dark energy, vacuum energy, inflation and early universe cosmology, singularities and the hierarchy problem---all involve gravity as a key component. And of all gravitational phenomena, black holes stand out in their elegant simplicity, while harbouring some of the most remarkable predictions of General Relativity: event horizons, singularities and ergoregions. The hitherto invisible landscape of the gravitational Universe is being unveiled before our eyes: the historical direct detection of gravitational waves by the LIGO-Virgo collaboration marks the dawn of a new era of scientific exploration. Gravitational-wave astronomy will allow us to test models of black hole formation, growth and evolution, as well as models of gravitational-wave generation and propagation. It will provide evidence for event horizons and ergoregions, test the theory of General Relativity itself, and may reveal the existence of new fundamental fields. The synthesis of these results has the potential to radically reshape our understanding of the cosmos and of the laws of Nature. The purpose of this work is to present a concise, yet comprehensive overview of the state of the art in the relevant fields of research, summarize important open problems, and lay out a roadmap for future progress.

Pulsars in the Galactic Center (GC) are important probes of General Relativity, star formation, stellar dynamics, stellar evolution, and the interstellar medium. Despite years of searching, only a handful of pulsars in the central 0.5 deg are known. The high-frequency sensitivity of ngVLA will open a new window for discovery and characterization of pulsars in the GC. A pulsar in orbit around the GC black hole, Sgr A*, will provide an unprecedented probe of black hole physics and General Relativity.

The violation of the weak equivalence principle (EP) in the gravitational field of the Earth, described by the E\"otv\"os parameter $\eta_\oplus$, was recently constrained to the level $\left|\eta_\oplus\right| \lesssim 10^{-14}$ by the MICROSCOPE space mission. The E\"otv\"os parameter $\eta_{\rm DM}$, pertaining to the differential couplings of dark matter (DM) and ordinary matter, was only tested to the level $\left| \eta_{\rm DM} \right| \lesssim 10^{-5}$ by the E\"ot-Wash group and lunar laser ranging. This test is limited by the EP-violating driving force in the Solar neighborhood that is determined by the Galactic distribution of DM. Here we propose a novel celestial experiment using the orbital dynamics from radio timing of binary pulsars, and obtain a competing limit on $\eta_{\rm DM}$ from a neutron star (NS) -- white dwarf (WD) system, PSR J1713+0747. The result benefits from the large material difference between the NS and the WD, and the large gravitational binding energy of the NS. If we can discover a binary pulsar within $\sim 10$ parsecs of the Galactic center, where the driving force is much larger in the expected DM spike, precision timing will improve the test of the universality of free fall towards DM and constrain various proposed couplings of DM to the Standard Model by several orders of magnitude. Such a test probes the hypothesis that gravity is the only long-range interaction between DM and ordinary matter.

The MeerKAT telescope represents an outstanding opportunity for radio pulsar timing science with its unique combination of a large collecting area and aperture efficiency (effective area $\sim$7500 m$^2$), system temperature ($T<20$K), high slew speeds (1-2 deg/s), large bandwidths (770 MHz at 20cm wavelengths), southern hemisphere location (latitude $\sim -30^\circ$) and ability to form up to four sub-arrays. The MeerTime project is a five-year program on the MeerKAT array by an international consortium that will regularly time over 1000 radio pulsars to perform tests of relativistic gravity, search for the gravitational wave signature induced by supermassive black hole binaries in the timing residuals of millisecond pulsars, explore the interiors of neutron stars through a pulsar glitch monitoring programme, explore the origin and evolution of binary pulsars, monitor the swarms of pulsars that inhabit globular clusters and monitor radio magnetars. In addition to these primary programmes, over 1000 pulsars will have their arrival times monitored and the data made immediately public. The MeerTime pulsar backend comprises two server-class machines each of which possess four Graphics Processing Units. Up to four pulsars can be coherently dedispersed simultaneously up to dispersion measures of over 1000 pc cm$^{-3}$. All data will be provided in psrfits format. The MeerTime backend will be capable of producing coherently dedispersed filterbank data for timing multiple pulsars in the cores of globular clusters that is useful for pulsar searches of tied array beams. All MeerTime data will ultimately be made available for public use, and any published results will include the arrival times and profiles used in the results.

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.

Pulsars in relativistic binary systems have emerged as fantastic natural laboratories for testing theories of gravity, the most prominent example being the double pulsar, PSR J0737$-$3039. The HTRU-South Low Latitude pulsar survey represents one of the most sensitive blind pulsar surveys taken of the southern Galactic plane to date, and its primary aim has been the discovery of new relativistic binary pulsars. Here we present our binary pulsar searching strategy and report on the survey's flagship discovery, PSR J1757$-$1854. A 21.5-ms pulsar in a relativistic binary with an orbital period of 4.4 hours and an eccentricity of 0.61, this double neutron star (DNS) system is the most accelerated pulsar binary known, and probes a relativistic parameter space not yet explored by previous pulsar binaries.

We report the discovery of PSR J1757$-$1854, a 21.5-ms pulsar in a highly-eccentric, 4.4-h orbit around a neutron star (NS) companion. PSR J1757$-$1854 exhibits some of the most extreme relativistic parameters of any known pulsar, including the strongest relativistic effects due to gravitational-wave (GW) damping, with a merger time of 76 Myr. Following a 1.6-yr timing campaign, we have measured five post-Keplerian (PK) parameters, yielding the two component masses ($m_\text{p}=1.3384(9)\,\text{M}_\odot$ and $m_\text{c}=1.3946(9)\,\text{M}_\odot$) plus three tests of general relativity (GR), which the theory passes. The larger mass of the NS companion provides important clues regarding the binary formation of PSR J1757$-$1854. With simulations suggesting 3-$\sigma$ measurements of both the contribution of Lense-Thirring precession to the rate of change of the semi-major axis and the relativistic deformation of the orbit within $\sim7-9$ years, PSR J1757$-$1854 stands out as a unique laboratory for new tests of gravitational theories.

In this work, we report the discovery and characterization of PSR J1411+2551, a new binary pulsar discovered in the Arecibo 327 MHz Drift Pulsar Survey. Our timing observations of the radio pulsar in the system span a period of about 2.5 years. This timing campaign allowed a precise measurement of its spin period (62.4 ms) and its derivative (9.6 $\pm$ 0.7) $\times 10^{-20}\, \rm s\, s^{-1}$; from these, we derive a characteristic age of $\sim 10\,$Gyr and a surface magnetic field strength of 2.5 $\times 10^{9}$ G. These numbers indicate that this pulsar was mildly recycled by accretion of matter from the progenitor of the companion star. The system has an eccentric ($e\, = \, 0.17$) 2.61 day orbit. This eccentricity allows a highly significant measurement of the rate of advance of periastron, $\dot{\omega} = 0.07686 \pm 0.00046 ^{\circ}~{\rm yr}^{-1}$. Assuming general relativity accurately models the orbital motion, this implies a total system mass M = $2.538 \pm 0.022 M_{\odot}$. The minimum companion mass is $0.92\, M_{\odot}$ and the maximum pulsar mass is $1.62\, M_{\odot}$. The large companion mass and the orbital eccentricity suggest that PSR J1411+2551 is a double neutron star system; the lightest known to date including the DNS merger GW 170817. Furthermore, the relatively low orbital eccentricity and small proper motion limits suggest that the second supernova had a relatively small associated kick; this and the low system mass suggest that it was an ultra-stripped supernova.

Multi-decade observing campaigns of the globular clusters 47 Tucanae and M15 have led to an outstanding number of discoveries. Here, we report on the latest results of the long-term observations of the pulsars in these two clusters. For most of the pulsars in 47 Tucanae we have measured, among other things, their higher-order spin period derivatives, which have in turn provided stringent constraints on the physical parameters of the cluster, such as its distance and gravitational potential. For M15, we have studied the relativistic spin precession effect in PSR B2127+11C. We have used full-Stokes observations to model the precession effect, and to constrain the system geometry. We find that the visible beam of the pulsar is swiftly moving away from our line of sight and may very soon become undetectable. On the other hand, we expect to see the opposite emission beam sometime between 2041 and 2053.

Double neutron star (DNS) systems represent extreme physical objects and the endpoint of an exotic journey of stellar evolution and binary interactions. Large numbers of DNS systems and their mergers are anticipated to be discovered using the Square-Kilometre-Array searching for radio pulsars and high-frequency gravitational wave detectors (LIGO/VIRGO), respectively. Here we discuss all key properties of DNS systems, as well as selection effects, and combine the latest observational data with new theoretical progress on various physical processes with the aim of advancing our knowledge on their formation. We examine key interactions of their progenitor systems and evaluate their accretion history during the high-mass X-ray binary stage, the common envelope phase and the subsequent Case BB mass transfer, and argue that the first-formed NSs have accreted at most $\sim 0.02\;M_{\odot}$. We investigate DNS masses, spins and velocities, and in particular correlations between spin period, orbital period and eccentricity. Numerous Monte Carlo simulations of the second supernova (SN) events are performed to extrapolate pre-SN stellar properties and probe the explosions. All known close-orbit DNS systems are consistent with ultra-stripped exploding stars. Although their resulting NS kicks are often small, we demonstrate a large spread in kick magnitudes which may, in general, depend on the past interaction history of the exploding star and thus correlate with the NS mass. We analyze and discuss NS kick directions based on our SN simulations. Finally, we discuss the terminal evolution of close-orbit DNS systems until they merge and possibly produce a short $\gamma$-ray burst.

Pulsar timing and gravitational-wave (GW) detectors are superb laboratories to study gravity theories in the strong-field regime. Here we combine those tools to test the mono-scalar-tensor theory of Damour and Esposito-Far{\`e}se (DEF), which predicts nonperturbative scalarization phenomena for neutron stars (NSs). First, applying Markov-chain Monte Carlo techniques, we use the absence of dipolar radiation in the pulsar-timing observations of five binary systems composed of a NS and a white dwarf, and eleven equations of state (EOSs) for NSs, to derive the most stringent constraints on the two free parameters of the DEF scalar-tensor theory. Since the binary-pulsar bounds depend on the NS mass and the EOS, we find that current pulsar-timing observations leave scalarization windows, i.e., regions of parameter space where scalarization can still be prominent. Then, we investigate if these scalarization windows could be closed and if pulsar-timing constraints could be improved by laser-interferometer GW detectors, when spontaneous (or dynamical) scalarization sets in during the early (or late) stages of a binary NS (BNS) evolution. For the early inspiral of a BNS carrying constant scalar charge, we employ a Fisher matrix analysis to show that Advanced LIGO can improve pulsar-timing constraints for some EOSs, and next-generation detectors, such as the Cosmic Explorer and Einstein Telescope, will be able to improve those bounds for all eleven EOSs. Using the late inspiral of a BNS, we estimate that for some of the EOSs under consideration the onset of dynamical scalarization can happen early enough to improve the constraints on the DEF parameters obtained by combining the five binary pulsars. Thus, in the near future the complementarity of pulsar timing and direct observations of GWs on the ground will be extremely valuable in probing gravity theories in the strong-field regime.

We report on the results of a 4-year timing campaign of PSR~J2222$-0137$, a 2.44-day binary pulsar with a massive white dwarf (WD) companion, with the Nan\c{c}ay, Effelsberg and Lovell radio telescopes. Using the Shapiro delay for this system, we find a pulsar mass $m_{p}=1.76,\pm\,0.06,M_\odot$ and a WD mass $m_{c}\,=\,1.293\,\pm\,0.025\, M_\odot$. We also measure the rate of advance of periastron for this system, which is marginally consistent with the GR prediction for these masses. The short lifetime of the massive WD progenitor star led to a rapid X-ray binary phase with little ($< \, 10^{-2} \, M_\odot$) mass accretion onto the neutron star (NS); hence, the current pulsar mass is, within uncertainties, its birth mass; the largest measured to date. We discuss the discrepancy with previous mass measurements for this system; we conclude that the measurements presented here are likely to be more accurate. Finally, we highlight the usefulness of this system for testing alternative theories of gravity by tightly constraining the presence of dipolar radiation. This is of particular importance for certain aspects of strong-field gravity, like spontaneous scalarization, since the mass of PSR~J2222$-0137$ puts that system into a poorly tested parameter range.

Einstein's General Theory of Relativity (GR) successfully describes gravity. The most fundamental predictions of GR are black holes (BHs), but in spite of many convincing BH candidates in the Universe, there is no conclusive experimental proof of their existence using astronomical observations in the electromagnetic spectrum. Are BHs real astrophysical objects? Does GR hold in its most extreme limit or are alternatives needed? The prime target to address these fundamental questions is in the center of our own Galaxy, which hosts the closest and best-constrained supermassive BH candidate in the Universe, Sagittarius A* (Sgr A*). Three different types of experiments hold the promise to test GR in a strong-field regime using observations of Sgr A* with new-generation instruments. The first experiment aims to image the relativistic plasma emission which surrounds the event horizon and forms a "shadow" cast against the background, whose predicted size (~50 microarcseconds) can now be resolved by upcoming VLBI experiments at mm-waves such as the Event Horizon Telescope (EHT). The second experiment aims to monitor stars orbiting Sgr A* with the upcoming near-infrared interferometer GRAVITY at the Very Large Telescope (VLT). The third experiment aims to time a radio pulsar in tight orbit about Sgr A* using radio telescopes (including the Atacama Large Millimeter Array or ALMA). The BlackHoleCam project exploits the synergy between these three different techniques and aims to measure the main BH parameters with sufficient precision to provide fundamental tests of GR and probe the spacetime around a BH in any metric theory of gravity. Here, we review our current knowledge of the physical properties of Sgr A* as well as the current status of such experimental efforts towards imaging the event horizon, measuring stellar orbits, and timing pulsars around Sgr A*.

Symmetries play important roles in modern theories of physical laws. In this paper, we review several experimental tests of important symmetries associated with the gravitational interaction, including the universality of free fall for self-gravitating bodies, time-shift symmetry in the gravitational constant, local position invariance and local Lorentz invariance of gravity, and spacetime translational symmetries. Recent experimental explorations for post-Newtonian gravity are discussed, of which, those from pulsar astronomy are highlighted. All of these tests, of very different aspects of gravity theories, at very different length scales, favor to very high precision the predictions of the strong equivalence principle (SEP) and, in particular, general relativity which embodies SEP completely. As the founding principles of gravity, these symmetries are motivated to be promoted to even stricter tests in future.

Timing results for the black-widow pulsar J2051-0827 are presented, using a 21-year dataset from four European Pulsar Timing Array telescopes and the Parkes radio telescope. This dataset, which is the longest published to date for a black-widow system, allows for an improved analysis that addresses previously unknown biases. While secular variations, as identified in previous analyses, are recovered, short-term variations are detected for the first time. Concurrently, a significant decrease of approx. 2.5x10-3 cm-3 pc in the dispersion measure associated with PSR J2051-0827 is measured for the first time and improvements are also made to estimates of the proper motion. Finally, PSR J2051-0827 is shown to have entered a relatively stable state suggesting the possibility of its eventual inclusion in pulsar timing arrays.

The Double Pulsar system PSR J0737-3039A/B has proven to be an excellent laboratory for high precision tests of general relativity. With additional years of timing measurements and new telescopes like the Square Kilometre Array (SKA), the precision of these tests will increase and new effects like the Lense-Thirring precession of the orbit will become measurable. Here, we discuss the prospects of measuring the Lense-Thirring effect and thereby constraining the equations of state at supra-nuclear densities in neutron stars using the Double Pulsar.

One century after its formulation, Einstein's general relativity has made remarkable predictions and turned out to be compatible with all experimental tests. Most of these tests probe the theory in the weak-field regime, and there are theoretical and experimental reasons to believe that general relativity should be modified when gravitational fields are strong and spacetime curvature is large. The best astrophysical laboratories to probe strong-field gravity are black holes and neutron stars, whether isolated or in binary systems. We review the motivations to consider extensions of general relativity. We present a (necessarily incomplete) catalog of modified theories of gravity for which strong-field predictions have been computed and contrasted to Einstein's theory, and we summarize our current understanding of the structure and dynamics of compact objects in these theories. We discuss current bounds on modified gravity from binary pulsar and cosmological observations, and we highlight the potential of future gravitational wave measurements to inform us on the behavior of gravity in the strong-field regime.

The black hole in the center of the Milky Way, Sgr A*, has the largest mass-to-distance ratio among all known black holes in the Universe. This property makes Sgr A* the optimal target for testing the gravitational no-hair theorem. In the near future, major developments in instrumentation will provide the tools for high-precision studies of its spacetime via observations of relativistic effects in stellar orbits, in the timing of pulsars, and in horizon-scale images of its accretion flow. We explore here the prospect of measuring the properties of the black-hole spacetime using all these three types of observations. We show that the correlated uncertainties in the measurements of the black-hole spin and quadrupole moment using the orbits of stars and pulsars are nearly orthogonal to those obtained from measuring the shape and size of the shadow the black hole casts on the surrounding emission. Combining these three types of observations will, therefore, allow us to assess and quantify systematic biases and uncertainties in each measurement and lead to a highly accurate, quantitative test of the gravitational no-hair theorem.

The discovery and timing of radio pulsars within the Galactic centre is a fundamental aspect of the SKA Science Case, responding to the topic of "Strong Field Tests of Gravity with Pulsars and Black Holes" (Kramer et al. 2004; Cordes et al. 2004). Pulsars have in many ways proven to be excellent tools for testing the General theory of Relativity and alternative gravity theories (see Wex (2014) for a recent review). Timing a pulsar in orbit around a companion, provides a unique way of probing the relativistic dynamics and spacetime of such a system. The strictest tests of gravity, in strong field conditions, are expected to come from a pulsar orbiting a black hole. In this sense, a pulsar in a close orbit ($P_{\rm orb}$ < 1 yr) around our nearest supermassive black hole candidate, Sagittarius A* - at a distance of ~8.3 kpc in the Galactic centre (Gillessen et al. 2009a) - would be the ideal tool. Given the size of the orbit and the relativistic effects associated with it, even a slowly spinning pulsar would allow the black hole spacetime to be explored in great detail (Liu et al. 2012). For example, measurement of the frame dragging caused by the rotation of the supermassive black hole, would allow a test of the "cosmic censorship conjecture." The "no-hair theorem" can be tested by measuring the quadrupole moment of the black hole. These are two of the prime examples for the fundamental studies of gravity one could do with a pulsar around Sagittarius A*. As will be shown here, SKA1-MID and ultimately the SKA will provide the opportunity to begin to find and time the pulsars in this extreme environment.

The anticipated discovery of a pulsar in orbit with a black hole is expected to provide a unique laboratory for black hole physics and gravity. In this context, the next generation of radio telescopes, like the Five-hundred-metre Aperture Spherical radio Telescope (FAST) and the Square Kilometre Array (SKA), with their unprecedented sensitivity, will play a key role. In this paper, we investigate the capability of future radio telescopes to probe the spacetime of a black hole and test gravity theories, by timing a pulsar orbiting a stellar-mass-black-hole (SBH). Based on mock data simulations, we show that a few years of timing observations of a sufficiently compact pulsar-SBH (PSR-SBH) system with future radio telescopes would allow precise measurements of the black hole mass and spin. A measurement precision of one per cent can be expected for the spin. Measuring the quadrupole moment of the black hole, needed to test GR's no-hair theorem, requires extreme system configurations with compact orbits and a large SBH mass. Additionally, we show that a PSR-SBH system can lead to greatly improved constraints on alternative gravity theories even if they predict black holes (practically) identical to GR's. This is demonstrated for a specific class of scalar-tensor theories. Finally, we investigate the requirements for searching for PSR-SBH systems. It is shown that the high sensitivity of the next generation of radio telescopes is key for discovering compact PSR-SBH systems, as it will allow for sufficiently short survey integration times.

Very long baseline interferometry at millimetre/submillimetre wavelengths (mmVLBI) offers the highest achievable spatial resolution at any wavelength in astronomy. The anticipated inclusion of ALMA as a phased array into a global VLBI network will bring unprecedented sensitivity and a transformational leap in capabilities for mmVLBI. Building on years of pioneering efforts in the US and Europe the ongoing ALMA Phasing Project (APP), a US-led international collaboration with MPIfR-led European contributions, is expected to deliver a beamformer and VLBI capability to ALMA by the end of 2014 (APP: Fish et al. 2013, arXiv:1309.3519). This report focuses on the future use of mmVLBI by the international users community from a European viewpoint. Firstly, it highlights the intense science interest in Europe in future mmVLBI observations as compiled from the responses to a general call to the European community for future research projects. A wide range of research is presented that includes, amongst others: - Imaging the event horizon of the black hole at the centre of the Galaxy - Testing the theory of General Relativity an/or searching for alternative theories - Studying the origin of AGN jets and jet formation - Cosmological evolution of galaxies and BHs, AGN feedback - Masers in the Milky Way (in stars and star-forming regions) - Extragalactic emission lines and astro-chemistry - Redshifted absorption lines in distant galaxies and study of the ISM and circumnuclear gas - Pulsars, neutron stars, X-ray binaries - Testing cosmology - Testing fundamental physical constants

We report on the discovery of four millisecond pulsars (MSPs) in the High Time Resolution Universe (HTRU) pulsar survey being conducted at the Parkes 64-m radio telescope. All four MSPs are in binary systems and are likely to have white dwarf companions. In addition, we present updated timing solutions for 12 previously published HTRU MSPs, revealing new observational parameters such as five proper motion measurements and significant temporal dispersion measure variations in PSR J1017-7156. We discuss the case of PSR J1801-3210, which shows no significant period derivative after four years of timing data. Our best-fit solution shows a period derivative of the order of $10^{-23}$, an extremely small number compared to that of a typical MSP. However, it is likely that the pulsar lies beyond the Galactic Centre, and an unremarkable intrinsic period derivative is reduced to close to zero by the Galactic potential acceleration. Furthermore, we highlight the potential to employ PSR J1801-3210 in the strong equivalence principle test due to its wide and circular orbit. In a broader comparison with the known MSP population, we suggest a correlation between higher mass functions and the presence of eclipses in `very low-mass binary pulsars', implying that eclipses are observed in systems with high orbital inclinations. We also suggest that the distribution of the total mass of binary systems is inversely-related to the Galactic height distribution. Finally, we report on the first detection of PSRs J1543-5149 and J1811-2404 as gamma-ray pulsars.

The centre of our Milky Way harbours the closest candidate for a supermassive black hole. The source is thought to be powered by radiatively inefficient accretion of gas from its environment. This form of accretion is a standard mode of energy supply for most galactic nuclei. X-ray measurements have already resolved a tenuous hot gas component from which it can be fed. However, the magnetization of the gas, a crucial parameter determining the structure of the accretion flow, remains unknown. Strong magnetic fields can influence the dynamics of the accretion, remove angular momentum from the infalling gas, expel matter through relativistic jets and lead to the observed synchrotron emission. Here we report multi-frequency measurements with several radio telescopes of a newly discovered pulsar close to the Galactic Centre and show that its unusually large Faraday rotation indicates a dynamically relevant magnetic field near the black hole. If this field is accreted down to the event horizon it provides enough magnetic flux to explain the observed emission from the black hole, from radio to X-rays.

Gravitational preferred frame effects are generally predicted by alternative theories that exhibit an isotropic violation of local Lorentz invariance of gravity. They are described by three parameters in the parametrized post-Newtonian formalism. One of their strong-field generalizations, $\hat \alpha_2$, induces a precession of a pulsar's spin around its movement direction with respect to the preferred frame. We constrain $\hat \alpha_2$ by using the non-detection of such a precession using the characteristics of the pulse profile. In our analysis we use a large number of observations from the 100-m Effelsberg radio telescope, which cover a time span of approximately 15 years. By combining data from two solitary millisecond pulsars, PSRs B1937+21 and J1744-1134, we get a limit of $|\hat \alpha_2| < 1.6 \times 10^{-9}$ at 95% confidence level, which is more than two orders of magnitude better than its best weak-field counterpart from the Solar system.

Within the parameterized post-Newtonian (PPN) formalism, there could be an anisotropy of local gravity induced by an external matter distribution, even for a fully conservative metric theory of gravity. It reflects the breakdown of the local position invariance of gravity and, within the PPN formalism, is characterized by the Whitehead parameter $\xi$. We present three different kinds of observation, from the Solar system and radio pulsars, to constrain it. The most stringent limit comes from recent results on the extremely stable pulse profiles of solitary millisecond pulsars, that gives $|\hat \xi| < 3.9 \times 10^{-9}$ (95% CL), where the hat denotes the strong-field generalization of $\xi$. This limit is six orders of magnitude more constraining than the current best limit from superconducting gravimeter experiments. It can be converted into an upper limit of $\sim 4 \times 10^{-16}$ on the spatial anisotropy of the gravitational constant.

Many physically motivated extensions to general relativity (GR) predict significant deviations in the properties of spacetime surrounding massive neutron stars. We report the measurement of a 2.01 +/- 0.04 solar mass pulsar in a 2.46-hr orbit with a 0.172 +/- 0.003 solar mass white dwarf. The high pulsar mass and the compact orbit make this system a sensitive laboratory of a previously untested strong-field gravity regime. Thus far, the observed orbital decay agrees with GR, supporting its validity even for the extreme conditions present in the system. The resulting constraints on deviations support the use of GR-based templates for ground-based gravitational wave detectors. Additionally, the system strengthens recent constraints on the properties of dense matter and provides insight to binary stellar astrophysics and pulsar recycling.

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.

New tests are proposed to constrain possible deviations from local Lorentz invariance and local position invariance in the gravity sector. By using precise timing results of two binary pulsars, i.e., PSRs J1012+5307 and J1738+0333, we are able to constrain (strong-field) parametrized post-Newtonian parameters $\hat{\alpha}_1$, $\hat{\alpha}_2$, $\hat{\xi}$ to high precision, among which, $|\hat{\xi}| < 3.1\times10^{-4}$ (95% C.L.) is reported here for the first time.

The discovery of a pulsar (PSR) in orbit around a black hole (BH) is expected to provide a superb new probe of relativistic gravity and BH properties. Apart from a precise mass measurement for the BH, one could expect a clean verification of the dragging of space-time caused by the BH spin. In order to measure the quadrupole moment of the BH for testing the no-hair theorem of general relativity (GR), one has to hope for a sufficiently massive BH. In this respect, a PSR orbiting the super-massive BH in the center of our Galaxy would be the ultimate laboratory for gravity tests with PSRs. But even for gravity theories that predict the same properties for BHs as GR, a PSR-BH system would constitute an excellent test system, due to the high grade of asymmetry in the strong field properties of these two components. Here we highlight some of the potential gravity tests that one could expect from different PSR-BH systems, utilizing present and future radio telescopes, like FAST and SKA.

Radio pulsars in relativistic binary systems are unique tools to study the curved space-time around massive compact objects. The discovery of a pulsar closely orbiting the super-massive black hole at the centre of our Galaxy, Sgr A*, would provide a superb test-bed for gravitational physics. To date, the absence of any radio pulsar discoveries within a few arc minutes of Sgr A* has been explained by one principal factor: extreme scattering of radio waves caused by inhomogeneities in the ionized component of the interstellar medium in the central 100 pc around Sgr A*. Scattering, which causes temporal broadening of pulses, can only be mitigated by observing at higher frequencies. Here we describe recent searches of the Galactic centre region performed at a frequency of 18.95 GHz with the Effelsberg radio telescope.

Preferred frame effects (PFEs) are predicted by a number of alternative gravity theories which include vector or additional tensor fields, besides the canonical metric tensor. In the framework of parametrized post-Newtonian (PPN) formalism, we investigate PFEs in the orbital dynamics of binary pulsars, characterized by the two strong-field PPN parameters, \alpha_1 and \alpha_2. In the limit of a small orbital eccentricity, \alpha_1 and \alpha_2 contributions decouple. By utilizing recent radio timing results and optical observations of PSRs J1012+5307 and J1738+0333, we obtained the best limits of \alpha_1 and \alpha_2 in the strong-field regime. The constraint on \alpha_1 also surpasses its counterpart in the weak-field regime.

In the post-Newtonian parametrization of semi-conservative gravity theories, local Lorentz invariance (LLI) violation is characterized by two parameters, alpha_1 and alpha_2. In binary pulsars the isotropic violation of LLI in the gravitational sector leads to characteristic preferred frame effects (PFEs) in the orbital dynamics, if the barycenter of the binary is moving relative to the preferred frame with a velocity w. For small-eccentricity binaries, the effects induced by alpha_1 and alpha_2 decouple, and can therefore be tested independently. We use recent timing results of two compact pulsar-white dwarf binaries with known 3D velocity, PSRs J1012+5307 and J1738+0333, to constrain PFEs for strongly self-gravitating bodies. We derive a limit |alpha_2| < 1.8e-4 (95% CL), which is the most constraining limit for strongly self-gravitating systems up to now. Concerning alpha_1, we propose a new, robust method to constrain this parameter. Our most conservative result, alpha_1 = -0.4^{+3.7}_{-3.1} e-5 (95% CL) from PSR J1738+0333, constitutes a significant improvement compared to current most stringent limits obtained both in Solar system and binary pulsar tests. We also derive corresponding limits for alpha_1 and alpha_2 for a preferred frame that is at rest with respect to our Galaxy, and preferred frames that locally co-move with the rotation of our Galaxy. (Abridged)

PSR J1738+0333 is one of the four millisecond pulsars known to be orbited by a white dwarf companion bright enough for optical spectroscopy. Of these, it has the shortest orbital period, making it especially interesting for a range of astrophysical and gravity related questions. We present a spectroscopic and photometric study of the white dwarf companion and infer its radial velocity curve, effective temperature, surface gravity and luminosity. We find that the white dwarf has properties consistent with those of low-mass white dwarfs with thick hydrogen envelopes, and use the corresponding mass-radius relation to infer its mass; M_WD = 0.181 +/- +0.007/-0.005 solar masses. Combined with the mass ratio q=8.1 +/- 0.2 inferred from the radial velocities and the precise pulsar timing ephemeris, the neutron star mass is constrained to M_PSR = 1.47 +/- +0.07/-0.06 solar masses. Contrary to expectations, the latter is only slightly above the Chandrasekhar limit. We find that, even if the birth mass of the neutron star was only 1.20 solar masses, more than 60% of the matter that left the surface of the white dwarf progenitor escaped the system. The accurate determination of the component masses transforms this system in a laboratory for fundamental physics by constraining the orbital decay predicted by general relativity. Currently, the agreement is within 1 sigma of the observed decay. Further radio timing observations will allow precise tests of white dwarf models, assuming the validity of general relativity.

In this paper, we review tests of the strong equivalence principle (SEP) derived from pulsar-white dwarf data. The extreme difference in binding energy between both components and the precise measurement of the orbital motion provided by pulsar timing allow the only current precision SEP tests for strongly self-gravitating bodies. We start by highlighting why such tests are conceptually important. We then review previous work where limits on SEP violation are obtained with an ensemble of wide binary systems with small eccentricity orbits. Then we propose a new SEP violation test based on the measurement of the variation of the orbital eccentricity de/dt. This new method has the following advantages: a) unlike previous methods it is not based on probabilistic considerations, b) it can make a direct detection of SEP violation, c) the measurement of de/dt is not contaminated by any known external effects, which implies that this SEP test is only restricted by the measurement precision of de/dt. In the final part of the review, we conceptually compare the SEP test with the test for dipolar radiation damping, a phenomenon closely related to SEP violation, and speculate on future prospects by new types of tests in globular clusters and future triple systems.

(abridged) We report the results of a 10-year timing campaign on PSR J1738+0333, a 5.85-ms pulsar in a low-eccentricity 8.5-hour orbit with a low-mass white dwarf companion (...) The measurements of proper motion and parallax allow for a precise subtraction of the kinematic contribution to the observed orbital decay; this results in a significant measurement of the intrinsic orbital decay: (-25.9 +/- 3.2) \times 10^{-15} s/s. This is consistent with the orbital decay from the emission of gravitational waves predicted by general relativity, (-27.7 +1.5/-1.9) \times 10^{-15} s/s (...). This agreement introduces a tight upper limit on dipolar gravitational wave emission, a prediction of most alternative theories of gravity for asymmetric binary systems such as this. We use this limit to derive the most stringent constraints ever on a wide class of gravity theories, where gravity involves a scalar field contribution. When considering general scalar-tensor theories of gravity, our new bounds are more stringent than the best current solar-system limits over most of the parameter space, and constrain the matter-scalar coupling constant {\alpha}_0^2 to be below the 10^{-5} level. For the special case of the Jordan-Fierz-Brans-Dicke, we obtain the one-sigma bound {\alpha}_0^2 < 2 \times 10^{-5}, which is within a factor two of the Cassini limit. We also use our limit on dipolar gravitational wave emission to constrain a wide class of theories of gravity which are based on a generalization of Bekenstein's Tensor-Vector-Scalar gravity (TeVeS), a relativistic formulation of Modified Newtonian Dynamics (MOND).

The discovery of radio pulsars in compact orbits around Sgr A* would allow an unprecedented and detailed investigation of the spacetime of the supermassive black hole. This paper shows that pulsar timing, including that of a single pulsar, has the potential to provide novel tests of general relativity, in particular its cosmic censorship conjecture and no-hair theorem for rotating black holes. These experiments can be performed by timing observations with 100 micro-second precision, achievable with the Square Kilometre Array for a normal pulsar at frequency above 15 GHz. Based on the standard pulsar timing technique, we develop a method that allows the determination of the mass, spin, and quadrupole moment of Sgr A*, and provides a consistent covariance analysis of the measurement errors. Furthermore, we test this method in detailed mock data simulations. It seems likely that only for orbital periods below ~0.3 yr is there the possibility of having negligible external perturbations. For such orbits we expect a ~10^-3 test of the frame dragging and a ~10^-2 test of the no-hair theorem within 5 years, if Sgr A* is spinning rapidly. Our method is also capable of identifying perturbations caused by distributed mass around Sgr A*, thus providing high confidence in these gravity tests. Our analysis is not affected by uncertainties in our knowledge of the distance to the Galactic center, R0. A combination of pulsar timing with the astrometric results of stellar orbits would greatly improve the measurement precision of R0.

Abbreviated: We investigate the potential of detecting the gravitational wave from individual binary black hole systems using pulsar timing arrays (PTAs) and calculate the accuracy for determining the GW properties. This is done in a consistent analysis, which at the same time accounts for the measurement of the pulsar distances via the timing parallax. We find that, at low redshift, a PTA is able to detect the nano-Hertz GW from super massive black hole binary systems with masses of $\sim10^8 - 10^{10}\,M_{\sun}$ less than $\sim10^5$\,years before the final merger, and those with less than $\sim10^3 - 10^4$ years before merger may allow us to detect the evolution of binaries. We derive an analytical expression to describe the accuracy of a pulsar distance measurement via timing parallax. We consider five years of bi-weekly observations at a precision of 15\,ns for close-by ($\sim 0.5 - 1$\,kpc) pulsars. Timing twenty pulsars would allow us to detect a GW source with an amplitude larger than $5\times 10^{-17}$. We calculate the corresponding GW and binary orbital parameters and their measurement precision. The accuracy of measuring the binary orbital inclination angle, the sky position, and the GW frequency are calculated as functions of the GW amplitude. We note that the "pulsar term", which is commonly regarded as noise, is essential for obtaining an accurate measurement for the GW source location. We also show that utilizing the information encoded in the GW signal passing the Earth also increases the accuracy of pulsar distance measurements. If the gravitational wave is strong enough, one can achieve sub-parsec distance measurements for nearby pulsars with distance less than $\sim 0.5 - 1$\,kpc.

We have conducted radio timing observations of the eclipsing millisecond binary pulsar J2051-0827 with the European Pulsar Timing Array network of telescopes and the Parkes radio telescope, spanning over 13 years. The increased data span allows significant measurements of the orbital eccentricity, e = (6.2 {\pm} 1.3) {\times} 10^{-5} and composite proper motion, {\mu}_t = 7.3 {\pm} 0.4 mas/yr. Our timing observations have revealed secular variations of the projected semi-major axis of the pulsar orbit which are much more extreme than those previously published; and of the orbital period of the system. Investigation of the physical mechanisms producing such variations confirm that the variations of the semi-major axis are most probably caused by classical spin-orbit coupling in the binary system, while the variations in orbital period are most likely caused by tidal dissipation leading to changes in the gravitational quadrupole moment of the companion.

Parallax measurements of pulsars allow for accurate measurements of the interstellar electron density and contribute to accurate tests of general relativity using binary systems. The Square Kilometre Array (SKA) will be an ideal instrument for measuring the parallax of pulsars, because it has a very high sensitivity, as well as baselines extending up to several thousands of kilometres. We performed simulations to estimate the number of pulsars for which the parallax can be measured with the SKA and the distance to which a parallax can be measured. We compare two different methods. The first method measures the parallax directly by utilising the long baselines of the SKA to form high angular resolution images. The second method uses the arrival times of the radio signals of pulsars to fit a transformation between time coordinates in the terrestrial frame and the comoving pulsar frame directly yielding the parallax. We find that with the first method a parallax with an accuracy of 20% or less can be measured up to a maximum distance of 13 kpc, which would include 9,000 pulsars. By timing pulsars with the most stable arrival times for the radio emission, parallaxes can be measured for about 3,600 millisecond pulsars up to a distance of 9 kpc with an accuracy of 20%.

(abridged) PSR J1903+0327, a millisecond pulsar in an eccentric (e = 0.44) 95-day orbit with a (~ 1Msun) companion poses a challenge to our understanding of stellar evolution in binary and multiple-star systems. Here we describe optical and radio observations which rule out most of the scenarios proposed to explain formation of this system. Radio timing measurements of three post-Keplerian effects yield the most precise measurement of the mass of a millisecond pulsar to date: 1.667 +/- 0.021 solar masses (99.7% confidence limit) (...). Optical spectroscopy of a proposed main sequence counterpart show that its orbital motion mirrors the pulsar's 95-day orbit; being therefore its binary companion (...) The optical detection also provides a measurement of the systemic radial velocity of the binary; this and the proper motion measured from pulsar timing allow the determination of the systemic 3-D velocity in the Galaxy. We find that the system is always within 270 pc of the plane of the Galaxy, but always more than 3 kpc away from the Galactic centre. Thus an exchange interaction in a dense stellar environment (like a globular cluster or the Galactic centre) is not likely to be the origin of this system. We suggest that after the supernova that formed it, the neutron star was in a tight orbit with a main-sequence star, the present companion was a tertiary farther out. The neutron star then accreted matter from its evolving inner MS companion, forming a millisecond pulsar. The former donor star then disappears, either due to a chaotic 3-body interaction with the outer star (caused by the expansion of the inner orbit that necessarily results from mass transfer), or in the case of a very compact inner system, due to ablation/accretion by the newly formed millisecond pulsar.

Pulsars with compact companions in close eccentric orbits are unique laboratories for testing general relativity and alternative theories of gravity. Moreover, they are excellent targets for future gravitational wave experiments like LISA and they are also highly important for understanding the equation of state of super-dense matter and the evolution of massive binaries. Here we report on optical observations of the 1.02 M_sun companion to the pulsar PSR J1141-6545. We detect an optical counterpart with apparent magnitudes V=25.08(11) and R=24.38(14), consistent with the timing position of the pulsar. We demonstrate that our results are in agreement with a white dwarf companion. However the latter is redder than expected and the inferred values are not consistent with the theoretical cooling tracks, preventing us from deriving the exact age. Our results confirm the importance of the PSR J1141-6545 system for gravitational experiments.

Massive gravitons are features of some alternatives to general relativity. This has motivated experiments and observations that, so far, have been consistent with the zero mass graviton of general relativity, but further tests will be valuable. A basis for new tests may be the high sensitivity gravitational wave experiments that are now being performed, and the higher sensitivity experiments that are being planned. In these experiments it should be feasible to detect low levels of dispersion due to nonzero graviton mass. One of the most promising techniques for such a detection may be the pulsar timing program that is sensitive to nano-Hertz gravitational waves. Here we present some details of such a detection scheme. The pulsar timing response to a gravitational wave background with the massive graviton is calculated, and the algorithm to detect the massive graviton is presented. We conclude that, with 90% probability, massles gravitons can be distinguished from gravitons heavier than $3\times 10^{-22}$ eV (Compton wave length $\lambda_{\rm g}=4.1 \times 10^{12}$ km), if biweekly observation of 60 pulsars are performed for 5 years with pulsar RMS timing accuracy of 100 ns. If 60 pulsars are observed for 10 years with the same accuracy, the detectable graviton mass is reduced to $5\times 10^{-23}$ eV ($\lambda_{\rm g}=2.5 \times 10^{13}$ km); for 5-year observations of 100 or 300 pulsars, the sensitivity is respectively $2.5\times 10^{-22}$ ($\lambda_{\rm g}=5.0\times 10^{12}$ km) and $10^{-22}$ eV ($\lambda_{\rm g}=1.2\times 10^{13}$ km). Finally, a 10-year observation of 300 pulsars with 100 ns timing accuracy would probe graviton masses down to $3\times 10^{-23}$ eV ($\lambda_{\rm g}=4.1\times 10^{13}$ km).

(abridged) In this paper, we express the relativistic propagational delay of light in the space-time of a binary system (commonly known as the "Shapiro delay") as a sum of harmonics of the orbital period of the system. We do this first for near-circular orbits as a natural expansion of an existing orbital model for low-eccentricity binary systems. The amplitudes of the 3rd and higher harmonics can be described by two new post-Keplerian (PK) parameters proportional to the amplitudes of the third and fourth harmonics (h_3, h_4). For high orbital inclinations we use a PK parameter proportional to the ratio of amplitudes of successive harmonics (sigma) instead of h_4. The new PK parameters are much less correlated with each other than r and s and provide a superior description of the constraints introduced by the Shapiro delay on the orbital inclination and the masses of the components of the binary (...). We extend the h_3,sigma parameterisation to eccentric binaries with high orbital inclinations. For some such binaries we can measure extra PK parameters and test general relativity using the Shapiro delay parameters. In this case we can use the measurement of h_3 as a test of general relativity. We show that this new test is not only more stringent than the r test, but it is even more stringent than the previous s test. Until now this new parametric test could only be derived statistically from an analysis of a probabilistic chi2 map.

We present preliminary results of 7 years of Arecibo timing of the pulsar-white dwarf binary PSR J1738+0333. We can measure the proper motion, parallax with excellent precision and have detected the orbital decay. Furthermore, the companion has been detected at optical wavelengths and a mass ratio of 8.1 +/- 0.3 has been measured from the orbital variation of its Doppler shift. Once the companion mass is determined from the optical measurements, this system will provide strong limits for the radiation of dipolar gravitational waves. Assuming that general relativity holds, the fast-improving measurement of the orbital decay, combined with the measurement of the mass ratio, will provide an independent and precise measurement of the component masses.

Presently the double pulsar allows for the measurement of six post-Keplerian parameters. In addition, its double-line nature gives access to the projected semi-major axes of both orbits. We use this wealth of information to pose some very general restrictions on a wide class of conservative and semi-conservative theories of gravity.

We report on the high precision timing analysis of the pulsar-white dwarf binary PSR J1012+5307. Using 15 years of multi-telescope data from the European Pulsar Timing Array (EPTA) network, a significant measurement of the variation of the orbital period is obtained. Using this ideal strong-field gravity laboratory we derive theory independent limits for both the dipole radiation and the variation of the gravitational constant.

We present results from the high precision timing analysis of the pulsar-white dwarf (WD) binary PSR J1012+5307 using 15 years of multi-telescope data. Observations were performed regularly by the European Pulsar Timing Array (EPTA) network, consisting of Effelsberg, Jodrell Bank, Westerbork and Nan\c{c}ay. All the timing parameters have been improved from the previously published values, most by an order of magnitude. In addition, a parallax measurement of $\pi = 1.2(3)$ mas is obtained for the first time for PSR J1012+5307, being consistent with the optical estimation from the WD companion. Combining improved 3D velocity information and models for the Galactic potential the complete evolutionary Galactic path of the system is obtained. A new intrinsic eccentricity upper limit of $e<8.4\times 10^{-7}$ is acquired, one of the smallest calculated for a binary system and a measurement of the variation of the projected semi-major axis also constrains the system's orbital orientation for the first time. It is shown that PSR J1012+5307 is an ideal laboratory for testing alternative theories of gravity. The measurement of the change of the orbital period of the system of $\dot{P}_{b} = 5(1)\times 10^{-14}$ is used to set an upper limit on the dipole gravitational wave emission that is valid for a wide class of alternative theories of gravity. Moreover, it is shown that in combination with other binary pulsars PSR J1012+5307 is an ideal system to provide self-consistent, generic limits, based only on millisecond pulsar data, for the dipole radiation and the variation of the gravitational constant $\dot{G}$.

In 2004, McLaughlin et al. discovered a phenomenon in the radio emission of PSR J0737-3039B (B) that resembles drifting sub-pulses. The repeat rate of the sub-pulses is equal to the spin frequency of PSR J0737-3039A (A); this led to the suggestion that they are caused by incidence upon B's magnetosphere of electromagnetic radiation from A. Here we describe a geometrical model which predicts the delay of B's sub-pulses relative to A's radio pulses. We show that measuring these delays is equivalent to tracking A's rotation from the point of view of an hypothetical observer located near B. This has three main astrophysical applications: (a) to determine the sense of rotation of A relative to its orbital plane; (b) to estimate where in B's magnetosphere the radio sub-pulses are modulated and (c) to provide an independent estimate of the mass ratio of A and B. The latter might improve existing tests of gravitational theories using this system.

In this paper we develop a consistent, phenomenological methodology to measure preferred-frame effects (PFEs) in binary pulsars that exhibit a high rate of periastron advance. We show that in these systems the existence of a preferred frame for gravity leads to an observable characteristic `signature' in the timing data, which uniquely identifies this effect. We expand the standard Damour-Deruelle timing formula to incorporate this `signature' and show how this new PFE timing model can be used to either measure or constrain the parameters related to a violation of the local Lorentz invariance of gravity in the strong internal fields of neutron stars. In particular, we demonstrate that in the presence of PFEs we expect a set of the new timing parameters to have a unique relationship that can be measured and tested incontrovertibly. This new methodology is applied to the Double Pulsar, which turns out to be the ideal test system for this kind of experiments.The currently available dataset allows us only to study the impact of PFEs on the orbital precession rate, d omega/dt, providing limits that are, at the moment, clearly less stringent than existing limits on PFE strong-field parameters. However, simulations show that the constraints improve fast in the coming years, allowing us to study all new PFE timing parameters and to check for the unique relationship between them. Finally, we show how a combination of several suitable systems in a "PFE antenna array", expected to be availabe for instance with the Square-Kilometre-Array (SKA), provides full sensitivity to possible violations of local Lorentz invariance in strong gravitational fields in all directions of the sky. This PFE antenna array may eventually allow us to determine the direction of a preferred frame should it exist.

We present results and applications of high precision timing measurements of the millisecond pulsar J1012+5307. Combining our radio observations with results based on optical observations, we derive complete 3-D velocity information for this system. Correcting for Doppler effects, we derive the intrinsic spin parameters of this pulsar and a characteristic age of 8.6 +/- 1.9 Gyr. Our upper limit for the orbital eccentricity of only 8 * 10^-7 (68% C.L.) is the smallest ever measured for a binary system. We demonstrate that this makes the pulsar an ideal laboratory to test certain aspects of alternative theories of gravitation. Our precise measurements suggest deviations from a simple pulsar spin-down timing model, which are consistent with timing noise and the extrapolation of the known behaviour of slowly rotating pulsars.

We report here on multifrequency radio observations of the pulsed emission from PSR B1259-63 around the time of the closest approach to its B2e companion star. There was a general increase in the pulsar's dispersion measure and scatter broadening, and a decrease in the flux density towards periastron although changes in these parameters were seen on timescales as short as minutes. The pulsed emission disappeared 16 days prior to periastron and remained undetectable until 16 days after periastron. The observations are used to determine the parameters of the wind from the Be star. We show that a simple model, in which the wind density varies with radius as r^-2, provides a good fit to the data. The wind is clumpy with size scales <= 10^10 cm, densities of 10^6 cm^-3 and a velocity of 2000 km/s at a distance of 20 - 50 stellar radii. We find a correlation between dispersion measure variations and the pulse scattering times, suggesting that the same electrons are responsible for both effects.

We have used the Effelsberg 100-m radio telescope to conduct a 1400-MHz ($\lambda$ 21-cm) search for young and rapidly rotating radio pulsars along a 2 deg$^2$ strip of the northern Galactic plane defined by $28^{\circ} \leq l \leq 30^{\circ}$ and $|b|\leq0.5^{\circ}$. This region lies close to the Scutum spiral arm which is already known to contain a number of radio and X-ray pulsars. The search was nominally sensitive to pulsars with 1400-MHz flux densities above 0.3 mJy; this represents a threefold improvement in sensitivity over all previous searches of this region of the Galaxy. Four new long-period pulsars were discovered as a result of this survey. All three previously known pulsars in this region were also detected. The four new pulsars are relatively young ($<$ 0.5 Myr), weak ($<1$ mJy) sources with dispersion measures in the range 180--910 cm$^{-3}$ pc. None of the newly-discovered pulsars are associated with catalogued supernova remnants.

A change of the component separation in the profiles of the binary pulsar PSR B1913+16 has been observed for the first time (Kramer 1998) as expected by geodetic precession. In this work we extend the previous work by accounting for recent data from the Effelsberg 100-m telescope and Arecibo Observatory and testing model predictions. We demonstrate how the new information will provide additional information on the solutions of the system geometry.

We use recent information on geodetic precession of the binary pulsar B1913+16 along with measurements of its orbital parameters and proper motion to derive new constraints on the immediate progenitor of this double neutron star system. As part of our analysis we model the motion of the binary in the Galaxy after the second supernova explosion, and we derive constraints on the unknown radial velocity. We also obtain limits on the magnitude and direction of the kick velocity imparted to the pulsar companion during the second supernova explosion. We consider the complete set of possible cases, depending on the kinematic age of the system and the 180 degree ambiguity in the pulsar spin orientation. Most interestingly, we find that the natal kick must have been directed almost perpendicular to the spin axis of the neutron star progenitor, independent of the specific pre-supernova configuration. Such a tight constraint on the kick direction has important implications for the physical mechanism responsible for the kick.

The discovery of frame-dragging effects in binary pulsar timing experiments requires a compact companion with sufficiently large spin. A pulsar orbiting a fast rotating black hole could provide an appropriate test system. In this paper we address questions concerning the identification of a black hole companion in such a system, the measurability of the frame dragging caused by the rotation of the black hole, and the measurability of the quadrupole moment, which would prove the presence of a Kerr black hole.

For radio pulsars in orbit with a compact companion, pulsar timing observations have proved to be a powerful tool for identifying the physical nature of the companion. Unfortunately, perhaps the most intriguing system where such a tool could be used, a pulsar in orbit with a black hole, has yet to be discovered. In this paper we give a detailed investigation of what one can learn about the black hole companion via timing observations of the pulsar. We present an analytic calculation for the propagation delay caused by the frame-dragging effect and show that it has the same functional behavior as the modulation of the observed rotational phase of the pulsar caused by the deflection of the radio signals in the gravitational field of the companion (bending delay). Thus, contrary to statements of other authors, the frame-dragging delay is unlikely to be separately measurable in pulsar binaries where the companion is a stellar mass black hole. We demonstrate, however, that the precession of the binary orbit caused by the relativistic spin-orbit coupling can lead to observable effects which can be used to set a lower limit to the black-hole spin, or possibly allow the determination of its magnitude and orientation. We give parameter estimates for two possible systems, a 10 solar mass black hole in orbit either with a young (~0.1 s) pulsar or with a millisecond pulsar. Finally, we discuss the measurability of the quadrupole moment of the rotating black hole companion which would test the presence of a Kerr black hole. As an interesting side result of our calculations, we can give a further argument why the companion of PSR J0045-7319 cannot be a Kerr black hole.

The pulsar PSR B1259-63 is in a highly eccentric 3.4-yr orbit with the Be star SS 2883. Timing observations of this pulsar, made over a 7-yr period using the Parkes 64-m radio telescope, cover two periastron passages, in 1990 August and 1994 January. The timing data cannot be fitted by the normal pulsar and Keplerian binary parameters. A timing solution including a (non-precessing) Keplerian orbit and timing noise (represented as a polynomial of fifth order in time) provide a satisfactory fit to the data. However, because the Be star probably has a significant quadrupole moment, we prefer to interpret the data by a combination of timing noise, dominated by a cubic phase term, and $\dot\omega$ and $\dot x$ terms. We show that the $\dot\omega$ and $\dot x$ are likely to be a result of a precessing orbit caused by the quadrupole moment of the tilted companion star. We further rule out a number of possible physical effects which could contribute to the timing data of PSR B1259-63 on a measurable level.

In binary radio pulsars with a main-sequence star companion, the spin-induced quadrupole moment of the companion gives rise to a precession of the binary orbit. As a first approximation one can model the secular evolution caused by this classical spin-orbit coupling by linear-in-time changes of the longitude of periastron and the projected semi-major axis of the pulsar orbit. This simple representation of the precession of the orbit neglects two important aspects of the orbital dynamics of a binary pulsar with an oblate companion. First, the quasiperiodic effects along the orbit, due to the anisotropic $1/r^3$ nature of the quadrupole potential. Secondly, the long-term secular evolution of the binary orbit which leads to an evolution of the longitude of periastron and the projected semi-major axis which is non-linear in time. In this paper a simple timing formula for binary radio pulsars with a main-sequence star companion is presented which models the short-term secular and most of the short-term periodic effects caused by the classical spin-orbit coupling. I also give extensions of the timing formula which account for long-term secular changes in the binary pulsar motion. It is shown that the short-term periodic effects are important for the timing observations of the binary pulsar PSR B1259--63. The long-term secular effects are likely to become important in the next few years of timing observations of the binary pulsar PSR J0045--7319. They could help to restrict or even determine the moments of inertia of the companion star and thus probe its internal structure. Finally, I reinvestigate the spin-orbit precession of the binary pulsar PSR J0045--7319 since the analysis given in the literature is based on an incorrect expression for the precession of the longitude of periastron.

The Strong Equivalence Principle (SEP) demands, besides the validity of the Einstein Equivalence Principle, that all self-gravitating bodies feel the same acceleration in an external gravitational field. It has been found that metric theories of gravity other that than general relativity typically predict a violation of the SEP. In case of the Earth-Moon system (weak field system) this violation is called the Nordtvedt effect. It has been shown by Damour and Sch\"afer, that small-eccentricity long-orbital-period binary pulsars with a white dwarf companion provide excellent conditions to test the SEP in strong field regimes. Based on newly discovered binary pulsars this paper investigates a possible violation of the SEP in strong field regimes. New limits with an improved confidence level are presented. The results of this paper lead to constrains on the combination $\epsilon/2-\zeta$ of the only two (post)$^2$-Newtonian parameters $\epsilon$ and $\zeta$ that arise from the (post)$^2$-Newtonian approximation of the tensor-multi-scalar theory of Damour and Esposito-Far{\`e}se.

We investigate the possibility of identifying massive objects (lenses) in the Galactic Center region (GC) by means of pulsar timing. The well known intensity change due to microlensing is found to be less important. For typical stellar masses, the flux magnification can be significant only if the lens passes very close to the pulsar--Earth axis. We show that in the case of a pulsar the time varying travel-time delay, which is observable because of the pulsating nature of pulsar radiation, is a much more powerful tool to investigate mass distributions in the GC. We find that a travel-time delay is measurable even for rather large distances between the lens and the pulsar--Earth axis. The time varying travel-time delay can be used to determine the mass of the lens and the ratio of the transverse velocity to the minimum impact parameter. Based on Monte Carlo simulations of synthetic data sets we give the expected accuracies in the parameter determination for various masses. We argue that pulsars found behind the very center of our Galaxy would provide an excellent opportunity to test the mass distribution in the Galactic nucleus, and therefore to distinguish between a super massive black hole (SMBH) and a super-dense star cluster (SDSC) within the central 0.1 pc.