Search

The observation of neutron stars enables the otherwise impossible study of fundamental physical processes. The timing of binary radio pulsars is particularly powerful, as it enables precise characterization of their (three-dimensional) positions and orbits. PSR J0437–4715 is an important millisecond pulsar for timing array experiments and is also a primary target for the Neutron Star Interior Composition Explorer (NICER). The main aim of the NICER mission is to constrain the neutron star equation of state by inferring the compactness ( M _p / R ) of the star. Direct measurements of the mass M _p from pulsar timing therefore substantially improve constraints on the radius R and the equation of state. Here we use observations spanning 26 yr from Murriyang, the 64 m Parkes radio telescope, to improve the timing model for this pulsar. Among the new precise measurements are the pulsar mass M _p = 1.418 ± 0.044 M _⊙ , distance D = 156.96 ± 0.11 pc, and orbital inclination angle i = 137.°506 ± 0.°016, which can be used to inform the X-ray pulse profile models inferred from NICER observations. We demonstrate that these results are consistent between multiple data sets from the Parkes Pulsar Timing Array (PPTA), each modeled with different noise assumptions. Using the longest available PPTA data set, we measure an apparent second derivative of the pulsar spin frequency and discuss how this can be explained either by kinematic effects due to the proper motion and radial velocity of the pulsar or excess low-frequency noise such as a gravitational-wave background.

Using the upGREAT instrument on board the Stratospheric Observatory for Infrared Astronomy, we imaged [C ii ] 157.74 and [O i ] 63.18 μ m line emission from a bright photodissociation region (PDR) associated with an ionized bubble located in the Nessie Nebula, a filamentary infrared dark cloud. A comparison with Australia Telescope Compact Array data reveals a classic photodissociation region (PDR) structure, with a uniform progression from ionized gas, to photodissociated gas, and to molecular gas from the bubble’s interior to its exterior. [O i ] line emission from the bubble’s PDR reveals self-absorption features. Toward a far-IR bright protostar, both [O i ] and [C ii ] show an absorption feature at a velocity of −18 km s ^−1 , the same velocity as an unrelated foreground molecular cloud. Since the gas density in typical molecular clouds is well below the [O i ] and [C ii ] critical densities, the excitation temperatures for both lines are low (∼20 K). The Meudon models demonstrate that the surface of a molecular cloud, externally illuminated by a standard G _0 = 1 interstellar radiation field, can produce absorption features in both transitions. Thus, the commonly observed [O i ] and [C ii ] self-absorption and absorption features plausibly arise from the subthermally excited, externally illuminated photodissociated envelopes of molecular clouds. The luminous young stellar object AGAL337.916-00.477, located precisely where the expanding bubble strikes the Nessie filament, is associated with two shock tracers: NH _3 (3,3) maser emission and SiO 2−1 emission, indicating an interaction between the bubble and the filament. The interaction of the expanding bubble with its parental dense filament has triggered star formation.

The noise in millisecond pulsar (MSP) timing data can include contributions from observing instruments, the interstellar medium, the solar wind, solar system ephemeris errors, and the pulsars themselves. The noise environment must be accurately characterized in order to form the null hypothesis from which signal models can be compared, including the signature induced by nanohertz-frequency gravitational waves (GWs). Here we describe the noise models developed for each of the MSPs in the Parkes Pulsar Timing Array (PPTA) third data release, which have been used as the basis of a search for the isotropic stochastic GW background. We model pulsar spin noise, dispersion measure variations, scattering variations, events in the pulsar magnetospheres, solar wind variability, and instrumental effects. We also search for new timing model parameters and detected Shapiro delays in PSR J0614−3329 and PSR J1902−5105. The noise and timing models are validated by testing the normalized and whitened timing residuals for Gaussianity and residual correlations with time. We demonstrate that the choice of noise models significantly affects the inferred properties of a common-spectrum process. Using our detailed models, the recovered common-spectrum noise in the PPTA is consistent with a power law with a spectral index of γ = 13/3, the value predicted for a stochastic GW background from a population of supermassive black hole binaries driven solely by GW emission.

Pulsar timing arrays aim to detect nanohertz-frequency gravitational waves (GWs). A background of GWs modulates pulsar arrival times and manifests as a stochastic process, common to all pulsars, with a signature spatial correlation. Here we describe a search for an isotropic stochastic gravitational-wave background (GWB) using observations of 30 millisecond pulsars from the third data release of the Parkes Pulsar Timing Array (PPTA), which spans 18 yr. Using current Bayesian inference techniques we recover and characterize a common-spectrum noise process. Represented as a strain spectrum ${h}_{c}=A{(f/1{\mathrm{yr}}^{-1})}^{\alpha }$ , we measure $A={3.1}_{-0.9}^{+1.3}\times {10}^{-15}$ and α = −0.45 ± 0.20, respectively (median and 68% credible interval). For a spectral index of α = −2/3, corresponding to an isotropic background of GWs radiated by inspiraling supermassive black hole binaries, we recover an amplitude of $A={2.04}_{-0.22}^{+0.25}\times {10}^{-15}$ . However, we demonstrate that the apparent signal strength is time-dependent, as the first half of our data set can be used to place an upper limit on A that is in tension with the inferred common-spectrum amplitude using the complete data set. We search for spatial correlations in the observations by hierarchically analyzing individual pulsar pairs, which also allows for significance validation through randomizing pulsar positions on the sky. For a process with α = −2/3, we measure spatial correlations consistent with a GWB, with an estimated false-alarm probability of p ≲ 0.02 (approx. 2 σ ). The long timing baselines of the PPTA and the access to southern pulsars will continue to play an important role in the International Pulsar Timing Array.