Probing the Interstellar Medium and Dark Matter with Pulsars

Pulsars are rapidly rotating, highly magnetised neutron stars which emit electromagnetic radiation from their magnetic poles in the form of highly collimated beams. Pulsars are known as a powerful tool to probe the interstellar medium (ISM) and its constituents in the Miky Way. In this thesis we focus on probing the non-baryonic entities in the Milky Way, namely interstellar magnetic fields and dark matter. The first part of the thesis is dedicated to the investigation of Galactic magnetic fields, which are a major agent in the dynamics and energy balance of the ISM, and general evolution of the Galaxy. Small-scale turbulent magnetic fields in the Milky Way can be probed by monitoring variations in the Faraday rotation of linearly polarised radiation of pulsars. Following this idea, we use high-cadence, low-frequency observations from a set of selected pulsars carried out with German LOw-Frequency ARray (LOFAR) stations. The method that is used to determine the Faraday rotation measures (RMs) of pulsar signals is the Bayesian generalised Lomb-Scargle periodogram technique, developed in this thesis. We find that measured RMs are strongly affected by the highly time-variable terrestrial ionosphere. We have mitigated the ionospheric contribution assuming a thin-layer model of the ionosphere. We conclude that within this approximation the ionospheric RM corrections are accurate to ~ 0.06 - 0.07 rad/m2, which defines our sensitivity towards long-term astrophysical RM variations. Following these results, we investigate the sensitivity to the turbulence in the magnetised ISM between the pulsar and observer. No astrophysically credible signal has been detected. We discuss implications of the non-detection and analyse the possibilities for future investigations. The second part of this thesis deals with dark matter - a matter which accounts for about a quarter of the energy density of the Universe, and the nature of which is still under debate. The ultralight scalar field dark matter is one of the compelling dark matter candidates, which leaves characteristic imprints in the times of arrival of radio pulses from pulsars. We search for traces of ultralight scalar-field dark matter in the Galaxy using the latest Parkes Pulsar Timing Array dataset that contains the times of arrival of 26 pulsars. No statistically significant signal has been detected. Therefore, we set an upper limit on the local dark matter density. The most stringent constraints are still one order of magnitude above the local dark-matter density inferred from kinematics of stars in the Milky Way. We conclude by discussing the prospects of detecting the fuzzy dark matter with future radio astronomical facilities

Parkes Pulsar Timing Array constraints on ultralight scalar-field dark matter

It is widely accepted that dark matter contributes about a quarter of the critical mass-energy density in our Universe. The nature of dark matter is currently unknown, with the mass of possible constituents spanning nearly one hundred orders of magnitude. The ultralight scalar field dark matter, consisting of extremely light bosons with m ∼ 10^(−22)  eV and often called “fuzzy” dark matter, provides intriguing solutions to some challenges at sub-Galactic scales for the standard cold dark matter model. As shown by Khmelnitsky and Rubakov, such a scalar field in the Galaxy would produce an oscillating gravitational potential with nanohertz frequencies, resulting in periodic variations in the times of arrival of radio pulses from pulsars. The Parkes Pulsar Timing Array (PPTA) has been monitoring 20 millisecond pulsars at two- to three-week intervals for more than a decade. In addition to the detection of nanohertz gravitational waves, PPTA offers the opportunity for direct searches for fuzzy dark matter in an astrophysically feasible range of masses. We analyze the latest PPTA data set which includes timing observations for 26 pulsars made between 2004 and 2016. We perform a search in this data set for evidence of ultralight dark matter in the Galaxy using Bayesian and Frequentist methods. No statistically significant detection has been made. We, therefore, place upper limits on the local dark matter density. Our limits, improving on previous searches by a factor of 2 to 5, constrain the dark matter density of ultralight bosons with m ≤ 10^(−23)  eV to be below 6  GeV cm^(−3) with 95% confidence in the Earth neighborhood. Finally, we discuss the prospect of probing the astrophysically favored mass range m ≳ 10^(−22)  eV with next-generation pulsar timing facilities

