High Energy Astrophysical Phenomena (astro-ph.HE)

Abstract: We present the spectroscopy of the neutron star low-mass X-ray binary 4U 1636-53 using six simultaneous XMM-Newton and Rossi X-ray Timing Explorer observations. We applied different self-consistent reflection models to explore the features when the disk is illuminated by either the corona or the neutron star surface. We found that the spectra could be well fitted by these two types of models, with the derived emissivity index below a typical value of 3. The relative low emissivity can be explained if the neutron star and the corona, working together as an extended illuminator, simultaneously illuminate and ionize the disk. Additionally, the derived ionization parameter in the lamppost geometry is larger than the theoretical prediction. This inconsistency likely suggests that the corona does not emit isotropically in a realistic context. Furthermore, we also found that there is a possible trend between the height of the corona and the normalization of the disk emission. This could be understood either as a variation in the reflected radiation pressure or in the context of a jet base. Finally, we found that the disk is less ionized if it is illuminated by the neutron star, indicating that the illuminating source has significant influence on the physical properties of the disk.

Abstract: Star-forming galaxies (SFGs) have been established as an important source population in the extra-galactic $\gamma$-ray background (EGB). Their intensive star-formation creates an abundance of environments able to accelerate particles, and these build-up a rich sea of cosmic rays (CRs). Above GeV energies, CR protons can undergo hadronic interactions with their environment to produce $\gamma$-rays. SFGs can operate as CR proton "calorimeters", where a large fraction of the CR energy is converted to $\gamma$-rays. However, CRs also deposit energy and momentum to modify the thermal and hydrodynamic conditions of the gas in SFGs, and can become a powerful driver of outflows. Such outflows are ubiquitous among some types of SFGs, and have the potential to severely degrade their CR proton calorimetry. This diminishes their contribution to the EGB. In this work, we adopt a self-consistent treatment of particle transport in outflows from SFGs to assess their calorimetry. We use 1D numerical treatments of galactic outflows driven by CRs and thermal gas pressure, accounting for the dynamical effects and interactions of CRs. We show the impact CR-driven flows have on the relative contribution of SFG populations to the EGB, and investigate the properties of SFGs that contribute most strongly.

Abstract: Dark matter could decay into Standard Model particles producing neutrinos directly or indirectly. The resulting flux of neutrinos from these decays could be detectable at neutrino telescopes and would be associated with massive celestial objects where dark matter is expected to be accumulated. Recent observations of high-energy astrophysical neutrinos at IceCube might hint at a signal produced by the decay of TeV to PeV scale dark matter. This analysis searches for neutrinos from decaying dark matter in nearby galaxy clusters and galaxies. We focus on dark matter masses from 10 TeV to 1 EeV and four decay channels: $\nu\bar{\nu}$, $\tau^{+}\tau^{-}$, $W^{+}W^{-}$, $b\bar{b}$. Three galaxy clusters, seven dwarf galaxies, and the Andromeda galaxy are chosen as targets and stacked within the same source class. A well-established IceCube data sample is used, which contains 11 years of upward-going track-like events. In this contribution, we present preliminary results of the analysis.

Abstract: We present the results of photometric and spectroscopic monitoring campaigns of the changing look AGN NGC~2617 carried out from 2016 until 2022 and covering the wavelength range from the X-ray to the near-IR. The facilities included the telescopes of the SAI MSU, MASTER Global Robotic Net, the 2.3-m WIRO telescope, Swift, and others. We found significant variability at all wavelengths and, specifically, in the intensities and profiles of the broad Balmer lines. We measured time delays of ~ 6 days (~ 8 days) in the responses of the H-beta (H-alpha) line to continuum variations. We found the X-ray variations to correlate well with the UV and optical (with a small time delay of a few days for longer wavelengths). The K-band lagged the B band by 14 +- 4 days during the last 3 seasons, which is significantly shorter than the delays reported previously by the 2016 and 2017--2019 campaigns. Near-IR variability arises from two different emission regions: the outer part of the accretion disc and a more distant dust component. The HK-band variability is governed primarily by dust. The Balmer decrement of the broad-line components is inversely correlated with the UV flux. The change of the object's type, from Sy1 to Sy1.8, was recorded over a period of ~ 8 years. We interpret these changes as a combination of two factors: changes in the accretion rate and dust recovery along the line of sight.

Abstract: By comparing the physical properties of pulsars hosted by core-collapsed (CCed) and non-core-collapsed (Non-CCed) globular clusters (GCs), we find that pulsars in CCed GCs rotate significantly slower than their counterparts in Non-CCed GCs. Additionally, radio luminosities at 1.4 GHz in CCed GCs are higher. These findings are consistent with the scenario that dynamical interactions in GCs can interrupt angular momentum transfer processes and surface magnetic field decay during the recycling phase. Our results suggest that such effects in CCed GCs are stronger due to more frequent disruptions of compact binaries. This is further supported by the observation that both estimated disruption rates and the fraction of isolated pulsars are predominantly higher in CCed GCs.

