Search

Context: Timing of pulsar PSR J0337+1715 provides a unique opportunity to test the strong equivalence principle (SEP) with a strongly self-gravitating object. This is due to its unique situation in a triple stellar system with two white dwarfs. Aims: Our previous study suggested the presence of a strong low-frequency signal in the timing residuals. We set out to model it on a longer dataset in order to determine its nature and improve accuracy. Methods: Three models are considered: chromatic or achromatic red-noise, and a small planet in a hierarchical orbit with the triple stellar system. These models are implemented in our numerical timing model. We perform Bayesian inference of posterior distributions. Best fits are compared using information-theoretic criteria. Results: Chromatic red noise from dispersion-measure variations is ruled out. Achromatic red noise or a planet in keplerian orbit provide the best fits. If it is red noise then it appears exceptionally strong. Assuming the presence of a planet, we obtain a marginal detection of mutual interactions which allows us to constrain its mass to $\sim 0.5 M_{\rm Moon}$ as well as its inclination. The latter is intriguingly coincident with a Kozai resonance. We show that a longer observation span will ultimately lead to a clear signature of the planet model due to its mutual interactions with the triple system. We produce new limits on SEP violation: $|\Delta| < 1.5\cdot 10^{-6}$ or $|\Delta| < 2.3\cdot 10^{-6}$ at 95\% confidence level under the planet or red-noise hypothesis, respectively. This model dependence emphasises the need for additional data and model selection. As a by-product, we estimate a rather low supernova kick velocity of $\sim 110-125 \rm km/s$, strengthening the idea that it is a necessary condition for the formation of pulsar triple systems., Comment: Astronomy and Astrophysics - A\&A, In press

The discovery of the first pulsar in a binary star system, the Hulse--Taylor pulsar, 50 years ago opened up an entirely new field of experimental gravity. For the first time it was possible to investigate strong-field and radiative aspects of the gravitational interaction. Continued observations of the Hulse--Taylor pulsar eventually led, among other confirmations of the predictions of general relativity (GR), to the first evidence for the reality of gravitational waves. In the meantime, many more radio pulsars have been discovered that are suitable for testing GR and its alternatives. One particularly remarkable binary system is the Double Pulsar, which has far surpassed the Hulse--Taylor pulsar in several respects. In addition, binary pulsar-white dwarf systems have been shown to be particularly suitable for testing alternative gravitational theories, as they often predict strong dipolar gravitational radiation for such asymmetric systems. A rather unique pulsar laboratory is the pulsar in a hierarchical stellar triple, that led to by far the most precise confirmation of the strong-field version of the universality of free fall. Using radio pulsars, it could be shown that additional aspects of the Strong Equivalence Principle apply to the dynamics of strongly self-gravitating bodies, like the local position and local Lorentz invariance of the gravitational interaction. So far, GR has passed all pulsar tests with flying colours, while at the same time many alternative gravity theories have either been strongly constrained or even falsified. New telescopes, instrumentation, timing and search algorithms promise a significant improvement of the existing tests and the discovery of (qualitatively) new, more relativistic binary systems., Comment: Invited review article for Living Reviews in Relativity

The Event Horizon Telescope (EHT) observed in 2017 the supermassive black hole at the center of the Milky Way, Sagittarius A* (Sgr A*), at a frequency of 228.1 GHz ($\lambda$=1.3 mm). The fundamental physics tests that even a single pulsar orbiting Sgr A* would enable motivate searching for pulsars in EHT datasets. The high observing frequency means that pulsars - which typically exhibit steep emission spectra - are expected to be very faint. However, it also negates pulse scattering, an effect that could hinder pulsar detections in the Galactic Center. Additionally, magnetars or a secondary inverse Compton emission could be stronger at millimeter wavelengths than at lower frequencies. We present a search for pulsars close to Sgr A* using the data from the three most-sensitive stations in the EHT 2017 campaign: the Atacama Large Millimeter/submillimeter Array, the Large Millimeter Telescope and the IRAM 30 m Telescope. We apply three detection methods based on Fourier-domain analysis, the Fast-Folding-Algorithm and single pulse search targeting both pulsars and burst-like transient emission; using the simultaneity of the observations to confirm potential candidates. No new pulsars or significant bursts were found. Being the first pulsar search ever carried out at such high radio frequencies, we detail our analysis methods and give a detailed estimation of the sensitivity of the search. We conclude that the EHT 2017 observations are only sensitive to a small fraction ($\lesssim$2.2%) of the pulsars that may exist close to Sgr A*, motivating further searches for fainter pulsars in the region., Comment: 33 pages, 7 figures, 6 Tables. Accepted for publication in ApJ

