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 40\%$. 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.2$; (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 $< {\rm few}\times 10^7\,$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.

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 (PTAs) 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 pulsar timing arrays, in providing an independent statistical test of the nature of SMBHB candidates discovered in electromagnetic searches.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

We report on an optical (SDSS) and X-ray (XMM) study of an optically selected sample of four dual AGN systems at projected separations of 30--60~kpc. All sources are detected in the X-ray band (0.3-10 keV); seven objects are optically identified as Seyfert, while one source, optically classified as a LINER, is likely powered by accretion in virtue of its relatively high X-ray luminosity (1.2$\times10^{41}$ erg/s). Six of the eight objects are obscured in X-rays with N$_{\rm H} \geq$ 10$^{23}$ cm$^{-2}$; three of these, whose X-ray spectrum is dominated by a reflection component, are likely Compton-thick (N$_{\rm H} \geq$ 10$^{24}$ cm$^{-2}$). This finding is in agreement with the hypothesis that galaxy encounters are effective in driving gas inflow toward the nuclear region, thus increasing the obscuration. We compare the absorption properties in our dual AGN with those in larger samples observed in X-rays but selected in different ways (optical, IR and hard X-rays). We find that the obscured (N$_{\rm H} \geq$ 10$^{22}$ cm$^{-2}$) AGN fraction within the larger sample is 84$\pm$4 per cent (taking into account the 90 per cent error on the N$_{\rm H}$ measure) up to large pair separations ($\sim$100~kpc). This is statistically higher than the fraction of obscured AGN in isolated galaxies found in X-ray surveys.

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.

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.

In this paper we present an overview of the MAGNA (Multiple AGN Activity) project aiming at a comprehensive study of multiple supemassive black hole systems. With the main goal to characterize the sources in merging systems at different stages of evolution, we selected a sample of objects optically classified as multiple systems on the basis of emission line diagnostics and started a massive multiband observational campaign. Here we report on the discovery of the exceptionally high AGN density compact group SDSS~J0959+1259. A multiband study 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. We also present a preliminary analysis of the multiple AGN system SDSS~J1038+3921.}

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.

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.

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.

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.

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.

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.

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.

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.

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)

Observations of the Galactic Center (GC) have accumulated a multitude of "forensic" evidence indicating that several million years ago the center of the Milky Way galaxy was teaming with starforming and accretion-powered activity -- this paints a rather different picture from the GC as we understand it today. We examine a possibility that this epoch of activity could have been triggered by the infall of a satellite galaxy into the Milky Way which began at the redshift of 10 and ended few million years ago with a merger of the Galactic supermassive black hole with an intermediate mass black hole brought in by the inspiralling satellite.

The intracluster medium of galaxy clusters is a weakly collisional, high-beta plasma in which the transport of heat and momentum occurs primarily along magnetic-field lines. Anisotropic heat conduction allows convective instabilities to be driven by temperature gradients of either sign, the magnetothermal instability (MTI) in the outskirts of non-isothermal clusters and the heat-flux buoyancy-driven instability (HBI) in their cooling cores. We employ the Athena MHD code to investigate the nonlinear evolution of these instabilities, self-consistently including the effects of anisotropic viscosity (i.e. Braginskii pressure anisotropy), anisotropic conduction, and radiative cooling. We highlight the importance of the microscale instabilities that inevitably accompany and regulate the pressure anisotropies generated by the HBI and MTI. We find that, in all but the innermost regions of cool-core clusters, anisotropic viscosity significantly impairs the ability of the HBI to reorient magnetic-field lines orthogonal to the temperature gradient. Thus, while radio-mode feedback appears necessary in the central few tens of kpc, conduction may be capable of offsetting radiative losses throughout most of a cool core over a significant fraction of the Hubble time. Magnetically-aligned cold filaments are then able to form by local thermal instability. Viscous dissipation during the formation of a cold filament produces accompanying hot filaments, which can be searched for in deep Chandra observations of nearby cool-core clusters. In the case of the MTI, anisotropic viscosity maintains the coherence of magnetic-field lines over larger distances than in the inviscid case, providing a natural lower limit for the scale on which the field can fluctuate freely. In the nonlinear state, the magnetic field exhibits a folded structure in which the field-line curvature and field strength are anti-correlated.

