Search

In the standard two-component crust-superfluid model of a neutron star, timing noise can arise when the two components are perturbed by stochastic torques. Here it is demonstrated how to analyse fluctuations in radio pulse times of arrival with a Kalman filter to measure physical properties of the two-component model, including the crust-superfluid coupling time-scale and the variances of the crust and superfluid torques. The analysis technique, validated previously on synthetic data, is applied to observations with the Molonglo Observatory Synthesis Telescope of the representative pulsar PSR J1359−6038. It is shown that the two-component model is preferred to a one-component model, with log Bayes factor 6.81 ± 0.02. The coupling time-scale and the torque variances on the crust and superfluid are measured with 90 per cent confidence to be |$10^{7.1^{+0.8}_{-0.5}}$| |$\rm {s}$| and |$10^{-24.0^{+0.4}_{-5.6}}$| |$\rm {rad^2~s^{-3}}$| and |$10^{-21.7^{+3.5}_{-0.9}}$| |$\rm {rad^2~s^{-3}}$|⁠ , respectively. [ABSTRACT FROM AUTHOR]

The S-shaped swing of the linear polarization position angle (PPA) observed in many pulsars can be interpreted by the rotating vector model (RVM). However, efforts to fit the RVM for a large sample of pulsars observed with the MeerKAT telescope as a part of the Thousand-Pulsar-Array (TPA) programme, only succeeded for about half the cases. High time-resolution studies suggest that the failed cases arise due to the presence of orthogonal polarization modes, or highly disordered distribution of PPA points. One such example is PSR J1645−0317. Recently it has been shown that the RVM can be recovered in this pulsar by using only time samples which are greater than 80 per cent linearly polarized. In this work, we test this novel approach on the brightest 249 pulsars from the TPA sample, of which 177 yield sufficient highly polarized samples to be amenable to our method. Remarkably, only nine of these pulsars (5 per cent) now fail to fit the RVM as opposed to 59 per cent from the original analysis. This result favours the paradigm that the underlying mechanism is coherent curvature radiation. [ABSTRACT FROM AUTHOR]

We analyse 13 yr of Fermi Large Area Telescope Pass 8 events from 127 gamma-ray emitting millisecond pulsars (MSPs) in the energy range 0.1–100 GeV and significantly detect 118 MSPs. We fit the stacked emission with a log parabola (LP) spectral model that we show is preferred to two previously published models. We consider the influence of pulsar properties and observed geometric effects on spectral features by defining energy flux colours for both the individual MSPs, and our stacked model as a baseline. There is no correlation of colours with pulsar luminosity, |$\dot{E}$|⁠ , surface magnetic field, or magnetic impact angle. We also find that pulsar geometry has little effect on the observed gamma-ray spectrum, which is in tension with previous modelling of gamma-ray emission with respect to pulsar geometry. Our LP MSP model is applicable to problems where an ensemble of gamma-ray MSPs is considered, such as that of the Galactic Centre excess or in the case of emission from globular clusters. [ABSTRACT FROM AUTHOR]

The common-spectrum process, characterized by the Hellings–Downs angular correlation and observed by pulsar timing array collaborations, such as NANOGrav, PPTA, EPTA, and CPTA, can be explained by the scalar-induced gravitational waves (SIGWs). The energy density of SIGWs exhibits universal behavior in the infrared regions. Utilizing a broken power law parameterization for the primordial curvature power spectrum, we clarify the PTA signal through the infrared characteristics of the SIGWs, using Bayesian analysis to provide posterior distributions. Bayesian factors emphasize the statistical preference for the SIGW model over explanations involving supermassive black hole binaries. [ABSTRACT FROM AUTHOR]

Radio observations from normal pulsars indicate that the coherent radio emission is excited by curvature radiation from charge bunches. In this review, we provide a systematic description of the various observational constraints on the radio emission mechanism. We have discussed the presence of highly polarized time samples where the polarization position angle follow two orthogonal well-defined tracks across the profile that closely match the rotating vector model in an identical manner. The observations also show the presence of circular polarization, with both the right and left handed circular polarization seen across the profile. Other constraints on the emission mechanism are provided by the detailed measurements of the spectral index variation across the profile window, where the central part of the profile, corresponding to the core component, has a steeper spectrum than the surrounding cones. Finally, the detailed measurements of the subpulse drifting behaviour can be explained by considering the presence of non-dipolar field on the stellar surface and the formation of the partially screened gap (PSG) above the polar cap region. The PSG gives rise to a non-stationary plasma flow that has a multi-component nature, consisting of highly energetic primary particles, secondary pair plasma, and iron ions discharged from the surface, with large fragmentation resulting in dense plasma clouds and lower-density inter-cloud regions. The physical properties of the outflowing plasma and the observational constraints lead us to consider coherent curvature radiation as the most viable explanation for the emission mechanism in normal pulsars, where propagation effects due to adiabatic walking and refraction are largely inconsequential. [ABSTRACT FROM AUTHOR]

In this study, we present an analysis of the fine-tuning required in various inflationary models in order to explain the production of Primordial Black Holes (PBHs). We specifically examine the degree of fine-tuning necessary in two prominent single-field inflationary models: those with an inflection point and those with step-like features in the potential. Our findings indicate that models with step-like features generally require less fine-tuning compared to those with an inflection point, making them more viable for consistent PBH production. An interesting outcome of these models is that, in addition to improved fine-tuning, they may also predict low-frequency signals that can be detected by pulsar timing array (PTA) collaborations. Additionally, we extend our analysis to multifield inflationary models to assess whether the integration of additional fields can further alleviate the fine-tuning demands. The study also explores the role of a spectator field and its impact on the fine-tuning process. Our results indicate that although mechanisms involving a spectator field can circumvent the issue of fine-tuning parameters for PBH production, both multifield models and models with step-like features present promising alternatives. While fine-tuning involves multiple considerations, our primary objective is to evaluate various inflationary models to identify the one that most naturally explains the formation of PBHs. Hence, this study introduces a novel approach by categorizing existing PBH mechanisms, paving the way for subsequent research to prioritize models that minimize the need for extensive fine-tuning. [ABSTRACT FROM AUTHOR]

