Search

Radio pulsars allow the study of the ionised interstellar medium and its dispersive effects, a major noise source in gravitational wave searches using pulsars. In this paper, we compare the functionality and reliability of three commonly used schemes to measure temporal variations in interstellar propagation effects in pulsar-timing data. We carry out extensive simulations at low observing frequencies (100-200 MHz) by injecting long-term correlated noise processes with power-law spectra and white noise, to evaluate the robustness, accuracy and precision of the following three mitigation methods: epoch-wise (EW) measurements of interstellar dispersion; the DMX method of simultaneous, piece-wise fits to interstellar dispersion; and DMGP, which models dispersion variations through Gaussian processes using a Bayesian analysis method. We then evaluate how reliably the input signals are reconstructed and how the various methods react to the presence of achromatic long-period noise. All the methods perform well, provided the achromatic long-period noise is modeled for DMX and DMGP. The most precise method is DMGP, followed by DMX and EW, while the most accurate is EW, followed by DMX and DMGP. We also test different scenarios including simulations of L-band ToAs and realistic DM injection, with no significant variation in the obtained results. Given the nature of our simulations and our scope, we deem that EW is the most reliable method to study the Galactic ionized media. Future works should be conducted to confirm this result via more realistic simulations. We note that DM GP and DMX seem to be the most performing techniques in removing long-term correlated noise, and hence for gravitational wave studies. However, full simulations of pulsar timing array experiments are needed to support this interpretation.

Synchrotron X-ray emission has been detected from nearly a dozen young supernova remnants (SNRs). X-rays of synchrotron origin exhibit linear polarization in a regular, non-randomly oriented magnetic field. The significant polarized X-ray emission from four such SNRs has already been reported on the basis of observations with the Imaging X-ray Polarimetry Explorer (IXPE). The magnetic-field structure as derived from IXPE observations is radial for Cassiopeia A, Tycho's SNR, and SN 1006, and tangential for RX J1713.7-3946. The latter together with the recent detection of a tangential magnetic field in SNR 1E 0102.2-7219 by the Australia Telescope Compact Array in the radio band shows that tangential magnetic fields can also be present in young SNRs. Thus, the dichotomy in polarization between young and middle-aged SNRs (radial magnetic fields in young SNRs, but tangential magnetic fields in middle-aged SNRs), previously noticed in the radio band, deserves additional attention. The present analysis of IXPE observations determines, for the first time, a magnetic-field structure in the northwestern rim of Vela Jr, also known as RX J0852.0-4622, and provides a new example of a young SNR with a tangential magnetic field., Comment: Accepted for publication in Astronomy and Astrophysics

X-ray polarization is a unique new probe of the particle acceleration in astrophysical jets made possible through the Imaging X-ray Polarimetry Explorer. Here we report on the first dense X-ray polarization monitoring campaign on the blazar Mrk 421. Our observations were accompanied by an even denser radio and optical polarization campaign. We find significant short-timescale variability in both X-ray polarization degree and angle, including a $\sim90^\circ$ angle rotation about the jet axis. We attribute this to random variations of the magnetic field, consistent with the presence of turbulence but also unlikely to be explained by turbulence alone. At the same time, the degree of lower-energy polarization is significantly lower and shows no more than mild variability. Our campaign provides further evidence for a scenario in which energy-stratified shock-acceleration of relativistic electrons, combined with a turbulent magnetic field, is responsible for optical to X-ray synchrotron emission in blazar jets., Comment: 23 pages, including 8 pages of appendices. 12 figures, 3 tables. Submitted to ApJ

The Small and Large Magellanic Clouds are the only galaxies outside our own in which radio pulsars have been discovered to date. The sensitivity of the MeerKAT radio interferometer offers an opportunity to search for a population of more distant extragalactic pulsars. The TRAPUM (TRansients And PUlsars with MeerKAT) collaboration has performed a radio-domain search for pulsars and transients in the dwarf star-forming galaxies Sextans A and B, situated at the edge of the local group 1.4 Mpc away. We conducted three 2-hour multi-beam observations at L-band (856-1712 MHz) with the full array of MeerKAT. No pulsars were found down to a radio pseudo-luminosity upper limit of 7.9$\pm$0.4 Jy kpc$^{2}$ at 1400 MHz, which is 28 times more sensitive than the previous limit from the Murriyang telescope. This luminosity is 30 per cent greater than that of the brightest known radio pulsar and sets a cut-off on the luminosity distributions of the entire Sextans A and B galaxies for unobscured radio pulsars beamed in our direction. A Fast Radio Burst was detected in one of the Sextans A observations at a Dispersion Measure (DM) of 737 pc cm$^{-3}$. We believe this is a background event not associated with the dwarf galaxy due to its large DM and its S/N being strongest in the wide-field incoherent beam of MeerKAT., Comment: 11 pages, 9 figures, 5 tables. Accepted for publication in Monthly Notices of the Royal Astronomical Society

PSR 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., Comment: 17 pages, 11 figures, 3 tables

Current observational evidence reveals that fast radio bursts (FRBs) exhibit bandwidths ranging from a few dozen MHz to several GHz. Traditional FRB searches primarily employ matched filter methods on time series collapsed across the entire observational bandwidth. However, with modern ultra-wideband receivers featuring GHz-scale observational bandwidths, this approach may overlook a significant number of events. We investigate the efficacy of sub-banded searches for FRBs, a technique seeking bursts within limited portions of the bandwidth. These searches aim to enhance the significance of FRB detections by mitigating the impact of noise outside the targeted frequency range, thereby improving signal-to-noise ratios. We conducted a series of Monte Carlo simulations, for the $400$-MHz bandwidth Parkes 21-cm multi-beam (PMB) receiver system and the Parkes Ultra-Wideband Low (UWL) receiver, simulating bursts down to frequency widths of about $100$\,MHz. Additionally, we performed a complete reprocessing of the high-latitude segment of the High Time Resolution Universe South survey (HTRU-S) of the Parkes-Murriyang telescope using sub-banded search techniques. Simulations reveal that a sub-banded search can enhance the burst search efficiency by $67_{-42}^{+133}$ % for the PMB system and $1433_{-126}^{+143}$ % for the UWL receiver. Furthermore, the reprocessing of HTRU led to the confident detection of eighteen new bursts, nearly tripling the count of FRBs found in this survey. These results underscore the importance of employing sub-banded search methodologies to effectively address the often modest spectral occupancy of these signals., Comment: Accepted for publication (A&A)

The Large Magellanic Cloud (LMC) presents a unique environment for pulsar population studies due to its distinct star formation characteristics and proximity to the Milky Way. As part of the TRAPUM (TRAnsients and PUlsars with MeerKAT) Large Survey Project, we are using the core array of the MeerKAT radio telescope (MeerKAT) to conduct a targeted search of the LMC for radio pulsars at L-band frequencies, 856-1712$\,$MHz. The excellent sensitivity of MeerKAT, coupled with a 2-hour integration time, makes the survey 3 times more sensitive than previous LMC radio pulsar surveys. We report the results from the initial four survey pointings which has resulted in the discovery of seven new radio pulsars, increasing the LMC radio pulsar population by 30 per cent. The pulse periods of these new pulsars range from 278 to 1690$\,$ms, and the highest dispersion measure is 254.20$\,$pc$\,$cm$^{-3}$. We searched for, but did not find any significant pulsed radio emission in a beam centred on the SN$\,$1987A remnant, establishing an upper limit of 6.3$\,{\mu}$Jy on its minimum flux density at 1400$\,$MHz., Comment: 12 pages, 4 figures, 4 tables. Accepted for publication in Monthly Notices of the Royal Astronomical Society