Abstract: We present explicit expressions for Rayleigh and Raman scattering cross-sections and phase matrices of the ground $1s$ state hydrogen atom based on the Kramers-Heisenberg dispersion formula. The Rayleigh scattering leaves the hydrogen atom in the ground-state while the Raman scattering leaves the hydrogen atom in either $ns$ ($n\geq2$; $s$-branch) or $nd$ ($n\geq3$; $d$-branch) excited state, and the Raman scattering converts incident ultraviolet (UV) photons around the Lyman resonance lines into optical-infrared (IR) photons. We show that this Raman wavelength conversion of incident flat UV continuum in dense hydrogen gas with a column density of $N_{\text{H}} > 10^{21}~\text{cm}^{-2}$ can produce broad emission features centred at Balmer, Paschen, and higher-level lines, which would mimic Doppler-broadened hydrogen lines with the velocity width of $\gtrsim 1,000~\text{km}~\text{s}^{-1}$ that could be misinterpreted as signatures of Active Galactic Nuclei, supernovae, or fast stellar winds. We show that the phase matrix of the Rayleigh and Raman $s$-branch scatterings is identical to that of the Thomson scattering while the Raman $d$-branch scattering is more isotropic, thus the Paschen and higher-level Raman features are depolarized compared to the Balmer features due to the flux contribution from the Raman $d$-branch. We argue that observations of the line widths, line flux ratios, and linear polarization of multiple optical/IR hydrogen lines are crucial to discriminate between the Raman-scattered broad emission features and Doppler-broadened emission lines.

Abstract: The spin orientations of spinning binary black hole (BBH) mergers detected by ground-based gravitational wave detectors such as LIGO and Virgo can provide important clues about the formation of such binaries. However, these spin tilts, i.e., the angles between the spin vector of each black hole and the binary's orbital angular momentum vector, can change due to precessional effects as the black holes evolve from a large separation to their merger. The tilts inferred at a frequency in the sensitive band of the detectors by comparing the signal with theoretical waveforms can thus be significantly different from the tilts when the binary originally formed. These tilts at the binary's formation are well approximated in many scenarios by evolving the BBH backwards in time to a formally infinite separation. Using the tilts at infinite separation also places all binaries on an equal footing in analyzing their population properties. In this paper, we perform parameter estimation for simulated BBHs and investigate the differences between the tilts one infers directly close to merger and those obtained by evolving back to infinite separation. We select simulated observations such that their configurations show particularly large differences in their orientations close to merger and at infinity. While these differences may be buried in the statistical noise for current detections, we show that in future plus-era (A$+$ and Virgo$+$) detectors, they can be easily distinguished in some cases. We also consider the tilts at infinity for BBHs in various spin morphologies and at the endpoint of the up-down instability. In particular, we find that we are able to easily identify the up-down instability cases as such from the tilts at infinity.

Abstract: When modeling the population of merging binary black holes, analyses generally focus on characterizing the distribution of primary (i.e. more massive) black holes in the binary, while simplistic prescriptions are used for the distribution of secondary masses. However, the secondary mass distribution provides a fundamental observational constraint on the formation history of coalescing binary black holes. If both black holes experience similar stellar evolutionary processes prior to collapse, as might be expected in dynamical formation channels, the primary and secondary mass distributions would show similar features. If they follow distinct evolutionary pathways (for example, due to binary interactions that break symmetry between the initially more massive and less massive star), their mass distributions may differ. We explicitly fit the secondary mass distribution, finding that if the primary and secondary mass distributions are different, the previously-identified peak in the primary mass distribution may be driven by an even larger peak in the secondary mass distribution. Alternatively, if we assume that the two masses are drawn from the same underlying distribution, they both show a peak at $31.4_{-2.6}^{+2.3} \, M_{\odot}$. This value is shifted lower than that obtained when assuming the peak only exists in the marginal primary mass distribution, placing this feature in further tension with expectations from a pulsational pair-instability supernova pileup. We anticipate that by the end of the fifth LIGO-Virgo-KAGRA observing run, we will be able to determine whether the data prefer distinct or identical component mass distributions to $>4\sigma$, providing important clues to the formation history of coalescing binary black holes.

