Explore the vast range of titles available.

Because of their extremely high densities, massive neutron stars can be used to test gravity. Based on spectroscopy of its white dwarf companion, Antoniadis et al. (p. 448 ) identified a millisecond pulsar as a neutron star twice as heavy as the Sun. The observed binary's orbital decay is consistent with that predicted by general relativity, ruling out previously untested strong-field phenomena predicted by alternative theories. The binary system has a peculiar combination of properties and poses a challenge to our understanding of stellar evolution. Observations of a pulsar confirm general relativity in the strong-field regime and reveal a perplexing stellar binary. Many physically motivated extensions to general relativity (GR) predict substantial deviations in the properties of spacetime surrounding massive neutron stars. We report the measurement of a 2.01 ± 0.04 solar mass ( M ☉ ) pulsar in a 2.46-hour orbit with a 0.172 ± 0.003 M ☉ 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.

Pulsar timing experiments require high fidelity template profiles in order to minimize the biases in pulse time-of-arrival (TOA) measurements and their uncertainties. Efforts to acquire more precise TOAs given fixed effective area of telescopes, finite receiver noise, and limited integration time have led pulsar astronomers to the solution of implementing ultra-wideband receivers. This solution, however, has run up against the problem that pulse profile shapes evolve with frequency, which raises the question of how to properly measure and analyze TOAs obtained using template-matching methods. This paper proposes a new method for one facet of this problem, that of template profile generation, and demonstrates it on the well-timed millisecond pulsar J1713+0747. Specifically, we decompose pulse profile evolution into a linear combination of basis eigenvectors, the coefficients of which change slowly with frequency such that their evolution is modeled simply by a sum of low degree piecewise polynomial spline functions. These noise-free, high fidelity, frequency-dependent templates can be used to make measurements of so-called \"wideband TOAs\" simultaneously with an estimate of the instantaneous dispersion measure. The use of wideband TOAs is becoming important for pulsar timing array experiments, as the volume of datasets comprised of conventional, subbanded TOAs are quickly becoming unwieldly for the Bayesian analyses needed to uncover latent gravitational wave signals. Although motivated by high precision timing experiments, our technique is applicable in more general pulsar observations.

Modeling of frequency-dependent effects, contributed by the turbulence in the free electron density of interstellar plasma, is required to enable the detection of the expected imprints from the stochastic gravitational-wave (GW) background in pulsar timing data. In this work, we present an investigation of temporal variations of interstellar medium for a set of millisecond pulsars (MSPs) with the upgraded GMRT aided by large fractional bandwidth at lower observing frequencies. Contrary to the conventional narrow-band analysis using a frequency invariant template profile, we applied \\(PulsePortraiture\\) based wide-band timing analysis while correcting for the evolution of the pulsar profile with frequency. Implementation of \\(PulsePortraiture\\) based wide-band timing method for the GMRT discovered MSPs to probe the DM variations resulted in a DM precision of \\(10^-4\\,pc~cm^-3\\). In general, we achieve similar DM and timing precision from wide-band timing compared to the narrow-band timing with matching temporal variations of DMs. This wide-band timing study of newly discovered MSPs over a wide frequency range highlights the effectiveness of profile-modeling at low frequencies and probes the potential of using them in pulsar timing array.

We present an algorithm for the simultaneous measurement of a pulse time-of-arrival (TOA) and dispersion measure (DM) from folded wideband pulsar data. We extend the prescription from Taylor (1992) to accommodate a general two-dimensional template \"portrait\", the alignment of which can be used to measure a pulse phase and DM. We show that there is a dedispersion reference frequency that removes the covariance between these two quantities, and note that the recovered pulse profile scaling amplitudes can provide useful information. We experiment with pulse modeling by using a Gaussian-component scheme that allows for independent component evolution with frequency, a \"fiducial component\", and the inclusion of scattering. We showcase the algorithm using our publicly available code on three years of wideband data from the bright millisecond pulsar J1824-2452A (M28A) from the Green Bank Telescope, and a suite of Monte Carlo analyses validates the algorithm. By using a simple model portrait of M28A we obtain DM trends comparable to those measured by more standard methods, with improved TOA and DM precisions by factors of a few. Measurements from our algorithm will yield precisions at least as good as those from traditional techniques, but is prone to fewer systematic effects and is without ad hoc parameters. A broad application of this new method for dispersion measure tracking with modern large-bandwidth observing systems should improve the timing residuals for pulsar timing array experiments, like the North American Nanohertz Observatory for Gravitational Waves.