We report new radio timing solutions from a three-year observing campaign conducted with the MeerKAT and Murriyang telescopes for nine Small Magellanic Cloud pulsars, increasing the number of characterised rotation-powered extragalactic pulsars by 40 per cent. We can infer from our determined parameters that the pulsars are seemingly all isolated, that six are ordinary pulsars, and that three of the recent MeerKAT discoveries have a young characteristic age of under 100 kyr and have undergone a spin-up glitch. Two of the sources, PSRs J0040$-$7337 and J0048$-$7317, are energetic young pulsars with spin-down luminosities of the order of 10$^{36}$ erg s$^{-1}$. They both experienced a large glitch, with a change in frequency of about 30 $\mu$Hz, and a frequency derivative change of order $-10^{-14}$ Hz s$^{-1}$. These glitches, the inferred glitch rate, and the properties of these pulsars (including potentially high inter-glitch braking indices) suggest these neutron stars might be Vela-like repeating glitchers and should be closely monitored in the future. The position and energetics of PSR J0048$-$7317 confirm it is powering a new Pulsar Wind Nebula (PWN) detected as a radio continuum source; and similarly the association of PSR J0040$-$7337 with the PWN of Supernova Remnant (SNR) DEM S5 (for which we present a new Chandra image) is strengthened. Finally, PSR J0040$-$7335 is also contained within the same SNR but is a chance superposition. It has also been seen to glitch with a change of frequency of $10^{-2}$ $\mu$Hz. This work more than doubles the characterised population of SMC radio pulsars., Comment: 20 pages, 13 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society

The Sculptor Galaxy (NGC 253), located in the Southern Hemisphere, far off the Galactic Plane, has a relatively high star-formation rate of about 7 M$_{\odot}$ yr$^{-1}$ and hosts a young and bright stellar population, including several super star clusters and supernova remnants. It is also the first galaxy, apart from the Milky Way Galaxy to be associated with two giant magnetar flares. As such, it is a potential host of pulsars and/or fast radio bursts in the nearby Universe. The instantaneous sensitivity and multibeam sky coverage offered by MeerKAT therefore make it a favourable target. We searched for pulsars, radio-emitting magnetars, and fast radio bursts in NGC 253 as part of the TRAPUM large survey project with MeerKAT. We did not find any pulsars during a four-hour observation, and derive a flux density limit of 4.4 $\mu$Jy at 1400 MHz, limiting the pseudo-luminosity of the brightest putative pulsar in this galaxy to 54 Jy kpc$^2$. Assuming universality of pulsar populations between galaxies, we estimate that detecting a pulsar as bright as this limit requires NGC 253 to contain a pulsar population of $\gtrsim$20 000. We also did not detect any single pulses and our single pulse search flux density limit is 62 mJy at 1284 MHz. Our search is sensitive enough to have detected any fast radio bursts and radio emission similar to the brighter pulses seen from the magnetar SGR J1935+2154 if they had occurred during our observation., Comment: 6 pages, 3 figures, added citation in section 3, corrected values in section 4, published in MNRAS

The eclipses seen in the radio emission of some pulsars can be invaluable to study the properties of the material from the companion stripped away by the pulsar. We present a study of six consecutive eclipses of PSR J0024-7204O in the globular cluster 47 Tucanae as seen by the MeerKAT radio telescope in the UHF (544-1088 MHz) band. A high scintillation state boosted the signal during one of the orbits and allowed a detailed study of the eclipse properties. We measure significant dispersion measure (DM) variations and detect strong scattering that seems to be the dominating mechanism of the eclipses at these frequencies. A complete drop in the linear polarization together with a small increase in the rotation measure suggests the presence of a magnetic field of $\sim 2$ mG. The study of multiple eclipses allowed us to measure difference in the lengths of the eclipses and DM differences of $\sim 0.01$ pc cm$^{-3}$ in consecutive orbits. One orbit in particular shows a delay in recovery of the linear polarization and a visible delay in the arrival of the pulses caused by a stronger scattering event. We suggest that these are caused by a higher variance of density fluctuations during the event., Comment: 14 pages, 8 figures, 9 tables. Accepted for publication on Monthly Notices of the Royal Astronomical Society

We present Crab X-ray polarization measurements using IXPE data with a total exposure of 300ks, three times more than the initial 2022 discovery paper. Polarization is detected in three times more pulsar phase bins, revealing an S-shaped $+40^\circ$ polarization angle sweep in the main pulse and ${>}1\sigma$ departures from the OPTIMA optical polarization in both pulses, suggesting different radiation mechanisms or sites for the polarized emission at the two wavebands. Our polarization map of the inner nebula reveals a toroidal magnetic field, as seen in prior IXPE analyses. Along the southern jet, the magnetic field orientation relative to the jet axis changes from perpendicular to parallel and the polarization degree decreases by ${\sim}6\%$. These observations may be explained by kink instabilities along the jet or a collision with a dense, jet-deflecting medium at the tip. Using spectropolarimetric analysis, we find asymmetric polarization in the four quadrants of the inner nebula, as expected for a toroidal field geometry, and a spatial correlation between polarization degree and photon index., Comment: 14 pages, 7 figures, 3 tables. Accepted by The Astrophysical Journal

We present multiwavelength polarization measurements of the luminous blazar Mrk~501 over a 14-month period. The 2--8 keV X-ray polarization was measured with the Imaging X-ray Polarimetry Explorer (IXPE) with six 100-ks observations spanning from 2022 March to 2023 April. Each IXPE observation was accompanied by simultaneous X-ray data from NuSTAR, Swift/XRT, and/or XMM-Newton. Complementary optical-infrared polarization measurements were also available in the B, V, R, I, and J bands, as were radio polarization measurements from 4.85 GHz to 225.5 GHz. Among the first five IXPE observations, we did not find significant variability in the X-ray polarization degree and angle with IXPE. However, the most recent sixth observation found an elevated polarization degree at $>3\sigma$ above the average of the other five observations. The optical and radio measurements show no apparent correlations with the X-ray polarization properties. Throughout the six IXPE observations, the X-ray polarization degree remained higher than, or similar to, the R-band optical polarization degree, which remained higher than the radio value. This is consistent with the energy-stratified shock scenario proposed to explain the first two IXPE observations, in which the polarized X-ray, optical, and radio emission arises from different regions., Comment: Accepted for publication in ApJ

Fast Radio Bursts (FRBs) are cosmological sub-second bursts of coherent radio emission, whose source is still unknown. To date, the galactic magnetar SGR 1935+2154 is the only astrophysical object known to emit radio bursts akin to FRBs, albeit less powerful, supporting suggestions that FRBs originate from magnetars. Many remarkable properties of FRBs, e.g. the dichotomy between repeaters and one-off sources, and their power-law energy distributions (with typical index $\sim 2-3$), are not well understood yet. Moreover, the huge radio power released by the most active repeaters is challenging even for the magnetic energy reservoir of magnetars. Here we assume that FRBs originate from co-rotating hot-spots anchored in neutron star magnetospheres and get occasionally amplified by large factors via gravitational self-lensing in the strong NS field. We evaluate the probability of amplification and show that (i) a power-law energy distribution of events $\propto E^{-(2- 3)}$ is generally expected, (ii) all FRB sources may be regarded as repeating, their appearance as one-off sources or repeaters being determined by the critical dependence of the amplification probability on the emission geometry and source orientation relative to Earth and (iii) the most active repeaters, in particular, correspond to extremely rare and finely-tuned orientations ($\sim$ one in $10^6$), leading to large probabilities of amplification which make their bursts frequently detectable. At the same time, their power release appears enhanced, typically by factors $\gtrsim 10$, easing their energy budget problem., Comment: 23 pages (17 main text + Appendices), 13 Figures, accepted for publication on The Astrophysical Journal

