High Energy Astrophysical Phenomena (astro-ph.HE)

Abstract: A catalog of the Galactic population of X-ray pulsars in high-mass X-ray binary (HMXB) systems is presented. It contains information about 82 confirmed sources: 18 persistent and 64 transient pulsars. Their basic parameters include spin period, spin evolution with global and local spin-up/spin-down and duration, orbital period, X-ray luminosity, magnetic field strength measured by cyclotron line analysis, distance, spectral and luminosity class, observable parameters of massive companions, which are shown in the tables provided, with corresponding references. Candidates of the HMXB pulsars are also listed for further careful consideration.

Abstract: We present an analysis of the effects of uncertainties in the atmosphere models on the radius, mass, and other neutron star parameter constraints for the NICER observations of rotation-powered millisecond pulsars. To date, NICER has applied the X-ray pulse profile modeling technique to two millisecond-period pulsars: PSR J0030+0451 and the high-mass pulsar PSR J0740+6620. These studies have commonly assumed a deep-heated fully-ionized hydrogen atmosphere model, although they have explored the effects of partial-ionization and helium composition in some cases. Here we extend that exploration and also include new models with partially ionized carbon composition, externally heated hydrogen, and an empirical atmospheric beaming parametrization to explore deviations in the expected anisotropy of the emitted radiation. None of the studied atmosphere cases have any significant influence on the inferred radius of PSR J0740+6620, possibly due to its X-ray faintness, tighter external constraints, and/or viewing geometry. In the case of PSR J0030+0451 both the composition and ionization state could significantly alter the inferred radius. However, based on the evidence (prior predictive probability of the data), partially ionized hydrogen and carbon atmospheres are disfavored. The difference in the evidence for ionized hydrogen and helium atmospheres is too small to be decisive for most cases, but the inferred radius for helium models trends to larger sizes around or above 14-15 km. External heating or deviations in the beaming that are less than $5\,\%$ at emission angles smaller than 60 degrees, on the other hand, have no significant effect on the inferred radius.

Abstract: The IceCube Neutrino Observatory, a cubic-kilometer-scale neutrino detector at the geographic South Pole, has reached a number of milestones in the field of neutrino astrophysics: the discovery of a high-energy astrophysical neutrino flux, the temporal and directional correlation of neutrinos with a flaring blazar, and a steady emission of neutrinos from the direction of an active galaxy of a Seyfert II type and the Milky Way. The next generation neutrino telescope, IceCube-Gen2, currently under development, will consist of three essential components: an array of about 10,000 optical sensors, embedded within approximately 8 cubic kilometers of ice, for detecting neutrinos with energies of TeV and above, with a sensitivity five times greater than that of IceCube; a surface array with scintillation panels and radio antennas targeting air showers; and buried radio antennas distributed over an area of more than 400 square kilometers to significantly enhance the sensitivity of detecting neutrino sources beyond EeV. This contribution describes the design and status of IceCube-Gen2 and discusses the expected sensitivity from the simulations of the optical, surface, and radio components.

Abstract: In 2019 the NICER collaboration published the first mass and radius inferred for PSR J0030+0451, thanks to NICER observations, and consequent constraints on the equation of state characterising dense matter. Two independent analyses found a mass of $\sim 1.3-1.4\,\mathrm{M_\odot}$ and a radius of $\sim 13\,$km. They also both found that the hot spots were all located on the same hemisphere, opposite to the observer, and that at least one of them had a significantly elongated shape. Here we reanalyse, in greater detail, the same NICER data set, incorporating the effects of an updated NICER response matrix and using an upgraded analysis framework. We expand the adopted models and jointly analyse also XMM-Newton data, which enables us to better constrain the fraction of observed counts coming from PSR J0030+0451. Adopting the same models used in previous publications, we find consistent results, although with more stringent inference requirements. We also find a multi-modal structure in the posterior surface. This becomes crucial when XMM-Newton data is accounted for. Including the corresponding constraints disfavors the main solutions found previously, in favor of the new and more complex models. These have inferred masses and radii of $\sim [1.4 \mathrm{M_\odot}, 11.5$ km] and $\sim [1.7 \mathrm{M_\odot}, 14.5$ km], depending on the assumed model. They display configurations that do not require the two hot spots generating the observed X-rays to be on the same hemisphere, nor to show very elongated features, and point instead to the presence of temperature gradients and the need to account for them.