Trojan asteroids are found in the equilateral triangle Lagrange points of the Sun-Jupiter system in great number, though they also exist less prolifically in other parts of the Solar System. Despite up to planetary mass Trojans being predicted in extrasolar systems (i.e. exotrojans), they remain unconfirmed, although strong candidate evidence has emerged recently. For the first time, we extend the search for exotrojans to radio pulsars with low-mass (\\(\\sim0.01\\,\\rm{M}_\\odot\\)) companions using accurately measured pulse times of arrival. With techniques developed for detecting the reflex motion of a star due to a librating Trojan, we place \\(\\sim 1\\,\\rm{M}_\\oplus\\) upper mass constraints on potential exotrojans around eight pulsars observed in the NANOGrav 15-year data set. We find weak evidence consistent with \\(\\sim2\\)--4\\(\\,\\rm{M}_{\\rm J}\\) exotrojans in the PSR~J0023+0923 and PSR~J1705\\(-\\)1903 binary systems, though the signals likely have a different, unknown source. We also place a libration-independent upper mass constraint of \\(\\sim8\\)\\,M\\(_{\\rm J}\\) on exotrojans in the PSR~J1641+8049 system by looking for an inconsistency between the times of superior conjunction as measured by optical light curves and those predicted by radio timing. These results offer initial observational constraints on the existence of exotrojans around pulsars, while their possible formation mechanisms remain unexplored.

Free electrons in the interstellar medium refract and diffract radio waves along multiple paths, resulting in angular and temporal broadening of radio pulses that limits pulsar timing precision. We determine multifrequency, multi-epoch scattering times for the large dispersion measure millisecond pulsar J1903+0327 by developing a three component model for the emitted pulse shape that is convolved with a best fit pulse broadening function (PBF) identified from a family of thin-screen and extended-media PBFs. We show that the scattering time, \\(\\tau\\), at a fiducial frequency of 1500 MHz changes by approximately 10% over a 5.5yr span with a characteristic timescale of approximately 100 days. We also constrain the spectral index and inner scale of the wavenumber spectrum of electron density variations along this line of sight. We find that the scaling law for \\(\\tau\\) vs. radio frequency is strongly affected by any mismatch between the true and assumed PBF or between the true and assumed intrinsic pulse shape. We show using simulations that refraction is a plausible cause of the epoch dependence of \\(\\tau\\), manifesting as changes in the PBF shape and \\(1/e\\) time scale. Finally, we discuss the implications of our scattering results on pulsar timing including time of arrival delays and dispersion measure misestimation.

Pulsar timing array experiments have recently uncovered evidence for a nanohertz gravitational wave background by precisely timing an ensemble of millisecond pulsars. The next significant milestones for these experiments include characterizing the detected background with greater precision, identifying its source(s), and detecting continuous gravitational waves from individual supermassive black hole binaries. To achieve these objectives, generating accurate and precise times of arrival of pulses from pulsar observations is crucial. Incorrect polarization calibration of the observed pulsar profiles may introduce errors in the measured times of arrival. Further, previous studies (e.g., van Straten 2013; Manchester et al. 2013) have demonstrated that robust polarization calibration of pulsar profiles can reduce noise in the pulsar timing data and improve timing solutions. In this paper, we investigate and compare the impact of different polarization calibration methods on pulsar timing precision using three distinct calibration techniques: the Ideal Feed Assumption (IFA), Measurement Equation Modeling (MEM), and Measurement Equation Template Matching (METM). Three NANOGrav pulsars-PSRs J1643\\(-\\)1224, J1744\\(-\\)1134, and J1909\\(-\\)3744-observed with the 800 MHz and 1.5 GHz receivers at the Green Bank Telescope (GBT) are utilized for our analysis. Our findings reveal that all three calibration methods enhance timing precision compared to scenarios where no polarization calibration is performed. Additionally, among the three calibration methods, the IFA approach generally provides the best results for timing analysis of pulsars observed with the GBT receiver system. We attribute the comparatively poorer performance of the MEM and METM methods to potential instabilities in the reference noise diode coupled to the receiver and temporal variations in the profile of the reference pulsar, respectively.