Millisecond pulsars are subject to accelerations in globular clusters (GCs) that manifest themselves in both the first and second spin period time derivatives, and can be used to explore the mass distribution of the potentials they inhabit. Here we report on over 20 yr of pulsar timing observations of five millisecond radio pulsars in the core of the core-collapse GC NGC 6752 with the Parkes (Murriyang) and MeerKAT radio telescopes, which have allowed us to measure the proper motions, positions, and first and second time derivatives of the pulsars. The pulsar timing parameters indicate that all the pulsars in the core experience accelerations and jerks that can be explained only if an amount of nonluminous mass of at least 2.56x10^3 M_SUN is present in the core of NGC 6752. On the other hand, our studies highly disfavor the presence of an intermediate-mass black hole at the center of the cluster, with a mass equal to or greater than ~3000M_SUN., Comment: Published in ApJ. Replaced with the final published version

We present the first X-ray spectropolarimetric results for Cygnus X-1 in its soft state from a campaign of five IXPE observations conducted during 2023 May-June. Companion multiwavelength data during the campaign are likewise shown. The 2-8 keV X-rays exhibit a net polarization degree PD=1.99%+/-0.13% (68% confidence). The polarization signal is found to increase with energy across IXPE's 2-8 keV bandpass. The polarized X-rays exhibit an energy-independent polarization angle of PA=-25.7+/-1.8 deg. East of North (68% confidence). This is consistent with being aligned to Cyg X-1's AU-scale compact radio jet and its pc-scale radio lobes. In comparison to earlier hard-state observations, the soft state exhibits a factor of 2 lower polarization degree, but a similar trend with energy and a similar (also energy-independent) position angle. When scaling by the natural unit of the disk temperature, we find the appearance of a consistent trendline in the polarization degree between soft and hard states. Our favored polarimetric model indicates Cyg X-1's spin is likely high (a* above ~0.96). The substantial X-ray polarization in Cyg X-1's soft state is most readily explained as resulting from a large portion of X-rays emitted from the disk returning and reflecting off the disk surface, generating a high polarization degree and a polarization direction parallel to the black hole spin axis and radio jet. In IXPE's bandpass, the polarization signal is dominated by the returning reflection emission. This constitutes polarimetric evidence for strong gravitational lensing of X-rays close to the black hole., Comment: 20 pages, accepted for publication in ApJL

Recent observations of X-ray pulsars (XRPs) performed by the Imaging X-ray Polarimetry Explorer (IXPE) have made it possible to investigate the intricate details of these objects in a new way, thanks to the added value of X-ray polarimetry. Here we present the results of the IXPE observations of SMC X-1, a member of the small group of XRPs displaying super-orbital variability. SMC X-1 was observed by IXPE three separate times during the high state of its super-orbital period. The observed luminosity in the 2-8 keV energy band of $L=2\times10^{38}$ erg/s makes SMC X-1 the brightest XRP ever observed by IXPE. We detect significant polarization in all three observations, with values of the phase-averaged polarization degree (PD) and polarization angle (PA) of $3.2\pm0.8$% and $97\deg\pm8\deg$ for Observation 1, $3.0\pm0.9$% and $90\deg\pm8\deg$ for Observation 2, and $5.5\pm1.1$% and $80\deg\pm6\deg$ for Observation 3, for the spectro-polarimetric analysis. The observed PD shows an increase over time with decreasing luminosity, while the PA decreases in decrements of 10\deg. The phase-resolved spectro-polarimetric analysis reveals significant detection of polarization in three out of seven phase bins, with the PD ranging between 2% and 10%, and a corresponding range in the PA from $\sim$70\deg\ to $\sim$100\deg. The pulse-phase resolved PD displays an apparent anti-correlation with the flux. Using the rotating vector model, we obtain constraints on the pulsar's geometrical properties for the individual observations. The position angle of the pulsar displays an evolution over time supporting the idea that we observe changes related to different super-orbital phases. Scattering in the wind of the precessing accretion disk may be responsible for the behavior of the polarimetric properties observed during the high-state of SMC X-1's super-orbital period., Comment: 11 pages, 11 figures, submitted to A&A

We report the X-ray polarization properties of the high-synchrotron-peaked (HSP) blazar PKS 2155$-$304 based on observations with the Imaging X-ray Polarimetry Explorer (IXPE). We observed the source between Oct 27 and Nov 7, 2023. We also conducted an extensive contemporaneous multiwavelength (MW) campaign. We find that during the first half ($T_1$) of the IXPE pointing, the source exhibited the highest X-ray polarization degree detected for an HSP blazar thus far, (30.7$\pm$2.0)%, which dropped to (15.3$\pm$2.1)% during the second half ($T_2$). The X-ray polarization angle remained stable during the IXPE pointing at 129.4$^\circ$$\pm$1.8$^\circ$ and 125.4$^\circ$$\pm$3.9$^\circ$ during $T_1$ and $T_2$, respectively. Meanwhile, the optical polarization degree remained stable during the IXPE pointing, with average host-galaxy-corrected values of (4.3$\pm$0.7)% and (3.8$\pm$0.9)% during the $T_1$ and $T_2$, respectively. During the IXPE pointing, the optical polarization angle changed achromatically from $\sim$140$^\circ$ to $\sim$90$^\circ$ and back to $\sim$130$^\circ$. Despite several attempts, we only detected (99.7% conf.) the radio polarization once (during $T_2$, at 225.5 GHz): with degree (1.7$\pm$0.4)% and angle 112.5$^\circ$$\pm$5.5$^\circ$. The direction of the broad pc-scale jet is rather ambiguous and has been found to point to the east and south at different epochs; however, on larger scales (> 1.5 pc) the jet points toward the southeast ($\sim$135$^\circ$), similar to all of the MW polarization angles. Moreover, the X-ray to optical polarization degree ratios of $\sim$7 and $\sim$4 during $T_1$ and $T_2$, respectively, are similar to previous IXPE results for several HSP blazars. These findings, combined with the lack of correlation of temporal variability between the MW polarization properties, agree with an energy-stratified shock-acceleration scenario in HSP blazars., Comment: 17 pages, 10 figures, 4 tables, Accepted for publication in A&A

Pulsar timing arrays (PTAs) are designed to detect low-frequency gravitational waves (GWs). GWs induce achromatic signals in PTA data, meaning that the timing delays do not depend on radio-frequency. However, pulse arrival times are also affected by radio-frequency dependent "chromatic" noise from sources such as dispersion measure (DM) and scattering delay variations. Furthermore, the characterization of GW signals may be influenced by the choice of chromatic noise model for each pulsar. To better understand this effect, we assess if and how different chromatic noise models affect achromatic noise properties in each pulsar. The models we compare include existing DM models used by NANOGrav and noise models used for the European PTA Data Release 2 (EPTA DR2). We perform this comparison using a subsample of six pulsars from the NANOGrav 15 yr data set, selecting the same six pulsars as from the EPTA DR2 six-pulsar dataset. We find that the choice of chromatic noise model noticeably affects the achromatic noise properties of several pulsars. This is most dramatic for PSR J1713+0747, where the amplitude of its achromatic red noise lowers from $\log_{10}A_{\text{RN}} = -14.1^{+0.1}_{-0.1}$ to $-14.7^{+0.3}_{-0.5}$, and the spectral index broadens from $\gamma_{\text{RN}} = 2.6^{+0.5}_{-0.4}$ to $\gamma_{\text{RN}} = 3.5^{+1.2}_{-0.9}$. We also compare each pulsar's noise properties with those inferred from the EPTA DR2, using the same models. From the discrepancies, we identify potential areas where the noise models could be improved. These results highlight the potential for custom chromatic noise models to improve PTA sensitivity to GWs.

