High Energy Astrophysical Phenomena (astro-ph.HE)

Abstract: Swift J1749.4-2807 is the only known eclipsing accreting millisecond X-ray pulsar. In this paper, we report on 7 thermonuclear (Type-I) X-ray bursts observed by NICER during its 2021 outburst. The first 6 bursts show slow rises and long decays, indicative of mixed H/He fuel, whereas the last burst shows fast rise and decay, suggesting He-rich fuel. Time-resolved spectroscopy of the bursts revealed typical phenomenology (i.e., an increase in black body temperature during the burst rise, and steady decrease in the decay), however they required a variable $N_\mathrm{H}$. We found that the values of $N_\mathrm{H}$ during the bursts were roughly double those found in the fits of the persistent emission prior to each burst. We interpret this change in absorption as evidence of burst-disc interaction, which we observe due to the high inclination of the system. We searched for burst oscillations during each burst and detected a signal in the first burst at the known spin frequency of the neutron star (517.92 Hz). This is the first time burst oscillations have been detected from Swift J1749.4-2807. We further find that each X-ray burst occurs on top of an elevated persistent count rate. We performed time-resolved spectroscopy on the combined data of the bursts with sufficient statistics (i.e., the clearest examples of this phenomenon) and found that the black body parameters evolve to hotter temperatures closer to the onset of the bursts. We interpret this as a consequence of an unusual marginally stable burning process similar to that seen through mHz QPOs.

Abstract: Neutron stars are known to be efficient accelerators producing particles with ultra-relavistic energies. As a by product they also emit copiously photons from radio wavelengths up to gamma-rays. As a follow up of our previous work on particle acceleration simulation near neutron stars, in this paper, we discuss the impact of radiation reaction on test particles injected into their magnetosphere. We therefore neglect the interaction between particles through the electromagnetic field as well as gravitation. We found that while due solely to the Lorentz force electrons reach Lorentz factors up to $\gamma=10^{14}$ and protons up to $\gamma=10^{10.7}$, when radiation reaction is enabled electrons reach energies up to $\gamma=10^{10.5}$ and protons up to $\gamma=10^{8.3}$. The latter values are more realistic since the radiation reaction feedback is predominant within the magnetosphere. Moreover, as expected, symmetrical behaviours between the north and the south hemispheres is highlighted, either with respect to the location around the neutron star or with respect to particles of opposite charge to mass ratio~$q/m$. Consequently, it is useless to simulate the full set of geometrical parameters to get an overview of all possibilities. The study of the influence of the magnetic dipolar moment inclination shows similar behaviours whether or not radiation reaction is switch on: protons (respectively electrons) impact less the surface of the neutron star as the inclination angle increases (decreases for electrons) while if the rotation and magnetic axes are aligned, all the protons impact onto the neutron star, and all the electrons impact the surface if the rotation and magnetic axes are anti-aligned. Similarly, we still find that particles can be ejected away from the neutron star, with preferred directions, were it for the number of particles or for their Lorentz factors

Abstract: In this work, we analyzed the long term gamma-ray data by a Fermi Large Area Telescope (Fermi-LAT) of blazar S2 0109+22, ranging from 2008 to 2023. The quasi-periodic oscillations (QPOs) of blazars aided in investigating the physical properties of internal supermassive black holes, the nature of variability, and the underlying radiation mechanism. We employed four different methods--Weighted Wavelet Z-transform, Lomb-Scargle periodogram, REDFIT and phase folded light curve analysis, for searching QPO signals. Our analysis identified a possible QPO behavior with a periodicity of $\sim$600 days in November 2013 to January 2023 at a significance level of 3.5 $\sigma$. This QPO signal sustained $\sim$9 years, corresponding to 5.6 cycles, which was in good agreement with the previously observed of periodicity $\sim$657 days in radio. We explained this phenomenon based on the accretion model and the lighthouse effect, in a binary black hole system.

Abstract: Using the Hubble Ultraviolet Globular Cluster Survey (HUGS) and additional HST archival data, we have carried out a search for optical counterparts to the low-luminosity Chandra X-ray sources in the globular cluster M4 (NGC 6121). We have also searched for optical or X-ray counterparts to radio sources detected by the VLA. We find 24 new confident optical counterparts to Chandra sources for a total of 40, including the 16 previously identified. Of the 24 new identifications, 18 are stellar coronal X-ray sources (active binaries, ABs), the majority located along the binary sequence in a V-I colour-magnitude diagram and generally showing an H-alpha excess. In addition to confirming the previously detected cataclysmic variable (CV, CX4), we identify one confident new CV (CX76), and two candidates (CX81 and CX101). One MSP is known in M4 (CX12), and another strong candidate has been suggested (CX1); we identify some possible MSP candidates among optical and radio sources, such as VLA20, which appears to have a white dwarf counterpart. One X-ray source with a sub-subgiant optical counterpart and a flat radio spectrum (CX8, VLA31) is particularly mysterious. The radial distribution of X-ray sources suggests a relaxed population of average mass ~ 1.2 - 1.5 Msun. Comparing the numbers of ABs, MSPs, and CVs in M4 with other clusters indicates that AB numbers are proportional to cluster mass (primordial population), MSPs to stellar encounter rate (dynamically formed population), while CVs seem to be produced both primordially and dynamically.

