by Fabian, A. C. and Wilkins, D. R. and Miller, J. M. and Reis, R. C. and Reynolds, C. S. and Cackett, E. M. and Nowak, M. A. and Pooley, G. G. and Pottschmidt, K. and Sanders, J. S. and Ross, R. R. and Wilms, J.
7 pages, 10 figures, MNRAS in press

The spin of Cygnus X-1 is measured by fitting reflection models to Suzaku data covering the energy band 0.9-400 keV. The inner radius of the accretion disc is found to lie within 2 gravitational radii (r_g=GM/c^2) and a value for the dimensionless black hole spin is obtained of 0.97^{+0.014}_{-0.02}. This agrees with recent measurements using the continuum fitting method by Gou et al. and of the broad iron line by Duro et al. The disc inclination is measured at 23.7^{+6.7}_{-5.4} deg, which is consistent with the recent optical measurement of the binary system inclination by Orosz et al of 27+/-0.8 deg. We pay special attention to the emissivity profile caused by irradiation of the inner disc by the hard power-law source. The X-ray observations and simulations show that the index q of that profile deviates from the commonly used, Newtonian, value of 3 within 3r_g, steepening considerably within 2r_g, as expected in the strong gravity regime.

by Hlavacek-Larrondo, J. and Fabian, A. C. and Edge, A. C. and Hogan, M. T.
9 pages, 3 figures, accepted for publication in MNRAS

We investigate where brightest cluster galaxies (BCGs) sit on the fundamental plane of black hole (BH) activity, an established relation between the X-ray luminosity, the radio luminosity and the mass of a BH. Our sample mostly consists of BCGs that lie at the centres of massive, strong cooling flow clusters, therefore requiring extreme mechanical feedback from their central active galactic nucleus (AGN) to offset cooling of the intracluster plasma (L_mech>10^44-45 erg/s). Based on the BH masses derived from the M_BH-sigma and M_BH-M_K correlations, we find that all of our objects are offset from the plane such that they appear to be less massive than predicted from their X-ray and radio luminosities (to more than a 99 per cent confidence level). For these objects to be consistent with the fundamental plane, the M_BH-sigma and M_BH-M_K correlations therefore seem to underestimate the BH masses of BCGs, on average by a factor of 10. Our results suggest that the standard relationships between BH mass and host galaxy properties no longer hold for these extreme galaxies. Furthermore, our results imply that if these BHs follow the fundamental plane, then many of those that lie in massive, strong cool core clusters must be ultramassive with M_BH>10^10M_sun. This exceeds the largest BH masses known and has important ramifications for our understanding of the formation and evolution of BHs.

by Reynolds, C. S. and Brenneman, L. W. and Lohfink, A. M. and Trippe, M. L. and Miller, J. M. and Fabian, A. C. and Nowak, M. A.
Submitted to the Astrophysical Journal

The analysis of relativistically broadened X-ray spectral features from the inner accretion disk provides a powerful tool for measuring the spin of supermassive black holes (SMBH) in active galactic nuclei (AGN). However, AGN spectra are often complex and careful analysis employing appropriate and self-consistent models are required if one is to obtain robust results. In this paper, we revisit the deep July-2009 Suzaku observation of the Seyfert galaxy NGC3783 in order to study in a rigorous manner the robustness of the inferred black hole spin parameter. Using Monte Carlo Markov Chain (MCMC) techniques, we identify a (partial) modeling degeneracy between the iron abundance of the disk and the black hole spin parameter. We show that the data for NGC3783 strongly require both supersolar iron abundance (Z_Fe=2-4Zsun) and a rapidly spinning black hole (a>0.88). We discuss various astrophysical considerations that can affect the measured abundance. We note that, while the abundance enhancement inferred in NGC3783 is modest, the X-ray analysis of some other objects has found extreme iron abundances. We introduce the hypothesis that the radiative levitation of iron ions in the innermost regions of radiation-dominated AGN disks can enhance the photospheric abundance of iron. We show that radiative levitation is a plausible mechanism in the very inner regions of high accretion rate AGN disks.

by Zorzi, Alejandra F. and Muzzio, Juan C.
10 pages, 5 figures. The Article has been submitted for publication in Monthly Notices of the RAS by the Royal Astronomical Society and Blackwell Publishing. At this time has been accepted but not yet published

We used the N-body code of Hernquist and Ostriker (1992) to build a dozen cuspy ({\gamma}\approx 1) triaxial models of stellar systems through dissipationless collapses of initially spherical distributions of 10^6 particles. We chose four sets of initial conditions that resulted in models morphologically resembling E2, E3, E4 and E5 galaxies, respectively. Within each set, three different seed numbers were selected for the random number generator used to create the initial conditions, so that the three models of each set are statistically equivalent. We checked the stability of our models using the values of their central densities and of their moments of inertia, which turned out to be very constant indeed. The changes of those values were all less than 3 per cent over one Hubble time and, moreover, we show that the most likely cause of those changes are relaxation effects in the numerical code. We computed the six Lyapunov exponents of nearly 5,000 orbits in each model in order to recognize regular, partially and fully chaotic orbits. All the models turned out to be highly chaotic, with less than 25 per cent of their orbits being regular. We conclude that it is quite possible to obtain cuspy triaxial stellar models that contain large fractions of chaotic orbits and are highly stable. The difficulty to build such models with the method of Schwarzschild (1979) should be attributed to the method itself and not to physical causes.

by Williams, Liliya L. R. and Barnes, Eric I. and Hjorth, Jens
27 pages, 18 figures; accepted to MNRAS

