Pulsar In Orbit Research Articles

Magnetic Field Geometry of the Gamma-Ray Binary PSR B1259–63 Revealed via X-Ray Polarization

Some X-ray binaries containing an energetic pulsar in orbit around a normal star accelerate particles to high energies in the shock cone formed where the pulsar and stellar winds collide. The magnetic field geometry in the acceleration region in such binaries is unknown. We performed the first measurement of the polarization of the X-ray synchrotron emission from a gamma-ray-emitting binary system. We observed PSR B1259–63 with the Imaging X-ray Polarimetry Explorer during an X-ray bright phase following the periastron passage in 2024 June. X-ray polarization is detected with a polarization degree of 8.3% ± 1.5% at a significance of 5.3σ. The X-ray polarization angle is aligned with the axis of the shock cone at the time of the observation. This indicates that the predominant component of the magnetic field in the acceleration region is oriented perpendicular to the shock-cone axis.

Reciprocating Magnetic Fields in the Pulsar Wind Observed from the Black Widow Pulsar J1720-0534

We report the radio observations of the eclipsing black widow pulsar J1720−0534, a 3.26 ms pulsar in orbit with a low mass companion of mass 0.029 to 0.034 M ⊙. We obtain the phase-connected timing ephemeris and polarization profile of this millisecond pulsar (MSP) using the Five-hundred-meter Aperture Spherical radio Telescope (FAST), the Green Bank Telescope (GBT), and the Parkes Telescope. For the first time from such a system, an oscillatory polarization angle change was observed from a particular eclipse egress with partial depolarization, indicating 10-milliGauss-level reciprocating magnetic fields oscillating in a length scale of 5 × 103 km (assuming an orbital inclination angle of 90°) outside the companion’s magnetosphere. The dispersion measure variation observed during the ingresses and egresses shows the rapid raising of the electron density in the shock boundary between the companion’s magnetosphere and the surrounding pulsar wind. We suggest that the observed oscillatory magnetic fields originate from the pulsar wind outside the companion’s magnetosphere.

Green Bank Telescope Discovery of the Redback Binary Millisecond Pulsar PSR J0212+5321

We report the discovery of a 2.11 ms binary millisecond pulsar during a targeted search of the redback optical candidate coincident with the γ-ray source 3FGL J0212.5+5320 using the Robert C. Byrd Green Bank Telescope (GBT) with the Breakthrough Listen backend at L band. Over a seven month period, five pointings were made near inferior conjunction of the pulsar in its 20.9 hr orbit, resulting in two detections, lasting 12 and 42 minutes. The pulsar dispersion measure (DM) of 25.7 pc cm−3 corresponds to a distance of 1.15 kpc in the NE2001 Galactic electron density model, consistent with the Gaia parallax distance of 1.16 ± 0.03 kpc for the companion star. We suspect the pulsar experiences wide orbital eclipses, similar to other redbacks, as well as scintillation and DM delays caused by its interaction with its companion and surroundings. Although the pulsar was only detected over ≈3.7% of the orbit, its measured acceleration is consistent with published binary parameters from optical radial velocity spectroscopy and light-curve modeling of the companion star, and it provides a more precise mass ratio and a projected semimajor axis for the pulsar orbit. We also obtained a refined optical photometric orbit ephemeris, and observed variability of the tidally distorted companion over 7 yr. A hard X-ray light curve from NuSTAR shows expected orbit-modulated emission from the intrabinary shock. The pulsar parameters and photometric ephemeris greatly restrict the parameter space required to search for a coherent timing solution including pulsar spin-down rate, either using Fermi γ-rays or further radio pulse detections.

Testing space–time geometries and theories of gravity at the Galactic centre with pulsar’s time delay

