Search

Secure confirmation that a gravitational wave (GW) has been gravitationally lensed would bring together these two pillars of General Relativity for the first time. This breakthrough is challenging for many reasons, including: GW sky localization uncertainties dwarf the angular scale of gravitational lensing, the mass and structure of gravitational lenses is diverse, the mass function of stellar remnant compact objects is not yet well constrained, and GW detectors do not operate continuously. We introduce a new approach that is agnostic to the mass and structure of the lenses, compare the efficiency of different methods for lensed GW discovery, and explore detection of lensed kilonova counterparts as a direct method for localising candidates. Our main conclusions are: (1) lensed neutron star mergers (NS-NS) are magnified into the "mass gap" between NS and black holes, therefore selecting candidates from public GW alerts with high mass gap probability is efficient, (2) the rate of detectable lensed NS-NS will approach one per year in the mid-2020s, (3) the arrival time difference between lensed NS-NS images is $1\,\rm sec\lesssim\Delta t\lesssim1\,year$, and thus well-matched to the operations of GW detectors and optical telescopes, (4) lensed kilonova counterparts are faint at peak (e.g.\ $r_{\rm AB}\simeq24-26$ in the mid-2020s), fade quickly ($d, Comment: Accepted by MNRAS. 20 pages, 10 figures

Mergers of black hole-neutron star (BHNS) binaries have now been observed by GW detectors with the recent announcement of GW200105 and GW200115. Such observations not only provide confirmation that these systems exist, but will also give unique insights into the death of massive stars, the evolution of binary systems and their possible association with gamma-ray bursts, $r$-process enrichment and kilonovae. Here we perform binary population synthesis of isolated BHNS systems in order to present their merger rate and characteristics for ground-based GW observatories. We present the results for 420 different model permutations that explore key uncertainties in our assumptions about massive binary star evolution (e.g. mass transfer, common-envelope evolution, supernovae), and the metallicity-specific star formation rate density, and characterize their relative impacts on our predictions. We find intrinsic local BHNS merger rates spanning $\mathcal{R}_{\rm{m}}^0 \approx 4$-$830\,\rm{Gpc}^{-3}\,\rm{yr}^{-1}$ for our full range of assumptions. This encompasses the rate inferred from recent BHNS GW detections, and would yield detection rates of $\mathcal{R}_{\rm{det}} \approx 1$-$180\, \rm{yr}^{-1}$ for a GW network consisting of LIGO, Virgo and KAGRA at design sensitivity. We find that the binary evolution and metallicity-specific star formation rate density each impact the predicted merger rates by order $\mathcal{O}(10)$. We also present predictions for the GW detected BHNS merger properties and find that all 420 model variations predict that $\lesssim 5\%$ of the BHNS mergers have BH masses $\gtrsim 18\,M_{\odot}$, total masses $ \gtrsim 20\,M_{\odot}$, chirp masses $\gtrsim 5.5\,M_{\odot}$, mass ratios $ \gtrsim 12$ or $\lesssim 2$. Moreover, we find that massive NSs $\gtrsim 2\,M_{\odot}$ are expected to be commonly detected in BHNS mergers in almost all our model variations., 38 pages, 18 figures, accepted to MNRAS. The authors welcome suggestions and feedback. All data and code to reproduce the results in this paper are publicly available

Direct collapse black holes~(BH) are promising candidates for producing massive $z\gtrsim 6$ quasars, but their formation requires fine-tuned conditions. In this work, we use cosmological zoom simulations to study systematically the impact of requiring: 1) low gas angular momentum, and 2) a minimum incident Lyman-Werner~(LW) flux in order to form BH seeds. We probe the formation of seeds (with initial masses of $M_{\rm seed} \sim 10^4$ - $10^6 M_{\odot}/h)$ in halos with a total mass $> 3000\times M_{\mathrm{seed}}$ and a dense, metal poor gas mass $> 5\times M_{\mathrm{seed}}$. We find that the seed-forming halos have a prior history of star formation and metal enrichment, but contain pockets of dense, metal poor gas. When seeding is further restricted to halos with low gas spins, the number of seeds formed is suppressed by factors of $\sim6$ compared to the baseline model, regardless of the seed mass. Seed formation is much more strongly impacted if the dense, metal poor gas is required to have a critical LW flux ($J_{\mathrm{crit}}$). Even for $J_{\mathrm{crit}}$ values as low as $50J_{21}$, no $8\times10^{5}M_{\odot}/h$ seeds are formed. While lower mass ($1.25\times10^{4},1\times10^{5} M_{\odot}/h$) seeds do form, they are strongly suppressed~(by factors of $\sim10-100$) compared to the baseline model at gas mass resolutions of $\sim10^4~M_{\odot}/h$ (with even stronger suppression at higher resolutions). As a result, BH merger rates are also similarly suppressed. Since early BH growth is dominated by mergers in our models, no seeds are able to grow to the supermassive regime~($\gtrsim10^6 M_{\odot}/h$) by $z=7$. Our results hint that producing the bulk of the $z\gtrsim6$ supermassive BH population may require alternate seeding scenarios that do not depend on the LW flux, early BH growth dominated by rapid or super-Eddington accretion, or a combination of these possibilities., 21 pages, 13 figures

Using bona-fide black hole (BH) mass estimates from reverberation mapping and the line ratio [SiVI] 1.963$\mu$m/Br$\gamma_{\rm broad}$ as tracer of the AGN ionising continuum, a novel BH-mass scaling relation of the form log($M_{\rm BH}) = (6.40\pm 0.17) - (1.99\pm 0.37) \times$ log ([SiVI]/Br$\gamma_{\rm broad})$, dispersion 0.47 dex, over the BH mass interval, $10^6 - 10^8$ M$_{\odot}$ is found. Following on the geometrically thin accretion disc approximation and after surveying a basic parameter space for coronal lines production, we believe one of main drivers of the relation is the effective temperature of the disc, which is effectively sampled by the [SiVI] 1.963$\mu$m coronal line for the range of BH masses considered. By means of CLOUDY photoionisation models, the observed anti-correlation appears to be formally in line with the thin disc prediction T_disc $\propto {M_{\rm BH}}^{-1/4}$., Comment: 22 pages, 16 figures. Accepted for publication in MNRAS. arXiv admin note: text overlap with arXiv:2010.00075

Supermassive black hole masses (MBH) can dynamically be estimated with various methods and using different kinematic tracers. Different methods have only been cross-checked for a small number of galaxies and often show discrepancies. To understand these discrepancies, detailed cross-comparisons of additional galaxies are needed. We present the first part of our cross-comparison between stellar- and gas-based MBH estimates in the nearby fast-rotating early-type galaxy NGC 6958. The measurements presented here are based on ground-layer adaptive optics-assisted Multi-Unit Spectroscopic Explorer (MUSE) science verification data at around 0.6 arcsec spatial resolution. The spatial resolution is a key ingredient for the measurement and we provide a Gaussian parametrisation of the adaptive optics-assisted point spread function (PSF) for various wavelengths. From the MUSE data, we extracted the stellar kinematics and constructed dynamical models. Using an axisymmetric Schwarzschild technique, we measured an MBH of (3.6+2.7-2.4)\times 10^8 Msun at 3\sigma significance taking kinematical and dynamical systematics (e.g.,radially-varying mass-to-light ratio) into account. We also added a dark halo, but our data does not allow to constrain the dark matter fraction. Adding dark matter with an abundance matching prior results in a 25 per cent more massive black hole. Jeans anisotropic models return MBH of (4.6+2.5-2.7) \times 10^8 Msun and (8.6+0.8-0.8) \times 10^8 Msun at 3\sigma confidence for spherical and cylindrical alignment of the velocity ellipsoid, respectively. In a follow-up study, we will compare the stellar-based MBH with those from cold and warm gas tracers, which will provide additional constraints for the MBH for NGC 6958, and insights into assumptions that lead to potential systematic uncertainty., Comment: Accepted for publication in MNRAS, 21 pages, 18 Figures (Including Appendix). Supplementary information is added at the end of the file (3 pages)