Since globular clusters (GCs) are old, low-N systems their dynamics is widely believed to be fully dominated by collisional two-body processes, and their surface brightness profiles are fit by King models. However, for many GCs, especially those with HST-resolved central regions, and `extra-tidal’ features, King models provide poor fits. We suggest that this is partly because collisionless dynamics is also important and contribute to shaping the cluster properties. We show using time-scale and length-scale arguments that except for the very centers of clusters, collisionless dynamics should be more important than collisional. We then fit 38 GCs analyzed by Noyola and Gebhardt (2006) with (collisional) King and (collisionless) DARKexp models over the full available radial range, and find that the latter provide a better fit to 29 GCs; for six of these the fit is at least ~5x better in term of rms. DARKexp models are theoretically derived maximum entropy equilibrium states of self-gravitating collisionless systems and have already been shown to fit the results of dark matter N-body simulations. (We do not attempt fits with ad hoc fitting functions.)

by Kirsten, F. and Vlemmings, W. H. T.
6 pages, 2 figures, accepted for publication in A&A

Intermediate mass black holes (IMBHs) with expected masses M_BH ~ 10^4 M_sun are thought to bridge the gap between stellar mass black holes (M_BH ~ 3 – 100 M_sun) and supermassive black holes found at the centre of galaxies (M_BH > 10^6 M_sun). Until today, no IMBH has been confirmed observationally. The most promising objects to host an IMBH as their central mass are globular clusters. Here, we present high sensitivity multi-epoch 1.6 GHz very long baseline interferometry observations of the globular cluster M15 that has been suggested to host an IMBH. Assuming the IMBH to be accreting matter from its surrounding we expect to detect it as a point source moving with the global motion of the cluster. However, we do not detect any such object within a radius of 6000 AU of the cluster centre in any of the five observations spread over more than one year. This rules out any variability of the putative IMBH on the time scale of one to two months. To get the most stringent upper limit for the flux density of the putative IMBH we concatenate the data of all five epochs. In this data we measure a 3{\sigma} upper flux limit of 10 {\mu}Jy for a central source. We employ the fundamental plane of black hole activity to estimate the mass of the central IMBH candidate. Based on previous X-ray observations of M15 our measurements indicate a 3{\sigma} upper mass limit of ~500 M_sun.

by Dexter, Jason and Fragile, P. Chris
17 pages, 18 figures, submitted to MNRAS; for movies and version with high-res figures see http://astro.berkeley.edu/~jdexter/tiltedsgra

High-resolution, multi-wavelength, and time-domain observations of the Galactic centre black hole candidate, Sgr A*, allow for a direct test of contemporary accretion theory. To date, all models have assumed alignment between the accretion disc and black hole angular momentum axes, but this is unjustified for geometrically thick accretion flows like that onto Sgr A*. Instead, we calculate images and spectra from a set of simulations of accretion flows misaligned (’tilted’) by 15 degrees from the black hole spin axis and compare them with millimetre (mm) to near-infrared (NIR) observations. Non-axisymmetric standing shocks from eccentric fluid orbits dominate the emission, leading to a wide range of possible image morphologies. These effects invalidate previous parameter estimates from model fitting, including estimates of the dimensionless black hole spin, except possibly at low values of spin or tilt. At 1.3mm, the images have crescent morphologies, and the black hole shadow may still be accessible to future mm-VLBI observations. Shock heating leads to high energy electrons (T > 10^12 K), which can naturally produce the observed NIR flux, spectral index, and rapid variability (’flaring’). This NIR emission is uncorrelated with that in the mm, which also agrees with observations.

These are the first models to self-consistently explain the time-variable mm to NIR emission of Sgr A*. Predictions of the model include significant structural changes observable with mm-VLBI on both the dynamical (hour) and Lense-Thirring precession (day-year) timescales; and ~30-50 microarcsecond changes in centroid position from extreme gravitational lensing events during NIR flares, detectable with the future VLT instrument GRAVITY. If the observed NIR emission is caused by shock heating in a tilted accretion disc, then the Galactic centre black hole has a positive, non-zero spin parameter (a > 0).

by Congedo, Giuseppe
PhD thesis defended at University of Trento on 26th March 2012. Advisors: Stefano Vitale, Mauro Hueller. Committee: Eugenio Coccia (Univ. of Rome, Tor Vergata), Philippe Jetzer (Univ. of Z\”urich), Eric Plagnol (APC-CNRS, Paris), Rita Dolesi (Univ. Of Trento)

LISA is the proposed ESA-NASA gravitational wave detector in the 0.1 mHz – 0.1 Hz band. LISA Pathfinder is the down-scaled version of a single LISA arm. The arm — named Doppler link — can be treated as a differential accelerometer, measuring the relative acceleration between test masses. LISA Pathfinder — the in-flight test of the LISA instrumentation — is currently in the final implementation and planned to be launched in 2014. It will set stringent constraints on the ability to put test masses in geodesic motion to within the required differential acceleration of 3\times10^{-14} m s^{-2} Hz^{-1/2} and track their relative motion to within the required differential displacement measurement noise of 9\times10^{-12} m Hz^{-1/2}, around 1 mHz. Given the scientific objectives, it will carry out — for the first time with such high accuracy required for gravitational wave detection — the science of spacetime metrology, in which the Doppler link between two free-falling test masses measures the curvature. This thesis contains a novel approach to the calculation of the Doppler response to gravitational waves. It shows that the parallel transport of 4-vectors records the history of gravitational wave signals. In practice, the Doppler link is implemented with 4 bodies in LISA and 3 bodies in LISA Pathfinder. To compensate for noise sources a control logic is implemented during the measurement. The closed-loop dynamics of LISA Pathfinder can be condensed into operators acting on the motion coordinates, handling the couplings, as well as the cross-talks. The scope of system identification is the optimal calibration of the instrument. This thesis describes some data analysis procedures applied to synthetic experiments and shows the relevance of system identification for the success of LISA Pathfinder in demonstrating the principles of spacetime metrology for all future space-based missions.

by Berti, Emanuele and Gualtieri, Leonardo and Horbatsch, Michael and Alsing, Justin
9 pages, 4 figures

Scalar-tensor theories are among the simplest extensions of general relativity. In theories with light scalars, deviations from Einstein’s theory of gravity are determined by the scalar mass m_s and by a Brans-Dicke-like coupling parameter \omega_{BD}. We show that gravitational-wave observations of nonspinning neutron star-black hole binary inspirals can be used to set upper bounds on the combination m_s/\sqrt{\omega_{BD}}. We estimate via a Fisher matrix analysis that individual observations with signal-to-noise ratio \rho would yield (m_s/\sqrt{\omega_{\rm BD}})(\rho/10)\lesssim 10^{-15}, 10^{-16} and 10^{-19} eV for Advanced LIGO, ET and eLISA, respectively. A statistical combination of multiple observations may further improve this bound.

by Treuthardt, Patrick and Seigar, Marc S. and Sierra, Amber D. and Al-Baidhany, Ismaeel and Salo, Heikki and Kennefick, Daniel and Kennefick, Julia and Lacy, Claud H. S.
17 pages, 1 table, 11 figures, accepted for publication in MNRAS

The discovery of a relationship between supermassive black hole (SMBH) mass and spiral arm pitch angle (P) is evidence that SMBHs are tied to the overall secular evolution of a galaxy. The discovery of SMBHs in late-type galaxies with little or no bulge suggests that an underlying correlation between the dark matter halo concentration and SMBH mass (MBH) exists, rather than between the bulge mass and MBH. In this paper we measure P using a two-dimensional fast fourier transform and estimate the bar pattern speeds of 40 barred spiral galaxies from the Carnegie-Irvine Galaxy Survey. The pattern speeds were derived by estimating the gravitational potentials of our galaxies from Ks-band images and using them to produce dynamical simulation models. The pattern speeds allow us to identify those galaxies with low central dark halo densities, or fast rotating bars, while P provides an estimate of MBH. We find that a wide range of MBH exists in galaxies with low central dark matter halo densities, which appears to support other theoretical results. We also find that galaxies with low central dark halo densities appear to follow more predictable trends in P versus de Vaucouleurs morphological type (T) and bar strength versus T than barred galaxies in general. The empirical relationship between MBH and total gravitational mass of a galaxy (Mtot) allows us to predict the minimum Mtot that will be observationally measured of our fast bar galaxies. These predictions will be investigated in a subsequent paper.

by Nixon, Chris
5 pages, 4 figures. Accepted for publication in MNRAS

In general, when gas accretes on to a supermassive black hole binary it is likely to have no prior knowledge of the binary angular momentum. Therefore a circumbinary disc forms with a random inclination angle, theta, to the binary. It is known that for theta 90 degrees the disc wholly counteraligns if it satisfies cos(theta) 90 degrees and this criterion is not satisfied the same disc may counteralign its inner regions and, on longer timescales, coalign its outer regions. I show that for typical disc parameters, describing an accretion event on to a supermassive black hole binary, a misaligned circumbinary disc is likely to wholly co– or counter–align with the binary plane. This is because the binary angular momentum dominates the disc angular momentum. However with extreme parameters (binary mass ratio M_2/M_1 << 1 or binary eccentricity e ~ 1) the same disc may simultaneously co- and counter-align. It is known that coplanar prograde circumbinary discs are stable. I show that coplanar retrograde circumbinary discs are also stable. A chaotic accretion event on to an SMBH binary will therefore result in a coplanar circumbinary disc that is either prograde or retrograde with respect to the binary plane.

by Lützgendorf, Nora and Kissler-Patig, Markus and Gebhardt, Karl and Baumgardt, Holger and Noyola, Eva and Jalali, Behrang and de Zeeuw, P. Tim and Neumayer, Nadine
12 pages, 9 figures, 2 tables, accepted for publication in A&A

Globular clusters are an excellent laboratory for stellar population and dynamical research. Recent studies have shown that these stellar systems are not as simple as previously assumed. With multiple stellar populations as well as outer rotation and mass segregation they turn out to exhibit high complexity. This includes intermediate-mass black holes which are proposed to sit at the centers of some massive globular clusters. Today’s high angular resolution ground based spectrographs allow velocity-dispersion measurements at a spatial resolution comparable to the radius of influence for plausible IMBH masses, and to detect changes in the inner velocity-dispersion profile. Together with high quality photometric data from HST, it is possible to constrain black-hole masses by their kinematic signatures. We determine the central velocity-dispersion profile of the globular cluster NGC 2808 using VLT/FLAMES spectroscopy. In combination with HST/ACS data our goal is to probe whether this massive cluster hosts an intermediate-mass black hole at its center and constrain the cluster mass to light ratio as well as its total mass. We derive a velocity-dispersion profile from integral field spectroscopy in the center and Fabry Perot data for larger radii. High resolution HST data are used to obtain the surface brightness profile. Together, these data sets are compared to dynamical models with varying parameters such as mass to light ratio profiles and black-hole masses. Using analytical Jeans models in combination with variable M/L profiles from N-body simulations we find that the best fit model is a no black hole solution. After applying various Monte Carlo simulations to estimate the uncertainties, we derive an upper limit of the back hole mass of M_BH < 1 x 10^4 M_SUN (with 95 % confidence limits) and a global mass-to-light ratio of M/L_V = (2.1 +- 0.2) M_SUN/L_SUN.

by Bédorf, Jeroen and Zwart, Simon Portegies
To appear in: European Physical Journal “Special Topics” : “Computer Simulations on Graphics Processing Units” . 18 pages, 8 figures

In this short review we present the developments over the last 5 decades that have led to the use of Graphics Processing Units (GPUs) for astrophysical simulations. Since the introduction of NVIDIA’s Compute Unified Device Architecture (CUDA) in 2007 the GPU has become a valuable tool for N-body simulations and is so popular these days that almost all papers about high precision N-body simulations use methods that are accelerated by GPUs. With the GPU hardware becoming more advanced and being used for more advanced algorithms like gravitational tree-codes we see a bright future for GPU like hardware in computational astrophysics.

by Nyland, Kristina and Marvil, Josh and Wrobel, Joan and Young, Lisa and Zauderer, B. Ashley
Accepted for publication to ApJ

We present the results of deep, high-resolution, 5 GHz Expanded Very Large Array (EVLA) observations of the nearby, dwarf lenticular galaxy and intermediate mass black hole candidate (M ~4.5 x 10^5 M_sun), NGC 404. For the first time, radio emission at frequencies above 1.4 GHz has been detected in this galaxy. We found a modestly resolved source in the NGC 404 nucleus with a total radio luminosity of 7.6 +/- 0.7 x 10^17 W/Hz at 5 GHz and a spectral index from 5 to 7.45 GHz of alpha = -0.88 +/- 0.30. NGC 404 is only the third central intermediate mass black hole candidate detected in the radio regime with subarcsecond resolution. The position of the radio source is consistent with the optical center of the galaxy and the location of a known, hard X-ray point source (Lx ~1.2 x 10^37 erg/s). The faint radio and X-ray emission could conceivably be produced by an X-ray binary, star formation, a supernova remnant or a low-luminosity AGN powered by an intermediate mass black hole. In light of our new EVLA observations, we find that the most likely scenario is an accreting intermediate mass black hole, with other explanations incompatible with the observed X-ray and/or radio luminosities or statistically unlikely.

by Mullaney, J. R. and Daddi, E. and Béthermin, M. and Elbaz, D. and Juneau, S. and Pannella, M. and Sargent, M. T. and Alexander, D. M. and Hickox, R. C.
6 Pages, 3 Figures, Submitted to ApJL

Using X-ray stacking analyses we estimate the average amounts of supermassive black hole (SMBH) growth taking place in star-forming galaxies (SFGs) at z~1 and z~2 as a function of galaxy stellar mass (M*). We find the average rate of SMBH growth taking place in SFGs follows remarkably similar trends with both M* and redshift as the average star-formation rates (SFRs) of these galaxies (i.e., dM_BH/dt ~ M*^(0.86+/-0.39) for the z~1 sample and dM_BH/dt ~ M*^(1.05+/-0.36) for the z~2 sample). It follows that the ratio of SMBH growth rate to SFR is (a) flat with respect to galaxy stellar mass (b) not evolving with redshift and (c) close to the ratio required to maintain/establish a SMBH to stellar mass ratio of ~10^(-3) as also inferred from today’s M_BH-M_Bulge relationship. We interpret this as evidence that SMBHs have, on average, grown in-step with their host galaxies since at least z~2, irrespective of host galaxy mass and AGN triggering mechanism and that the relative growth rates are more important in establishing inferred M_BH-M* relationships than the seed SMBH masses or merger history. Based on these results we speculate that it is the availability of gas reservoirs that regulates both cosmological SMBH growth and star-formation.

by Bédorf, Jeroen and Gaburov, Evghenii and Zwart, Simon Portegies
5 pages, 2 figures. Proceedings of “Advances in Computational Astrophysics: methods, tools and outcomes”, June 13-17, 2011, Cefalu, Sicily, Italy, eds. Capuzzo Dolcetta, Limongi, Tornambe and Giobbi

We present a gravitational hierarchical N-body code that is designed to run efficiently on Graphics Processing Units (GPUs). All parts of the algorithm are executed on the GPU which eliminates the need for data transfer between the Central Processing Unit (CPU) and the GPU. Our tests indicate that the gravitational tree-code outperforms tuned CPU code for all parts of the algorithm and show an overall performance improvement of more than a factor 20, resulting in a processing rate of more than 2.8 million particles per second.

by Doroshenko, V. T. and Sergeev, S. G. and Klimanov, S. A. and Pronik, V. I. and Efimov, Yu. S.
17 pages, 10 figures, accepted for publication in MNRAS

We present results of the optical spectral and photometric observations of the nucleus of Markarian 6 made with the 2.6-m Shajn telescope at the Crimean Astrophysical Observatory. The continuum and emission Balmer line intensities varied more than by a factor of two during 1992-2008. The lag between the continuum and Hbeta emission line flux variations is 21.1+-1.9 days. For the Halpha line the lag is about 27 days but its uncertainty is much larger. We use Monte-Carlo simulation of the random time series to check the effect of our data sampling on the lag uncertainties and we compare our simulation results with those obtained by random subset selection (RSS) method of Peterson et al. (1998). The lag in the high-velocity wings are shorter than in the line core in accordance with the virial motions. However, the lag is slightly larger in the blue wing than in the red wing. This is a signature of the infall gas motion. Probably the BLR kinematic in the Mrk 6 nucleus is a combination of the Keplerian and infall motions. The velocity-delay dependence is similar for individual observational seasons. The measurements of the Hbeta line width in combination with the reverberation lag permits us to determine the black hole mass, M_BH=(1.8+-0.2)x10^8 M_sun. This result is consistent with the AGN scaling relationships between the BLR radius and the optical continuum luminosity (R_BLR is proportional to L^0.5) as well as with the black-hole mass-luminosity relationship (M_BH-L) under the Eddington luminosity ratio for Mrk 6 to be L_bol/L_Edd ~ 0.01.

by Caranicolas, Nicolaos D. and Zotos, Euaggelos E.
Published in Mechanics Research Communications journal

A 3D-dynamical model is constructed for the study of motion in the central regions of a disk galaxy with a double nucleus. Using the results of the 2D-model, we find the regions of initial conditions in the (x,px,z,py)=EJ, (y=pz=0) phase space producing regular or chaotic orbits. The majority of stars are on chaotic orbits. All chaotic orbits come arbitrary close to one or to both nuclei. Regular orbits are those starting near the stable periodic orbits of the 2D-system with small values of z0. All regular orbits circulate around only one of the two nuclei.

by Roth, Nathaniel and Kasen, Daniel and Hopkins, Philip F. and Quataert, Eliot
15 pages, 14 figures. Submitted to ApJ; v2 Corrected spelling and other small typos

Observational and theoretical arguments suggest that the momentum carried in mass outflows from AGN can reach several times L / c, corresponding to outflow rates of hundreds of solar masses per year. Radiation pressure on lines alone may not be sufficient to provide this momentum deposition, and the transfer of reprocessed IR radiation in dusty nuclear gas has been postulated to provide the extra enhancement. The efficacy of this mechanism, however, will be sensitive to multi-dimensional effects such as the tendency for the reprocessed radiation to preferentially escape along sight-lines of lower column density. We use Monte Carlo radiative transfer calculations to determine the radiation force on dusty gas residing within approximately 10 parsecs from an accreting super-massive black hole. We calculate the net rate of momentum deposition in the surrounding gas and estimate the mass-loss rate in the resulting outflow as a function of solid angle for different black hole luminosities, sightline-averaged column densities, clumping parameters, and opening angles of the dusty gas. We find that these dust-driven winds carry momentum fluxes of 1-5 times L / c and correspond to mass-loss rates of 10-100 solar masses per year for a 10^8 solar mass black hole radiating at or near its Eddington limit. These results help to explain the origin of high velocity molecular and atomic outflows in local ULIRGs, and can inform numerical simulations of galaxy evolution including AGN feedback.

by Armijo, Matias Montesinos
21 pages, 9 figures

In this paper I review and discuss the basic concepts of accretion disks, focused especially on the case of accretion disks around black holes. The well known alpha-model is revisited, showing the strengths and weaknesses of the model. Other turbulent viscosity prescription, based on the Reynolds number, that may improve our understanding of the accretion paradigm is discussed. A simple but efficient mathematical model of a self-gravitating accretion disk, as well as observational evidence of these objects, are also included.

by Bromley, Benjamin C. and Kenyon, Scott J. and Geller, Margaret J. and Brown, Warren R.
ApJ Letters, accepted; 13 pages, including 3 figures

We examine whether disrupted binary stars can fuel black hole growth. In this mechanism, tidal disruption produces a single hypervelocity star (HVS) ejected at high velocity and a former companion star bound to the black hole. After a cluster of bound stars forms, orbital diffusion allows the black hole to accrete stars by tidal disruption at a rate comparable to the capture rate. In the Milky Way, HVSs and the S star cluster imply similar rates of 10^{-5}–10^{-3} yr^{-1} for binary disruption. These rates are consistent with estimates for the tidal disruption rate in nearby galaxies and imply significant black hole growth from disrupted binaries on 10 Gyr time scales.

by Strader, Jay and Chomiuk, Laura and Maccarone, Thomas and Miller-Jones, James and Seth, Anil and Heinke, Craig and Sivakoff, Gregory
6 pages, 2 figures. ApJL in press

With a goal of searching for accreting intermediate-mass black holes (IMBHs), we report the results of ultra-deep Jansky VLA radio continuum observations of the cores of three Galactic globular clusters: M15, M19, and M22. We reach rms noise levels of 1.5-2.1 uJy/beam at an average frequency of 6 GHz. No sources are observed at the center of any of the clusters. For a conservative set of assumptions about the properties of the accretion, we set 3-sigma upper limits on IMBHs from 360-980 M_sun. These limits are among the most stringent obtained for any globular cluster. They add to a growing body of work that suggests either (a) IMBHs ~> 1000 M_sun are rare in globular clusters, or (b) when present, IMBHs accrete in an extraordinarily inefficient manner.

by Giacomazzo, Bruno and Baker, John G. and Miller, M. Coleman and Reynolds, Christopher S. and van Meter, James R.
6 pages, 5 figures

Coalescing supermassive black hole binaries are produced by the mergers of galaxies and are the most powerful sources of gravitational waves accessible to space-based gravitational observatories. Some such mergers may occur in the presence of matter and magnetic fields and hence generate an electromagnetic counterpart. In this paper we present the first general relativistic simulations of magnetized plasma around merging supermassive black holes using the general relativistic magnetohydrodynamic code Whisky. By considering different magnetic field strengths, going from non-magnetically dominated to magnetically dominated regimes, we explore how magnetic fields affect the dynamics of the plasma and the possible emission of electromagnetic signals. In particular we observe a total amplification of the magnetic field of ~2 orders of magnitude which is driven by the accretion onto the binary and that leads to much stronger electromagnetic signals, more than a factor of 10^4 larger than comparable calculations done in the force-free regime where such amplifications are not possible.

by Haiman, Zoltán
47 pages with 306 references; this review is a chapter in “The First Galaxies – Theoretical Predictions and Observational Clues”, Springer Astrophysics and Space Science Library, Eds. T. Wiklind, V. Bromm & B. Mobasher, in press

Supermassive black holes (SMBHs) are common in local galactic nuclei, and SMBHs as massive as several billion solar masses already exist at redshift z=6. These earliest SMBHs may grow by the combination of radiation-pressure-limited accretion and mergers of stellar-mass seed BHs, left behind by the first generation of metal-free stars, or may be formed by more rapid direct collapse of gas in rare special environments where dense gas can accumulate without first fragmenting into stars. This chapter offers a review of these two competing scenarios, as well as some more exotic alternative ideas. It also briefly discusses how the different models may be distinguished in the future by observations with JWST, (e)LISA and other instruments.

by Steier, Frank and Cervantes, Felipe Guzmán and Marín, Antonio F. García and Gerardi, Domenico and Heinzel, Gerhard and Danzmann, Karsten

LISA Pathfinder is a technology demonstration mission for the Laser Interferometer Space Antenna (LISA). The main experiment on-board LISA Pathfinder is the so-called LISA Technology Package (LTP) which has the aim to measure the differential acceleration between two free-falling test masses with an accuracy of 3×10^(-14) ms^(-2)/sqrt[Hz] between 1 mHz and 30 mHz. This measurement is performed interferometrically by the Optical Metrology System (OMS) on-board LISA Pathfinder. In this paper we present the development of an experimental end-to-end testbed of the entire OMS. It includes the interferometer and its sub-units, the interferometer back-end which is a phasemeter and the processing of the phasemeter output data. Furthermore, 3-axes piezo actuated mirrors are used instead of the free-falling test masses for the characterisation of the dynamic behaviour of the system and some parts of the Drag-free and Attitude Control System (DFACS) which controls the test masses and the satellite. The end-to-end testbed includes all parts of the LTP that can reasonably be tested on earth without free-falling test masses. At its present status it consists mainly of breadboard components. Some of those have already been replaced by Engineering Models of the LTP experiment. In the next steps, further Engineering Models and Flight Models will also be inserted in this testbed and tested against well characterised breadboard components. The presented testbed is an important reference for the unit tests and can also be used for validation of the on-board experiment during the mission.

by Cervantes, Felipe Guzmán and Steier, Frank and Wanner, Gudrun and Heinzel, Gerhard and Danzmann, Karsten
Applied Physics B – Lasers and Optics (2008)

We present recent sensitivity measurements of the LISA Technology Package interferometer with articulated mirrors as test masses, actuated by piezo-electric transducers. The required longitudinal displacement resolution of 9 pm/sqrt[Hz] above 3 mHz has been demonstrated with an angular noise that corresponds to the expected in on-orbit operation. The excess noise contribution of this test mass jitter onto the sensitive displacement readout was completely subtracted by fitting the angular interferometric data streams to the longitudinal displacement measurement. Thus, this cross-coupling constitutes no limitation to the required performance of the LISA Technology Package interferometry.

by Sutton, Andrew D. and Roberts, Timothy P. and Walton, Dominic J. and Gladstone, Jeanette C. and Scott, Amy E.
Accepted for publication in MNRAS, 26 pages, 9 figures, 10 tables

We present the results from an X-ray and optical study of a new sample of eight extreme luminosity ultraluminous X-ray source (ULX) candidates, which were selected as the brightest ULXs (with L_X > 5×10^40 erg/s) located within 100 Mpc identified in a cross correlation of the 2XMM-DR1 and RC3 catalogues. These objects are so luminous that they are difficult to describe with current models of super-Eddington accretion onto all but the most massive stellar remnants; hence they are amongst the most plausible candidates to host larger, intermediate-mass black holes (IMBHs). Two objects are luminous enough in at least one observation to be classed as hyperluminous X-ray source (HLX) candidates, including one persistent HLX in an S0 galaxy that (at 3×10^41 erg/s) is the second most luminous HLX yet detected. The remaining seven sources are located in spiral galaxies, and several appear to be closely associated with regions of star formation as is common for many less luminous ULXs. However, the X-ray characteristics of these extreme ULXs appear to diverge from the less luminous objects. They are typically harder, possessing absorbed power-law continuum spectra with photon indexes ~ 1.7, and are potentially more variable on short timescales, with data consistent with ~ 10-20 per cent rms variability on timescales of 0.2-2 ks. These properties appear consistent with the sub-Eddington hard state, which given the observed luminosities of these objects suggests the presence of IMBHs with masses in the range 10^3-10^4 M_Sun. As such, this strengthens the case for these brightest ULXs as good candidates for the eventual conclusive detection of the highly elusive IMBHs. However, we caution that a combination of the highest plausible super-Eddington accretion rates and the largest permitted stellar black hole remnants cannot be ruled out without future, improved observations.

by Baruteau, Clément and Ramirez-Ruiz, Enrico and Masset, Frédéric
5 pages, 5 figures, accepted for publication in MNRAS

We present two-dimensional hydrodynamical simulations of the tidal interaction between a supermassive black hole binary with moderate mass ratio, and the fossil gas disc where it is embedded. Our study extends previous one-dimensional height-integrated disc models, which predicted that the density of the gas disc between the primary and the secondary black holes should rise significantly during the ultimate stages of the binary’s hardening driven by the gravitational radiation torque. This snow-plough mechanism, as we call it, would lead to an increase in the bolometric luminosity of the system prior to the binary merger, which could be detected in conjunction with the gravitational wave signal. We argue here that the snow-plough mechanism is unlikely to occur. In two-dimensions, when the binary’s hardening timescale driven by gravitational radiation becomes shorter than the disc’s viscous drift timescale, fluid elements in the inner disc get funneled to the outer disc through horseshoe trajectories with respect to the secondary. Mass leakage across the secondary’s gap is thus found to be effective and, as a result, the predicted accretion disc luminosity will remain at roughly the same level prior to merger.

by Nayakshin, Sergei and Power, Chris and King, Andrew R.
8 pages; submitted to ApJ; version that takes referee’s comments into account

We argue that current observations of $latex M – \sigma$ relations for galaxies can be used to constrain theories of super-massive black holes (SMBH) feeding. In particular, assuming that SMBH mass is limited only by the feedback on the gas that feeds it, we show that SMBHs fed via a planar galaxy scale gas flow, such as a disc or a bar, should be much more massive than their counterparts fed by quasi-spherical inflows. This follows from the relative inefficiency of AGN feedback on a flattened inflow. We find that even under the most optimistic conditions for SMBH feedback on flattened inflows, the mass at which the SMBH expels the gas disc and terminates its own growth is a factor of several higher than the one established for quasi-spherical inflows. Any beaming of feedback away from the disc and any disc self-shadowing strengthens this result further. Contrary to this theoretical expectation, recent observations have shown that SMBH in pseudobulge galaxies (which are associated with barred galaxies) are typically under- rather than over-massive when compared with their classical bulge counterparts at a fixed value of $latex \sigma$. We conclude from this that SMBHs are not fed by large (100 pc to many kpc) scale gas discs or bars, most likely because such planar flows are turned into stars too efficiently to allow any SMBH growth. Based on this and other related observational evidence, we argue that most SMBHs grow by chaotic accretion of gas clouds with a small and nearly randomly distributed direction of angular momentum.

by Berti, Emanuele and Kesden, Michael and Sperhake, Ulrich
12 pages, 4 figures, 1 table

Recent numerical relativity simulations have shown that the final black hole produced in a binary merger can recoil with a velocity as large as 5,000 km/s. Because of enhanced gravitational-wave emission in the so-called “hang-up” configurations, this maximum recoil occurs when the black-hole spins are partially aligned with the orbital angular momentum. We revisit our previous statistical analysis of post-Newtonian evolutions of black-hole binaries in the light of these new findings. We demonstrate that despite these new configurations with enhanced recoil velocities, spin alignment during the post-Newtonian stage of the inspiral will still significantly suppress (or enhance) kick magnitudes when the initial spin of the more massive black hole is more (or less) closely aligned with the orbital angular momentum than that of the smaller hole. We present a preliminary study of how this post-Newtonian spin alignment affects the ejection probabilities of supermassive black holes from their host galaxies with astrophysically motivated mass ratio and initial spin distributions. We find that spin alignment suppresses (enhances) ejection probabilities by ~ 40% (20%) for an observationally motivated mass-dependent galactic escape velocity, and by an even greater amount for a constant escape velocity of 1,000 km/s. Kick suppression is thus at least a factor two more efficient than enhancement.

by Cervantes, Felipe Guzmán and Livas, Jeffrey and Silverberg, Robert and Buchanan, Ernest and Stebbins, Robin

LISA will use quadrant photoreceivers as front-end devices for the phasemeter measuring the motion of drag-free test masses in both angular orientation and separation. We have set up a laboratory testbed for the characterization of photoreceivers. Some of the limiting noise sources have been identified and their contribution has been either measured or derived from the measured data. We have built a photoreceiver with a 0.5 mm diameter quadrant photodiode with an equivalent input current noise of better than 1.8 pA/sqrt[Hz] below 20 MHz and a 3dB bandwidth of 34 MHz.

by Sabha, Nadeen and Eckart, Andreas and Merritt, David and Zamaninasab, Mohammad and Witzel, Gunther and García-Marín, Macarena and Jalali, Behrang and Valencia-S., Monica and Yazici, Senol and Buchholz, Rainer and Shahzamanian, Banafsheh and Straubmeier, Christian
13 pages, 11 figures, submitted to A&A

Sagittarius A*, the super-massive black hole at the center of the Milky Way, is surrounded by a small cluster of high velocity stars, known as the S-stars. We aim at constraining the amount and nature of the stellar and dark mass associated with the cluster in the immediate vicinity of Sagittarius A*. We use near-infrared imaging to determine the $latex K_s$-band luminosity function of the S-star cluster members, and use the distribution of the diffuse background emission and the stellar number density counts around the central black hole. This allows us to determine the stellar light and mass contribution that we can expect from the faint members of the cluster. We then use post-Newtonian N-body techniques to investigate the effect of stellar perturbations on the motion of S2, as a means of detecting the number and masses of the perturbers. We find that the stellar mass derived from the $latex K_s$-band luminosity extrapolation is much smaller than the amount of mass that might be present considering the uncertainties in the orbital motion of the star S2. Also the amount of light from the fainter S-cluster members is below the amount of residual light at the position of the S-star cluster after removing the bright cluster members. If the distribution of stars and stellar remnants is strongly enough peaked near Sagittarius A*, observed changes in the orbital elements of S2 can be used to constrain both their masses and numbers. Based on simulations of the cluster of high velocity stars we find that at a wavelength of 2.2 $latex \mu$m close to the confusion level for 8 m class telescopes blend stars will occur (preferentially near the position of Sagittarius A*) that last for typically 3 years before they dissolve due to proper motions.

by Straub, Odele and Vincent, Frederic H. and Abramowicz, Marek A. and Gourgoulhon, Eric and Paumard, Thibaut
11 pages, 10 figures, submitted to A & A

We calculate the “observed at infinity” image and spectrum of the accretion structure in Sgr A*, by modelling it as an optically thin, constant angular momentum ion torus in hydrodynamic equilibrium. The physics we consider includes a two-temperature plasma, a toroidal magnetic field, as well as radiative cooling by bremsstrahlung, synchrotron and inverse Compton processes. Our relativistic model has the virtue of being fully analytic and very simple, depending only on eight tunable parameters: the black hole spin and the inclination of the spin axis to our line of sight, the torus angular momentum, the polytropic index, the magnetic to total pressure ratio, the central values of density and electron temperature and the ratio of electron to ion temperatures. The observed image and spectrum are calculated numerically using the ray-tracing code GYOTO. Our results demonstrate that the ion torus model is able to account for the main features of the accretion structure surrounding Sgr A*.

by Hosokawa, Takashi and Omukai, Kazuyuki and Yorke, Harold W.
11 pages, 10 figures

Direct collapse of supermassive stars (SMSs) is a possible pathway for generating supermassive black holes in the early universe. It is expected that an SMS could form via very rapid mass accretion with Mdot ~ 0.1 – 1 Msun/yr during the gravitational collapse of an atomic-cooling primordial gas cloud. In this paper we study how stars would evolve under such extreme rapid mass accretion, focusing on the early evolution until the stellar mass reaches 1000 Msun. To this end we numerically calculate the detailed interior structure of accreting stars with primordial element abundances. Our results show that for accretion rates higher than 0.01 Msun/yr, stellar evolution is qualitatively different from that expected at lower rates. While accreting at these high rates the star always has a radius exceeding 100 Rsun, which increases monotonically with the stellar mass. The mass-radius relation for stellar masses exceeding ~ 100 Msun follows the same track with R_* \propto M_*^0.5 in all cases with accretion rates > 0.01 Msun/yr; at a stellar mass of 1000 Msun the radius is about 7000 Rsun (~= 30 AU). With higher accretion rates the onset of hydrogen burning is shifted towards higher stellar masses. In particular, for accretion rates exceeding Mdot > 0.1 Msun/yr, there is no significant hydrogen burning even after 1000 Msun have accreted onto the protostar. Such “supergiant” protostars have effective temperatures as low as Teff ~= 5000 K throughout their evolution and because they hardly emit ionizing photons, they do not create an HII region or significantly heat their immediate surroundings. Thus, radiative feedback is unable to hinder the growth of rapidly accreting stars to masses in excess of 1000 Msun, as long as material is accreted at rates Mdot > 0.01 Msun/yr.

by Mishra, Chandra Kant
17 pages, 2 figures

The loss rate of linear momentum from a binary system composed of compact objects (radially falling towards each other under mutual gravitational influence) has been investigated using the multipolar post-Minkowskian approach. The 2.5PN accurate analytical formula for the linear momentum flux is provided, in terms of the separation of the two objects, in harmonic coordinates, both for a finite and infinite initial separation. The 2.5PN formulas for the linear momentum flux are finally used to estimate the recoil velocity accumulated during a premerger phase of the binary evolution.

by Roedig, Constanze and Sesana, Alberto and Dotti, Massimo and Cuadra, Jorge and Amaro-Seoane, Pau and Haardt, Francesco
14 pages, 12 Figures, submitted to MNRAS

We analyse 3D SPH simulations of the evolution of initially quasi-circular massive black hole binaries (BHBs) residing in the central hollow (cavity) of self-gravitating circumbinary discs. We perform a set of simulations adopting different thermodynamics for the gas within the cavity and for the ‘numerical size’ of the black holes. We study the interplay between gas accretion and gravity torques in changing the binary elements (semi-major axis and eccentricity) and its total angular momentum budget. We pay special attention to the gravity torques, by analysing their physical origin and location. We show that (i) the BHB eccentricity grows due to gravity torques from the inner edge of the disc, independently of the accretion and the adopted thermodynamics; (ii) the semi-major axis decay depends not only on the gravity torques but also on their subtle interplay with the disc-binary angular momentum transfer due to accretion; (iii) the spectral structure of the gravity torques is predominately caused by disc edge overdensities and spiral arms developing in the body of the disc; (iv) the net gravity torque changes sign across the BHB corotation radius: gas inside this radius exerts a net positive torque, while streams located outside this radius (but within the cavity) exert a net negative torque. The relative importance of the two might depend on the thermodynamical properties of the instreaming gas and is crucial in assessing the disc–binary angular momentum transfer; (v) the net torque manifests as a purely kinematic effect as it stems from the low density cavity, where the material flows in and out in highly eccentric orbits. Thus both accretion onto the black holes and the interaction with gas streams inside the cavity must be taken into account to assess the fate of the binary.

by Saxton, R. D. and Read, A. M. and Esquej, P. and Komossa, S. and Dougherty, S. and Rodriguez-Pascual, P. and Barrado, D.
Accepted for publication in Astronomy and Astrophysics. 11 pages, 11 figures

SDSS J120136.02+300305.5 was detected in an XMM-Newton slew from June 2010 with a flux 56 times higher than an upper limit from ROSAT, corresponding to Lx~3×10^44 ergs/s. It has the optical spectrum of a quiescent galaxy (z=0.146). Overall the X-ray flux has evolved consistently with the canonical t^-5/3 model, expected for returning stellar debris from a tidal disruption event, fading by a factor ~300 over 300 days. In detail the source is very variable and became invisible to Swift between 27 and 48 days after discovery, perhaps due to self-absorption. The X-ray spectrum is soft but is not the expected tail of optically thick thermal emission. It may be fit with a Bremsstrahlung or double-power-law model and is seen to soften with time and declining flux. Optical spectra taken 12 days and 11 months after discovery indicate a deficit of material in the broad line and coronal line regions of this galaxy, while a deep radio non-detection implies that a jet was not launched during this event.

by Inayoshi, Kohei and Omukai, Kazuyuki
9 pages, 4 figures, accepted for publication in MNRAS

We propose a new scenario for supermassive star (SMS;>10^5Msun) formation in shocked regions of colliding cold accretion flows near the centers of first galaxies. Recent numerical simulations indicate that assembly of a typical first galaxy with virial temperature (~10^4K) proceeds via cold and dense flows penetrating deep to the center, where the supersonic streams collide each other to develop a hot and dense (~10^4K, ~10^3/cc) shocked gas. The post-shock layer first cools by efficient Ly alpha emission and contracts isobarically until 8000K. Whether the layer continues the isobaric contraction depends on the density at this moment: if the density is high enough for collisionally exciting H2 rovibrational levels (>10^4/cc), enhanced H2 collisional dissociation suppresses the gas to cool further. In this case, the layer fragments into massive (>10^5Msun) clouds, which collapse isothermally (~8000K) by the Ly alpha cooling without subsequent fragmentation. As an outcome, SMSs are expected to form and evolve eventually to seeds of supermassive black holes (SMBH). By calculating thermal evolution of the post-shock gas, we delimit the range of post-shock conditions for the SMS formation, which can be expressed as: T>6000K/(n/10^4/cc) for n5000-6000K for n>10^4/cc, depending somewhat on initial ionization degree. We found that metal enrichment does not affect the above condition for metallicity below 10^-3Zsun if metals are in the gas phase, while condensation of several percent of metals into dust decreases this critical value of metallicity by an order of magnitude. Unlike the previously proposed scenario for SMS formation, which postulates extremely strong ultraviolet radiation to quench H2 cooling, our scenario naturally explains the SMBH seed formation in the assembly process of the first galaxies, even without such a strong radiation.

by Wegg, Christopher
4 pages, 1 figure. Accepted by ApJ

We describe a pseudo-Newtonian potential which, to within 1% error at all angular momenta, reproduces the precession due to general relativity of particles whose specific orbital energy is small compared to c^2 in the Schwarzschild metric. For bound orbits the constraint of low energy is equivalent to requiring the apoapsis of a particle to be large compared to the Schwarzschild radius. Such low energy orbits are ubiquitous close to supermassive black holes in galactic nuclei, but the potential is relevant in any context containing particles on low energy orbits. Like the more complex post-Newtonian expressions, the potential correctly reproduces the precession in the far-field, but also correctly reproduces the position and magnitude of the logarithmic divergence in precession for low angular momentum orbits. An additional advantage lies in its simplicity, both in computation and implementation. We also provide two simpler, but less accurate potentials, for cases where orbits always remain at large angular momenta, or when the extra accuracy is not needed. In all of the presented cases the accuracy in precession in low energy orbits exceeds that of the well known potential of Paczynski & Wiita (1980), which has ~30% error in the precession at all angular momenta.

by Miralda-Escude, Jordi
submitted to ApJ

We propose that the cloud of gas moving on a highly eccentric orbit around the central black hole in our Galaxy, reported by Gillessen et al., is produced by a wind from photoevaporating debris orbiting around a star with a small circumstellar disk. The disk is tidally truncated to less than 1 AU at the peribothron passage, and a cloud like the observed one is recreated by the wind at every orbit. The star-disk system, which may have been producing the cloud for hundreds of orbits in the past, is proposed to have formed when the star flew by a stellar black hole and was tidally disrupted and deflected to its present orbit. Encounters of low-mass stars with stellar black holes are likely to occur at the location of this cloud, because of the high density of stellar black holes expected to have migrated to the Galactic center by mass segregation. The rate of these encounters at a small enough impact parameter to disrupt the star may reasonably be ~ 10^{-6} per year. The flyby should have spun up the star and pulled out a substantial fraction of its mass as tidal debris, part of which fell back onto the star and created a small disk. Since then, the disk may have expanded by absorbing angular momentum from the star up to the tidal truncation radius. Thereafter, the strong tidal perturbation of the outer disk edge at every peribothron may create gas streams moving out to larger radius that can photoevaporate and generate the wind that produces the cloud at every orbit. The model predicts that when the cloud is disrupted at the next peribothron passage in 2013, a smaller unresolved cloud will follow the star on the same orbit that will gradually grow. An increased infrared luminosity from the disk may also become detectable during the peribothron passage.

by Hobbs, Alexander and Power, Chris and Nayakshin, Sergei and King, Andrew
7 pages, 3 figures. Accepted for publication in MNRAS

Accretion onto supermassive black holes (SMBHs) in galaxy formation simulations is frequently modelled by the Bondi-Hoyle formalism. Here we examine the validity of this approach analytically and numerically. We argue that the character of the flow where one evaluates the gas properties is unlikely to satisfy the simple Bondi-Hoyle model. Only in the specific case of hot virialised gas with zero angular momentum and negligible radiative cooling is the Bondi-Hoyle solution relevant. In the opposite extreme, where the gas is in a state of free-fall at the evaluation radius due to efficient cooling and the dominant gravity of the surrounding halo, the Bondi-Hoyle formalism can be erroneous by orders of magnitude in either direction. This may impose artificial trends with halo mass in cosmological simulations by being wrong by different factors for different halo masses. We propose an expression for the sub-grid accretion rate which interpolates between the free-fall regime and the Bondi-Hoyle regime, therefore taking account of the contribution of the halo to the gas dynamics.

by Aarseth, Sverre
8 pages, 9 figures, accepted MNRAS

We report on results of fully consistent N-body simulations of globular cluster models with N = 100 000 members containing neutron stars and black holes. Using the improved `algorithmic regularization’ method of Hellstrom and Mikkola for compact subsystems, the new code NBODY7 enables for the first time general relativistic coalescence to be achieved for post-Newtonian terms and realistic parameters. Following an early stage of mass segregation, a few black holes form a small dense core which usually leads to the formation of one dominant binary. The subsequent evolution by dynamical shrinkage involves the competing processes of ejection and mergers by radiation energy loss. Unless the binary is ejected, long-lived triple systems often exhibit Kozai cycles with extremely high inner eccentricity (e > 0.999) which may terminate in coalescence at a few Schwarzschild radii. A characteristic feature is that ordinary stars as well as black holes and even BH binaries are ejected with high velocities. On the basis of the models studied so far, the results suggest a limited growth of a few remaining stellar mass black holes in globular clusters.

