Search

With its capacity to observe $\sim 10^{5-6}$ faint active galactic nuclei (AGN) out to redshift $z\approx 6$, Roman is poised to reveal a population of $10^{4-6}\, {\rm M_\odot}$ black holes during an epoch of vigorous galaxy assembly. By measuring the light curves of a subset of these AGN and looking for periodicity, Roman can identify several hundred massive black hole binaries (MBHBs) with 5-12 day orbital periods, which emit copious gravitational radiation and will inevitably merge on timescales of $10^{3-5}$ years. During the last few months of their merger, such binaries are observable with the Laser Interferometer Space Antenna (LISA), a joint ESA/NASA gravitational wave mission set to launch in the mid-2030s. Roman can thus find LISA precursors, provide uniquely robust constraints on the LISA source population, help identify the host galaxies of LISA mergers, and unlock the potential of multi-messenger astrophysics with massive black hole binaries., Comment: White Paper for the Nancy Grace Roman Space Telescope's Core Community Surveys (https://roman.gsfc.nasa.gov/science/ccs_white_papers.html)

We present a study of optically selected dual AGN with projected separations of 3--97~kpc. Using multi-wavelength (MWL) information (optical, X-rays, mid-IR), we characterized the intrinsic nuclear properties of this sample and compared them with those of isolated systems. Among the 124 X-ray detected AGN candidates, 52 appear in pairs and 72 as single X-ray sources. Through MWL analysis, we confirmed the presence of the AGN in a fraction >80\% of the detected targets in pairs (42 over 52). X-ray spectral analysis confirms the trend of increasing AGN luminosity with decreasing separation, suggesting that mergers may have contributed in triggering more luminous AGN. Through X/mid-IR ratio $vs$ X-ray colors, we estimated a fraction of Compton-thin AGN (with 10$^{22}$ cm$^{-2}$ $<$ N$_{\rm H} <$10$^{24}$ cm$^{-2}$) of about 80\%, while about 16\% are Compton thick (CT, with N$_{\rm H}>$10$^{24}$ cm$^{-2}$) sources. These fractions of obscured sources are larger than those found in samples of isolated AGN, confirming that pairs of AGN show higher obscuration. This trend is further confirmed by comparing the de-reddened [O\ III] emission with the observed X-ray luminosity. However, the derived fraction of Compton-thick sources in this sample at early stage of merging is lower than reported for late-merging dual-AGN samples. Comparing N$_{\rm H}$ from X-rays with that derived from E(B-V) from Narrow Line Regions, we find that the absorbing material is likely associated with the torus or the Broad Line Regions. We also explored the X-ray detection efficiency of dual-AGN candidates, finding that, when observed properly (at on-axis positions and with long exposures), X-ray data represent a powerful way to confirm and investigate dual-AGN systems., Comment: 12 pages, 4 figures. Accepted for publication in MNRAS

Dual active galaxy nuclei (dAGNs) trace the population of post-merger galaxies and are the precursors to massive black hole (MBH) mergers, an important source of gravitational waves that may be observed by LISA. In Paper I of this series, we used the population of nearly 2000 galaxy mergers predicted by the TNG50-3 simulation to seed semi-analytic models of the orbital evolution and coalescence of MBH pairs with initial separations of about 1 kpc. Here, we calculate the dAGN luminosities and separation of these pairs as they evolve in post-merger galaxies, and show how the coalescence fraction of dAGNs changes with redshift. We find that because of the several Gyr long dynamical friction timescale for orbital evolution, the fraction of dAGNs that eventually end in a MBH merger grows with redshift and does not pass 50% until a redshift of 1. However, dAGNs in galaxies with bulge masses >10^10 solar masses, or comprised of near-equal mass MBHs, evolve more quickly and have higher than average coalescence fractions. At any redshift, dAGNs observed with small separations (> 0.7 kpc) have a higher probability of merging in a Hubble time than more widely separated systems. As found in Paper I, radiation feedback effects can significantly reduce the number of MBH mergers, and this could be manifested as a larger than expected number of widely separated dAGNs. We present a method to estimate the MBH coalescence rate as well as the potential LISA detection rate given a survey of dAGNs. Comparing these rates to the eventual LISA measurements will help determine the efficiency of dynamical friction in post-merger galaxies., Comment: Submitted to ApJ, 13 pages, 8 figures

The Laser Interferometer Space Antenna (LISA) will be a transformative experiment for gravitational wave astronomy, and, as such, it will offer unique opportunities to address many key astrophysical questions in a completely novel way. The synergy with ground-based and space-born instruments in the electromagnetic domain, by enabling multi-messenger observations, will add further to the discovery potential of LISA. The next decade is crucial to prepare the astrophysical community for LISA's first observations. This review outlines the extensive landscape of astrophysical theory, numerical simulations, and astronomical observations that are instrumental for modeling and interpreting the upcoming LISA datastream. To this aim, the current knowledge in three main source classes for LISA is reviewed; ultracompact stellar-mass binaries, massive black hole binaries, and extreme or intermediate mass ratio inspirals. The relevant astrophysical processes and the established modeling techniques are summarized. Likewise, open issues and gaps in our understanding of these sources are highlighted, along with an indication of how LISA could help making progress in the different areas. New research avenues that LISA itself, or its joint exploitation with upcoming studies in the electromagnetic domain, will enable, are also illustrated. Improvements in modeling and analysis approaches, such as the combination of numerical simulations and modern data science techniques, are discussed. This review is intended to be a starting point for using LISA as a new discovery tool for understanding our Universe.

The hot intracluster medium (ICM) is thought to be quiescent with low observed velocity dispersions. Surface brightness fluctuations of the ICM also suggest that its turbulence is subsonic with a Kolmogorov scaling relation, indicating that the viscosity is suppressed and the kinetic energy cascades to small scales unscathed. However, recent observations of the cold gas filaments in galaxy clusters find that the scaling relations are steeper than that of the hot plasma, signaling kinetic energy losses and the presence of supersonic flows. In this work we use high-resolution simulations to explore the turbulent velocity structure of the cold filaments at the cores of galaxy clusters. Our results indicate that supersonic turbulent structures can be "frozen" in the cold gas that cools and fragments out of a fast, $10^7$ K outflow driven by the central active galactic nucleus (AGN), when the radiative cooling time is shorter than the dynamical sound-crossing time. After the cold gas formation, however, the slope of the velocity structure function (VSF) flattens significantly over short, 10 Myr timescales. The lack of flattened VSF in observations of H$\alpha$ filaments indicates that the H$\alpha$-emitting phase is short-lived for the cold gas in galaxy clusters. On the other hand, the ubiquity of supersonic turbulence revealed by observed filaments strongly suggests that supersonic outflows are an integral part of AGN-ICM interaction, and that AGN activity plays a crucial role at driving turbulence in galaxy clusters., Comment: 9 pages, 4 figures, accepted for publication in ApJL

