The MPIfR-MeerKAT Galactic Plane Survey: II. The eccentric double neutron star system PSR J1208a-5936 and a neutron star merger rate update

The MPIfR-MeerKAT Galactic Plane survey at L-band (MMGPS-L) is the most sensitive pulsar survey in the Southern Hemisphere, providing 78 discoveries in an area of 900 sq. deg. Here, we present a follow-up study of one of these new discoveries, PSR J1208a-5936, a 28.71-ms recycled pulsar in a double neutron star system with an orbital period of Pb=0.632 days and an eccentricity of e=0.348, merging within the Hubble time. Through timing of almost one year of observations, we detected the relativistic advance of periastron ( Ì = 0.918(1) deg yra-1), resulting in a total system mass of Mt=2.586(5) M·. We also achieved low-significance constraints on the amplitude of the Einstein delay and Shapiro delay, in turn yielding constraints on the pulsar mass (Mp = 1.26a-0.25+0.13 M·), the companion mass (Mc = 1.32a-0.13+0.25 M·), and the inclination angle (i=57 ± 12). This system is highly eccentric compared to other Galactic field double neutron stars with similar periods, possibly hinting at a larger-than-usual supernova kick during the formation of the second-born neutron star. The binary will merge within 7.2(2) Gyr due to the emission of gravitational waves, making it a progenitor of the neutron star merger events seen by ground-based gravitational wave observatories. With the improved sensitivity of the MMGPS-L, we updated the Milky Way neutron star merger rate to be RMWnew = 25a-9+19 Myra-1 within 90% credible intervals, which is lower than previous studies based on known Galactic binaries owing to the lack of further detections despite the highly sensitive nature of the survey. This implies a local cosmic neutron star merger rate of Rlocalnew = 293a-103+222 Gpca-3 yra-1, which is consistent with LIGO and Virgo O3 observations. With this, we also predict the observation of 10a-4+8 neutron star merger events during the LIGO-Virgo-KAGRA O4 run. We predict the uncertainties on the component masses and the inclination angle will be reduced to 5- 10a-3 M· and 0.4 after two decades of timing, and that in at least a decade from now the detection of b and the sky proper motion will serve to make an independent constraint of the distance to the system

The MPIfR-MeerKAT Galactic Plane Survey II. The eccentric double neutron star system PSR J1208-5936 and a neutron star merger rate update