Millisecond pulsars (MSPs) are believed to be very old neutron stars (NSs) whose age may exceed significantly 108$$ {10}^8 $$ years. Although cooling scenarios of isolated NSs predict for that age a surface temperature Ts~104$$ {T}_s\sim {10}^4 $$ K, observations of the nearest MSP J0437‐4715 indicate Ts$$ {T}_s $$ well above that value. Besides the heating of the polar cap surface by backflowing charged particles, Joule heating in the crust can contribute to the overall heat budget of MSPs. Because the dipolar field component, derived from P$$ P $$ and P˙$$ \dot{P} $$ measurements, is much too weak for remarkable heating, smaller‐scale structures should be analyzed to understand whether they can supply the demanded heat. For this purpose, we study the small‐scale field structure of radio pulsars. Magnetic field components, significantly stronger than the dipolar one, may exist especially at the surface of MSPs. We assign upper limits to the strength of single‐field components up to a multipolarity of l=10$$ l=10 $$ and the corresponding deviations from axial symmetry m≤l$$ m\le l $$. Arguments are provided that the decay of the small‐scale components with l=3$$ l=3 $$ or l=4$$ l=4 $$ of the crustal magnetic field may cause the relatively high surface temperature of isolated MSPs. [ABSTRACT FROM AUTHOR]

The Vela pulsar is a young neutron star with spin period of P = 89 ms and a measured low braking index (˜1.4) that is much less than the standard value of 3 predicted by the magnetic dipole radiation (MDR) model; however, its spin period evolution has been a mystery. In this article, we assume that the spin‐down of the Vela pulsar is attributed to both MDR and wind flow (hereafter MDRW), and find that the ratio of wind flow to the magnetic dipole radiation is about 80%, which is higher than that of the Crab pulsar (25%). In other words, the spin‐down torque of the Vela pulsar is dominated by the wind flow. The spin period (P) evolution of the Vela pulsar depends on its real age, where its supernova remnant age is assumed to be an indicator of its true age, estimated from 10 to 20 kyr, and then we obtain their initial spin periods of ˜53.89 and ˜20.90 ms, respectively, which are consistent with the observed initial spin period ranges of young pulsars. Furthermore, we find that the Vela‐like pulsar by MDRW can evolve to the long spin period of a thousand of seconds in less than million years, which can conveniently help the astronomers understand the recently observed ultra‐long period radio pulsars like GPM J1839‐10 (P = 21 min), GLEAM‐X J1627 (P = 18 min), as well as PSR J0901+4046 (P = 76 s). [ABSTRACT FROM AUTHOR]

The broad-spectrum polarization and spectral characteristics of pulsars contain crucial information about the origin of their radio emission. These properties, together with pulsar flux density variations, can also be used to guide future surveys of radio pulsars and probe the Galactic interstellar medium. In this paper, we present studies of 18 pulsars at high Galactic latitudes using the Ultra-Wideband Low (UWL) receiver of the Parkes radio telescope. For these pulsars, we measured their wideband flux densities, spectral indices, and polarization fractions. We obtain seven new rotation measures (RMs) and refine the RMs of another ten pulsars. In this sample of pulsars, we observed significant variations in their flux densities, suggesting that previous shallow surveys were likely to miss a population of pulsars at high galactic latitude. In addition, we identified a previously reported isolated pulsar (PSR J1947−18) as a potential binary system. [ABSTRACT FROM AUTHOR]

Traditional X-ray pulsar ranging and velocity measurement methods only estimate the radial position and velocity information of the pulsar. For non-linear orbits, errors in the transverse position and velocity of the pulsar lead to errors in the radial velocity of the pulsar, leading to distortion of the X-ray pulsar profile. Based on this, we propose using the distortion of the pulsar profile to infer the transverse position and velocity information of the pulsar. First, a model of the distortion of the pulsar profile due to errors in the transverse position and velocity is established, and the observable directions of the transverse position and velocity are given separately. Then, considering that the distortions in the pulsar profile caused by errors in the transverse position and velocity are indistinguishable, we establish a reactive motion state measure related to the observable directions for the transverse position and velocity errors as a new observable measure in X-ray pulsar navigation. The experimental results show that the precision of the reactive motion state measure reaches 0.57, equivalent to a position error of 284.50 m or a velocity error of 0.57 m/s. [ABSTRACT FROM AUTHOR]

It is expected that the extreme mass accretion rate onto strongly magnetized neutron stars results in the appearance of accretion columns above the stellar surface. For a distant observer, rotation of a star results in periodic variations of X-ray flux. Because the mass accretion rate fluctuates around the average value, the pulse profiles are not stable and demonstrate fluctuations as well. In the case of bright X-ray pulsars, however, pulse fluctuations are not solely attributed to variations in the mass accretion rate. They are also influenced by the variable height of the columns, which is dependent on the mass accretion rate. This study delves into the process of pulse profile formation in bright X-ray pulsars, taking into account stochastic fluctuations in the mass accretion rate, the corresponding variations in accretion column geometry, and gravitational bending. Our analysis reveals that potential eclipses of accretion columns by a neutron star during their spin period should manifest specific features in pulse profile variability. Applying a novel pulse profile analysis technique, we successfully detect these features in the bright X-ray transient V 0332+53 at luminosities |$\gtrsim 2\times 10^{38}\, {\rm erg\ \rm s^{-1}}$|⁠. This detection serves as compelling evidence for the eclipse of an accretion column by a neutron star. Detection of the eclipse places constraints on the relation between neutron star mass, radius, and accretion column height. Specifically, we can establish an upper limit on the accretion column height, which is crucial for refining theoretical models of extreme accretion. [ABSTRACT FROM AUTHOR]

Electron densities in different locations of our galaxy are obtained in pulsar astronomy by dividing the dispersion measure (DM) by the distance of the pulsar to Earth. The properties of the interstellar plasma are related to its heating by cosmic rays. Following the present analysis, DM's measurements are obtained with different properties of the corresponding plasmas at different temperatures. (1) For relatively low temperatures (around and below 10,000 (K)), the state of molecular, atomic and ionized hydrogen is analyzed by the interstellar medium model with partially ionized plasma. In this region, various spectroscopic effects are obtained. (2) For temperatures above this limit, the interstellar gas is found to be a completely ionized medium, and this plasma is defined as warm ionized medium (WIM) plasma. We show in the first part of the paper that this plasma is transparent, as obtained from the solution of Saha's equation. The index of refraction of the WIM plasma is real, but the plasma can be observed by dispersion measurements. (3) For very high temperatures (around 1 million (K)), the plasma is defined as hot ionized medium, where X-rays are obtained. We concentrate in the present work on the analysis of WIM plasma. We calculate the mass densities of this plasma and compare it with dark matter mass densities, which are found to be larger. But some factors which may reduce this difference are discussed. [ABSTRACT FROM AUTHOR]