Abstract: The availability of simultaneous X-ray and Very High Energy (VHE) observations of blazars helps to identify the plausible radiative contributors to the VHE emission. Under leptonic scenario, the VHE emission from BL Lacs are attributed to the synchrotron self Compton (SSC) emission. However, many BL Lacerate (BL Lacs) have shown significant hardening at VHE after correction for the Extra Galactic Background Light (EBL) attenuation. We study the spectral hardening of two nearby BL Lac objects, Mkn 421 and Mkn 501 having most number of simultaneous X-ray and VHE observations available among all the blazars. These BL Lacs are relatively close and the effect of EBL attenuation is relatively minimal/negligible. We study the scatter plot between the X-ray spectral indices and intrinsic VHE indices to identify the plausible origin of the VHE emission. For Mkn 501, the VHE spectral indices are steeper than X-ray spectra, suggesting the scattering process happening at extreme Klein-Nishina regime. On the other hand, for Mkn 421, the VHE spectra is remarkably harder than the X-ray spectra, which suggests an additional emission mechanism other than the SSC process. We show this hard VHE spectrum of Mkn 421 can be explained by considering the inverse Compton (IC) emission from a broken power law electron distribution with Maxwellian pileup. The possibility of the hadronic contribution at VHE {\gamma}-rays is also explored by modeling the hard spectrum under photomeson process.

Abstract: Aims. The main goal of the present paper is to provide the first systematic numerical study of the propagation of astrophysical relativistic jets, in the context of high-resolution shock-capturing resistive relativistic magnetohydrodynamics (RRMHD) simulations. We aim at investigating different values and models for the plasma resistivity coefficient, and at assessing their impact on the level of turbulence, the formation of current sheets and reconnection plasmoids, the electromagnetic energy content, and the dissipated power. Methods. We use the PLUTO code for simulations and we assume an axisymmetric setup for jets, endowed with both poloidal and toroidal magnetic fields, and propagating in a uniform magnetized medium. The gas is assumed to be characterized by a realistic Synge-like equation of state (Taub equation), appropriate for such type of astrophysical jets. The Taub equation is combined here for the first time with the Implicit-Explicit Runge-Kutta time-stepping procedure, as required in RRMHD simulations. Results. The main result is that turbulence is clearly suppressed for the highest values of resistivity (low Lundquist numbers), current sheets are broader, and plasmoids are barely present, while for low values of resistivity results are very similar to ideal runs, where dissipation is purely numerical. We find that recipes employing a variable resistivity based on the advection of a jet tracer or on the assumption of a uniform Lundquist number improve on the use of a constant coefficient and are probably more realistic, preserving the development of turbulence and of sharp current sheets, possible sites for the acceleration of the non-thermal particles producing the observed high-energy emission.

Abstract: We present a new energy transport code that models the time dependent and non-linear evolution of spectra of cosmic-ray nuclei, their secondaries, and photon target fields. The software can inject an arbitrary chemical composition including heavy elements up to iron nuclei. Energy losses and secondary production due to interactions of cosmic ray nuclei, secondary mesons, leptons, or gamma-rays with a target photon field are available for all relevant processes, e.g., photo-meson production, photo disintegration, synchrotron radiation, Inverse Compton scattering, and more. The resulting x-ray fluxes can be fed back into the simulation chain to correct the initial photon targets, resulting in a non-linear treatment of the energy transport. The modular structure of the code facilitates simple extension of interaction or target field models. We will show how the software can be used to improve predictions of observables in various astrophysical sources such as jetted active galactic nuclei (AGN). Since the software can model the propagation of heavy ultrahigh-energy cosmic rays inside the source it can precisely predict the chemical composition at the source. This will also refine predictions of neutrino emissions - they strongly depend on the chemical composition. This helps in the future to optimize the selection and analyses of data from the IceCube neutrino observatory with the aim to enhance the sensitivity of IceCube and reduce the number of trial factors.

Abstract: CRPropa is a Monte Carlo framework for simulating the propagation of (ultra-) high-energy particles in the Universe, including cosmic rays, gamma rays, electrons, and neutrinos. It covers energies from ZeV down to GeV for gamma rays and electrons, and TeV for cosmic rays and neutrinos, supporting various astrophysical environments such as the surroundings of astrophysical sources, galactic, and extragalactic environments. The newest version, CRPropa 3.2, represents a significant leap forward towards a universal multi-messenger framework, opening up the possibility for many more astrophysical applications. This includes extensions to simulate cosmic-ray acceleration and particle interactions within astrophysical source environments, a full Monte Carlo treatment of electromagnetic cascades, improved ensemble-averaged Galactic propagation, significant performance improvements for cosmic-ray tracking through magnetic fields, and a user-friendly implementation of custom photon fields, among many more enhancements. This contribution will give an overview of the new features and present several applications to cosmic-ray and gamma-ray propagation.