Abstract: Throughout the past three decades, only a few tens of observations have been made of optical flashes contemporaneous with gamma-ray bursts (GRBs), despite the thousands of GRBs that have been detected during that same timeframe. In this work, we present light curves from the Transiting Exoplanet Survey Satellite (TESS) for a sample of 7 GRBs that were localized to within 10" by the Swift X-ray Telescope. For each burst, we characterize both the prompt emission, if it exists, and the afterglow, and conduct searches for late-time emission from a supernova or kilonova component. We also constrain the physical parameters of the burst based on the TESS light curve, and present a novel method to account for the effects of TESS's cosmic ray mitigation strategy on the observed flux from these GRBs. This allows us to establish upper limits on the true magnitude of any GRB-associated optical flash. Finally, we discuss how TESS's continuous monitoring and new weekly downlink schedule are proving invaluable in the rapid follow-up and characterization of short-duration transients, including GRBs; these could potentially enable TESS to detect electromagnetic counterparts to gravitational-wave events.

Abstract: Eclipsing of the X-ray emitting region in active galactic nuclei (AGN) is a potentially powerful probe to examine the AGN environment and absorber properties. Here we study the eclipse data from the 2016 XMM-Newton observation of NGC 6814 using a colour-colour analysis. Colours (i.e. hardness ratios) can provide the advantage of better time-resolution over spectral analysis alone. Colour-colour grids are constructed to examine the effects of different parameters on the observed spectral variability during the eclipse. Consistent with previous spectral analysis, the variations are dominated by changes in the column density and covering fraction of the absorber. However, during maximum eclipse the behaviour of the absorber changes. Just after ingress, the eclipse is described by changes in column density and covering fraction, but prior to egress, the variations are dominated by changes in column density alone. Simulations are carried out to consider possible absorber geometries that might produce this behaviour. The behaviour is inconsistent with a single, homogeneous cloud, but simulations suggest that multiple clouds, perhaps embedded in a highly ionised halo, could reproduce the results. In addition, we determine the orbital covering factor (fraction of orbital path-length) based on evidence of several eclipses in the 2016, 64-day Swift light curve. We estimate that ~ 2-4 per cent of the orbit is covered by obscuring clouds and that the distribution of clouds is not isotropic.

Abstract: Only a tiny fraction ~ 1% of stellar tidal disruption events (TDE) generate powerful relativistic jets evidenced by luminous hard X-ray and radio emissions. We propose that a key property responsible for both this surprisingly low rate and a variety of other observations is the typically large misalignment {\psi} between the orbital plane of the star and the spin axis of the supermassive black hole (SMBH). Such misaligned disk/jet systems undergo Lense-Thirring precession together about the SMBH spin axis. We find that TDE disks precess sufficiently rapidly that winds from the accretion disk will encase the system on large scales in a quasi-spherical outflow. We derive the critical jet efficiency {\eta} > {\eta}crit for both aligned and misaligned precessing jets to successfully escape from the disk-wind ejecta. As {\eta}crit is higher for precessing jets, less powerful jets only escape after alignment with the SMBH spin. Alignment can occur through magneto-spin or hydrodynamic mechanisms, which we estimate occur on typical timescales of weeks and years, respectively. The dominant mechanism depends on {\eta} and the orbital penetration factor \b{eta}. Hence depending only on intrinsic parameters of the event {{\psi},{\eta},\b{eta}}, we propose that each TDE jet can either escape prior to alignment, thus exhibiting erratic X-ray light curve and two-component radio afterglow (e.g., Swift J1644+57) or escape after alignment. Relatively rapid magneto-spin alignments produce relativistic jets exhibiting X-ray power-law decay and bright afterglows (e.g., AT2022cmc), while long hydrodynamic alignments give rise to late jet escape and delayed radio flares (e.g., AT2018hyz).

Abstract: Tidal disruption events (TDEs) are exotic transients that can lead to temporary super-Eddington accretion onto a supermassive black hole. Such accretion mode is naturally expected to result in powerful outflows of ionized matter. However, to date such an outflow has only been directly detected in the X-ray band in a single TDE, ASASSN-14li. This outflow has a low velocity of just a few 100 km/s, although there is also evidence for a second, ultra-fast phase. Here we present the detection of a low-velocity outflow in a second TDE, ASASSN-20qc. The high-resolution X-ray spectrum reveals an array of narrow absorption lines, each blueshifted by a few 100 km/s, which cannot be described by a single photo-ionization phase. For the first time, we confirm the multiphase nature of a TDE outflow, with at least two phases and two distinct velocity components. One highly ionized phase is outflowing at $910^{+90}_{-80}$ km/s, while a lower ionization component is blueshifted by $400_{-120}^{+100}$ km/s. We perform time-resolved analysis of the X-ray spectrum and detect that, surprisingly, the mildly ionized absorber strongly varies in ionization parameter over the course of a single 60 ks observation, indicating that its distance from the black hole may be as low as 400 gravitational radii. We discuss these findings in the context of TDEs and compare this newly detected outflow with that of ASASSN-14li.