The sensitivity of the MeerKAT radio interferometer is an opportunity to probe deeper into the population of rare and faint extragalactic pulsars. The TRAPUM (TRAnsients and PUlsars with MeerKAT) collaboration has conducted a radio-domain search for accelerated pulsars and transients in the Small Magellanic Cloud (SMC). This partially targeted survey, performed at L-band (856-1712 MHz) with the core array of the MeerKAT telescope in 2-h integrations, is twice as sensitive as the latest SMC radio pulsar survey. We report the discovery of seven new SMC pulsars, doubling this galaxy's radio pulsar population and increasing the total extragalactic population by nearly a quarter. We also carried out a search for accelerated millisecond pulsars in the SMC Globular Cluster NGC 121 using the full array of MeerKAT. This improved the previous upper limit on pulsed radio emission from this cluster by a factor of six. Our discoveries reveal the first radio pulsar-PWN systems in the SMC, with only one such system previously known outside our galaxy (the "Crab pulsar twin" in the Large Magellanic Cloud, PSR J0540$-$6919). We associate the 59 ms pulsar discovery PSR J0040$-$7337, now the fastest spinning radio pulsar in the SMC, with the bow-shock Pulsar Wind Nebula (PWN) of Supernova Remnant DEM S5. We also present a new young pulsar with a 79 ms period, PSR J0048$-$7317, in a PWN recently discovered in a MeerKAT radio continuum image. Using the multi-beam capability of MeerKAT, we localised our pulsar discoveries, and two previous Murriyang discoveries, to a positional uncertainty of a few arcseconds., Comment: 32 pages, 14 figures, 10 tables. Accepted for publication in Monthly Notices of the Royal Astronomical Society

The measurement of pulsar pulse times-of-arrival (ToAs) is a crucial step in detecting low-frequency gravitational waves. To determine ToAs, we can use template-matching to compare each observed pulse profile with a standard template. However, using different combinations of templates and template-matching methods (TMMs) without careful consideration may lead to inconsistent results. In pulsar timing array (PTA) experiments, distinct ToAs from the same observations can be obtained, due to the use of diverse templates and TMMs. In other words, employing diverse approaches can yield different timing results and would thus have a significant impact on subsequent gravitational wave searches. In this paper, we examine several commonly used combinations to analyze their effect on pulse ToAs. we evaluate the potential impact of template and TMM selection on thirteen typical millisecond pulsars within the European PTA. We employ pulsar timing methods, specifically the root mean square and reduced chi-square $\chi_r^2$ of the residuals of the best timing solution to assess the outcomes. Additionally, we evaluate the system-limited noise floor (SLNF) for each pulsar at various telescopes operating around 1.4~GHz using frequency-resolved templates. Our findings suggest that utilizing data-derived and smoothed templates in conjunction with the Fourier-domain with Markov-chain Monte Carlo (FDM) TMM is generally the most effective approach, though there may be exceptions that require further attention. Furthermore, we determine that pulse phase jitter noise does not significantly limit the current precision of the European PTA's timing, as jitter levels derived from other studies are much smaller than the SLNF., Comment: Accepted by A&A

Discovery of pulsations from a number of ultra-luminous X-ray (ULX) sources proved that accretion onto neutron stars can produce luminosities exceeding the Eddington limit by a couple of orders of magnitude. The conditions necessary to achieve such high luminosities as well as the exact geometry of the accretion flow in the neutron star vicinity are, however, a matter of debate. The pulse phase-resolved polarization measurements that became possible with the launch of the IXPE can be used to determine the pulsar geometry and its orientation relative to the orbital plane. They provide an avenue to test different theoretical models of ULX pulsars. In this paper we present the results of three IXPE observations of the first Galactic ULX pulsar Swift J0243.6+6124 during its 2023 outburst. We find strong variations of the polarization characteristics with the pulsar phase. The average polarization degree increases from about 5% to 15% as the flux dropped by a factor of three in the course of the outburst. The polarization angle (PA) as function of the pulsar phase shows two peaks in the first two observations, but changes to a characteristic sawtooth pattern in the remaining data set. This is not consistent with a simple rotating vector model. Assuming the existence of an additional constant polarized component, we were able to fit the three observations with a common rotating vector model and obtain constraints on the pulsar geometry. In particular, we find the pulsar angular momentum inclination with respect to the line-of-sight of 15-40 deg, the magnetic obliquity of 60-80 deg, and the pulsar spin position angle of -50 deg, which differs from the constant component PA of about 10 deg. Combining these X-ray measurements with the optical PA, we find evidence for a 30 deg misalignment between the pulsar spin and the binary orbital axis., Comment: 13 pages, 10 figures, submitted to A&A

Supernova remnants (SNRs) provide insights into cosmic-ray acceleration and magnetic field dynamics at shock fronts. Recent X-ray polarimetric measurements by the Imaging X-ray Polarimetry Explorer (IXPE) have revealed radial magnetic fields near particle acceleration sites in young SNRs, including Cassiopeia A, Tycho, and SN 1006. We present here the spatially-resolved IXPE X-ray polarimetric observation of the northwestern rim of SNR RX J1713.7-3946. For the first time, our analysis shows that the magnetic field in particle acceleration sites of this SNR is oriented tangentially with respect to the shock front. Because of the lack of precise Faraday-rotation measurements in the radio band, this was not possible before. The average measured polarization degree (PD) of the synchtrotron emission is 12.5 {\pm} 3.3%, lower than the one measured by IXPE in SN 1006, comparable to the Tycho one, but notably higher than the one in Cassiopeia A. On sub-parsec scales, localized patches within RX J1713.7-3946 display PD up to 41.5 {\pm} 9.5%. These results are compatible with a shock-compressed magnetic field. However, in order to explain the observed PD, either the presence of a radial net magnetic field upstream of the shock, or partial reisotropization of the turbulence downstream by radial magneto-hydrodynamical instabilities, can be invoked. From comparison of PD and magnetic field distribution with {\gamma}-rays and 12 CO data, our results provide new inputs in favor of a leptonic origin of the {\gamma}-ray emission., Comment: 18 pages, 6 figures, 2 tables, published in ApJ Letters

We report the discovery of ten new pulsars in the globular cluster Terzan 5 as part of the Transients and Pulsars with MeerKAT (TRAPUM) Large Survey Project. We observed Terzan 5 at L-band (856--1712 MHz) with the MeerKAT radio telescope for four hours on two epochs, and performed acceleration searches of 45 out of 288 tied-array beams covering the core of the cluster. We obtained phase-connected timing solutions for nine discoveries, covering nearly two decades of archival observations from the Green Bank Telescope for all but one. Highlights include PSR J1748$-$2446ao which is an eccentric ($e = 0.32$) wide-orbit (orbital period $P_{\rm b} = 57.55$ d) system. We were able to measure the rate of advance of periastron ($\dot{\omega}$) for this system allowing us to determine a total mass of $3.17 \pm \, 0.02\, \rm M_{\odot}$. With a minimum companion mass ($M_{\rm c}$) of $\sim 0.8\, \rm M_{\odot}$, PSR J1748$-$2446ao is a candidate double neutron star (DNS) system. If confirmed to be a DNS, it would be the fastest spinning pulsar ($P = 2.27$ ms) and the longest orbital period measured for any known DNS system. PSR J1748$-$2446ap has the second highest eccentricity for any recycled pulsar ($e \sim 0.905$) and for this system we can measure the total mass ($1.997 \pm 0.006\, \rm M_{\odot}$) and also estimate the individual pulsar and companion masses. PSR J1748$-$2446ar is an eclipsing redback (minimum $M_{\rm c} \sim 0.34\, \rm M_{\odot}$) system whose properties confirm it to be the counterpart to a previously published source identified in radio and X-ray imaging. With these discoveries, the total number of confirmed pulsars in Terzan 5 is 49, the highest for any globular cluster so far. These discoveries further enhance the rich set of pulsars known in Terzan 5 and provide scope for a deeper understanding of binary stellar evolution, cluster dynamics and ensemble population studies., Comment: 23 pages, 11 figures, 5 tables, published in A&A