Modeling the late inspiral and merger of supermassive black holes is central to understanding accretion processes and the conditions under which electromagnetic emission accompanies gravitational waves. We use fully general relativistic, hydrodynamics simulations to investigate how electromagnetic signatures correlate with black hole spins, mass ratios, and the gaseous environment in this final phase of binary evolution. In all scenarios, we find some form of characteristic electromagnetic variability whose pattern depends on the spins and binary mass ratios. Binaries in hot accretion flows exhibit a flare followed by a sudden drop in luminosity associated with the plunge and merger, as well as quasi-periodic oscillations correlated with the gravitational waves during the inspiral. Conversely, circumbinary disk systems are characterized by a low luminosity of variable emission, suggesting challenging prospects for their detection.

It has recently been suggested that conduction-driven magnetohydrodynamic (MHD) instabilities may operate at all radii within an intracluster medium (ICM), and profoundly affect the structure of a cluster's magnetic field. Where MHD instabilities dominate the dynamics of an ICM, they will re-orient magnetic field lines perpendicular to the temperature gradient inside a cooling core, or parallel to the temperature gradient outside it. This characteristic structure of magnetic field could be probed by measurements of polarized radio emission from background sources. Motivated by this possibility we have constructed 3-d models of a magnetized cooling core cluster and calculated Faraday rotation measure (RM) maps in the plane of the sky under realistic observing conditions. We compare a scenario in which magnetic field geometry is characterized by conduction driven MHD instabilities to that where it is determined by isotropic turbulent motions. We find that future high-sensitivity spectro-polarimetric measurements of RM, such as will be enabled by the Expanded Very Large Array and Square Kilometer Array can distinguish between these two cases with plausible exposure times. Such observations will test the existence of conduction-driven MHD instabilities in dynamically relaxed cooling core clusters. More generally, our findings imply that observations of Faraday RM should be able to discern physical mechanisms that result in qualitatively different magnetic field topologies, without a priori knowledge about the nature of the processes.

What are the properties of accretion flows in the vicinity of coalescing supermassive black holes (SBHs)? The answer to this question has direct implications for the feasibility of coincident detections of electromagnetic (EM) and gravitational wave (GW) signals from coalescences. Such detections are considered to be the next observational grand challenge that will enable testing general relativity in the strong, nonlinear regime and improve our understanding of evolution and growth of these massive compact objects. In this paper we review the properties of the environment of coalescing binaries in the context of the circumbinary disk and hot, radiatively inefficient accretion flow models and use them to mark the extent of the parameter space spanned by this problem. We report the results from an ongoing, general relativistic, hydrodynamical study of the inspiral and merger of black holes, motivated by the latter scenario. We find that correlated EM+GW oscillations can arise during the inspiral phase followed by the gradual rise and subsequent drop-off in the light curve at the time of coalescence. While there are indications that the latter EM signature is a more robust one, a detection of either signal coincidentally with GWs would be a convincing evidence for an impending SBH binary coalescence. The observability of an EM counterpart in the hot accretion flow scenario depends on the details of a model. In the case of the most massive binaries observable by the Laser Interferometer Space Antenna, upper limits on luminosity imply that they may be identified by EM searches out to z~0.1-1. However, given the radiatively inefficient nature of the gas flow, we speculate that a majority of massive binaries may appear as low luminosity AGN in the local universe.

