Astrophysics

In Jiang et al. (2020), we reported a possible bright flash (hereafter GN-z11-flash) from a galaxy GN-z11 at z ~ 11. Recently, Steinhardt et al. (2021; arXiv:2101.12738) found 27 images with transient signals in Keck MOSFIRE archival data and claimed that GN-z11-flash was more likely from a moving object in our Solar system. We show that the Steinhardt et al.'s definition of the chance probability and their methodology of finding GN-z11-flash-like transients are problematic in several aspects. In particular, none of their transients is analogous to GN-z11-flash, and none of them is positionally coincident with a known object in their imaging data. In Jiang et al., we performed a comprehensive analysis of the origin of GN-z11-flash and ruled out, to the best of our knowledge, the possibility of known man-made objects or moving objects in the Solar system, based on all available information and our current understanding of these objects. Steinhardt et al. did not use such information and did not analyse the GN-z11-flash event itself. The majority of their transients are apparently low-Earth orbit satellites or aircrafts. Therefore, their analysis can neither prove nor disprove our results. Finally, we present a method to estimate the chance probability of finding GN-z11-flash-like transients in archival data. Based on this method and the archival data used by Steinhardt et al., we obtain a loose upper limit of the probability that actually support the original results of Jiang et al. (2020).

Superluminous supernovae (SLSNe) are massive star explosions too luminous to be powered by traditional energy sources, such as radioactive 56Ni. These transients may instead be powered by a central engine, such as a millisecond pulsar or magnetar, whose relativistic wind inflates a nebula of high energy particles and radiation behind the expanding ejecta. We present 3D Monte Carlo radiative transfer calculations which follow the production and thermalization of high energy radiation from the nebula into optical radiation and, conversely, determine the gamma-ray emission that escapes the ejecta without thermalizing. We track the evolution of photons and matter in a coupled two-zone ("wind/nebula" and "ejecta") model, accounting for the range of radiative processes. We identify a novel mechanism by which gamma-gamma pair creation in the upstream pulsar wind regulates the mean energy of particles entering the nebula over the first several years after the explosion, rendering our results on this timescale insensitive to the (uncertain) intrinsic wind pair multiplicity. To explain the observed late-time steepening of SLSNe optical light curves as being the result of gamma-ray leakage, the nebular magnetization must be very low, epsB <~ 1e-6-1e-4. For higher epsB, synchrotron emission quickly comes to dominate the thermalized nebula radiation, and being readily absorbed because of its lower photon energies, results in the SN optical light curve tracking the spin-down power even to late times >~ 1 yr, inconsistent with observations. For magnetars to remain viable contenders for powering SLSNe, we conclude that either magnetic dissipation in the wind/nebula is extremely efficient, or that the spin-down luminosity decays significantly faster than the canonical dipole rate ~1/t^2 in a way that coincidentally mimicks gamma-ray escape.

Several large neutrino telescopes, operating at various sites around the world, have as their main objective the first detection of neutrinos emitted by a gravitational collapse in the Milky Way. The success of these observation programs depends on the rate of supernova core collapse in the Milky Way, R . In this work, standard statistical techniques are used to combine several independent results. Their consistency is discussed and the most critical input data are identified. The inference on R is further tested and refined by including direct information on the occurrence rate of gravitational collapse events in the Milky Way and in the Local Group, obtained from neutrino telescopes and electromagnetic surveys. A conservative treatment of the errors yields a combined rate R=1.63±0.46 (100 yr) −1 ; the corresponding time between core collapse supernova events turns out to be T= 61 +24 −14 ~yr. The importance to update the analysis of the stellar birthrate method is emphasized.