Observations found that star clusters contain a large fraction of binaries. Tight binaries are an important heating source that influences the long-term dynamical evolution of star clusters. However, due to the limitation of $N$-body tool, previous theoretical modelling for globular clusters (GCs) by using direct $N$-body simulations have not investigated how a large fraction of primordial binaries affect their long-term evolution. In this work, by using the high-performance $N$-body code, PeTar, we carry out star-by-star models for intermediate massive GCs ($N=100000$) with the primordial binary fraction varying from 0 to 1. We find that when a stellar-mass black hole (BH) subsystem exists, the structural evolution of GCs (core and half-mass radii) only depends on the properties of massive primordial binaries, because they affect the number of BH binaries (BBHs), which dominate the binary heating process. Low-mass binaries including double white dwarf binaries (BWDs) have almost no influence on the dynamics. Meanwhile, only gravitational wave (GW) mergers from BBHs are strongly affected by dynamical interactions, while low-mass mergers from BWDs show no difference in the isolated environment (field) and in GCs. Low-mass binaries become important only after most BHs escape and the core collapse of light stars occurs. Our result suggests that for $N$-body modelling of GCs with a black hole subsystem dominating binary heating, it is not necessary to include low-mass binaries. These binaries can be studied separately by using standalone binary stellar evolution codes. This way can significantly reduce the computing cost., Comment: 10 pages, 6 figures, submitted to MNRAS

We present results from the B-type Binaries Characterisation (BBC) programme, a multi-epoch spectroscopic study of 88 early B-type binary candidates in the 30 Doradus region of the Large Magellanic Cloud (LMC). From radial-velocity analysis of 29 observational epochs we confirm the binary status of 64 of our targets, comprising 50 SB1 and 14 SB2 B-type binaries. A further 20 systems (classified as SB1*) show clear signs of periodicity but with more tentative periods. Orbital solutions are presented for these 84 systems, providing the largest homogeneous sample to date of the binary properties of early B-type stars. Our derived orbital-period distribution is generally similar to those for samples of more massive (O-type) binaries in both the LMC and the Galaxy. This similarity with the properties of the more massive O-type binaries is important as early B-type stars are expected to account for the majority of core-collapse supernovae. Differences in the period distributions of the different samples start to increase above 4 d, and are also present between the earliest (B0-0.7) and later-type (B1-2.5) systems within the BBC sample, although further study is required to understand if this is an observational bias or a real physical effect. We have examined the semi-amplitude velocities and orbital periods of our sample to identify potential candidates that could hide compact companions. Comparing with probability distributions of finding black hole companions to OB-type stars from a recent theoretical study, we have found 16 binaries in the higher probability region that warrant further study., 42 pages (main body: 28 pages, supplementary material: 14 pages), 20 figures. Accepted for publication in MNRAS

Dual-Active Galactic Nuclei (AGN) are a natural consequence of the hierarchical structure formation scenario, and can provide an important test of various models for black hole growth. However, due to their rarity and difficulty to find at high redshift, very few confirmed dual-AGN are known at the epoch where galaxy formation peaks. Here we report the discovery of a gravitationally lensed dual-AGN system at redshift 2.37 comprising two optical/IR quasars separated by 6.5+/-0.6 kpc, and a third compact (R_eff = 0.45+/-0.02 kpc) red galaxy that is offset from one of the quasars by 1.7+/-0.1 kpc. From Very Large Array imaging at 3 GHz, we detect 600 and 340 pc-scale radio emission that is associated with both quasars. The 1.4 GHz luminosity densities of the radio sources are about 10^24.35 W / Hz, which is consistent with weak jets. However, the low brightness temperature of the emission is also consistent with star-formation at the level of 850 to 1150 M_sun / yr. Although this supports the scenario where the AGN and/or star-formation is being triggered through an ongoing triple-merger, a post-merger scenario where two black holes are recoiling is also possible, given that neither has a detected host galaxy., Comment: 5 pages, 3 figures, accepted for publication in MNRAS Letters

We explore, using numerical simulations, the influence of mass and distance on the evolution of triple black hole systems. Following in the direction of Burrau's famous 3,4,5 problem, black holes are initially placed at the vertices of Pythagorean triangles. Numerical integration of orbits was conducted using relativistic corrections (post-Newtonian) up to the 2.5$^{th}$ order with ARCcode. As a descriptor of the evolution of the systems, the lifetimes, the number of two-body encounters and the number of mergers were all analysed. We found that as the mass unit of the black holes increased there was strong positive correlation with the fraction of mergers (0.9868), strong negative correlation with the average number of two-body encounters (-0.9016) and the average lifetimes of the triple systems decayed exponentially (determination coefficient of 0.9986). Around the mass unit range of 10$^{5.5}$M$_{\odot}$-10$^{5.6}$M$_{\odot}$, there was a transition from escape dominated dynamics to merger dominated dynamics. However, in the mass unit range of 10$^6$ M$_{\odot}$ $-$ 10$^{9}$M$_{\odot}$ with 1 pc distance unit, we find that 25\% of cases resulted in the escape of a supermassive black hole (SMBH) which may be a cause for wandering SMBH's found in some galaxies/galactic merger remnants., Comment: 11 pages, 14 figures

3C 84 (NGC 1275) is the radio source at the centre of the Perseus cluster and exhibits a bright radio jet. We observed the source with the Global Millimeter VLBI Array (GMVA) between 2008 and 2015, with a typical angular resolution of ∼50 μas. The observations revealed a consistent double nuclear structure separated by ∼770 gravitational radii assuming a black hole mass of 3.2 × 108 M⊙. The region is likely too broad and bright to be the true jet base anchored in the accretion disc or black hole ergosphere. A cone and parabola were fit to the stacked (time averaged) image of the nuclear region. The data did not strongly prefer either fit, but combined with a jet/counter-jet ratio analysis, an upper limit on the viewing angle to the inner jet region of ≤35° was found. This provides evidence for a variation of the viewing angle along the jet (and therefore a bent jet) within ∼0.5 pc of the jet launching region. In the case of a conical jet, the apex is located ∼2400 gravitational radii upstream of the bright nuclear region and up to ∼600 gravitational radii upstream in the parabolic case. We found a possible correlation between the brightness temperature and relative position angle of the double nuclear components, which may indicate rotation within the jet. © 2021 The Author(s)., This research was supported by Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education (NRF-2021R1A6A3A01086420). JAH was supported by Korea Research Fellowship Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Science and ICT, South Korea (2018H1D3A1A02032824) and the research grant (2021R1C1C1009973). S-SL was supported by the National Research Foundation of Korea (NRF) grant funded by the Korea government (MIST) (2020R1A2C2009003). J-YK is supported for this research by the International Max-Planck Research School (IMPRS) for Astronomy and Astrophysics at the University of Bonn and Cologne. This research was supported by an appointment to the NASA Postdoctoral Program at the Goddard Space Flight Center, administered by Universities Space Research Association through a contract with NASA. RL acknowledges the support of the Spanish Ministerio de Economía y Competitividad (grant PID2019-108995GB-C21), the Consejería de Economía, Conocimiento, Empresas y Universidad, Junta de Andalucía (grant P18-FR-1769), and the Consejo Superior de Investigaciones Científicas (grant 2019AEP112). GFP is supported for this research by the International Max-Planck Research School (IMPRS) for Astronomy and Astrophysics at the University of Bonn and Cologne. ST and MK acknowledge support via NRF grant 2019R1F1A1059721. This research has made use of data obtained with the Global Millimeter VLBI Array (GMVA), which consists of telescopes operated by the MPIfR, IRAM, Onsala, Metsahovi, Yebes, the Korean VLBI Network, the Green Bank Observatory, and the Very Long Baseline Array (VLBA). The VLBA is an instrument of the NRAO, which is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc. The data were correlated at the correlator of the Max Planck Institute for Radioastronomy (MPIfR) in Bonn, Germany. This work made use of the Swinburne University of Technology software correlator, developed as part of the Australian Major National Research Facilities Programme and operated under licence., With funding from the Spanish government through the Severo Ochoa Centre of Excellence accreditation SEV-2017-0709.

