Search

With its extreme density of stars and stellar remnants, dense young massive clusters, high specific star formation rate, intense radiation field, high magnetic field strength, and properties of the interstellar medium that resemble those in high redshift galaxies and starbursts, the Galactic Centre is the most extreme environment that we can observe in detail. It is also the only nucleus of a galaxy that we can observe with a resolution of just a few milli parsecs. This makes it a crucial target to understand the physics of galactic nuclei and star formation, as well as the connection between them. It enables studies of a large number of otherwise rare objects, such as extremely massive stars and stellar remnants, at a well-defined distance, thus facilitating the interpretation of their properties. The Galactic Centre has been and is being studied intensively with the most advanced facilities. In this White Paper, we advocate for a large-area, multi-wavelength survey with the Square Kilometre Array of an area of about 1.25x0.3 deg**2 (180x40 pc**2), centered on the massive black hole Sagittarius A* and for repeated deep observations of the nuclear star cluster over a decade, which will allow the community to address multiple science problems with a single data set.

We observe systematic profile changes in the visible pulsar of the compact double neutron star system PSR~J1946+2052 using observations with the Five-hundred-meter Aperture Spherical radio Telescope (FAST). The interpulse of PSR~J1946+2052 changed from single-peak to double-peak shape from 2018 to 2021. We attribute this evolution as the result of the relativistic spin precession of the pulsar. With the high sensitivity of FAST, we also measure significant polarization for the first time, allowing us to model this with the precessional rotating vector model. Assuming, to the first order, a circular hollow-cone-like emission beam pattern and taking the validity of general relativity, we derive the binary's orbital inclination angle (${63^\circ}^{+5^\circ}_{-3^\circ}$) and pulsar's spin geometry. Pulsar's spin vector and the orbital angular momentum vector are found to be only slightly misaligned (${0.21^\circ}^{+0.28^\circ}_{-0.10^\circ}$).The quoted uncertainties do not reflect the systematic uncertainties introduced by our model assumptions. By simulating future observations of profile and polarization evolution, we estimate that we could constrain the precession rate within a $43\%$ uncertainty in 9 years. Hence, we suggest that the system's profile evolution could be combined with precise pulsar timing to test general relativity in the future., Comment: 12 pages, 9 figures, accepted for publication in ApJ

The shape and polarisation properties of the radio pulse profiles of radio-loud magnetars provide a unique opportunity to investigate their magnetospheric properties. Gaussian Process Regression analysis was used to investigate the variation in the total intensity shape of the radio pulse profiles of the magnetar Swift J1818.0-1607. The observed profile shape was found to evolve through three modes between MJDs 59104 and 59365. The times at which these transitions occurred coincided with changes in the amplitude of modulations in the spin-down rate. The amount of linear and circular polarisation was also found to vary significantly with time. Lomb-Scargle periodogram analysis of the spin-down rate revealed three possibly harmonically related frequencies. This could point to the magnetar experiencing seismic activity. However, no profile features exhibited significant periodicity, suggesting no simple correlations between the profile variability and fluctuations of the spin-down on shorter timescales within the modes. Overall, this implies the mode changes seen are a result of local magnetospheric changes, with other theories, such as precession, less able to explain these observations., Comment: 11 pages, 6 figures, accepted for publication in MNRAS

Magnetars are highly-magnetised rotating neutron stars that are predominantly observed as high-energy sources. Six of this class of neutron star are known to also emit radio emission, and magnetars are, thus, a favoured model for the origin for at least some of the Fast Radio Bursts (FRBs). If magnetars, or neutron stars in general, are indeed responsible, sharp empirical constraints on the mechanism producing radio emission are required. Here we report on the detection of polarised quasi-periodic sub-structure in the emission of all well-studied radio-detected magnetars. A correlation previously seen, relating sub-structure in pulsed emission of radio emitting neutron stars to their rotational period, is extended, and shown to now span more than six of orders of magnitude in pulse period. This behaviour is not only seen in magnetars but in members of all classes of radio-emitting rotating neutron stars, regardless of their evolutionary history, their power source or their inferred magnetic field strength. If magnetars are responsible for FRBs, it supports the idea of being able to infer underlying periods from sub-burst timescales in FRBs., Comment: published in Nature Astronomy, 48 pages (main text, methods, extended data figures & tables and supplementary information), in total 17 figures, 4 tables

