Astrophysics

The Cadmium Zinc Telluride Imager (CZTI) is an imaging instrument onboard AstroSat. This instrument operates as a nearly open all-sky detector above ~60 keV, making possible long integrations irrespective of the spacecraft pointing. We present a technique based on the AstroSat-CZTI data to explore the hard X-ray characteristics of the γ -ray pulsar population. We report highly significant ( ??0? ) detection of hard X-ray (60--380 keV) pulse profile of the Crab pulsar using ??5000 ks of CZTI observations within 5 to 70 degrees of Crab position in the sky, using a custom algorithm developed by us. Using Crab as our test source, we estimate the off-axis sensitivity of the instrument and establish AstroSat-CZTI as a prospective tool in investigating hard X-ray characteristics of γ -ray pulsars as faint as 10 mCrab.

In this work, at first we present a model of studying astrophysical flows of binary systems and microquasars based on the laws of relativistic magnetohydrodynamics. Then, by solving the time independent transfer equation, we estimate the primary and secondary particle distributions within the hadronic astrophysical jets as well as the emissivities of high energy neutrinos and γ -rays. One of our main goals is, by taking into consideration the various energy-losses of particles into the hadronic jets, to determine through the transport equation the respective particle distributions focusing on relativistic hadronic jets of binary systems. As a concrete example we examine the extragalactic binary system LMC X-1 located in the Large Magellanic Cloud, a satellite galaxy of our Milky Way Galaxy.

Binary millisecond pulsars (MSPs) are believed to have descended from low-mass X-ray binaries (LMXBs), which have experienced substantial mass transfer and tidal circularization. Therefore, they should have very circular orbits. However, the discovery of several eccentric binary MSPs (with eccentricity e??.01??.1 ) challenges this standard picture. Three models have been proposed thus far based on accretion-induced collapse of massive white dwarfs (WDs), neutron star-strange star transition, and formation of circumbinary disks. All of them are subject to various uncertainties, and are not entirely consistent with observations. Here we propose an alternative model taking into account the influence of thermonuclear flashes on proto-WDs. We assume that the flashes lead to asymmetrical mass ejection, which imparts a mild kick on the proto-WDs. By simulating orbital changes of binary MSPs with multiple shell flashes, we show that it is possible to reproduce the observed eccentricities, provided that the kick velocities are around a few kms ?? .

We introduce a family of equations of state (EoS) for hybrid neutron star (NS) matter that is obtained by a two-zone parabolic interpolation between a soft hadronic EoS at low densities and a set of stiff quark matter EoS at high densities within a finite region of chemical potentials μ H <μ< μ Q . Fixing the hadronic EoS as the APR one and chosing the color-superconductiong, nonlocal NJL model with two free parameters for the quark phase, we perform Bayesian analyses with this two-parameter family of hybrid EoS. Using three different sets of observational constraints that include the mass of PSR J0740+6620, the tidal deformability for GW170817 and the mass-radius relation for PSR J0030+0451 from NICER as obligatory (set 1), while set 2 uses the possible upper limit on the maximum mass from GW170817 as additional constraint and set 3 instead the possibility that the lighter object in the asymmetric binary merger GW190814 is a neutron star. We confirm that in any case the quark matter phase has to be color superconducting with the dimensionless diquark coupling approximately fulfilling the Fierz relation η D =0.75 and the most probable solutions exhibiting a proportionality between η D and η V , the coupling of the repulsive vector interaction that is required for a sufficiently large maximum mass. We anticipate the outcome of the NICER radius measurement on PSR J0740+6220 as a fictitious constraint and find evidence for claiming that GW190814 was a binary black hole merger if the radius will be 11 km or less.

The increase in the sensitivity of gravitational wave interferometers will bring additional detections of binary black hole and double neutron star mergers. It will also very likely add many merger events of black hole - neutron star binaries. Distinguishing mixed binaries from binary black holes mergers for high mass ratios could be challenging because in this situation the neutron star coalesces with the black hole without experiencing significant disruption. To investigate the transition of mixed binary mergers into those behaving more like binary black hole coalescences, we present results from merger simulations for different mass ratios. We show how the degree of deformation and disruption of the neutron star impacts the inspiral and merger dynamics, the properties of the final black hole, the accretion disk formed from the circularization of the tidal debris, the gravitational waves, and the strain spectrum and mismatches. The results also show the effectiveness of the initial data method that generalizes the Bowen-York initial data for black hole punctures to the case of binaries with neutron star companions.

In this pilot study, we investigate the use of a deep learning (DL) model to temporally evolve the dynamics of gas accreting onto a black hole in the form of a radiatively inefficient accretion flow (RIAF). We have trained a machine to forecast such a spatiotemporally chaotic system -- i.e. black hole weather forecasting -- using a convolutional neural network (CNN) and a training dataset which consists of numerical solutions of the hydrodynamical equations, for a range of initial conditions. We find that deep neural networks seem to learn well black hole accretion physics and evolve the accretion flow orders of magnitude faster than traditional numerical solvers, while maintaining a reasonable accuracy for a long time. For instance, CNNs predict well the temporal evolution of a RIAF over a long duration of 8? 10 4 GM/ c 3 , which corresponds to 80 dynamical times at r=100GM/ c 2 . The DL model is able to evolve flows from initial conditions not present in the training dataset with good accuracy. Our approach thus seems to generalize well. Once trained, the DL model evolves a turbulent RIAF on a single GPU four orders of magnitude faster than usual fluid dynamics integrators running in parallel on 200 CPU cores. We speculate that a data-driven machine learning approach should be very promising for accelerating not only fluid dynamics simulations, but also general relativistic magnetohydrodynamic ones.