The quest for high-redshift quasars has led to a series of record-breaking sources, with the current record holder at $z=7.642$. Here, we show how future detections of $z>8$ quasars impact the constraints on the parameters for black hole growth and seed models. Using broad flat priors on the growth parameters (Eddington ratio $\,f_{\rm Edd}$, duty cycle ${\cal D}$, seed mass $M_{\rm \bullet, seed}$ and radiative efficiency $\epsilon$), we show that the large uncertainties in their determination decrease by a factor $\sim 5$ when a quasar's detection redshift goes from $z=9$ to $z=12$. In this high-redshift regime, $\epsilon$ tends to the lowest value allowed, and the distribution for $M_{\rm \bullet, seed}$ peaks well inside the heavy seed domain. Remarkably, two quasars detected at $z > 7$ with low accretion rates (J1243+0100 and J0313-1806) already tighten the available parameter space, requiring $M_{\rm \bullet, seed} > 10^{3.5} \,{\rm M_\odot}$ and $\epsilon < 0.1$. The radiative efficiency is a crucial unknown, with factor $\sim 2$ changes able to modify the predicted mass by $\sim 3$ orders of magnitude already at $z\sim 9$. The competing roles of inefficient accretion (decreasing $\epsilon$) and black hole spin-up (increasing $\epsilon$) significantly impact growth models. Finally, we suggest that yields currently predicted by upcoming quasar surveys (e.g., Euclid) will be instrumental for determining the most-likely seed mass regime. For example, assuming thin-disk accretion, a detection of a quasar with $M_\bullet \sim 10^{10} \,{\rm M_\odot}$ by $z\sim 9-10$ would exclude the entire parameter space available for light seeds and dramatically reduce the one for heavy seeds., Comment: Accepted for publication in MNRAS. This is the final version of the manuscript. 8 pages, 5 figures

The orbital dynamics of fast spinning neutron stars encountering a massive black hole (BH) with unbounded orbits are investigated using the quadratic-in-spin Mathisson–Papapetrou–Dixon (MPD) formulation. We consider the motion of the spinning neutron stars with astrophysically relevant speed in the gravity field of the BH. For such slow-speed scattering, the hyperbolic orbits followed by these neutron stars all have near the e = 1 eccentricity, and have distinct properties compared with those of e ≫ 1. We have found that, compared with geodesic motion, the spin–orbit and spin–spin coupling will lead to a variation of scattering angles at spatial infinity, and this variation is more prominent for slow-speed scattering than fast-speed scattering. Such a variation leads to an observable difference in pulse-arrival-time within a few hours of observation, and up to a few days or months for larger BH masses or longer spinning periods. Such a relativistic pulsar-BH system also emits a burst of gravitational waves (GWs) in the sensitivity band of Laser Interferometric Space Antenna, and for optimal settings, can be seen up to $100\, {\rm Mpc}$ away. A radio follow up of such a GW burst with SKA or FAST will allow for measuring the orbital parameters with high accuracy and testing the predictions of general relativity.

We present a near-infrared spectroscopic survey of narrow-line Seyfert 1 galaxies in the southern hemisphere (using the SOFI instrument on the ESO-NTT telescope), sampled from optical surveys. We examine the kinematics of the broad-line region, probed by the emission line width of hydrogen (Pa$\alpha$ and H$\beta$). We observed 57 objects, of which we could firmly measure Pa$\alpha$ in 49 cases. We find that a single Lorentzian fit (preferred on theoretical grounds) is preferred over multi-component Gaussian fits to the line profiles; a lack of narrow-line region emission, overwhelmed by the pole-on view of the broad line region (BLR) light, supports this. We recompute the catalog black hole (BH) mass estimates, using the values of FWHM and luminosity of H$\beta$, both from catalog values and re-fitted Lorentzian values. We find a relationship slope greater than unity compared to the catalog values. We ascribe this to contamination by galactic light or difficulties with line flux measurements. However, the comparison of masses computed by the fitted Lorentzian and Gaussian measurements show a slope close to unity. Comparing the BH masses estimated from both Pa$\alpha$ and H$\beta$, the line widths and fluxes shows deviations from expected; in general, however, the computed BH masses are comparable. We posit a scenario where an intermixture of dusty and dust-free clouds (or alternately a structured atmosphere) differentially absorbs the line radiation of the BLR, due to dust absorption and hydrogen bound-free absorption., Comment: Accepted for publication by MNRAS. 15 pages, 15 figures

Black-hole (BH) accretion disks formed in compact-object mergers or collapsars may be major sites of the rapid-neutron-capture (r-)process, but the conditions determining the electron fraction (Y_e) remain uncertain given the complexity of neutrino transfer and angular-momentum transport. After discussing relevant weak-interaction regimes, we study the role of neutrino absorption for shaping Y_e using an extensive set of simulations performed with two-moment neutrino transport and again without neutrino absorption. We vary the torus mass, BH mass and spin, and examine the impact of rest-mass and weak-magnetism corrections in the neutrino rates. We also test the dependence on the angular-momentum transport treatment by comparing axisymmetric models using the standard alpha-viscosity with viscous models assuming constant viscous length scales (l_t) and three-dimensional magnetohydrodynamic (MHD) simulations. Finally, we discuss the nucleosynthesis yields and basic kilonova properties. We find that absorption pushes Y_e towards ~0.5 outside the torus, while inside increasing the equilibrium value Y_e^eq by ~0.05-0.2. Correspondingly, a substantial ejecta fraction is pushed above Y_e=0.25, leading to a reduced lanthanide fraction and a brighter, earlier, and bluer kilonova than without absorption. More compact tori with higher neutrino optical depth, tau, tend to have lower Y_e^eq up to tau~1-10, above which absorption becomes strong enough to reverse this trend. Disk ejecta are less (more) neutron-rich when employing an l_t=const. viscosity (MHD treatment). The solar-like abundance pattern found for our MHD model marginally supports collapsar disks as major r-process sites, although a strong r-process may be limited to phases of high mass-infall rates, Mdot>~ 2 x 10^(-2) Msun/s., 37 pages, 22 figures, 5 tables, accepted to MNRAS, minor changes compared to previous version

Luminous active galactic nuclei (AGN) and X-Ray binaries (XRBs) often contain geometrically thin, radiatively cooled accretion discs. According to theory, these are -- in many cases -- initially highly misaligned with the black hole equator. In this work, we present the first general relativistic magnetohydrodynamic simulations of very thin (h/r~0.015-0.05) accretion discs around rapidly spinning (a~0.9) black holes and tilted by 45-65 degrees. We show that the inner regions of the discs with h/r, 8 pages, 5 figures, accompanying animations included in YouTube playlist: https://www.youtube.com/playlist?list=PLDO1oeU33GwlaPSME1TdCto1Y3P6yG91L