The Event Horizon Telescope (EHT) observed in 2017 the supermassive black hole at the center of the Milky Way, Sagittarius A* (Sgr A*), at a frequency of 228.1 GHz ($\lambda$=1.3 mm). The fundamental physics tests that even a single pulsar orbiting Sgr A* would enable motivate searching for pulsars in EHT datasets. The high observing frequency means that pulsars - which typically exhibit steep emission spectra - are expected to be very faint. However, it also negates pulse scattering, an effect that could hinder pulsar detections in the Galactic Center. Additionally, magnetars or a secondary inverse Compton emission could be stronger at millimeter wavelengths than at lower frequencies. We present a search for pulsars close to Sgr A* using the data from the three most-sensitive stations in the EHT 2017 campaign: the Atacama Large Millimeter/submillimeter Array, the Large Millimeter Telescope and the IRAM 30 m Telescope. We apply three detection methods based on Fourier-domain analysis, the Fast-Folding-Algorithm and single pulse search targeting both pulsars and burst-like transient emission; using the simultaneity of the observations to confirm potential candidates. No new pulsars or significant bursts were found. Being the first pulsar search ever carried out at such high radio frequencies, we detail our analysis methods and give a detailed estimation of the sensitivity of the search. We conclude that the EHT 2017 observations are only sensitive to a small fraction ($\lesssim$2.2%) of the pulsars that may exist close to Sgr A*, motivating further searches for fainter pulsars in the region., Comment: 33 pages, 7 figures, 6 Tables. Accepted for publication in ApJ

We report detections of scintillation arcs for pulsars in globular clusters M5, M13 and M15 for the first time using the Five-hundred-meter Aperture Spherical radio Telescope (FAST). From observations of these arcs at multiple epochs, we infer that screen-like scattering medium exists at distances $4.1_{-0.3}^{+0.2}$ kpc, $6.7_{-0.2}^{+0.2}$ kpc and $1.3_{-1.0}^{+0.7}$ kpc from Earth in the directions of M5, M13 and M15, respectively. This means M5's and M13's scattering screens are located at $3.0_{-0.2}^{+0.1}$ kpc and $4.4_{-0.1}^{+0.1}$ kpc above the galactic plane, whereas, M15's is at $0.6_{-0.5}^{+0.3}$ kpc below the plane. We estimate the scintillation timescale and decorrelation bandwidth for each pulsar at each epoch using the one-dimensional auto-correlation in frequency and time of the dynamic spectra. We found that the boundary of the Local Bubble may have caused the scattering of M15, and detected the most distant off-plane scattering screens to date through pulsar scintillation, which provides evidence for understanding the medium circulation in the Milky Way., Comment: 9 pages, 3 figures, matches the published version

We present 294 pulsars found in GeV data from the Large Area Telescope (LAT) on the Fermi Gamma-ray Space Telescope. Another 33 millisecond pulsars (MSPs) discovered in deep radio searches of LAT sources will likely reveal pulsations once phase-connected rotation ephemerides are achieved. A further dozen optical and/or X-ray binary systems co-located with LAT sources also likely harbor gamma-ray MSPs. This catalog thus reports roughly 340 gamma-ray pulsars and candidates, 10% of all known pulsars, compared to $\leq 11$ known before Fermi. Half of the gamma-ray pulsars are young. Of these, the half that are undetected in radio have a broader Galactic latitude distribution than the young radio-loud pulsars. The others are MSPs, with 6 undetected in radio. Overall, >235 are bright enough above 50 MeV to fit the pulse profile, the energy spectrum, or both. For the common two-peaked profiles, the gamma-ray peak closest to the magnetic pole crossing generally has a softer spectrum. The spectral energy distributions tend to narrow as the spindown power $\dot E$ decreases to its observed minimum near $10^{33}$ erg s$^{-1}$, approaching the shape for synchrotron radiation from monoenergetic electrons. We calculate gamma-ray luminosities when distances are available. Our all-sky gamma-ray sensitivity map is useful for population syntheses. The electronic catalog version provides gamma-ray pulsar ephemerides, properties and fit results to guide and be compared with modeling results., Comment: 142 pages. Accepted by the Astrophysical Journal Supplement