Abstract: We study the motion of neutron superfluid vortices in a spinning-down neutron star, assuming axisymmetry of the flow and ignoring motion of vortices about the rotation axis. We find that the vortex array, if initially rectilinear, is soon substantially deformed as the star spins down; vortices are swept outward by the Magnus force, accumulating in regions of the inner crust where they pin, accompanied by significant bending of the vortex array. As the star spins down to below a spin rate of ~20 Hz (twice the spin rate of the Vela pulsar), the Magnus and pinning forces gradually compress the vortex array into dense sheets that follow spherical shells. In some cases, the vortex array bends on itself and reconnects, forming one or more tori of vortex rings that contain superfluid ``rivers" with significant angular momentum. Vortex sheets are likely to form near the base of the inner crust, in the regime of nuclear pasta.

Abstract: We present nebular spectroscopy of SN 2020hvf, a Type Ia supernova (SN Ia) with an early bump in its light curve. SN 2020hvf shares many spectroscopic and photometric similarities to the carbon-rich high-luminosity "03fg-like" SNe Ia. At $>$240 days after peak brightness, we detect unambiguous emission from [Ca II] $\lambda\lambda$7291, 7324 which is never observed in normal-SNe Ia and only seen in peculiar subclasses. SN 2020hvf displays "saw-tooth" emission profiles near 7300 A that cannot be explained with single symmetric velocity components of [Fe II], [Ni II], and [Ca II], indicating an asymmetric explosion. The broad [Ca II] emission is best modeled by two velocity components offset by 1,220 km s$^{-1}$, which could be caused by ejecta associated with each star in the progenitor system, separated by their orbital velocity. For the first time in a SN Ia, we identify narrow (${\rm FWHM} = 180\pm40$ km s$^{-1}$) [Ca II] emission, which we associate with a wind from a surviving, puffed-up companion star. Few published spectra have sufficient resolution and signal-to-noise ratio necessary to detect similar narrow [Ca II] emission, however, we have detected similar line profiles in other 03fg-like SNe Ia. The extremely narrow velocity width of [Ca II] has only otherwise been observed in SNe Iax at late times. Since this event likely had a double-degenerate "super-Chandrasekhar" mass progenitor system, we suggest that a single white dwarf (WD) was fully disrupted and a wind from a surviving companion WD is producing the observed narrow emission. It is unclear if this unique progenitor and explosion scenario can explain the diversity of 03fg-like SNe Ia, potentially indicating that multiple progenitor channels contribute to this subclass.

Abstract: The origin of astrophysical neutrinos is one of the most debated topics today. Perhaps the most robust evidence of neutrino counterpart comes from supermassive black holes in active galactic nuclei associated with strongly collimated outflows, or jets, that can accelerate particles to relativistic energies and produce neutrinos through hadronic interactions. Similar outflows can also be found from X-ray binaries, or `microquasars', that consist of a neutron star or a stellar-mass black hole accreting matter from a non-degenerate companion star. In some cases, these systems can accelerate particles up to GeV energies implying an efficient acceleration mechanism in their jets. Neutrino production in microquasar jets can be expected with suitable conditions and a hadronic particle population. Microquasar Cyg X-3 is a unique, short orbital period X-ray binary hosting a Wolf-Rayet companion star with a strong stellar wind. The interaction of the dense stellar wind with a relativistic jet leads to particle collisions followed by high-energy gamma-ray and potentially neutrino emission. Here, using the 10-year neutrino candidate sample of the IceCube neutrino observatory, we find that the events with the highest spatial association with Cyg X-3 occur during short-lived high-energy gamma-ray flaring periods indicating the possible astrophysical nature of these events.

Abstract: To determine the nature of the PeVatron's emission (hadronic or leptonic), it is essential to characterize the physical parameters of the environment from where it originates. We unambiguously confirm the association of molecular gas with the PeVatron candidate LHAASO J2108+5157 using unprecedented high angular-resolution (17$^{\prime \prime}$) $^{12,13}$CO($J$=1$\rightarrow$0) observations carried out with the Nobeyama 45m radio telescope. We characterize a molecular cloud in the vicinity of the PeVatron candidate LHAASO J2108+5157 by determining its physical parameters from our $^{12,13}$CO($J$=1$\rightarrow$0) line observations. We use an updated estimation of the distance to the cloud, which allows us to obtain a more reliable result. The molecular emission is compared with excess gamma-ray images obtained with Fermi--LAT at energies above 2 GeV to search for spatial correlations and test a possible hadronic ($\pi^0$ decay) origin for the gamma-ray emission. We find that the morphology of the spatial distribution of the CO emission is strikingly similar to that of the Fermi--LAT excess gamma-ray. By combining our observations with archival 21cm HI line data, the nucleons (HI + H$_2$) number density of the target molecular cloud is found to be 133.0 $\pm$ 45.0 cm$^{-3}$, for the measured angular size of 0.55 $\pm$ 0.02$^\circ$ at a distance of 1.6 $\pm$ 0.1 kpc. The resulting total mass of the cloud is M(HI +H$_2$) = 7.5$\pm$2.9$\times$10$^3$ M$_{\odot}$. Under a hadronic scenario, we obtain a total energy of protons of W$_p$ = 4.3$\pm$1.5 $\times$ 10$^{46}$ erg with a cutoff of 700$\pm$300 TeV, which reproduces the sub-PeV gamma-ray emission. We identified a molecular cloud in the vicinity of LHAASO J2107+5157 as the main target where cosmic rays from an unknown PeVatron produce the observed gamma-ray emission via $\pi^0$ decay.