Using MeerKAT, we have discovered three new millisecond pulsars (MSPs) in the bulge globular cluster M62: M62H, M62I, and M62J. All three are in binary systems, which means all ten known pulsars in the cluster are in binaries. M62H has a planetary-mass companion with a median mass $M_{\rm c,med} \sim 3$ M$_{\rm J}$ and a mean density of $\rho \sim 11$ g cm$^{-3}$. M62I has an orbital period of 0.51 days and a $M_{\rm c,med} \sim 0.15$ M$_{\odot}$. Neither of these low-mass systems exhibit eclipses. M62J has only been detected in the two UHF band (816 MHz) observations with a flux density $S_{816} = 0.08$ mJy. The non-detection in the L-band (1284 MHz) indicates it has a relatively steep spectrum ($\beta < -3.1$). We also present 23-yr-long timing solutions obtained using data from the Parkes "Murriyang", Effelsberg and MeerKAT telescopes for the six previously known pulsars. For all these pulsars, we measured the second spin-period derivatives and the rate of change of orbital period caused by the gravitational field of the cluster, and their proper motions. From these measurements, we conclude that the pulsars' maximum accelerations are consistent with the maximum cluster acceleration assuming a core-collapsed mass distribution. Studies of the eclipses of the redback M62B and the black widow M62E at four and two different frequency bands, respectively, reveal a frequency dependence with longer and asymmetric eclipses at lower frequencies. The presence of only binary MSPs in this cluster challenges models which suggest that the MSP population of core-collapsed clusters should be dominated by isolated MSPs., Comment: 22 pages, 8 figures, 6 tables. Accepted for publication in MNRAS

Redbacks are millisecond pulsar binaries with low mass, irradiated companions. These systems have a rich phenomenology that can be used to probe binary evolution models, pulsar wind physics, and the neutron star mass distribution. A number of high-confidence redback candidates have been identified through searches for variable optical and X-ray sources within the localisation regions of unidentified but pulsar-like Fermi-LAT gamma-ray sources. However, these candidates remain unconfirmed until pulsations are detected. As part of the TRAPUM project, we searched for radio pulsations from six of these redback candidates with MeerKAT. We discovered three new radio millisecond pulsars, PSRs J0838$-$2527, J0955$-$3947 and J2333$-$5526, confirming their redback nature. PSR J0838$-$2827 remained undetected for two years after our discovery despite repeated observations, likely due to evaporated material absorbing the radio emission for long periods of time. While, to our knowledge, this system has not undergone a transition to an accreting state, the disappearance, likely caused by extreme eclipses, illustrates the transient nature of spider pulsars and the heavy selection bias in uncovering their radio population. Radio timing enabled the detection of gamma-ray pulsations from all three pulsars, from which we obtained 15-year timing solutions. All of these sources exhibit complex orbital period variations consistent with gravitational quadrupole moment variations in the companion stars. These timing solutions also constrain the binary mass ratios, allowing us to narrow down the pulsar masses. We find that PSR J2333$-$5526 may have a neutron star mass in excess of 2 M$_{\odot}$., Comment: 19 pages, 7 figures, accepted for publication in MNRAS

We used the Imaging X-ray Polarimetry Explorer (IXPE) satellite to measure, for the first time, the 2-8 keV polarization of NGC 1068. We pointed IXPE for a net exposure time of 1.15 Ms on the target, in addition to two ~ 10 ks each Chandra snapshots in order to account for the potential impact of several ultraluminous X-ray source (ULXs) within IXPE's field-of-view. We measured a 2 - 8 keV polarization degree of 12.4% +/- 3.6% and an electric vector polarization angle of 101{\deg} +/- 8{\deg} at 68% confidence level. If we exclude the spectral region containing the bright Fe K lines and other soft X-ray lines where depolarization occurs, the polarization fraction rises up to 21.3% +/- 6.7% in the 3.5 - 6.0 keV band, with a similar polarization angle. The observed polarization angle is found to be perpendicular to the parsec scale radio jet. Using a combined Chandra and IXPE analysis plus multi-wavelength constraints, we estimated that the circumnuclear "torus" may sustain a half-opening angle of 50{\deg} - 55{\deg} (from the vertical axis of the system). Thanks to IXPE, we have measured the X-ray polarization of NGC 1068 and found comparable results, both in terms of polarization angle orientation with respect to the radio-jet and torus half-opening angle, to the X-ray polarimetric measurement achieved for the other archetypal Compton-thick AGN : the Circinus galaxy. Probing the geometric arrangement of parsec-scale matter in extragalactic object is now feasible thanks to X-ray polarimetry., Comment: 9 pages, 5 figures, 2 tables, submited to A&A

Recently, the Galactic magnetar SGR J1935+2154 has garnered attention due to its emission of an extremely luminous radio burst, reminiscent of Fast Radio Bursts (FRBs). SGR J1935+2154 is one of the most active magnetars, displaying flaring events nearly every year, including outbursts as well as short and intermediate bursts. Here, we present our results on the properties of the persistent and bursting X-ray emission from SGR J1935+2154, during the initial weeks following its outburst on October 10, 2022. The source was observed with XMM-Newton and NuSTAR (quasi-)simultaneously during two epochs, separated by $\sim$5 days. The persistent emission spectrum is well described by an absorbed blackbody plus power-law model up to an energy of $\sim$25 keV. No significant changes were observed in the blackbody temperature ($kT_{\rm BB}\sim$ 0.4 keV) and emitting radius ($R_{\rm BB}\sim$ 1.9 km) between the two epochs. However, we observed a slight variation in the power-law parameters. Moreover, we detected X-ray pulsations in all the datasets and derived a spin period derivative of $\dot{P} = 5.52(5) \times 10^{-11}$ ss. This is 3.8 times larger than the value measured after the first recorded outburst in 2014. Additionally, we performed quasi-simultaneous radio observations using three 25--32-m class radio telescopes for a total of 92.5 hr to search for FRB-like radio bursts and pulsed emission. However, our analysis did not reveal any radio bursts or periodic emission., Comment: 16 pages, 6 figures, accepted for publication on ApJ

We report on X-ray polarization measurements of the extra-galactic Crab-like PSR B0540-69 and its Pulsar Wind Nebula (PWN) in the Large Magellanic Cloud (LMC), using a ~850 ks Imaging X-ray Polarimetry Explorer (IXPE) exposure. The PWN is unresolved by IXPE. No statistically significant polarization is detected for the image-averaged data, giving a 99% confidence polarization upper limit (MDP99) of 5.3% in 2-8 keV energy range. However, a phase-resolved analysis detects polarization for both the nebula and pulsar in the 4-6 keV energy range. For the PWN defined as the off-pulse phases, the polarization degree (PD) of (24.5 ${\pm}$ 5.3)% and polarization angle (PA) of (78.1 ${\pm}$ 6.2){\deg} is detected at 4.6${\sigma}$ significance level, consistent with the PA observed in the optical band. In a single on-pulse window, a hint of polarization is measured at 3.8${\sigma}$ with polarization degree of (50.0 ${\pm}$ 13.1)% and polarization angle of (6.2 ${\pm}$ 7.4){\deg}. A 'simultaneous' PSR/PWN analysis finds two bins at the edges of the pulse exceeding 3${\sigma}$ PD significance, with PD of (68 ${\pm}$ 20)% and (62 ${\pm}$ 20)%; intervening bins at 2-3${\sigma}$ significance have lower PD, hinting at additional polarization structure., Comment: 10 pages, 7 figures, 4 tables, author's version of the paper accepted for publication in ApJ