In 2018 an ultra-wide-bandwidth low-frequency (UWL) receiver was installed on the 64-m Parkes Radio Telescope enabling observations with an instantaneous frequency coverage from 704 to 4032 MHz. Here, we present the analysis of a three-year data set of 35 millisecond pulsars observed with the UWL by the Parkes Pulsar Timing Array (PPTA), using wideband timing methods. The two key differences compared to typical narrow-band methods are, firstly, generation of two-dimensional templates accounting for pulse shape evolution with frequency and, secondly, simultaneous measurements of the pulse time-of-arrival (ToA) and dispersion measure (DM). This is the first time that wideband timing has been applied to a uniform data set collected with a single large-fractional bandwidth receiver, for which such techniques were originally developed. As a result of our study, we present a set of profile evolution models and new timing solutions including initial noise analysis. Precision of our ToA and DM measurements is in the range of 0.005 \\(-\\) 2.08 \\(\\mu\\)s and (0.043\\(-\\)14.24)\\(\\times10^{-4}\\) cm\\(^{-3}\\) pc, respectively, with 94% of the pulsars achieving a median ToA uncertainty of less than 1 \\(\\mu\\)s.

In this work, we present polarization profiles for 23 millisecond pulsars observed at 820 MHz and 1500 MHz with the Green Bank Telescope as part of the NANOGrav pulsar timing array. We calibrate the data using Mueller matrix solutions calculated from observations of PSRs B1929+10 and J1022+1001. We discuss the polarization profiles, which can be used to constrain pulsar emission geometry, and present both the first published radio polarization profiles for nine pulsars and the discovery of very low intensity average profile components (\"microcomponents\") in four pulsars. Using the Faraday rotation measures, we measure for each pulsar and use it to calculate the Galactic magnetic field parallel to the line of sight for different lines of sight through the interstellar medium. We fit for linear and sinusoidal trends in time in the dispersion measure and Galactic magnetic field and detect magnetic field variations with a period of one year in some pulsars, but overall find that the variations in these parameters are more consistent with a stochastic origin.

The radio galaxy 3C 66B has been hypothesized to host a supermassive black hole binary (SMBHB) at its center based on electromagnetic observations. Its apparent 1.05-year period and low redshift (\\(\\sim0.02\\)) make it an interesting testbed to search for low-frequency gravitational waves (GWs) using Pulsar Timing Array (PTA) experiments. This source has been subjected to multiple searches for continuous GWs from a circular SMBHB, resulting in progressively more stringent constraints on its GW amplitude and chirp mass. In this paper, we develop a pipeline for performing Bayesian targeted searches for eccentric SMBHBs in PTA data sets, and test its efficacy by applying it on simulated data sets with varying injected signal strengths. We also search for a realistic eccentric SMBHB source in 3C 66B using the NANOGrav 12.5-year data set employing PTA signal models containing Earth term-only as well as Earth+Pulsar term contributions using this pipeline. Due to limitations in our PTA signal model, we get meaningful results only when the initial eccentricity \\(e_0<0.5\\) and the symmetric mass ratio \\(\\eta>0.1\\). We find no evidence for an eccentric SMBHB signal in our data, and therefore place 95% upper limits on the PTA signal amplitude of \\(88.1\\pm3.7\\) ns for the Earth term-only and \\(81.74\\pm0.86\\) ns for the Earth+Pulsar term searches for \\(e_0<0.5\\) and \\(\\eta>0.1\\). Similar 95% upper limits on the chirp mass are \\((1.98 \\pm 0.05) \\times 10^9\\,M_{\\odot}\\) and \\((1.81 \\pm 0.01) \\times 10^9\\,M_{\\odot}\\). These upper limits, while less stringent than those calculated from a circular binary search in the NANOGrav 12.5-year data set, are consistent with the SMBHB model of 3C 66B developed from electromagnetic observations.