High Energy Astrophysical Phenomena (astro-ph.HE)

Abstract: We report the detection of millihertz quasi-periodic oscillations ($\mathrm{mHz}$ QPOs) from the low-mass X-ray binary XTE 1701$-$462. The discovery came from a search of the legacy data set of the Rossi X-ray Timing Explorer, in order to detect the periodic signals in all observations of sources exhibiting thermonuclear bursts. We found that $47$ out of $860$ observations of XTE 1701$-$462; covering the 2006--7 outburst exhibits signals with a significance above the detection threshold, which was determined separately for each observation via a Monte Carlo approach. We chose the four strongest candidates, each with maximum power exceeding $4\sigma$ of the simulated wavelet noise power distribution, to demonstrate the properties of the QPOs. The frequencies of the signals in the four observations are $\sim 3.5\;\text{to}\;5.6\; \mathrm{mHz}$, and the fractional R.M.S. amplitudes vary between $0.74 \pm 0.05\,\%$ and $3.54 \pm 0.04\,\%$. Although previously reported signals in other sources typically disappear immediately before a burst, we do not observe this behaviour in XTE 1701$-$462. Instead, we found that the QPOs and bursts occurred in separate accretion regimes. When the persistent luminosity dropped near the end of the outburst, the source showed bursts and no QPOs were detected, which is the behaviour predicted by theory for the transition from stable to unstable burning. On the basis of this new detection, we reassess the cases for identifying these $\mathrm{mHz}$ QPOs in this and other sources as arising from marginally stable burning.

Abstract: We present an investigation into the spectroscopic properties of non-equilibrium photoionization processes operating in a time-evolving mode. Through a quantitative comparison between equilibrium and time-evolving models, we find that the time-evolving model exhibits a broader distribution of charge states compared to the equilibrium model, accompanied by a slight shift in the peak ionization state depending on the source variability and gas density. The time-evolving code, tpho in SPEX, has been successfully employed to analyze the spectral properties of warm absorbers in the Seyfert galaxy NGC 3783. The incorporation of variability in the tpho model improves the fits of the time-integrated spectra, providing more accurate descriptions to the average charge states of several elements, in particular for Fe which is peaked around Fe XIX. The inferred densities and distances of the relevant X-ray absorber components are estimated to be approximately a few 1E11 per cubic meter and less than 1 pc, respectively. Furthermore, the updated fit suggests a potential scenario in which the observed absorbers are being expelled from the central AGN at the escape velocities. This implies that these absorbers might not play a significant role in the AGN feedback mechanism.

Abstract: Some supernovae (SNe), such as Type IIn SNe, are powered by collision of the SN ejecta with a dense circumstellar matter (CSM). Their emission spectra show characteristic line shapes of combined broad emission and narrow P-Cyg lines, which should closely relate to the CSM structure and the mass-loss mechanism that creates the dense CSM. We quantitatively investigate the relationship between the line shape and the CSM structure by Monte Carlo radiative transfer simulations, considering two representative cases of dense CSM formed by steady and eruptive mass loss. Comparing the H$\alpha$ emission between the two cases, we find that a narrow P-Cyg line appears in the eruptive case while it does not appear in the steady case, due to the difference in the velocity gradient in the dense CSM. We also reproduce the blue-shifted photon excess observed in some SNe IIn, which is formed by photon transport across the shock wave and find the relationship between the velocity of the shocked matter and the amount of the blue shift of the photon excess. We conclude that the presence or absence of narrow P-Cyg lines can distinguish the mass loss mechanism, and suggest high-resolution spectroscopic observations with $\lambda/ \Delta \lambda \gtrsim 10^4$ after the light curve peak for applying this diagnostic method.

Abstract: Current simulations of air showers produced by ultra-high energy cosmic rays (UHECRs) do not satisfactorily describe recent experimental data, particularly when looking at the muonic shower component relative to the electromagnetic one. Discrepancies can be seen in both average values and on an individual shower-by-shower basis. It is thought that the muonic part of the air showers isn't accurately represented in simulations, despite various attempts to boost the number of muons within standard hadronic interaction physics. In this study, we investigate whether modifying the final state of events created with Sibyll~2.3d in air shower simulations can achieve a more consistent description of the muon content observed in experimental data. We create several scenarios where we separately increase the production of baryons, $\rho^0$, and strange particles to examine their impact on realistic air shower simulations. Our results suggest that these ad-hoc modifications can improve the simulations, providing a closer match to the observed muon content in air showers. One side-effect of the increased muon production in the considered model versions is a smaller difference in the predicted total muon numbers for proton and iron showers. However, more research is needed to find out whether any of these adjustments offers a realistic solution to the mismatches seen in data, and to identify the precise physical process causing these changes in the model. We hope that these modified model versions will also help to develop improved machine-learning analyses of air shower data and to estimate sys.{} uncertainties related to shortcomings of hadronic interaction models.

Abstract: In this paper, we present an updated version of our model (KYNXiltr) which considers thermal reverberation of a standard Novikov-Thorne accretion disc illuminated by an X-ray point-like source. Previously, the model considered only two cases of black hole spins, and assumed a colour correction factor $f_{\rm col} = 2.4$. Now, we extend the model to any spin value and colour correction. In addition, we consider two scenarios of powering the X-ray corona, either via accretion, or external to the accretion disc. We use KYNXiltr to fit the observed time lags obtained from intense monitoring of four local Seyfert galaxies (NGC 5548, NGC 4395, Mrk 817, and Fairall 9). We consider various combinations of black hole spin, colour correction, corona height, and fraction of accretion power transferred to the corona. The model fits well the overall time-lags spectrum in these sources (for a large parameter space). For NGC 4593 only, we detect a significant excess of delays in the U-band. The contribution of the diffuse BLR emission in the time-lags spectrum of this source is significant. It is possible to reduce the large best-fitting parameter space by combining the results with additional information, such as the observed Eddington ratio and average X-ray luminosity. We also provide an update to the analytic expression provided by Kammoun et al., for an X-ray source that is not powered by the accretion process, which can be used for any value of colour correction, and for two values of the black hole spin (0 and 0.998).