Weakly magnetized neutron stars in X-ray binaries show complex phenomenology with several spectral components that can be associated with the accretion disk, boundary and/or spreading layer, a corona, and a wind. Spectroscopic information alone is, however, not enough to disentangle these components. Additional information about the nature of the spectral components and in particular the geometry of the emission region can be provided by X-ray polarimetry. One of the objects of the class, a bright, persistent, and rather peculiar galactic Type I X-ray burster was observed with the Imaging X-ray Polarimetry Explorer (IXPE) and the X-ray Multi-Mirror Mission Newton (XMM-Newton). Using the XMM-Newton data we estimated the current state of the source as well as detected strong absorption lines associated with the accretion disk wind. IXPE data showed the source to be significantly polarized in the 2-8 keV energy band with the overall polarization degree (PD) of 1.4% at a polarization angle (PA) of -2 degrees (errors at 68% confidence level). During the two-day long observation, we detected rotation of the PA by about 70 degrees with the corresponding changes in the PD from 2% to non-detectable and then up to 5%. These variations in polarization properties are not accompanied by visible changes in spectroscopic characteristics. The energy-resolved polarimetric analysis showed a significant change in polarization, from being strongly dependent on energy at the beginning of the observation to being almost constant with energy in the later parts of the observation. As a possible interpretation, we suggest the presence of a constant component of polarization, strong wind scattering, or different polarization of the two main spectral components with individually peculiar behavior. The rotation of the PA suggests a 30-degree misalignment of the neutron star spin from the orbital axis., Comment: 12 pages, 12 figures

Magnetar-like activity has been observed in a large variety of neutron stars. PSR J1846-0258 is a young 327 ms radio-quiet pulsar with a large rotational power ($\sim 8 \times 10^{36}$ erg s$^{-1}$), and resides at the center of the supernova remnant Kes 75. It is one of the rare examples of a high magnetic field pulsar showing characteristics both of magnetars and radio pulsars, and can thus provide important clues on the differences in the emission mechanisms between these two classes. In 2006, PSR J1846-0258 was detected to undergo an outburst for the first time, accompanied by a large flux increase, millisecond X-ray bursts, significant spectral changes and a large timing glitch. In the period between May-June 2020, after fourteen years of quiescent stable emission, the source underwent a second magnetar-like outburst, which was followed up with several observations by NICER, XMM-Newton, NuSTAR and Swift. In this work, we report on the long-term timing and X-ray spectral properties of the source following the 2020 outburst, and place upper limits on any source activity at radio wavelengths. We demonstrate that the pulsed flux increased by a factor $> 6$ during the outburst, followed by non-trivial variability in the spin-down rate. Our timing analysis shows that the spin frequency and its derivative are clearly affected by magnetospheric activity due to the outburst. We find hints for an oscillation in the frequency derivative with a timescale of 50-60 days, recovering later on to stable quiescence., Comment: 22 pages, 10 figures; accepted for publication by ApJ; comments are welcome

Observations of linear polarization in the 2-8 keV energy range with the Imaging X-ray Polarimetry Explorer (IXPE) explore the magnetic field geometry and dynamics of the regions generating non-thermal radiation in relativistic jets of blazars. These jets, particularly in blazars whose spectral energy distribution peaks at X-ray energies, emit X-rays via synchrotron radiation from high-energy particles within the jet. IXPE observations of the X-ray selected BL Lac-type blazar 1ES 1959+650 in 2022 May 3-4 showed a significant linear polarization degree of $\Pi_\mathrm{x} = 8.0\% \pm 2.3\%$ at an electric-vector position angle $\psi_\mathrm{x} = 123^\circ \pm 8^\circ$. However, in 2022 June 9-12, only an upper limit of $\Pi_\mathrm{x} \leq 5.1\%$ could be derived (at the 99% confidence level). The degree of optical polarization at that time $\Pi_\mathrm{O} \sim 5\%$ is comparable to the X-ray measurement. We investigate possible scenarios for these findings, including temporal and geometrical depolarization effects. Unlike some other X-ray selected BL Lac objects, there is no significant chromatic dependence of the measured polarization in 1ES 1959+650, and its low X-ray polarization may be attributed to turbulence in the jet flow with dynamical timescales shorter than 1 day., Comment: 16 pages, 4 figures, accepted for publication in The Astrophysical Journal

We present the first X-ray polarimetric study of the dipping accreting neutron star 4U 1624$-$49 with the Imaging X-ray Polarimetry Explorer (IXPE). We report a detection of polarization in the non-dip time intervals with a confidence level of 99.99%. We find an average polarization degree (PD) of $3.1\pm0.7$% and a polarization angle of $81\pm6$ degrees east of north in the 2-8 keV band. We report an upper limit on the PD of 22% during the X-ray dips with 95% confidence. The PD increases with energy, reaching from $3.0\pm0.9$% in the 4-6 keV band to $6\pm2$% in the 6-8 keV band. This indicates the polarization likely arises from Comptonization. The high PD observed is unlikely to be produced by Comptonization in the boundary layer or spreading layer alone. It can be produced by the addition of an extended geometrically thin slab corona covering part of the accretion disk, as assumed in previous models of dippers, and/or a reflection component from the accretion disk.

We report the discovery of X-ray polarization from the X-ray-bright filament. G0.13-0.11 in the Galactic center (GC) region. This filament features a bright, hard X-ray source that is most plausibly a pulsar wind nebula (PWN) and an extended and structured diffuse component. Combining the polarization signal from IXPE with the imaging/spectroscopic data from Chandra, we find that X-ray emission of G0.13-0.11 is highly polarized PD=$57(\pm18)$% in the 3-6 keV band, while the polarization angle is PA=$21^\circ(\pm9^\circ)$. This high degree of polarization proves the synchrotron origin of the X-ray emission from G0.13-0.11. In turn, the measured polarization angle implies that the X-ray emission is polarized approximately perpendicular to a sequence of nonthermal radio filaments that may be part of the GC Radio Arc. The magnetic field on the order of $100\,{\rm\mu G}$ appears to be preferentially ordered along the filaments. The above field strength is the fiducial value that makes our model self-consistent, while the other conclusions are largely model independent., Comment: Accepted for publication in A&A

The Imaging X-ray Polarimetry Explorer (IXPE) measured with high significance the X-ray polarization of the brightest Z-source Scorpius X-1, resulting in the nominal 2-8 keV energy band in a polarization degree of 1.0(0.2)% and a polarization angle of 8(6){\deg} at 90% of confidence level. This observation was strictly simultaneous with observations performed by NICER, NuSTAR, and Insight-HXMT, which allowed for a precise characterization of its broad-band spectrum from soft to hard X-rays. The source has been observed mainly in its soft state, with short periods of flaring. We also observed low-frequency quasi-periodic oscillations. From a spectro-polarimetric analysis, we associate a polarization to the accretion disk at <3.2% at 90% of confidence level, compatible with expectations for an electron-scattering dominated optically thick atmosphere at the Sco X-1 inclination of 44{\deg}; for the higher-energy Comptonized component, we obtain a polarization of 1.3(0.4)%, in agreement with expectations for a slab of Thomson optical depth of ~7 and an electron temperature of ~3 keV. A polarization rotation with respect to previous observations by OSO-8 and PolarLight, and also with respect to the radio-jet position angle, is observed. This result may indicate a variation of the polarization with the source state that can be related to relativistic precession or to a change in the corona geometry with the accretion flow.