ABSTRACTWe developed a numerical methodology to compute the fully relativistic propagation time of photons emitted by a pulsar in orbit around a massive compact object, like the supermassive black hole Sagittarius A* in the Galactic centre, whose gravitational field is described by a generic spherically symmetric space–time. Pulsars at the Galactic centre are usually regarded as the next major precision probe for theories of gravity, filling the current experimental gap between horizon-scale gravity tests and those at larger scales. We retain a completely general approach, which allows us to apply our code to the Schwarzschild space–time (by which we successfully validate our methodology) and to three different well-motivated alternatives to the standard black hole paradigm. The results of our calculations highlight departures spanning several orders of magnitudes in timing residuals, that are supposed to be detectable with future observing facilities like the Square Kilometer Array.

Variable scintillation arcs of millisecond pulsars observed with the Large European Array for Pulsars

ABSTRACT We present the first large sample of scintillation arcs in millisecond pulsars (MSPs), analysing 12 sources observed with the Large European Array for Pulsars (LEAP), and the Effelsberg 100-m telescope. We estimate the delays from multipath propagation, measuring significant correlated changes in scattering time-scales over a 10 yr time span. Many sources show compact concentrations of power in the secondary spectrum, which in PSRs J0613−0200 and J1600−3053 can be tracked between observations, and are consistent with compact scattering at fixed angular positions. Other sources such as PSRs J1643−1224 and J0621+1002 show diffuse, asymmetric arcs which are likely related to phase-gradients across the scattering screen. PSR B1937+21 shows at least three distinct screens which dominate at different times and evidence of varying screen axes or multiscreen interactions. We model annual and orbital arc curvature variations in PSR J0613−0200, providing a measurement of the longitude of ascending node, resolving the sense of the orbital inclination, where our best-fit model is of a screen with variable axis of anisotropy over time, corresponding to changes in the scattering of the source. Unmodelled variations of the screen’s axis of anisotropy are likely to be a limiting factor in determining orbital parameters with scintillation, requiring careful consideration of variable screen properties, or independent very long baseline interferometry (VLBI) measurements. Long-term scintillation studies such as this serve as a complementary tool to pulsar timing, to measure a source of correlated noise for pulsar timing arrays, solve pulsar orbits, and to understand the astrophysical origin of scattering screens.

Characterizing the Regular Orbits of Binary Pulsars: An Initial Prospection Study

The dynamical interaction between binary systems is crucial in understanding the nature of orbital motion under the influence of gravitational potential. In our study, we focused on investigating the effects of dynamical forces on the regularity of binary pulsar orbits, which represent a pure two-body system. To incorporate the necessary time dependence and have a regular 3-D axisymmetric potential, we utilized the Rebound package as a numerical integrator. This package integrates the motion of particles under the influence of gravity, allowing for changing orbital parameters at a given instant, and providing a variety of integrators to be used. By analyzing the regularity properties of binaries and their sensitivity to initial conditions, we gained insight into the importance of considering even small perturbations to the system, as they can lead to significant changes in its dynamics.

Relativistic effects in a mildly recycled pulsar binary: PSR J1952+2630

We report the results of timing observations of PSR J1952+2630, a 20.7 ms pulsar in orbit with a massive white dwarf companion. We performed six months of timing observations with the Arecibo radio telescope in 2020 and used data from FAST from 2021. Together with previously published data, this represents a total timing baseline of 11 yr since the discovery of the pulsar in 2010. For the first time, we present a polarimetric profile of the pulsar and determine its rotation measure (RM), − 145.79 ± 0.15 rad m−2. With the increased timing baseline, we obtain improved estimates for astrometric, spin, and binary parameters for this system. In particular, we obtain an imporvement of an order of magnitude on the proper motion, and, for the first time, we detect three post-Keplerian parameters in this system: the advance of periastron ω̇, the orbital decay Ṗb, and the Shapiro delay (measured in the form of the h3 parameter). With the detection of these relativistic effects, we constrain the pulsar mass to 1.20−0.29+0.28 M⊙ and the mass of its companion to 0.97−0.13+0.16 M⊙. The current value of Ṗb is consistent with the General Relativity expectation for the masses obtained using ω̇ and h3. The excess (4.2−73.1+70.2 fs s−1) represents a limit on the emission of dipolar gravitational waves (GWs) from this system. This results in a limit on the difference in effective scalar couplings for the pulsar and companion (predicted by scalar-tensor theories of gravity; STTs) of |αp − αc|< 4.8 × 10−3 (68% C.L.), which does not yield a competitive test for STTs. However, our simulations of future timing campaigns of this system, based on the timing precision we have achieved with FAST, show that by 2032, the precision of Ṗb and ω̇ will allow for much more precise masses and much tighter constraints on the orbital decay contribution from dipolar GWs, resulting in |αp − αc|< 1.3 × 10−3 (68% C.L.). For comparison, we obtain |αp − αc|< 1.9 × 10−3 and < 3.3 × 10−3 from PSR J1738+0333 and PSR J2222−0137, respectively. We also present the constraints this system will place on the {α0, β0} parameters of Damour-Esposito-Farèse (DEF) gravity by 2032. They are comparable to those of PSR J1738+0333. Unlike PSR J1738+0333, PSR J1952+2630 will not be limited in its mass measurement and has the potential to place even more restrictive limits on DEF gravity in the future. Further improvements to this test will likely be limited by uncertainties in the kinematic contributions to Ṗb because of the lack of precise distance measurements.

