Binary Pulsar System Research Articles

In what sense a neutron star-black hole binary is the holy grail for testing gravity?

Pulsars in binary systems have been very successful to test the validity of general relativity in the strong field regime [1-4]. So far, such binaries include neutron star-white dwarf (NS-WD) and neutron star-neutron star (NS-NS) systems. It is commonly believed that a neutron star-black hole (NS-BH) binary will be much superior for this purpose. But in what sense is this true? Does it apply to all possible deviations?

Growth of cosmic structure: Probing dark energy beyond expansion

The quantity and quality of cosmic structure observations have greatly accelerated in recent years, and further leaps forward will be facilitated by imminent projects. These will enable us to map the evolution of dark and baryonic matter density fluctuations over cosmic history. The way that these fluctuations vary over space and time is sensitive to several pieces of fundamental physics: the primordial perturbations generated by GUT-scale physics; neutrino masses and interactions; the nature of dark matter and dark energy. We focus on the last of these here: the ways that combining probes of growth with those of the cosmic expansion such as distance-redshift relations will pin down the mechanism driving the acceleration of the Universe.One way to explain the acceleration of the Universe is invoke dark energy parameterized by an equation of state w. Distance measurements provide one set of constraints on w, but dark energy also affects how rapidly structure grows; the greater the acceleration, the more suppressed the growth of structure. Upcoming surveys are therefore designed to probe w with direct observations of the distance scale and the growth of structure, each complementing the other on systematic errors and constraints on dark energy. A consistent set of results will greatly increase the reliability of the final answer.Another possibility is that there is no dark energy, but that General Relativity does not describe the laws of physics accurately on large scales. While the properties of gravity have been measured with exquisite precision at stellar system scales and densities, within our solar system and by binary pulsar systems, its properties in different environments are poorly constrained. To fully understand if General Relativity is the complete theory of gravity we must test gravity across a spectrum of scales and densities. Rapid developments in gravitational wave astronomy and numerical relativity are directed at testing gravity in the high curvature, high density regime. Cosmological evolution provides a polar opposite test bed, probing how gravity behaves in the lowest curvature, low density environments.There are a number of different implementations of astrophysically relevant modifications of gravity. Generically, the models are able to reproduce the distance measurements while at the same time altering the growth of structure. In particular, as detailed below, the Poisson equation relating over-densities to gravitational potentials is altered, and the potential that determines the geodesics of relativistic particles (such as photons) differs from the potential that determines the motion of non-relativistic particles. Upcoming surveys will exploit these differences to determine whether the acceleration of the Universe is due to dark energy or to modified gravity.To realize this potential, both wide field imaging and spectroscopic redshift surveys play crucial roles. Projects including DES, eBOSS, DESI, PFS, LSST, Euclid, and WFIRST are in line to map more than a 1000cubic-billion-light-year volume of the Universe. These will map the cosmic structure growth rate to 1% in the redshift range 0<z<2, over the last 3/4 of the age of the Universe.

The Confrontation between General Relativity and Experiment.

The status of experimental tests of general relativity and of theoretical frameworks for analyzing them is reviewed and updated. Einstein’s equivalence principle (EEP) is well supported by experiments such as the Eötvös experiment, tests of local Lorentz invariance and clock experiments. Ongoing tests of EEP and of the inverse square law are searching for new interactions arising from unification or quantum gravity. Tests of general relativity at the post-Newtonian level have reached high precision, including the light deflection, the Shapiro time delay, the perihelion advance of Mercury, the Nordtvedt effect in lunar motion, and frame-dragging. Gravitational wave damping has been detected in an amount that agrees with general relativity to better than half a percent using the Hulse-Taylor binary pulsar, and a growing family of other binary pulsar systems is yielding new tests, especially of strong-field effects. Current and future tests of relativity will center on strong gravity and gravitational waves.

A COMPREHENSIVE STUDY OF RELATIVISTIC GRAVITY USING PSR B1534+12