by Kawasaki, Masahiro and Kusenko, Alexander and Yanagida, Tsutomu T.
12 pages, 3 figures

Supermassive black holes exist in the centers of galaxies, including Milky Way, but there is no compelling theory of their formation. Furthermore, observations of quasars imply that supermassive black holes have already existed at some very high redshifts, suggesting the possibility of their primordial origin. In a class of well-motivated models, inflationary epoch could include two or more periods of inflation dominated by different scalar fields. The transition between such periods of inflation could enhance the spectrum of density perturbations on some specific scale, which could lead to formation of primordial black holes with a very narrow range of masses of the order of 0.1 to 1 million solar masses. These primordial black holes could have provided the requisite seeds for the observed population of supermassive black holes.

by Li, Yan-Rong and Wang, Jian-Min and Ho, Luis C.
To appear in the ApJ, 11 pages and 9 figures

The spin is an important but poorly constrained parameter for describing supermassive black holes (SMBHs). Using the continuity equation of SMBH number density, we explicitly obtain the mass-dependent cosmological evolution of the radiative efficiency for accretion, which serves as a proxy for SMBH spin. Our calculations make use of the SMBH mass function of active and inactive galaxies (derived in the first paper of this series), the bolometric luminosity function of active galactic nuclei (AGNs), corrected for the contribution from Compton-thick sources, and the observed Eddington ratio distribution. We find that the radiative efficiency generally increases with increasing black hole mass at high redshifts (z>~1), roughly as \eta \propto M_bh^0.5, while the trend reverses at lower redshifts, such that the highest efficiencies are attained by the lowest mass black holes. Black holes with M_bh>~10^8.5M_sun maintain radiative efficiencies as high as \eta~0.3-0.4 at high redshifts, near the maximum for rapidly spinning systems, but their efficiencies drop dramatically (by an order of magnitude) by z~0. The pattern for lower mass holes is somewhat more complicated but qualitatively similar. Assuming that the standard accretion disk model applies, we suggest that the accretion history of SMBHs and their accompanying spins evolve in two distinct regimes: an early phase of prolonged accretion, plausibly driven by major mergers, during which the black hole spins up, then switching to a period of random, episodic accretion, governed by minor mergers and internal secular processes, during which the hole spins down. The transition epoch depends on mass, mirroring other evidence for “cosmic downsizing” in the AGN population; it occurs at z~2 for high-mass black holes, and somewhat later, at z~1, for lower-mass systems.