Recently, two ultralong-period radio transients, GLEAM-X J162759.5-523504.3 (J1627) and GPM J1839-10 (J1839), were discovered with spin periods longer than 1000 s. The origin of these two ultralong-period radio transients is intriguing in understanding the spin evolution of neutron stars (NSs). In this work, we examine whether the interaction between strong magnetized NSs and fallback disks can spin NSs down to the observed ultralong period. Our simulations found that the magnetar + fallback disk model can account for the observed period, period derivative, and X-ray luminosity of J1627 in the quasi-spin-equilibrium stage. To evolve to the current state of J1627, the initial mass-accretion rate of the fallback disk and the magnetic field of the NS are in the range of (1.1–30) × 1024 g s−1 and (2–5) × 1014 G, respectively. In the active lifetime of the fallback disk, it is impossible for J1839 to achieve the observed upper limit of the period derivative. Therefore, we propose that J1839 may be in the second ejector phase after the fallback disk becomes inactive. Those NSs with a magnetic field of (2–6) × 1014 G and a fallback disk with an initial mass-accretion rate of ∼1024–1026 g s−1 are possible progenitors of J1839. [ABSTRACT FROM AUTHOR]

The formation of e± plasmas within pulsar magnetospheres through quantum electrodynamics (QED) cascades in vacuum gaps is widely acknowledged. This paper aims to investigate the effect of photon polarization during the QED cascade occurring over the polar cap of a pulsar. We employ a Monte Carlo-based QED algorithm that accurately accounts for both spin and polarization effects during photon emission and pair production in both single-particle and particle-in-cell (PIC) simulations. Our findings reveal distinctive properties in the photon polarization of curvature radiation (CR) and synchrotron radiation (SR). CR photons exhibit high linear polarization parallel to the plane of the curved magnetic field lines, whereas SR photons, on average, demonstrate weak polarization. As the QED cascade progresses, SR photons gradually dominate over CR photons, thus reducing the average degree of photon polarization. Additionally, our study highlights an intriguing observation: the polarization of CR photons enhances e± pair production by approximately 5 per cent, in contrast to the inhibition observed in laser–plasma interactions. Our self-consistent QED-PIC simulations in the corotating frame reproduce the essential results obtained from single-particle simulations. [ABSTRACT FROM AUTHOR]

Polarized radio emission of RRAT J1854+0306 is investigated with single pulses using Five-hundred-meter Aperture Spherical Telescope. Its emission is characterized by nulls, narrow and weak pulses, and occasional wide and intense bursts with a nulling fraction of 53.2 per cent. Its burst emission is typically of one rotation, and occasionally of two or three or even five rotations at the most, but without significant periodicity. The integrated pulse profile has an 'S'-shaped position angle curve that is superposed with orthogonal modes, from which geometry parameters are obtained. Individual pulses exhibit diverse profile morphology with single, double, or multiple peaks. The intensity and width of these pulses are highly correlated, and bright pulses generally have wide profiles with multiple peaks. These nulling behaviours, profile morphology, and polarization demonstrate that a rotating radio transient has the same physical origins as the normal pulsars. [ABSTRACT FROM AUTHOR]

This study investigates the influence of different braking mechanisms on the formation of three long-period radio pulsars (PSRs J0250+5854, J2251−3711, and J0901−4046): plasma-filled magnetosphere in combination with magnetic field decay, fall-back disc, and r-mode instability. These braking mechanisms can also affect the thermal evolution of pulsars. By comparing the model-predicted values with observational data such as spin periods, period derivatives, and upper limits of the bolometric luminosity of these pulsars, we find that these three braking mechanisms can reasonably explain the spin period and the period derivative within a certain range of parameters for these sources. The model-predicted bolometric luminosity associated with magnetic field dissipation exceeds the upper limit for PSR J0901−4046 but falls below the upper limits for PSR J0250+5854 and PSR J2251−3711. The model-predicted bolometric luminosity within the fall-back disc model exceeds the observed results, whereas the luminosity within the r-mode instability falls below the observed upper limit for these three pulsars. However, a conflict arises when we consider the pulsar radio activity and the accretion phases within the fall-back disc model simultaneously. By combining data from radio and X-ray observations, along with precise measurements of surface thermal emissions, we can enhance our understanding of the braking mechanisms involved in the formation of long-period radio pulsars or constrain key model parameters. Finding more long-period pulsars in the future and conducting multiband observations will enhance our understanding of the formation and nature of long-period radio pulsars. [ABSTRACT FROM AUTHOR]

We report here on the timing of 597 pulsars over the last four years with the MeerKAT telescope. We provide times of arrival, pulsar ephemeris files, and per-epoch measurements of the flux density, dispersion measure (DM), and rotation measure (RM) for each pulsar. In addition, we use a Gaussian process to model the timing residuals to measure the spin frequency derivative at each epoch. We also report the detection of 11 glitches in nine individual pulsars. We find significant DM and RM variations in 87 and 76 pulsars, respectively. We find that the DM variations scale approximately linearly with DM, which is broadly in agreement with models of the ionized interstellar medium. The observed RM variations seem largely independent of DM, which may suggest that the RM variations are dominated by variations in the interstellar magnetic field on the line of sight, rather than varying electron density. We also find that normal pulsars have around 5 times greater amplitude of DM variability compared to millisecond pulsars, and surmise that this is due to the known difference in their velocity distributions. [ABSTRACT FROM AUTHOR]

This work presents a methodological approach to generate realistic γ-ray light curves of pulsars, resembling reasonably well the observational ones observed by the Fermi -Large Area Telescope (Fermi -LAT) instrument, fitting at the same time their high-energy spectra. The theoretical light curves are obtained from a spectral and geometrical model of the synchro-curvature emission. Despite our model relying on a few effective physical parameters, the synthetic light curves present the same main features observed in the observational γ-ray light-curve zoo, such as the different shapes, variety in the number of peaks, and a diversity of peak widths. The morphological features of the light curves allow us to statistically compare the observed properties. In particular, we find that the proportion of the number of peaks found in our synthetic light curves is in agreement with the observational one provided by the third Fermi -LAT pulsar catalogue. We also found that the detection probability due to beaming is much higher for orthogonal rotators (approaching 100  per cent) than for small inclination angles (less than 20 per cent). The small variation in the synthetic skymaps generated for different pulsars indicates that the geometry dominates over timing and spectral properties in shaping the γ-ray light curves. This means that geometrical parameters such as the inclination angle can be in principle constrained by γ-ray data alone independently of the specific properties of a pulsar. At the same time, we find that γ-ray spectra seen by different observers can slightly differ, opening the door to constraining the viewing angle of a particular pulsar. [ABSTRACT FROM AUTHOR]