We present updated analyses of pulse profiles and their arrival-times from PSR B1534+12, a 37.9-ms radio pulsar in orbit with another neutron star. A high-precision timing model is derived from twenty-two years of timing data, and accounts for all astrophysical processes that systematically affect pulse arrival-times. Five "post-Keplerian" parameters are measured that represent relativistic corrections to the standard Keplerian quantities of the pulsar's binary orbit. These relativistic parameters are then used to test general relativity by comparing the measurements with their predicted values. We conclude that relativity theory is confirmed to within 0.17% of its predictions. Furthermore, we derive the following astrophysical results from our timing analysis: a distance of $d_{\rm GR} = 1.051 \pm 0.005$ kpc to the pulsar-binary system, by relating the "excess" orbital decay to Galactic parameters; evidence for pulse "jitter" in PSR B1534+12 due to short-term magnetospheric activity; and evolution in pulse-dispersion properties. As a secondary study, we also present several analyses on pulse-structure evolution and its connection to relativistic precession of the pulsar's spin axis. The precession-rate measurement yields a value of $\Omega_1^{\rm spin}$ = $0.59^{+0.12}_{-0.08}$ $^{\circ}$/year (68% confidence) that is consistent with expectations, and represents an additional test of relativistic gravity.

BINARY PULSAR SYSTEM PSR B1259-63/LS2883 AS A GAMMA-RAY EMITTER

Observations of the binary pulsar PSR B1259-63/LS2883 in the high energy and very high energy domains have revealed a few quite unusual features. One of the most puzzling phenomena is the bright GeV flare detected with Fermi/LAT in 2011 January, approximately one month after periastron passage. Since the maximum luminosity in the high energy band during the flare nearly achieved the level of the pulsar spin-down energy losses, it is likely that the particles, responsible for this emission component, had a strongly anisotropic distribution, which resulted in the emission enhancement. One of the most prolific scenarios for such an emission enhancement is the Doppler boosting, which is realized in sources with relativistic motions. Interestingly, a number of hydrodynamical simulations have predicted a formation of highly relativistic outflows in binary pulsar systems, therefore scenarios, involving relativistic boosting, are very natural for these systems. However a more detailed analysis of such a possibility, presented in this study, reveals certain limitations which put strict constraints on the maximum luminosity achievable in this scenario. These constraints render the "Doppler boosting" scenario to be less feasible, especially for the synchrotron models.

Modulated gamma-ray emission from compact millisecond pulsar binary systems

A significant amount of the millisecond pulsars has been discovered within binary systems. In several such binary systems the masses of the companion stars have been derived allowing to distinguish two classes of objects, called the Black Widow and the Redback binaries. Pulsars in these binary systems are expected to produce winds which, colliding with stellar winds, create conditions for acceleration of electrons. These electrons should interact with the anisotropic radiation from the companion stars producing gamma-ray emission modulated with the orbital period of the binary system. We consider the interaction of a millisecond pulsar (MSP) wind with a very inhomogeneous stellar wind from the companion star within binary systems of the Black Widow and Redback types. It is expected that the pulsar wind should mix efficiently with the inhomogeneous stellar wind. Electrons accelerated in such mixed, turbulent winds can interact with the magnetic field and also strong radiation from the companion star producing not only synchrotron radiation but also gamma-rays in the the Inverse Compton process. Applying numerical methods, we calculated the GeV-TeV gamma-ray spectra and the light curves expected from some millisecond pulsar binary systems. It is concluded that energetic millisecond pulsar binary systems create a new class of TeV gamma-ray sources which could be detectable by the future Cherenkov arrays (e.g. CTA) and possibly also by the extensive campaigns with the present arrays (HESS, MAGIC, VERITAS). However, gamma-ray emission from the millisecond pulsar binary systems is predicted to have different features than those observed in the case of massive TeV gamma-ray binaries such as LS I 303 61 or LS 5039. The maximum in the TeV gamma-ray orbital light curve should appear when the MSP is behind the companion star.

Accretion, ablation and propeller evolution in close millisecond pulsar binary systems

A model for the formation and evolution of binary millisecond radio pulsars in systems with low mass companions (< 0.1 Msun) is investigated using a binary population synthesis technique. Taking into account the non conservative evolution of the system due to mass loss from an accretion disk as a result of propeller action and from the companion via ablation by the pulsar, the transition from the accretion powered to rotation powered phase is investigated. It is shown that the operation of the propeller and ablation mechanisms can be responsible for the formation and evolution of black widow millisecond pulsar systems from the low mass X-ray binary phase at an orbital period of ~0.1 day. For a range of population synthesis input parameters, the results reveal that a population of black widow millisecond pulsars characterized by orbital periods as long as ~0.4 days and companion masses as low as ~0.005 Msun can be produced. The orbital periods and minimum companion mass of this radio millisecond pulsar population critically depend on the thermal bloating of the semi-degenerate hydrogen mass losing component, with longer orbital periods for a greater degree of bloating. Provided that the radius of the companion is increased by about a factor of 2 relative to a fully degenerate, zero temperature configuration, an approximate agreement between observed long orbital periods and theoretical modeling of hydrogen rich donors can be achieved. We find no discrepancy between the estimated birth rates for LMXBs and black widow systems, which on average are $\sim~1.3\times10^{-5} {\rm yr}^{-1}$ and $1.3\times10^{-7} {\rm yr}^{-1}$ respectively.