Analysis of the ionized interstellar medium and orbital dynamics of PSR J1909-3744 using scintillation arcs

ABSTRACT Long-term studies of binary millisecond pulsars (MSPs) provide precise tests of strong-field gravity, and can be used to measure neutron-star masses. PSR J1909-3744, a binary MSP has been the subject of several pulsar timing analyses. The edge-on orbit enables measurement of its mass using the Shapiro delay; however, there is degeneracy in the sense of the inclination angle, i, and multiple solutions for the longitude of ascending node, Ω. Radio pulsars scintillate due to inhomogeneities in the ionized interstellar medium (IISM). This can result in scintillation arcs in the power spectrum of the dynamic spectrum that can use these to study the interstellar medium and constrain binary pulsar orbits. Here, we study the scintillation of PSR J1909-3744 using observations from the 64-m Parkes Radio Telescope (Murriyang) over ≈13 yr, using techniques to study scintillation in a lower signal-to-noise regime. By monitoring annual and orbital variations of the arc-curvature measurements, we are able to characterize the velocity of the IISM. We find that the statistics of the IISM remained stationary over this time and a slightly anisotropic model (axial ratio ≳1.2) is preferred. We measure the relative distance to a single dominant thin scattering screen at s = 0.49 ± 0.04, or Ds = 590 ± 50 pc, with an angle of anisotropy ζ = 85 ± 6° (East of North) and velocity in the direction of anisotropy VIISM,ζ = 14 ± 10 km s−1. By combining a physical model of the IISM and current pulsar timing results, we also constrain Ω = 225 ± 3° and i = 86.46 ± 0.05°.

Photon Theory of Gravity – An Advance from Einstein’s Relativity

Based on a postulate that photons of low frequencies (undetectable by current technology) are the gravity force carrier, the paper derives quantitative results that are the same as or very similar to those derived in the special and general relativity theories and explains experiments and observations better. These quantitative results include the mass-energy formula, the energy momentum equation, and those for relative mass, the transverse Doppler effect, gravitational red shift, planetary precession, the deflection angle of light in gravitational lensing, the orbits around a black hole, and the strength and direction of gravitational waves (orbit decay of pulsars). Moreover, the explanations are different from those in Einstein’s relativity theory, such as the explanation of the null Doppler effect of electromagnetic waves reflected from a transversely moving surface, the reason for gravitational red shift, and the size of the light sphere around a black hole. The paper claims that both the high-order Doppler effect and the gravitational red shift occur only at the point of photon emission. The paper also explains why the predicted pulsar orbit decay is close but differs from calculations based on observations.

Binary Star System Decay by Graviton Interaction

The action of gravitons in a binary star system is modelled as the locus of points on an ellipse synchronous to the elliptic orbit of the binary star. In their interaction between the masses in the system the rotational energy of the gravitons is reduced by gravitational redshift, which accounts for the decay of the binary star orbital period. This model is able to fit a broad range of eccentricities of binary pulsar orbits and orbital period decay comparable to the General Relativistic gravitational wave model.

