A Dynamical Potential-Density Pair for Star Clusters With Nearly Isothermal Interiors

We present a potential-density pair designed to model nearly isothermal star clusters (and similar self-gravitating systems) with a central core and an outer turnover radius, beyond which density falls off as $r^{-4}$. In the intermediate zone, the profile is similar to that of an isothermal sphere (density $\rho \propto r^{-2}$), somewhat less steep than the King 62 profile, and with the advantage that many dynamical quantities can be written in a simple closed form. We derive new analytic expressions for the cluster binding energy and velocity dispersion, and apply these to create toy models for cluster core collapse and evaporation. We fit our projected surface brightness profiles to observed globular and open clusters, and find that the quality of the fit is generally at least as good as that for the surface brightness profiles of King 62. This model can be used for convenient computation of the dynamics and evolution of globular and nuclear star clusters.Comment: 6 pages, 5 figures. Published in ApJL; changes to match published versio

The Tree-Particle-Mesh N-body Gravity Solver

The Tree-Particle-Mesh (TPM) N-body algorithm couples the tree algorithm for directly computing forces on particles in an hierarchical grouping scheme with the extremely efficient mesh based PM structured approach. The combined TPM algorithm takes advantage of the fact that gravitational forces are linear functions of the density field. Thus one can use domain decomposition to break down the density field into many separate high density regions containing a significant fraction of the mass but residing in a very small fraction of the total volume. In each of these high density regions the gravitational potential is computed via the tree algorithm supplemented by tidal forces from the external density distribution. For the bulk of the volume, forces are computed via the PM algorithm; timesteps in this PM component are large compared to individually determined timesteps in the tree regions. Since each tree region can be treated independently, the algorithm lends itself to very efficient parallelization using message passing. We have tested the new TPM algorithm (a refinement of that originated by Xu 1995) by comparison with results from Ferrell & Bertschinger's P^3M code and find that, except in small clusters, the TPM results are at least as accurate as those obtained with the well-established P^3M algorithm, while taking significantly less computing time. Production runs of 10^9 particles indicate that the new code has great scientific potential when used with distributed computing resources.Comment: 24 pages including 9 figures, uses aaspp4.sty; revised to match published versio

Reasoning From Fossils: Learning From the Local Black Hole Population About the Evolution of Quasars

We discuss a simple model for the growth of supermassive black holes (BHs) at the center of spheroidal stellar systems. In particular, we assess the hypotheses that (1) star formation in spheroids and BH fueling are proportional to one another, and (2) the BH accretion luminosity stays near the Eddington limit during luminous quasar phases. With the aid of this simple model, we are able to interpret many properties of the QSO luminosity function, including the puzzling steep decline of the characteristic luminosity from redshift z=2 to to z=0: indeed the residual star formation in spheroidal systems is today limited to a small number of bulges, characterized by stellar velocity dispersions a factor of 2-3 smaller those of the elliptical galaxies hosting QSOs at z > 2. A simple consequence of our hypotheses is that the redshift evolution of the QSO emissivity and of the star formation history in spheroids should be roughly parallel. We find this result to be broadly consistent with our knowledge of the evolution of both the global star formation rate, and of the evolution of the QSO emissivity, but we identify interesting discrepancies at both low and high redshifts, to which we offer tentative solutions. Finally, our hypotheses allow us to present a robust method to derive the duty cycle of QSO activity, based on the observed QSO luminosity function, and on the present-day relation between the masses of supermassive BHs and those of their spheroidal host stellar systems. The duty cycle is found to be substantially less than unity, with characteristic values in the range (3-6)x10^(-3), and we compute that the average bolometric radiative efficiency is epsilon=0.07. Finally, we find that the growth in mass of individual black holes at high redshift (z>2) can be dominated by mergers, and is therefore not necessarily limited by accretion.Comment: Submitted to ApJ, 26 preprint pages with 3 figure