Coincident detections of electromagnetic (EM) and gravitational wave (GW) signatures from coalescence events of supermassive black holes are the next observational grand challenge. Such detections will provide the means to study cosmological evolution and accretion processes associated with these gargantuan compact objects. More generally, the observations will enable testing general relativity in the strong, nonlinear regime and will provide independent cosmological measurements to high precision. Understanding the conditions under which coincidences of EM and GW signatures arise during supermassive black hole mergers is therefore of paramount importance. As an essential step towards this goal, we present results from the first fully general relativistic, hydrodynamical study of the late inspiral and merger of equal-mass, spinning supermassive black hole binaries in a gas cloud. We find that variable EM signatures correlated with GWs can arise in merging systems as a consequence of shocks and accretion combined with the effect of relativistic beaming. The most striking EM variability is observed for systems where spins are aligned with the orbital axis and where orbiting black holes form a stable set of density wakes, but all systems exhibit some characteristic signatures that can be utilized in searches for EM counterparts. In the case of the most massive binaries observable by the Laser Interferometer Space Antenna, calculated luminosities imply that they may be identified by EM searches to z = 1, while lower mass systems and binaries immersed in low density ambient gas can only be detected in the local universe.

Detection of electromagnetic (EM) counterparts of pre-coalescence binaries has very important implications for our understanding of the evolution of these systems as well as the associated accretion physics. In addition, a combination of EM and gravitational wave signatures observed from coalescing supermassive black hole binaries (SBHBs) would provide independent measurements of redshift and luminosity distance, thus allowing for high precision cosmological measurements. However, a statistically significant sample of these objects is yet to be attained and finding them observationally has proven to be a difficult task. Here we discuss existing observational evidence and how further advancements in the theoretical understanding of observational signatures of SBHBs before and after the coalescence can help in future searches.

We perform a suite of simulations of cooling cores in clusters of galaxies in order to investigate the effect of the recently discovered heat flux buoyancy instability (HBI) on the evolution of cores. Our models follow the 3-dimensional magnetohydrodynamics (MHD) of cooling cluster cores and capture the effects of anisotropic heat conduction along the lines of magnetic field, but do not account for the cosmological setting of clusters or the presence of AGN. Our model clusters can be divided into three groups according to their final thermodynamical state: catastrophically collapsing cores, isothermal cores, and an intermediate group whose final state is determined by the initial configuration of magnetic field. Modeled cores that are reminiscent of real cluster cores show evolution towards thermal collapse on a time scale which is prolonged by a factor of ~2-10 compared with the zero-conduction cases. The principal effect of the HBI is to re-orient field lines to be perpendicular to the temperature gradient. Once the field has been wrapped up onto spherical surfaces surrounding the core, the core is insulated from further conductive heating (with the effective thermal conduction suppressed to less than 1/100th of the Spitzer value) and proceeds to collapse. We speculate that, in real clusters, the central AGN and possibly mergers play the role of "stirrers," periodically disrupting the azimuthal field structure and allowing thermal conduction to sporadically heat the core.

The association of an electromagnetic signal with the merger of a pair of supermassive black holes would have many important implications. For example, it would provide new information about gas and magnetic field interactions in dynamical spacetimes as well as a combination of redshift and luminosity distance that would enable precise cosmological tests. A proposal first made by Bode & Phinney (2007) is that because radiation of gravitational waves during the final inspiral and merger of the holes is abrupt and decreases the mass of the central object by a few percent, there will be waves in the disk that can steepen into shocks and thus increase the disk luminosity in a characteristic way. We evaluate this process analytically and numerically. We find that shocks only occur when the fractional mass loss exceeds the half-thickness (h/r) of the disk, hence significant energy release only occurs for geometrically thin disks which are thus at low Eddington ratios. This strongly limits the effective energy release, and in fact our simulations show that the natural variations in disk luminosity are likely to obscure this effect entirely. However, we demonstrate that the reduction of luminosity caused by the retreat of the inner edge of the disk following mass loss is potentially detectable. This decrease occurs even if the disk is geometrically thick, and lasts for a duration on the order of the viscous time of the modified disk. Observationally, the best prospect for detection would be a sensitive future X-ray instrument with a field of view of on the order of a square degree, or possibly a wide-field radio array such as the Square Kilometer Array, if the disk changes produce or interrupt radio emission from a jet.