Radio absorption in high-mass gamma-ray binaries

High-mass γ-ray binaries consist of a presumptive pulsar in orbit with a massive star. The intense outflows from the star can absorb radio emission from the pulsar, making the detection of pulsation difficult. In this work, we present the basic geometry and formulae that describe the absorption process of a pulsar in binary with an O/B star and apply our model to two typical and well-studied binaries: PSR B1259−63/LS 2883 and LS 5039. We investigate the influences of the equatorial disc of LS 2883 with different orientations on the dispersion measure and free-free absorption of the radio pulsation from PSR B1259−63. The observed data are consistent with the disc inserted on the orbital plane with a relatively large inclination angle. For LS 5039, due to its tight orbit, it was believed that the strong wind absorption makes detecting radio emissions from the putative pulsar unlikely. However, considering the wind interaction and orbital motion, a bow shock cavity and a Coriolis shock would be formed, thereby allowing the pulsations to partially avoid stellar outflow absorption. We investigate the dependence of the radio optical depth on the observing frequencies, the orbital inclination angle, and the wind parameters. We suppose that the presumptive pulsar in LS 5039 is similar to PSR B1259−63 with pulsed emission extending to several tens of gigahertz. In that case, there could be a transparent window for radio pulsations when the pulsar is moving around the inferior conjunction. The following deep monitoring of LS 5039 and other systems by radio telescopes at high radio frequencies might reveal the nature of compact objects in the future. Alternatively, even a null detection could still provide further constraints on the properties of the putative pulsar and stellar outflows.

Refined Mass and Geometric Measurements of the High-mass PSR J0740+6620

Abstract We report results from continued timing observations of PSR J0740+6620, a high-mass, 2.8 ms radio pulsar in orbit with a likely ultracool white dwarf companion. Our data set consists of combined pulse arrival-time measurements made with the 100 m Green Bank Telescope and the Canadian Hydrogen Intensity Mapping Experiment telescope. We explore the significance of timing-based phenomena arising from general relativistic dynamics and variations in pulse dispersion. When using various statistical methods, we find that combining ∼1.5 yr of additional, high-cadence timing data with previous measurements confirms and improves on previous estimates of relativistic effects within the PSR J0740+6620 system, with the pulsar mass m p = 2.08 − 0.07 + 0.07 M ⊙ (68.3% credibility) determined by the relativistic Shapiro time delay. For the first time, we measure secular variation in the orbital period and argue that this effect arises from apparent acceleration due to significant transverse motion. After incorporating contributions from Galactic differential rotation and off-plane acceleration in the Galactic potential, we obtain a model-dependent distance of d = 1.14 − 0.15 + 0.17 kpc (68.3% credibility). This improved distance confirms the ultracool nature of the white dwarf companion determined from recent optical observations. We discuss the prospects for future observations with next-generation facilities, which will likely improve the precision on m p for J0740+6620 by an order of magnitude within the next few years.

Testing the black hole no-hair theorem with Galactic Center stellar orbits

Theoretical investigations have provided proof-of-principle calculations suggesting measurements of stellar or pulsar orbits near the Galactic Center could strongly constrain the properties of the Galactic Center black hole, local matter, and even the theory of gravity itself. In this work, we develop both a Markov chain Monte Carlo and an analytic model (Fisher matrix) to understand what properties are well-constrained and why. We conclude that existing astrometric measurements cannot constrain the spin of the Galactic Center black hole. Extrapolating to the precision and cadence of future experiments, we anticipate that the black hole spin can be measured with the known star S2. Our calculations show that we can measure the dimensionless black hole spin to a precision of $\sim$0.1 with weekly measurements of the orbit of S2 for 40 years using the GRAVITY telescope's best resolution at the Galactic Center. An analytic expression is derived for the measurement uncertainty of the black hole spin using the Fisher matrix in terms of observation strategy, star's orbital parameters, and instrument resolution. From it, we conclude that highly eccentric orbits can provide better constraints on the spin and that an orbit with a higher eccentricity is more favorable even when the orbital period is longer. We also apply it to S62, S4711, and S4714 and show whether they can constrain the black hole spin sooner than S2. If in addition, future measurements include the discovery of a new, tighter stellar orbit, then future data could conceivably enable tests of strong field gravity, by directly measuring the black hole quadrupole moment. Our simulations show that with a stellar orbit similar to that of S2 but at one-fifth the distance to the Galactic Center and GRAVITY's resolution limits, we can start to test the no-hair theorem with 20 years of weekly orbital measurements.