We evaluate the cosmological coalescence and detection rates for massive black hole (MBH) binaries targeted by the gravitational wave observatory Laser Interferometer Space Antenna (LISA). Our calculation starts with a population of gravitationally unbound MBH pairs, drawn from the TNG50-3 cosmological simulation, and follows their orbital evolution from kpc scales all the way to coalescence using a semi-analytic model developed in our previous work. We find that for a majority of MBH pairs that coalesce within a Hubble time dynamical friction is the most important mechanism that determines their coalescence rate. Our model predicts a MBH coalescence rate < 0.45/ yr and a LISA detection rate < 0.34/ yr. Most LISA detections should originate from 10^6 - 10^6.8 solar masses MBHs in gas-rich galaxies at redshifts 1.6 < z < 2.4, and have a characteristic signal to noise ratio SNR ~ 100. We however find a dramatic reduction in the coalescence and detection rates, as well as the average SNR, if the effects of radiative feedback from accreting MBHs are taken into account. In this case, the MBH coalescence rate is reduced by 78% (to < 0.1/ yr), and the LISA detection rate is reduced by 94% (to 0.02/ yr), whereas the average SNR is ~ 10. We emphasize that our model provides a lower limit on the LISA detection rate, consistent with other works in the literature that draw their MBH pairs from cosmological simulations., Comment: Accepted to ApJ, 20 pages, 10 figures

Supermassive black hole binaries (SMBHBs) are a natural outcome of galaxy mergers and should form frequently in galactic nuclei. Sub-parsec binaries can be identified from their bright electromagnetic emission, e.g., Active Galactic Nuclei (AGN) with Doppler shifted broad emission lines or AGN with periodic variability, as well as from the emission of strong gravitational radiation. The most massive binaries (with total mass >10^8 M_sol) emit in the nanohertz band and are targeted by Pulsar Timing Arrays (PTAs). Here we examine the synergy between electromagnetic and gravitational wave signatures of SMBHBs. We connect both signals to the orbital dynamics of the binary and examine the common link between them, laying the foundation for joint multi-messenger observations. We find that periodic variability arising from relativistic Doppler boost is the most promising electromagnetic signature to connect with GWs. We delineate the parameter space (binary total mass/chirp mass versus binary period/GW frequency) for which joint observations are feasible. Currently multi-messenger detections are possible only for the most massive and nearby galaxies, limited by the sensitivity of PTAs. However, we demonstrate that as PTAs collect more data in the upcoming years, the overlapping parameter space is expected to expand significantly., Comment: submitted to MNRAS

The next two decades are expected to open the door to the first coincident detections of electromagnetic (EM) and gravitational wave (GW) signatures associated with massive black hole (MBH) binaries heading for coalescence. These detections will launch a new era of multimessenger astrophysics by expanding this growing field to the low-frequency GW regime and will provide an unprecedented understanding of the evolution of MBHs and galaxies. They will also constitute fundamentally new probes of cosmology and would enable unique tests of gravity. The aim of this Living Review is to provide an introduction to this research topic by presenting a summary of key findings, physical processes and ideas pertaining to EM counterparts to MBH mergers as they are known at the time of this writing. We review current observational evidence for close MBH binaries, discuss relevant physical processes and timescales, and summarize the possible EM counterparts to GWs in the precursor, coalescence, and afterglow stages of a MBH merger. We also describe open questions and discuss future prospects in this dynamic and quick-paced research area., Comment: Invited review article accepted for publication in Living Reviews in Relativity. 127 pages, 31 figures. The revised version includes reviewers' comments, updated references, new Table 1, and new Section 4.4

Warm ionized and cold neutral outflows with velocities exceeding $100\,{\rm km\,s}^{-1}$ are commonly observed in galaxies and clusters. Theoretical studies however indicate that ram pressure from a hot wind, driven either by the central active galactic nucleus (AGN) or a starburst, cannot accelerate existing cold gas to such high speeds without destroying it. In this work we explore a different scenario, where cold gas forms in a fast, radiatively cooling outflow with temperature $T\lesssim 10^7\,{\rm K}$. Using 3D hydrodynamic simulations, we demonstrate that cold gas continuously fragments out of the cooling outflow, forming elongated filamentary structures extending tens of kiloparsecs. For a range of physically relevant temperature and velocity configurations, a ring of cold gas perpendicular to the direction of motion forms in the outflow. This naturally explains the formation of transverse cold gas filaments such as the blue loop and the horseshoe filament in the Perseus cluster. Based on our results, we estimate that the AGN outburst responsible for the formation of these two features drove bipolar outflows with velocity $>2,000\,{\rm km\,s}^{-1}$ and total kinetic energy $>8\times10^{57}\,{\rm erg}$ about $\sim10$ Myr ago. We also examine the continuous cooling in the mixing layer between hot and cold gas, and find that radiative cooling only accounts for $\sim10\%$ of the total mass cooling rate, indicating that observations of soft X-ray and FUV emission may significantly underestimate the growth of cold gas in the cooling flow of galaxy clusters., Comment: 8 pages, 3 figures, accepted for publication in ApJL

The science objectives of the LISA mission have been defined under the implicit assumption of a 4 yr continuous data stream. Based on the performance of LISA Pathfinder, it is now expected that LISA will have a duty cycle of $\approx 0.75$, which would reduce the effective span of usable data to 3 yr. This paper reports the results of a study by the LISA Science Group, which was charged with assessing the additional science return of increasing the mission lifetime. We explore various observational scenarios to assess the impact of mission duration on the main science objectives of the mission. We find that the science investigations most affected by mission duration concern the search for seed black holes at cosmic dawn, as well as the study of stellar-origin black holes and of their formation channels via multi-band and multi-messenger observations. We conclude that an extension to 6 yr of mission operations is recommended., Comment: 50 pages, 19 figures, 5 tables. Matches version published in GERG

In several previous studies, quasars exhibiting broad emission lines with >1000 km/s velocity offsets with respect to the host galaxy rest frame have been discovered. One leading hypothesis for the origin of these velocity-offset broad lines is the dynamics of a binary supermassive black hole (SMBH). We present high-resolution radio imaging of 34 quasars showing these velocity-offset broad lines with the Very Long Baseline Array (VLBA), aimed at finding evidence for the putative binary SMBHs (such as dual radio cores), and testing the competing physical models. We detect exactly half of the target sample from our VLBA imaging, after implementing a 5 detection limit. While we do not resolve double radio sources in any of the targets, we obtain limits on the instantaneous projected separations of a radio-emitting binary for all of the detected sources under the assumption that a binary still exists within our VLBA angular resolution limits. We also assess the likelihood that a radio-emitting companion SMBH exists outside of our angular resolution limits, but its radio luminosity is too weak to produce a detectable signal in the VLBA data. Additionally, we compare the precise sky positions afforded by these data to optical positions from both the SDSS and Gaia DR2 source catalogs. We find projected radio/optical separations on the order of 10 pc for three quasars. Finally, we explore how future multi-wavelength campaigns with optical, radio, and X-ray observatories can help discriminate further between the competing physical models., Comment: Accepted for publishing by ApJ