We discover three new weak pulse components in two known pulsars, one in PSR J0304+1932 and two in PSR J1518+4904. These components are emitted about half way between the main emission beam and the interpulse beam (beam from the opposite pole). They are separated from their main pulse peak by $99^{\circ}\pm{3}^{\circ}$ for J0304+1932, $123^{\circ}.6\pm{0^{\circ}.7}$ (leading) and $93^{^{\circ}}\pm 0^{\circ}.4$ (trailing) for J1518+4904, respectively. Their peak-intensity ratios to main pulses are: $\sim$ 0.06% for J0304+1932, $\sim$ 0.17% and $\sim$ 0.83% for J1518+4904. We also analyzed flux fluctuation and profile variation of the emissions for two pulsars. The results show correlations between the weak pulses and their main pulses, indicating that these emissions come from the same pole. We estimated the emission altitude of these weak pulses and derived a height of about half of the pulsar's light-cylinder radius. These pulse components are a unique sample of high-altitude emissions from pulsars, and challenge the current pulsar emission models., Comment: 10 pages, 13 figures, Accepted for publication in ApJ

BlackHoleCam is a project funded by a European Research Council Synergy Grant to build a complete astrophysical description of nearby supermassive black holes by using a combination of radio imaging, pulsar observations, stellar astrometry and general relativistic magneto-hydrodynamic models. BlackHoleCam scientists are active partners of the Event Horizon Telescope Consortium. In this talk I will discuss the use of pulsars orbiting Sagittarius A* for tests of General Relativity, the current difficulties in detecting such sources, recent results from the Galactic Centre magnetar PSR J1745-2900 and how BlackHoleCam aims to search for undiscovered pulsars in the Galactic Centre., Comment: 4 pages, 3 figures, Proceedings of The Fifteenth Marcel Grossmann Meeting on General Relativity, University of Rome La Sapienza, Italy, 1 - 7 July 2018

Interpretation of resolved polarized images of black holes by the Event Horizon Telescope (EHT) requires predictions of the polarized emission observable by an Earth-based instrument for a particular model of the black hole accretion system. Such predictions are generated by general relativistic radiative transfer (GRRT) codes, which integrate the equations of polarized radiative transfer in curved spacetime. A selection of ray-tracing GRRT codes used within the EHT collaboration is evaluated for accuracy and consistency in producing a selection of test images, demonstrating that the various methods and implementations of radiative transfer calculations are highly consistent. When imaging an analytic accretion model, we find that all codes produce images similar within a pixel-wise normalized mean squared error (NMSE) of 0.012 in the worst case. When imaging a snapshot from a cell-based magnetohydrodynamic simulation, we find all test images to be similar within NMSEs of 0.02, 0.04, 0.04, and 0.12 in Stokes I, Q, U , and V respectively. We additionally find the values of several image metrics relevant to published EHT results to be in agreement to much better precision than measurement uncertainties., Comment: Accepted for publication in ApJ

We present the first sub-arcsecond localised Fast Radio Burst (FRB) detected using MeerKAT. FRB 20210405I was detected in the incoherent beam using the MeerTRAP pipeline on 2021 April 05 with a signal to noise ratio of 140.8 and a dispersion measure of 565.17 pc cm$^{-3}$. It was detected while MeerTRAP was observing commensally with the ThunderKAT large survey project, and was sufficiently bright that we could use the ThunderKAT 8s images to localise the FRB. Two different models of the dispersion measure in the Milky Way and halo suggest that the source is either right at the edge of the Galaxy, or outside. This highlights the uncertainty in the Milky Way dispersion measure models, particularly in the Galactic Plane, and the uncertainty of Milky Way halo models. Further investigation and modelling of these uncertainties will be facilitated by future detections and localisations of nearby FRBs. We use the combined localisation, dispersion measure, scattering, specific luminosity and chance coincidence probability information to find that the origin is most likely extra-galactic and identify the likely host galaxy of the FRB: 2MASS J1701249$-$4932475. Using SALT spectroscopy and archival observations of the field, we find that the host is a disk/spiral galaxy at a redshift of $z=0.066$., Comment: 15 pages, 4 tables, 10 figures. Accepted to MNRAS