THEEINSTEIN@HOMESEARCH FOR RADIO PULSARS AND PSR J2007+2722 DISCOVERY

Einstein@Home aggregates the computer power of hundreds of thousands of volunteers from 193 countries, to search for new neutron stars using data from electromagnetic and gravitational-wave detectors. This paper presents a detailed description of the search for new radio pulsars using Pulsar ALFA survey data from the Arecibo Observatory. The enormous computing power allows this search to cover a new region of parameter space; it can detect pulsars in binary systems with orbital periods as short as 11 minutes. We also describe the first Einstein@Home discovery, the 40.8 Hz isolated pulsar PSR J2007+2722, and provide a full timing model. PSR J2007+2722's pulse profile is remarkably wide with emission over almost the entire spin period. This neutron star is most likely a disrupted recycled pulsar, about as old as its characteristic spin-down age of 404 Myr. However there is a small chance that it was born recently, with a low magnetic field. If so, upper limits on the X-ray flux suggest but can not prove that PSR J2007+2722 is at least ~ 100 kyr old. In the future, we expect that the massive computing power provided by volunteers should enable many additional radio pulsar discoveries.

PROPAGATION AND STABILITY OF SUPERLUMINAL WAVES IN PULSAR WINDS

Nonlinear electromagnetic waves with superluminal phase velocity can propagate in the winds around isolated pulsars, and around some pulsars in binary systems. Using a short-wavelength approximation, we find and analyze an integrable system of equations that govern their evolution in spherical geometry. A confined mode is identified that stagnates to finite pressure at large radius and can form a precursor to the termination shock. Using a simplified criterion, we find this mode is stable for most isolated pulsars, but may be unstable if the external pressure is high, such as in the pulsar wind nebulae in starburst galaxies and in W44. Pulsar winds in eccentric binary systems, such as PSR 1259-63, may go through phases with stable and unstable electromagnetic precursors, as well as phases in which the density is too high for these modes to propagate.

The High Time Resolution Universe Pulsar Survey – VII. Discovery of five millisecond pulsars and the different luminosity properties of binary and isolated recycled pulsars

This paper presents the discovery and timing parameters for five millisecond pulsars (MSPs), four in binary systems with probable white dwarf companions and one isolated, found in ongoing processing of the High Time Resolution Universe Pulsar Survey (HTRU). We also present high quality polarimetric data on four of them. These further discoveries confirm the high potential of our survey in finding pulsars with very short spin periods. At least two of these five MSPs are excellent candidates to be included in the Pulsar Timing Array projects. Thanks to the wealth of MSP discoveries in the HTRU survey, we revisit the question of whether the luminosity distributions of isolated and binary MSPs are different. Using the Cordes and Lazio distance model and our new and catalogue flux density measurements, we find that 41 of the 42 most luminous MSPs in the Galactic disk are in binaries and a statistical analysis suggests that the luminosity functions differ with 99.9% significance. We conclude that the formation process that leads to solitary MSPs affects their luminosities, despite their period and period derivatives being similar to those of pulsars in binary systems.

Vainshtein mechanism in binary pulsars

We compute the scalar gravitational radiation from a binary pulsar system in the simplest model that exhibits the Vainshtein mechanism. The mechanism is successful in screening the effect from scalar fields conformally coupled to matter, although gravitational radiation is less suppressed relative to its general relativity predictions than static fifth forces effects within the pulsar system. This is due to a combination of two effects: firstly the existence of monopole and dipole radiation; secondly the Vainshtein suppression comes from the hierarchy of scales between the inverse frequency scale and the Vainshtein radius, rather than the orbital radius of the pulsar system. Extensions of these results will have direct relevance to infrared modifications of gravity, such as massive gravity theories, which are known to exhibit a Vainshtein mechanism. Generalization to Galileon models with higher order interactions are likely to provide stronger constraints.