While theory and simulations indicate that galaxy mergers play an important role in the cosmological evolution of accreting black holes and their host galaxies, samples of Active Galactic Nuclei (AGN) in galaxies at close separations are still small. In order to increase the sample of AGN pairs, we undertook an archival project to investigate the X-ray properties of a SDSS-selected sample of 32 galaxy pairs with separations $\le$150 kpc containing one optically-identified AGN, that were serendipitously observed by XMM-Newton. We discovered only one X-ray counterpart among the optically classified non-active galaxies, with a weak X-ray luminosity ($\simeq$5$\times$10$^{41}$ erg s$^{-1}$). 59% (19 out of 32) of the AGN in our galaxy pair sample exhibit an X-ray counterpart, covering a wide range in absorption corrected X-ray luminosity (5$\times$10$^{41}$-2$\times$10$^{44}$ erg s$^{-1}$). More than 79% of these AGN are obscured (column density $N_H>$10$^{22}$ cm$^{-2}$), with more than half thereof ({\it i.e.}, about 47% of the total AGN sample) being Compton-thick. AGN/no-AGN pairs are therefore more frequently X-ray obscured (by a factor $\simeq$1.5) than isolated AGN. When compared to a luminosity and redshift-matched sample of {\it bona fide} dual AGN, AGN/no-AGN pairs exhibit one order-of-magnitude lower X-ray column density in the same separation range ($>$10 kpc). A small sample (4 objects) of AGN/no-AGN pairs with sub-pc separation are all heavily obscured, driving a formal anti-correlation between the X-ray column density and the galaxy pair separation in these systems. These findings suggest that the galactic environment has a key influence on the triggering of nuclear activity in merging galaxies., Comment: Accepted on MNRAS, 17 pages, 9 figures

Outflows driven by active galactic nuclei (AGN) are an important channel for accreting supermassive black holes (SMBHs) to interact with their host galaxies and clusters. Properties of the outflows are however poorly constrained due to the lack of kinetically resolved data of the hot plasma that permeates the circumgalactic and intracluster space. In this work, we use a single parameter, outflow-to-accretion mass-loading factor $m=\dot{M}_{\rm out}/\dot{M}_{\rm BH}$, to characterize the outflows that mediate the interaction between SMBHs and their hosts. By modeling both M87 and Perseus, and comparing the simulated thermal profiles with the X-ray observations of these two systems, we demonstrate that $m$ can be constrained between $200-500$. This parameter corresponds to a bulk flow speed between $4,000-7,000\,{\rm km\,s}^{-1}$ at around 1 kpc, and a thermalized outflow temperature between $10^{8.7}-10^{9}\,{\rm K}$. Our results indicate that the dominant outflow speeds in giant elliptical galaxies and clusters are much lower than in the close vicinity of the SMBH, signaling an efficient coupling with and deceleration by the surrounding medium on length scales below 1 kpc. Consequently, AGNs may be efficient at launching outflows $\sim10$ times more massive than previously uncovered by measurements of cold, obscuring material. We also examine the mass and velocity distribution of the cold gas, which ultimately forms a rotationally supported disk in simulated clusters. The rarity of such disks in observations indicates that further investigations are needed to understand the evolution of the cold gas after it forms., Comment: 23 pages, 12 figures, 1 table, accepted for publication in ApJ

Observational searches for dual active galactic nuclei (dAGNs) at kiloparsec separations are crucial for understanding the role of galaxy mergers in the evolution of galaxies. In addition, kpc-scale dAGNs may serve as the parent population of merging massive black hole (MBH) binaries, an important source of gravitational waves. We use a semi-analytical model to describe the orbital evolution of unequal mass MBH pairs under the influence of stellar and gaseous dynamical friction in post-merger galaxies. We quantify how the detectability of approximately 40,000 kpc-scale dAGNs depends on the structure of their host galaxies and the orbital properties of the MBH pair. Our models indicate that kpc-scale dAGNs are most likely to be detected in gas-rich post-merger galaxies with smaller stellar bulges and relatively massive, rapidly rotating gas disks. The detectability is also increased in systems with MBHs of comparable masses following low eccentricity prograde orbits. In contrast, dAGNs with retrograde, low eccentricity orbits are some of the least detectable systems among our models. The dAGNs in models in which the accreting MBHs are allowed to exhibit radiative feedback are characterized by a significantly lower overall detectability. The suppression in detectability is most pronounced in gas-rich merger remnant galaxies, where radiation feedback is more likely to arise. If so, then large, relatively gas poor galaxies may be the best candidates for detecting dAGNs., Comment: Accepted to ApJ, 16 pages, 13 figures

This white paper describes the science case for Very Long Baseline Interferometry (VLBI) and provides suggestions towards upgrade paths for the European VLBI Network (EVN). The EVN is a distributed long-baseline radio interferometric array, that operates at the very forefront of astronomical research. Recent results, together with the new science possibilities outlined in this vision document, demonstrate the EVN's potential to generate new and exciting results that will transform our view of the cosmos. Together with e-MERLIN, the EVN provides a range of baseline lengths that permit unique studies of faint radio sources to be made over a wide range of spatial scales. The science cases are reviewed in six chapters that cover the following broad areas: cosmology, galaxy formation and evolution, innermost regions of active galactic nuclei, explosive phenomena and transients, stars and stellar masers in the Milky Way, celestial reference frames and space applications. The document concludes with identifying the synergies with other radio, as well as multi-band/multi-messenger instruments, and provide the recommendations for future improvements. The appendices briefly describe other radio VLBI arrays, the technological framework for EVN developments, and a selection of spectral lines of astrophysical interest below 100 GHz. The document includes a glossary for non-specialists, and a list of acronyms at the end., Comment: Editors: Tiziana Venturi, Zsolt Paragi, Michael Lindqvist. "EVN Vision Document (2020)", 7 chapters, 3 appendices, 191 pages, 64 figures. Chapter coordinators and all the contributors are listed in the document

Dynamical friction (DF) against stars and gas is thought to be an important mechanism for orbital evolution of massive black holes (MBHs) in merger remnant galaxies. Recent theoretical investigations however show that DF does not always lead to MBH inspiral. For MBHs evolving in gas-rich backgrounds, the ionizing radiation that emerges from the innermost parts of their accretion flow can affect the surrounding gas in such a way to cause the MBHs to accelerate and gain orbital energy. This effect was dubbed "negative DF". We use a semi-analytic model to study the impact of negative DF on pairs of MBHs in merger remnant galaxies evolving under the combined influence of stellar and gaseous DF. Our results show that for a wide range of merger galaxy and MBH properties negative DF reduces the MBH pairing probability by $\sim 46\%$. The suppression of MBH pairing is most severe in galaxies with one or more of these properties: (1) a gas fraction of $f_g \geq 0.1$; (2) a galactic gas disk rotating close to the circular velocity; (3) MBH pairs in prograde, low eccentricity orbits, and (4) MBH pairs with mass $< 10^8\,$M$_\odot$. The last point is of importance because MBH pairs in this mass range are direct progenitors of merging binaries targeted by the future space-based gravitational wave observatory LISA., Comment: Accepted for publication to ApJ, 8 pages, 4 figures