An empirical model is presented that links the demographics of AGN to their ensemble X-ray variability properties. Observations on the incidence of AGN in galaxies are combined with (i) models of the Power Spectrum Density (PSD) of the flux variations of AGN and (ii) parameterisations of the black-hole mass vs stellar-mass scaling relation, to predict the mean excess variance of active black-hole populations in cosmological volumes. We show that the comparison of the model with observational measurements of the ensemble excess variance as a function of X-ray luminosity provides a handle on both the PSD models and the black-hole mass vs stellar mass relation. We find strong evidence against a PSD model that is described by a broken power-law and a constant overall normalisation. Instead our analysis indicates that the amplitude of the PSD depends on the physical properties of the accretion events, such as the Eddington ratio and/or the black hole mass. We also find that current observational measurements of the ensemble excess variance are consistent with the black-hole mass vs stellar mass relation of local spheroids based on dynamically determined black-hole masses. We also discuss future prospects of the proposed approach to jointly constrain the PSD of AGN and the black-hole mass vs stellar mass relation as a function of redshift. The mock AGN catalogue used in the paper is made available. The contents are described in the README file. The code used in the paper is available at https://github.com/ageorgakakis/EnsembleVariability. Full details can be found in the associated paper at https://arxiv.org/abs/2105.08752., The code that uses this data is available at https://github.com/ageorgakakis/EnsembleVariability The paper describing the modeling approach is available at https://arxiv.org/abs/2105.08752, {"references":["https://arxiv.org/abs/2105.08752"]}

The mergers of two neutron stars are typically accompanied by broad-band electromagnetic emission from either a relativistic jet or a kilonova. It has also been long predicted that coherent radio emission will occur during the merger phase or from a newly formed neutron star remnant, however this emission has not been seen to date. This paper presents the deepest limits for this emission from a neutron star merger folowing triggered LOFAR observations of the short gamma-ray burst (SGRB) 181123B, starting 4.4 minutes after the GRB occurred. During the X-ray plateau phase, a signature of ongoing energy injection, we detect no radio emission to a 3$\sigma$ limit of 153 mJy at 144 MHz (image integration time of 136 seconds), which is significantly fainter than the predicted emission from a standard neutron star. At a redshift of 1.8, this corresponds to a luminosity of $2.5 \times 10^{44}$ erg s$^{-1}$. Snapshot images were made of the radio observation on a range of timescales, targeting short duration radio flashes similar to fast radio bursts (FRBs). No emission was detected in the snapshot images at the location of GRB 181123B enabling constraints to be placed on the prompt coherent radio emission model and emission predicted to occur when a neutron star collapses to form a black hole. At the putative host redshift of 1.8 for GRB 181123B, the non detection of the prompt radio emission is two orders of magnitude lower than expected for magnetic reconnection models for prompt GRB emission and no magnetar emission is expected., Comment: MNRAS Accepted, 11 pages, 6 figures

Active galactic nuclei (AGN) forming in the early universe are thought to be the primary source of hard ionizing photons contributing to the reionization of intergalactic helium. However, the number density and spectral properties of high-redshift AGN remain largely unconstrained. In this work, we make use of physically-informed models calibrated with a wide variety of available observations to provide estimates for the role of AGN throughout the Epoch of Reionization. We present AGN luminosity functions in various bands between z = 2 to 7 predicted by the well-established Santa Cruz semi-analytic model, which includes modelling of black hole accretion and AGN feedback. We then combine the predicted AGN populations with a physical spectral model for self-consistent estimates of ionizing photon production rates, which depend on the mass and accretion rate of the accreting supermassive black hole. We then couple the predicted comoving ionizing emissivity with an analytic model to compute the subsequent reionization history of intergalactic helium and hydrogen. This work demonstrates the potential of coupling physically motivated analytic or semi-analytic techniques to capture multi-scale physical processes across a vast range of scales (here, from AGN accretion disks to cosmological scales). Our physical model predicts an intrinsic ionizing photon budget well above many of the estimates in the literature, meaning that helium reionization can comfortably be accomplished even with a relatively low escape fraction. We also make predictions for the AGN populations that are expected to be detected in future \emph{James Webb Space Telescope} surveys., Comment: 28 pages, 23 figures, accepted for publication in MNRAS

We investigate systematically four outbursts of black hole system GX 339-4 observed by the Rossi X-ray Timing Explorer (RXTE) in both spectral and timing domains and find that these outbursts have some common properties although they experience different 'q' tracks in the hardness-intensity diagram (HID). While the spectral indices are around 1.5 in low hard state (LHS), 2.4 in soft intermediate state (SIMS) and high soft state (HSS), the spectral parameters of thermal, non-thermal and reflection components vary significantly in transitions from LHS to HIMS. Also the quasi periodic oscillation (QPO) shows a peculiar behavior during the state transition between LHS and HIMS: the RMS drop of type C fundamental QPO is accompanied with showing-up of the second harmonic. Interestingly, the QPO RMS is found to have a similar linear relationship with the non-thermal fraction of emission in different outbursts. These findings provide more clues to our understanding the outburst of the black hole X-ray binary system., 13 papes, 12 figures

Gravitational Wave (GW) astronomy promises to observe different kinds of astrophysical sources. Here we explore the possibility of detection of GWs from hyperbolic interactions of compact stars with ground-based interferometric detectors. It is believed that a bound compact cluster, such as a globular cluster, can be a primary environment for these interactions. We estimate the detection rates for such events by considering local geometry within the cluster, accounting for scattering probability of compact stars at finite distances, and assuming realistic cluster properties guided by available numerical models, their formation times, and evolution of stars inside them. We find that, even in the conservative limit, it may be possible to detect such black hole encounters in the next few years by the present network of observatories with the ongoing sensitivity upgrades and one to few events per year with the next generation observatories. In practice, actual detection rates can significantly surpass the estimated average rates, since the chances of finding outliers in a very large population can be high. Such observations (or, no observation) may provide crucial constraints to estimate the number of isolated compact stars in the universe. These detections will be exciting discoveries on their own and will be complimentary to observations of binary mergers bringing us one step closer to address a fundamental question, how many black holes are there in the observable universe., 10 pages, 4 figures

Possible formation scenarios of supermassive black holes in the early universe include rapid growth from less massive seed black holes (BHs) via super-Eddington accretion or runaway mergers, yet both of these scenarios would require seed BHs to efficiently sink to and be trapped in the galactic center via dynamical friction. This may not be true for their complicated dynamics in clumpy high-$z$ galaxies. In this work we study this "sinking problem" with state-of-the-art high-resolution cosmological simulations, combined with both direct $N$-body integration of seed BH trajectories and post-processing of randomly generated test particles with a newly developed dynamical friction estimator. We find that seed BHs less massive than $10^8\,M_\odot$ (i.e., all but the already-supermassive seeds) cannot efficiently sink in typical high-$z$ galaxies. We also discuss two possible solutions: dramatically increasing the number of seeds such that one seed can end up trapped in the galactic center by chance, or seed BHs being embedded in dense structures (e.g. star clusters) with effective masses above the mass threshold. We discuss the limitations of both solutions., Comment: accepted to MNRAS