Double neutron star (DNS) systems offer excellent opportunities to test gravity theories. We report the timing results of PSR J1901+0658, the first pulsar discovered in the FAST Galactic Plane Pulsar Snapshot (GPPS) Survey. Based on timing observations by FAST over 5 yr, we obtain the phase-coherent timing solution and derive the precise measurements of its position, spin parameters, orbital parameters, and dispersion measure. It has a period of 75.7 ms, a period derivative of 2.169(6) × 10−19 s s−1, and a characteristic age of 5.5 Gyr. This pulsar is in an orbit with a period of 14.45 d and an eccentricity of 0.366. One post-Keplerian parameter, periastron advance, has been well-measured as being 0.00531(9) deg yr−1, from which the total mass of this system is derived to be 2.79(7) M⊙. The pulsar has the mass upper limit of 1.68 M⊙, so the lower limit for the companion mass is 1.11 M⊙. Because PSR J1901+0658 is a partially recycled pulsar in an eccentric binary orbit with such a large companion mass, it should be in a DNS system according to the evolution history of the binary system. [ABSTRACT FROM AUTHOR]

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 c, med ∼ 3 MJ and a mean density of ρ ∼ 11 g cm−3. M62I has an orbital period of 0.51 d and a M c, med ∼ 0.15 M⊙. Neither of these low-mass systems exhibit eclipses. M62J has only been detected in the two Ultra High Frequency 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 (β < −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. [ABSTRACT FROM AUTHOR]

It is believed that cores of neutron stars provide a natural laboratory where exotic high baryon density phases of quantum chromo dynamics (QCD) may exist. In fact, the theoretically well-established neutron superfluid phase is also believed to be found only inside neutron stars. Focus on neutrons stars has tremendously intensified in recent years with the direct detection of gravitational waves (GWs) by LIGO/Virgo from binary neutron star (BNS) merger events which has allowed the possibility of directly probing the properties of the interior of a neutron star. A truly remarkable phenomenon manifested by rapidly rotating neutron stars is in their avatar as Pulsars. The accuracy of pulsar timing can reach the level of one part in 10 1 5 , comparable to that of atomic clocks. Indeed, it was such a great accuracy which had allowed the first indirect detection of GWs from a BNS system. Such an incredible accuracy of pulse timings points to a very interesting possibility. Any deformation of the pulsar, even if it is extremely tiny, has the potential of leaving its imprints on the pulses through introduction of tiny perturbations in the entire moment of inertia (MI) tensor. While, the diagonal components of perturbed MI tensor affect the pulse timings, the off-diagonal components lead to wobbling of pulsar, directly affecting the pulse profile. This opens up a new window of opportunity for exploring various phase transitions occurring inside a pulsar core, through induced density fluctuations, which may be observable as perturbations in the pulse timing as well as its profile. Such perturbations also naturally induce a rapidly changing quadrupole moment of the star, thereby providing a new source of GW emission. Another remarkable possibility arises when we consider the effect of an external GW on neutron star. With the possibility of detecting any minute changes in its configuration through pulse observations, the neutron star has the potential of performing as a Weber detector of GW. This brief review will focus on these specific aspects of a pulsar. Specifically, the focus will be on the type of physics which can be probed by utilizing the effect of changes in the MI tensor of the pulsar on pulse properties. [ABSTRACT FROM AUTHOR]

We have searched for radio pulsations toward 49 Fermi Large Area Telescope (LAT) 1FGL Catalog γ -ray sources using the Green Bank Telescope at 350 MHz. We detected 18 millisecond pulsars (MSPs) in blind searches of the data; 10 of these were discoveries unique to our survey. 16 are binaries, with eight having short orbital periods P B < 1 day. No radio pulsations from young pulsars were detected, although three targets are coincident with apparently radio-quiet γ -ray pulsars discovered in LAT data. Here, we give an overview of the survey and present radio and γ -ray timing results for the 10 MSPs discovered. These include the only isolated MSP discovered in our survey and six short- P B binary MSPs. Of these, three have very-low-mass companions (M c ≪ 0.1 M ⊙) and hence belong to the class of black widow pulsars. Two have more massive, nondegenerate companions with extensive radio eclipses and orbitally modulated X-ray emission consistent with the redback class. Significant γ -ray pulsations have been detected from nine of the discoveries. This survey and similar efforts suggest that the majority of Galactic γ -ray sources at high Galactic latitudes are either MSPs or relatively nearby nonrecycled pulsars, with the latter having on average a much smaller radio/ γ -ray beaming ratio as compared to MSPs. It also confirms that past surveys suffered from an observational bias against finding short- P B MSP systems. [ABSTRACT FROM AUTHOR]

Context. The radio-quiet γ-ray pulsar PSR J2021+4026 is a peculiar Fermi-LAT pulsar showing repeated and quasi-periodic mode changes. Its γ-ray flux shows repeated variations between two states at intervals of ∼3.5 years. These events occur over timescales < 100 days and are correlated with sudden changes in the spin-down rate. Multiwavelength observations also revealed an X-ray phase shift relative to the γ-ray profile for one of the events. PSR J2021+4026 is currently the only known isolated γ-ray pulsar showing significant variability, and thus it has been the object of thorough investigations. Aims. The goal of our work is to study the mode changes of PSR J2021+4026 with improved detail. By accurately characterizing variations in the γ-ray spectrum and pulse profile, we aim to relate the Fermi-LAT observations to theoretical models. We also aim to interpret the mode changes in terms of variations in the structure of a multipolar dissipative magnetosphere. Methods. We continually monitored the rotational evolution and the γ-ray flux of PSR J2021+4026 using more than 13 years of Fermi-LAT data with a binned likelihood approach. We investigated the features of the phase-resolved spectrum and pulse profile, and from these we inferred the macroscopic conductivity, the electric field parallel to the magnetic field, and the curvature radiation cutoff energy. These physical quantities are related to the spin-down rate and the γ-ray flux and therefore are relevant to the theoretical interpretation of the mode changes. We introduced a simple magnetosphere model that combines a dipole field with a strong quadrupole component. We simulated magnetic field configurations to determine the positions of the polar caps for different sets of parameters. Results. We clearly detect the previous mode changes and confirm a more recent mode change that occurred around June 2020. We provide a full set of best-fit parameters for the phase-resolved γ-ray spectrum and the pulse profile obtained in five distinct time intervals. We computed the relative variations in the best-fit parameters, finding typical flux changes between 13% and 20%. Correlations appear between the γ-ray flux and the spectral parameters, as the peak of the spectrum shifts by ∼10% toward lower energies when the flux decreases. The analysis of the pulse profile reveals that the pulsed fraction of the light curve is larger when the flux is low. Finally, the magnetosphere simulations show that some configurations could explain the observed multiwavelength variability. However, self-consistent models are required to reproduce the observed magnitudes of the mode changes. [ABSTRACT FROM AUTHOR]

Recently, there has been renewed interest in the Be X-ray binary (Be/XRB) SAX J1324.4−6200 because of its spatial coincidence with a variable γ-ray source detected by Fermi/LAT. To explore more thoroughly its properties, new X-ray observations were carried out in 2023 by NuSTAR, XMM-Newton, and Swift satellites, jointly covering the energy range from 0.2 − 79 keV. SAX J1324.4−6200 was caught at an X-ray flux of ∼10−11 erg cm−2 s−1. The X-ray spectrum fits well with an absorbed power law with a high energy cutoff. Other acceptable fits require an additional blackbody component (kTbb ≈ 1.1 keV) or a Gaussian in absorption (Egabs ≈ 6.9 keV). We measured a NuSTAR spin period of 175.8127 ± 0.0036 s and an XMM-Newton spin period of 175.862 ± 0.025 s. We show that all the available spin period measurements of SAX J1324.4−6200, spanning 29 yr, are highly correlated with time, resulting in a remarkably stable spin-down of Ṗ = 6.09 ± 0.06 × 10−9 s s−1. We find that if SAX J1324.4−6200 hosts an accretion-powered pulsar, accretion torque models indicate a surface magnetic field of ∼1012 − 13 G. The X-ray properties emerging from our analysis strenghten the hypothesis that SAX J1324.4−6200 belongs to the small group of persistent Be/XRBs. We also performed radio observations with the Parkes Murriyang telescope, to search for radio pulsations. However, no radio pulsations compatible with the rotational ephemeris of SAX J1324.4−6200 were detected. We rule out the hypothesis that SAX J1324.4−6200 is a γ-ray binary where the emission is produced by interactions between the pulsar and the companion winds. Other models commonly used to account for the production of γ-rays in accreting pulsars cannot reproduce the bright emission from SAX J1324.4−6200. We examined other possible mechanisms behind the γ-ray emission and note that there is a ∼0.5% chance probability that an unknown extragalactic active galactic nucleus (AGN) observed through the Galactic plane may coincidentally fall within the Fermi/LAT error circle of the source and be responsible for the γ-ray emission. [ABSTRACT FROM AUTHOR]

Context. Spiral galaxies M31 and M33 are among the γ-ray sources detected by the Fermi Large Area Telescope (LAT). Aims. We aim to model the broadband non-thermal emission of the central region of M31 (a LAT point source) and of the disk of M33 (a LAT extended source), as part of our continued survey of non-thermal properties of local galaxies that includes the Magellanic Clouds. Methods. We analysed the observed emission from the central region of M31 (R < 5.5 kpc) and the disk-sized emission from M33 (R ∼ 9 kpc). For each galaxy, we self-consistently modelled the broadband spectral energy distribution of the diffuse non-thermal emission based on published radio and γ-ray data. All relevant radiative processes involving relativistic and thermal electrons (synchrotron, Compton scattering, bremsstrahlung, and free–free emission and absorption), along with relativistic protons (π0 decay following interaction with thermal protons), were considered, using exact emissivity formulae. We also used the Fermi-LAT-validated γ-ray emissivities for pulsars. Results. Joint spectral analyses of the emission from the central region of M31 and the extended disk of M33 indicate that the radio emission is composed of both primary and secondary electron synchrotron and thermal bremsstrahlung, whereas the γ-ray emission may be explained as a combination of diffuse pionic, pulsar, and nuclear-BH-related emissions in M31 and plain diffuse pionic emission (with an average proton energy density of 0.5 eV cm−3) in M33. Conclusions. The observed γ-ray emission from M33 appears to be mainly hadronic. This situation is similar to other local galaxies, namely, the Magellanic Clouds. In contrast, we have found suggestions of a more complex situation in the central region of M31, whose emission could be an admixture of pulsar emission and hadronic emission, with the latter possibly originating from both the disk and the vicinity of the nuclear black hole. The alternative modelling of the spectra of M31 and M33 is motivated by the different hydrogen distribution in the two galaxies: The hydrogen deficiency in the central region of M31 partially unveils emissions from the nuclear BH and the pulsar population in the bulge and inner disk. If this were to be the case in M33 as well, these emissions would be outshined by diffuse pionic emission originating within the flat central-peak gas distribution in M33. [ABSTRACT FROM AUTHOR]

Since the discovery of pulsars, the rotating-lighthouse model has been the choice of model to explain the radiation pulsation of pulsars. After discovering that some main sequence stars (e.g. CU Virginis) and ultracool dwarfs (e.g. TVLM 513−46546) also behave like pulsars, the lighthouse model was again adopted to explain their pulse signals. Our research found that if we use the magnetic-field oscillation (MO) model, we can explain the pulse radiation results better regardless of whether its source is a neutron star, a main-sequence star, or an ultracool dwarf. We propose a verifiable prediction that can be used to evaluate the MO model. Our prediction is that there is a 90° phase lag between the magnetic field and radio signal of TVLM 513−46546, and the zero-crossing point of the magnetic field is the moment when the direction of the light's circular polarization is reversed. No new observations are needed to check this prediction, but certain existing data needs to be re-mined. [ABSTRACT FROM AUTHOR]

Within the magnetospheric radius, the geometry of accretion flow in X-ray pulsars is shaped by a strong magnetic field of a neutron star. Starting at the magnetospheric radius, accretion flow follows field lines and reaches the stellar surface in small regions located close to the magnetic poles of a star. At low mass accretion rates, the dynamics of the flow is determined by gravitational attraction and rotation of the magnetosphere due to the centrifugal force. At the luminosity range close to the Eddington limit and above it, the flow is additionally affected by the radiative force. We construct a model simulating accretion flow dynamics over the magnetosphere, assuming that the flow strictly follows field lines and is affected by gravity, radiative, and centrifugal forces only. The magnetic field of a neutron star is taken to be dominated by the dipole component of arbitrary inclination with respect to the accretion disc plane. We show that accretion flow becomes unstable at high mass accretion rates and tends to fluctuate quasi-periodically with a typical period comparable to the free-fall time from the inner disc radius. The inclination of a magnetic dipole with respect to the disc plane and strong anisotropy of X-ray radiation stabilize the mass accretion rate at the poles of a star, but the surface density of material covering the magnetosphere fluctuates even in this case. [ABSTRACT FROM AUTHOR]

The Snake is a remarkable Galactic Centre radio filament with a morphology characterized by two kinks along its ∼20 arcmin extent. The major and minor kinks are located where the filament is most distorted from a linear magnetized structure running perpendicular to the Galactic plane. We present Chandra , VLA, and MeerKAT data and report the detection of an X-ray and radio source at the location of the major kink. High-resolution radio images of the major kink reveal a compact source with a steep spectrum with spectral index α ∼ −2.7 surrounded by extended emission. The radio luminosity and steep spectrum of the compact source are consistent with a pulsar. We also show flattening of the spectrum and enhanced synchrotron emissivity away from the position of the major kink along the Snake, which suggests injection of relativistic particles along the Snake. We argue that the major kink is created by a fast-moving (∼500–1000 km s−1) object punching into the Snake, distorting its magnetic structure, and producing X-ray emission. X-ray emission pinpoints an active acceleration site where the interaction is taking place. A secondary kink is argued to be induced by the impact of the high-velocity object producing the major kink. [ABSTRACT FROM AUTHOR]

We have determined positions, proper motions, and parallaxes of 77 millisecond pulsars (MSPs) from ∼3 yr of MeerKAT radio telescope observations. Our timing and noise analyses enable us to measure 35 significant parallaxes (12 of them for the first time) and 69 significant proper motions. Eight pulsars near the ecliptic have an accurate proper motion in ecliptic longitude only. PSR J0955−6150 has a good upper limit on its very small proper motion (<0.4 mas yr−1). We used pulsars with accurate parallaxes to study the MSP velocities. This yields 39 MSP transverse velocities, and combined with MSPs in the literature (excluding those in Globular Clusters) we analyse 66 MSPs in total. We find that MSPs have, on average, much lower velocities than normal pulsars, with a mean transverse velocity of only 78(8) km s−1 (MSPs) compared with 246(21) km s−1 (normal pulsars). We found no statistical differences between the velocity distributions of isolated and binary MSPs. From Galactocentric cylindrical velocities of the MSPs, we derive 3D velocity dispersions of σρ, σϕ, σ z  = 63(11), 48(8), 19(3) km s−1. We measure a mean asymmetric drift with amplitude 38(11) km s−1, consistent with expectation for MSPs, given their velocity dispersions and ages. The MSP velocity distribution is consistent with binary evolution models that predict very few MSPs with velocities >300 km s−1 and a mild anticorrelation of transverse velocity with orbital period. [ABSTRACT FROM AUTHOR]

To the first post-Newtonian order, the orbital angular momentum of the fast-revolving inner binary of the triple system PSR J0337+1715, made of a millisecond pulsar and a white dwarf, induces an annular gravitomagnetic field which displaces the line of apsides of the slower orbit of the other, distant white dwarf by − 1.2 milliarcseconds per year. The current accuracy in determining the periastron of the outer orbit is 63.9 milliarcseconds after 1.38 years of data collection. By hypothesizing a constant rate of measurement of the pulsar's times of arrivals over the next 10 years, assumed equal to the present one, it can be argued that the periastron will be finally known to a ≃ 0.15 milliarcseconds level, while its cumulative gravitomagnetic retrograde shift will be as large as − 12 milliarcseconds. The competing post-Newtonian gravitolectric periastron advance due to the inner binary's masses, nominally amounting to 74.3 milliarcseconds per year, can be presently modelled to an accuracy level as good as ≃ 0.04 milliarcseconds per year. The mismodeling in the much larger Newtonian periastron rate due to the quadrupolar term of the multipolar expansion of the gravitational potential of a massive ring representing the inner binary, whose nominal size for PSR J0337+1715 is 0.17 degrees per year, might be reduced down to the ≃ 0.5 milliarcseconds per year level over the next 10 years. Thus, a first measurement of such a novel form of gravitomagnetism, although undoubtedly challenging, might be, perhaps, feasible in a not too distant future. [ABSTRACT FROM AUTHOR]

We analyzed the timing parameters (the rotational frequency ν , the first (ν ˙) and second (ν ¨) time-derivatives of frequency) and the derived parameters of a sample of pulsars for which ν ¨ (470 pulsars) were recorded in the Australian Telescope National Facility (ATNF) pulsar catalog. We formed various subsamples, those with braking indices n < 0 and n > 0 , and glitching and non-glitching pulsars. Our statistical analyses of the timing and derived parameters indicated some level of differences and similarities among the parameters analyzed. Glitching pulsars appear to have a higher rotational frequency than non-glitching pulsars, and pulsars with n > 0 appear to rotate faster than those with n < 0 . Our results also suggest that glitching pulsars have lower values of | n | (where | n | is the absolute value of the braking index), and it is lower for the subsample with n > 0 than for the subsample with n < 0 . We believe that the results obtained could be useful in understanding the evolution of pulsar spin. [ABSTRACT FROM AUTHOR]

We discuss spherically symmetric dynamical systems in the framework of a general model of f (R) gravity, i.e. f (R) = R e ζ R , where ζ is a dimensional quantity in squared length units [L 2 ]. We initially assume that the internal structure of such systems is governed by the Krori–Barua ansatz, alongside the presence of fluid anisotropy. By employing astrophysical observations obtained from the pulsar SAX J1748.9-2021, derived from bursting X-ray binaries located within globular clusters, we determine that ζ is approximately equal to ± 5 km 2 . In particular, the model is capable of producing stable configurations for SAX J1748.9-2021, encompassing both its geometric and physical characteristics. We show that, within the framework of f (R) gravity, the Krori–Barua ansatz establishes semi-analytical connections between the radial ( p r ) and tangential ( p t ) pressures, and the density (ρ ). These relations are described as p r ≈ v r 2 (ρ - ρ I) and p t ≈ v t 2 (ρ - ρ II) . In this context, v r and v t denote the sound speeds in the radial and tangential directions, respectively. Meanwhile, ρ I pertains to the surface density, and ρ II is derived from the model parameters. These connections are consistent with the equations of state derived from the best-fit solutions identified in the ongoing investigation. Notably, within the framework of f (R) gravity where ζ is negative, the maximum compactness, denoted as C, is inherently limited to values that do not exceed the Buchdahl limit. This contrasts with general relativity or f (R) gravity with ζ positive, where the compactness has the potential to asymptotically reach the black hole threshold ( C → 1 ). The model predictions suggest a central density that largely exceeds the saturation nuclear density, which is ρ nuc = 3 × 10 14 g/cm 3 . Also the surface density ρ I surpasses ρ nuc . We obtain a mass-radius diagram, corresponding to the boundary density, which is consistent with other observational data. [ABSTRACT FROM AUTHOR]

Pulsar wind nebulae (PWNe) are clouds of magnetized relativistic electron/positron plasma supplied from the central pulsar. However, the number of radio-emitting particles inside a PWN is larger than the expectation from the study of pulsar magnetospheres and their origin is still unclear. A stochastic acceleration of externally injected particles by a turbulence inside the PWN is proposed by our previous studies. In this paper, the previous stochastic acceleration model of the PWN broadband spectra is improved by taking into account the time evolution of the turbulent energy and then the total energy balance inside a PWN is maintained. The turbulent energy supplied from the central pulsar is wasted by the backreaction from the stochastic particle acceleration and the adiabatic cooling according to the PWN expansion. The model is applied to the Crab Nebula and reproduces the current broadband emission spectrum, especially the flat radio spectrum, although the time evolution of the turbulent energy (diffusion coefficient) is a bit complicated compared with our previous studies, where we assumed an exponential behavior of the diffusion coefficient. [ABSTRACT FROM AUTHOR]

We report here results from pulse arrival time delay analysis of the eclipsing high-mass X-ray binary (HMXB) pulsar LMC X-4 using observations made with the Rossi X-ray Timing Explorer, XMM –Newton, NuSTAR (Nuclear Spectroscopic Telescope ARray), and AstroSat. Combining the orbital parameters determined from these observations with the historical measurements dating back to 1998, we have extended the T π/2 epoch history of LMC X-4 by about 4600 binary orbits spanning about 18 yr. We also report mid-eclipse time measurements (T ecl) using data obtained from wide-field X-ray monitors of MAXI -GSC (Monitor of All-sky X-ray Image – Gas Slit Camera) and Swift - BAT (Burst Alert Telescope). Combining the new T π/2 and T ecl estimates with all the previously reported values, we have significantly improved the orbital evolution measurement, which indicates that the orbital period is evolving at a time-scale (⁠|$P_{\rm orb}/\dot{P}_{\rm orb}$|⁠) of about 0.8 Myr. For the first time in an accreting X-ray pulsar system, we confirm the existence of a second derivative of the orbital period, having an evolution time-scale (⁠|$\dot{P}_{\mathrm{ orb}}/\ddot{P}_{\mathrm{ orb}}$|⁠) of about 55 yr. Detection of a second derivative of the orbital period in LMC X-4 makes its orbital evolution time-scale more uncertain, which may also be true for other HMXBs. Independent solutions for the orbital evolution measurement using the mid-eclipse data and the pulse timing data are consistent with each other, and help us put an upper limit of 0.009 on the eccentricity of the binary system. [ABSTRACT FROM AUTHOR]

In their catalogue of pulsars' radio spectra, Swainston et al. distinguish between five different forms of these spectra: those that can be fitted with (i) a simple power law, (ii) a broken power law, (iii) a low-frequency turn-over, (iv) a high-frequency turn-over or (v) a double turn-over spectrum. Here, we choose two examples from each of these categories and fit them with the spectral distribution function of the caustics that are generated by the superluminally moving current sheet in the magnetosphere of a non-aligned neutron star. In contrast to the prevailing view that the curved features of pulsars' radio spectra arise from the absorption of the observed radiation in high-density environments, our results imply that these features are intrinsic to the emission mechanism. We find that all observed features of pulsar spectra (including those that are normally fitted with simple or broken power laws) can be described by a single spectral distribution function and regarded as manifestations of a single emission mechanism. From the results of an earlier analysis of the emission from a pulsar's current sheet and the values of the fit parameters for each spectrum, we also determine the physical characteristics of the central neutron star of each considered example and its magnetosphere. [ABSTRACT FROM AUTHOR]

We report a detection of GeV γ-ray emission potentially originating from the pulsar wind nebula in CTA 1 by analysing about 15 yr of Fermi Large Area Telescope data. By selecting an energy range from 50 GeV to 1 TeV to remove contamination from the γ-ray pulsar PSR J0007+7303, we have discovered an extended γ-ray source with a test statistic value of ∼44.94 in the region of CTA 1. The obtained flux is measured to be 6.71 ± 2.60 × 10−12 erg cm−2 s−1 with a spectral index of 1.61 ± 0.36, which allows for a smooth connection with the flux in the TeV band. CTA 1 is also considered to be associated with 1LHAASO J0007+7303u, which is an ultra-high-energy source listed in the recently published catalogue of the Large High Altitude Air Shower Observatory. We assume that the radiation originates from the pulsar wind nebula and that its multiwavelength spectral energy distribution can be explained well with a time-dependent one-zone model. [ABSTRACT FROM AUTHOR]

The origin of fast radio bursts (FRBs), the brightest cosmic explosion in radio bands, remains unknown. We introduce here a novel method for a comprehensive analysis of active FRBs' behaviors in the time-energy domain. Using "Pincus Index" and "Maximum Lyapunov Exponent", we were able to quantify the randomness and chaoticity, respectively, of the bursting events and put FRBs in the context of common transient physical phenomena, such as pulsar, earthquakes, and solar flares. In the bivariate time-energy domain, repeated FRB bursts' behaviors deviate significantly (more random, less chaotic) from pulsars, earthquakes, and solar flares. The waiting times between FRB bursts and the corresponding energy changes exhibit no correlation and remain unpredictable, suggesting that the emission of FRBs does not exhibit the time and energy clustering observed in seismic events. The pronounced stochasticity may arise from a singular source with high entropy or the combination of diverse emission mechanisms/sites. Consequently, our methodology serves as a pragmatic tool for illustrating the congruities and distinctions among diverse physical processes. [ABSTRACT FROM AUTHOR]

The magnetic dipole moment μ of an accretion-powered pulsar in magnetocentrifugal equilibrium cannot be inferred uniquely from time-averaged pulse period and aperiodic X-ray flux data, because the radiative efficiency η 0 of the accretion is unknown, as are the mass, radius, and distance of the star. The degeneracy associated with the radiative efficiency is circumvented if fluctuations of the pulse period and aperiodic X-ray flux are tracked with a Kalman filter, whereupon μ can be measured uniquely up to the uncertainties in the mass, radius, and distance. Here, the Kalman filter analysis is demonstrated successfully in practice for the first time on Rossi X-ray Timing Explorer observations of the X-ray transient SXP 18.3 in the Small Magellanic Cloud (SMC), which is monitored regularly. The analysis yields μ = 8.0 − 1.2 + 1.3 × 10 30 G cm 3 and η 0 = 0.04 − 0.01 + 0.02 , compared to μ = 5.0 − 1.0 + 1.0 × 10 30 G cm 3 as inferred traditionally from time-averaged data assuming η 0 = 1. The analysis also yields time-resolved estimates of two hidden state variables, the mass accretion rate and the Maxwell stress at the disk–magnetosphere boundary. The success of the demonstration confirms that the Kalman filter analysis can be applied in the future to study the magnetic moments and disk–magnetosphere physics of accretion-powered pulsar populations in the SMC and elsewhere. [ABSTRACT FROM AUTHOR]

We report the phase-resolved spectral results of the first Galactic pulsating ultraluminous X-ray source (PULX) Swift J0243.6+6124, modeling its 2017–2018 outburst peak using data collected by the Hard X-ray Modulation Telescope (Insight-HXMT). The broad energy coverage of Insight-HXMT allows us to obtain a more accurate spectral continuum to reduce the coupling of broad iron line profiles with other components. We use three different continuum spectrum models but obtain similar iron line results. For the first time, we detect the pulse characteristics of the broad iron line in a PULX. The variation in the width and intensity of this iron line with σ ∼ 1.2–1.5 keV has a phase offset of about 0.25 from the pulse phase. We suggest that the uneven irradiation of the thick inner disk by the accretion column produces the modulated variation of the broad iron line. In addition, the nonpulsed narrow line is suggested to come from the outer disk region. [ABSTRACT FROM AUTHOR]

The accretion flow within the magnetospheric radius of bright X-ray pulsars can form an optically thick envelope, concealing the central neutron star from the distant observer. Most photons are emitted at the surface of a neutron star and leave the system after multiple reflections by the accretion material covering the magnetosphere. Reflections cause momentum to be transferred between photons and the accretion flow, which contributes to the radiative force and should thus influence the dynamics of accretion. We employ Monte Carlo simulations and estimate the acceleration along magnetic field lines due to the radiative force as well as the radiation pressure across magnetic field lines. We demonstrate that the radiative acceleration can exceed gravitational acceleration along the field lines, and similarly, radiation pressure can exceed magnetic field pressure. Multiple reflections of X-ray photons back into the envelope tend to amplify both radiative force along the field lines and radiative pressure. We analyse the average photon escape time from the magnetosphere of a star and show that its absolute value is weakly dependent on the magnetic field strength of a star and roughly linearly dependent on the mass accretion rate being |$\sim 0.1\, {\rm s}$| at |$\dot{M}\sim 10^{20}\, {\rm g\, s^{-1}}$|⁠. At high mass accretion rates, the escape time can be longer than free-fall time from the inner disc radius. [ABSTRACT FROM AUTHOR]

We present a comprehensive timing analysis on the Be/X-ray binary pulsar RX J0440.9+4431 using observations from NICER and Insight - HXMT during the 2022–2023 outburst. The power density spectrum (PDS) of RX J0440.9+4431 exhibits typical aperiodic variability in X-ray flux across a wide frequency range. During a supercritical accretion state, we detect quasi-periodic oscillations (QPOs) at 0.2–0.5 Hz in the light curves of five pulses for RX J0440.9+4431. The observed QPOs manifest during flares, while the flares appear at the peaks of the pulse profiles on a time-scale of seconds and are primarily caused by an increase in hard photons. These flares can be explained by increased material ingestion in the accretion column at a fixed phase, primarily generating hard photons. Alternatively, an increase in accretion rate, independent of phase, may result in highly beamed hard photons within the accretion column, causing the flares. We argue the origin of QPOs to instabilities within the accretion flow. Additionally, we find that the break frequencies in the noise power spectra align well with |$\propto L_{\mathrm{x}}^{3 / 7}$| across three orders of magnitude in the luminosity, which points to a relatively strong magnetic field in RX J0440.9+4431, estimated to be ~1013 G. [ABSTRACT FROM AUTHOR]

Pulsar time has different physical mechanisms from atomic time, and these two time-scales are complementary in terms of stability performance. The joint timekeeping of pulsar and atomic time is currently an important application direction. Our research utilizes a wavelet analysis method to perform multiscale analysis of time signals in both the time and frequency domains, and assigns corresponding weights to different frequency bands for synthesis, which solves the problem of balancing long- and short-term stability in unified weighted synthesis. In this paper, a joint time-scale APT (atomic pulsar time) is obtained by a fusion of pulsar ensemble time (EPT) and atomic time TA(NTSC) (AT) based on wavelet analysis. The results show that APT effectively integrates the stability advantages of both EPT and AT, and the short-term stability of APT at 0.4 yr is consistent with that of AT, reaching 5.6 × 10−15 as evaluated by σ z variance. The stability of APT remains consistent with EPT after 3.2 yr, and the long-term stability at 13.1 yr is close to 1.3 × 10−16, perfectly retaining the excellent long-term stability characteristics of pulsar time. The root mean square of APT's clock difference relative to TT(BIPM) is much lower than that of TA(NTSC). Moreover, we propose ways to improve the stability of APT through experimental comparisons. The results indicate that APT has the potential to serve as a reference for steering UTC(NTSC) in the future, improving the long-term timekeeping capability of local atomic time in China, thereby comprehensively improving the stability and reliability of China's time reference. [ABSTRACT FROM AUTHOR]

The γ -ray sky above a few tens of megaelectronvolts (MeV) reveals some of the most powerful and energetic phenomena of our Universe. The Astrorivelatore Gamma ad Immagini LEggero (AGILE) Gamma-ray Mission was launched in 2007 with the aim of observing celestial sources by means of three instruments covering a wide range of energies, from hard X-rays up to 30 GeV. Thanks to its wide field of view, AGILE set to observe and detect emission from pulsars, pulsar wind nebulae, gamma-ray bursts, active galactic nuclei, fast radio bursts, terrestrial gamma-ray flashes, and the electromagnetic counterparts of neutrinos and gravitational waves. In particular, the fast on-ground processing and analysis chain allowed the AGILE team to promptly respond to transient events, and activate or participate in multiwavelength observing campaigns. Eventually, after 17 years of operations, the AGILE Italian scientific satellite re-entered the atmosphere on 14 February 2024, ending its intense activity as a hunter of some of the most energetic cosmic sources in the Universe that emit X and γ -rays. We will review the most relevant AGILE results to date and their impact on the advancements of theoretical models. [ABSTRACT FROM AUTHOR]