Interpretation of resolved polarized images of black holes by the Event Horizon Telescope (EHT) requires predictions of the polarized emission observable by an Earth-based instrument for a particular model of the black hole accretion system. Such predictions are generated by general relativistic radiative transfer (GRRT) codes, which integrate the equations of polarized radiative transfer in curved spacetime. A selection of ray-tracing GRRT codes used within the EHT collaboration is evaluated for accuracy and consistency in producing a selection of test images, demonstrating that the various methods and implementations of radiative transfer calculations are highly consistent. When imaging an analytic accretion model, we find that all codes produce images similar within a pixel-wise normalized mean squared error (NMSE) of 0.012 in the worst case. When imaging a snapshot from a cell-based magnetohydrodynamic simulation, we find all test images to be similar within NMSEs of 0.02, 0.04, 0.04, and 0.12 in Stokes I, Q, U , and V respectively. We additionally find the values of several image metrics relevant to published EHT results to be in agreement to much better precision than measurement uncertainties., Comment: Accepted for publication in ApJ

Benefitting from the unequaled precision of the pulsar timing technique, binary pulsars are important testbeds of gravity theories, providing some of the tightest bounds on alternative theories of gravity. One class of well-motivated alternative gravity theories, the scalar-tensor gravity, predict large deviations from general relativity for neutron stars through a nonperturbative phenomenon known as spontaneous scalarization. This effect, which cannot be tested in the Solar System, can now be tightly constrained using the latest results from the timing of a set of 7 binary pulsars (PSRs J0348+0432, J0737$-$3039A, J1012+5307, J1738+0333, J1909$-$3744, J1913+1102, and J2222$-$0137), especially with the updated parameters of PSRs J0737$-$3039A, J1913+1102, and J2222$-$0137. Using new timing results, we constrain the neutron star's effective scalar coupling, which describes how strongly neutron stars couple to the scalar field, to a level of $\left|\alpha_{\rm A}\right| \lesssim 6 \times 10^{-3}$ in a Bayesian analysis. Our analysis is thorough, in the sense that our results apply to all neutron star masses and all reasonable equations of state of dense matters, in the full relevant parameter space. It excludes the possibility of spontaneous scalarization of neutron stars, at least within a class of scalar-tensor gravity theories., Comment: 14 pages, 5 figures; accepted by Class. Quantum Grav. Letters

Very-long-baseline interferometry (VLBI) observations of active galactic nuclei at millimeter wavelengths have the power to reveal the launching and initial collimation region of extragalactic radio jets, down to $10-100$ gravitational radii ($r_g=GM/c^2$) scales in nearby sources. Centaurus A is the closest radio-loud source to Earth. It bridges the gap in mass and accretion rate between the supermassive black holes (SMBHs) in Messier 87 and our galactic center. A large southern declination of $-43^{\circ}$ has however prevented VLBI imaging of Centaurus A below ${\lambda}1$cm thus far. Here, we show the millimeter VLBI image of the source, which we obtained with the Event Horizon Telescope at $228$GHz. Compared to previous observations, we image Centaurus A's jet at a tenfold higher frequency and sixteen times sharper resolution and thereby probe sub-lightday structures. We reveal a highly-collimated, asymmetrically edge-brightened jet as well as the fainter counterjet. We find that Centaurus A's source structure resembles the jet in Messier 87 on ${\sim}500r_g$ scales remarkably well. Furthermore, we identify the location of Centaurus A's SMBH with respect to its resolved jet core at ${\lambda}1.3$mm and conclude that the source's event horizon shadow should be visible at THz frequencies. This location further supports the universal scale invariance of black holes over a wide range of masses., Comment: 27 pages, 9 figures. This is a post-peer-review, pre-copyedit version of an article published in Nature Astronomy. The final authenticated version is available online at: https://doi.org/10.1038/s41550-021-01417-w