Young stellar clusters across nearly five orders of magnitude in mass appear to follow a power-law mass-radius relationship (MRR), $R_{\star} \propto M_{\star}^{\alpha}$, with $\alpha \approx 0.2 - 0.33$. We develop a simple analytic model for the cluster mass-radius relation. We consider a galaxy disc in hydrostatic equilibrium, which hosts a population of molecular clouds that fragment into clumps undergoing cluster formation and feedback-driven expansion. The model predicts a mass-radius relation of $R_{\star} \propto M_{\star}^{1/2}$ and a dependence on the kpc-scale gas surface density $R_{\star} \propto \Sigma_{\rm g}^{-1/2}$, which results from the formation of more compact clouds (and cluster-forming clumps within) at higher gas surface densities. This environmental dependence implies that the high-pressure environments in which the most massive clusters can form also induce the formation of clusters with the smallest radii, thereby shallowing the observed MRR at high-masses towards the observed $R_{\star} \propto M_{\star}^{1/3}$. At low cluster masses, relaxation-driven expansion induces a similar shallowing of the MRR. We combine our predicted MRR with a simple population synthesis model and apply it to a variety of star-forming environments, finding good agreement. Our model predicts that the high-pressure formation environments of globular clusters at high redshift naturally led to the formation of clusters that are considerably more compact than those in the local Universe, thereby increasing their resilience to tidal shock-driven disruption and contributing to their survival until the present day., Comment: Accepted to MNRAS. Revised model incorporates a new dependence on Toomre Q

We calculate the vorticity generation rate in the accretion disc near a slowly rotating black hole in the low-velocity, weak-field limit of general relativity. Specifically, we find that the frame-dragging effect due to the black hole’s rotation – manifested through the gravitomagnetic field – can generate vorticity in a moving plasma in the accretion disc. The mechanism remains operational as long as the accretion disc has non-negligible vertical height and is independent of the exact thermodynamical profile of the disc. The enstrophy density generation rate, as a measure of turbulence and dissipation, is presented, which indicates that the frame-dragging effect can disrupt the stability of the disc away from the z = 0 plane.

Analysis of timing variation of the eclipsing binary RR Dra implies the existence of an additional object with mass of no less than 3.5(±0.2) M⊙. This third object is orbiting around the central eclipsing pair once every 98(±1) yr, making the whole a hierarchical triple system. However, reliable photometric solutions based on light curves from Transiting Exoplanet Survey Satellite and WASP show that contribution of the third-light takes only about 2 per cent of the total luminosity. It could photometrically be attributed to other unrelated star located within the aperture. The tertiary component is thus a massive object but invisible in optical wavelengths. Besides, evidence of misalignment between the barycenter and the centre of light of the system is also found. This strengthen the existence of a hidden black hole candidate in the form of the third body. The potential black hole may play an essential role in extracting angular momentum from the central binary pair, forming the current state. As a good laboratory, RR Dra is an interesting system that can help to study stellar-mass black hole under the circumbinary case.

X-ray flares in gamma-ray bursts (GRBs) are believed to be generated by the late activities of central engine, and thus provide an useful tool to diagnose the properties of central objects. In this paper, we work on a GRB X-ray flare sample whose bulk Lorentz factors are constrained by two different methods and the jet opening angles are determined by the jet breaks in afterglow lightcurves. Considering a hyperaccreting stellar-mass black hole (BH) as the central engine of GRBs and the Blandford \& Znajek process (BZ) as the jet production mechanism, we constrain the parameters of central engine by using the X-ray flare data. We find that the BZ mechanism is so powerful making it possible to interpret both GRB prompt emissions and bright X-ray flares. The wind parameter ($p$) and accreted mass ($M_d$) fall into reasonable ranges. Our result is also applied to GRB 170817A. The late X-ray flare in GRB 170817A, if it is true, might not be a BH origin., 9 pages, 4 figures, 1 table. Accepted for publication in MNRAS

Relativistic magnetic reconnection is a powerful agent through which magnetic energy can be tapped in astrophysics, energizing particles that then produce observed radiation. In some systems, the highest energy photons come from particles Comptonizing an ambient radiation bath supplied by an external source. If the emitting particle energies are high enough, this inverse Compton (IC) scattering enters the Klein-Nishina regime, which differs from the low-energy Thomson IC limit in two significant ways. First, radiative losses become inherently discrete, with particles delivering an order-unity fraction of their energies to single photons. Second, Comptonized photons may pair-produce with the ambient radiation, opening up another channel for radiative feedback on magnetic reconnection. We analytically study externally illuminated highly magnetized reconnecting systems for which both of these effects are important. We identify a universal (initial magnetization-independent) quasi-steady state in which gamma-rays emitted from the reconnection layer are absorbed in the upstream region, and the resulting hot pairs dominate the energy density of the inflow plasma. However, a true pair cascade is unlikely, and the number density of created pairs remains subdominant to that of the original plasma for a wide parameter range. Future particle-in-cell simulation studies may test various aspects. Pair-regulated Klein-Nishina reconnection may explain steep spectra (quiescent and flaring) from flat-spectrum radio quasars and black hole accretion disc coronae., 43 pages, 22 figures, submitted to MNRAS

Recent detections of gravitational waves from mergers of neutron stars (NSs) and black holes (BHs) in the low and high-end mass gap regimes pose a puzzle to standard stellar and binary evolution theory. Mass-gap mergers may originate from successive mergers in hierarchical systems such as quadruples. Here, we consider repeated mergers of NSs and BHs in stellar 2+2 quadruple systems, in which secular evolution can accelerate the merger of one of the inner binaries. Subsequently, the merger remnant may interact with the companion binary, yielding a second-generation merger. We model the initial stellar and binary evolution of the inner binaries as isolated systems. In the case of successful compact object formation, we subsequently follow the secular dynamical evolution of the quadruple system. When a merger occurs, we take into account merger recoil, and model subsequent evolution using direct N-body integration. With different assumptions on the initial properties, we find that the majority of first-generation mergers are not much affected by secular evolution, with their observational properties mostly consistent with isolated binaries. A small subset shows imprints of secular evolution through residual eccentricity in the LIGO band, and retrograde spin-orbit orientations. Second-generation mergers are ~10^7 times less common than first-generation mergers, and can be strongly affected by scattering (i.e., three-body interactions) induced by the first-generation merger. In particular, scattering can account for mergers within the low-end mass gap, although not the high-end mass gap. Also, in a few cases, scattering could explain highly eccentric LIGO sources and negative effective spin parameters., Accepted for publication in MNRAS (updated references). 38 pages, 91 figures. Includes supplementary material appended to the main manuscript

Local variations in the intergalactic medium (IGM) neutral hydrogen fraction will affect the Ly-$\alpha$ absorption signature of protoclusters identified in tomographic surveys. Using the IllustrisTNG simulations, we investigate how the AGN proximity effect and hot, collisionally ionised gas arising from gravitational infall and black hole feedback changes the Ly-$\alpha$ absorption associated with $M_{z=0}\simeq10^{14}\,M_\odot$ protoclusters at $z\simeq2.4$. We find that protocluster galaxy overdensities exhibit a weak anti-correlation with Ly-$\alpha$ transmission in IGM transmission maps, but local HI ionisation enhancements due to hot $T>10^{6}\rm\,K$ gas or nearby AGN can disrupt this relationship within individual protoclusters. On average, however, we find that strong reductions in the IGM neutral fraction are limited to within $\lesssim 5h^{-1}\,\textrm{cMpc}$ of the dark matter haloes. Local ionisation enhancements will therefore have a minimal impact on the completeness of protocluster identification in tomographic surveys if smoothing Ly-$\alpha$ transmission maps over scales of $\sim4 h^{-1}\,\textrm{cMpc}$, as is typically done in observations. However, if calibrating the relationship between the matter density and Ly-$\alpha$ transmission in tomographic maps using simple analytical models for the Ly-$\alpha$ forest opacity, the presence of hot gas around haloes can still result in systematically lower estimates of $M_{z=0}$ for the most massive protoclusters., Comment: 14 pages (including 1 appendix), 5 figures (plus 2 in the appendix), 1 table. Accepted for publication in MNRAS

