Search

Pulsar glitch is a phenomenon characterized by abrupt changes in the spin period over less than a minute. We present a comprehensive analysis of glitches in four gamma-ray pulsars by combining the timing observation data from Fermi -Large Area Telescope (Fermi -LAT) and Parkes 64 m radio telescope. The timing data of five pulsars, namely PSRs J1028–5819, J1420–6048, J1509–5850, J1709–4429 (B1706–44), and J1718–3825, are examined over 14 yr of observations for each. A total of 12 glitches are identified in four pulsars, including a previously unreported glitch. That is, a new small glitch is identified for PSR J1718–3825 in MJD |$\sim$| 59121(8), with a fractional glitch size of |$\Delta \nu /\nu \sim 1.9(2) \times 10^{-9}$|⁠. For PSR J1420–6048, our investigation confirms the presence of two linear recovery terms during the evolution of |$\dot{\nu }$| following glitches 4, 6, and 8. Moreover, an exponential recovery process was identified after glitch 8, with a recovery fraction (Q) of |$Q = 0.0131(5)$| and a corresponding time-scale of |$\tau _{\rm d} = 100(6)$| d. Regarding the fourth glitch of PSR J1709–4429, our analysis reveals the presence of two exponential recovery terms with degree of recovery and decay time-scales Q 1 = 0.0104(5), |$\tau _{\rm d1}=72(4)$| d, and Q 2 = 0.006(1), |$\tau _{\rm d2}=4.2(6)$|  d, respectively. For the remaining previously reported glitches, we also refine the glitch epochs and recovery process through precise fitting of the timing data. We discuss how multiband data of glitches can help better characterize the glitch recoveries and constrain the underlying physics of glitch events. Our findings demonstrate that the accumulation of observational data reveals the rich complexity of the glitch phenomenon, aiding in the search for a well-established interpretation. [ABSTRACT FROM AUTHOR]

Radio searches for single pulses provide the opportunity to discover one-off events, fast transients, and some pulsars that might otherwise be missed by conventional periodicity searches. The MeerTRAP real-time search pipeline operates commensally to observations with the MeerKAT telescope. Here, we report on 26 new Galactic radio transients, mostly rotating radio transients (RRATs) and also the detection of one RRAT and two pulsars that were independently discovered by other surveys. The dispersion measures of two of the new sources marginally exceed the Galactic contribution depending on the electron density model used. Using a simple method of fitting a Gaussian function to individual pulses, and obtaining positions of arcsec accuracy from image-based localizations using channelized voltage data from our transient buffer, we have derived timing solutions spanning multiple years for five sources. The timing parameters imply ages of several Myr and low surface magnetic field strengths which is characteristic of RRATs. We were able to measure spin periods for eight more transients, including one source which appears to rotate every 17.5 s. A majority of the sources have only been seen in one observation, sometimes despite multiple return visits to the field. Some sources exhibit complex emission features like component switching and periodic microstructure. [ABSTRACT FROM AUTHOR]

We perform microscopic simulations of the thermal relaxation of warm neutral plasmas of astrophysical importance. Using molecular dynamics we study the thermal relaxation of a hot neutral fluid of finite-size neutron-rich ions kept in a fixed-wall cool box. We show how the interplay among particle size, Yukawa interaction range, and density are key to understand the features of the time-dependent thermal relaxation curve |$T(t)$|⁠. We show that, under certain conditions, these systems exhibit faster cooling from increasingly larger initial temperature values. They also display non-Newtonian thermal behaviour, including oscillations, that can be effectively interpreted as the consequence of the existence of a non-trivial system memory function. Finally, we consider the impact of multicomponent admixtures in the simulated system. We discuss these results and their possible extensions to astrophysical scenarios, where screened plasmas with Yukawa-like electrostatic potentials are usually involved. [ABSTRACT FROM AUTHOR]

r -mode oscillations in rotating neutron stars are a source of continuous gravitational radiation. We investigate the excitation of r -modes by the mechanical impact on the neutron star surface of stochastically accreted clumps of matter, assuming that the Chandrasekhar–Friedman–Schutz instability is not triggered. The star is idealized as a slowly rotating, unmagnetized, one-component fluid with a barotropic equation of state in Newtonian gravity. It is found that the r -mode amplitude depends weakly on the equation of state but sensitively on the rotation frequency |$\nu _{\rm s}$|⁠. The gravitational wave strain implicitly depends on the equation of state through the damping time-scale. The root-mean-square strain is |$h_{\rm rms} \approx 10^{-35} (\nu _{\rm s}/ 10\, {\rm Hz})^{2} (R_*/10\, {\rm km})^2 (\Delta t_{\rm acc}/1\, {\rm yr})^{1/2} (f_{\rm acc}/1\, {\rm kHz})^{-1/2} (\dot{M}/10^{-8}\text{M}_{\odot } \, \text{yr}^{-1}) (v/0.4c) (d/1\, {\rm kpc})^{-1}$|⁠ , which is comparable to the strain from g -, p -, and f -modes excited by stochastic accretion, where |$R_*$| is the radius of the star, |$\Delta t_{\rm acc}$| is the uninterrupted duration of an accretion episode, |$f_{\rm acc}$| is the mean clump impact frequency, |$\dot{M}$| is the accretion rate, v is the impact speed, and d is the distance of the star from the Earth. An observational test is proposed, based on the temporal autocorrelation function of the gravitational wave signal, to discern whether the Chandrasekhar–Friedman–Schutz instability switches on and coexists with impact-excited r -modes before or during a gravitational wave observation. [ABSTRACT FROM AUTHOR]

We present a detailed spectral study of an intermittent-AMXP Aql X-1 during the pulse-on and pulse-off stages by using the archival Rossi X-ray Timing Explorer (RXTE) data. We first perform temporal analysis by using Z |$_n^2$| technique in three different energy bands, 3.0–13.0, 13.0–23.0, and 23.0–33.0 keV, for the last 128 s time segment of the RXTE data including pulse-on region. We show that the pulse is the most significant in the softest band. We, then, show that the spectrum is represented the best via combination of absorbed blackbody, disc blackbody, and a Gaussian line. We modelled the last four segments of the data 30188-03-05-00 to better compare pulse-on and pulse-off stages. We found a vague residual in the spectral fit of the pulse-on segment between |$\sim$| 3.0 and 13.0 keV, which agrees with the result of temporal analysis. We show that the residual may be represented with an extra blackbody component with the temperature of 1.75 keV and the radius of 0.75 |$\pm$| 0.49 km. For deeper analysis, we performed phase-resolved spectroscopy to the last 128 s, pulse-on , segment. We obtain two separate spectra for the spin phase ranges of 0.75–0.25 (pulse-high) and 0.25–0.75 (pulse-low), and followed the same procedure. We display that the residual becomes more clear for pulse-high compared to the pulse-low. We report that the additional blackbody component, which models the residual, indicates a hotspot from the surface of the neutron star with the radius of 1.65 |$\pm$| 0.74 km whose temperature is 1.65 keV. [ABSTRACT FROM AUTHOR]