Interstellar scattering (ISS) of radio pulsar emission can be used as a probe of the ionised interstellar medium (IISM) and causes corruptions in pulsar timing experiments. Two types of ISS phenomena (intensity scintillation and pulse broadening) are caused by electron density fluctuations on small scales (< 0.01 AU). Theory predicts that these are related, and both have been widely employed to study the properties of the IISM. Larger scales ($\sim$1-100\,AU) cause measurable changes in dispersion and these can be correlated with ISS observations to estimate the fluctuation spectrum over a very wide scale range. IISM measurements can often be modeled by a homogeneous power-law spatial spectrum of electron density with the Kolmogorov ($-11/3$) spectral exponent. Here we aim to test the validity of using the Kolmogorov exponent with PSR~J0826+2637. We do so using observations of intensity scintillation, pulse broadening and dispersion variations across a wide fractional bandwidth (20 -- 180\,MHz). We present that the frequency dependence of the intensity scintillation in the high frequency band matches the expectations of a Kolmogorov spectral exponent but the pulse broadening in the low frequency band does not change as rapidly as predicted with this assumption. We show that this behavior is due to an inhomogeneity in the scattering region, specifically that the scattering is dominated by a region of transverse size $\sim$40\,AU. The power spectrum of the electron density, however, maintains the Kolmogorov spectral exponent from spatial scales of 5$\times10^{-6}$\,AU to $\sim$100\,AU., Comment: Accepted for publication in MNRAS, typo fixed

PSR J2140$-$2311B is a 13-ms pulsar discovered in 2001 in a 7.8-hour Green Bank Telescope (GBT) observation of the core-collapsed globular cluster M30 and predicted to be in a highly eccentric binary orbit. This pulsar has eluded detection since then, therefore its precise orbital parameters have remained a mystery until now. In this work, we present the confirmation of this pulsar using observations taken with the UHF receivers of the MeerKAT telescope as part of the TRAPUM Large Survey Project. Taking advantage of the beamforming capability of our backends, we have localized it, placing it $1.2(1)^\prime$ from the cluster centre. Our observations have enabled the determination of its orbit: it is highly eccentric ($e = 0.879$) with an orbital period of $6.2$ days. We also measured the rate of periastron advance, $\dot{\omega} = 0.078 \pm 0.002\, \rm deg \, yr^{-1}$. Assuming that this effect is fully relativistic, general relativity provides an estimate of the total mass of the system, $M_{\rm TOT} = 2.53 \pm 0.08$ M$_{\odot}$, consistent with the lightest double neutron star systems known. Combining this with the mass function of the system gives the pulsar and companion masses of $m_p < 1.43 \, \rm M_{\odot}$ and $m_c > 1.10 \, \rm M_{\odot}$ respectively. The massive, undetected companion could either be a massive WD or a NS. M30B likely formed as a result of a secondary exchange encounter. Future timing observations will allow the determination of a phase-coherent timing solution, vastly improving our uncertainty in $\dot{\omega}$ and likely enabling the detection of additional relativistic effects which will determine $m_p$ and $m_c$., Comment: Accepted for publication in the Astrophysical Journal Letters (ApJL)

We present the first 2.5 years of data from the MeerKAT Pulsar Timing Array (MPTA), part of MeerTime, a MeerKAT Large Survey Project. The MPTA aims to precisely measure pulse arrival times from an ensemble of 88 pulsars visible from the Southern Hemisphere, with the goal of contributing to the search, detection and study of nanohertz-frequency gravitational waves as part of the International Pulsar Timing Array. This project makes use of the MeerKAT telescope, and operates with a typical observing cadence of two weeks using the L-band receiver that records data from 856-1712 MHz. We provide a comprehensive description of the observing system, software, and pipelines used and developed for the MeerTime project. The data products made available as part of this data release are from the 78 pulsars that had at least $30$ observations between the start of the MeerTime programme in February 2019 and October 2021. These include both sub-banded and band-averaged arrival times, as well as the initial timing ephemerides, noise models, and the frequency-dependent standard templates (portraits) used to derive pulse arrival times. After accounting for detected noise processes in the data, the frequency-averaged residuals of $67$ of the pulsars achieved a root-mean-square residual precision of $< 1 \mu \rm{s}$. We also present a novel recovery of the clock correction waveform solely from pulsar timing residuals, and an exploration into preliminary findings of interest to the international pulsar timing community. The arrival times, standards and full Stokes parameter calibrated pulsar timing archives are publicly available., Comment: 18 pages, 6 figures