We report on an observational campaign on the bright black hole X-ray binary Swift J1727.8$-$1613 centered around five observations by the Imaging X-ray Polarimetry Explorer (IXPE). These observations track for the first time the evolution of the X-ray polarization of a black hole X-ray binary across a hard to soft state transition. The 2--8 keV polarization degree decreased from $\sim$4\% to $\sim$3\% across the five observations, but the polarization angle remained oriented in the North-South direction throughout. Based on observations with the Australia Telescope Compact Array (ATCA), we find that the intrinsic 7.25 GHz radio polarization aligns with the X-ray polarization. Assuming the radio polarization aligns with the jet direction (which can be tested in the future with higher spatial resolution images of the jet), our results imply that the X-ray corona is extended in the disk plane, rather than along the jet axis, for the entire hard intermediate state. This in turn implies that the long ($\gtrsim$10 ms) soft lags that we measure with the Neutron star Interior Composition ExploreR (NICER) are dominated by processes other than pure light-crossing delays. Moreover, we find that the evolution of the soft lag amplitude with spectral state does not follow the trend seen for other sources, implying that Swift J1727.8$-$1613 is a member of a hitherto under-sampled sub-population., Comment: Accepted for publication in ApJ. 21 pages, 8 figures

Cir X-1 is a neutron star X-ray binary characterized by strong variations in flux during its eccentric $\sim$16.6 days orbit. There are also strong variations in the spectral state, and historically it has shown both atoll and Z state properties. We observed the source with the Imaging X-ray Polarimetry Explorer during two orbital segments, 6 days apart, for a total of 263~ks. We find an X-ray polarization degree in these segments of $1.6\%\pm0.3\%$ and $1.4\%\pm0.3\%$ at polarization angles of $37^\circ\pm5^\circ$ and $-12^\circ\pm7^\circ$, respectively. Thus we observed a rotation of the polarization angle by $49^\circ\pm8^\circ$ along the orbit. Because variations of accretion flow, and then of the hardness ratio, are expected during the orbit, we also studied the polarization binned in hardness ratio, and found the polarization angle differing by $67^\circ\pm11^\circ$ between the lowest and highest values of the hardness ratio. We discuss possible interpretations of this result that could indicate a possible misalignment between the symmetry axes of the accretion disk and the Comptonizing region caused by the misalignment of the neutron star's angular momentum with respect to the orbital one., Comment: 15 pages, 6 figures, Accepted for publication in ApJL

In an accreting X-ray pulsar, a neutron star accretes matter from a stellar companion through an accretion disk. The high magnetic field of the rotating neutron star disrupts the inner edge of the disc, funneling the gas to flow onto the magnetic poles on its surface. Hercules X-1 is in many ways the prototypical X-ray pulsar; it shows persistent X-ray emission and it resides with its companion HZ Her, a two-solar-mass star, at about 7~kpc from Earth. Its emission varies on three distinct timescales: the neutron star rotates every 1.2~seconds, it is eclipsed by its companion each 1.7~days, and the system exhibits a superorbital period of 35~days which has remained remarkably stable since its discovery. Several lines of evidence point to the source of this variation as the precession of the accretion disc, the precession of the neutron star or both. Despite the many hints over the past fifty years, the precession of the neutron star itself has yet not been confirmed or refuted. We here present X-ray polarization measurements with the Imaging X-ray Polarimetry Explorer (IXPE) which probe the spin geometry of the neutron star. These observations provide direct evidence that the 35-day-period is set by the free precession of the neutron star crust, which has the important implication that its crust is somewhat asymmetric fractionally by a few parts per ten million. Furthermore, we find indications that the basic spin geometry of the neutron star is altered by torques on timescale of a few hundred days., Comment: 37 pages

We report on IXPE, NICER and XMM-Newton observations of the magnetar 1E 2259+586. We find that the source is significantly polarized at about or above 20% for all phases except for the secondary peak where it is more weakly polarized. The polarization degree is strongest during the primary minimum which is also the phase where an absorption feature has been identified previously (Pizzocaro et al. 2019). The polarization angle of the photons are consistent with a rotating vector model with a mode switch between the primary minimum and the rest of the rotation of the neutron star. We propose a scenario in which the emission at the source is weakly polarized (as in a condensed surface) and, as the radiation passes through a plasma arch, resonant cyclotron scattering off of protons produces the observed polarized radiation. This confirms the magnetar nature of the source with a surface field greater than about 1015 G, Comment: 13 pages, 13 figures, 5 tables

The HTRU-S Low Latitude survey data within 1$^{\circ}$of the Galactic Centre (GC) were searched for pulsars using the Fast Folding Algorithm (FFA). Unlike traditional Fast Fourier Transform (FFT) pipelines, the FFA optimally folds the data for all possible periods over a given range, which is particularly advantageous for pulsars with low-duty cycle. For the first time, a search over acceleration was included in the FFA to improve its sensitivity to binary pulsars. The steps in dispersion measure (DM) and acceleration were optimised, resulting in a reduction of the number of trials by 86 per cent. This was achieved over a search period range from 0.6-s to 432-s, i.e. 10 per cent of the observation time (4320s), with a maximum DM of 4000 pc cm$^{-3}$ and an acceleration range of $\pm 128$m s$^{-2}$. The search resulted in the re-detections of four known pulsars, including a pulsar which was missed in previous FFT processing of this survey. This result indicates that the FFA pipeline is more sensitive than the FFT pipeline used in the previous processing of the survey within our parameter range. Additionally, we discovered a 1.89-s pulsar, PSR J1746-2829, with a large DM, located~0.5 from the GC. Follow-up observations revealed that this pulsar has a relatively flat spectrum($\alpha=-0.9\pm0.1$) and has a period derivative of $\sim1.3\times10^{-12}$ s s$^{-1}$, implying a surface magnetic field of $\sim5.2\times10^{13}$ G and a characteristic age of $\sim23000$ yr. While the period, spectral index, and surface magnetic field strength are similar to many radio magnetars, other characteristics such as high linear polarization are absent., Comment: 12 pages, 6 figures, 4 tables, Accepted for publication on Monthly Notices of the Royal Astronomical Society

We present X-ray polarimetry observations from the Imaging X-ray Polarimetry Explorer (IXPE) of three low spectral peak and one intermediate spectral peak blazars, namely 3C 273, 3C 279, 3C 454.3, and S5 0716+714. For none of these objects was IXPE able to detect X-ray polarization at the 3$\sigma$ level. However, we placed upper limits on the polarization degree at $\sim$10-30\%. The undetected polarizations favor models where the X-ray band is dominated by unpolarized photons upscattered by relativistic electrons in the jets of blazars, although hadronic models are not completely eliminated. We discuss the X-ray polarization upper limits in the context of our contemporaneous multiwavelength polarization campaigns., Comment: 19 pages, 10 figures, 6 tables; submitted to the Astrophysical Journal

We present the results of the first X-ray polarimetric observation of the low-mass X-ray binary 4U 1957+115, performed with the Imaging X-ray Polarimetry Explorer in May 2023. The binary system has been in a high-soft spectral state since its discovery and is thought to host a black hole. The $\sim$571 ks observation reveals a linear polarisation degree of $1.9\% \pm 0.6\%$ and a polarisation angle of $-41^\circ.8 \pm 7^\circ.9$ in the 2-8 keV energy range. Spectral modelling is consistent with the dominant contribution coming from the standard accretion disc, while polarimetric data suggest a significant role of returning radiation: photons that are bent by strong gravity effects and forced to return to the disc surface, where they can be reflected before eventually reaching the observer. In this setting, we find that models with a black hole spin lower than 0.96 and an inclination lower than $50^\circ$ are disfavoured., Comment: 12 pages, 10 figures, 2 tables, accepted for publication in A&A