The study of high-energy gamma rays from passive Giant Molecular Clouds (GMCs) in our Galaxy is an indirect way to characterize and probe the paradigm of the "sea" of cosmic rays in distant parts of the Galaxy. By using data from the High Altitude Water Cherenkov (HAWC) observatory, we measure the gamma-ray flux above 1 TeV of a set of these clouds to test the paradigm. We selected high-galactic latitude clouds that are in HAWC's field-of-view and which are within 1~kpc distance from the Sun. We find no significant excess emission in the cloud regions, nor when we perform a stacked log-likelihood analysis of GMCs. Using a Bayesian approach, we calculate 95\% credible intervals upper limits of the gamma-ray flux and estimate limits on the cosmic-ray energy density of these regions. These are the first limits to constrain gamma-ray emission in the multi-TeV energy range ( > 1 TeV) using passive high-galactic latitude GMCs. Assuming that the main gamma-ray production mechanism is due to proton-proton interaction, the upper limits are consistent with a cosmic-ray flux and energy density similar to that measured at Earth.

The INTEGRAL hard X-ray surveys have proven to be of fundamental importance. INTEGRAL has mapped the Galactic plane with its large field of view and excellent sensitivity. Such hard X-ray snapshots of the whole Milky Way on a time scale of a year are beyond the capabilities of past and current narrow-FOV grazing incidence X-ray telescopes. By expanding the INTEGRAL X-ray survey into shorter timescales, a productive search for transient X-ray emitters was made possible. In more than fifteen years of operation, the INTEGRAL observatory has given us a sharper view of the hard X-ray sky, and provided the triggers for many follow-up campaigns from radio frequencies to gamma-rays. In addition to conducting a census of hard X-ray sources across the entire sky, INTEGRAL has carried out, through Earth occultation maneuvers, unique observations of the large-scale cosmic X-ray background, which will without question be included in the annals of X-ray astronomy as one of the mission's most salient contribution to our understanding of the hard X-ray sky.

Central compact objects are young neutron stars emitting thermal X-rays with bolometric luminosities L X in the range 10 32 - 10 34 erg/s. Gourgouliatos, Hollerbach and Igoshev recently suggested that peculiar emission properties of central compact objects can be explained by tangled magnetic field configurations formed in a stochastic dynamo during the proto-neutron star stage. In this case the magnetic field consists of multiple small-scale components with negligible contribution of global dipolar field. We study numerically three-dimensional magneto-thermal evolution of tangled crustal magnetic fields in neutron stars. We find that all configurations produce complicated surface thermal patterns which consist of multiple small hot regions located at significant separations from each other. The configurations with initial magnetic energy of 2.5??0? 10 47 erg have temperatures of hot regions that reach ??.2 keV, to be compared with the bulk temperature of ??.1 keV in our simulations with no cooling. A factor of two in temperature is also seen in observations of central compact objects. The hot spots produce periodic modulations in light curve with typical amplitudes of ????1 %. Therefore, the tangled magnetic field configuration can explain thermal emission properties of some central compact objects.

We present a hyperspectral cube of the Crab Nebula obtained with the imaging Fourier transform spectrometer SITELLE on the Canada-France-Hawaii telescope. We describe our techniques used to deconvolve the 310 000 individual spectra (R = 9 600) containing Halpha, [NII]6548,6583, and [SII]6716,6731 emission lines and create a detailed three-dimensional reconstruction of the supernova remnant assuming uniform global expansion. We find that the general boundaries of the 3D volume occupied by the Crab are not strictly ellipsoidal as commonly assumed, and instead appear to follow a "heart-shaped" distribution that is symmetrical about the plane of the pulsar wind torus. Conspicuous restrictions in the bulk distribution of gas consistent with constrained expansion coincide with positions of the dark bays and east-west band of He-rich filaments, which may be associated with interaction with a pre-existing circumstellar disk. The distribution of filaments follows an intricate honeycomb-like arrangement with straight and rounded boundaries at large and small scales that are anti-correlated with distance from the center of expansion. The distribution is not unlike the large-scale rings observed in supernova remnants 3C 58 and Cassiopeia A, where it has been attributed to turbulent mixing processes that encouraged outwardly expanding plumes of radioactive 56Ni-rich ejecta. These characteristics reflect critical details of the original supernova of 1054 CE and its progenitor star, and may favour a low-energy explosion of an iron-core progenitor. We demonstrate that our main findings are robust despite regions of non-homologous expansion driven by acceleration of material by the pulsar wind nebula.