Rapidly rotating neutron star progenitors

Rotating proto-neutron stars can be important sources of gravitational waves to be searched for by present-day and future interferometric detectors. It was demonstrated by Imshennik that in extreme cases the rapid rotation of a collapsing stellar core may lead to fission and formation of a binary proto-neutron star which subsequently merges due to gravitational wave emission. In the present paper, we show that such dynamically unstable collapsing stellar cores may be the product of a former merger process of two stellar cores in a common envelope. We applied population synthesis calculations to assess the expected fraction of such rapidly rotating stellar cores which may lead to fission and formation of a pair of proto-neutron stars. We have used the BSE population synthesis code supplemented with a new treatment of stellar core rotation during the evolution via effective core-envelope coupling, characterized by the coupling time, $\tau_c$. The validity of this approach is checked by direct MESA calculations of the evolution of a rotating 15 $M_\odot$ star. From comparison of the calculated spin distribution of young neutron stars with the observed one, reported by Popov and Turolla, we infer the value $\tau_c \simeq 5 \times 10^5$ years. We show that merging of stellar cores in common envelopes can lead to collapses with dynamically unstable proto-neutron stars, with their formation rate being $\sim 0.1-1\%$ of the total core collapses, depending on the common envelope efficiency.Comment: 10 pages, 4 figures, accepted for publication in MNRA

Parkes Pulsar Timing Array constraints on ultralight scalar-field dark matter

It is widely accepted that dark matter contributes about a quarter of the critical mass-energy density in our Universe. The nature of dark matter is currently unknown, with the mass of possible constituents spanning nearly one hundred orders of magnitude. The ultralight scalar field dark matter, consisting of extremely light bosons with m ∼ 10^(−22)  eV and often called “fuzzy” dark matter, provides intriguing solutions to some challenges at sub-Galactic scales for the standard cold dark matter model. As shown by Khmelnitsky and Rubakov, such a scalar field in the Galaxy would produce an oscillating gravitational potential with nanohertz frequencies, resulting in periodic variations in the times of arrival of radio pulses from pulsars. The Parkes Pulsar Timing Array (PPTA) has been monitoring 20 millisecond pulsars at two- to three-week intervals for more than a decade. In addition to the detection of nanohertz gravitational waves, PPTA offers the opportunity for direct searches for fuzzy dark matter in an astrophysically feasible range of masses. We analyze the latest PPTA data set which includes timing observations for 26 pulsars made between 2004 and 2016. We perform a search in this data set for evidence of ultralight dark matter in the Galaxy using Bayesian and Frequentist methods. No statistically significant detection has been made. We, therefore, place upper limits on the local dark matter density. Our limits, improving on previous searches by a factor of 2 to 5, constrain the dark matter density of ultralight bosons with m ≤ 10^(−23)  eV to be below 6  GeV cm^(−3) with 95% confidence in the Earth neighborhood. Finally, we discuss the prospect of probing the astrophysically favored mass range m ≳ 10^(−22)  eV with next-generation pulsar timing facilities

The second data release from the European Pulsar Timing Array:IV. Implications for massive black holes, dark matter, and the early Universe