We estimate the merger timescale of spectroscopically-selected, subparsec supermassive black hole binary (SMBHB) candidates by comparing their expected contribution to the gravitational wave background (GWB) with the sensitivity of current pulsar timing array (PTA) experiments and in particular, with the latest upper limit placed by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). We find that the average timescale to coalescence of such SMBHBs is $\langle t_{\rm evol} \rangle > 6\times 10^4\,$yr, assuming that their orbital evolution in the PTA frequency band is driven by emission of gravitational waves. If some fraction of SMBHBs do not reside in spectroscopically detected active galaxies, and their incidence in active and inactive galaxies is similar, then the merger timescale could be $\sim 10$ times longer, $\langle t_{\rm evol} \rangle > 6\times 10^5\,$yr. These limits are consistent with the range of timescales predicted by theoretical models and imply that all the SMBHB candidates in our spectroscopic sample could be binaries without violating the observational constraints on the GWB. This result illustrates the power of the multi-messenger approach, facilitated by the PTAs, in providing an independent statistical test of the nature of SMBHB candidates discovered in electromagnetic searches., Comment: Published in ApJL (9 pages, 3 figures)

Motivated by observational searches for massive black hole (MBH) pairs at kiloparsec separations we develop a semi-analytic model to describe their orbital evolution under the influence of stellar and gaseous dynamical friction (DF). The goal of this study is to determine how the properties of the merger remnant galaxy and the MBHs affect the likelihood and timescale for formation of a close MBH pair with separation of < 1 pc. We compute approximately 40,000 configurations that cover a wide range of host galaxy properties and investigate their impact on the orbital evolution of unequal mass MBH pairs. We find that the percentage for MBH pairing within a Hubble time is larger than 80% in remnant galaxies with a gas fraction < 20% and in galaxies hosting MBH pairs with total mass > 10^6 solar mass and mass ratios > 1/4. Among these, the remnant galaxies characterized by the fastest formation of close, gravitationally bound MBHs have one or more of the following properties: (1) large stellar bulge, (2) comparable mass MBHs and (3) a galactic gas disk rotating close to the circular speed. In such galaxies, the MBHs with the shortest inspiral times, which are likely progenitors of coalescing MBHs, are either on circular prograde orbits or on very eccentric retrograde orbits. Our model also indicates that remnant galaxies with opposite properties, that host slowly evolving MBH pairs, are the most likely hosts of dual AGNs at kiloparsec separations., Comment: Accepted to ApJ, 16 pages, 11 figures

Hyperaccretion occurs when the gas inflow rate onto a black hole (BH) is so high that the radiative feedback cannot reverse the accretion flow. This extreme process is a promising mechanism for the rapid growth of seed BHs in the early universe, which can explain high-redshift quasars powered by billion solar mass BHs. In theoretical models, spherical symmetry is commonly adopted for hyperaccretion flows; however, the sustainability of such structures on timescales corresponding to the BH growth has not been addressed yet. Here we show that stochastic interactions between the ionizing radiation from the BH and nonuniform accretion flow can lead to the formation of a rotating gas disk around the BH. Once the disk forms, the supply of gas to the BH preferentially occurs via biconical-dominated accretion flow perpendicular to the disk, avoiding the centrifugal barrier of the disk. Biconical-dominated accretion flows from opposite directions collide in the vicinity of the BH supplying high-density, low angular momentum gas to the BH, whereas most of the gas with nonnegligible angular momentum is deflected to the rotationally supported outflowing decretion disk. The disk becomes reinforced progressively as more mass from the biconical flow transfers to the disk and some of the outflowing gas from the disk is redirected to the biconical accretion funnels through a meridional structure. This axisymmetric hydrodynamic structure of a biconical-dominated accretion flow and decretion disk continues to provide uninterrupted flow of high-density gas to the BH., Comment: 14 pages, 10 figures, accepted for publication in ApJ

Galaxy clusters are the most massive collapsed structures in the universe whose potential wells are filled with hot, X-ray emitting intracluster medium. Observations however show that a significant number of clusters (the so-called cool-core clusters) also contain large amounts of cold gas in their centres, some of which is in the form of spatially extended filaments spanning scales of tens of kiloparsecs. These findings have raised questions about the origin of the cold gas, as well as its relationship with the central active galactic nucleus (AGN), whose feedback has been established as a ubiquitous feature in such galaxy clusters. Here we report a radiation hydrodynamic simulation of AGN feedback in a galaxy cluster, in which cold filaments form from the warm, AGN-driven outflows with temperatures between $10^4$ and $10^7$ K as they rise in the cluster core. Our analysis reveals a new mechanism, which, through the combination of radiative cooling and ram pressure, naturally promotes outflows whose cooling time is shorter than their rising time, giving birth to spatially extended cold gas filaments. Our results strongly suggest that the formation of cold gas and AGN feedback in galaxy clusters are inextricably linked and shed light on how AGN feedback couples to the intracluster medium., Comment: 24 pages, 8 figures. This is the authors' final accepted version including the supplementary information and references therein. The published version is available at http://doi.org/10.1038/s41550-020-1090-7

The quest for binary and dual supermassive black holes (SMBHs) at the dawn of the multi-messenger era is compelling. Detecting dual active galactic nuclei (AGN) -- active SMBHs at projected separations larger than several parsecs -- and binary AGN -- probing the scale where SMBHs are bound in a Keplerian binary -- is an observational challenge. The study of AGN pairs (either dual or binary) also represents an overarching theoretical problem in cosmology and astrophysics. The AGN triggering calls for detailed knowledge of the hydrodynamical conditions of gas in the imminent surroundings of the SMBHs and, at the same time, their duality calls for detailed knowledge on how galaxies assemble through major and minor mergers and grow fed by matter along the filaments of the cosmic web. This review describes the techniques used across the electromagnetic spectrum to detect dual and binary AGN candidates and proposes new avenues for their search. The current observational status is compared with the state-of-the-art numerical simulations and models for formation of dual and binary AGN. Binary SMBHs are among the loudest sources of gravitational waves (GWs) in the Universe. The search for a background of GWs at nHz frequencies from inspiralling SMBHs at low redshifts, and the direct detection of signals from their coalescence by the Laser Interferometer Space Antenna in the next decade, make this a theme of major interest for multi-messenger astrophysics. This review discusses the future facilities and observational strategies that are likely to significantly advance this fascinating field., Comment: 69 pages, 21 figures. Accepted for publication in New Astronomy Reviews

The past two decades have seen a tremendous investment in observational facilities that promise to reveal new and unprecedented discoveries about the universe. In comparison, the investment in theoretical work is completely dwarfed, even though theory plays a crucial role in the interpretation of these observations, predicting new types of phenomena, and informing observing strategies. In this white paper, we argue that in order to reach the promised critical breakthroughs in astrophysics over the next decade and well beyond, the national agencies must take a serious approach to investment in theoretical astrophysics research. We discuss the role of theory in shaping our understanding of the universe, and then we provide a multi-level strategy, from the grassroots to the national, to address the current underinvestment in theory relative to observational work., Comment: 7 pages, Submitted to the Astro2020 call for APC white papers

