2025-04-04 | | Total: 11
We investigate the X-ray variability of the supergiant X-ray binary 4U 1907+09 using the new NuSTAR observation of 2024. The source had a relatively stable flux level during previous NuSTAR observations, but the flux varied significantly during the current one. The light curve exhibits dips (off-state) and flares (on-state). The phase-coherent timing analysis during the on-state yields a pulse period of 443.99(4) s, showing the pulsar's continued spin-down. The pulse profiles show an asymmetric double-peaked structure with a phase separation of 0.47 between the two peaks. A cyclotron resonance scattering feature (CRSF) is also detected at ∼17.6 keV, along with its harmonic at ∼38 keV, persisting across all flux states. Flux-resolved spectroscopy reveals that the CRSF remains constant despite a 25-fold change in flux. The spectral parameters like photon index and e-fold energy are out of phase with the pulse shape, whereas cutoff energy is in phase with the pulse shape. The source's luminosity during the on-state is 2.85×1035 erg s−1, consistent with a "pencil" beam radiation pattern expected at this flux level from a collisionless gas-mediated shock. These results offer further insights into the accretion dynamics and magnetic field geometry of this system.
Based on well-grounded Galactic neutron star populations formed from radio pulsar population syntheses of canonical pulsars (CPs) and millisecond pulsars (MSPs), we use the latest Fermi-LAT catalog (4FGL-DR4) to investigate the implications of proposed γ−ray luminosity models. Using Monte Carlo techniques, we calculate the number of CPs and MSPs that would comprise the sample of pulsar-like unidentified sources (PLUIDs) in 4FGL-DR4. While radio beaming fractions were used to scale the sizes of the populations, when forming the mock 4FGL-DR4 samples, we make the simplifying assumption that all γ−ray pulsars are beaming towards the Earth. We then explore the observable outcomes of seven different γ−ray luminosity models. Four of the models provide a good match to the observed number of PLUIDs, while three others significantly over-predict the number of PLUIDs. For these latter models, either the average beaming fraction of γ−ray pulsars is more like 25--50\%, or a revision in the luminosity scaling is required. Most of the radio detectable MSPs that our models predict as part of the PLUIDs within 4FGL-DR4 are, unsurprisingly, fainter than the currently observed sample and at larger dispersion measures. For CPs, in spite of an excellent match to the observed radio population, none of the γ−ray models we investigated could replicate the observed sample of 150 γ−ray CPs. Further work is required to understand this discrepancy. For both MSPs and CPs, we provide encouraging forecasts for targeted radio searches of PLUIDs from 4FGL-DR4 to elucidate the issues raised in this study.
4U 1820-30 is a ultracompact X-ray binary located in the globular cluster NGC 6624, consisting of a neutron star accreting material from a helium white dwarf companion characterized by the shortest known orbital period for this type of star (11.4 minutes). Despite extensive studies, the detection of the relativistic Fe K emission line, has been inconsistently reported and no measurement of the system inclination has been achieved. In this work, we investigate the broadband spectral and polarimetric properties of 4U 1820-30, exploring the presence of a reflection component and its role in shaping the polarization signal. We analyzed simultaneous X-ray observations from NICER, NuSTAR, and IXPE. The spectral continuum was modeled with a disk blackbody, a power law, and a Comptonization component with seed photons originating from a boundary layer. We detected a strong reflection component, described, for the first time, with two self-consistent models ({\tt Relxillns} and {\tt Rfxconv}), allowing us to provide a measurement of the system inclination angle (about 31 degrees), supporting the low-inclination hypothesis. Subsolar iron abundances were detected in the accretion disk and interstellar medium, probably related to the source location in a metal-poor globular cluster. The polarization increases from an upper limit of 1.2% in the 2--4 keV band up to about 8% in the 7--8 keV range. The disk is expected to be orthogonally polarized to these components, which may help to explain the decreasing of the observed polarization at low energies. However, the high polarization degree we found challenges the current models, also taking into consideration the relatively low inclination angle derived from the spectral analysis.
GRB 170519A was discovered by \emph{Swift}/BAT, and then observed by \emph{Swift}/XRT, \emph{Swift}/UVOT, and ground-based telescopes. We report Lick/KAIT observations of GRB 170519A, and make temporal analysis and spectral joint fits of its multiwavelength light curves. The observations present a relatively complete afterglow structure, including two X-ray flares (Flares I and II), optical onset (Slice 1), normal decay (Slices 2 and 3), and a possible jet break. The spectrum of the bright X-ray flare (Flare II) indicates that a thermal component exists at t=190--240~s. \textbf{The blackbody emits in the photospheric radius Rph∼1011 cm,} and its temperature (kT) decreases with time from \textbf{1.08 to 0.37 keV, its Lorentz factor of blackbody (ΓBB) decreases with time from 67.71 to 46.70. The luminosity of the blackbody (LBB), kT and ΓBB follow the relations LBB∝kT2.49±0.03 and ΓBB∝L0.27BB (estimated from \cite{fan2012}).} In the optical light curves, there is an onset bump in the early-time afterglow, rising with an index αO,1≈−0.43 and peaking ∼1174.9 s since the BAT trigger. The bump then decays with αO,2≈0.88 in the normal decay phase, and the X-ray flux decays with a similar index of αX,1≈0.95. There is no obvious spectral evolution in the normal decay phases, with photon index ˆΓ=1.86 and 1.92 in Slices 2 and 3, respectively. We find that the multiwavelength light curves of the GRB 170519A afterglow can be well fitted by an external shock with time-dependent ϵB. In the early afterglow, the value of ϵB decays rapidly from 4.29×10−2 to 8.23×10−3.
We conducted targeted fast radio burst (FRB) and pulsar searches on eight pulsing ultraluminous X-ray sources (PULXs) using the Five-hundred-meter Aperture Spherical Radio Telescope (FAST) and the Parkes 64-meter Radio Telescope (Murriyang) to investigate whether PULXs could be progenitors of FRBs. FAST carried out 12 observations of four PULXs, totaling 8 hours, while Parkes conducted 12 observations of the remaining four PULXs, totaling 11 hours. No significant signals were detected through single-pulse and periodic searches, covering a dispersion measure (DM) range of 0-5000 pc cm−3, placing stringent upper limits on the radio flux density from these sources. The results imply that accretion processes and dense stellar winds in PULXs likely suppress or attenuate potential coherent emission in radio band. Additionally, the beaming factor and luminosity of FRBs associated with PULXs, as well as the highly relativistic and magnetized nature of their outflows, may limit detectability. Non-detection yielded from the observations covering the full orbital phases of PULXs can also constrain the theoretical models that link FRB emission to highly magnetized neutron stars in binary systems.
Pulsar halos are regions around middle-aged pulsars extending out to tens of parsecs. The large extent of the halos and well-defined central cosmic-ray accelerators make this new class of Galactic sources an ideal laboratory for studying cosmic-ray transport. LHAASO J0621+3755 is a candidate pulsar halo associated with the middle-aged gamma-ray pulsar PSR J0622+3749. We observed LHAASO J0621+3755 with VERITAS and XMM-Newton in the TeV and X-ray bands, respectively. For this work, we developed a novel background estimation technique for imaging atmospheric Cherenkov telescope observations of such extended sources. No halo emission was detected with VERITAS (0.3--10 TeV) or XMM-Newton (2--7 keV) within 1 degree and 10 arcmin around PSR J0622+3749, respectively. Combined with the LHAASO-KM2A and Fermi-LAT data, VERITAS flux upper limits establish a spectral break at ~1--10 TeV, a unique feature compared with Geminga, the most studied pulsar halo. We model the gamma-ray spectrum and LHAASO-KM2A surface brightness as inverse Compton emission and find suppressed diffusion around the pulsar, similar to Geminga. A smaller diffusion suppression zone and harder electron injection spectrum than Geminga are necessary to reproduce the spectral cutoff. A magnetic field <= 1 uG is required by our XMM-Newton observation and synchrotron spectral modeling, consistent with Geminga. Our findings support slower diffusion and lower magnetic field around pulsar halos than the Galactic averages, hinting at magnetohydrodynamic turbulence around pulsars. Additionally, we report the detection of an X-ray point source spatially coincident with PSR J0622+3749, whose periodicity is consistent with the gamma-ray spin period of 333.2 ms. The soft spectrum of this source suggests a thermal origin.
Energetic bursts from strongly magnetized neutron stars, known as magnetars, are typically detected in clusters. Once an active episode begins, anywhere from a few to thousands of hard X-ray bursts can occur over durations ranging from days to months. The temporal clustering of these recurrent bursts during an active episode suggests an underlying mechanism that triggers multiple bursts in rapid succession. These burst clusters are likely crucial for understanding the processes driving magnetar activity. In this study, we investigate the repetitive short X-ray burst behavior of magnetars through crustal interactions, employing the cellular automaton model for the magnetar crust proposed by Lander (2023). Our simulations, based on physically motivated criteria, successfully reproduce burst clustering. Additionally, the durations and energetics of active episodes in our simulations agree well with observational data. We discuss the potential physical mechanisms underlying burst clusters observed in numerous magnetars, as well as the reactivations of an individual magnetar.
In 2021, RS Ophiuchi was the first nova to be detected in the very-high-energy (TeV) gamma-ray domain, directly testifying of efficient acceleration of charged particles up to at least the TeV range at the nova shock. Surprisingly, the TeV gamma-ray signal peaks ∼2 days after the GeV signal and the origin of this delay has still not been clearly understood. We investigate the possibility that this delay is due to the effect of gamma-ray absorption resulted from interactions between gamma rays and optical photons copiously emitted during the outburst. We model particle acceleration at a nova shock to obtain the gamma-ray emission produced in interactions between the accelerated particles and the shocked gas. The effect of gamma-ray absorption is then included in details using the radiative transfer equation. We find that this can naturally account for the delay between the peaks of GeV and TeV gamma-ray lightcurves. This result emphasizes the importance of gamma-ray absorption for interpreting gamma-ray observations of novae in the TeV range which, in turn, demonstrates the necessity of a multi-wavelength view for unraveling the underlying physics of particle acceleration in these systems.
We present a modified Richardson-Lucy (RL) algorithm tailored for image reconstruction in MeV gamma-ray observations, focusing on its application to the upcoming Compton Spectrometer and Imager (COSI) mission. Our method addresses key challenges in MeV gamma-ray astronomy by incorporating Bayesian priors for sparseness and smoothness while optimizing background components simultaneously. We introduce a novel sparsity term suitable for Poisson-sampled data in addition to a smoothness prior, allowing for flexible reconstruction of both point sources and extended emission. The performance of the algorithm is evaluated using simulated three-month COSI observations of gamma-ray lines of 44Ti (1.157 MeV), 26Al (1.809 MeV), and positron annihilation (0.511 MeV), respectively, representing various spatial features. Our results demonstrate significant improvements over conventional RL methods, particularly in suppressing artificial structures in point source reconstructions and retaining diffuse spatial structures. This work represents an important step towards establishing a robust data analysis for studying nucleosynthesis, positron annihilation, and other high-energy phenomena in our Galaxy.
Axion-like particles (ALPs) can be copiously produced in binary neutron star (BNS) mergers through nucleon-nucleon bremsstrahlung if the ALP-nucleon couplings gaN are sizable. Furthermore, the ALP-photon coupling gaγ may trigger conversions of ultralight ALPs into photons in the magnetic fields of the merger remnant and of the Milky Way. This effect would lead to a potentially observable short gamma-ray signal, in coincidence with the gravitational-wave signal produced during the merging process. This event could be detected through multi-messenger observation of BNS mergers employing the synergy between gravitational-wave detectors and gamma-ray telescopes. Here, we study the sensitivity of current and proposed MeV gamma-ray experiments to detect such a signal. We find that the proposed instruments can reach a sensitivity down to gaγ≳few×10−13GeV−1 for ma≲10−9 eV, comparable with the SN 1987A limit.
In light-shining-through-walls experiments, axions and axion-like particles (ALPs) are searched for by exposing an optically thick barrier to a laser beam. In a magnetic field, photons could convert into ALPs in front of the barrier and reconvert behind it, giving rise to a signal that can occur only in the presence of such hidden particles. In this work, we utilize the light-shining-through-walls concept and apply it to astrophysical scales. Namely, we consider eclipsing binary systems, consisting of a neutron star, which is a bright source of X-rays, and a companion star with a much larger radius. Space observatories such as XMM-Newton and NuSTAR have performed extensive measurements of such systems, obtaining data on both out-of-eclipse photon rates and those during eclipses. The latter are typically O(102−103) times smaller, due to the fact that X-rays propagating along the line of sight from the neutron star to the X-ray observatory do not pass through the barrier that is the companion star. Using this attenuation, we derive a constraint on ALP-photon coupling of gaγ≃10−10GeV−1 for the LMC X-4 eclipsing binary system, surpassing current bounds from light-shining-through-walls experiments by several orders of magnitude. We also present future prospects that could realistically improve this limit by an order of magnitude in gaγ, making it competitive with some of the strongest limits derived to date.