We investigate whether measurements of the neutron star mass and radius or the tidal deformability can provide information about the presence of hyperons inside a neutron star. This is achieved by considering two inference models, with and without hyperons, based on a field theoretical approach. While current observations do not distinguish between the two scenarios, we have shown that data simulating expected observations from future large-area X-ray timing telescopes could provide some information through Bayes' factors. Inference using simulated data generated from an equation of state containing hyperons decisively favours the hyperonic model over the nucleonic model. However, a 2 per cent uncertainty in the mass and radius determination may not be sufficient to constrain the parameters of the model when only six neutron star mass–radius measurements are considered. [ABSTRACT FROM AUTHOR]

There has been a rapid increase in the known fast radio burst (FRB) population, yet the progenitor(s) of these events have remained an enigma. A small number of FRBs have displayed some level of quasi-periodicity in their burst profile, which can be used to constrain their plausible progenitors. However, these studies suffer from the lack of polarization data which can greatly assist in constraining possible FRB progenitors and environments. Here, we report on the detection and characterisation of FRB 20230708A by the Australian Square Kilometre Array Pathfinder (ASKAP), a burst which displays a rich temporal and polarimetric morphology. We model the burst time series to test for the presence of periodicity, scattering and scintillation. We find a potential period of T  = 7.267 ms within the burst, but with a low statistical significance of 1.77 |$\sigma$|⁠. Additionally, we model the burst's time- and frequency-dependent polarization to search for the presence of (relativistic and non-relativistic) propagation effects. We find no evidence to suggest that the high circular polarization seen in FRB 20230708A is generated by Faraday conversion. The majority of the properties of FRB 20230708A are broadly consistent with a (non-millisecond) magnetar model in which the quasi-periodic morphology results from microstructure in the beamed emission, but other explanations are not excluded. [ABSTRACT FROM AUTHOR]

We conducted a GPU-accelerated reprocessing of |$\sim 87~{{\ \rm per\ cent}}$| of the archival data from the High Time Resolution Universe South Low Latitude (HTRU-S LowLat) pulsar survey by implementing a pulsar search pipeline that was previously used to reprocess the Parkes Multibeam Pulsar Survey (PMPS). We coherently searched the full 72-min observations of the survey with an acceleration search range up to |$|50|\, \rm m\, s^{-2}$|⁠ , which is most sensitive to binary pulsars experiencing nearly constant acceleration during 72 min of their orbital period. Here we report the discovery of 71 pulsars, including six millisecond pulsars, of which five are in binary systems, and seven pulsars with very high dispersion measures (DM |$\gt 800 \, \rm pc \, cm^{-3}$|⁠). These pulsar discoveries largely arose by folding candidates to a much lower spectral signal-to-noise ratio than in previous surveys and by exploiting the coherence of folding over the incoherent summing of the Fourier components to discover new pulsars as well as candidate classification techniques. We show that these pulsars could be fainter and on average more distant as compared with both the previously reported 100 HTRU-S LowLat pulsars and the background pulsar population in the survey region. We have assessed the effectiveness of our search method and the overall pulsar yield of the survey. We show that through this reprocessing we have achieved the expected survey goals, including the predicted number of pulsars in the survey region, and discuss the major causes why these pulsars were missed in previous processing of the survey. [ABSTRACT FROM AUTHOR]

We present GTC-OSIRIS phase-resolved optical spectroscopy of three compact binary millisecond pulsars, or 'spiders': PSR J1048+2339, PSR J1810+1744, and (for the first time) PSR J1908+2105. For the companion in each system, the temperature is traced throughout its orbit, and radial velocities are measured. The radial velocities are found to vary with the absorption features used when measuring them, resulting in different radial velocity curve semi-amplitudes: for J1048 (⁠|$K_\mathrm{metals, red} = 344 \pm 4{\mathrm{\, km\, s^{-1}}}{}$| and |$K_\mathrm{metals, blue} = 372 \pm 3{\mathrm{\, km\, s^{-1}}}{}$|⁠) and, tentatively, for J1810 (⁠|$K_\mathrm{Balmer} = 448 \pm 19{\mathrm{\, km\, s^{-1}}}{}$| and |$K_\mathrm{metals} = 491 \pm 32{\mathrm{\, km\, s^{-1}}}{}$|⁠). With existing inclination constraints, this gives the neutron star (NS) and companion masses |$M_\mathrm{NS} = 1.50$| – |$2.04{\, {\rm M}_{\odot }}{}$| and |$M_2 = 0.32$| – |$0.40{\, {\rm M}_{\odot }}{}$| for J1048, and |$M_\mathrm{NS} \gt 1.7{\, {\rm M}_{\odot }}{}$| and |$M_2 = 0.05$| – |$0.08{\, {\rm M}_{\odot }}{}$| for J1810. For J1908, we find an upper limit of |$K_2 \lt 32{\mathrm{\, km\, s^{-1}}}{}$|⁠ , which constrains its mass ratio |$q = M_2 / M_\mathrm{NS} \gt 0.55$| and inclination |$i \lt 6.0^\circ$|⁠ , revealing the previously misunderstood system to be the highest mass ratio, lowest inclination redback yet. This raises questions for the origins of its substantial radio eclipses. Additionally, we find evidence of asymmetric heating in J1048 and J1810, and signs of metal enrichment in J1908. We also explore the impact of inclination on spectroscopic temperatures, and demonstrate that the temperature measured at quadrature (⁠|$\phi = 0.25\,\mathrm{ and}\,0.75$|⁠) is essentially independent of inclination, and thus can provide additional constraints on photometric modelling. [ABSTRACT FROM AUTHOR]

Mature neutron stars are expected to exhibit gravity g modes due to stratification caused by a varying matter composition in the high-density core. By employing the BSk equation-of-state family, and working within the relativistic Cowling approximation, we examine how subtle differences in the nuclear matter assumptions impact on the g-mode spectrum. We investigate the possibility of detecting individual g-mode resonances during a binary inspiral with current and next-generation ground-based detectors, like Cosmic Explorer and the Einstein Telescope. Our results suggest that these resonances may be within the reach of future detectors, especially for low-mass stars with |$M\lesssim 1.4\,\mathrm{ M}_\odot$|⁠. [ABSTRACT FROM AUTHOR]

IGR J17591–2342, a transient accretion-powered millisecond X-ray pulsar, was discovered during its 2018 outburst. Here, we present a timing and spectral analysis of the source using AstroSat data of the same outburst. From the timing analysis, we obtain updated values of binary orbital parameters, which reveal an average pulsar spin frequency of 527.425 6984(8) Hz. The pulse profiles can be fit well with four harmonically related sinusoidal components with fractional amplitudes of fundamental and second, third, and fourth harmonics as ~13 per cent, ~6 per cent, ~0.9 per cent, and ~0.2 per cent, respectively. The energy-dependent study of pulse profiles in the range of 3–20 keV shows that the fractional amplitude of both the fundamental and first overtone is consistent with being constant across the considered energy band. Besides, a decaying trend has been observed for both the fundamental and first overtone in the phase-delay versus energy relation, resulting in soft X-ray (2.8–3.3 keV) phase lags of |$\sim$| 0.05 and |$\sim$| 0.13 with respect to |$\le 15$|  keV photons, for the fundamental and first overtone, respectively. The combined spectra from the Large Area X-ray Proportional Counters and the Soft X-ray Telescope aboard AstroSat in the 1–18 keV range can be fit well with an absorbed model consisting of a Comptonization, a blackbody and a Gaussian emission-line component yielding as best-fitting parameters a blackbody seed photon temperature |$kT_{\rm bb}$| |$\sim 0.95 \pm 0.03$| keV, and an electron temperature |$kT_{\rm e}$| |$\sim 1.54 \pm 0.03$| keV. The spectral aspects suggest the scattering of photons from the accretion disc or the neutron star's surface. [ABSTRACT FROM AUTHOR]