We present a method for comparing the H$\beta$ emission-line profiles of observed supermassive black hole (SBHB) candidates and models of sub-parsec SBHBs in circumbinary disks. Using the approach based on principal component analysis we infer the values of the binary parameters for the spectroscopic SBHB candidates and evaluate the parameter degeneracies, representative of the uncertainties intrinsic to such measurements. We find that as a population, the SBHB candidates favor the average value of the semimajor axis corresponding to $\log(a/M) \approx 4.20\pm 0.42$ and comparable mass ratios, $q>0.5$. If the SBHB candidates considered are true binaries, this result would suggest that there is a physical process that allows initially unequal mass systems to evolve toward comparable mass ratios (e.g., accretion that occurs preferentially onto the smaller of the black holes) or point to some, yet unspecified, selection bias. Our method also indicates that the SBHB candidates equally favor configurations in which the mini-disks are coplanar or misaligned with the binary orbital plane. If confirmed for true SBHBs, this finding would indicate the presence of a physical mechanism that maintains misalignment of the mini-disks down to sub-parsec binary separations (e.g., precession driven by gravitational torques). The probability distributions of the SBHB parameters inferred for the observed SBHB candidates and our control group of AGNs are statistically indistinguishable, implying that this method can in principle be used to interpret the observed emission-line profiles once a sample of confirmed SBHBs is available but cannot be used as a conclusive test of binarity., Comment: Accepted to ApJ (31 pages, 19 figures), includes referee's comments

Pulsar timing arrays (PTAs) are on the verge of detecting low-frequency gravitational waves (GWs) from supermassive black hole binaries (SMBHBs). With continued observations of a large sample of millisecond pulsars, PTAs will reach this major milestone within the next decade. Already, SMBHB candidates are being identified by electromagnetic surveys in ever-increasing numbers; upcoming surveys will enhance our ability to detect and verify candidates, and will be instrumental in identifying the host galaxies of GW sources. Multi-messenger (GW and electromagnetic) observations of SMBHBs will revolutionize our understanding of the co-evolution of SMBHs with their host galaxies, the dynamical interactions between binaries and their galactic environments, and the fundamental physics of accretion. Multi-messenger observations can also make SMBHBs 'standard sirens' for cosmological distance measurements out to $z\simeq0.5$. LIGO has already ushered in breakthrough insights in our knowledge of black holes. The multi-messenger detection of SMBHBs with PTAs will be a breakthrough in the years $2020-2030$ and beyond, and prepare us for LISA to help complete our views of black hole demographics and evolution at higher redshifts., Comment: Submitted to the Astro2020 decadal review. This is one of five core white papers written by members of the NANOGrav Collaboration, see also: J.Cordes et al., S.R.Taylor et al., X.Siemens et al., and E.Fonseca et al

Coalescing, massive black-hole (MBH) binaries are the most powerful sources of gravitational waves (GWs) in the Universe, which makes MBH science a prime focus for ongoing and upcoming GW observatories. The Laser Interferometer Space Antenna (LISA) -- a gigameter scale space-based GW observatory -- will grant us access to an immense cosmological volume, revealing MBHs merging when the first cosmic structures assembled in the Dark Ages. LISA will unveil the yet unknown origin of the first quasars, and detect the teeming population of MBHs of $10^4 - 10^7$ solar masses. forming within protogalactic halos. The Pulsar Timing Array, a galactic-scale GW survey, can access the largest MBHs the Universe, detecting the cosmic GW foreground from inspiraling MBH binaries of about 10^9 solar masses. LISA can measure MBH spins and masses with precision far exceeding that from electromagnetic (EM) probes, and together, both GW observatories will provide the first full census of binary MBHs, and their orbital dynamics, across cosmic time. Detecting the loud gravitational signal of these MBH binaries will also trigger alerts for EM counterpart searches, from decades (PTAs) to hours (LISA) prior to the final merger. By witnessing both the GW and EM signals of MBH mergers, precious information will be gathered about the rich and complex environment in the aftermath of a galaxy collision. The unique GW characterization of MBHs will shed light on the deep link between MBHs of $10^4-10^{10}$ solar masses and the grand design of galaxy assembly, as well as on the complex dynamics that drive MBHs to coalescence., Comment: 5 pages, 2 figures submitted to the Astro2020 Decadal

The tidal disruption of stars by (super-)massive black holes in galactic nuclei has been discussed in theoretical terms for about 30 years but only in the past decade have we been able to detect such events in substantial numbers. Thus, we are now starting to carry out observational tests of models for the disruption. We are also formulating expectations for the inspiral and disruption of white dwarfs by intermediate-mass black holes with masses $< 10^5\;{\rm M}_\odot$. Such events are very rich with information and open a new window to intermediate-mass black holes, thought to live in dwarf galaxies and star clusters. They can inform us of the demographics of intermediate-mass black holes, stellar populations and dynamics in their immediate vicinity, and the physics of accretion of hydrogen-deficient material. The combination of upcoming transient surveys using ground-based, electromagnetic observatories and low-frequency gravitational wave observations is ideal for exploiting tidal disruptions of white dwarfs. The detection rate of gravitational wave signals, optimistically, may reach a few dozen per year in a volume up to $z\approx 0.1$. Gravitational wave observations are particularly useful because they yield the masses of the objects involved and allow determination of the spin of the black hole, affording tests of physical models for black hole formation and growth. They also give us advance warning of the electromagnetic flares by weeks or more. The right computing infrastructure for modern models for the disruption process and event rates will allow us to make the most of the upcoming observing facilities., Comment: White Paper no. 16, submitted on Feb 18, 2019 to Astro2020 (2020 Decadal Survey on Astronomy and Astrophyscs)

Recent observations of giant ellipticals and brightest cluster galaxies (BCGs) provide tentative evidence for a correlation between the luminosity of the H$\alpha$ emitting gas filaments and the strength of feedback associated with the active galactic nucleus (AGN). Motivated by this, we use 3D radiation-hydrodynamic simulations with the code Enzo to examine and quantify the relationship between the observable properties of the H$\alpha$ filaments and the kinetic and radiative feedback from supermassive black holes in BCGs. We find that the spatial extent and total mass of the filaments show positive correlations with AGN feedback power and can therefore be used as probes of the AGN activity. We also examine the relationship between the AGN feedback power and velocity dispersion of the H$\alpha$ filaments and find that the kinetic luminosity shows statistically significant correlation with the component of the velocity dispersion along the jet axis but not the components perpendicular to it., Comment: Accepted to ApJL (7 pages, 4 figures). Revised version, includes referee's comments