In this letter, we note that the observed in the LIGO / Virgo experiment ratio of the detection rate of black holes to the rate of detection of binary neutron stars requires the assumption of a "conservative" collapse of massive stars into a black hole: almost all the mass of the collapsing star goes under the horizon. This is consistent with the large masses of black holes detected by LIGO/Virgo. On the other hand, the assumption of a small loss of matter during the collapse into a black hole is in good agreement with the small eccentricity of Single-lined Binaries. At the same time, the absence of X-rays from most black holes in binary systems with blue stars is explained. We argue that three sets of LIGO / Virgo observations and data on the Single-lined Binary with a Candidate Black Hole Component confirm the scenario of the evolution of massive field binaries.

We present a novel technique to study Type Ia supernovae by constraining surviving companions of historical extragalactic SN by combining archival photographic plates and Hubble Space Telescope imaging. We demonstrate this technique for Supernova 1972E, the nearest known SN Ia in over 125 years. Some models of Type Ia supernovae describe a white dwarf with a non-degenerate companion that donates enough mass to trigger thermonuclear detonation. Hydrodynamic simulations and stellar evolution models show that these donor stars survive the explosion, and show increased luminosity for at least a thousand years. Thus, late-time observations of the exact location of a supernova after its ejecta have faded can constrain the presence of a surviving donor star and progenitor models. We find the explosion site of SN 1972E by analyzing 17 digitized photographic plates taken with the European Southern Observatory 1m Schmidt and 1 plate taken with the Cerro Tololo Inter-American Observatory 1.5m telescope from 1972-1974. Using the \textit{Gaia} eDR3 catalog to determine Supernova 1972E's equatorial coordinates yields: α = 13 h 39 m 52.708 s ± 0.004 s and δ = ??31\degree 40' 8\farcs97 ± 0\farcs04 (ICRS). In 2005, HST/ACS imaged NGC 5253, the host galaxy of SN 1972E, with the F435W , F555W , and F814W filters covering the explosion site. The nearest source detected is offset by 3.0 times our positional precision, and is inconsistent with the colors expected of a surviving donor star. Thus, the 2005 HST observation rules out all Helium-star companion models, and the most luminous main-sequence companion model currently in the literature. The remaining main-sequence companion models could be tested with deeper HST imaging.

The most extreme active galactic nuclei (AGN) are the radio active ones whose relativistic jet propagates close to our line of sight. These objects were first classified according to their emission line features into flat-spectrum radio quasars (FSRQs) and BL Lacertae objects (BL Lacs). More recently, observations revealed a trend between these objects known as the \emph{blazar sequence}, along with an anti-correlation between the observed power and the frequency of the synchrotron peak. In the present work, we propose a fairly simple idea that could account for the whole blazar population: all jets are launched with similar energy per baryon, independently of their power. In the case of FSRQs, the most powerful jets, manage to accelerate to high bulk Lorentz factors, as observed in the radio. As a result, they have a rather modest magnetization in the emission region, resulting in magnetic reconnection injecting a steep particle energy distribution and, consequently, steep emission spectra in the γ -rays. For the weaker jets, namely BL Lacs, the opposite holds true; i.e., the jet does not achieve a very high bulk Lorentz factor, leading to more magnetic energy available for non-thermal particle acceleration, and harder emission spectra at frequencies ≳ GeV. In this scenario, we recover all observable properties of blazars with our simulations, including the \emph{blazar sequence} for models with mild baryon loading ( 50≲μ≲80 ). This interpretation of the blazar population, therefore, tightly constrains the energy per baryon of blazar jets regardless of their accretion rate.

The accreting millisecond X-ray pulsar Swift J1756.9 ??2508 went into outburst in April 2018 and June 2019, 8.7 yr after the previous activity period. We investigated the temporal, timing and spectral properties of these two outbursts using data from NICER, XMM-Newton, NuSTAR, INTEGRAL, Swift and Insight-HXMT. The two outbursts exhibited similar broad-band spectra and X-ray pulse profiles. For the first time, we report the detection of the pulsed emission up to ??00 keV observed by Insight-HXMT during the 2018 outburst. We also found the pulsation up to ??0 keV observed by NICER and NuSTAR during the 2019 outburst. We performed a coherent timing analysis combining the data from two outbursts. The binary system is well described by a constant orbital period over a time span of ??2 years. The time-averaged broad-band spectra are well fitted by an absorbed thermal Comptonization model in a slab geometry with the electron temperature k T e =40 -50 keV, Thomson optical depth ???.3 , blackbody seed photon temperature k T bb,seed ??0.7-0.8 keV and hydrogen column density of N H ??.2? 10 22 cm ?? . We searched the available data for type-I (thermonuclear) X-ray bursts, but found none, which is unsurprising given the estimated low peak accretion rate ( ??.05 of the Eddington rate) and generally low expected burst rates for hydrogen-poor fuel. Based on the history of four outbursts to date, we estimate the long-term average accretion rate at roughly 5? 10 ??2 M ??yr ?? for an assumed distance of 8 kpc. The expected mass transfer rate driven by gravitational radiation in the binary implies the source can be no closer than 4 kpc.