The black-hole images obtained with the Event Horizon Telescope (EHT) are expected to be variable at the dynamical timescale near their horizons. For the black hole at the center of the M87 galaxy, this timescale (5-61 days) is comparable to the 6-day extent of the 2017 EHT observations. Closure phases along baseline triangles are robust interferometric observables that are sensitive to the expected structural changes of the images but are free of station-based atmospheric and instrumental errors. We explored the day-to-day variability in closure phase measurements on all six linearly independent non-trivial baseline triangles that can be formed from the 2017 observations. We showed that three triangles exhibit very low day-to-day variability, with a dispersion of $\sim3-5^\circ$. The only triangles that exhibit substantially higher variability ($\sim90-180^\circ$) are the ones with baselines that cross visibility amplitude minima on the $u-v$ plane, as expected from theoretical modeling. We used two sets of General Relativistic magnetohydrodynamic simulations to explore the dependence of the predicted variability on various black-hole and accretion-flow parameters. We found that changing the magnetic field configuration, electron temperature model, or black-hole spin has a marginal effect on the model consistency with the observed level of variability. On the other hand, the most discriminating image characteristic of models is the fractional width of the bright ring of emission. Models that best reproduce the observed small level of variability are characterized by thin ring-like images with structures dominated by gravitational lensing effects and thus least affected by turbulence in the accreting plasmas., Comment: Accepted for Publication in ApJ

Our understanding of strong gravity near supermassive compact objects has recently improved thanks to the measurements made by the Event Horizon Telescope (EHT). We use here the M87* shadow size to infer constraints on the physical charges of a large variety of nonrotating or rotating black holes. For example, we show that the quality of the measurements is already sufficient to rule out that M87* is a highly charged dilaton black hole. Similarly, when considering black holes with two physical and independent charges, we are able to exclude considerable regions of the space of parameters for the doubly-charged dilaton and the Sen black holes., Comment: 15 pages, 4 figures, published in PRD on May 19

Synchrotron radiation from hot gas near a black hole results in a polarized image. The image polarization is determined by effects including the orientation of the magnetic field in the emitting region, relativistic motion of the gas, strong gravitational lensing by the black hole, and parallel transport in the curved spacetime. We explore these effects using a simple model of an axisymmetric, equatorial accretion disk around a Schwarzschild black hole. By using an approximate expression for the null geodesics derived by Beloborodov (2002) and conservation of the Walker-Penrose constant, we provide analytic estimates for the image polarization. We test this model using currently favored general relativistic magnetohydrodynamic simulations of M87*, using ring parameters given by the simulations. For a subset of these with modest Faraday effects, we show that the ring model broadly reproduces the polarimetric image morphology. Our model also predicts the polarization evolution for compact flaring regions, such as those observed from Sgr A* with GRAVITY. With suitably chosen parameters, our simple model can reproduce the EVPA pattern and relative polarized intensity in Event Horizon Telescope images of M87*. Under the physically motivated assumption that the magnetic field trails the fluid velocity, this comparison is consistent with the clockwise rotation inferred from total intensity images., Comment: 29 pages, 11 figures, published in ApJ on May 3