Recent observations provide evidence that some cool-core clusters (CCCs) host quasars in their brightest cluster galaxies (BCGs). Motivated by these findings we use 3D radiation-hydrodynamic simulations with the code Enzo to explore the joint role of the kinetic and radiative feedback from supermassive black holes (SMBHs) in BCGs. We implement kinetic feedback as sub-relativistic plasma outflows and model radiative feedback using the ray-tracing radiative transfer or thermal energy injection. In our simulations the central SMBH transitions between the radiatively efficient and radiatively inefficient states on timescales of a few Gyr, as a function of its accretion rate. The timescale for this transition depends primarily on the fraction of power allocated to each feedback mode, and to a lesser degree on the overall feedback luminosity of the active galactic nucleus (AGN). Specifically, we find that (a) kinetic feedback must be present at both low and high accretion rates in order to prevent the cooling catastrophe, and (b) its contribution likely accounts for > 10% of the total AGN feedback power, since below this threshold simulated BCGs tend to host radio-loud quasars most of the time, in apparent contrast with observations. We also find a positive correlation between the AGN feedback power and the mass of the cold gas filaments in the cluster core, indicating that observations of H$\alpha$ filaments can be used as a measure of AGN feedback., Comment: Published in ApJ (28 pages, 14 figures)

The Next-Generation Very Large Array (ngVLA) has the potential to be a workhorse for the discovery and study of paired supermassive black holes either at large separations (dual) or in tightly bound systems (binary). In this chapter, we outline the science case for the study of these supermassive pairs, and summarize discovery methods that can be used at radio wavelengths to discover them: including morphological, spectral, and time-domain identifications. One critical aspect of this work is that multi-messenger binary black hole studies may be possible with the ngVLA when combined with gravitational-wave searches using pulsar timing array techniques. However, long-baseline interferometery (>>1000 km) will make this possibility more likely by expanding the redshift range at which radio emission arising from two separate black holes may be resolved and studied., Comment: 12 pages, 3 figures; Chapter in the ngVLA science book

We present an improved semi-analytic model for calculation of the broad optical emission-line signatures from sub-parsec supermassive black hole binaries (SBHBs) in circumbinary disks. The second-generation model improves upon the treatment of radiative transfer by taking into account the effect of the radiation driven accretion disk wind on the properties of the emission-line profiles. Analysis of 42.5 million modeled emission-line profiles shows that correlations between the profile properties and SBHB parameters identified in the first-generation model are preserved, indicating that their diagnostic power is not diminished. The profile shapes are a more sensitive measure of the binary orbital separation and the degree of alignment of the black hole mini-disks, and are less sensitive to the SBHB mass ratio and orbital eccentricity. We also find that modeled profile shapes are more compatible with the observed sample of SBHB candidates than with our control sample of regular AGNs. Furthermore, if the observed sample of SBHBs is made up of genuine binaries, it must include compact systems with comparable masses, and misaligned mini-disks. We note that the model described in this paper can be used to interpret the observed emission-line profiles once a sample of confirmed SBHBs is available but cannot be used to prove that the observed SBHB candidates are true binaries., Comment: Accepted to ApJ (25 pages, 16 figures), includes referee's comments

Motivated by observational searches for sub-parsec supermassive black hole binaries (SBHBs) we develop a modular analytic model to determine the likelihood for detection of SBHBs by ongoing spectroscopic surveys. The model combines the parametrized rate of orbital evolution of SBHBs in circumbinary disks with the selection effects of spectroscopic surveys and returns a multivariate likelihood for SBHB detection. Based on this model we find that in order to evolve into the detection window of the spectroscopic searches from larger separations in less than a Hubble time, $10^8M_\odot$ SBHBs must, on average, experience angular momentum transport faster than that provided by a disk with accretion rate $0.06\,\dot{M}_E$. Spectroscopic searches with yearly cadence of observations are in principle sensitive to binaries with orbital separations $< {\rm few}\times 10^4\, r_g$ ($r_g = GM/c^2$ and $M$ is the binary mass), and for every one SBHB in this range there should be over 200 more gravitationally bound systems with similar properties, at larger separations. Furthermore, if spectra of all SBHBs in this separation range exhibit the AGN-like emission lines utilized by spectroscopic searches, the projection factors imply five undetected binaries for each observed $10^8M_\odot$ SBHB with mass ratio $0.3$ and orbital separation $10^4\,r_g$ (and more if some fraction of SBHBs is inactive). This model can be used to infer the most likely orbital parameters for observed SBHB candidates and to provide constraints on the rate of orbital evolution of SBHBs, if observed candidates are shown to be genuine binaries., Comment: Accepted to ApJ (16 pages, 11 figures). Revised version includes referee's comments

Accretion of gas and interaction of matter and radiation are at the heart of many questions pertaining to black hole (BH) growth and coevolution of massive BHs and their host galaxies. To answer them it is critical to quantify how the ionizing radiation that emanates from the innermost regions of the BH accretion flow couples to the surrounding medium and how it regulates the BH fueling. In this work we use high resolution 3-dimensional (3D) radiation-hydrodynamic simulations with the code Enzo, equipped with adaptive ray tracing module Moray, to investigate radiation-regulated BH accretion of cold gas. Our simulations reproduce findings from an earlier generation of 1D/2D simulations: the accretion powered UV and X-ray radiation forms a highly ionized bubble, which leads to suppression of BH accretion rate characterized by quasi-periodic outbursts. A new feature revealed by the 3D simulations is the highly turbulent nature of the gas flow in vicinity of the ionization front. During quiescent periods between accretion outbursts, the ionized bubble shrinks in size and the gas density that precedes the ionization front increases. Consequently, the 3D simulations show oscillations in the accretion rate of only ~2-3 orders of magnitude, significantly smaller than 1D/2D models. We calculate the energy budget of the gas flow and find that turbulence is the main contributor to the kinetic energy of the gas but corresponds to less than 10% of its thermal energy and thus does not contribute significantly to the pressure support of the gas., Comment: 13 pages, 11 figures, accepted for publication in ApJ

We have been spectroscopically monitoring 88 quasars selected to have broad H$\beta$ emission lines offset from their systemic redshift by thousands of km s$^{-1}$. By analogy with single-lined spectroscopic binary stars, we consider these quasars to be candidates for hosting supermassive black hole binaries (SBHBs). In this work we present new radial velocity measurements, typically 3-4 per object over a time period of up to 12 years in the observer's frame. In 29/88 of the SBHB candidates no variability of the shape of the broad H$\beta$ profile is observed, which allows us to make reliable measurements of radial velocity changes. Among these, we identify three objects that have displayed systematic and monotonic velocity changes by several hundred km s$^{-1}$ and are prime targets for further monitoring. Because the periods of the hypothetical binaries are expected to be long, we cannot hope to observe many orbital cycles during our lifetimes. Instead, we seek to evaluate the credentials of the SBHB candidates by attempting to rule out the SBHB hypothesis. In this spirit, we present a method for placing a lower limit on the period, and thus the mass, of the SBHBs under the assumption that the velocity changes we observe are due to orbital motion. Given the duration of our monitoring campaign and the uncertainties in the radial velocities, we were able to place a lower limit on the total mass in the range $4.7\times10^4-3.8\times10^8$ $M_{\scriptscriptstyle \odot}$, which does not yet allow us to rule out the SBHB hypothesis for any candidates., Comment: Accepted for publication in MNRAS, full version of Table 1 will be available in MNRAS online version of the paper