Abstract: The radioactive isotopes of 44Ti and 56Ni are important products of explosive nucleosynthesis, which play a key role for supernova (SN) diagnostics and were detected in several nearby young SN remnants. However, most SN models based on non-rotating single stars predict yields of 44Ti that are much lower than the values inferred from observations. We present, for the first time, the nucleosynthesis yields from a self-consistent three-dimensional (3D) SN simulation of an approximately 19 Msun progenitor star that reaches an explosion energy comparable to that of SN 1987A and that covers the evolution of the neutrino-driven explosion until more than 7 seconds after core bounce. We find a significant enhancement of the Ti/Fe yield compared to recent spherically symmetric (1D) models and demonstrate that the long-time evolution is crucial to understand the efficient production of 44Ti due to the non-monotonic temperature and density histories of ejected mass elements. Additionally, we identify characteristic signatures of the nucleosynthesis in proton-rich ejecta, in particular high yields of 45Sc and 64Zn.

Abstract: The detection of high-energy astrophysical neutrinos by IceCube has opened a new window on our Universe. While IceCube has measured the flux of these neutrinos at energies up to several PeV, much remains to be discovered regarding their origin and nature. Currently, measurements are limited by the small sample size of astrophysical neutrinos and by the difficulty of discriminating between electron and tau neutrinos. TAMBO is a next-generation neutrino observatory specifically designed to detect tau neutrinos in the 1-100 PeV energy range, enabling tests of neutrino physics at high energies and the characterization of astrophysical neutrino sources. The observatory will comprise an array of water Cherenkov and plastic scintillator detectors deployed on the face of the Colca Canyon in the Peruvian Andes. This unique geometry will facilitate a high-purity measurement of astrophysical tau neutrino properties. In this talk, I will present the prospects of TAMBO in the context of next-generation neutrino observatories and provide an overview of its current status.

Abstract: Origin of wide varieties of quasi-periodic oscillation (QPO) observed in compact sources is still not well established. Its frequencies range from mHz to kHz spanning all compact objects. Are different QPOs, with different frequencies, originating from different Physics? We propose that the emergence of QPOs is the result of nonlinear resonance of fundamental modes present in accretion disks forced by external modes including that of the spin of the underlying compact object. Depending on the properties of accreting flow, e.g. its velocity and gradient, resonances, and any mode locking, take place at different frequencies, exhibiting low to high frequency QPOs. We explicitly demonstrate the origin of higher frequency QPOs for black holes and neutron stars by a unified model and outline how the same physics could be responsible to produce lower frequency QPOs. The model also predicts the spin of black holes, and constrains the radius of neutron stars and the mass of both.

Abstract: AT2022cmc was recently reported as the first on-axis jetted tidal disruption event (TDE) discovered in the last decade, and the fourth on-axis jetted TDE candidate known so far. In this work, we present NuSTAR hard X-ray (3--30 keV) observations of AT2022cmc, as well as soft X-ray (0.3--6 keV) observations obtained by NICER, Swift, and XMM-Newton. Our analysis reveals that the broadband X-ray spectra can be well described by a broken power-law with $f_\nu \propto \nu^{-0.5}$ ($f_\nu \propto \nu^{-1}$) below (above) the rest-frame break energy of $E_{\rm bk}\sim 10$ keV at observer-frame $t_{\rm obs}=7.8$ and 17.6 days since discovery. At $t_{\rm obs} = 36.2$ days, the X-ray spectrum is consistent with either a single power-law of $f_\nu \propto \nu^{-0.8}$ or a broken power-law with spectral slopes similar to the first two epochs. By modeling the spectral energy distribution evolution from radio to hard X-ray across the three NuSTAR observing epochs, we find that the sub-millimeter/radio emission originates from external shocks at large distances $\gtrsim\! 10^{17}$ cm from the black hole, the UV/optical light comes from a thermal envelope with radius $\sim\!10^{15}$ cm, and the X-ray emission is consistent with synchrotron radiation powered by energy dissipation at intermediate radii within the (likely magnetically dominated) jet. Our interpretation differs from the model proposed by Pasham et al. (2023) where both the radio and X-rays come from the same emitting zone in a matter-dominated jet. Our model for the jet X-ray emission has broad implications on the nature and radiation mechanism of relativistic jets in other sources such as gamma-ray bursts.