Astrophysics

With the black hole mass function (BHMF; assuming an exponential cutoff at a mass of ∼40 M ⊙ ) of coalescing binary black hole systems constructed with the events detected in the O1 run of the advanced LIGO/Virgo network, Liang et al.(2017) predicted that the birth of the lightest intermediate mass black holes (LIMBHs; with a final mass of ≳100 M ⊙ ) is very likely to be caught by the advanced LIGO/Virgo detectors in their O3 run. The O1 and O2 observation run data, however, strongly favor a cutoff of the BHMF much sharper than the exponential one. In this work we show that a power-law function followed by a sudden drop at ∼40 M ⊙ by a factor of ∼ a few tens and then a new power-law component extending to ≥100 M ⊙ are consistent with the O1 and O2 observation run data. With this new BHMF, quite a few LIMBH events can be detected in the O3 observation run of advanced LIGO/Virgo. The first LIMBH born in GW190521, an event detected in the early stage of the O3 run of advanced LIGO/Virgo network, provides additional motivation for our hypothesis.

The recent gravitational-wave transient GW190521 has been interpreted by the LIGO-Virgo collaboration (LVC) as sourced by a binary black hole (BH) merger. According to the LVC parameter estimation, at least one of these progenitors falls into the so-called pair-instability supernova mass gap. This raises the important question of how and when these progenitors \textit{formed}. In this paper, we use a relativistic accretion model obtained from General Relativity hydrodynamics simulations to analyze the scenario wherein the GW190521 original progenitors (OPs) formed at lower masses (and spins) and grew to their estimated LVC parameters by relativistic accretion. We consider that the environment wherein the binary is immersed has density gradients as well as a dependence on the Mach number of the gas. Taking the LVC parameter estimation at z=0.82 as the endpoint of the accretion evolution, we estimate the initial masses and spins of the OPs at three different red shifts z=100, 50 , and 20 . We found three distinct possible types of OPs: (i) 10 ?? M ???? M ??almost non-rotating (with Kerr spin parameter a ??< 10 ?? ) primordial BHs; (ii) 3 M ????0 M ??slowly rotating ( 10 ?? < a ??<0.5 ) stellar-mass BHs; (iii) 40 M ????0 M ??BHs with a moderate spin parameter a ????.5 , which could originate from the collapse of high mass Pop III stars. The mass spread is due to varying the density gradient and the relativistic Mach number of the cosmic plasma; the variation of the masses due to the origin at different red shifts, on the other hand, is negligible, ??% . For high Mach number scenarios, the BHs have low mass and spin accretion rates, leading to OPs with masses and spins close to the GW190521 LVC estimated values. . .

We present the Galactic Radio Explorer (GReX), an all-sky monitor to probe the brightest bursts in the radio sky. Building on the success of STARE2, we will search for fast radio bursts (FRBs) emitted from Galactic magnetars as well as bursts from nearby galaxies. GReX will search down to ~ten microseconds time resolution, allowing us to find new super giant radio pulses from Milky Way pulsars and study their broadband emission. The proposed instrument will employ ultra-wide band (0.7--2 GHz) feeds coupled to a high performance (receiver temperature <10 K) low noise amplifier (LNA) originally developed for the DSA-110 and DSA-2000 projects. In GReX Phase I (GReX-I), unit systems will be deployed at Owens Valley Radio Observatory (OVRO), NASA's Goldstone station, and at Telescope Array, Delta Utah. Phase II will expand the array, placing feeds in India, Australia, and elsewhere in order to build up to continuous coverage of nearly 4 ? steradians and to increase our exposure to the Galactic plane. We model the local magnetar population to forecast for GReX, finding the improved sensitivity and increased exposure to the Galactic plane could lead to dozens of FRB-like bursts per year.