Dynamical friction is thought to be a principal mechanism responsible for orbital evolution of massive black holes (MBHs) in the aftermath of galactic mergers and an important channel for formation of gravitationally bound MBH binaries. We use 2D radiative hydrodynamic simulations to investigate the efficiency of dynamical friction in the presence of radiative feedback from an MBH moving through a uniform density gas. We find that ionizing radiation that emerges from the innermost parts of the MBH's accretion flow strongly affects the dynamical friction wake and renders dynamical friction inefficient for a range of physical scenarios. MBHs in this regime tend to experience positive net acceleration, meaning that they speed up, contrary to the expectations for gaseous dynamical friction in absence of radiative feedback. The magnitude of this acceleration is however negligibly small and should not significantly alter the velocity of MBHs over relevant physical timescales. Our results suggest that suppression of dynamical friction is more severe at the lower mass end of the MBH spectrum which, compounded with inefficiency of the gas drag for lower mass objects in general, implies that $< 10^7$ solar mass MBHs have fewer means to reach the centers of merged galaxies. These findings provide formulation for a sub-resolution model of dynamical friction in presence of MBH radiative feedback that can be easily implemented in large scale simulations., Comment: 8 pages of published article in ApJ, 4 pages of Erratum attached

Motivated by advances in observational searches for sub-parsec supermassive black hole binaries (SBHBs) made in the past few years we develop a semi-analytic model to describe spectral emission line signatures of these systems. The goal of this study is to aid the interpretation of spectroscopic searches for binaries and help test one of the leading models of binary accretion flows in the literature: SBHB in a circumbinary disk. In this work we present the methodology and a comparison of the preliminary model with the data. We model SBHB accretion flows as a set of three accretion disks: two mini-disks that are gravitationally bound to the individual black holes and a circumbinary disk. Given a physically motivated parameter space occupied by sub-parsec SBHBs, we calculate a synthetic database of nearly 15 million broad optical emission line profiles and explore the dependence of the profile shapes on characteristic properties of SBHBs. We find that the modeled profiles show distinct statistical properties as a function of the semi-major axis, mass ratio, eccentricity of the binary, and the degree of alignment of the triple disk system. This suggests that the broad emission line profiles from SBHB systems can in principle be used to infer the distribution of these parameters and as such merit further investigation. Calculated profiles are more morphologically heterogeneous than the broad emission lines in observed SBHB candidates and we discuss improved treatment of radiative transfer effects which will allow direct statistical comparison of the two groups., Comment: 27 pages, 20 figures, accepted to ApJ

Observations have revealed a relative paucity of red giant (RG) stars within the central 0.5pc in the Galactic Center (GC). Motivated by this finding we investigate the hypothesis that collisions of stars with a fragmenting accretion disk are responsible for the observed dearth of evolved stars. We use 3-dimensional hydrodynamic simulations to model a star with radius $10 R_{\odot}$ and mass $1 M_{\odot}$, representative of the missing population of RGs, colliding with high density clumps. We find that multiple collisions with clumps of column density $\gtrsim10^{8}\, {\rm g\,cm^{-2}}$ can strip a substantial fraction of the star's envelope and in principle render it invisible to observations. Simulations confirm that repeated impacts are particularly efficient in driving mass loss as partially stripped RGs expand and have increased cross sections for subsequent collisions. Because the envelope is unbound on account of the kinetic energy of the star, any significant amount of stripping of the RG population in the GC should be mirrored by a systematic decay of their orbits and possibly by their enhanced rotational velocity. To be viable, this scenario requires that the total mass of the fragmenting disk has been several orders of magnitude higher than that of the early type stars which now form the stellar disk in the GC., Comment: 15 pages, 11 figures, accepted for publication in ApJ

We examine radiation-regulated accretion onto intermediate-mass and massive black holes (BHs) embedded in a bulge component. Using spherically symmetric one-dimensional radiation-hydrodynamics simulations, we track the growth of BHs accreting from a cold, neutral gas reservoir with temperature T=10^4 K. We find that the accretion rate of BHs embedded in bulges is proportional to r_{B,eff}/r_B, where r_{B,eff} is the increased effective Bondi radius that includes the gravitational potential of the bulge, and r_B is the Bondi radius of the BH. The radiative feedback from the BH suppresses the cold accretion rate to ~1 percent of the Bondi rate when a bulge is not considered. However, we find that the BH fueling rate increases rapidly when the bulge mass M_bulge is greater than the critical value of 10^6 M_sun and is proportional to M_bulge. Since the critical bulge mass is independent of the central BH mass M_{BH}, the growth rate of BHs with masses of 10^2, 10^4, and 10^6 M_sun exhibits distinct dependencies on the bulge-to-BH mass ratio. Our results imply that light seed BHs (<= 10^2 M_sun) which might be the remnants of the Pop III stars, cannot grow through accretion coevally with the early assembly of the bulge of the host galaxies until the bulge reaches the critical mass. However, massive BH seeds (>= 10^5 M_sun) that may form via direct collapse, are more likely to be embedded in a supercritical bulge and thus can grow efficiently coupling to the host galaxies and driving the early evolution of the M_{BH}-$\sigma$ relationship., Comment: 10 pages, 10 Figures, 1 Table, submitted to ApJ

The dynamics of massive black holes (BHs) in galaxy mergers is a rich field of research that has seen much progress in recent years. In this contribution we briefly review the processes describing the journey of BHs during mergers, from the cosmic context all the way to when BHs coalesce. If two galaxies each hosting a central BH merge, the BHs would be dragged towards the center of the newly formed galaxy. If/when the holes get sufficiently close, they coalesce via the emission of gravitational waves. How often two BHs are involved in galaxy mergers depends crucially on how many galaxies host BHs and on the galaxy merger history. It is therefore necessary to start with full cosmological models including BH physics and a careful dynamical treatment. After galaxies have merged, however, the BHs still have a long journey until they touch and coalesce. Their dynamical evolution is radically different in gas-rich and gas-poor galaxies, leading to a sort of "dichotomy" between high-redshift and low-redshift galaxies, and late-type and early-type, typically more massive galaxies., Comment: Proceedings of FM14, The Gravitational Wave Symphony of Structure Formation, XXIXth IAU General Assembly, to be published in Astronomy in Focus, Volume 1

We present new spectroscopic observations that are part of our continuing monitoring campaign of 88 quasars at z<0.7 whose broad H$\beta$ lines are offset from their systemic redshifts by a few thousand km/s. These quasars have been considered candidates for hosting supermassive black hole binaries (SBHBs) by analogy with single-lined spectroscopic binary stars. We present the data and describe our improved analysis techniques, which include an extensive evaluation of uncertainties. We also present a variety of measurements from the spectra that are of general interest and will be useful in later stages of our analysis. Additionally, we take this opportunity to study the variability of the optical continuum and integrated flux of the broad H$\beta$ line. We compare the variability properties of the SBHB candidates to those of a sample of typical quasars with similar redshifts and luminosities observed multiple times during the Sloan Digital Sky Survey. We find that the variability properties of the two samples are similar (variability amplitudes of 10-30% on time scales of approximately 1-7 years) and that their structure functions can be described by a common model with parameters characteristic of typical quasars. These results suggest that the broad-line regions of SBHB candidates have a similar extent as those of typical quasars. We discuss the implications of this result for the SBHB scenario and ensuing constraints on the orbital parameters., Comment: 19 pages in ApJ format. Accepted for publication in the Astrophysical Journal Supplements