by Yagi, Kent
32 pages, 13 figures, accepted to CQG

DECIGO Path Finder (DPF) is a space-borne gravitational wave (GW) detector with sensitivity in the frequency band 0.1–100Hz. As a first step mission to DECIGO, it is aiming for launching in 2016–2017. Although its main objective is to demonstrate technology for GW observation in space, DPF still has a chance of detecting GW signals and performing astrophysical observations. With an observable range up to 50 kpc, its main targets are GW signals from galactic intermediate mass black hole (IMBH) binaries. By using inspiral-merger-ringdown phenomenological waveforms, we perform both pattern-averaged analysis and Monte Carlo simulations including the effect of detector motion to find that the masses and (effective) spins of the IMBHs could be determined with errors of a few percent, should the signals be detected. Since GW signals from IMBH binaries with masses above $latex 10^4 M_\odot$ cannot be detected by ground-based detectors, these objects can be unique sources for DPF. If the inspiral signal of a $latex 10^3M_\odot$ IMBH binary is detected with DPF, it can give alert to the ringdown signal for the ground-based detectors $latex 10^2$–$latex 10^3$s before coalescence. We also estimate the possible bound on the graviton Compton wavelength from a possible IMBH binary in $latex \omega$ Centauri. We obtain a slightly weaker constraint than the solar system experiment and an about 2 orders of magnitude stronger constraint than the one from binary pulsar tests. Unfortunately, the detection rate of IMBH binaries is rather small.