The presence of relativistic electrons within the diffuse gas phase of galaxy clusters is now well established, but their detailed origin remains unclear. Cosmic ray protons are also expected to accumulate during the formation of clusters and would lead to gamma-ray emission through hadronic interactions within the thermal gas. Recently, the detection of gamma-ray emission has been reported toward the Coma cluster with Fermi-LAT. Assuming that this gamma-ray emission arises from hadronic interactions in the ICM, we aim at exploring the implication of this signal on the cosmic ray content of the Coma cluster. We use the MINOT software to build a physical model of the cluster and apply it to the Fermi-LAT data. We also consider contamination from compact sources and the impact of various systematic effects. We confirm that a significant gamma-ray signal is observed within the characteristic radius θ 500 of the Coma cluster, with a test statistic TS~27 for our baseline model. The presence of a possible point source may account for most of the observed signal. However, this source could also correspond to the peak of the diffuse emission of the cluster itself and extended models match the data better. We constrain the cosmic ray to thermal energy ratio within R 500 to X CRp = 1.79 +1.11 ??.30 \% and the slope of the energy spectrum of cosmic rays to α= 2.80 +0.67 ??.13 . Finally, we compute the synchrotron emission associated with the secondary electrons produced in hadronic interactions assuming steady state. This emission is about four times lower than the overall observed radio signal, so that primary cosmic ray electrons or reacceleration of secondary electrons is necessary to explain the total emission. Assuming an hadronic origin of the signal, our results provide the first quantitative measurement of the cosmic ray proton content in a cluster.[Abridged]

We analyze about 12 years of Fermi-LAT data in the direction of the Andromeda galaxy (M31). We robustly characterize its spectral and morphological properties against systematic uncertainties related to the modeling of the Galactic diffuse emission. We perform this work by adapting and exploiting the potential of the skyFACT adaptive template fitting algorithm. We reconstruct the gamma-ray image of M31 in a template-independent way, and we show that flat spatial models are preferred by data, indicating an extension of the γ -ray emission of about 0.3-0.4 degree for the bulge of M31. This study also suggests that a second component, extending to at least 1 degree, contributes to the observed total emission. We quantify systematic uncertainties related to mis-modeling of Galactic foreground emission at the level of 2.9%.

Using data from the \emph{Fermi} Large Area Telescope we perform a detailed study of the GeV emission in the direction of G51.26+0.09 to constrain its origin, its possible relation to this SNR, the star-forming region G051.010+00.060 seen nearby in the sky, or the pulsars known in the region, and to derive the properties of the underlying cosmic ray particles producing the non-thermal radiation. We also study properties of the environment which could shed light on the nature of the source of the gamma rays. We compare the morphology of the gamma-ray radiation to that of the emission detected at radio wavelengths in previous observations. Modeling the measured spectrum and fluxes of the high-energy radiation allows us to derive the properties of the particle populations that could produce this emission in several possible scenarios. We use existing data from 13 CO emission and neutral hydrogen emission to probe the environment, searching for possible morphological features associated to the gamma rays and SNR. We rule out the star-forming region G051.010+00.060 as the origin of the GeV emission. The correspondence seen between the gamma-ray and radio morphologies support an SNR scenario, where the object responsible is more extended than G51.26+0.09, or is made up of more than one unresolved SNR. Given the flat spectral energy distribution observed at GeV energies and the radio flux upper limits, we also rule out bremsstrahlung emission as the origin of the gamma rays. A pulsar wind nebula origin of the high-energy photons, associated to the pulsar PSR J1926+1613, cannot be ruled out or confirmed, due to its unknown parameters such as spin-down power and age, while the pulsars PSR J1924+1628 and PSR J1924+1631 are too far away to be the source of gamma rays.

Convection and turbulence in core-collapse supernovae (CCSNe) are inherently three-dimensional in nature. However, 3D simulations of CCSNe are computationally demanding. Thus, it is valuable to modify simulations in spherical symmetry to incorporate 3D effects using some parametric model. In this paper, we report on the formulation and implementation of general relativistic (GR) neutrino-driven turbulent convection in the spherically symmetric core-collapse supernova code \texttt{GR1D}. This is based upon the recently proposed method of Supernova Turbulence in Reduced-dimensionality (\textit{STIR}) in Newtonian simulations from \cite{Couch2020_STIR}. When the parameters of this model are calibrated to 3D simulations, we find that our GR formulation of \textit{STIR} requires larger turbulent eddies to achieve a shock evolution similar to the original \textit{STIR} model. We also find that general relativity may alter the correspondence between progenitor mass and successful vs.~failed explosions.