We report on the first pulsar and transient survey of the Galactic Center (GC) with the Atacama Large Millimeter/submillimeter Array (ALMA). The observations were conducted during the Global Millimeter VLBI Array campaign in 2017 and 2018. We carry out searches using timeseries of both total intensity and other polarization components in the form of Stokes parameters. We incorporate acceleration and its derivative in the pulsar search, and also search in segments of the entire observation to compensate for potential orbital motion of the pulsar. While no new pulsar is found, our observations yield the polarization profile of the GC magnetar PSR J1745-2900 at mm-wavelength for the first time, which turns out to be nearly 100 % linearly polarized. Additionally, we estimate the survey sensitivity placed by both system and red noise, and evaluate its capability of finding pulsars in orbital motion with either Sgr A* or a binary companion. We show that the survey is sensitive to only the most luminous pulsars in the known population, and future observations with ALMA in Band-1 will deliver significantly deeper survey sensitivity on the GC pulsar population., Comment: 16 pages, 11 figures, accepted for publication in ApJ

The 2017 Event Horizon Telescope (EHT) observations of the central source in M87 have led to the first measurement of the size of a black-hole shadow. This observation offers a new and clean gravitational test of the black-hole metric in the strong-field regime. We show analytically that spacetimes that deviate from the Kerr metric but satisfy weak-field tests can lead to large deviations in the predicted black-hole shadows that are inconsistent with even the current EHT measurements. We use numerical calculations of regular, parametric, non-Kerr metrics to identify the common characteristic among these different parametrizations that control the predicted shadow size. We show that the shadow-size measurements place significant constraints on deviation parameters that control the second post-Newtonian and higher orders of each metric and are, therefore, inaccessible to weak-field tests. The new constraints are complementary to those imposed by observations of gravitational waves from stellar-mass sources., Comment: Physical Review Letters

Radio pulsars provide some of the most important constraints for our understanding of matter at supranuclear densities. So far, these constraints are mostly given by precision mass measurements of neutron stars (NS). By combining single measurements of the two most massive pulsars, J0348$+$0432 and J0740$+$6620, the resulting lower limit of 1.98 $M_\odot$ (99% confidence) of the maximum NS mass, excludes a large number of equations of state (EOSs). Further EOS constraints, complementary to other methods, are likely to come from the measurement of the moment of inertia (MOI) of binary pulsars in relativistic orbits. The Double Pulsar, PSR J0737$-$3039A/B, is the most promising system for the first measurement of the MOI via pulsar timing. Reviewing this method, based in particular on the first MeerKAT observations of the Double Pulsar, we provide well-founded projections into the future by simulating timing observations with MeerKAT and the SKA. For the first time, we account for the spin-down mass loss in the analysis. Our results suggest that an MOI measurement with 11% accuracy (68% confidence) is possible by 2030. If by 2030 the EOS is sufficiently well known, however, we find that the Double Pulsar will allow for a 7% test of Lense-Thirring precession, or alternatively provide a $\sim3\sigma$-measurement of the next-to-leading order gravitational wave damping in GR. Finally, we demonstrate that potential new discoveries of double NS systems with orbital periods shorter than that of the Double Pulsar promise significant improvements in these measurements and the constraints on NS matter., Comment: 13 pages, 8 figures. Accepted by MNRAS

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., Comment: 14 pages, 10 figures; accepted by PRD

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}$., Comment: 14 pages, 5 figures, 4 tables; accepted by ApJ

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., Comment: Accepted in A&A

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., Comment: To be published in the ASP Monograph Series, "Science with a Next-Generation VLA", ed. E. J. Murphy (ASP, San Francisco, CA)

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., Comment: White Paper for the COST action "Gravitational Waves, Black Holes, and Fundamental Physics", 272 pages, 12 figures; v4: updated references and author list. Overall improvements and corrections. To appear in Classical and Quantum Gravity

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., Comment: 5 pages, 2 figures; accepted by PRL

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., Comment: 17 pages in journal format, 3 figures, 2 tables, accepted for publication in the Astrophysical Journal. Some references fixed and updated

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., Comment: 19 pages, 11 figures; accepted by Physical Review X

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., Comment: 6 pages, 3 figures; Contribution to the proceedings of "The Fourteenth Marcel Grossmann Meeting", University of Rome "La Sapienza", Rome, July 12-18, 2015

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., Comment: 23 pages, 9 figures. Invited review article, minor modifications to match published version

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., Comment: to appear in the Astrophysical Journal

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., Comment: 188 pages, 46 figures, 6 tables, 903 references. Matches version published in Classical and Quantum Gravity. Supplementary data files available at http://www.phy.olemiss.edu/~berti/research/ and http://centra.tecnico.ulisboa.pt/network/grit/files/