Transitional millisecond pulsars (tMSPs) are neutron star systems that alternate between a rotation-powered radio millisecond pulsar state and an accretion disc-dominated low-mass X-ray binary (LMXB)-like state on multi-year time-scales. During the LMXB-like state, the X-ray emission from tMSPs switches between 'low' and 'high' X-ray brightness modes on a time-scale of seconds to minutes (or longer), while the radio emission shows variability on time-scales of roughly minutes. Coordinated Very Large Array (VLA) and Chandra observations of the nearby tMSP PSR J1023+0038 uncovered a clear anticorrelation between radio and X-ray luminosities such that the radio emission consistently peaks during the X-ray low modes. In addition, there are sometimes also radio/X-ray flares that show no obvious correlation. In this paper, we present simultaneous radio and X-ray observations of a promising tMSP candidate system, 3FGL J1544.6 |$-$| 1125, which shows optical, |$\gamma$| -ray, and X-ray phenomena similar to PSR J1023+0038, but which is challenging to study because of its greater distance. Using simultaneous VLA and Chandra observations, we find that the radio and X-ray emission are consistent with being anticorrelated in a manner similar to PSR J1023+0038. We discuss how our results help in understanding the origin of bright radio emission from tMSPs. The greater sensitivity of upcoming telescopes such as the Square Kilometre Array will be crucial for studying the correlated radio/X-ray phenomena of tMSP systems. [ABSTRACT FROM AUTHOR]

A theoretical framework for calculating the mass–radius curve of compact stars in the Kaluza–Klein space–time is introduced, with one additional compact spatial dimension. Static, spherically symmetric solutions are considered, with the equation of state provided by a zero temperature, interacting multidimensional Fermi gas. To model the strong force between baryons, a repulsive potential is introduced, which is linear in the particle number density. The maximal mass of compact stars is calculated for different model parameters, and with a physical parameter choice, it satisfies observational data, meaning that it is possible to model simple, realistic objects within this framework. Based on this comparison, a limiting size for the observational regime of extra dimensions in compact stars is provided, with |$r_\mathrm{c} \gtrsim 0.2$|  fm. [ABSTRACT FROM AUTHOR]

We present results from three-dimensional, magnetohydrodynamic, core-collapse simulations of 16 progenitors following until 0.5 s after bounce. We use non-rotating solar-metallicity progenitor models with zero-age main-sequence mass between 9 and 24 |${\rm M}_{\odot }$|⁠. The examined progenitors cover a wide range of the compactness parameter including a peak around |$23 \, {\rm M}_{\odot }$|⁠. We find that neutrino-driven explosions occur for all models within 0.3 s after bounce. We also find that the properties of the explosions and the central remnants are well correlated with the compactness. Early shock evolution is sensitive to the mass accretion rate on to the central core, reflecting the density profile of the progenitor stars. The most powerful explosions with diagnostic explosion energy |$E_{\rm dia} \sim 0.75 \times 10^{51}$| erg are obtained by 23 and 24 |${\rm M}_{\odot }$| models, which have the highest compactness among the examined models. These two models exhibit spiral standing-accretion-shock-instability motions during 150–230 ms after bounce preceding a runaway shock expansion and leave a rapidly rotating neutron star with spin periods |$\sim 50$| ms. Our models predict the gravitational masses of the neutron star ranging between |$1.22$| and |$1.67 {\rm M}_{\odot }$| and their spin periods 0.04 – 4 s. The number distribution of these values roughly matches observation. On the other hand, our models predict small hydrodynamic kick velocity (15–260  |${\rm km \, s}^{-1}$|⁠), although they are still growing at the end of our simulations. Further systematic studies, including rotation and binary effects, as well as long-term simulations up to several seconds, will enable us to explore the origin of various core-collapse supernova explosions. [ABSTRACT FROM AUTHOR]

A fraction of the ultraluminous X-ray (ULX) sources are known to be accreting neutron stars as they show coherent X-ray pulsations with pulse periods ranging from |${\sim} 1{\!-\!}30$|  s. While initially thought to host intermediate-mass black holes, ULXs have since been recognized as a diverse class of objects, including ULX pulsars. These pulsars require models specifically tailored to account for their unique accretion physics, distinct from those used for Galactic black hole binaries. The X-ray spectra of all Galactic accreting X-ray pulsars (including sources in the Magellanic Clouds) are dominated by a high-energy cut-off power law and some of the sources show a soft excess, some emission lines, cyclotron absorption features, etc. In this work, we undertake a comprehensive analysis of the broad-band X-ray spectra of five ULX pulsars using simultaneous XMM – Newton and NuSTAR observations and show that their X-ray spectra can be effectively described by spectral models, similar to those used for the local accretion-powered X-ray pulsars. A soft excess is detected in all the sources which is also consistent with the local X-ray pulsars that have low absorption column density. We have marginal detection or low upper limit on the presence of the iron K-alpha emission line from these sources, which is a key difference of the ULX pulsars with the local accreting X-ray pulsars. We discuss the implication of this on the nature of the binary companion and the accretion mechanism in the ULX pulsars. [ABSTRACT FROM AUTHOR]

We study long-term evolution of the matter ejected in a black hole neutron star (BH–NS) merger employing the results of a long-term numerical-relativity simulation and nucleosynthesis calculation, in which both dynamical and post-merger ejecta formation is consistently followed. In particular, we employ the results for the merger of a |$1.35\, \mathrm{ M}_\odot$| NS and a |$5.4\, \mathrm{ M}_\odot$| BH with the dimensionless spin of 0.75. We confirm the finding in the previous studies that thermal pressure induced by radioactive heating in the ejecta significantly modifies the morphology of the ejecta. We then compute the kilonova (KN) light curves employing the ejecta profile obtained by the long-term evolution. We find that our present BH–NS model results in a KN light curve that is fainter yet more enduring than that observed in AT2017gfo. This is due to the fact that the emission is primarily powered by the lanthanide-rich dynamical ejecta, in which a long photon diffusion time-scale is realized by the large mass and high opacity. While the peak brightness of the KN emission in both the optical and near-infrared bands is fainter than or comparable to those of binary NS models, the time-scale maintaining the peak brightness is much longer in the near-infrared band for the BH–NS KN model. Our result indicates that a BH–NS merger with massive ejecta can observationally be identified by the long lasting (> two weeks) near-infrared emission. [ABSTRACT FROM AUTHOR]