Precision timing of millisecond pulsars in binary systems enables observers to detect the relativistic Shapiro delay induced by space time curvature. When favourably aligned, this enables constraints to be placed on the component masses and system orientation. Here we present the results of timing campaigns on seven binary millisecond pulsars observed with the 64-antenna MeerKAT radio telescope that show evidence of Shapiro delay: PSRs~J0101$-$6422, J1101$-$6424, J1125$-$6014, J1514$-$4946, J1614$-$2230, J1732$-$5049, and J1909$-$3744. Evidence for Shapiro delay was found in all of the systems, and for three the orientations and data quality enabled strong constraints on their orbital inclinations and component masses. For PSRs~J1125$-$6014, J1614$-$2230 and J1909$-$3744, we determined pulsar masses to be $M_{\rm p} = 1.68\pm 0.17 \, {\rm M_{\odot}} $, $1.94\pm 0.03 \, {\rm M_{\odot}} $ and $1.45 \pm 0.03 \, {\rm M_{\odot}}$, and companion masses to be $M_{\rm c} = 0.33\pm 0.02 \, {\rm M_{\odot}} $, $0.495\pm 0.005 \, {\rm M_{\odot}} $ and $0.205 \pm 0.003 \, {\rm M_{\odot}}$, respectively. This provides the first independent confirmation of PSR~J1614$-$2230's mass, one of the highest known. The Shapiro delays measured for PSRs~J0101$-$6422, J1101$-$6424, J1514$-$4946, and J1732$-$5049 were only weak, and could not provide interesting component mass limits. Despite a large number of millisecond pulsars being routinely timed, relatively few have accurate masses via Shapiro delays. We use simulations to show that this is expected, and provide a formula for observers to assess how accurately a pulsar mass can be determined. We also discuss the observed correlation between pulsar companion masses and spin period, and the anti-correlation between recycled pulsar mass and their companion masses., Comment: 19 pages, 8 figures

In spite of the rich phenomenology of the polarization properties of radio pulsars, the rotating vector model (RVM) created 50 years ago remains the best method to determine the beam geometry of a pulsar. We apply the RVM to a sample of 854 radio pulsars observed with the MeerKAT telescope in order to draw conclusions about the population of pulsars as a whole. The main results are that (i) the geometrical interpretation of the position angle traverse is valid in the majority of the population, (ii) the pulsars for which the RVM fails tend to have a high fraction of circular polarization compared to linear polarization, (iii) emission heights obtained through both geometrical and relativistic methods show that the majority of pulsars must have emission heights less than 1000~km independent of spin period, (iv) orthogonal mode jumps are seen in the position angle traverse in about one third of the population. All these results are weakly dependent on the pulsar spin-down energy., Comment: accepted by MNRAS

Magnetars are conjectured to be highly magnetized neutron stars (NSs). Strong internal magnetic field and elasticity in the crust may deform the stars and lead to free precession. We study the precession dynamics of triaxially-deformed NSs incorporating the near-field and the far-field electromagnetic torques. We obtain timing residuals for different NS geometries and torques. We also investigate the polarized X-ray and radio signals from precessing magnetars. The modulations on the Stokes parameters are obtained for thermal X-rays emitted from the surface of magnetars. For radio signals, we apply the simple rotating vector model (RVM) to give the modulations on the position angle (PA) of the polarization. Our results are comprehensive, ready to be used to search for magnetar precession with timing data and polarizations of X-ray and radio emissions. Future observations of precessing magnetars will give us valuable information on the geometry and the strength of the strong magnetic fields, the emission geometry, as well as the equation of state (EoS) of NSs., Comment: 18 pages, 18 figures; accepted by MNRAS