Before the 1970s, precision tests for gravity theories were constrained to the weak gravitational fields of the Solar system. Hence, only the weak-field slow-motion aspects of relativistic celestial mechanics could be investigated. Testing gravity beyond the first post-Newtonian contributions was for a long time out of reach. The discovery of the first binary pulsar by Russell Hulse and Joseph Taylor in the summer of 1974 initiated a completely new field for testing the relativistic dynamics of gravitationally interacting bodies. For the first time the back reaction of gravitational wave emission on the binary motion could be studied. Furthermore, the Hulse-Taylor pulsar provided the first test bed for the orbital dynamics of strongly self-gravitating bodies. To date there are a number of pulsars known, which can be utilized for precision test of gravity. Depending on their orbital properties and their companion, these pulsars provide tests for various different aspects of relativistic dynamics. Besides tests of specific gravity theories, like general relativity or scalar-tensor gravity, there are pulsars that allow for generic constraints on potential deviations of gravity from general relativity in the quasi-stationary strong-field and the radiative regime. This article presents a brief overview of this modern field of relativistic celestial mechanics, reviews some of the highlights of gravity tests with radio pulsars, and discusses their implications for gravitational physics and astronomy, including the upcoming gravitational wave astronomy., Comment: 80 pages, 22 figures, to appear in the Brumberg Festschrift, edited by S. M. Kopeikein, and to be published by de Gruyter, Berlin, 2014. Submitted to arXiv after final formatting for de Gruyter

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., Comment: 13 pages, 4 figures; accepted by Classical and Quantum Gravity

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., Comment: 25 pages, 9 figures; accepted by Classical and Quantum Gravity

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., Comment: This is the authors' version of the work. The definite version is published in Science Online, 26 April 2013, Vol: 340, Issue: 6131 doi: 10.1126/science.1233232. 54 pages, 18 figures, 4 tables

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., Comment: 3 pages, 1 figure; submitted to Proceedings of the Thirteenth Marcel Grossmann Meeting, (eds.) Kjell Rosquist, Robert T Jantzen, Remo Ruffini

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., Comment: 3 pages, 1 figure; submitted to Proceedings of IAU Symposium 291 -- Neutron Stars and Pulsars: Challenges and Opportunities after 80 years, (ed.) J. van Leeuwen

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), Comment: 34 pages, 8 figures, 2 tables; accepted by Classical and Quantum Gravity

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., Comment: 19 pages, two figures; in press in Classical and Quantum Gravity as an invited review on the "focus issue" on Tests of the Weak Equivalence Principle, organized by Clive Speake and Clifford Will, proof corrections added

(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)., Comment: Accepted for publication in MNRAS. 18 pages in emulate MNRAS format, 9 figures and 1 table

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)., Comment: 13 pages, 5 figures, Accepted by ApJ

(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., Comment: 15 pages in emulate MNRAS format, 4 tables, 8 figures. Accepted for publication in MNRAS

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., Comment: 3 pages, 1 figure, proceedings of the 12th Marcel Grossmann meeting

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., Comment: 3 pages, Proceedings of the 12th Marcel Grossmann Meeting on General Relativity (MG 12)

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., Comment: 6 pages, PostScript, 2 figures included. Submitted to Astronomy & Astrophysics. replaced by accepted version (minor corrections)