by Tang, Alf and Sumner, Timothy J.
10 pages, 7 figures

Although cosmic string scenario for galaxy formation is disfavored by CMB data, it is of great interest in the generation of cosmic gravitational-wave background. This research aims to develop an algorithm to extract cosmic gravitational-wave background produced by cosmic strings from the LISA data stream, and apply the algorithm to the simulated data stream containing the background produced by cosmic strings with various strength to study the detection threshold for this source. For 1-yr observation, It is found that the detection threshold of G{\mu} is 3.12 \times 10^-16 in the standard scenario. In the case that p and {\epsilon} are adjustable, the detectable region in parameter space is defined by (G{\mu})^2/3 {\epsilon}^-1/3 / p> 4.6 \times 10-11.

by Binney, James
43 pages. Lectures given at the XXIII Canary islands winter school of astrophysics, November 2011 To appear in “Secular Evolution of Galaxies” ed. J. Falcon-Barroso & J.H. Knapen

The text of lectures to the 2011 Tenerife Winter School. The School’s theme was “Secular Evolution of Galaxies” and my task was to present the underlying stellar-dynamical theory. Other lecturers were speaking on the role of bars and chemical evolution, so these topics are avoided here. We start with an account of the connections between isolating integrals, quasiperiodicity and angle-action variables – these variables played a unifying role throughout the lectures. This leads on to the phenomenon of resonant trapping and how this can lead to chaos in cuspy potentials and phase-space mixing in slowly evolving potentials. Surfaces of section and frequency analysis are introduced as diagnostics of phase-space structure. Real galactic potentials include a fluctuating part that drives the system towards unattainable thermal equilibrium. Two-body encounters are only one source of fluctuations, and all fluctuations will drive similar evolution. We derive the orbit-averaged Fokker-Planck equation and relations that hold between the second-order diffusion coefficients and both the power spectrum of the fluctuations and the first-order diffusion coefficients. From the observed heating of the solar neighbourhood we show that the second-order diffusion coefficients must scale as J^{1/2}. We show that periodic spiral structure shifts angular momentum outwards, heating at the Lindblad resonances and mixing at corotation. The equation that would yield the normal modes of a stellar disc is first derived and then used to discuss the propagation of tightly-wound spiral waves. The winding up of such waves is explains why cool stellar discs are responsive systems that amplify ambient noise. An explanation is offered of why the Lin-Shu-Kalnajs dispersion relation and even global normal-mode calculations provide a very incomplete understanding of the dynamics of stellar discs.

