2024-12-04 | | Total: 19
GRB 221009A is one of the brightest transients ever observed with the highest peak gamma-ray flux for a gamma-ray burst (GRB). A type Ic-BL supernova (SN), SN 2022xiw, was definitively detected in late-time JWST spectroscopy (t = 195 days, observer-frame). However, photometric studies have found SN 2022xiw to be less luminous (10-70%) than the canonical GRB-SN, SN 1998bw. We present late-time Hubble Space Telescope (HST)/WFC3 and JWST/NIRCam imaging of the afterglow and host galaxy of GRB 221009A at t ~ 185, 277, and 345 days post-trigger. Our joint archival ground, HST, and JWST light curve fits show strong support for a break in the light curve decay slope at t = 50 +/- 10 days (observer-frame) and a supernova at $1.4^{+0.37}_{-0.40} \times$ the optical/NIR flux of SN 1998bw. This break is consistent with an interpretation as a jet break when requiring slow-cooling electrons in a wind medium with the electron energy spectral index, p > 2, and $\nu_m < \nu_c$. Our light curve and joint HST/JWST spectral energy distribution (SED) also show evidence for the late-time emergence of a bluer component in addition to the fading afterglow and supernova. We find consistency with the interpretations that this source is either a young, massive, low-metallicity star cluster or a scattered light echo of the afterglow with a SED shape of $f_{\nu} \propto \nu^{2.0\pm1.0}$.
The supersonic flow motions associated with infall of baryonic gas toward sheets and filaments, as well as cluster mergers, produces large-scale shock waves. The shocks associated with galaxy clusters can be classified mainly into two categories: internal shocks appear in the hot intracluster medium within the viral radius, and external accretion shocks form in the outer cold region well outside of the virial radius. Cosmic-ray (CR) electrons and/or protons accelerated by these shocks are expected to produce gamma-rays through inverse-Compton scatterings (ICS) or inelastic $pp$ collisions respectively. Recent studies have found a spatially extended GeV source within the virial radius, consistent with the internal shock origin. Here we report the detection of a new GeV source at a distance of about 2.8$^\circ$ from the center of the Coma cluster through the analysis of 16.2 years of Fermi-LAT data. The hard spectrum of the source, in agreement with the ICS origin, and its location in a large-scale filament of galaxies points to the external accretion shock origin. The gamma-ray ($0.1-10^3$ GeV) luminosity of the source, $1.4\times 10^{42}~ {\rm erg~s^{-1}}$, suggests that a fraction $\sim 10^{-3}$ of the kinetic energy flux through the shock-surface is transferred to relativistic CR electrons.
Context. Radio galaxies with visible two-sided jet structures, such as NGC 1052, are sources of particular interest to study the collimation and shock structure of active galactic nuclei jets. High-resolution very-long-baseline interferometry observations of such sources can resolve and study the jet collimation profile and probe different physical mechanisms. Aims. In this paper, we study the physics of double-sided radio sources at parsec scales, and in particular investigate whether propagating shocks can give rise to the observed asymmetry between jet and counterjet. Methods. We carry out special relativistic hydrodynamic simulations and perform radiative transfer calculations of an over-pressured perturbed jet. During the radiative transfer calculations we incorporate both thermal and nonthermal emission while taking the finite speed of light into account. To further compare our results to observations, we create more realistic synthetic data including the properties of the observing array as well as the image reconstruction via multifrequency regularized maximum likelihood methods. We finally introduce a semi-automatized method for tracking jet components and extracting jet kinematics. Results. We show that propagating shocks in an inherently symmetric double-sided jet can lead to partially asymmetric jet collimation profiles due to time delay effects and relativistic beaming. These asymmetries may appear on specific epochs, with one jet evolving near conically and the other one parabolically (width profile evolving with a slope of 1 and 0.5, respectively). However, these spurious asymmetries are not significant when observing the source evolve for an extended amount of time. Conclusions. Purely observational effects are not enough to explain a persisting asymmetry in the jet collimation profile of double-sided jet sources and hint at evidence for asymmetrically launched jets.
Rotation-powered pulsars represent the main class of identified gamma-ray sources in the Galaxy. The wealth of observational data collected by the AGILE and Fermi gamma-ray space telescopes in the GeV range, and by ground-based Cherenkov telescopes in the TeV band provide invaluable insights into how relativistic plasmas dissipate and accelerate particles. Decoding the information contained in the gamma-ray pulses profile is an important step to understand how pulsars work. In this study, we aim at putting an ab initio plasma model of pulsar magnetospheres to the test, in light of the most recent gamma-ray observations in the GeV and TeV bands. To this end, we present of a new series of global particle-in-cell simulations of an inclined pulsar magnetosphere. High-quality synthetic pulse profiles in the synchrotron and inverse Compton channels are reconstructed to study in greater details their morphology and their energy dependence. We also perform a fit of observed lightcurves with the model, using the third Fermi-LAT gamma-ray pulsar catalog. Reconnection in the wind current sheet powers synchrotron and inverse Compton emission. The modeled pulse profiles reproduce some of the salient features of observed gamma-ray pulsars, including the mysterious Vela-like lightcurves, such as: the generic double-peaked structure, the presence of a bridge or third peak in between the main pulses, the pulse narrowing with increasing energy. The bolometric synchrotron radiative efficiency is strictly limited by the reconnection rate. Our global kinetic simulations are able to match observed pulse profiles. Such direct comparisons will help drive and focus future simulation developments.
We conduct the first-ever Pulsar Polarization Array (PPA) analysis to detect the ultralight Axion-Like Dark Matter (ALDM) using the polarization data of 22 millisecond pulsars from the third data release of Parkes Pulsar Timing Array. As one of the major dark matter candidates, the ultralight ALDM exhibits a pronounced wave nature on astronomical scales and offers a promising solution to small-scale structure issues within local galaxies. While the linearly polarized pulsar light travels through the ALDM galactic halo, its position angle (PA) can be subject to an oscillation induced by the ALDM Chern-Simons coupling with electromagnetic field. The PPA is thus especially suited for detecting the ultralight ALDM by correlating polarization data across the arrayed pulsars. To accomplish this task, we develop an advanced Bayesian analysis framework that allows us to construct pulsar PA residual time series, model noise contributions properly and search for pulsar cross-correlations. We find that for an ALDM density of $\rho_0=0.4\,\textrm{GeV}/\textrm{cm}^3$, the Parkes PPA offers the best global limits on the ALDM Chern-Simons coupling, namely $\lesssim 10^{-13.5}-10^{-12.2}~{\rm GeV}^{-1}$, for the mass range of $10^{-22} - 10^{-21}~{\rm eV}$. The crucial role of pulsar cross-correlation in recognizing the nature of the derived limits is also highlighted.
We present the statistic results of GeV spectral breaks of bright gamma-ray flat-spectrum radio quasars (FSRQs) in the energy range of 0.1-10 GeV based on New Pass 8 date of the Large Area Telescope abroad on Fermi Gamma-ray Space Telescope. We have fitted the 15-year average gamma-ray spectra of 755 FSRQs by using both a broken power-law (BPL) and the Logarithmic Parabolas (LP) models, and obtained 87 bright gamma-ray FSRQs with their integrated photon fluxes greater than $2.16\times 10^{-8}$ cm$^{-2}$ s$^{-1}$. From our results, the FSRQ population shows similar preferences for both the BPL and LP models in gamma-ray spectral fitting and clustering analysis suggests that BPL-preferred and LP-preferred FSRQs belong to the same category. Our results indicate that GeV spectral breaks in bright gamma-ray FSRQs are located at $\rm 2.90\pm1.92$ GeV in the rest frame, and the observed change in photon index is $\rm \Delta \gamma =0.45 \pm 0.19$, which is consistent with the expected value for a cooling break of electrons scattering seed photons.
Data derived from general relativistic MHD simulations of accretion onto black holes can be used as input to a post-processing scheme that predicts the radiated spectrum. Combining a relativistic Compton-scattering radiation transfer solution in the corona with detailed local atmosphere solutions incorporating local ionization and thermal balance within the disk photosphere, it is possible to study both spectral formation and intrinsic spectral variability in the radiation from relativistic accretion disks. With this method, we find that radiatively-efficient systems with black holes of $10M_\odot$ accreting at $\approx 0.01$ in Eddington units produce spectra very similar to those observed in the hard states of X-ray binaries. The spectral shape above 10keV is well-described by a power-law with an exponential cutoff. Intrinsic turbulent variations lead to order-unity fluctuations in bolometric luminosity, logarithmic spectral slope fluctuations over a range $\approx 0.15$, and factor of 2 changes in the cut-off energy on timescales $\sim 50\, (M_{\rm BH}/10 M_\odot)$ms. Within the corona, the range of gas temperature spans more than an order of magnitude. The wide distribution of temperatures is central to defining the spectrum: the logarithmic spectral slope is harder by $\sim 0.3$ and the cut-off energy larger by a factor $\sim 10 - 30$ than if the coronal temperature everywhere were its mass-weighted mean.
Binary black hole (BBH) evolution in the discs of active galactic nuclei (AGN) is a promising channel for gravitational wave (GW)-driven mergers. It is however unclear whether binaries interacting with the surrounding disc undergo orbital contraction or expansion. We develop a simple analytic model of accreting BBHs in AGN discs to follow the orbital evolution from the disc-dominated regime at large separations into the GW-driven regime at small separations (the coupled `disc+GW'-driven evolution). We obtain that accreting binaries expand in thick discs with aspect ratio greater than a critical value ($> h_\mathrm{crit}$); whereas accreting binaries contract and eventually merge in thin discs ($< h_\mathrm{crit}$). Interestingly, accreting BBHs can experience faster mergers compared to non-accreting counterparts, with a non-monotonic dependence on the disc aspect ratio. The orbital contraction is usually coupled with eccentricity growth in the disc-dominated regime, which lead to accelerated inspirals in the GW-driven regime. We quantify the resulting BBH merger timescales in AGN discs ($\tau_\mathrm{merger} \sim 10^5 - 10^7$ yr) and estimate the associated GW merger rates ($\mathcal{R} \sim (0.2 - 5) \, \text{Gpc}^{-3} \text{yr}^{-1}$). Overall, accreting binaries may efficiently contract and merge in thin discs, hence this particular BBH-in-AGN channel may provide a non-negligible contribution to the observed GW merger event rate.
Observations of tidal disruption events (TDEs) on a timescale of years after the main flare show evidence of continued activity in the form of optical/UV emission, quasi-periodic eruptions, and delayed radio flares. Motivated by this, we explore the time evolution of these disks using semi-analytic models to follow the changing disk properties and feeding rate to the central black hole (BH). We find that thermal instabilities typically begin $\sim150-250\,{\rm days}$ after the TDE, causing the disk to cycle between high and low accretion states for up to $\sim10-20\,{\rm yrs}$. The high state is super-Eddington, which may be associated with outflows that eject $\sim10^{-3}-10^{-1}\,M_\odot$ with a range of velocities of $\sim0.03-0.3c$ over a span of a couple of days and produce radio flares. In the low state, the accretion rate slowly grows over many months to years as continued fallback accretion builds the mass of the disk. In this phase, the disk may reach luminosities of $\sim10^{41}-10^{42}\,{\rm erg\,s^{-1}}$ in the UV as seen in some late-time observations. We highlight the importance of the iron-opacity "bump" at $\approx2\times10^5\,{\rm K}$ in generating sufficiently high luminosities. This work suggests that joint optical/UV observations with radio monitoring could be key for following the disk state as the radio flares are produced.
Metallicity plays a crucial role in the evolution of massive stars and their final core-collapse supernova (CCSN) explosions. Integral-field-unit (IFU) spectroscopy can provide a spatially resolved view of SN host galaxies and serve as a powerful tool to study SN metallicities. Early transient surveys targeted bright galaxies with high star formation and SN rates; as a result, the discovered SNe are significantly biased toward high metallicities. More recently, the untargeted, wide-field transient surveys, such as ASAS-SN and ZTF, have discovered a large number of SNe without such a bias. In this work, we construct a large and unbiased sample of SNe discovered by (quasi-) untargted searches, consisting of 209 SNe of Types II (with unknown subtypes), IIP, IIn, IIb, Ib and Ic at z $\leq$ 0.02 with VLT/MUSE observations. This is currently the largest CCSN sample with IFU observations. With the strong-line method, we reveal the spatially-resolved metallicity maps of the SN host galaxies and acquire accurate metallicity measurements for the SN sites. Our results show that the SN metallicities range from 12 + log(O/H) = 8.1 to 8.7 dex, and the metallicity distributions for different SN types are very close to each other, with mean and median values of 8.4 - 8.5 dex. We carefully analysed the stochastic sampling effect, showing that our large sample size narrows the 1$\sigma$ uncertainty down to only 0.05 dex. The apparent metallicity differences among SN types are all within 3$\sigma$ uncertainties and the metallicity distributions for different SN types are all consistent with being randomly drawn from the same reference distribution. This suggests that metallicity plays a minor role in the origin of different CCSN types and some other metallicity-insensitive processes, such as binary interaction, dominate the distinction of CCSN types.
Many questions posed in the Astro2020 Decadal survey in both the New Messengers and New Physics and the Cosmic Ecosystems science themes require a gamma-ray mission with capabilities exceeding those of existing (e.g. Fermi, Swift) and planned (e.g. COSI) observatories. ComPair, the Compton Pair telescope, is a prototype of such a next-generation gamma-ray mission. It had its inaugural balloon flight from Ft. Sumner, New Mexico in August 2023. To continue the goals of the ComPair project to develop technologies that will enable a future gamma-ray mission, the next generation of ComPair (ComPair-2) will be upgraded to increase the sensitivity and low-energy transient capabilities of the instrument. These advancements are enabled by AstroPix, a silicon monolithic active pixel sensor, in the tracker and custom dual-gain silicon photomultipliers and front-end electronics in the calorimeter. This effort builds on design work for the All-sky Medium Energy Gamma-ray Observatory eXplorer (AMEGO-X) concept that was submitted the 2021 MIDEX Announcement of Opportunity. Here we describe the ComPair-2 prototype design and integration and testing plans to advance the readiness level of these novel technologies.
Ultra-light axion-like particles (ALPs) can be a viable solution to the dark matter problem. The scalar field associated with ALPs, coupled to the electromagnetic field, acts as an active birefringent medium, altering the polarisation properties of light through which it propagates. In particular, oscillations of the axionic field induce monochromatic variations of the plane of linearly polarised radiation of astrophysical signals. The radio emission of millisecond pulsars provides an excellent tool to search for such manifestations, given their high fractional linear polarisation and negligible fluctuations of their polarisation properties. We have searched for the evidence of ALPs in the polarimetry measurements of pulsars collected and preprocessed for the European Pulsar Timing Array (EPTA) campaign. Focusing on the twelve brightest sources in linear polarisation, we searched for an astrophysical signal from axions using both frequentist and Bayesian statistical frameworks. For the frequentist analysis, which uses Lomb-Scargle periodograms at its core, no statistically significant signal has been found. The model used for the Bayesian analysis has been adjusted to accommodate multiple deterministic systematics that may be present in the data. A statistically significant signal has been found in the dataset of multiple pulsars with common frequency between $10^{-8}$ Hz and $2\times10^{-8}$ Hz, which can most likely be explained by the residual Faraday rotation in the terrestrial ionosphere. Strong bounds on the coupling constant $g_{a\gamma}$, in the same ballpark as other searches, have been obtained in the mass range between $6\times10^{-24}$ eV and $5\times10^{-21}$ eV. We conclude by discussing problems that can limit the sensitivity of our search for ultra-light axions in the polarimetry data of pulsars, and possible ways to resolve them.
Accurate $N$-body simulations of multiple systems such as binaries and triples are essential for understanding the formation and evolution of interacting binaries and binary mergers, including gravitational wave sources, blue stragglers and X-ray binaries. The logarithmic time-transformed explicit symplectic integrator (LogH), also known as algorithmic regularization, is a state-of-the-art method for this purpose.However, we show that this method is accurate for isolated Kepler orbits because of its ability to trace Keplerian trajectories, but much less accurate for hierarichal triple systems. The method can lead to an unphysical secular evolution of inner eccentricity in Kozal-Lidov triples, despite a small energy error. We demonstrate that hybrid methods, which apply LogH to the inner binary and alternative methods to the outer bodies, are significantly more effective, though not symplectic. Additionally, we introduce a more efficient hybrid method, BlogH, which eliminates the need for time synchronization and is time symmetric. The method is implemented in the few-body code SDAR. We explore suitable criteria for switching between the LogH and BlogH methods for general triple systems. These hybrid methods have the potential to enhance the integration performance of hierarchial triples.
Stellar-mass black holes descend from high-mass stars, most of which had stellar binary companions. However, the number of those binary systems that survive the binary evolution and black hole formation is uncertain by multiple orders of magnitude. The survival rate is particularly uncertain for massive stars with low-mass companions, which are thought to be the progenitors of most black hole X-ray binaries. We present a search for close black hole companions (separations less than 20 solar radii) to AFGK-type stars in TESS, i.e. the non-accreting counterparts to and progenitors of low-mass X-ray binaries. Such black holes can be detected by the tidally induced ellipsoidal deformation of the visible star, and the ensuing photometric light-curve variations. From an initial sample of 4.7 million TESS stars, we have selected 457 candidates for such variations. However, spectroscopic followup of 251 of them shows that none are consistent with a close black hole companion. On the basis of this non-detection, we determine (2 $\sigma$ confidence) that fewer than one in $10^5$ Solar-type stars in the Solar neighbourhood host a short-period black hole companion. This upper limit is in tension with a number of ``optimistic'' population models in the literature that predict short-period black hole companions around one in $10^{4-5}$ stars. Our limits are still consistent with other models that predict only a few in $10^{7-8}$.
The James Webb Space Telescope (JWST) has revealed an unexpectedly high abundance of UV luminous galaxies at redshifts $z\gtrsim 10$, challenging `standard' galaxy formation models. This study investigates the role of rapidly rotating (massive) stars undergoing chemically homogeneous evolution (CHE) in reconciling this potential tension. These stars are more compact, hotter, and exhibit enhanced UV emission. We find that the rest-frame UV luminosity of star-forming galaxies can be significantly enhanced by a factor of $\sim 3-6$ when CHE stars above a minimum initial mass of $m_{\star,\min}^{\rm CHE}\sim 2-10\ \rm M_\odot$ account for more than half of the total stellar mass following a Salpeter initial mass function. As a result, the UV luminosity functions observed at $z\sim 12-16$ can be reproduced with less extreme values of star formation efficiency and UV luminosity stochastic variability. Our results highlight the potential of CHE in explaining the UV-bright galaxy populations detected by JWST and call for future work to explore the broader astrophysical implications of CHE and its associated phenomena in the early universe, such as gamma-ray bursts, compact object binaries, and metal enrichment.
Horizon-scale imaging with the Event Horizon Telescope (EHT) has provided transformative insights into supermassive black holes but its resolution and scope are limited by ground-based constraints such as the size of the Earth, its relatively slow rotation, and atmospheric delays. Space-based very long baseline interferometry (VLBI) offers the capability for studying a larger and more diverse sample of black holes. We identify a number of nearby supermassive black holes as prime candidates for horizon-scale imaging at millimeter wavelengths, and use source characteristics such as angular size, sky distribution, and variability timescales to shape the design of a space-based array. We identify specific metrics that serve as key predictors of image fidelity and scientific potential, providing a quantitative basis for optimizing mission design parameters. Our analysis demonstrates that the optimal configuration requires two space-based elements in high Earth orbits (HEO) that are not coplanar and are apparently counter-rotating. Our results delineate the key requirements for a space-based VLBI mission, enabling detailed studies of black hole shadows, plasma dynamics, and jet formation, advancing black hole astrophysics beyond the current capabilities of the EHT.
A large population of binary systems in the Universe emitting gravitational waves (GW) would produce a stochastic noise, known as the gravitational wave background (GWB). The properties of the GWB directly depend on the attributes of its constituents. If the binary systems are in eccentric orbits, it is well established that the GW power they radiate strongly depends on their instantaneous orbital phase. Consequently, their power spectrum varies over time, and the resulting GWB can appear nonstationary. In this work, we estimate the amplitude of time-dependent fluctuations in the GWB power spectrum as a function of the eccentricity of the binaries. Specifically, we focus on the GWB produced by a population of supermassive black hole binaries (SMBHB) that should be observable by pulsar timing arrays (PTA). We show that a large population of homogeneously distributed equal SMBHBs produces nonstationary features that are undetectable by current PTA datasets. However, using more realistic and astrophysically motivated populations of SMBHBs, we show that the nonstationarity might become very large and detectable, especially in the case of more massive and eccentric populations. In particular, when one binary is slightly brighter than the GWB, we demonstrate that time fluctuations can become significant. This is also true for individual binary systems with a low signal-to-noise ratio (SNR) relative to the GWB (SNR $\approx$ 1), which standard data analysis methods would struggle to detect. The detection of nonstationary features in the GWB could indicate the presence of some relatively bright GW sources in eccentric orbits, offering new insights into the origins of the signal.
Superradiance, the process by which waves are amplified through energy and angular momentum transfer, can also occur in horizonless objects like boson stars, due to both the real space and internal field space rotations. In this work, we study superradiance in the frequency and time domains for static and spinning boson stars, constructed within general relativity and with a self-interacting complex scalar field as a matter source. Using linear perturbation theory and three dimensional nonlinear simulations, we calculate amplification factors and analyze energy and angular momentum transfer in scattering processes, with results showing consistency between approaches. Wave scattering inside a cavity containing a boson star is also examined, demonstrating the effects of confinement on amplification.
In this work, we analyze the characteristics of electromagnetic (EM) radiation associated with scalar and axion field oscillations in different background field setups. Because the scalar field and axion field have different parity and couple with the EM field in different forms, the EM signals generated by the scalar and axion can be used to distinguish them. More interestingly, resonance effects amplify the difference between the two fields and consequent EM signal strength, which helps us distinguish and detect them in future observations.