'Z' type neutron star low-mass X-ray binaries typically show a 'Z'-like three-branched track in their hardness intensity diagram. However, a few such 'Z' sources show an additional branch known as the extended flaring branch (EFB). EFB has been poorly studied, and its origin is not known. It is thought to be an extension of the flaring branch (FB) or associated with Fe K  |$\alpha$| complex or an additional continuum due to the radiative recombination continuum (RRC) process. Using AstroSat observations, we have detected the EFB from two 'Z' sources, GX 340+0 and GX 5–1, and performed a broad-band spectral analysis in the 0.5–22 keV energy range. During EFB, both sources show the presence of a significant RRC component with absorption edges at |$7.91^{+0.16}_{-0.15}$| and |$8.10^{+0.16}_{-0.17}$|  keV, respectively along with blackbody radiation and thermal Comptonization. No signature of RRC was detected during the FB, which is adjoint to the EFB. No Fe K  |$\alpha$| complex is detected. Interestingly, inside EFB dips of GX 5-1, for the first time, we have detected flaring events of 30–60 s, which can be modelled with a single blackbody radiation. During the FB to EFB transition, an increase in the blackbody radius by a factor of 1.5–2 is observed in both sources. Our analysis strongly suggests that EFB is not an extension of FB or caused by the Fe K  |$\alpha$| complex. Rather, it is caused by a sudden expansion of the hot, thermalized boundary layer and subsequent rapid cooling. [ABSTRACT FROM AUTHOR]

We present the first model aimed at understanding how the Meissner effect in a young neutron star affects its macroscopic magnetic field. In this model, field expulsion occurs on a dynamical time-scale, and is realized through two processes that occur at the onset of superconductivity: fluid motions causing the dragging of field lines, followed by magnetic reconnection. Focusing on magnetic fields weaker than the superconducting critical field, we show that complete Meissner expulsion is but one of four possible generic scenarios for the magnetic-field geometry, and can never expel magnetic flux from the centre of the star. Reconnection causes the release of up to |$\sim 5\times 10^{46}\, \mathrm{erg}$| of energy at the onset of superconductivity, and is only possible for certain favourable early-phase dynamics and for pre-condensation fields |$10^{12}\, \mathrm{G}\lesssim B\lesssim 5\times 10^{14}\, \mathrm{G}$|⁠. Fields weaker or stronger than this are predicted to thread the whole star. [ABSTRACT FROM AUTHOR]

Evidence has accumulated for an as-yet unaccounted for source of heat located at shallow depths within the accreted neutron star crust. However, the nature of this heat source is unknown. I demonstrate that the inferred depth of carbon ignition in X-ray superbursts can be used as an additional constraint for the magnitude and depth of shallow heating. The inferred shallow heating properties are relatively insensitive to the assumed crust composition and carbon fusion reaction rate. For low-accretion rates, the results are weakly dependent on the duration of the accretion outburst, so long as accretion has ensued for enough time to replace the ocean down to the superburst ignition depth. For accretion rates at the Eddington rate, results show a stronger dependence on the outburst duration. Consistent with earlier work, it is shown that urca cooling does not impact the calculated superburst ignition depth unless there is some proximity in depth between the heating and cooling sources. [ABSTRACT FROM AUTHOR]

This paper examines the oscillatory behaviour of relativistic, non-self-gravitating, charged-fluid toroidal structures within the context of the Kerr metric. The primary objective is to explore how thick accretion discs influence the mass and spin measurements of black holes and neutron stars in low-mass X-ray binaries (LMXBs) through quasi-periodic oscillations (QPOs) models. To achieve this, we conduct a local analysis within a general relativistic framework, determining the radial epicyclic and orbital frequencies in a perfect fluid disc. The tori are modelled using a non-Keplerian distribution of specific angular momentum, and we analyse how the oscillation properties depend on the model's angular momentum distribution parameters. Subsequently, we connect these oscillatory frequencies to high-frequency QPOs observed in LMXBs, enabling us to calculate the optimal mass and spin values for each studied source. [ABSTRACT FROM AUTHOR]

Several sources of repeating coherent bursts of radio emission with periods of many minutes have now been reported in the literature. These 'ultralong period' (ULP) sources have no clear multiwavelength counterparts and challenge canonical pulsar emission models, leading to debate regarding their nature. In this work, we report the discovery of a bright, highly polarized burst of radio emission at low Galactic latitude as part of a wide-field survey for transient and variable radio sources. ASKAP J175534.9 |$-$| 252749.1 does not appear to repeat, with only a single intense two-minute |$\sim$| 200-mJy burst detected from 60 h of observations. The burst morphology and polarization properties are comparable to those of classical pulsars but the duration is more than one hundred times longer, analogous to ULPs. Combined with the existing ULP population, this suggests that these sources have a strong Galactic latitude dependence and hints at an unexplored population of transient and variable radio sources in the thin disc of the Milky Way. The resemblance of this burst with both ULPs and pulsars calls for a unified coherent emission model for objects with spin periods from milliseconds to tens of minutes. However, whether or not these are all neutron stars or have the same underlying power source remains open for debate. [ABSTRACT FROM AUTHOR]

A neutron star (NS) accreting matter from a companion star in a low-mass X-ray binary (LMXB) system can spin up to become a millisecond pulsar (MSP). Properties of many such MSP systems are known, which is excellent for probing fundamental aspects of NS physics when modelled using the theoretical computation of NS LMXB evolution. Here, we systematically compute the long-term evolution of NS, binary, and companion parameters for NS LMXBs using the stellar evolution code mesa. We consider the baryonic to gravitational mass conversion to calculate the NS mass evolution and show its cruciality for the realistic computation of some parameters. With computations using many combinations of parameter values, we find the general nature of the complex NS spin frequency (⁠|$\nu$|⁠) evolution, which depends on various parameters, including accretion rate, fractional mass-loss from the system, and companion star magnetic braking. Further, we utilize our results to precisely match some main observed parameters, such as |$\nu$|⁠ , orbital period (⁠|$P_{\rm orb}$|⁠), etc. of four accreting millisecond X-ray pulsars (AMXPs). By providing the |$\nu$|⁠ , |$P_{\rm orb}$|⁠ , and the companion mass spaces for NS LMXB evolution, we indicate the distribution and plausible evolution of a few other AMXPs. We also discuss the current challenges in explaining the parameters of AMXP sources with brown dwarf companions and indicate the importance of modelling the transient accretion in LMXBs as a possible solution. [ABSTRACT FROM AUTHOR]

The inner and outer cores of neutron stars are believed to contain type-I and -II proton superconductors, respectively. The type-I superconductor exists in an intermediate state, comprising macroscopic flux-free and flux-containing regions, while the type-II superconductor is flux-free, except for microscopic, quantized flux tubes. Here, we show that the inner and outer cores are coupled magnetically, when the macroscopic flux tubes subdivide dendritically into quantized flux tubes, a phenomenon called flux branching. An important implication is that up to |${\sim} 10^{12} (r_1/10^6 \, {\rm cm}) \, {\rm erg}$| of energy are required to separate a quantized flux tube from its progenitor macroscopic flux tube, where |$r_1$| is the length of the macroscopic flux tube. Approximating the normal-superconducting boundary as sharp, we calculate the magnetic coupling energy between a quantized and macroscopic flux tube due to flux branching as a function of, |$f_1$|⁠ , the radius of the type-I inner core divided by the radius of the type-II outer core. Strong coupling delays magnetic field decay in the type-II superconductor. For an idealized inner core containing only a type-I proton superconductor and poloidal flux, and in the absence of ambipolar diffusion and diamagnetic screening, the low magnetic moments (⁠|${\lesssim} 10^{27} \, {\rm G \, cm^3}$|⁠) of recycled pulsars imply |$f_1 \lesssim 10^{-1.5}$|⁠. [ABSTRACT FROM AUTHOR]

