High Energy Astrophysical Phenomena (astro-ph.HE)

Abstract: We have studied the spectro-temporal properties of the neutron star low mass X-ray binary GX 9$+$1 using data from \textit{NuSTAR/FPM} and \textit{AstroSat/SXT} and \textit{LAXPC}. The hardness-intensity diagram of the source showed it to be in the soft spectral state during both observations. \textit{NuSTAR} spectral analysis yielded an inclination angle ($\theta$) $=$ 29$\substack{+3\\-4}^{\circ}$ and inner disk radius ($R_{in}$) $\leq$ 19.01 km. Assuming that the accretion disk was truncated at the Alfv\'en radius during the observation, the upper limit of the magnetic dipole moment ($\mu$) and the magnetic field strength ($B$) at the poles of the neutron star in GX 9$+$1 were calculated to be 1.45$\times$$10^{26}$ G cm$^3$ and 2.08$\times$$10^8$ G, respectively (for $k_A$ $=$ 1). Flux resolved spectral analysis with \textit{AstroSat} data showed the source to be in the soft spectral state ($F_{disk}$/$F_{total}$ $\sim$0.9) with a monotonic increase in mass accretion rate ($\dot{m}$) along the banana branch. The analysis also showed the presence of absorption edges at $\sim$1.9 and $\sim$2.4 keV, likely due to Si XIII and S XV, respectively. Temporal analysis with \textit{LAXPC-20} data in the 0.02 $-$ 100 Hz range revealed the presence of noise components, which could be characterized with broad Lorentzian components.

Abstract: Observations of long duration gamma-ray bursts (GRBs) with TeV emission during their afterglow have been on the rise. Recently, GRB 221009A, the most energetic GRB ever observed, was detected by the {LHAASO} experiment in the energy band 0.2 - 7 TeV. Here, we interpret its afterglow in the context of a hybrid model in which the TeV spectral component is explained by the proton-synchrotron process while the low energy emission from optical to X-ray is due to synchrotron radiation from electrons. We constrained the model parameters using the observed optical, X-ray and TeV data. By comparing the parameters of this burst and of GRB 190114C, we deduce that the VHE emission at energies $\geq$ 1 TeV in the GRB afterglow requires large explosion kinetic energy, $E \gtrsim 10^{54}$~erg and a reasonable circumburst density, $n\gtrsim 10$~cm$^{-3}$. This results in a small injection fractions of particles accelerated to a power-law, $\sim 10^{-2}$. {A significant fraction of shock energy must be allocated to a near equipartition magnetic field, $\epsilon_B \sim 10^{-1}$, while electrons should only carry a small fraction of this energy, $\epsilon_e \sim 10^{-3}$. Under these conditions required for a proton synchrotron model, namely $\epsilon_B \gg \epsilon_e$, the SSC component is substantially sub-dominant over proton-synchrotron as a source of TeV photons.} These results lead us to suggest that proton-synchrotron process is a strong contender for the radiative mechanisms explaining GRB afterglows in the TeV band.

Abstract: Black hole X-ray binaries routinely exhibit Quasi Periodic Oscillations (QPOs) in their Power density spectrum. Studies of QPOs have demonstrated immense ability to understand these dynamical systems although their unambiguous origin still remains a challenge. We investigate the energy-dependent properties of the Type-C QPOs detected for H 1743-322 as observed with AstroSat in its two X-ray outbursts of 2016 and 2017. The combined broadband LAXPC and SXT spectrum is well modelled with a soft thermal and a hard Comptonization component. The QPO exhibits soft/negative lags i.e. variation in soft band lags the variation in hard band, although the upper harmonic shows opposite behaviour i.e. hard/positive lags. Here, we model energy-dependent properties (fractional root mean square and time-lag variation with energy) of the QPO and its upper harmonic individually with a general scheme that fits these properties by utilizing the spectral information and consequently allows to identify the radiative component responsible for producing the variability. Considering the truncated disk picture of accretion flow, a simple model with variation in inner disk temperature, heating rate and fractional scattering with time delays is able to describe the fractional RMS and time-lag spectra. In this work, we show that this technique can successfully describe the energy-dependent features and identify the spectral parameters generating the variability.