Abstract: De-excitation gamma-ray lines, produced by nuclei colliding with protons, provide information about astrophysical environments where particles have kinetic energies of $10-100$ MeV per nucleon. In general, such environments can be categorized into two types: the interaction between non-thermal MeV cosmic rays and ambient gas, and the other is thermal plasma with a temperature above a few MeV. In this paper, we focus on the latter type and investigate the production of de-excitation gamma-ray lines in very hot thermal plasma, especially the dependence of the line profile on the plasma temperature. We have calculated the line profile of prompt gamma rays from $^{12}$C and $^{16}$O and found that when nuclei have a higher temperature than protons, gamma-ray line profiles can have a complex shape unique to each nucleus species. This is caused by anisotropic gamma-ray emission in the nucleus rest frame. We propose that the spectroscopy of nuclear de-excitation gamma-ray lines may enable to probe energy distribution in very hot astrophysical plasmas. This diagnostics can be a new and powerful technique to investigate the physical state of a two-temperature accretion flows onto a black hole, especially the energy distributions of the protons and nuclei, which are difficult to access for any other diagnostics.

Abstract: Late-time spectra of Type Ia supernovae (SNe Ia) are important in clarifying the physics of their explosions, as they provide key clues to the inner structure of the exploding white dwarfs. We examined late-time optical spectra of 36 SNe Ia, including five from our own project (SNe 2019np, 2019ein, 2021hpr, 2021wuf, and 2022hrs), with phase coverage of $\sim 200$ to $\sim 400$ days after maximum light. At this late phase, the outer ejecta have become transparent and the features of inner iron-group elements emerge in the spectra. Based on multicomponent Gaussian fits and reasonable choices for the pseudocontinuum around Ni and Fe emission features, we get reliable estimates of the Ni to Fe ratio, which is sensitive to the explosion models of SNe Ia. Our results show that the majority (about 67%) of our SNe Ia are more consistent with the sub-Chandrasekhar-mass (i.e., double-detonation) model, although they could be affected by evolutionary or ionisation effects. Moreover, we find that the Si II $\lambda$6355 velocity measured around the time of maximum light tends to increase with the Ni to Fe ratio for the subsample with either redshifted or blueshifted nebular velocities, suggesting that progenitor metallicity might play an important role in accounting for the observed velocity diversity of SNe Ia.

Abstract: Recent analysis of the kilonova, AT2017gfo, has indicated that this event was highly spherical. This may challenge hydrodynamics simulations of binary neutron star mergers, which usually predict a range of asymmetries, and radiative transfer simulations show a strong direction dependence. Here we investigate whether the synthetic spectra from a 3D kilonova simulation of asymmetric ejecta from a hydrodynamical merger simulation can be compatible with the observational constraints suggesting a high degree of sphericity in AT2017gfo. Specifically, we determine whether fitting a simple P-Cygni line profile model leads to a value for the photospheric velocity that is consistent with the value obtained from the expanding photosphere method. We would infer that our kilonova simulation is highly spherical at early times, when the spectra resemble a blackbody distribution. The two independently inferred photospheric velocities can be very similar, implying a high degree of sphericity, which can be as spherical as inferred for AT2017gfo, demonstrating that the photosphere can appear spherical even for asymmetrical ejecta. The last-interaction velocities of radiation escaping the simulation show a high degree of sphericity, supporting the inferred symmetry of the photosphere. We find that when the synthetic spectra resemble a blackbody the expanding photosphere method can be used to obtain an accurate luminosity distance (within 4-7 per cent).

Abstract: Horizon-scale observations of the jetted active galactic nucleus M87 are compared with simulations spanning a broad range of dissipation mechanisms and plasma content in three-dimensional general relativistic flows around spinning black holes. Observations of synchrotron radiation from radio to X-ray frequencies can be compared with simulations by adding prescriptions specifying the relativistic electron-plus-positron distribution function and associated radiative transfer coefficients. A suite of time-varying simulations with various spins and plasma magnetizations is chosen to represent distinct possibilities for the M87 jet/accretion flow/black hole (JAB) system. We then input turbulent heating and equipartition-based emission prescriptions (and piecewise combinations thereof) in the time-dependent 3D simulations, in which jet morphology, polarization and variation are "observed" and compared with real observations so as to try to infer the rules that govern the polarized emissivity. The models in this paper support a magnetically arrested disk (MAD) with several possible spin/emission model combinations supplying the jet in M87, whose inner jet and black hole shadow have been observed down to the photon ring at 230 GHz by the Event Horizon Telescope (EHT). We also show that some MAD cases that are dominated by intrinsic circular polarization have near-linear V/I dependence on unpaired electron or positron content while SANE polarization exhibits markedly greater positron-dependent Faraday effects -- future probes of the SANE/MAD dichotomy and plasma content with the EHT. This is the second work in a series also applying the "observing" simulations methodology to near-horizon regions of supermassive black holes in Sgr A* and 3C 279.