The MMGPS-L is the most sensitive pulsar survey in the Southern Hemisphere. We present a follow-up study of one of these new discoveries, PSR J1208-5936, a 28.71-ms recycled pulsar in a double neutron star system with an orbital period of Pb=0.632 days and an eccentricity of e=0.348. Through timing of almost one year of observations, we detected the relativistic advance of periastron (0.918(1) deg/yr), resulting in a total system mass of Mt=2.586(5) Mo. We also achieved low-significance constraints on the amplitude of the Einstein delay and Shapiro delay, in turn yielding constraints on the pulsar mass (Mp=1.26(+0.13/-0.25) Mo), the companion mass (Mc=1.32(+0.25/-0.13) Mo, and the inclination angle (i=57(2) degrees). This system is highly eccentric compared to other Galactic field double neutron stars with similar periods, possibly hinting at a larger-than-usual supernova kick during the formation of the second-born neutron star. The binary will merge within 7.2(2) Gyr due to the emission of gravitational waves. With the improved sensitivity of the MMGPS-L, we updated the Milky Way neutron star merger rate to be 25(+19/-9) Myr$^{-1}$ within 90% credible intervals, which is lower than previous studies based on known Galactic binaries owing to the lack of further detections despite the highly sensitive nature of the survey. This implies a local cosmic neutron star merger rate of 293(+222/-103} Gpc/yr, consistent with LIGO and Virgo O3 observations. With this, we predict the observation of 10(+8/-4) neutron star merger events during the LIGO-Virgo-KAGRA O4 run. We predict the uncertainties on the component masses and the inclination angle will be reduced to 5x10$^{-3}$ Mo and 0.4 degrees after two decades of timing, and that in at least a decade from now the detection of the shift in Pb and the sky proper motion will serve to make an independent constraint of the distance to the system

Comparing recent PTA results on the nanohertz stochastic gravitational wave background

The Australian, Chinese, European, Indian, and North American pulsar timing array (PTA) collaborations recently reported, at varying levels, evidence for the presence of a nanohertz gravitational wave background (GWB). Given that each PTA made different choices in modeling their data, we perform a comparison of the GWB and individual pulsar noise parameters across the results reported from the PTAs that constitute the International Pulsar Timing Array (IPTA). We show that despite making different modeling choices, there is no significant difference in the GWB parameters that are measured by the different PTAs, agreeing within $1\sigma$. The pulsar noise parameters are also consistent between different PTAs for the majority of the pulsars included in these analyses. We bridge the differences in modeling choices by adopting a standardized noise model for all pulsars and PTAs, finding that under this model there is a reduction in the tension in the pulsar noise parameters. As part of this reanalysis, we "extended" each PTA's data set by adding extra pulsars that were not timed by that PTA. Under these extensions, we find better constraints on the GWB amplitude and a higher signal-to-noise ratio for the Hellings and Downs correlations. These extensions serve as a prelude to the benefits offered by a full combination of data across all pulsars in the IPTA, i.e., the IPTA's Data Release 3, which will involve not just adding in additional pulsars, but also including data from all three PTAs where any given pulsar is timed by more than as single PTA.Comment: 21 pages, 9 figures, submitted to Ap

Comparing Recent Pulsar Timing Array Results on the Nanohertz Stochastic Gravitational-wave Background

The Australian, Chinese, European, Indian, and North American pulsar timing array (PTA) collaborations recently reported, at varying levels, evidence for the presence of a nanohertz gravitational-wave background (GWB). Given that each PTA made different choices in modeling their data, we perform a comparison of the GWB and individual pulsar noise parameters across the results reported from the PTAs that constitute the International Pulsar Timing Array (IPTA). We show that despite making different modeling choices, there is no significant difference in the GWB parameters that are measured by the different PTAs, agreeing within 1σ. The pulsar noise parameters are also consistent between different PTAs for the majority of the pulsars included in these analyses. We bridge the differences in modeling choices by adopting a standardized noise model for all pulsars and PTAs, finding that under this model there is a reduction in the tension in the pulsar noise parameters. As part of this reanalysis, we "extended" each PTA's data set by adding extra pulsars that were not timed by that PTA. Under these extensions, we find better constraints on the GWB amplitude and a higher signal-to-noise ratio for the Hellings–Downs correlations. These extensions serve as a prelude to the benefits offered by a full combination of data across all pulsars in the IPTA, i.e., the IPTA's Data Release 3, which will involve not just adding in additional pulsars but also including data from all three PTAs where any given pulsar is timed by more than a single PTA

Seismic hazard assessment in Central Asia: Outcomes from a site approach

In the process of updating existing PSHA maps in Central Asia, a first step is the evaluation of the seismic hazard in terms of macroseismic intensity by applying a data driven method. Following the Site Approach to Seismic Hazard Assessment (SASHA) [11], the evaluation of the probability of exceedance of any given intensity value over a fixed exposure time, is mainly based on the seismic histories available at different locations without requiring any a-priori assumption about seismic zonation. The effects of earthquakes not included in the seismic history can be accounted by propagating the epicentral information through a Intensity Prediction Equation developed for the analyzed area. In order to comply with existing building codes in the region that use macroseismic intensity instead of PGA, we evaluated the seismic hazard at 2911 localities using a macroseismic catalog composed by 5322 intensity data points relevant to 75 earthquakes in the magnitude range 4.6-8.3. The results show that for most of the investigated area the intensity having a probability of at least 10% to be exceeded in 50 years is VIII. The intensity rises to IX for some area struck by strong earthquakes in the past, like the Chou-Kemin-Chilik fault zone in northern Tien-Shan, between Kyrgyzstan and Kazakhstan, or in Gissar range between Tajikistan and Uzbekistan. These values are about one intensity unit less than those evaluated in the Global Seismic Hazard Assessment Program (GSHAP: Ulomov, The GSHAP Region 7 working group [29]). Moreover, hazard curves have been extracted for the main towns of Central Asia and the results compared with the estimates previously obtained. A good agreement has been found for Bishkek (Kyrgyzstan) and Dushanbe (Tajikistan), while a lower probability of occurrence of I=VIII has been obtained for Tashkent (Uzbekistan) and a larger one for I=IX in Almaty (Kazakhstan). (C) 2012 Elsevier Ltd. All rights reserved

Triple trouble with PSR J1618-3921: Mass measurements and orbital dynamics of an eccentric millisecond pulsar

International audiencePSR J1618-3921 is one of five known millisecond pulsars (MSPs) in eccentric orbits (eMPSs) located in the Galactic plane, whose formation is poorly understood. Earlier studies of these objects revealed significant discrepancies between observation and predictions from standard binary evolution scenarios of pulsar-Helium white dwarf binaries. We conducted observations with the L-band receiver of the MeerKAT radio telescope and the UWL receiver of the Parkes Murriyang radio telescope between 2019 and 2021. These data were added to archival observations. We perform an analysis of this joint 23-year-dataset. We use the recent observations to give a brief account of the emission properties of J1618-3921, including a Rotating Vector model fit of the linear polarisation position angle of the pulsar. The long timing baseline allowed for a highly significant measurement of the rate of advance of periastron of $\dot{\omega}$. We can only report a low significance detection of the orthometric Shapiro delay parameters $h_3$ and $\varsigma$, leading to mass estimates of the total and individual binary masses. We detect an unexpected change in the orbital period of, which is an order of magnitude larger and carries an opposite sign to what is expected from Galactic acceleration and the Shklovskii effect. We also detect a significant second derivative of the spin frequency. Furthermore, we report an unexpected, abrupt change of the mean pulse profile in June 2021 with unknown origin. We propose that the anomalous $\dot{P_b}$ and $\ddot{f}$ indicate an additional varying acceleration due to a nearby mass, i.e., the J1618-3921 binary system is likely part of a hierarchical triple. This finding suggests that at least some eMSPs might have formed in triple star systems. Although the uncertainties are large, the binary companion mass is consistent with the $P_b$ - $M_{WD}$ relation