We present a formulation of special relativistic dissipative hydrodynamics (SRDHD) derived from the well-established Müller–Israel–Stewart (MIS) formalism using an expansion in deviations from ideal behaviour. By re-summing the non-ideal terms, our approach extends the Euler equations of motion for an ideal fluid through a series of additional source terms that capture the effects of bulk viscosity, shear viscosity, and heat flux. For efficiency these additional terms are built from purely spatial derivatives of the primitive fluid variables. The series expansion is parametrized by the dissipation strength and time-scale coefficients, and is therefore rapidly convergent near the ideal limit. We show, using numerical simulations, that our model reproduces the dissipative fluid behaviour of other formulations. As our formulation is designed to avoid the numerical stiffness issues that arise in the traditional MIS formalism for fast relaxation time-scales, it is roughly an order of magnitude faster than standard methods near the ideal limit. [ABSTRACT FROM AUTHOR]

PSR J0737−3039A/B is unique among double neutron star systems. Its near-perfect edge-on orbit causes the fast spinning pulsar A to be eclipsed by the magnetic field of the slow spinning pulsar B. Using high-sensitivity MeerKAT radio observations combined with updated constraints on the system geometry, we studied the impact of these eclipses on the incident polarization properties of pulsar A. Averaging light curves together after correcting for the rotation of pulsar B revealed enormous amounts of circular polarization and rapid changes in the linear polarization position angle, which occur at phases where emission from pulsar A is partially transmitted through the magnetosphere of pulsar B. These behaviours confirm that the eclipse mechanism is the result of synchrotron absorption in a relativistic pair-plasma confined to the closed-field region of pulsar B's truncated dipolar magnetic field. We demonstrate that changes in circular polarization handedness throughout the eclipses are directly tied to the average line of sight magnetic field direction of pulsar B, from which we unambiguously determine the complete magnetic and viewing geometry of the pulsar. [ABSTRACT FROM AUTHOR]

Tellurium is a primary candidate for the identification of the 2.1  |$\, \mu$| m emission line in kilonovae (KNe) spectra AT2017gfo and GRB230307A. Despite this, there is currently an insufficient amount of atomic data available for this species. We calculate the required atomic structure and collisional data, particularly the data required for accurate non-local-thermodynamic-equilibrium (NLTE) modelling of the low temperatures and densities in KNe. We use a multiconfigurational Dirac–Hartree–Fock method to produce optimized one-electron orbitals for Te i - iii. As a result energy levels and Einstein A-coefficients for Te i - iii have been calculated. These orbitals are then employed within Dirac R -matrix collision calculations to provide electron-impact-excitation collision strengths that were subsequently averaged according to a thermal Maxwellian distribution. Subsequent tardis simulations using this new atomic data reveal no significant changes to the synthetic spectra due to the very minor contribution of Te at early epochs. NLTE simulations with the colradpy package reveal optically thin spectra consistent with the increasing prominence of the Te iii 2.1 |$\, \mu$| m line as the KNe ejecta cools. This is reinforced by the estimation of luminosities at nebular KNe conditions. New line ratios for both observation and laboratory benchmarks of the atomic data are proposed. [ABSTRACT FROM AUTHOR]

Strongly magnetized neutron stars are popular candidates for producing detectable electromagnetic and gravitational-wave signals. Gravitational collapses of neutron stars triggered by a phase transition from hadrons to deconfined quarks in the cores could also release a considerable amount of energy in the form of gravitational waves and neutrinos. Hence, the formation of a magnetized hybrid star from such a phase-transition-induced collapse is an interesting scenario for detecting all these signals. These detections may provide essential probes for the magnetic field and composition of such stars. Thus far, a dynamical study of the formation of a magnetized hybrid star from a phase-transition-induced collapse has yet to be realized. Here, we investigate the formation of a magnetized hybrid star with a purely toroidal field and its properties through dynamical simulations. We find that the maximum values of rest-mass density and magnetic field strength increase slightly and these two quantities are coupled in phase during the formation. We then demonstrate that all microscopic and macroscopic quantities of the resulting hybrid star vary drastically when the maximum magnetic field strength goes beyond a threshold of |$\sim 5 \times 10^{17}$| G, but they are insensitive to the magnetic field below this threshold. Specifically, the magnetic deformation makes the rest-mass density drop significantly, suppressing the matter fraction in the mixed phase. These behaviours agree with those in the equilibrium models of previous studies. Therefore, this work provides a solid support for the magnetic effects on a hybrid star. [ABSTRACT FROM AUTHOR]

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-h multibeam 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 signal-to-noise ratio being strongest in the wide-field incoherent beam of MeerKAT. [ABSTRACT FROM AUTHOR]

Comprehensive spectrotemporal analyses of the Z-type neutron star low-mass X-ray binary GX 5−1 were performed using 10 broad-band observations from AstroSat /Soft X-ray Telescope and Large Area X-ray Proportional Counter (LAXPC) instruments. The LAXPC-20 hardness–intensity diagram showed horizontal and normal branches (HBs and NBs) of the Z track which exhibited secular motion. The time-averaged spectra in the energy range 0.7–25.0 keV could be fitted with the model combination – |$\tt {constant}\, \times \, \tt {tbabs}\, \times \, \tt {edge}\, \times \, \tt {edge}\, \times \, \tt {thcomp}\, \times \, \tt {diskbb}$|⁠. This yielded |$\Gamma \, \sim$| 2, |$kT_{\mathrm{ e}}\, \sim$| 3.3 keV, and |$F_{\mathrm{ disc}}$| / |$F_{\mathrm{ total}}\, \sim$| 0.8 indicating the soft/intermediate spectral state of the source during the observations. Flux-resolved spectral analysis revealed a positive correlation between |$kT_{\mathrm{ in}}$| and |$F_{\mathrm{ bol}}$|⁠. However, a negative correlation was observed between them in one of the NBs. Time-averaged temporal analysis revealed multiple HB oscillations (HBOs) and NB oscillations (NBOs), and peaked noise components in the |$\sim$| 5–50 Hz range. Furthermore, flux-resolved temporal analysis showed that the frequency of the HBOs correlates positively whereas the strength of HBOs correlates negatively with |$F_{\mathrm{ bol}}$|⁠ , indicating their probable origin from the accretion disc. In contrast, the frequency and strength of NBOs remain fairly constant with |$F_{\mathrm{ bol}}$|⁠ , suggesting that they originate from a different region in the system. Using the relativistic precession model along with highest frequency of the HBO, the upper limits of the magnetic dipole moment (⁠|$\mu$|⁠) and field strength (B) at the poles of the neutron star in the system were found to be 25.60 |$\times \, 10^{25}$| G cm3 and 3.64 |$\times \, 10^{8}$| G, respectively, for |$k_{\mathrm{ A}}$| = 1. [ABSTRACT FROM AUTHOR]