A search for recoiling supermassive black hole candidates recently yielded the best candidate thus far, SDSS J092712.65+294344.0 reported by Komossa et al. Here we propose the alternative hypothesis that this object is a supermassive black hole binary. From the velocity shift imprinted in the emission-line spectrum we infer an orbital period of ~190 years for a binary mass ratio of 0.1, a secondary black hole mass of 100 million solar masses, and assuming inclination and orbital phase angles of 45 degrees. In this model the origin of the blueshifted narrow emission lines is naturally explained in the context of an accretion flow within the inner rim of the circumbinary disk. We attribute the blueshifted broad emission lines to gas associated with a disk around the accreting secondary black hole. We show that, within the uncertainties, this binary system can be long lived and thus, is not observed in a special moment in time. The orbital motion of the binary can potentially be observed with the VLBA if at least the secondary black hole is a radio emitter. In addition, for the parameters quoted above, the orbital motion will result in a ~100 km/s velocity shift of the emission lines on a time scale of about a year, providing a direct observational test for the binary hypothesis.

We model the electromagnetic signatures of massive black hole binaries (MBHBs) with an associated gas component. The method comprises numerical simulations of relativistic binaries and gas coupled with calculations of the physical properties of the emitting gas. We calculate the UV/X-ray and the Halpha light curves and the Halpha emission profiles. The simulations are carried out with a modified version of the parallel tree SPH code Gadget. The heating, cooling, and radiative processes are calculated for two different physical scenarios, where the gas is approximated as a black-body or a solar metallicity gas. The calculation for the solar metallicity scenario is carried out with the photoionization code Cloudy. We focus on sub-parsec binaries which have not yet entered the gravitational radiation phase. The results from the first set of calculations, carried out for a coplanar binary and gas disk, suggest that there are pronounced outbursts in the X-ray light curve during pericentric passages. If such outbursts persist for a large fraction of the lifetime of the system, they can serve as an indicator of this type of binary. The predicted Halpha emission line profiles may be used as a criterion for selection of MBHB candidates from existing archival data. The orbital period and mass ratio of a binary may be inferred after carefully monitoring the evolution of the Halpha profiles of the candidates. The discovery of sub-parsec binaries is an important step in understanding of the merger rates of MBHBs and their evolution towards the detectable gravitational wave window.

Recent numerical relativistic simulations of black hole coalescence suggest that in certain alignments the emission of gravitational radiation can produce a kick of several thousand kilometers per second. This exceeds galactic escape speeds, hence unless there a mechanism to prevent this, one would expect many galaxies that had merged to be without a central black hole. Here we show that in most galactic mergers, torques from accreting gas suffice to align the orbit and spins of both black holes with the large-scale gas flow. Such a configuration has a maximum kick speed <200 km/s, safely below galactic escape speeds. We predict, however, that in mergers of galaxies without much gas, the remnant will be kicked out several percent of the time. We also discuss other predictions of our scenario, including implications for jet alignment angles and X-type radio sources.

We model the electromagnetic emission signatures of massive black hole binaries (MBHBs) with an associated gas component. The method comprises numerical simulations of relativistic binaries and gas coupled with calculations of the physical properties of the emitting gas. We calculate the accretion powered UV/X-ray and Halpha light curves and the Halpha emission line profiles. The simulations have been carried out with a modified version of the parallel tree SPH code Gadget. The heating, cooling, and radiative processes for the solar metallicity gas have been calculated with the photoionization code Cloudy. We investigate gravitationally bound, sub-parsec binaries which have not yet entered the gravitational radiation phase. The results from the first set of calculations, carried out for a coplanar binary and gas disk, suggest that the outbursts in the X-ray light curve are pronounced during pericentric passages and can serve as a fingerprint for this type of binaries if periodic outbursts are a long lived signature of the binary. The Halpha emission-line profiles also offer strong indications of a binary presence and may be used as a criterion for selection of MBHB candidates for further monitoring from existing archival data. The orbital period and mass ratio of a binary could be determined from the Halpha light curves and profiles of carefully monitored candidates. Although systems with the orbital periods studied here are not within the frequency band of the Laser Interferometer Space Antenna (LISA), their discovery is important for understanding of the merger rates of MBHBs and the evolution of such binaries through the last parsec and towards the detectable gravitational wave window.