We present a multi-wavelength study of a newly discovered compact group (CG), SDSS J0959+1259, based data from XMM-Newton, SDSS and the Calar Alto optical imager BUSCA. With a maximum velocity offset of 500 km s$^{-1}$, a mean redshift of 0.035, and a mean spatial extension of 480 kpc, this CG is exceptional in having the highest concentration of nuclear activity in the local Universe, established with a sensitivity limit L$_{X}>4\times $10$^{40}$ erg s$^{-1}$ in 2--10 keV band and R-band magnitude $M_R < -19$. The group is composed of two type-2 Seyferts, one type-1 Seyfert, two LINERs and three star forming galaxies. Given the high X-ray luminosity of LINERs which reaches $\sim 10^{41}$ erg s$^{-1}$, it is likely that they are also accretion driven, bringing the number of active nuclei in this group to to 5 out of 8 (AGN fraction of 60\%). The distorted shape of one member of the CG suggests that strong interactions are taking place among its galaxies through tidal forces. Therefore, this system represents a case study for physical mechanisms that trigger nuclear activity and star formation in CGs., Comment: 8 pages, 4 figures and 4 tables. Accepted for publication in MNRAS

Motivated by an improved multi-wavelength observational coverage of the transient sky, we investigate the importance of relativistic effects in disruptions of stars by non-spinning black holes (BHs). This paper focuses on calculating the ballistic rate of return of debris to the black hole as this rate is commonly assumed to be proportional to the light curve of the event. We simulate the disruption of a low mass main sequence star by BHs of varying masses ($10^5,10^6,10^7 M_\odot$) and of a white dwarf by a $10^5 M_\odot$ BH. Based on the orbital energy as well as angular momentum of the debris, we infer the orbital distribution and estimate the return rate of the debris following the disruption. We find two signatures of relativistic disruptions: a gradual rise as well as a delayed peak in the return rate curves relative to their Newtonian analogs. Assuming that the return rates are proportional to the light curves, we find that relativistic effects are in principle measurable given the cadence and sensitivity of the current transient sky surveys. Accordingly, using a simple model of a relativistic encounter with a Newtonian parametric fit of the peak leads to an overestimate in the BH mass by a factor of $\sim {\rm few}\times0.1$ and $\sim {\rm few}$ in the case of the main sequence star and white dwarf tidal disruptions, respectively., Comment: Accepted by Phys. Rev. D. Some changes in response to referee reports

Gravitationally bound supermassive black hole binaries (SBHBs) are thought to be a natural product of galactic mergers and growth of the large scale structure in the universe. They however remain observationally elusive, thus raising a question about characteristic observational signatures associated with these systems. In this conference proceeding I discuss current theoretical understanding and latest advances and prospects in observational searches for SBHBs., Comment: 17 pages, 4 figures. To appear in the Proceedings of 2014 Sant Cugat Forum on Astrophysics. Astrophysics and Space Science Proceedings, ed. C.Sopuerta (Berlin: Springer-Verlag)

The development of a new generation of theoretical models for tidal disruptions is timely, as increasingly diverse events are being captured in surveys of the transient sky. Recently, Gezari et al. reported a discovery of a new class of tidal disruption events: the disruption of a helium-rich stellar core, thought to be a remnant of a red giant (RG) star. Motivated by this discovery and in anticipation of others, we consider tidal interaction of an RG star with a supermassive black hole (SMBH) which leads to the stripping of the stellar envelope and subsequent inspiral of the compact core toward the black hole. Once the stellar envelope is removed the inspiral of the core is driven by tidal heating as well as the emission of gravitational radiation until the core either falls into the SMBH or is tidally disrupted. In the case of tidal disruption candidate PS1-10jh we find that there is a set of orbital solutions at high eccentricities in which the tidally stripped hydrogen envelope is accreted by the SMBH before the helium core is disrupted. This places the RG core in a portion of parameter space where strong tidal heating can lift the degeneracy of the compact remnant and disrupt it before it reaches the tidal radius. We consider how this sequence of events explains the puzzling absence of the hydrogen emission lines from the spectrum of PS1-10jh and gives rise to its other observational features., Comment: Revised version accepted for publication in ApJ, 788, 99 (12 pages, 2 figures)

The role played by magnetic fields in the intracluster medium (ICM) of galaxy clusters is complex. The weakly collisional nature of the ICM leads to thermal conduction that is channelled along field lines. This anisotropic heat conduction profoundly changes the stability of the ICM atmosphere, with convective stabilities being driven by temperature gradients of either sign. Here, we employ the Athena magnetohydrodynamic code to investigate the local non-linear behavior of the heat-flux driven buoyancy instability (HBI), relevant in the cores of cooling-core clusters where the temperature increases with radius. We study a grid of 2-d simulations that span a large range of initial magnetic field strengths and numerical resolutions. For very weak initial fields, we recover the previously known result that the HBI wraps the field in the horizontal direction thereby shutting off the heat flux. However, we find that simulations which begin with intermediate initial field strengths have a qualitatively different behavior, forming HBI-stable filaments that resist field-line wrapping and enable sustained vertical conductive heat flux at a level of 10--25% of the Spitzer value. While astrophysical conclusions regarding the role of conduction in cooling cores require detailed global models, our local study proves that systems dominated by HBI do not necessarily quench the conductive heat flux., Comment: 11 pages, 9 figures, Accepted for Publication in ApJ

Observations of low-metallicity halo stars have revealed a puzzling result: the abundance of \li7 in these stars is at least three times lower than their predicted primordial abundance. It is unclear whether the cause of this disagreement is a lack of understanding of lithium destruction mechanisms in stars or the non-standard physics behind the Big Bang Nucleosynthesis (BBN). Uncertainties related to the destruction of lithium in stars can be circumvented if lithium abundance is measured in the "pristine" gas of the low metallicity systems. The first measurement in one such system, the Small Magellanic Cloud (SMC), was found to be at the level of the pure expected primordial value, but is on the other hand, just barely consistent with the expected galactic abundance for the system at the SMC metallicity, where important lithium quantity was also produced in interactions of galactic cosmic rays (GCRs) and presents an addition to the already present primordial abundance. Due to the importance of the SMC lithium measurement for the resolution of the lithium problem, we here draw attention to the possibility of another post-BBN production channel of lithium, which could present an important addition to the observed SMC lithium abundance. Besides standard galactic cosmic rays, additional post-BBN production of lithium might come from cosmic rays accelerated in galaxy-galaxy interactions. This might be important for a system such is the SMC, which has experienced galaxy harassment in its history...(abridged), Comment: 9 pages. Accepted to PRD. Version update for some revisions