Date: Wed, 8 May 2024 | Total: 17
A wide variety of Galactic sources show transient emission at soft and hard X-ray energies: low-mass and high-mass X-ray binaries containing compact objects (e.g., novae, microquasars, transitional millisecond pulsars, supergiant fast X-ray transients), isolated neutron stars exhibiting extreme variability as magnetars as well as pulsar wind nebulae. Although most of them can show emission up to MeV and/or GeV energies, many have not yet been detected in the TeV domain by Imaging Atmospheric Cherenkov Telescopes. In this paper, we explore the feasibility of detecting new Galactic transients with the Cherenkov Telescope Array (CTA) and the prospects for studying them with Target of Opportunity observations. We show that CTA will likely detect new sources in the TeV regime, such as the massive microquasars in the Cygnus region, low-mass X-ray binaries with low- viewing angle, flaring emission from the Crab pulsar-wind nebula or other novae explosions, among others. We also discuss the multi-wavelength synergies with other instruments and large astronomical facilities.
We examine expected effective spin ($\chi_{{\rm eff},1\rm yr}$) and orbital eccentricity ($e_{1\rm yr}$) correlations for a population of observable equal-mass massive black hole binaries (MBHBs) with total redshifted mass $M_z\sim[10^{4.5},10^{7.5}]~{\rm M}_\odot$ embedded in a circumbinary disc (CBD), one-year before merging in the LISA band. We find a strong correlation between measurable eccentricity and negative effective spin for MBHBs that are carried to merger by retrograde accretion. This is due to the well-established eccentricity pumping of retrograde accretion and the formation of retrograde CBD-aligned mini-discs, as observed in hydrodynamical simulations. Conversely, prograde accretion channels result in positive $\chi_{{\rm eff},1\rm yr}$ and non-measurable $e_{1\rm yr}$. This clear contrast between the two CBD orientations -- and particularly the unique signature of retrograde configurations -- provides a promising way to unlock the mysteries of MBHB formation channels in the LISA era.
Pulsar wind nebulae (PWNe) are clouds of the magnetized relativistic electron/positron plasma supplied from the central pulsar. However, the number of radio-emitting particles inside a PWN is larger than the expectation from the study of pulsar magnetospheres and then their origin is still unclear. A stochastic acceleration of externally injected particles by a turbulence inside the PWN is proposed by our previous studies. In this paper, the previous stochastic acceleration model of the PWN broadband spectra is improved by taking into account the time evolution of the turbulent energy and then the total energy balance inside a PWN is maintained. The turbulent energy supplied from the central pulsar is wasted by the backreaction from the stochastic particle acceleration and the adiabatic cooling according the PWN expansion. The model is applied to the Crab Nebula and reproduce the current broadband emission spectrum, especially the flat radio spectrum although time evolution of the turbulent energy (diffusion coefficient) is a bit complicated compared with our previous studies, where we assumed an exponential behavior of the diffusion coefficient.
The study of temporal properties of variable sources can elucidate their physical processes. In this context, we present a critical study comparing three approaches to periodic or quasiperiodic behavior: Gaussian process, power spectrum, and wavelet analysis, using celerite, Lomb-Scargle periodograms, and weighted wavelet-Z transforms, respectively. We use 15 Swift-X-ray Telescope light curves of short gamma-ray bursts (sGRBs) as examples. A comprehensive analysis of two sGRB X-ray light curves is performed. The results reveal the importance of artifacts, largely in the form of false quasiperiodic oscillation signals, possibly introduced by preprocessing (such as detrending) or other aspects of the analysis. The exploration described in this paper can be helpful for future studies of variability in GRBs, active galactic nuclei, and other astronomical sources.
JWST observed the magnetar 4U 0142+61 with the MIRI and NIRCam instruments within a 77 min time interval on 2022 September 20-21. The low-resolution MIRI spectrum and NIRCam photometry show that the spectrum in the wavelength range 1.4-11 $\mu$m range can be satisfactorily described by an absorbed power-law model, $f_{\nu}\propto \nu^{-\alpha}$, with a spectral slope $\alpha =0.96\pm0.02$, interstellar extinction $A_V= 3.9\pm0.2$, and normalization $f_0 = 59.4\pm 0.5$ $\mu$Jy at $\lambda = 8$ $\mu$m. These observations do not support the passive disk model proposed by Wang et al. (2006), based on the Spitzer photometry, which was interpreted as evidence for a fallback disk from debris formed during the supernova explosion. We suggest a nonthermal origin for this emission and source variability as the most likely cause of discrepancies between the JWST data and other IR-optical observing campaigns. However, we cannot firmly exclude the presence of a large disk with a different dependence of the effective disk temperature on distance from the magnetar. Comparison with the power-law fit to the hard X-ray spectrum above 10 keV, measured by NuSTAR contemporaneously with JWST, shows that the X-ray spectrum is significantly harder. This may imply that the X-ray and IR nonthermal emission come from different sites in the magnetosphere of the magnetar.
The mergers of double neutron stars (DNSs) systems are believed to drive the majority of short $\gamma$-ray bursts (SGRBs), while also serving as production sites of heavy r-process elements. Despite being key to i) confirming the nature of the extragalactic SGRBs, ii) addressing the poorly-understood r-process enrichment in the ultra-faint dwarf galaxies (UFDGs), and iii) probing the formation process of DNS systems, the space velocity distribution of DNSs is still poorly constrained due to the small number of DNSs with well-determined astrometry. In this work, we determine new proper motions and parallaxes of two Galactic DNSs -- PSR J0509+3801 and PSR J1930-1852, using the Very Long Baseline Array, and estimate the transverse velocities $v_\perp$ of all the 11 isolated Galactic DNSs having proper motion measurements in a consistent manner. Our correlation analysis reveals that the DNS $v_\perp$ is tentatively correlated with three parameters: spin period, orbital eccentricity, and companion mass. With the preliminary $v_\perp$ distribution, we obtain the following findings. Firstly, the refined $v_\perp$ distribution is confirmed to agree with the observed displacements of the localized SGRBs from their host galaxy birth sites. Secondly, we estimate that around 11% and 25% of DNSs remain gravitationally bound to UFDGs with escape velocities of 15$\mathrm{~km~s^{-1}}$ and 25$\mathrm{~km~s^{-1}}$, respectively. Hence, the retained DNSs might indeed be responsible for the r-process enrichment confirmed so far in a few UFDGs. Finally, we discuss how a future ensemble of astrometrically determined DNSs may probe the multimodality of the $v_\perp$ distribution.
We present the MOST Hosts survey (Multi-Object Spectroscopy of Transient Hosts). The survey is planned to run throughout the five years of operation of the Dark Energy Spectroscopic Instrument (DESI) and will generate a spectroscopic catalog of the hosts of most transients observed to date, in particular all the supernovae observed by most public, untargeted, wide-field, optical surveys (PTF/iPTF, SDSS II, ZTF, DECAT, DESIRT). Scientific questions for which the MOST Hosts survey will be useful include Type Ia supernova cosmology, fundamental plane and peculiar velocity measurements, and the understanding of the correlations between transients and their host galaxy properties. Here, we present the first release of the MOST Hosts survey: 21,931 hosts of 20,235 transients. These numbers represent 36% of the final MOST Hosts sample, consisting of 60,212 potential host galaxies of 38,603 transients (a transient can be assigned multiple potential hosts). Of these galaxies, 40% do not appear in the DESI primary target list and therefore require a specific program like MOST Hosts. Of all the transients in the MOST Hosts list, only 26.7% have existing classifications, and so the survey will provide redshifts (and luminosities) for nearly 30,000 transients. A preliminary Hubble diagram and a transient luminosity-duration diagram are shown as examples of future potential uses of the MOST Hosts survey. The survey will also provide a training sample of spectroscopically observed transients for photometry-only classifiers, as we enter an era when most newly observed transients will lack spectroscopic classification. The MOST Hosts DESI survey data will be released through the Wiserep platform on a rolling cadence and updated to match the DESI releases. Dates of future releases and updates are available through the https://mosthosts.desi.lbl.gov website.
This study investigates the origins of GW230529, delving into its formation from massive stars within isolated binary systems. Utilizing population synthesis models, we present compelling evidence that the neutron star component forms first. However, the event's low signal-to-noise ratio introduces complexities in identifying the underlying physical mechanisms driving its formation. Augmenting our analysis with insights from numerical relativity, we estimate the final black hole mass and spin to be approximately $5.3 M_\odot$ and $0.53$, respectively. Furthermore, we employ the obtained posterior samples to calculate the ejecta mass and kilonova light curves resulting from r-process nucleosynthesis. We find the ejecta mass to range within $0-0.06 M_{\odot}$, contingent on the neutron star equation of state. The peak brightness of the kilonovae light curves indicates that targeted follow-up observations with a Rubin-like observatory may have detected this emission.
Galactic PeVatrons are sources that can accelerate cosmic rays to PeV energies. The high-energy cosmic rays are expected to interact with the surrounding ambient material or radiation, resulting in the production of gamma rays and neutrinos. To optimize for the detection of such associated production of gamma rays and neutrinos for a given source morphology and spectrum, a multi-messenger analysis that combines gamma rays and neutrinos is required. In this study, we use the Multi-Mission Maximum Likelihood framework (3ML) with IceCube Maximum Likelihood Analysis software (i3mla) and HAWC Accelerated Likelihood (HAL) to search for a correlation between 22 known gamma-ray sources from the third HAWC gamma-ray catalog and 14 years of IceCube track-like data. No significant neutrino emission from the direction of the HAWC sources was found. We report the best-fit gamma-ray model and 90% CL neutrino flux limit from the 22 sources. From the neutrino flux limit, we conclude that the gamma-ray emission from five of the sources can not be produced purely from hadronic interactions. We report the limit for the fraction of gamma rays produced by hadronic interactions for these five sources.
In an effort to search for faint sources of emission over arbitrary timescales, we present a novel method for analyzing forced photometry light curves in difference imaging from optical surveys. Our method "ATLAS Clean'' or ATClean, utilizes the reported fluxes, uncertainties, and fits to the point-spread function from difference images to quantify the statistical significance of individual measurements. We apply this method to control light curves across the image to determine whether any source of flux is present in the data for a range of specific timescales. From ATLAS $o$-band imaging at the site of the Type II supernova (SN) 2023ixf in M101 from 2015--2023, we show that this method accurately reproduces the 3$\sigma$ flux limits produced from other, more computationally expensive methods. We derive limits for emission on timescales of 5~days and 80-300~days at the site of SN\,2023ixf, which are 19.8 and 21.3~mag, respectively. The latter limits rule out variability for unextinguished red supergiants (RSG) with initial masses $>$22~$M_{\odot}$, comparable to the most luminous predictions for the SN 2023ixf progenitor system. We also compare our limits to short timescale outbursts, similar to those expected for Type IIn SN progenitor stars or the Type II SN 2020tlf, and rule out outburst ejecta masses of $>$0.021~$M_{\odot}$, much lower than the inferred mass of circumstellar matter around SN 2023ixf in the literature. In the future, these methods can be applied to any forced point-spread function photometry on difference imaging from other surveys, such as Rubin optical imaging.
We report the ASKAP discovery of a new Galactic supernova remnant (SNR) candidate G308.73+1.38, which we name Raspberry. This new SNR candidate has an angular size of 20.7 arcmin $\times$ 16.7 arcmin, and we measure a total integrated flux of 407$\pm$50 mJy. We estimate Raspberry's most likely diameter of 10$-$30 pc which would place it at a distance of 3$-$5 kpc, on the near side of the Milky Way's Scutum$-$Centaurus Arm. We also find a Stokes$-$V point source close to the centre of Raspberry with a $\sim$5$\sigma$ significance. This point source may be the remaining compact source, a neutron star, or possibly a pulsar, formed during the initial supernova event.
Li {\it et al.\/} (2024) reported a 4.605 day period in the repeating FRB 20121102A in addition to the previously reported 157 day modulation of its activity. This note suggests that the shorter period is the orbital period of a mass-transferring star orbiting a black hole, possibly of intermediate mass, and that the 157 day period is the precession period of an accretion disc around the black hole. The mass-losing star must be evolved.
GW170817 and GRB 170817A provided direct evidence that binary neutron star (NSNS) mergers can produce short gamma-ray bursts (sGRBs). However, questions remain about the nature of the central engine. Depending on the mass, the remnant from a NSNS merger may promptly collapse to a black hole (BH), form a hypermassive neutron star (HMNS) which undergoes a delayed collapse to a BH, a supramassive neutron star (SMNS) with a much longer lifetime, or an indefinitely stable NS. There is strong evidence that a BH with an accretion disk can launch a sGRB-compatible jet via the Blandford-Znajek mechanism, but whether a supramassive star can do the same is less clear. We have performed general relativistic magnetohydrodynamics (GRMHD) simulations of the merger of both irrotational and spinning, equal-mass NSNSs constructed from a piecewise polytropic representation of the SLy equation of state, with a range of gravitational (ADM) masses that yield remnants with mass above and below the supramassive limit. Each NS is endowed with a dipolar magnetic field extending from the interior into the exterior, as in a radio pulsar. We examine cases with different initial binary masses, including a case which produces a HMNS which collapses to a BH, and lower mass binaries that produce SMNS remnants. We find similar jet-like structures for both the SMNS and HMNS remnants that meet our basic critera for an incipient jet. The outflow for the HMNS case is consistent with a Blandford-Znajek (BZ) jet. There is sufficient evidence that such BZ-powered outflows can break out and produce ulrarelativistic jets so that we can describe the HMNS system as a sGRB progenitor. However, the incipient jets from the SMNS remnants have much more baryon pollution and we see indications of inefficient acceleration and mixing with the surrounding debris. Therefore, we cannot conclude that SMNS outflows are the progenitors of sGRBs.
There are two contradictory views of the eROSITA bubbles: either a 10 kpc-scale pair of giant bubbles blown by the Galactic center (GC), or a 100 pc-scale local structure coincidentally located in the direction of GC. A key element of this controversy is the distance to the bubbles. Based on the 3D dust distribution in the Galactic plane, we found three isolated, distant (500-800 pc) clouds at intermediate Galactic latitudes. Their projected morphologies perfectly match the X-ray shadows on the defining features of the north eROSITA bubble, i.e., the North Polar Spur (NPS) and the Lotus Petal Cloud (LPC), indicating that both the NPS and LPC are distant with a distance lower limit of nearly 1kpc. In the X-ray dark region between the NPS and LPC, we found a few polarized radio arcs and attributed them to the bubble's shock front. These arcs match up perfectly with the outer border of the NPS and LPC and provide a way to define the bubble's border. The border defined in this way can be well described by the line-of-sight tangent of a 3D skewed cup model rooted in the GC. We conclude that, instead of being two independent, distant features, NPS and LPC compose a single, giant bubble, which, therefore, is most plausibly a 10-kpc scale bubble rooted at the GC.
Modern photometric surveys of the sky suggest that many, perhaps most supernovae (SNe) associated with the explosion of massive stars are influenced at an appreciable level by their interaction with circumstellar material (CSM). The photometric and spectroscopic diversity of these transients point to a wide range of CSM properties in terms of mass, extent, composition, and location relative to the exploding star, suggesting progenitors that cover from standard to the most extreme mass loss rates. Surveys at high-cadence catch massive stars at shock breakout and inform us on the immediate mass loss history before core collapse. In contrast, long-term monitoring of these SNe cover the transition to the birth of a SN remnant and document the progenitor mass loss that took place centuries to millennia before explosion. Interacting SNe are therefore not just extraordinary astrophysical laboratories to study radiation-dominated shocks and probe the distant Universe, they also open the path to novel and fundamental studies on stellar evolution, stellar stability, or mass loss in single and binary massive stars.
Millisecond pulsars with white dwarf companions have typical eccentricities $e\sim 10^{-6}-10^{-3}$. The eccentricities of helium white dwarfs are explained well by applying the fluctuation-dissipation theorem to convective eddies in their red giant progenitors. We extend this theory to more massive carbon-oxygen (CO) white dwarfs with asymptotic giant branch (AGB) progenitors. Due to the radiation pressure in AGB stars, the dominant factor in determining the remnant white dwarf's eccentricity is the critical residual hydrogen envelope mass $m_{\rm env}$ required to inflate the star to giant proportions. Using a suite of MESA stellar evolution simulations with $\Delta m_{\rm c}=10^{-3}\,{\rm M}_\odot$ core-mass intervals, we resolved the AGB thermal pulses and found that the critical $m_{\rm env}\propto m_{\rm c}^{-6}$. This steep dependence causes the $e(m_{\rm c})$ relation to turn over, such that $e\sim 3\times 10^{-3}$ almost independently of the remnant CO white dwarf's mass $m_{\rm c}$. Nearly all of the measured eccentricities lie below this robust theoretical limit, indicating that the eccentricity is damped during the common-envelope inspiral that follows the unstable Roche-lobe overflow of the AGB star. Specifically, we focused on white dwarfs with median masses $m_{\rm c}>0.6\,{\rm M}_\odot$. These massive white dwarfs begin their inspiral with practically identical orbital periods and eccentricities, eliminating any dependence on the initial conditions. For this sub-sample, we find an empirical relation $e\propto P^{3/2}$ between the final period and eccentricity that is much tighter than previous studies - motivating theoretical work on the eccentricity evolution during the common envelope phase.
We present a multi-wavelength analysis, from the radio to the X-ray band, of the redshift $z=6.44$ VIK J2318$-$31 radio-loud (RL) quasi stellar object (QSO), one of the most distant currently known in this class. The work is based on newly obtained (uGMRT, ATCA, Chandra) as well as archival (GNIRS and X-Shooter) dedicated observations that have not been published yet. Based on the observed X-ray and radio emission, its relativistic jets are likely young and misaligned from our line of sight. Moreover, we can confirm, with simultaneous observations, the presence of a turnover in the radio spectrum at $\nu_{\rm peak} \sim 650$ MHz which is unlikely to be associated with self-synchrotron absorption. From the NIR spectrum we derived the mass of the central black hole, M$_{\rm BH}=8.1^{+6.8}_{-5.6} \times 10^8 {\rm M_{\odot}}$, and the Eddington ratio, $\lambda_{\rm EDD} = 0.8^{+0.8}_{-0.6}$, using broad emission lines as well as an accretion disc model fit to the continuum emission. Given the high accretion rate, the presence of a $\sim$8$\times$10$^8$ M$_\odot$ black hole at $z=6.44$ can be explained by a seed black hole ($\sim$10$^{4}$ M$_\odot$) that formed at $z\sim25$, assuming a radiative efficiency $\eta_{\rm d}\sim0.1$. However, by assuming $\eta_{\rm d}\sim0.3$, as expected for jetted systems, the mass observed would challenge current theoretical models of black hole formation.