In modelling a relativistic disc around a compact object, the self-gravity of the disc is often neglected while it needs to be incorporated for more accurate descriptions in several circumstances. Extending the Komatsu–Eriguchi–Hachisu self-consistent field method, we present numerical models of a rapidly rotating neutron star with a self-gravitating disc in stationary equilibrium. In particular, our approach allows us to obtain numerical solutions involving a massive disc with the rest mass |$\mathcal {O}(10^{-1})-\mathcal {O}(10^0)\, \mathrm{ M}_\odot$| closely attached to a rotating neutron star, given that the disc is mainly supported by the relativistic electron degeneracy pressure. We also assess the impact of self-gravity on the internal structure of the disc and the neutron star. These axisymmetric, stationary solutions can be employed for simulations involving the neutron star–disc system in the context of high-energy transients and gravitational-wave emissions. [ABSTRACT FROM AUTHOR]

We present a comprehensive study of the spectral properties of the accreting millisecond X-ray pulsar IGR J17498−2921 during its 2023 outburst. Similar to other accreting millisecond X-ray pulsars, the broadband spectral emission observed quasi-simultaneously by NICER and NuSTAR is well described by an absorbed Comptonized emission with an electron temperature of ∼17 keV plus a disk reflection component. The broadening of the disk reflection spectral features, such as a prominent iron emission line at 6.4–6.7 keV, is consistent with the relativistic motion of matter in a disk truncated at ∼21 Rg from the source, near the Keplerian corotation radius. From the high-cadence monitoring data obtained with NICER, we observed that the evolution of the photon index and the temperature of seed photons tracks variations in the X-ray flux. This is particularly evident close to a sudden ∼–0.25 cycle jump in the pulse phase, which occurs immediately following an X-ray flux flare and a drop in the pulse amplitude below the 3σ detection threshold. We also report on the non-detection of optical pulsations with TNG/SiFAP2 from the highly absorbed optical counterpart. [ABSTRACT FROM AUTHOR]

Context. On April 25, 2019, the LIGO-Virgo Collaboration discovered a gravitational-wave (GW) signal from a binary neutron star (BNS) merger, that is, GW190425. Due to the inferred large total mass, the origin of GW190425 remains unclear. Aims. Assuming GW190425 originated from the standard isolated binary evolution channel, its immediate progenitor is considered to be a close binary system, consisting of a He-rich star and a NS just after the common envelope phase. We aim to study the formation of GW190425 in a solar-like environment by using the detailed binary evolution code MESA. Methods. We perform detailed stellar structure and binary evolution calculations that take into account mass loss, internal differential rotation, and tidal interactions between a He-rich star and a NS companion. We explore the parameter space of the initial binary properties, including initial NS and He-rich masses and initial orbital period. Results. We find that the immediate post-common-envelope progenitor system, consisting of a primary ∼2.0 M⊙ (∼1.7 M⊙) NS and a secondary He-rich star with an initial mass of ∼3.0 − 5.5 M⊙ (∼5.5 − 6.0 M⊙) in a close binary with an initial period of ∼0.08 − 0.5 days (∼0.08 − 0.4 days), that experiences stable Case BB/BC mass transfer (MT) during binary evolution, can reproduce the formation of GW190425-like BNS events. Our studies reveal that the secondary He-rich star of the GW190425's progenitor before its core collapse can be efficiently spun up through tidal interaction, finally remaining as a NS with rotational energy even reaching ∼1052 erg, which is always much higher than the neutrino-driven energy of the supernova (SN) explosion. If the newborn secondary NS is a magnetar, we expect that GW190425 can be the remnant of a magnetar-driven SN, namely a magnetar-driven ultra-stripped SN, a superluminous SN, or a broad-line Type Ic SN. Conclusions. Our results show that GW190425 could be formed through the isolated binary evolution, which involves a stable Case BB/BC MT just after the common envelope phase. On top of that, we show the He-rich star can be tidally spun up, potentially forming a spinning magnetized NS (magnetar) during the second SN explosion. [ABSTRACT FROM AUTHOR]

We report on the broadband spectral and timing properties of the accreting millisecond X-ray pulsar IGR J17498−2921 during its April 2023 outburst. We used data from NICER (1–10 keV), NuSTAR (3–79 keV), Insight-HXMT (2–150 keV), and INTEGRAL (30–150 keV). We detected significant 401 Hz pulsations across the 0.5–150 keV band. The pulse fraction increases from ∼2% at 1 keV to ∼13% at 66 keV. We detected five type-I X-ray bursts, including three photospheric radius expansion bursts, with a rise time of ∼2 s and an exponential decay time of ∼5 s. The recurrence time is ∼9.1 h, which can be explained by unstable thermonuclear burning of hydrogen-deficient material on the neutron star surface. The quasi-simultaneous 1–150 keV broadband spectra from NICER, NuSTAR and INTEGRAL can be aptly fitted by an absorbed reflection model, relxillCp, and a Gaussian line of instrumental origin. The Comptonized emission from the hot corona is characterized by a photon index Γ of ∼1.8 and an electron temperature kTe of ∼40 keV. We obtained a low inclination angle i ∼ 34°. The accretion disk shows properties of strong ionization, log(ξ/erg cm s−1)∼4.5, over-solar abundance, AFe ∼ 7.7, and high density, log(ne/cm−3)∼19.5. However, a lower disk density with normal abundance and ionization could also be possible. Based on the inner disk radius of Rin = 1.67RISCO and the long-term spin-down rate of −3.1(2)×10−15 Hz s−1, we were able to constrain the magnetic field of IGR J17498−2921 to the range of (0.9 − 2.4)×108 G. [ABSTRACT FROM AUTHOR]

We report a detailed timing and spectral analysis of the X-ray pulsar 2S 1417−624 using the data from Insight-HXMT during the 2018 outburst. The pulse profiles are highly variable with respect to both unabsorbed flux and energy. A double-peaked pulse profile from the low flux evolved to a multi-peaked shape in the high-flux state. The pulse fraction is negatively correlated to the source flux in the range of ∼(1 − 6)×10−9 erg cm−2 s−1, consistent with Rossi X-ray Timing Explorer (RXTE) studies during the 2009 giant outburst. The energy-resolved pulse profiles around the peak outburst showed a four-peak shape in the low-energy bands and gradually evolved to triple peaks at higher energies. The continuum spectrum is well described by typical phenomenological models, such as the cut-off power law and the power law with high-energy cut-off models. Notably, we discovered high-energy cyclotron resonant scattering features (CRSFs) for the first time, which are around 100 keV with a statistical significance of ∼7σ near the peak luminosity of the outburst. This CRSF line is significantly detected with different continuum models and provides very robust evidence for its presence. Furthermore, pulse-phase-resolved spectroscopy confirmed the presence of the line, whose energy varied from 97 to 107 keV over the pulse phase and appeared to have a maximum value at the narrow peak phase of the profiles. [ABSTRACT FROM AUTHOR]