The Legacy Survey of Space and Time (LSST) by the Vera C. Rubin Observatory is expected to discover tens of millions of quasars. A significant fraction of these could be powered by coalescing massive black hole (MBH) binaries, since many quasars are believed to be triggered by mergers. We show that under plausible assumptions about the luminosity functions, lifetimes, and binary fractions of quasars, we expect the full LSST quasar catalogue to contain between 20 and 100 million compact MBH binaries with masses M = 105–9M⊙, redshifts z = 0–6, and orbital periods P = 1–70 d. Their light-curves are expected to be distinctly periodic, which can be confidently distinguished from stochastic red-noise variability, because LSST will cover dozens, or even hundreds of cycles. A very small subset of 10–150 ultracompact (P ≲ 1 d) binary quasars among these will, over ∼5–15 yr, evolve into the mHz gravitational-wave frequency band and can be detected by LISA. They can therefore be regarded as ‘LISA verification binaries’, analogous to short-period Galactic compact-object binaries. The practical question is how to find these handful of ‘needles in the haystack’ among the large number of quasars: this will likely require a tailored co-adding analysis optimized for this purpose.

The search for extraterrestrial intelligence (SETI) has been conducted for nearly 60 years. A Dyson Sphere, a spherical structure that surrounds a star and transports its radiative energy outward as an energy source for an advanced civilisation, is one of the main targets of SETI. In this study, we discuss whether building a Dyson Sphere around a black hole is effective. We consider six energy sources: (i) the cosmic microwave background, (ii) the Hawking radiation, (iii) an accretion disk, (iv) Bondi accretion, (v) a corona, and (vi) relativistic jets. To develop future civilisations (for example, a Type II civilisation), $4\times10^{26}\,{\rm W}$($1\,{\rm L_{\odot}}$) is expected to be needed. Among (iii) to (vi), the largest luminosity can be collected from an accretion disk, reaching $10^{5}\,{\rm L_{\odot}}$, enough to maintain a Type II civilisation. Moreover, if a Dyson Sphere collects not only the electromagnetic radiation but also other types of energy (e.g., kinetic energy) from the jets, the total collected energy would be approximately 5 times larger. Considering the emission from a Dyson Sphere, our results show that the Dyson Sphere around a stellar-mass black hole in the Milky Way ($10\,\rm kpc$ away from us) is detectable in the ultraviolet$(\rm 10-400\,{\rm nm)}$, optical$(\rm 400-760\,{\rm nm)}$, near-infrared($\rm 760\,{\rm nm}-5\,{\rm \mu m}$), and mid-infrared($\rm 5-40\,{\rm \mu m}$) wavelengths via the waste heat radiation using current telescopes such as Galaxy Evolution Explorer Ultraviolet Sky Surveys. Performing model fitting to observed spectral energy distributions and measuring the variability of radial velocity may help us to identify these possible artificial structures., Comment: This paper has been accepted for publication in MNRAS

We present multi-wavelength fast timing observations of the black hole X-ray binary MAXI J1820+070 (ASASSN-18ey), taken with the Karl G. Jansky Very Large Array (VLA), Atacama Large Millimeter/Sub-Millimeter Array (ALMA), Very Large Telescope (VLT), New Technology Telescope (NTT), Neutron Star Interior Composition Explorer (NICER), and XMM-Newton. Our data set simultaneously samples ten different electromagnetic bands (radio - X-ray) over a 7-hour period during the hard state of the 2018-2019 outburst. The emission we observe is highly variable, displaying multiple rapid flaring episodes. To characterize the variability properties in our data, we implemented a combination of cross-correlation and Fourier analyses. We find that the emission is highly correlated between different bands, measuring time-lags ranging from hundreds of milliseconds between the X-ray/optical bands to minutes between the radio/sub-mm bands. Our Fourier analysis also revealed, for the first time in a black hole X-ray binary, an evolving power spectral shape with electromagnetic frequency. Through modelling these variability properties, we find that MAXI J1820+070 launches a highly relativistic ($\Gamma=6.81^{+1.06}_{-1.15}$) and confined ($\phi=0.45^{+0.13}_{-0.11}$ deg) jet, which is carrying a significant amount of power away from the system (equivalent to $\sim0.6 \, L_{1-100{\rm keV}}$). We additionally place constraints on the jet composition and magnetic field strength in the innermost jet base region. Overall, this work demonstrates that time-domain analysis is a powerful diagnostic tool for probing jet physics, where we can accurately measure jet properties with time-domain measurements alone., Comment: 23 pages, 16 figures, accepted to MNRAS

Jets launched by the supermassive black holes in the centers of cool-core clusters are the most likely heat source to solve the cooling flow problem. One way for this heating to occur is through generation of a turbulent cascade by jet-inflated bubbles. Measurements of the X-ray intensity power spectra show evidence of this cascade in different regions of the cluster, constraining the role of driving mechanisms. We analyze feedback simulations of the Perseus cluster to constrain the effect of the jet activity on the intensity fluctuations and kinematics of the cluster atmosphere. We find that, within the inner 60 kiloparsecs, the power spectra of the predicted surface brightness fluctuations are broadly consistent with those measured by Chandra and that even a single episode of jet activity can generate a long-lasting imprint on the intensity fluctuations in the innermost region of the cluster. AGN-driven motions within the same region approach the values reported by Hitomi during and right after the AGN episode. However, the line-of-sight velocity dispersion excited by the jet in simulations underpredicts the Hitomi measurement. This indicates that driving a volume-filling sustained level of turbulence requires several episodes of jet activity, and/or additional processes drive turbulence outside the 60-kpc sphere. This also suggests that sharp edges of the bubbles in the innermost region of the cluster contribute substantially to the intensity of fluctuations, consistent with the Perseus observations in the inner 30-kpc region. We discuss new diagnostics to decompose annular power spectra to constrain past episodes of jet activity., Accepted for publication in MNRAS

Relativistic jets and disc-winds are typically observed in BH-XRBs and AGNs. However, many physical details of jet launching and the driving of disc winds from the underlying accretion disc are still not fully understood. In this study, we further investigate the role of the magnetic field strength and structure in launching jets and disc winds. In particular, we explore the connection between jet, wind, and the accretion disc around the central black hole. We perform axisymmetric GRMHD simulations of the accretion-ejection system using adaptive mesh refinement. Essentially, our simulations are initiated with a thin accretion disc in equilibrium. An extensive parametric study by choosing different combinations of magnetic field strength and initial magnetic field inclination is also performed. Our study finds relativistic jets driven by the Blandford \& Znajek (BZ) mechanism and the disc-wind driven by the Blandford \& Payne (BP) mechanism. We also find that plasmoids are formed due to the reconnection events, and these plasmoids advect with disc-winds. As a result, the tension force due to the poloidal magnetic field is enhanced in the inner part of the accretion disc, resulting in disc truncation and oscillation. These oscillations result in flaring activities in the jet mass flow rates. We find simulation runs with a lower value of the plasma-$��$, and lower inclination angle parameters are more prone to the formation of plasmoids and subsequent inner disc oscillations. Our models provide a possible template to understand spectral state transition phenomena in BH-XRBs., 21 pages, 20 figures, accepted for publication in MNRAS

We show that the odds of the mass-gap (secondary) object in GW190814 being a neutron star (NS) improve if one allows for a stiff high-density equation of state (EoS) or a large spin. Since its mass is $\in (2.50,2.67) M_{\odot}$, establishing its true nature will make it either the heaviest NS or the lightest black hole (BH), and can have far-reaching implications on NS EoS and compact object formation channels. When limiting oneself to the NS hypothesis, we deduce the secondary's properties by using a Bayesian framework with a hybrid EoS formulation that employs a parabolic expansion-based nuclear empirical parameterization around the nuclear saturation density augmented by a generic 3-segment piecewise polytrope (PP) model at higher densities and combining a variety of astrophysical observations. For the slow-rotation scenario, GW190814 implies a very stiff EoS and a stringent constraint on the EoS specially in the high-density region. On the other hand assuming the secondary object is a rapidly rotating NS, we constrain its rotational frequency to be $f=1170^{+389}_{-495}$ Hz, within a $90\%$ confidence interval. In this scenario, the secondary object in GW190814 would qualify as the fastest rotating NS ever observed. However, for this scenario to be viable, rotational instabilities would have to be suppressed both during formation and the subsequent evolution until merger, otherwise the secondary of GW190814 is more likely to be a BH., Comment: Matches with the MNRAS published version