Neutron stars harbour extremely strong magnetic fields within their solid outer crust. The topology of this field strongly influences the surface temperature distribution, and hence the star's observational properties. In this work, we present the first realistic simulations of the coupled crustal magneto-thermal evolution of isolated neutron stars in three dimensions with account for neutrino emission, obtained with the pseudo-spectral code Parody. We investigate both the secular evolution, especially in connection with the onset of instabilities during the Hall phase, and the short-term evolution following episodes of localised energy injection. Simulations show that a resistive tearing instability develops in about a Hall time if the initial toroidal field exceeds ~ 10 15 G. This leads to crustal failures because of the huge magnetic stresses coupled with the local temperature enhancement produced by dissipation. Localised heat deposition in the crust results in the appearance of hot spots on the star surface which can exhibit a variety of patterns. Since the transport properties are strongly influenced by the magnetic field, the hot regions tend to drift away and get deformed following the magnetic field lines while cooling. The shapes obtained with our simulations are reminiscent of those recently derived from NICER X-ray observations of the millisecond pulsar PSR J0030+0451.

We present photometric and spectroscopic observations of the nearby Type Ia SN 2019yvq, from its discovery ??1 day after explosion to ??100 days after its peak brightness. This SN exhibits several unusual features, most notably an extremely bright UV excess seen within ??5 days of its explosion. As seen in Swift UV data, this early excess outshines its "peak" brightness, making this object more extreme than other SNe with early UV/blue excesses (e.g. iPTF14atg and SN 2017cbv). In addition, it was underluminous ( M B =??8.4 ), relatively quickly declining ( ? m 15 (B)=1.35 ), and shows red colors past its early blue bump. Unusual (although not unprecedented) spectral features include extremely broad-lined and high-velocity Si absorption. Despite obvious differences in peak spectra, we classify SN 2019yvq as a transitional member of the 02es-like subclass due to its similarities in several respects (e.g. color, peak luminosity, peak Ti, nebular [Ca II]). We model this dataset with a variety of published models, including SN ejecta - companion shock interaction and sub-Chandrasekhar mass WD double detonation models. Radio constraints from the VLA place an upper limit of (4.5??0)? 10 ?? M ??/yr on the mass-loss rate from a symbiotic progenitor, which does not exclude a red giant or main sequence companion. Ultimately we find that no one model can accurately replicate all aspects of the dataset, and further we find that the ubiquity of early excesses in 02es-like SNe Ia requires a progenitor system that is capable of producing isotropic UV flux, ruling out some models for this class of objects.

We investigated the prompt light curves (LCs) of long gamma-ray bursts (GRBs) from the Swift/BAT catalog. We aimed to characterize their power spectral densities (PSDs), search for quasiperiodic oscillations (QPOs), and conduct novel analyses directly in the time domain. We analyzed the PSDs using Lomb-Scargle periodograms, and searched for QPOs using wavelet scalograms. We also attempted to classify the GRBs using the Hurst exponent, H , and the A?�T plane. The PSDs fall into three categories: power law (PL; P(f)??/ f β ) with index β??0,2) , PL with a non-negligible Poisson noise level (PLC) with β??1,3) , and a smoothly broken PL (SBPL; with Poisson noise level) yielding high-frequency index β 2 ??2,6) . The latter yields break time scales on the order of 1--100\,seconds. The PL and PLC models are broadly consistent with a fully developed turbulence, β=5/3 . For an overwhelming majority of GRBs (93\%), H>0.5 , implying ubiquity of the long-term memory. We find no convincing substructure in the A?�T plane. Finally, we report on 34 new QPOs: with one or more constant leading periods, as well as several chirping signals. The presence of breaks and QPOs suggests the existence of characteristic time scales that in at least some GRBs might be related to the dynamical properties of plasma trajectories in the accretion disks powering the relativistic jets.