We report on the coordinated observations of the neutron star low-mass X-ray binary (NS-LMXB) \gx in X-rays (IXPE, NICER, Nustar and INTEGRAL), optical (REM and LCO), near-infrared (REM), mid-infrared (VLT VISIR), and radio (ATCA). This Z-source was observed by \IXPE twice in March-April 2023 (Obs. 1 and 2). In the radio band, the source was detected, but only upper-limits to the linear polarization were obtained at a $3\sigma$ level of $6.1\%$ at 5.5 GHz and $5.9\%$ at 9 GHz in Obs.~1 and $12.5\%$ at 5.5~GHz and $20\%$ at 9~GHz in Obs.~2. The mid-IR, near-IR and optical observations suggest the presence of a compact jet which peaks in the mid- or far-IR. The X-ray polarization degree was found to be $3.7\% \pm 0.4 \%$ (at $90\%$ confidence level) during Obs.~1 when the source was in the horizontal branch of the Z-track and $1.8\% \pm 0.4 \%$ during Obs.~2 when the source was in the normal-flaring branch. These results confirm the variation of polarization degree as a function of the position of the source in the color-color diagram as for previously observed Z-track sources (Cyg~X-2 and XTE~1701$-$462). Evidence for a variation of the polarization angle $\sim 20^\circ$ with energy is found in both observations, likely related to the different, non-orthogonal polarization angles of the disk and Comptonization components which peak at different energies., Comment: Submitted to Astronomy and Astrophysics on 06 July 2023. Accepted on 21 November 2023

We conducted a polarimetry campaign from radio to X-ray wavelengths of the high-synchrotron-peak (HSP) blazar Mrk 421, including Imaging X-ray Polarimetry Explorer (IXPE) measurements on 2022 December 6-8. We detected X-ray polarization of Mrk 421 with a degree of $\Pi_{\rm X}$=14$\pm$1$\%$ and an electric-vector position angle $\psi_{\rm X}$=107$\pm$3$^{\circ}$ in the 2-8 keV band. From the time variability analysis, we find a significant episodic variation in $\psi_{\rm X}$. During 7 months from the first IXPE pointing of Mrk 421 in 2022 May, $\psi_{\rm X}$ varied across the range of 0$^{\circ}$ to 180$^{\circ}$, while $\Pi_{\rm X}$ maintained similar values within $\sim$10-15$\%$. Furthermore, a swing in $\psi_{\rm X}$ in 2022 June was accompanied by simultaneous spectral variations. The results of the multiwavelength polarimetry show that the X-ray polarization degree was generally $\sim$2-3 times greater than that at longer wavelengths, while the polarization angle fluctuated. Additionally, based on radio, infrared, and optical polarimetry, we find that rotation of $\psi$ occurred in the opposite direction with respect to the rotation of $\psi_{\rm X}$ over longer timescales at similar epochs. The polarization behavior observed across multiple wavelengths is consistent with previous IXPE findings for HSP blazars. This result favors the energy-stratified shock model developed to explain variable emission in relativistic jets. The accompanying spectral variation during the $\psi_{\rm X}$ rotation can be explained by a fluctuation in the physical conditions, e.g., in the energy distribution of relativistic electrons. The opposite rotation direction of $\psi$ between the X-ray and longer-wavelength polarization accentuates the conclusion that the X-ray emitting region is spatially separated from that at longer wavelengths., Comment: 17 pages, 13 figures, 4 tables; Accepted for publication in A&A

We present the results of observations aimed at discovering and studying pulsars in the core-collapsed globular cluster (GC) NGC 6522 performed by the MeerTIME and TRAPUM Large Survey Project with the MeerKAT telescope. We have discovered two new isolated pulsars bringing the total number of known pulsars in the cluster to six. PSR J1803$-$3002E is a mildly recycled pulsar with spin period of 17.9 ms while pulsar PSR J1803$-$3002F is a slow pulsar with spin period of 148.1 ms. The presence of isolated and slow pulsars is expected in NGC 6522 and confirms the predictions of previous theories for clusters at this stage in evolution. We further present a tentative timing solution for the millisecond pulsar (MSP) PSR J1803$-$3002C combining older observations taken with the Parkes 64m radio telescope, Murriyang. This solution implies a relatively small characteristic age of the pulsar in contrast with the old age of the GC. The presence of a slow pulsar and an apparently young MSP, both rare in GCs, suggests that their formation might be linked to the evolutionary stage of the cluster., Comment: 8 pages, 5 figues, 3 tables. Accepted for publication in Astronomy and Astrophysics

We present polarization measurements in the $2-8 \thinspace \mathrm{keV}$ band from blazar 1ES 0229+200, the first extreme high synchrotron peaked source to be observed by the Imaging X-ray Polarimetry Explorer (IXPE). Combining two exposures separated by about two weeks, we find the degree of polarization to be $\Pi_{X} = 17.9 \pm 2.8 \%$ at an electric-vector position angle $\psi_X = 25.0 \pm 4.6^{\circ}$ using a spectropolarimetric fit from joint IXPE and XMM-Newton observations. There is no evidence for the polarization degree or angle varying significantly with energy or time on both short time scales (hours) or longer time scales (days). The contemporaneous polarization degree at optical wavelengths was $>$7$\times$ lower, making 1ES 0229+200 the most strongly chromatic blazar yet observed. This high X-ray polarization compared to the optical provides further support that X-ray emission in high-peaked blazars originates in shock-accelerated, energy-stratified electron populations, but is in tension with many recent modeling efforts attempting to reproduce the spectral energy distribution of 1ES 0229+200 which attribute the extremely high energy synchrotron and Compton peaks to Fermi acceleration in the vicinity of strongly turbulent magnetic fields., Comment: 17 Pages, 6 Figures, Resubmitted to ApJ after addressing referee comments

We describe IXPE polarization observations of the Pulsar Wind Nebula (PWN) MSH15-52, the `Cosmic Hand'. We find X-ray polarization across the PWN, with B field vectors generally aligned with filamentary X-ray structures. High significance polarization is seen in arcs surrounding the pulsar and toward the end of the `jet', with polarization degree PD>70%, thus approaching the maximum allowed synchrotron value. In contrast, the base of the jet has lower polarization, indicating a complex magnetic field at significant angle to the jet axis. We also detect significant polarization from PSR B1509-58 itself. Although only the central pulse-phase bin of the pulse has high individual significance, flanking bins provide lower significance detections and, in conjunction with the X-ray image and radio polarization, can be used to constrain rotating vector model solutions for the pulsar geometry., Comment: To appear in the Astrophysical Journal

We report the first detection of the X-ray polarization of the bright transient Swift J1727.8-1613 with the Imaging X-ray Polarimetry Explorer. The observation was performed at the beginning of the 2023 discovery outburst, when the source resided in the bright hard state. We find a time- and energy-averaged polarization degree of 4.1%+/-0.2% and a polarization angle of 2.2+/-1.3 degrees (errors at 68% confidence level; this translates to about 20-sigma significance of the polarization detection). This finding suggests that the hot corona emitting the bulk of the detected X-rays is elongated, rather than spherical. The X-ray polarization angle is consistent with that found in sub-mm wavelengths. Since the sub-mm polarization was found to be aligned with the jet direction in other X-ray binaries, this indicates that the corona is elongated orthogonal to the jet., Comment: 10 pages, 3 figures, submitted