Abstract: The air shower simulation code CORSIKA has served as a key part of the simulation chain for numerous astroparticle physics experiments over the past decades. Due to retirement of the original developers and the increasingly difficult maintenance of the monolithic Fortran code of CORSIKA, a new air shower simulation framework has been developed over the course of the last years in C++, called CORSIKA 8. Besides the hadronic and muonic component, the electromagnetic component is one of the key constituents of an air shower. The cascade producing the electromagnetic component of an air shower is driven by bremsstrahlung and photoproduction of electron-positron pairs. At ultrahigh energies or in media with high densities, the bremsstrahlung and pair production processes are suppressed by the Landau-Pomeranchuk-Migdal (LPM) effect, which leads to more elongated showers compared to showers without the LPM suppression. Furthermore, photons at higher energies can produce muon pairs or interact hadronically with nucleons in the target medium, producing a muon component in electromagnetic air showers. In this contribution, we compare electromagnetic showers simulated with the latest Fortran version of CORSIKA and CORSIKA 8, which uses the library PROPOSAL for the electromagnetic component. While earlier validations of CORSIKA 8 electromagnetic showers focused on showers of lower energy, the recent implementation of the LPM effect, photo pair production of muons, and of photohadronic interactions allows now to make a physics-complete comparison also at high energies.

Abstract: We use the light curve and spectral synthesis code JEKYLL to calculate a set of macroscopically mixed Type IIb supernova (SN) models, which are compared to both previously published and new late-phase observations of SN 2020acat. The models differ in the initial mass, the radial mixing and expansion of the radioactive material, and the properties of the hydrogen envelope. The best match to the photospheric and nebular spectra and lightcurves of SN 2020acat is found for a model with an initial mass of 17 solar masses, strong radial mixing and expansion of the radioactive material, and a 0.1 solar mass hydrogen envelope with a low hydrogen mass-fraction of 0.27. The most interesting result is that strong expansion of the clumps containing radioactive material seems to be required to fit the observations of SN 2020acat both in the diffusion phase and the nebular phase. These "Ni bubbles" are expected to expand due to heating from radioactive decays, but the degree of expansion is poorly constrained. Without strong expansion there is a tension between the diffusion phase and the subsequent evolution, and models that fit the nebular phase produce a diffusion peak that is too broad. The diffusion phase lightcurve is sensitive to the expansion of the "Ni bubbles", as the resulting Swiss-cheese-like geometry decreases the effective opacity and therefore the diffusion time. This effect has not been taken into account in previous lightcurve modelling of stripped-envelope SNe, which may lead to a systematic underestimate of their ejecta masses. It should be emphasized, though, that JEKYLL is limited to a geometry that is spherically symmetric on average, and large-scale asymmetries may also play a role. The relatively high initial mass found for the progenitor of SN 2020acat places it at the upper end of the mass distribution of Type IIb SN progenitors, and a single star origin can not be excluded.

Abstract: One of the most challenging open questions regarding the origin of ultrahigh energy cosmic rays (UHECRs) deals with the shape of the source emission spectra. A commonly-used simplifying assumption is that the source spectra of the highest energy cosmic rays trace a Peters cycle, in which the maximum cosmic-ray energy scales linearly with $Z$, i.e., with the charge of the UHECR in units of the proton charge. However, this would only be a natural assumption for models in which UHECRs escape the acceleration region without suffering significant energy losses. In most cases, however, UHECRs interact in the acceleration region and/or in the source environment changing the shape of the source emission spectra. Energy losses are typically parameterized in terms of $Z$ and the UHECR baryon number $A$, and therefore one would expect the source emission spectra to be a function of both $Z$ and $A$. Taking a pragmatic approach, we investigate whether existing data favor any region of the $(Z,A)$ parameter space. Using data from the Pierre Auger Observatory, we carry out a maximum likelihood analysis of the observed spectrum and nuclear composition to shape the source emission spectra for the various particle species. We also study the impact of possible systematic uncertainties driven by hadronic models describing interactions in the atmosphere.