Radio timing in a millisecond pulsar – extreme/intermediate mass ratio binary system

Radio timing observations of a millisecond pulsar in orbit around the Galactic centre black hole (BH) or a BH at the centre of globular clusters could answer foundational questions in astrophysics and fundamental physics. Pulsar radio astronomy typically employs the post-Keplerian approximation to determine the system parameters. However, in the strong gravitational field around the central BH, higher order relativistic effects may become important. We compare the pulsar timing delays given by the post-Keplerian approximation with those given by a relativistic timing model. We find significant discrepancies between the solutions derived for the Einstein delay and the propagation delay (i.e. Roemer and Sharpiro delay) compared to the fully relativistic solutions. Correcting for these higher order relativistic effects is essential in order to construct accurate radio timing models for pulsar systems at the Galactic centre and the centre of globular clusters and informing issues related to their detection.

A revisit of PSR J1909−3744 with 15-yr high-precision timing

ABSTRACT We report on a high-precision timing analysis and an astrophysical study of the binary millisecond pulsar, PSR J1909−3744, motivated by the accumulation of data with well improved quality over the past decade. Using 15 yr of observations with the Nançay Radio Telescope, we achieve a timing precision of approximately 100 ns. We verify our timing results by using both broad-band and sub-band template matching methods to create the pulse time-of-arrivals. Compared with previous studies, we improve the measurement precision of secular changes in orbital period and projected semimajor axis. We show that these variations are both dominated by the relative motion between the pulsar system and the Solar system barycentre. Additionally, we identified four possible solutions to the ascending node of the pulsar orbit, and measured a precise kinetic distance of the system. Using our timing measurements and published optical observations, we investigate the binary history of this system using the stellar evolution code mesa, and discuss solutions based on detailed WD cooling at the edge of the WD age dichotomy paradigm. We determine the 3D velocity of the system and show that it has been undergoing a highly eccentric orbit around the centre of our Galaxy. Furthermore, we set up a constraint over dipolar gravitational radiation with the system, which is complementary to previous studies given the mass of the pulsar. We also obtain a new limit on the parametrized post-Newtonian parameter, $\left| \hat\alpha_1 \right|$ < 2.1 × 10−5 at 95 per cent confidence level, which is fractionally better than previous best published value and achieved with a more concrete method.

New graviton mass bound from binary pulsars

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.

Very Long Baseline Astrometry of PSR J1012+5307 and its Implications on Alternative Theories of Gravity

Abstract PSR J1012+5307, a millisecond pulsar in orbit with a helium white dwarf (WD), has been timed with high precision for about 25 yr. One of the main objectives of this long-term timing is to use the large asymmetry in gravitational binding energy between the neutron star and the WD to test gravitational theories. Such tests, however, will be eventually limited by the accuracy of the distance to the pulsar. Here, we present very long baseline interferometry (VLBI) astrometry results spanning approximately 2.5 yr for PSR J1012+5307, obtained with the Very Long Baseline Array as part of the project. These provide the first proper motion and absolute position for PSR J1012+5307 measured in a quasi-inertial reference frame. From the VLBI results, we measure a distance of kpc (all the estimates presented in the abstract are at 68% confidence) for PSR J1012+5307, which is the most precise obtained to date. Using the new distance, we improve the uncertainty of measurements of the unmodeled contributions to orbital period decay, which, combined with three other pulsars, places new constraints on the coupling constant for dipole gravitational radiation and the fractional time derivative of Newton’s gravitational constant in the local universe. As the uncertainties of the observed decays of orbital period for the four leading pulsar-WD systems become negligible in ≈10 yr, the uncertainties for and κ D will be improved to ≤1.5 × 10−13 yr−1 and ≤1.0 × 10−4, respectively, predominantly limited by the distance uncertainties.