by Gultekin, Kayhan and Cackett, Edward M. and Miller, Jon M. and Di Matteo, Tiziana and Markoff, Sera and Richstone, Douglas O.
Accepted by ApJ. 16 pages, 8 figures, 5 tables

We present Chandra observations of 12 galaxies that contain supermassive black holes with dynamical mass measurements. Each galaxy was observed for 30 ksec and resulted in a total of 68 point source detections in the target galaxies including supermassive black hole sources, ultraluminous X-ray sources, and extragalactic X-ray binaries. Based on our fits of the X-ray spectra, we report fluxes, luminosities, Eddington ratios, and slope of the power-law spectrum. Normalized to the Eddington luminosity, the 2–10 keV band X-ray luminosities of the SMBH sources range from $latex 10^{-8}$ to $latex 10^{-6}$, and the power-law slopes are centered at $latex \sim2$ with a slight trend towards steeper (softer) slopes at smaller Eddington fractions, implying a change in the physical processes responsible for their emission at low accretion rates. We find 20 ULX candidates, of which six are likely ($latex >90%$ chance) to be true ULXs. The most promising ULX candidate has an isotropic luminosity in the 0.3–10 keV band of $latex 1.0_{-0.3}^{+0.6} \times 10^{40}$ erg/s.

by Petri, Andrea and Ferrara, Andrea and Salvaterra, Ruben
12 pages, 9 figures, 1 table, MNRAS in press