The circular polarization of black hole accretion flows can encode properties of the underlying magnetic field structure. Using general relativistic magnetohydrodynamics (GRMHD) simulations, we study the imprint of magnetic field geometry on circular polarization images potentially observable by the Event Horizon Telescope (EHT). We decompose images into the different mechanisms that generate circular polarization in these models, which are sensitive to both the line of sight direction and twist of the magnetic field. In these models, a stable sign of the circular polarization over time, as observed for several sources, can be attributed to a stability of these properties. We illustrate how different aspects of a generic helical magnetic field geometry become imprinted on a circular polarization image. We also identify novel effects of light bending that affect the circular polarization image on event horizon scales. One consequence is the sign flipping of successive photon rings in face-on systems, which if observable and uncorrupted by Faraday rotation, can directly encode the handedness of the approaching magnetic field., Comment: 20 pages, 19 figures, in review with MNRAS

I propose a new type of common envelope jets supernova (CEJSN) events where instead of a single neutron star (NS; or a black hole; BH) a tight binary system of a NS and a main sequence star enters a common envelope evolution (CEE) with a red supergiant. The NS and the main sequence star of the tight binary system merge inside the red supergiant envelope and enter a CEE of their own. The NS accretes some mass through an accretion disk and launches jets that explodes the main sequence star. I estimate that the two jets that the NS launches at this phase carry an energy of ~10^{52} erg, about the same order of magnitude as the energy that the jets will carry when the NS or its BH remnant will enter the core in a later phase. For that, I term the entire event a double CEJSN. The outcome of the double CEJSN is a very long, months to years, and very energetic event, a total energy of ~10^{52} - 10^{53} erg, that will be observationally classified as a peculiar super-energetic event. I crudely estimate that new transient surveys should detect about one CEJSN event from a triple star system per year., Accepted for publication in Monthly Notices of the Royal Astronomical Society

The recently obtained hairy Kerr black holes, due to additional sources or surrounding fluid, like dark matter, with conserved energy-momentum tensor, have a deviation $\alpha$ and primary hair $l_0$, apart from rotation parameter $a$ and mass $M$. In the wake of the \textit{Event Horizon Telescope} (\textit{EHT}) observations of the supermassive black hole M87*, a recent surge in interest in black hole shadows suggests comparing the black holes in general relativity (GR) and modified theories of gravity (MoG) to assess these models' differences. Motivated by this, we take on an extensive study of the rotating hairy Kerr black holes, which encompasses, in particular cases, the Kerr black hole ($\alpha=0$). We investigate ergosphere and shadows of the black holes to infer that their size and shape are affected due to the $l_0$ and are found to harbour a richer chaotic structure. In particular, the hairy Kerr black holes possess smaller size but more distorted shadows when compared with Kerr black holes. We also estimate the parameters $l_0$ and $a$ associated with hairy Kerr black holes using the shadow observables. The inferred circularity deviation $\Delta C \leq 0.1$ for the M87* black hole is satisfied, whereas shadow angular diameter $\theta_{d}=42 \pm 3 \mu as$, within $1 \sigma$ region, for a given choice of $\alpha$, places bounds on the parameters $a$ and $l_0$. Interestingly, the shadow axial ratio obeying $1< D_x \lesssim 4/3$ is in agreement with the \textit{EHT} results and thus eventuates in the hairy Kerr black holes being suitable candidates for astrophysical black holes., Comment: 15 pages, 13 figures, 2 tables. Matched to the published version

The gamma-ray emitting galaxy SBS 0846+513 has been classified as a Narrow-Line Seyfert 1 from its spectroscopy, and on that basis it was thought likely to have a small central black hole hosted in a spiral galaxy. But very few of the gamma-ray Narrow-Line Seyfert 1s have high-resolution imaging of their hosts, so it is unknown how those expectations hold up for the gamma-emitting class. We have observed this galaxy in the J-band with the Large Binocular Telescope's LUCI1 camera and the ARGOS adaptive optics system. We estimate its black hole mass to lie between $7.70 \leq \log \frac{\text{M}}{\text{M}_\odot} \leq 8.19$, using the correlation with bulge luminosity, or $7.96 \leq \log \frac{\text{M}}{\text{M}_\odot} \leq 8.16$ using the correlation with S\'{e}rsic index, putting its mass at the high end of the Narrow Line Seyfert 1 range. These estimates are independent of the Broad Line Region viewing geometry and avoid underestimates due to looking down the jet axis. Its host shows evidence of a bulge + disc structure, both from two-dimensional modeling and isophote shape, in keeping with the expectations. Mergers and interactions appear to be common among the gamma-ray Narrow-Line Seyfert 1s, and we see some circumstantial evidence for companion galaxies or disturbed features in the host., Comment: 10 pages, 6 figures. Submitted to MNRAS

We report on the Insight-HXMT observations of the new black hole X-ray binary MAXI J1820+070 during its 2018 outburst. Detailed spectral analysis via the continuum fitting method shows an evolution of the inferred spin during its high soft sate. Moreover, the hardness ratio, the non-thermal luminosity and the reflection fraction also undergo an evolution, exactly coincident to the period when the inferred spin transition takes place. The unphysical evolution of the spin is attributed to the evolution of the inner disc, which is caused by the collapse of a hot corona due to condensation mechanism or may be related to the deceleration of a jet-like corona. The studies of the inner disc radius and the relation between the disc luminosity and the inner disc radius suggest that, only at a particular epoch, did the inner edge of the disc reach the innermost stable circular orbit and the spin measurement is reliable. We then constrain the spin of MAXI J1820+070 to be a*=0.2^{+0.2}_{-0.3}. Such a slowly spinning black hole possessing a strong jet suggests that its jet activity is driven mainly by the accretion disc rather than by the black hole spin., Comment: 14 pages, 13 figures, 5 tables, accepted for publication in MNRAS

We present sub-arcsecond resolution ALMA imaging of the CO(3-2) emission in two $z\sim2.5$ heavily reddened quasars (HRQs) - ULASJ1234+0907 and ULASJ2315+0143 - and their companion galaxies. Dynamical modeling of the resolved velocity fields enables us to constrain the molecular gas morphologies and host galaxy masses. Combining the new data with extensive multi-wavelength observations, we are able to study the relative kinematics of different molecular emission lines, the molecular gas fractions and the locations of the quasars on the M$_{\rm{BH}}$-M$_{\rm{gal}}$ relation. Despite having similar black-hole properties, the two HRQs display markedly different host galaxy properties and local environments. J1234 has a very massive host, M$_{\rm{dyn}} \sim 5 \times 10^{11}$M$_\odot$ and two companion galaxies that are similarly massive located within 200 kpc of the quasar. The molecular gas fraction is low ($\sim$6%). The significant ongoing star formation in the host galaxy is entirely obscured at rest-frame UV and optical wavelengths. J2315 is resolved into a close-separation major-merger ($\Delta$r=15 kpc; $\Delta$v=170 km/s) with a $\sim$1:2 mass ratio. The total dynamical mass is estimated to be $\lesssim$10$^{11}$M$_\odot$ and the molecular gas fraction is high ($>$45%). A new HSC image of the galaxy shows unobscured UV-luminous star-forming regions co-incident with the extended reservoir of cold molecular gas in the merger. We use the outputs from the Illustris simulations to track the growth of such massive black holes from $z\sim6$ to the present day. While J1234 is consistent with the simulated $z\sim2$ relation, J2315 has a black hole that is over-massive relative to its host galaxy., Comment: 18 pages, 12 figures, Accepted for publication in MNRAS