Discovery of a 2.8 s Pulsar in a 2 Day Orbit High-mass X-Ray Binary Powering the Ultraluminous X-Ray Source ULX-7 in M51

Abstract We discovered 2.8 s pulsations in the X-ray emission of the ultraluminous X-ray source (ULX) M51 ULX-7 within the UNSEeN project, which was designed to hunt for new pulsating ULXs (PULXs) with XMM-Newton. The pulse shape is sinusoidal, and large variations of its amplitude were observed even within single exposures (pulsed fraction from less than 5% to 20%). Source M51 ULX-7 is variable, generally observed at an X-ray luminosity between 1039 and 1040 erg s−1, located in the outskirts of the spiral galaxy M51a at a distance of 8.6 Mpc. According to our analysis, the X-ray pulsar orbits in a 2 day binary with a projected semimajor axis 28 lt-s. For a neutron star (NS) of 1.4 M ⊙, this implies a lower limit on the companion mass of 8 M ⊙, placing the system hosting M51 ULX-7 in the high-mass X-ray binary class. The barycentric pulse period decreased by ≃0.4 ms in the 31 days spanned by our 2018 May–June observations, corresponding to a spin-up rate . In an archival 2005 XMM-Newton exposure, we measured a spin period of ∼3.3 s, indicating a secular spin-up of , a value in the range of other known PULXs. Our findings suggest that the system consists of a massive donor, possibly an OB giant or supergiant, and a moderately magnetic (dipole field component in the range 1012 G G) accreting NS with weakly beamed emission ( ).

Accretion‐induced collapse to third family compact stars as trigger for eccentric orbits of millisecond pulsars in binaries

AbstractA numerical rotating neutron star solver is used to study the temporal evolution of accreting neutron stars using a multi‐polytrope model for the nuclear equation of state named ACB5. The solver is based on a quadrupole expansion of the metric, but confirms the results of previous works, revealing the possibility of an abrupt transition of a neutron star from a purely hadronic branch to a third‐family branch of stable hybrid stars, passing through an unstable intermediate branch. The accretion is described through a sequence of stationary rotating stellar configurations which lose angular momentum through magnetic dipole emission, while, at the same time, gaining angular momentum through mass accretion. The model has several free parameters which are inferred from observations. The mass accretion scenario is studied in dependence on the effectiveness of angular momentum transfer which determines at which spin frequency the neutron star will become unstable against gravitational collapse to the corresponding hybrid star on the stable third‐family branch. It is conceivable that the neutrino burst which accompanies the deconfinement transition may trigger a pulsar kick which results in the eccentric orbit. A consequence of the present model is the prediction of a correlation between the spin frequency of the millisecond pulsar in the eccentric orbit and its mass at birth.

Phase-resolved gamma-ray spectroscopy of the Crab pulsar observed by POLAR

The POLAR detector is a space based Gamma-Ray Burst (GRB) polarimeter sensitive in the 15-500 keV energy range. Apart from its main scientific goal as a Gamma-Ray Burst polarimeter it is also able to detect photons from pulsars in orbit. By using the six-months in-orbit observation data, significant pulsation from the PSR B0531+21 (Crab pulsar) was obtained. In this work, we present the precise timing analysis of the Crab pulsar, together with a phase-resolved spectroscopic study using a joint-fitting method adapted for wide field of view instruments like POLAR. By using single power law fitting over the pulsed phase, we obtained spectral indices ranging from 1.718 to 2.315, and confirmed the spectral evolution in a reverse S shape which is homogenous with results from other missions over broadband. We will also show, based on the POLAR in-orbit performance and Geant4 Monte-Carlo simulation, the inferred capabilities of POLAR-2, the proposed follow-up mission of POLAR on board the China Space Station (CSS), for pulsars studies.