MONTE CARLO SIMULATION OF NEUTRON STAR MASSES

We perform a Monte Carlo simulation to investigate the distribution and evolution of 66 mass measured pulsars. We get the best fits to the distribution at 1.35 ± 0.27M⨀ (1σ confidence level). In addition, we notice bimodal distributions in 1.34 ± 0.15M⨀ and 1.48 ± 0.53M⨀, this can be led to the idea that radio pulsars in binary systems have recycled. Thus we divide the data according to the characteristic spin period into two groups, millisecond pulsars (MSPs), P Spin ≤ 20 ms and less recycled pulsars P Spin ≥ 20 ms , respectively. We show that the distributions of MSPs at 1.42 ± 0.36M⨀, and 1.32 ± 0.18M⨀ for less recycled pulsars. As such, the mass of MSPs are heavier than those in less recycled pulsars by ~ 0.1M⨀, since they accreting material from their companions. On the other hand, the formation of heavier pulsars from the accretion induced collapse of accreting white dwarfs, must be invoked.

Impact of the orbital uncertainties on the timing of pulsars in binary systems

The detection of pulsations from an X-ray binary is an unambiguous signature of the presence of a neutron star in the system. When the pulsations are missed in the radio band, their detection at other wavelengths, like X-ray or gamma-rays, requires orbital demodulation, since the length of the observations are often comparable to, or longer than the system orbital period. The detailed knowledge of the orbital parameters of binary systems plays a crucial role in the detection of the spin period of pulsars, since any uncertainty in their determination translates into a loss in the coherence of the signal during the demodulation process. In this paper, we present an analytical study aimed at unveiling how the uncertainties in the orbital parameters might impact on periodicity searches. We find a correlation between the power of the signal in the demodulated arrival time series and the uncertainty in each of the orbital parameters. This correlation is also a function of the pulsar frequency. We test our analytical results with numerical simulations, finding good agreement between them. Finally, we apply our study to the cases of LS 5039 and LS I +61 303 and consider the current level of uncertainties in the orbital parameters of these systems and their impact on a possible detection of a hosted pulsar. We also discuss the possible appearance of a sideband ambiguity in real data. The latter can occur when, due to the use of uncertain orbital parameters, the power of a putative pulsar is distributed in frequencies lying nearby the pulsar period. Even if the appearance of a sideband is already a signature of a pulsar component, it may introduce an ambiguity in the determination of its period. We present here a method to solve the sideband issue.

A SHAPIRO DELAY DETECTION IN THE BINARY SYSTEM HOSTING THE MILLISECOND PULSAR PSR J1910–5959A

PSR J1910-5959A is a binary pulsar with a helium white dwarf companion located about 6 arcmin from the center of the globular cluster NGC6752. Based on 12 years of observations at the Parkes radio telescope, the relativistic Shapiro delay has been detected in this system. We obtain a companion mass Mc = 0.180+/-0.018Msun (1sigma) implying that the pulsar mass lies in the range 1.1Msun <= Mp <= 1.5Msun. We compare our results with previous optical determinations of the companion mass, and examine prospects for using this new measurement for calibrating the mass-radius relation for helium white dwarfs and for investigating their evolution in a pulsar binary system. Finally we examine the set of binary systems hosting a millisecond pulsar and a low mass helium white dwarf for which the mass of both stars has been measured. We confirm that the correlation between the companion mass and the orbital period predicted by Tauris & Savonije reproduces the observed values but find that the predicted Mp - Pb correlation over-estimates the neutron star mass by about 0.5Msun in the orbital period range covered by the observations. Moreover, a few systems do not obey the observed Mp - Pb correlation. We discuss these results in the framework of the mechanisms that inhibit the accretion of matter by a neutron star during its evolution in a low-mass X-ray binary.

Testing General Relativity and gravitational physics using the LARES satellite