Quiescent black hole X-ray binaries (X-ray luminosities 25 kpc, while initial optical astrometric measurements from Gaia Data Release 2 suggested likely distances of just a few kpc. Here, we use the most recent measurements from Gaia Early Data Release 3 and find a distance d=7.1(+4.8/-3.9) kpc and a potential kick velocity PKV=165(+81/-17) km/s, with distances up to ~20 kpc possible based on its parallax and proper motion. Even though there is now less tension between the parallax and donor-star based distance measurements, it remains an unresolved matter, and we conclude with suggestions on how to reconcile the two measurements., 13 pages, 5 figures, 1 table. Accepted for publication in MNRAS

Binary supermassive black hole (SMBH) systems result from galaxy mergers, and will eventually coalesce due to gravitational wave (GW) emission if the binary separation can be reduced to $\lesssim 0.1$ pc by other mechanisms. Here, we explore a gravitational self-lensing binary SMBH model for the sharp (duration $\sim 1$ hr), quasi-regular X-ray flares -- dubbed quasi-periodic eruptions -- recently observed from two low mass active galactic nuclei: GSN 069 and RX J1301.9+2747. In our model, the binary is observed $\sim$edge-on, such that each SMBH gravitationally lenses light from the accretion disc surrounding the other SMBH twice per orbital period. The model can reproduce the flare spacings if the current eccentricity of RX J1301.9+2747 is $\epsilon_0 \gtrsim 0.16$, implying a merger within $\sim 1000$ yrs. However, we cannot reproduce the observed flare profiles with our current calculations. Model flares with the correct amplitude are $\sim 2/5$ the observed duration, and model flares with the correct duration are $\sim 2/5$ the observed amplitude. Our modelling yields three distinct behaviours of self-lensing binary systems that can be searched for in current and future X-ray and optical time-domain surveys: i) periodic lensing flares, ii) partial eclipses (caused by occultation of the background mini-disc by the foreground mini-disc), and iii) partial eclipses with a very sharp in-eclipse lensing flare. Discovery of such features would constitute very strong evidence for the presence of a supermassive binary, and monitoring of the flare spacings will provide a measurement of periastron precession., Comment: Accepted for publication in MNRAS

During outbursts, the observational properties of black hole X-ray binaries (BHXBs) vary on timescales of days to months. These relatively short timescales make these systems ideal laboratories to probe the coupling between accreting material and outflowing jets as a the accretion rate varies. In particular, the origin of the hard X-ray emission is poorly understood and highly debated. This spectral component, which has a power-law shape, is due to Comptonisation of photons near the black hole, but it is unclear whether it originates in the accretion flow itself, or at the base of the jet, or possibly the interface region between them. In this paper we explore the disk-jet connection by modelling the multi-wavelength emission of MAXI J1836-194 during its 2011 outburst. We combine radio through X-ray spectra, X-ray timing information, and a robust joint-fitting method to better isolate the jet's physical properties. Our results demonstrate that the jet base can produce power-law hard X-ray emission in this system/outburst, provided that its base is fairly compact and that the temperatures of the emitting electrons are sub-relativistic. Because of energetic considerations, our model favours mildly pair-loaded jets carrying at least 20 pairs per proton. Finally, we find that the properties of the X-ray power spectrum are correlated with the jet properties, suggesting that an underlying physical process regulates both., 17 pages, 10 figures, accepted for publication on MNRAS

The past decade has seen significant progress in understanding galaxy formation and evolution using large-scale cosmological simulations. While these simulations produce galaxies in overall good agreement with observations, they employ different sub-grid models for galaxies and supermassive black holes (BHs). We investigate the impact of the sub-grid models on the BH mass properties of the Illustris, TNG100, TNG300, Horizon-AGN, EAGLE, and SIMBA simulations, focusing on the M_BH-M_star relation and the BH mass function. All simulations predict tight M_BH-M_star relations, and struggle to produce the lowest (M_BH< 10^7.5 Msun) in galaxies of M_star~10^10.5-10^11.5 Msun. While the time evolution of the mean M_BH-M_star relation is mild (10^9 Msun in most of the simulations. The BH mass function is dominated by efficiently accreting BHs (log10 f_Edd >-2$) at high redshifts, and transitions progressively from the high-mass to the low-mass end to be governed by inactive BHs. The transition time and the contribution of active BHs are different among the simulations, and can be used to evaluate models against observations., Comment: Published version - MNRAS

Some ultra-compact dwarf galaxies (UCDs) have elevated observed dynamical V-band mass-to-light ($M/L_{\rm V}$) ratios with respect to what is expected from their stellar populations assuming a canonical initial mass function (IMF). Observations have also revealed the presence of a compact dark object in the centers of several UCDs, having a mass of a few to 15\% of the present-day stellar mass of the UCD. This central mass concentration has typically been interpreted as a super-massive black hole, but can in principle also be a sub-cluster of stellar remnants. We explore the following two formation scenarios of UCDs, i) monolithic collapse and ii) mergers of star clusters in cluster complexes as are observed in massively star-bursting regions. We explore the physical properties of the UCDs at different evolutionary stages assuming different initial stellar masses of the UCDs and the IMF being either universal or changing systematically with metallicity and density according to the Integrated Galactic IMF (IGIMF) theory. While the observed elevated $M/L_{\rm V}$ ratios of the UCDs cannot be reproduced if the IMF is invariant and universal, the empirically derived IMF which varies systematically with density and metallicity shows agreement with the observations. Incorporating the UCD-mass-dependent retention fraction of dark remnants improves this agreement. In addition we apply the results of N-body simulations to young UCDs and show that the same initial conditions describing the observed $M/L_{\rm V}$ ratios reproduce the observed relation between the half-mass radii and the present-day masses of the UCDs. The findings thus suggest that the majority of UCDs that have elevated $M/L_{\rm V}$ ratios could have formed monolithically with significant remnant-mass components that are centrally concentrated, while those with small $M/L_V$ values may be merged star-cluster complexes., Accepted for publication in MNRAS. 16 pages, 11 figures, 1 table

We present 3D general relativistic magnetohydrodynamic(GRMHD) simulations of zero angular momentum accretion around a rapidly rotating black hole, modified by the presence of initially uniform magnetic fields. We consider serveral angles between the magnetic field direction and the black hole spin. In the resulting flows, the midplane dynamics are governed by magnetic reconnection-driven turbulence in a magnetically arrested (or a nearly arrested) state. Electromagnetic jets with outflow efficiencies ~10-200% occupy the polar regions, reaching several hundred gravitational radii before they dissipate due to the kink instability. The jet directions fluctuate in time and can be tilted by as much as ~30 degrees with respect to black hole spin, but this tilt does not depend strongly on the tilt of the initial magnetic field. A jet forms even when there is no initial net vertical magnetic flux since turbulent, horizon-scale fluctuations can generate a net vertical field locally. Peak jet power is obtained for an initial magnetic field tilted by 40-80 degrees with respect to the black hole spin because this maximizes the amount of magnetic flux that can reach the black hole. These simulations may be a reasonable model for low luminosity black hole accretion flows such as Sgr A* or M87., Comment: Accepted by MNRAS. For animations see [https://smressle.bitbucket.io/animations.html]