The European Pulsar Timing Array (EPTA) and Indian Pulsar Timing Array (InPTA) collaborations have measured a low-frequency common signal in the combination of their second and first data releases, respectively, with the correlation properties of a gravitational wave background (GWB). Such a signal may have its origin in a number of physical processes including a cosmic population of inspiralling supermassive black hole binaries (SMBHBs); inflation, phase transitions, cosmic strings, and tensor mode generation by the non-linear evolution of scalar perturbations in the early Universe; and oscillations of the Galactic potential in the presence of ultra-light dark matter (ULDM). At the current stage of emerging evidence, it is impossible to discriminate among the different origins. Therefore, for this paper, we consider each process separately, and investigated the implications of the signal under the hypothesis that it is generated by that specific process. We find that the signal is consistent with a cosmic population of inspiralling SMBHBs, and its relatively high amplitude can be used to place constraints on binary merger timescales and the SMBH-host galaxy scaling relations. If this origin is confirmed, this would be the first direct evidence that SMBHBs merge in nature, adding an important observational piece to the puzzle of structure formation and galaxy evolution. As for early Universe processes, the measurement would place tight constraints on the cosmic string tension and on the level of turbulence developed by first-order phase transitions. Other processes would require non-standard scenarios, such as a blue-tilted inflationary spectrum or an excess in the primordial spectrum of scalar perturbations at large wavenumbers. Finally, a ULDM origin of the detected signal is disfavoured, which leads to direct constraints on the abundance of ULDM in our Galaxy

The second data release from the European Pulsar Timing Array:IV. Implications for massive black holes, dark matter, and the early Universe

The European Pulsar Timing Array (EPTA) and Indian Pulsar Timing Array (InPTA) collaborations have measured a low-frequency common signal in the combination of their second and first data releases, respectively, with the correlation properties of a gravitational wave background (GWB). Such a signal may have its origin in a number of physical processes including a cosmic population of inspiralling supermassive black hole binaries (SMBHBs); inflation, phase transitions, cosmic strings, and tensor mode generation by the non-linear evolution of scalar perturbations in the early Universe; and oscillations of the Galactic potential in the presence of ultra-light dark matter (ULDM). At the current stage of emerging evidence, it is impossible to discriminate among the different origins. Therefore, for this paper, we consider each process separately, and investigated the implications of the signal under the hypothesis that it is generated by that specific process. We find that the signal is consistent with a cosmic population of inspiralling SMBHBs, and its relatively high amplitude can be used to place constraints on binary merger timescales and the SMBH-host galaxy scaling relations. If this origin is confirmed, this would be the first direct evidence that SMBHBs merge in nature, adding an important observational piece to the puzzle of structure formation and galaxy evolution. As for early Universe processes, the measurement would place tight constraints on the cosmic string tension and on the level of turbulence developed by first-order phase transitions. Other processes would require non-standard scenarios, such as a blue-tilted inflationary spectrum or an excess in the primordial spectrum of scalar perturbations at large wavenumbers. Finally, a ULDM origin of the detected signal is disfavoured, which leads to direct constraints on the abundance of ULDM in our Galaxy

Quality over Quantity: Optimizing pulsar timing array analysis for stochastic and continuous gravitational wave signals

The search for gravitational waves using Pulsar Timing Arrays (PTAs) is acomputationally expensive complex analysis that involves source-specific noisestudies. As more pulsars are added to the arrays, this stage of PTA analysiswill become increasingly challenging. Therefore, optimizing the number ofincluded pulsars is crucial to reduce the computational burden of dataanalysis. Here, we present a suite of methods to rank pulsars for use withinthe scope of PTA analysis. First, we use the maximization of thesignal-to-noise ratio as a proxy to select pulsars. With this method, we targetthe detection of stochastic and continuous gravitational wave signals. Next, wepresent a ranking that minimizes the coupling between spatial correlationsignatures, namely monopolar, dipolar, and Hellings & Downs correlations.Finally, we also explore how to combine these two methods. We test theseapproaches against mock data using frequentist and Bayesian hypothesis testing.For equal-noise pulsars, we find that an optimal selection leads to an increasein the log-Bayes factor two times steeper than a random selection for thehypothesis test of a gravitational wave background versus a common uncorrelatedred noise process. For the same test but for a realistic EPTA dataset, a subsetof 25 pulsars selected out of 40 can provide a log-likelihood ratio that is$89\%$ of the total, implying that an optimally selected subset of pulsars canyield results comparable to those obtained from the whole array. We expectthese selection methods to play a crucial role in future PTA data combinations.<br