The discovery of the accelerating expansion of the Universe, thought to be driven by a mysterious form of `dark energy' constituting most of the Universe, has further revived the interest in testing Einstein's theory of General Relativity. At the very foundation of Einstein's theory is the geodesic motion of a small, structureless test-particle. Depending on the physical context, a star, planet or satellite can behave very nearly like a test-particle, so geodesic motion is used to calculate the advance of the perihelion of a planet's orbit, the dynamics of a binary pulsar system and of an Earth orbiting satellite. Verifying geodesic motion is then a test of paramount importance to General Relativity and other theories of fundamental physics. On the basis of the first few months of observations of the recently launched satellite LARES, its orbit shows the best agreement of any satellite with the test-particle motion predicted by General Relativity. That is, after modelling its known non-gravitational perturbations, the LARES orbit shows the smallest deviations from geodesic motion of any artificial satellite. LARES-type satellites can thus be used for accurate measurements and for tests of gravitational and fundamental physics. Already with only a few months of observation, LARES provides smaller scatter in the determination of several low-degree geopotential coefficients (Earth gravitational deviations from sphericity) than available from observations of any other satellite or combination of satellites.

Seven pulsars in binary systems above the spin-up line

AbstractUsing data from the ATNF pulsar catalogue, 186 binary pulsars are shown in the magnetic field versus spin period (B-P) diagram, and their relationship to the spin-up line is investigated. Generally speaking, pulsars in binary systems should be below the spin-up line when they get enough accretion mass from their companions. It is found that there are seven binary pulsars above the spin-up line. Based on the parameters of these seven binary systems, we describe possible reasons why they are above the spin-up line.

Minimum accretion rate for millisecond pulsar formation in binary system

Abstract186 binary pulsars are shown in the magnetic field versus spin period (B-P) diagram, and their relations to the millisecond pulsars can be clearly seen. We declaim a minimum accretion rate for the millisecond pulsar formation both from the observation and theory. If the accretion rate is lower than the minimum accretion rate, the pulsar in binary system will not become a millisecond pulsar after the evolution.

Gamma-ray emission from pulsar binaries

AbstractPulsar winds, containing charged particles, waves and a net (phase-averaged) magnetic field, are thought to fuel the high-energy emission from several gamma-ray binaries. They terminate where the ram pressure matches that of the surroundings - which, in binaries, is provided by the wind of the companion. Before termination, pulsed emission can be produced by inverse Compton scattering of photons from the companion by particles in the waves. After termination, both the bulk kinetic energy of the particles and the Poynting flux in the waves are dissipated into an energetic particle population embedded in the surviving phase-averaged magnetic field. Pulsed emission is no longer possible, but a substantial flux of unpulsed high-energy photons can be produced. I will present results showing that the physical conditions at the termination shock can be divided into two regimes: a high density one, where current sheets in the wind are first compressed by an MHD shock and subsequently dissipate by reconnection, and a low density one, where the wind can first convert into an electromagnetic wave in the shock precursor, which then damps and merges into the wind nebula. The shocks surrounding isolated pulsars fall into the low-density category, but those around pulsars in binary systems, may transit from one regime to the other according to binary phase. The implications of the shock-structure dichotomy for these objects will be discussed.

Neutron star masses

AbstractNeutron star masses can be inferred from observations of binary pulsar systems, particularly by the measurement of relativistic phenomena within these orbits. The observed distribution of masses can be used to infer or constrain the equation of state for nuclear matter and to study astrophysical processes such as supernovae and binary star evolution. In this talk, I will review our present understanding of the neutron star mass distribution with an emphasis on the observational data.

GAMMA-RAY SIGNAL FROM THE PULSAR WIND IN THE BINARY PULSAR SYSTEM PSR B1259–63/LS 2883

Binary pulsar systems emit potentially detectable components of gamma ray emission due to Comptonization of the optical radiation of the companion star by relativistic electrons of the pulsar wind, both before and after termination of the wind. The recent optical observations of binary pulsar system PSR B1259-63/LS 2883 revealed radiation properties of the companion star which differ significantly from previous measurements. In this paper we study the implications of these observations for the interaction rate of the unshocked pulsar wind with the stellar photons and the related consequences for fluxes of high energy (HE) and very high energy (VHE) gamma rays. We show that the signal should be strong enough to be detected with Fermi close to the periastron passage, unless the pulsar wind is strongly anisotropic or the Lorentz factor of the wind is smaller than $10^3$ or larger that $10^5$. The higher luminosity of the optical star also has two important implications: (i) attenuation of gamma rays due to photon-photon pair production, and (ii) Compton drag of the unshocked wind. While the first effect has an impact on the lightcurve of VHE gamma rays, the second effect may significantly decrease the energy available for particle acceleration after termination of the wind.