We present the first X-ray observation of the energetic millisecond pulsar binary PSR J1431−4715, performed with XMM-Newton and complemented with fast optical multi-band photometry acquired with the ULTRACAM instrument at ESO-NTT. It is found as a faint X-ray source without a significant orbital modulation. This contrasts with the majority of systems that instead display substantial X-ray orbital variability. The X-ray spectrum is dominated by non-thermal emission and, due to the lack of orbital modulation, does not favour an origin in an intrabinary shock between the pulsar and companion star wind. While thermal emission from the neutron star polar cap cannot be excluded in the soft X-rays, the dominance of synchrotron emission favours an origin in the pulsar magnetosphere that we describe at both X-ray and gamma-ray energies with a synchro-curvature model. The optical multi-colour light curve folded at the 10.8 h orbital period is double-humped and dominated by ellipsoidal effects, but also affected by irradiation. The ULTRACAM light curves are fit with several models encompassing direct heating and a cold spot, or heat redistribution after irradiation either through convection or convection plus diffusion. Despite the inability to constrain the best irradiation models, the fits provide consistent system parameters, giving an orbital inclination of 59 ± 6° and a distance of 3.1 ± 0.3 kpc. The companion is found to be an F-type star, underfilling its Roche lobe (fRL = 73 ± 4%) with a mass of 0.20 ± 0.04 M⊙, confirming the redback status, but hotter than the majority of redbacks. The stellar dayside and nightside temperatures of 7500 K and 7400 K, respectively, indicate a weak irradiation effect on the companion, likely due to its high intrinsic luminosity. Although the pulsar mass cannot be precisely derived, a heavy (1.8−2.2 M⊙) neutron star is favoured. [ABSTRACT FROM AUTHOR]