This White Paper summarises potential key science topics to be achieved with Thai National Radio Telescope (TNRT). The commissioning phase has started in mid 2022. The key science topics consist of "Pulsars and Fast Radio Bursts (FRBs)", "Star Forming Regions (SFRs)", "Galaxy and Active Galactic Nuclei (AGNs)", "Evolved Stars", "Radio Emission of Chemically Peculiar (CP) Stars", and "Geodesy", covering a wide range of observing frequencies in L/C/X/Ku/K/Q/W-bands (1-115 GHz). As a single-dish instrument, TNRT is a perfect tool to explore time domain astronomy with its agile observing systems and flexible operation. Due to its ideal geographical location, TNRT will significantly enhance Very Long Baseline Interferometry (VLBI) arrays, such as East Asian VLBI Network (EAVN), Australia Long Baseline Array (LBA), European VLBI Network (EVN), in particular via providing a unique coverage of the sky resulting in a better complete "uv" coverage, improving synthesized-beam and imaging quality with reducing side-lobes. This document highlights key science topics achievable with TNRT in single-dish mode and in collaboration with VLBI arrays., Comment: 75 pages, 16 figures, 15 tables. White Paper for Potential Key Sciences to be Achieved with Thai National Radio Telescope (TNRT)

High sensitivity interstellar scintillation and polarization observations of PSR~B0656+14 made at three epochs over a year using the Five-hundred-meter Aperture Spherical radio Telescope (FAST) show that the scattering is dominated by two different compact regions. We identify the one nearer to the pulsar with the shell of the Monogem Ring, thereby confirming the association. The other is probably associated with the Local Bubble. We find that the observed position angles of the pulsar spin axis and the space velocity are significantly different, with a separation of $19\fdg3\pm$0\fdg8, inconsistent with a previously published near-perfect alignment of $1\degr\pm 2\degr$. The two independent scattering regions are clearly defined in the secondary spectra which show two strong forward parabolic arcs. The arc curvatures imply that the scattering screens corresponding to the outer and inner arcs are located approximately 28~pc from PSR B0656+14 and 185~pc from the Earth, respectively. Comparison of the observed Doppler profiles with electromagnetic simulations shows that both scattering regions are mildly anisotropic. For the outer arc, we estimate the anisotropy $A_R$ to be approximately 1.3, with the scattering irregularities aligned parallel to the pulsar velocity. For the outer arc, we compare the observed delay profiles with delay profiles computed from a theoretical strong-scattering model. Our results suggest that the spatial spectrum of the scattering irregularities in the Monogem Ring is flatter than Kolmogorov, but further observations are required to confirm this., Comment: 21 pages, 11 figures, accepted for publication in the Astronomical Journal

The large number of ongoing surveys for pulsars and transients at various radio observatories is motivated by the science obtained from these sources. Timing and polarisation analysis of relativistic binaries can place strong constraints on theories of gravity. The observation of a growing number of millisecond pulsars (MSPs) spread over the celestial sphere may allow the detection of a stochastic gravitational wave background arising from supermassive black hole binaries. A more complete sample of young pulsars improves our knowledge of neutron star birth and evolution. Transients such as fast radio bursts can serve to probe the intergalactic medium. The SPAN512 pulsar survey covers intermediate Galactic latitudes using the L-band receiver of the Nan\c{c}ay Radio Telescope (NRT). The survey covers 224 sq. deg. of the sky for a total exposure time of 2200 h. Population syntheses predict the discovery of 3 to 19 new normal pulsars and a few MSPs. We present detailed modelling of the NRT beam with its L-band receiver and its sensitivity which we used to precisely assess the expected survey yield. We used the flexible Pulsar Arecibo L-band Feed Array data processing pipeline to search the 47 TB of SPAN512 data for pulsars and transients. The SPAN512 survey discovered two new MSPs and one new middle-aged pulsar. We focus on the analysis of the 2.4-ms spin period pulsar J2205+6012 for which we also report the detection of gamma-ray pulsations. Its narrow pulse width (35 $\mu$s at an observing frequency of 2.55 GHz) allows for sub-microsecond timing precision over 8 years, with exciting prospects for pulsar timing array programs., Comment: 13 pages, 11 figures, 5 tables. Accepted in A&A