We study a model in which supermassive black holes (SMBHs) can grow by the combined action of gas accretion on heavy seeds and mergers of both heavy (m_s^h=10^5 Msol) and light (m_s^l = 10^2 Msol) seeds. The former result from the direct collapse of gas in T_s^h >1.5×10^4K, H_2-free halos; the latter are the endproduct of a standard H_2-based star formation process. The H_2-free condition is attained by exposing halos to a strong (J_21 > 10^3) Lyman-Werner UV background produced by both accreting BHs and stars, thus establishing a self-regulated growth regime. We find that this condition is met already at z close to 18 in the highly biased regions in which quasars are born. The key parameter allowing the formation of SMBHs by z=6-7 is the fraction of halos that can form heavy seeds: the minimum requirement is that f_heavy>0.001; SMBH as large as 2×10^10 Msol can be obtained when f_heavy approaches unity. Independently of f_heavy, the model produces a high-z stellar bulge-black hole mass relation which is steeper than the local one, implying that SMBHs formed before their bulge was in place. The formation of heavy seeds, allowed by the Lyman-Werner radiative feedback in the quasar-forming environment, is crucial to achieve a fast growth of the SMBH by merger events in the early phases of its evolution, i.e. z>7. The UV photon production is largely dominated by stars in galaxies, i.e. black hole accretion radiation is sub-dominant. Interestingly, we find that the final mass of light BHs and of the SMBH in the quasar is roughly equal by z=6; by the same time only 19% of the initial baryon content has been converted into stars. The SMBH growth is dominated at all epochs z > 7.2 by mergers (exceeding accretion by a factor 2-50); at later times accretion becomes by far the most important growth channel. We finally discuss possible shortcomings of the model.

by Tang, Alf and Sumner, Timothy J.
12 pages, 6 figures

In this paper we demonstrate a methodology to remove the power of the drift induced from random acceleration on LISA proof mass in the frequency domain. The drift must be cleaned from LISA time series data in advance of any further analysis. The cleaning is usually performed in the time domain by using a quadratic function to fit the time series data, and then removing the fitted part from the data. Having Fourier transformed the residuals, and then convolved with LISA transfer function, LISA sensitivity curve can be obtained. However, cosmic gravitational-wave background cannot be retrieved with this approach due to its random nature. Here we provide a new representation of power spectrum given by discrete Fourier transform, which is applied to find the function of the drift power for the cleaning in the frequency domain. We also give the probability distribution used to analyze the data in the frequency domain. We combine several techniques, including Markov Chain Monte Carlo method, simulated annealing, and Gelman & Rubin’s method, with Baye’s theorem to build the algorithm. The algorithm is utilized to analyze 24 simulations of LISA instrumental noise. We prove that the LISA sensitivity can be recovered through this approach. It can help us to build algorithms for some tasks which are must accomplished in the frequency domain for LISA data analysis. This method can be applied to other space-borne interferometers if charges on their proof masses cannot be perfectly cancelled.