We have modeled the time-variable profiles of the Balmer alpha emission from the nonaxisymmetric disk and debris tail created in the tidal disruption of a solar-type star by a million solar mass black hole. Two tidal disruption events were simulated using a three dimensional relativistic smoothed particle hydrodynamic code to describe the early evolution of the debris during the first 50-90 days. We have calculated the physical conditions and radiative processes in the debris using the photoionization code CLOUDY. We model the emission-line profiles in the period immediately after the accretion rate onto the black hole becomes significant. We find that the line profiles at these very early stages of the evolution of the postdisruption debris do not resemble the double-peaked profiles expected from a rotating disk, since the debris has not yet settled into such a stable structure. As a result of the uneven distribution of the debris and the existence of a ``tidal tail'' (the stream of returning debris), the line profiles depend sensitively on the orientation of the tail relative to the line of sight. Moreover, the predicted line profiles vary on fairly short time scales (of order hours to days). Given the accretion rate onto the black hole we also model the Balmer alpha light curve from the debris.

We have modeled the time-variable profiles of the Halpha emission line from the non-axisymmetric disk and debris tail created in the tidal disruption of a solar-type star by a million solar mass black hole. Two tidal disruption event simulations were carried out using a three dimensional relativistic smooth-particle hydrodynamic code, to describe the early evolution of the debris during the first fifty to ninety days. We have calculated the physical conditions and radiative processes in the debris using the photoionization code CLOUDY. We model the emission line profiles in the period immediately after the accretion rate onto the black hole became significant. We find that the line profiles at these very early stages of the evolution of the post-disruption debris do not resemble the double peaked profiles expected from a rotating disk since the debris has not yet settled into such a stable structure. As a result of the uneven distribution of the debris and the existence of a ``tidal tail'' (the stream of returning debris), the line profiles depend sensitively on the orientation of the tail relative to the line of sight. Moreover, the predicted line profiles vary on fairly short time scales (of order hours to days). Given the accretion rate onto the black hole we also model the Halpha light curve from the debris and the evolution of the Halpha line profiles in time.

We have modeled the time-variable profiles of the Halpha emission line from the non-axisymmetric disk and debris tail created in the tidal disruption of a solar-type star by a million solar mass black hole. We find that the line profiles at these very early stages of the evolution of the post-disruption debris do not resemble the double peaked profiles expected from a rotating disk since the debris has not yet settled into such a stable structure. The predicted line profiles vary on fairly short time scales (of order hours to days). As a result of the uneven distribution of the debris and the existence of a ``tidal tail'' (the stream of returning debris), the line profiles depend sensitively on the orientation of the tail relative to the line of sight. Given the illuminating UV/X-ray light curve, we also model the Halpha light curve from the debris.

We have obtained adaptive optics, high spatial resolution (0.15 arcsecond) K-band spectra and images of the region around the two active nuclei in NGC 6240 which show the presence of circumnuclear shocks. The data are consistent with the thermal excitation mechanism being the dominant one in the nuclear region. UV fluorescence and associative detachment may also contribute to the fraction of the energy emitted through molecular hydrogen transitions. The near-IR continuum emission appears closely associated with the two active nuclei. The morphological similarities between the near-IR images and the Chandra X-ray images indicate the same mechanisms may be responsible for the emission in near-IR and X-ray band.