Abstract: We present the extent to which anisotropies in the ultrahigh energy neutrino sky can probe the distribution of extreme astrophysical accelerators in the universe. In this talk, we discuss the origin of an anisotropic neutrino sky and show how observers can use this anisotropy to measure the evolution of ultrahigh energy neutrino sources - and therefore, the sources of ultrahigh energy cosmic rays - for the very first time.

Abstract: The Radio Neutrino Observatory in Greenland (RNO-G) is the only ultrahigh energy (UHE, ${\gtrsim}30$~PeV) neutrino monitor of the Northern sky and will soon be the world's most sensitive high-uptime detector of UHE neutrinos. Because of this, RNO-G represents an important piece of the multimessenger landscape over the next decade. In this talk, we will highlight RNO-G's multimessenger capabilities and its potential to provide key information in the search for the most extreme astrophysical accelerators. In particular, we will highlight opportunities enabled by RNO-G's unique field-of-view, its potential to constrain the sources of UHE cosmic rays, and its complementarity with IceCube at lower energies.

Abstract: New X-ray polarisation results are challenging our understanding of the accretion flow geometry in black hole binary systems. Even spectra dominated by a standard disc can give unexpected results, such as the high inclination black hole binary 4U 1630- 472, where the observed X-ray polarisation is much higher than predicted. This system also shows a strong, highly ionised wind, consistent with thermal-radiative driving from the outer disc, leading to speculation that scattering in the wind is responsible for the unexpectedly high polarisation degree from a standard optically thick disk. Here we show that this is not the case. The optically thin(ish) wind polarises the scattered light in a direction orthogonal to that predicted from a standard optically thick disc, reducing about 2% rather than enhancing the predicted polarisation of the total emission. This value is consistent with the polarisation difference between the disc-dominated soft state, where absorption lines by the wind are clearly seen, and the steep power-law state, where no absorption lines are seen. If this difference is genuinely due to the presence or absence of wind, the total polarisation direction must be orthogonal to the disc plane rather than parallel as expected from optically thick material.

Abstract: We present an extensive $\textit{Hubble Space Telescope}$ ($\textit{HST}$) rest-frame ultraviolet (UV) imaging study of the locations of Type I superluminous supernovae (SLSNe) within their host galaxies. The sample includes 65 SLSNe with detected host galaxies in the redshift range $z\approx 0.05-2$. Using precise astrometric matching with SN images, we determine the distributions of physical and host-normalized offsets relative to the host centers, as well as the fractional flux distribution relative to the underlying UV light distribution. We find that the host-normalized offsets of SLSNe roughly track an exponential disk profile, but exhibit an overabundance of sources with large offsets of $1.5-4$ times their host half-light radius. The SLSNe normalized offsets are systematically larger than those of long gamma-ray bursts (LGRBs), and even Type Ib/c and II SNe. Furthermore, we find that about 40\% of all SLSNe occur in the dimmest regions of their host galaxies (fractional flux of 0), in stark contrast to LGRBs and Type Ib/c and II SNe. We do not detect any significant trends in the locations of SLSNe as a function of redshift, or as a function of explosion and magnetar engine parameters inferred from modeling of their optical lights curves. The significant difference in SLSN locations compared to LGRBs (and normal core-collapse SNe) suggests that at least some of their progenitors follow a different evolutionary path. We speculate that SLSNe arise from massive runaway stars from disrupted binary systems, with velocities of $\sim 10^2$ km s$^{-1}$.