Context. Pair cascades in polar cap regions of neutron stars are considered to be an essential process in various models of coherent radio emissions of pulsars. The cascades produce pair plasma bunch discharges in quasi-periodic spark events. The cascade properties, and therefore also the coherent radiation, depend strongly on the magnetospheric plasma properties and vary significantly across and along the polar cap. Importantly, where the radio emission emanates from in the polar cap region is still uncertain. Aims. We investigate the generation of electromagnetic waves by pair cascades and their propagation in the polar cap for three representative inclination angles of a magnetic dipole, 0°, 45°, and 90°. Methods. We use two-dimensional particle-in-cell simulations that include quantum-electrodynamic pair cascades in a charge-limited flow from the star surface. Results. We find that the discharge properties are strongly dependent on the magnetospheric current profile in the polar cap and that transport channels for high intensity Poynting flux are formed along magnetic field lines where the magnetospheric currents approach zero and where the plasma cannot carry the magnetospheric currents. There, the parallel Poynting flux component is efficiently transported away from the star and may eventually escape the magnetosphere as coherent radio waves. The Poynting flux decreases with increasing distance from the star in regions of high magnetospheric currents. Conclusions. Our model shows that no process of energy conversion from particles to waves is necessary for the coherent radio wave emission. Moreover, the pulsar radio beam does not have a cone structure; rather, the radiation generated by the oscillating electric gap fields directly escapes along open magnetic field lines in which no pair creation occurs. [ABSTRACT FROM AUTHOR]

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 × 1038 erg s−1 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 ± 0.8% and 97° ±8° for Observation 1, 3.0 ± 0.9% and 90° ±8° for Observation 2, and 5.5 ± 1.1% and 80° ±6° 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°. 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 ∼70° to ∼100°. 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. [ABSTRACT FROM AUTHOR]

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 several 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 Imaging X-ray Polarimetry Explorer (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 in 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 a 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 ip = 15°–40°, the magnetic obliquity of θp = 60°–80°, and the pulsar spin position angle of χp ≈ −50°, which significantly differs from the constant component PA of about 10°. Combining these X-ray measurements with the optical PA, we find evidence for at least a 30° misalignment between the pulsar angular momentum and the binary orbital axis. [ABSTRACT FROM AUTHOR]

The four characteristic oscillation frequencies of accretion flows (in addition to the Keplerian orbital frequency) are often discussed in the context of the time variability of black hole and neutron star (NS) low-mass X-ray binaries (LMXBs). These four frequencies are the frequencies of the axisymmetric radial and vertical epicyclic oscillations, and the frequencies of non-axisymmetric oscillations corresponding to the periastron (radial) and Lense-Thirring (vertical) precessions. In this context, we investigated the effect of the quadrupole moment of a slowly rotating NS and provide complete formulae for calculating these oscillation and precession frequencies, as well as convenient approximations. Simple formulae corresponding to the geodesic limit of a slender torus (and test-particle motion) and the limit of a marginally overflowing torus (a torus exhibiting a critical cusp) are presented, and more general approximate formulae are included to allow calculations for arbitrarily thick tori. We provide the Wolfram Mathematica code used for our calculations together with the C++ and PYTHON codes for calculating the frequencies. Our formulae can be used for various calculations regarding the astrophysical signatures of the NS super-dense matter equation of state. For instance, we demonstrate that even for a given fixed number of free parameters, a model that accounts for fluid flow precession matches the frequencies of twin-peak quasiperiodic oscillations observed in NS LMXBs better than a model that uses geodesic precession. [ABSTRACT FROM AUTHOR]

We present a full general relativistic analytic solution for a radiation-pressure-supported equilibrium fluid torus orbiting a rotating neutron star (NS). We applied previously developed analytical methods that include the effects of both the NS's angular momentum and quadrupole moment in the Hartle-Thorne geometry. The structure, size, and shape of the torus are explored, with a particular focus on the critically thick solution – the cusp tori. For the astrophysically relevant range of NS parameters, we examined how our findings differ from those obtained for the Schwarzschild space-time. The solutions for rotating stars display signatures of an interplay between relativistic and Newtonian effects where the impact of the NS angular momentum and quadrupole moment are almost counterbalanced at a given radius. Nevertheless, the space-time parameters still strongly influence the size of tori, which can be shown in a coordinate-independent way. Finally, we discuss the importance of the size of the central NS which determines whether or not a surrounding torus exists. We provide a set of tools in a Wolfram Mathematica code, which establishes a basis for further investigation of the impact of the NSs' super-dense matter equation of state on the spectral and temporal behaviour of accretion tori. [ABSTRACT FROM AUTHOR]

Context. The connection between quasi-periodic oscillations (QPOs) and magnetic fields has been investigated in various celestial bodies. Magnetohydrodynamic (MHD) waves have been employed to explain the simultaneous upper and lower kilohertz (kHz) QPOs. Nevertheless, the intricate and undefined formation pathways of twin kHz QPOs present a compelling avenue for exploration. This study area holds great interest as it provides an opportunity for deriving crucial parameters related to compact stars. Aims. We strive to develop a self-consistent model elucidating the radiation mechanism of twin kHz QPOs, which we then compare it with observations. Methods. A sample of 28 twin kHz QPOs detected from the X-ray binary 4U 1636–53 was used for a comparison with the results of the Markov chain Monte Carlo calculations based on our model of the radiation mechanism of twin kHz QPOs, which is related to twin MHD waves. Results. We obtained 28 groups of parameters of 4U 1636–53 and a tight exponential fit between the flux and the temperature of seed photons to Compton up-scattering and find that the electron temperature in the corona around the neutron star decreases with increasing seed photon temperature. Conclusions. The origin of twin kHz QPOs are dual disturbances that arise from twin MHD waves that are generated at the innermost radius of an accretion disc. The seed photons can be transported through a high temperature corona and are Compton up-scattered. The photons that vary the frequencies of twin MHD waves lead to the observed twin kHz QPOs. [ABSTRACT FROM AUTHOR]

Context. 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 observations and predictions from standard binary evolution scenarios of pulsar-helium white dwarf (HeWD) binaries, especially in the case of PSR J0955−6150, for which mass measurements ruled out most eMSP formation models. Aims. We aim to measure the masses of the pulsar and its companion, and constrain the orbital configuration of PSR J1618−3921. This facilitates understanding similarities among eMSPs and could offer hints as to their formation mechanism. Methods. 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 Parkes and Nançay observations. We performed a full analysis on this joint data set with a timing baseline of 23 years. We also used the data from recent observations to give a brief account of the emission properties of J1618−3921, including a rotating vector model (RVM) fit of the linear polarisation position angle of the pulsar. Results. From the timing analysis, we measure a small but significant proper motion of the pulsar. The long timing baseline allowed for a highly significant measurement of the rate of advance of periastron of ω̇ = (0.00145±0.00010)°yr−1. Despite the tenfold improvement in timing precision from MeerKAT observations, we can only report a low-significance detection of the orthometric Shapiro delay parameters, h3 = 2.70−1.47+2.07 μs and ς = 0.68−0.09+0.13. Under the assumption of the validity of general relativity (GR), the self-consistent combination of these three parameters leads to mass estimates of the total and individual masses in the binary of Mtot = 1.42−0.19+0.20 M⊙, Mc = 0.20−0.03+0.11 M⊙, and Mp = 1.20−0.20+0.19 M⊙. We detect an unexpected change in the orbital period of Ṗb = −2.26−0.33+0.35 × 10−12, that is an order of magnitude larger and carries an opposite sign to what is expected from the Galactic acceleration and the Shklovskii effect, which are a priori the only non-negligible contributions expected for Ṗb. We also detect a significant second derivative of the spin frequency, f̈. The RVM fit reveals a viewing angle of ζ = (111 ± 1)°. Furthermore, we report an unexpected, abrupt change in the mean pulse profile in June 2021 of unknown origin. Conclusions. We propose that the anomalous Ṗb and f̈ we measure for J1618−3921 indicate an additional varying acceleration due to a nearby mass. The J1618−3921 binary system is likely part of a hierarchical triple, but with the third component much farther away than the outer component of the MSP in a triple star system, PSR J0337+1715. 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 Pb − MWD relation, which has been verified for circular HeWD binaries and also for the two HeWDs in the PSR J0337+1715 system. Future regular observations with the MeerKAT telescope will, due to the further extension of the timing baseline, improve the measurement of Ṗb and f̈. This will help us further understand the nature of this system, and perhaps improve our understanding of eMSPs in general. [ABSTRACT FROM AUTHOR]

This study utilizes the simultaneous broad-band observations of 4U 1735−44 from AstroSat , offering enhanced spectral and temporal resolution, to investigate its spectral properties, temporal behaviour, and burst characteristics. Spectral, type I X-ray burst, and temporal analyses on 4U 1735−44 were performed using AstroSat /Soft X-ray Telescope and Large Area X-ray Proportional Counter (LAXPC) observations. The hardness–intensity diagram from LAXPC-20 showed a positive correlation between hardness and intensity, with a pattern resembling the banana branch typical of atoll sources. Spectral analysis carried out in the 0.7–20.0 keV energy range, using the model combination – |$\tt {constant}$| |$\times$| |$\tt {tbabs}$| (⁠|$\tt {nthcomp}$| |$+$| |$\tt {diskbb}$| |$+$| |$\tt {bbodyrad}$|⁠), suggested a cool accretion disc truncated at a large distance from the neutron star in the system. Time-resolved spectral studies of two type I X-ray bursts detected from the source revealed evidence of photospheric radius expansion, allowing for an estimation of the source distance. Temporal analysis showed the presence of low-frequency quasi-periodic oscillation at |$\sim$| 69 Hz (3.3 |$\sigma$| significance with more than 99 per cent confidence) and prominent noise features below 30 Hz. [ABSTRACT FROM AUTHOR]

The recent observation of the GW230529 event indicates that black hole–neutron star binaries can contain low-mass black holes. Since lower mass systems are more favourable for tidal disruption, such events are promising candidates for multimessenger observations. In this study, we employ five finite-temperature, composition-dependent matter equations of state and present results from ten 3D general relativistic hydrodynamic simulations for the mass ratios |$q = 2.6$| and 5. Two of these simulations target the chirp mass and effective spin parameter of the GW230529 event, while the remaining eight contain slightly higher mass black holes, including both spinning (⁠|$a_{\mathrm{ BH}} = 0.7$|⁠) and non-spinning (⁠|$a_{\mathrm{ BH}} = 0$|⁠) models. We discuss the impact of the equation of state, spin, and mass ratio on black hole–neutron star mergers by examining both gravitational-wave and ejected matter properties. For the low-mass ratio model we do not see fast-moving ejecta for the softest equation of state model, but the stiffer model produces on the order of |$10^{-6}\,\mathrm{ M}_\odot$| of fast-moving ejecta, expected to contribute to an electromagnetic counterpart. Notably, the high-mass ratio model produces nearly the same amount of total dynamical ejecta, but yields 52 times more fast-moving ejecta than the low-mass ratio system. In addition, we observe that the black-hole spin tends to decrease the amount of fast-moving ejecta while increasing significantly the total ejected mass. Finally, we note that the disc mass tends to increase as the neutron star compactness decreases. [ABSTRACT FROM AUTHOR]

Magnetic field evolution of neutron stars is a long-standing debate. The rate of magnetic field decay for isolated, non-accreting neutron stars can be quantified by measuring the negative second derivative of the spin period. Alternatively, this rate can be estimated by observing an excess of thermal emission with respect to the standard cooling without additional heating mechanisms involved. One of the nearby cooling isolated neutron stars – RX J0720.4−3125, – offers a unique opportunity to probe the field decay as for this source there are independent measurements of the surface X-ray luminosity, the second spin period derivative, and magnetic field. We demonstrate that the evolution rate of the spin period derivative is in correspondence with the rate of dissipation of magnetic energy of the dipolar field if a significant part of the released energy is emitted in X-rays. The instantaneous time-scale for the magnetic field decay is |${\sim} 10^4$|  yr. [ABSTRACT FROM AUTHOR]

We have investigated the evolutionary connections of the isolated neutron star (NS) populations including radio pulsars (RPs), anomalous X-ray pulsars (AXPs), soft gamma repeaters (SGRs), dim isolated NSs (XDINs), 'high-magnetic field' RPs ('HBRPs'), central compact objects (CCOs), rotating radio transients (RRATs), and long-period pulsars (LPPs) in the fallback disc model. The model can reproduce these NS families as a natural outcome of different initial conditions (initial period, disc mass, and dipole moment, μ) with a continuous μ distribution in the |$\sim 10^{27} - 5 \times 10^{30}$| G cm |$^3$| range. Results of our simulations can be summarized as follows: (1) A fraction of 'HBRPs' with relatively high μ evolve into the persistent AXP/SGR properties, and subsequently become LPPs. (2) Persistent AXP/SGRs do not have evolutionary links with CCOs, XDINs, and RRATs. (3) For a wide range of μ , most RRATs evolve passing through RP or 'HBRP' properties during their early evolutionary phases. (4) A fraction of RRATs which have the highest estimated birth rate seem to be the progenitors of XDINs. (5) LPPs, whose existence was predicted by the fallback disc model, are the sources evolving in the late stage of evolution before the discs become inactive. These results provide concrete support to the ideas proposing evolutionary connections between the NS families to account for the 'birth rate problem', the discrepancy between the cumulative birth rate estimated for these systems and the core-collapse supernova rate. [ABSTRACT FROM AUTHOR]