In 2019, the Event Horizon Telescope Collaboration (EHTC) has published the first image of a supermassive black hole (SMBH) obtained via the Very Large Baseline Interferometry (VLBI) technique. In the future, it is expected that additional and more sensitive VLBI observations will be pursued for other nearby Active Galactic Nuclei (AGN), and it is therefore important to understand which possible features can be expected in such images. In this paper, we post-process General Relativistic Magneto-Hydrodynamical (GR-MHD) simulations which include resistivity, thus providing a self-consistent jet formation model (with resistive mass loading) launched from a thin disc. The ray-tracing is done using the General Relativistic Ray-Tracing code GRTRANS assuming synchrotron emission. We study the appearance of the black hole environment including the accretion disc, winds and jets under a large range of condition, varying black hole mass, accretion rate, spin, inclination angle, disc parameters and observed frequency. When we adopt M87-like parameters, we show that we can reproduce a ring-like feature (similar as observed by the EHT) for some of our simulations. The latter suggests that such thin disc models are thus likely consistent with the observed results. Depending on their masses, accretion rates, spin and other parameters, we note that other SMBHs may show additional features like winds and jets depending on the sensitivity that can be reached in the observations.

GrailQuest (Gamma-ray Astronomy International Laboratory for Quantum Exploration of Space-Time) is an ambitious astrophysical mission concept that uses a fleet of small satellites whose main objective is to search for a dispersion law for light propagation in vacuo. Within Quantum Gravity theories, different models for space-time quantization predict relative discrepancies of the speed of photons w.r.t. the speed of light that depend on the ratio of the photon energy to the Planck energy. This ratio is as small as 1E-23 for photons in the gamma-ray band (100 keV). Therefore, to detect this effect, light must propagate over enormous distances and the experiment must have extraordinary sensitivity. Gamma-Ray Bursts, occurring at cosmological distances, could be used to detect this tiny signature of space-time granularity. This can be obtained by coherently combine a huge number of small instruments distributed in space to act as a single detector of unprecedented effective area. This is the first example of high-energy distributed astronomy: a new concept of modular observatory of huge overall collecting area consisting in a fleet of small satellites in low orbits, with sub-microsecond time resolution and wide energy band (keV-MeV). The enormous number of collected photons will allow to effectively search these energy dependent delays. Moreover, GrailQuest will allow to perform temporal triangulation of impulsive events with arc-second positional accuracies: an extraordinary sensitive X-ray/Gamma all-sky monitor crucial for hunting the elusive electromagnetic counterparts of Gravitational Waves, that will play a paramount role in the future of Multi-messenger Astronomy. A pathfinder of GrailQuest is already under development through the HERMES (High Energy Rapid Modular Ensemble of Satellites) project: a fleet of six 3U cube-sats to be launched by the end of 2022.

After the prediction of many sub- and super-Chandrasekhar (at least a dozen for the latter) limiting mass white dwarfs, hence apparently peculiar class of white dwarfs, from the observations of luminosity of type Ia supernovae, researchers have proposed various models to explain these two classes of white dwarfs separately. We earlier showed that these two peculiar classes of white dwarfs, along with the regular white dwarfs, can be explained by a single form of the f(R) gravity, whose effect is significant only in the high-density regime, and it almost vanishes in the low-density regime. However, since there is no direct detection of such white dwarfs, it is difficult to single out one specific theory from the zoo of modified theories of gravity. We discuss the possibility of direct detection of such white dwarfs in gravitational wave astronomy. It is well-known that in f(R) gravity, more than two polarization modes are present. We estimate the amplitudes of all the relevant modes for the peculiar as well as the regular white dwarfs. We further discuss the possibility of their detections through future-based gravitational wave detectors, such as LISA, ALIA, DECIGO, BBO, or Einstein Telescope, and thereby put constraints or rule out various modified theories of gravity. This exploration links the theory with possible observations through gravitational wave in f(R) gravity.