We present a method for fitting orbit-superposition models to the kinematics of discrete stellar systems when the available stellar sample has been filtered by a known selection function. The fitting method can be applied to any model in which the distribution function is represented as a linear superposition of basis elements with unknown weights. As an example, we apply it to Fritz et al.'s kinematics of the innermost regions of the Milky Way's nuclear stellar cluster. Assuming spherical symmetry, our models fit a black hole of mass $M_\bullet=(3.76\pm0.22)\times10^6\,M_\odot$, surrounded by an extended mass $M_\star=(6.57\pm0.54)\times10^6\,M_\odot$ within 4 pc. Within 1 pc the best-fitting mass models have an approximate power-law density cusp $\rho\propto r^{-\gamma}$ with $\gamma=1.3\pm0.3$. We carry out an extensive investigation of how our modelling assumptions might bias these estimates: $M_\bullet$ is the most robust parameter and $\gamma$ the least. Internally the best-fitting models have broadly isotropic orbit distributions, apart from a bias towards circular orbits between 0.1 and 0.3 parsec.

It has been known for more than thirty years that the distribution of molecular gas in the innermost 300 parsecs of the Milky Way, the Central Molecular Zone, is strongly asymmetric. Indeed, approximately three quarters of molecular emission comes from positive longitudes, and only one quarter from negative longitudes. However, despite much theoretical effort, the origin of this asymmetry has remained a mystery. Here we show that the asymmetry can be neatly explained by unsteady flow of gas in a barred potential. We use high-resolution 3D hydrodynamical simulations coupled to a state-of-the-art chemical network. Despite the initial conditions and the bar potential being point-symmetric with respect to the Galactic Centre, asymmetries develop spontaneously due to the combination of a hydrodynamical instability known as the "wiggle instability" and the thermal instability. The observed asymmetry must be transient: observations made tens of megayears in the past or in the future would often show an asymmetry in the opposite sense. Fluctuations of amplitude comparable to the observed asymmetry occur for a large fraction of the time in our simulations, and suggest that the present is not an exceptional moment in the life of our Galaxy.

We run hydrodynamical simulations of a 2D isothermal non self-gravitating inviscid gas flowing in a rigidly rotating externally imposed potential formed by only two components: a monopole and a quadrupole. We explore systematically the effects of varying the quadrupole while keeping fixed the monopole and discuss the consequences for the interpretation of longitude-velocity diagrams in the Milky Way. We find that the gas flow can constrain the quadrupole of the potential and the characteristics of the bar that generates it. The exponential scale length of the bar must be at least $1.5\rm\, kpc$. The strength of the bar is also constrained. Our global interpretation favours a pattern speed of $\Omega=40\,\rm km s^{-1} {kpc}^{-1}$. We find that for most observational features, there exist a value of the parameters that matches each individual feature well, but is difficult to reproduce all the important features at once. Due to the intractably high number of parameters involved in the general problem, quantitative fitting methods that can run automatic searches in parameter space are necessary.

Gaussian processes are the ideal tool for modelling the Galactic ISM, combining statistical flexibility with a good match to the underlying physics. In an earlier paper we outlined how they can be employed to construct three-dimensional maps of dust extinction from stellar surveys. Gaussian processes scale poorly to large datasets though, which put the analysis of realistic catalogues out of reach. Here we show how a novel combination of the Expectation Propagation method and certain sparse matrix approximations can be used to accelerate the dust mapping problem. We demonstrate, using simulated Gaia data, that the resultant algorithm is fast, accurate and precise. Critically, it can be scaled up to map the Gaia catalogue.

We use hydrodynamical simulations to construct a new coherent picture for the gas flow in the Central Molecular Zone (CMZ), the region of our Galaxy within $R\leq 500\, \mathrm{pc}$. We relate connected structures observed in $(l,b,v)$ data cubes of molecular tracers to nuclear spiral arms. These arise naturally in hydrodynamical simulations of barred galaxies, and are similar to those that can be seen in external galaxies such as NGC4303 or NGC1097. We discuss a face-on view of the CMZ including the position of several prominent molecular clouds, such as Sgr B2, the $20\,{\rm km\, s^{-1}}$ and $50\,{\rm km\, s^{-1}}$ clouds, the polar arc, Bania Clump 2 and Sgr C. Our model is also consistent with the larger scale gas flow, up to $R\simeq 3\,\rm kpc$, thus providing a consistent picture of the entire Galactic bar region.

We derive the kinetic equation that describes the secular evolution of a large set of particles orbiting a dominant massive object, such as stars bound to a supermassive black hole or a proto-planetary debris disc encircling a star. Because the particles move in a quasi-Keplerian potential, their orbits can be approximated by ellipses whose orientations remain fixed over many dynamical times. The kinetic equation is obtained by simply averaging the BBGKY equations over the fast angle that describes motion along these ellipses. This so-called Balescu-Lenard equation describes self-consistently the long-term evolution of the distribution of quasi-Keplerian orbits around the central object: it models the diffusion and drift of their actions, induced through their mutual resonant interaction. Hence, it is the master equation that describes the secular effects of resonant relaxation. We show how it captures the phenonema of mass segregation and of the relativistic Schwarzschild barrier recently discovered in $N$-body simulations.

We present our mass estimate of the central black hole in the isolated spiral galaxy NGC 4414. Using natural guide star adaptive optics assisted observations with the Gemini Near-Infrared Integral Field Spectrometer (NIFS) and the natural seeing Gemini Multi-Object Spectrographs-North (GMOS), we derived two-dimensional stellar kinematic maps of NGC 4414 covering the central 1.5 arcsec and 10 arcsec, respectively, at a NIFS spatial resolution of 0.13 arcsec. The kinematic maps reveal a regular rotation pattern and a central velocity dispersion dip down to around 105 km/s. We constructed dynamical methods using two different methods: Jeans anisotropic dynamical modeling and axisymmetric Schwarzschild modeling. Both modeling methods give consistent results, but we cannot constrain the lower mass limit and only measure an upper limit for the black hole mass of Mbh= 1.56 x 10^6 Msun(at 3 sigma level) which is at least 1 sigma below the recent Mbh-sigma_e relations. Further tests with dark matter, mass-to-light ratio variation and different light models confirm that our results are not dominated by uncertainties. The derived upper limit mass is not only below the Mbh-sigma_e relation, but is also five times lower than the lower limit black hole mass anticipated from the resolution limit of the sphere of influence. This proves that via high quality integral field data we are now able to push black hole measurements down to at least five times less than the resolution limit.

The secular evolution of an infinitely thin tepid isolated galactic disc made of a finite number of particles is investigated using the inhomogeneous Balescu-Lenard equation expressed in terms of angle-action variables. The matrix method is implemented numerically in order to model the induced gravitational polarization. Special care is taken to account for the amplification of potential fluctuations of mutually resonant orbits and the unwinding of the induced swing amplified transients. Quantitative comparisons with ${N-}$body simulations yield consistent scalings with the number of particles and with the self-gravity of the disc: the fewer particles and the colder the disc, the faster the secular evolution. Secular evolution is driven by resonances, but does not depend on the initial phases of the disc. For a Mestel disc with ${Q \sim 1.5}$, the polarization cloud around each star boosts up its secular effect by a factor of the order of a thousand or more, promoting accordingly the dynamical relevance of self-induced collisional secular evolution. The position and shape of the induced resonant ridge are found to be in very good agreement with the prediction of the Balescu-Lenard equation, which scales with the square of the susceptibility of the disc. In astrophysics, the inhomogeneous Balescu-Lenard equation may describe the secular diffusion of giant molecular clouds in galactic discs, the secular migration and segregation of planetesimals in proto-planetary discs, or even the long-term evolution of population of stars within the Galactic centre. It could be used as a valuable check of the accuracy of ${N-}$body integrators over secular timescales.

Spiral arms that emerge from the ends of a galactic bar are important in interpreting observations of our and external galaxies. It is therefore important to understand the physical mechanism that causes them. We find that these spiral arms can be understood as kinematic density waves generated by librations around underlying ballistic closed orbits. This is even true in the case of a strong bar, provided the librations are around the appropriate closed orbits and not around the circular orbits that form the basis of the epicycle approximation. An important consequence is that it is a potential's orbital structure that determines whether a bar should be classified as weak or strong, and not crude estimates of the potential's deviation from axisymmetry.

We use a Cartesian grid to simulate the flow of gas in a barred Galactic potential and investigate the effects of varying the sound speed in the gas and the resolution of the grid. For all sound speeds and resolutions, streamlines closely follow closed orbits at large and small radii. At intermediate radii shocks arise and the streamlines shift between two families of closed orbits. The point at which the shocks appear and the streamlines shift between orbit families depends strongly on sound speed and resolution. For sufficiently large values of these two parameters, the transfer happens at the cusped orbit as hypothesised by Binney et al. over two decades ago. For sufficiently high resolutions the flow downstream of the shocks becomes unsteady. If this unsteadiness is physical, as appears to be the case, it provides a promising explanation for the asymmetry in the observed distribution of CO.

We present a new method for fitting simple hydrodynamical models to the (l,v) distribution of atomic and molecular gas observed in the Milky Way. The method works by matching features found in models and observations. It is based on the assumption that the large-scale features seen in (l,v) plots, such as ridgelines and the terminal velocity curve, are influenced primarily by the underlying large-scale Galactic potential and are only weakly dependent on local ISM heating and cooling processes. In our scheme one first identifies by hand the features in the observations: this only has to be done once. We describe a procedure for automatically extracting similar features from simple hydrodynamical models and quantifying the "distance" between each model's features and the observations. Application to models of the Galactic Bar region (|l|<30deg) shows that our feature-fitting method performs better than \chi^2 or envelope distances at identifying the correct underlying galaxy model.

We present a method for obtaining the likelihood function of distance and extinction to a star given its photometry. The other properties of the star (its mass, age, metallicity and so on) are marginalised assuming a simple Galaxy model. We demonstrate that the resulting marginalised likelihood function can be described faithfully and compactly using a Gaussian mixture model. For dust mapping applications we strongly advocate using monochromatic over bandpass extinctions, and provide tables for converting from the former to the latter for different stellar types.

We present a scheme for using stellar catalogues to map the three-dimensional distributions of extinction and dust within our Galaxy. Extinction is modelled as a Gaussian random field, whose covariance function is set by a simple physical model of the ISM that assumes a Kolmogorov-like power spectrum of turbulent fluctuations. As extinction is modelled as a random field, the spatial resolution of the resulting maps is set naturally by the data available; there is no need to impose any spatial binning. We verify the validity of our scheme by testing it on simulated extinction fields and show that its precision is significantly improved over previous dust-mapping efforts. The approach we describe here can make use of any photometric, spectroscopic or astrometric data; it is not limited to any particular survey. Consequently, it can be applied to a wide range of data from both existing and future surveys.

We present a new framework for estimating a galaxy's gravitational potential, Phi(x), from its stellar kinematics by adopting a fully non-parametric model for the galaxy's unknown action-space distribution function, f(J). Having an expression for the joint likelihood of Phi and f, the likelihood of Phi is calculated by using a Dirichlet process mixture to represent the prior on f and marginalising. We demonstrate that modelling machinery constructed using this framework is successful at recovering the potentials of some simple systems from perfect discrete kinematical data, a situation handled effortlessly by traditional moment-based methods, such as the virial theorem, but in which other, more modern, methods are less than satisfactory. We show how to generalise the machinery to account for realistic observational errors and selection functions. A practical implementation is likely to raise some interesting algorithmic and computational challenges.

We present three-dimensional eccentric disc models of the nucleus of M31, modelling the disc as a linear combination of thick rings of massless stars orbiting in the potential of a central black hole. Our models are nonparametric generalisations of the parametric models of Peiris & Tremaine. The models reproduce well the observed WFPC2 photometry, the detailed line-of-sight velocity distributions from STIS observations along P1 and P2, together with the qualitative features of the OASIS kinematic maps. We confirm Peiris & Tremaine's finding that nuclear discs aligned with the larger disc of M31 are strongly ruled out. Our optimal model is inclined at 57 degrees with respect to the line of sight of M31 and has a position angle of 55 degrees. It has a central black hole of mass 10^8 solar masses, and, when viewed in three dimensions, shows a clear enhancement in the density of stars around the black hole. The distribution of orbit eccentricities in our models is similar to Peiris & Tremaine's model, but we find significantly different inclination distributions, which might provide valuable clues to the origin of the disc.

This work aims to study the distribution of luminous and dark matter in Coma early-type galaxies. Dynamical masses obtained under the assumption that mass follows light do not match with the masses of strong gravitational lens systems of similar velocity dispersions. Instead, dynamical fits with dark matter halos are in good agreement with lensing results. We derive mass-to-light ratios of the stellar populations from Lick absorption line indices, reproducing well the observed galaxy colours. Even in dynamical models with dark matter halos the amount of mass that follows the light increases more rapidly with galaxy velocity dispersion than expected for a constant stellar initial mass function (IMF). While galaxies around sigma ~ 200 km/s are consistent with a Kroupa IMF, the same IMF underpredicts luminous dynamical masses of galaxies with sigma ~ 300 km/s by a factor of two and more. A systematic variation of the stellar IMF with galaxy velocity dispersion could explain this trend with a Salpeter IMF for the most massive galaxies. If the IMF is instead constant, then some of the dark matter in high velocity dispersion galaxies must follow a spatial distribution very similar to that of the light. A combination of both, a varying IMF and a component of dark matter that follows the light is possible as well. For a subsample of galaxies with old stellar populations we show that the tilt in the fundamental plane can be explained by systematic variations of the total (stellar + dark) mass inside the effective radius. We tested commonly used mass estimator formulae, finding them accurate at the 20-30% level.

We derive improved versions of the relations between supermassive black hole mass (M_BH) and host-galaxy bulge velocity dispersion (sigma) and luminosity (L) (the M-sigma and M-L relations), based on 49 M_BH measurements and 19 upper limits. Particular attention is paid to recovery of the intrinsic scatter (epsilon_0) in both relations. We find log(M_BH / M_sun) = alpha + beta * log(sigma / 200 km/s) with (alpha, beta, epsilon_0) = (8.12 +/- 0.08, 4.24 +/- 0.41, 0.44 +/- 0.06) for all galaxies and (alpha, beta, epsilon_0) = (8.23 +/- 0.08, 3.96 +/- 0.42, 0.31 +/- 0.06) for ellipticals. The results for ellipticals are consistent with previous studies, but the intrinsic scatter recovered for spirals is significantly larger. The scatter inferred reinforces the need for its consideration when calculating local black hole mass function based on the M-sigma relation, and further implies that there may be substantial selection bias in studies of the evolution of the M-sigma relation. We estimate the M-L relationship as log(M_BH / M_sun) = alpha + beta * log(L_V / 10^11 L_sun,V) of (alpha, beta, epsilon_0) = (8.95 +/- 0.11, 1.11 +/- 0.18, 0.38 +/- 0.09); using only early-type galaxies. These results appear to be insensitive to a wide range of assumptions about the measurement errors and the distribution of intrinsic scatter. We show that culling the sample according to the resolution of the black hole's sphere of influence biases the relations to larger mean masses, larger slopes, and incorrect intrinsic residuals.

We use oblate axisymmetric dynamical models including dark halos to determine the orbital structure of intermediate mass to massive Coma early-type galaxies. We find a large variety of orbital compositions. Averaged over all sample galaxies the unordered stellar kinetic energy in the azimuthal and the radial direction are of the same order, but they can differ by up to 40 percent in individual systems. In contrast, both for rotating and non-rotating galaxies the vertical kinetic energy is on average smaller than in the other two directions. This implies that even most of the rotating ellipticals are flattened by an anisotropy in the stellar velocity dispersions. Using three-integral axisymmetric toy models we show that flattening by stellar anisotropy maximises the entropy for a given density distribution. Collisionless disk merger remnants are radially anisotropic. The apparent lack of strong radial anisotropy in observed early-type galaxies implies that they may not have formed from mergers of disks unless the influence of dissipational processes was significant.

Axisymmetric, orbit-based dynamical models are used to derive dark matter scaling relations for Coma early-type galaxies. From faint to bright galaxies halo core-radii and asymptotic circular velocities increase. Compared to spirals of the same brightness, the majority of Coma early-types -- those with old stellar populations -- have similar halo core-radii but more than 2 times larger asymptotic halo velocities. The average dark matter density inside 2 reff decreases with increasing luminosity and is 6.8 times larger than in disk galaxies of the same B-band luminosity. Compared at the same stellar mass, dark matter densities in ellipticals are 13.5 times higher than in spirals. Different baryon concentrations in ellipticals and spirals cannot explain the higher dark matter density in ellipticals. Instead, the assembly redshift (1+z) of Coma early-type halos is likely about two times larger than of comparably bright spirals. Assuming that local spirals typically assemble at a redshift of one, the majority of bright Coma early-type galaxy halos must have formed around z = 2-3. For about half of our Coma galaxies the assembly redshifts match with constraints derived from stellar populations. We find dark matter densities and estimated assembly redshifts of our observed Coma galaxies in reasonable agreement with recent semi-analytic galaxy formation models.

We determine the mass of the black hole at the center of the spiral galaxy NGC 4258 by constructing axisymmetric dynamical models of the galaxy. These models are constrained by high spatial resolution imaging and long-slit spectroscopy of the nuclear region obtained with the {\em Hubble Space Telescope}, complemented by ground-based observations extending to larger radii. Our best mass estimate is $\MBH = (3.3 \pm 0.2) \times 10^7 \MSun $ for a distance of 7.28 Mpc (statistical errors only). This is within 15% of $ (3.82\pm 0.01) \times 10^7 \MSun$, the mass determined from the kinematics of water masers (rescaled to the same distance) assuming they are in Keplerian rotation in a warped disk. The construction of accurate dynamical models of NGC 4258 is somewhat compromised by an unresolved active nucleus and color gradients, the latter caused by variations in the stellar population and/or obscuring dust. These problems are not present in the $\sim 30$ other black hole mass determinations from stellar dynamics that have been published by us and other groups; thus, the relatively close agreement between the stellar dynamical mass and the maser mass in NGC 4258 enhances our confidence in the black hole masses determined in other galaxies from stellar dynamics using similar methods and data of comparable quality.

We present a general scheme for constructing Monte Carlo realizations of equilibrium, collisionless galaxy models with known distribution function (DF) f_0. Our method uses importance sampling to find the sampling DF f_s that minimizes the mean-square formal errors in a given set of projections of the DF f_0. The result is a multi-mass N-body realization of the galaxy model in which ``interesting'' regions of phase-space are densely populated by lots of low-mass particles, increasing the effective N there, and less interesting regions by fewer, higher-mass particles. As a simple application, we consider the case of minimizing the shot noise in estimates of the acceleration field for an N-body model of a spherical Hernquist model. Models constructed using our scheme easily yield a factor ~100 reduction in the variance in the central acceleration field when compared to a traditional equal-mass model with the same number of particles. When evolving both models with a real N-body code, the diffusion coefficients in our model are reduced by a similar factor. Therefore, for certain types of problems, our scheme is a practical method for reducing the two-body relaxation effects, thereby bringing the N-body simulations closer to the collisionless ideal.

(abridged) Photometry and long-slit spectroscopy are presented for a sample of 6 galaxies with a low surface brightness stellar disc and a bulge. The stellar and ionised-gas kinematics were measured along the major and minor axis in half of the sample galaxies, whereas the other half was observed only along two diagonal axes. Spectra along two diagonal axes were obtained also for one of the objects with major and minor axis spectra. The kinematic measurements extend in the disc region out to a surface-brightness level mu_R~24mag/arcsec^2 reaching in all cases the flat part of the rotation curve. The stellar kinematics turns out to be more regular and symmetric than the ionised-gas kinematics, which often shows the presence of non-circular, off-plane, and non-ordered motions. This raises the question about the reliability of the use of the ionised gas as the tracer of the circular velocity in the modeling of the mass distribution, in particular in the central regions of low surface brightness galaxies.

Supermassive black hole binaries (SMBHBs) are expected by the hierarchical galaxy formation model in $\Lambda$CDM cosmology. There is some evidence in the literature for SMBHBs in AGNs, but there are few observational constraints on the evolution of SMBHBs in inactive galaxies and gas-poor mergers. On the theoretical front, it is unclear how long is needed for a SMBHB in a typical galaxy to coalesce. In this paper we investigate the tidal interaction between stars and binary BHs and calculate the tidal disruption rates of stellar objects by the BH components of binary. We derive the interaction cross sections between SMBHBs and stars from intensive numerical scattering experiments with particle number $\sim10^7$ and calculate the tidal disruption rates by both single and binary BHs for a sample of realistic galaxy models, taking into account the general relativistic effect and the loss cone refilling because of two-body interaction. We estimate the frequency of tidal flares for different types of galaxies using the BH mass function in the literature. We find that because of the three-body slingshot effect, the tidal disruption rate in SMBHB system is more than one order of magnitude smaller than that in single SMBH system. The difference is more significant in less massive galaxies and does not depend on detailed stellar dynamical processes. Our calculations suggest that comparisons of the calculated tidal disruption rates for both single and binary BHs and the surveys of X-ray or UV flares at galactic centers could tell us whether most SMBHs in nearby galaxies are single and whether the SMBHBs formed in gas-poor galaxy mergers coalesce rapidly.

This paper presents a fast, economical particle-multiple-mesh N-body code optimized for large-N modelling of collisionless dynamical processes, such as black-hole wandering or bar-halo interactions, occurring within isolated galaxies. The code has been specially designed to conserve linear momentum. Despite this, it also has variable softening and an efficient block-timestep scheme: the force between any pair of particles is calculated using the finest mesh that encloses them both (respecting Newton's third law) and is updated only on the longest timestep of the two (which conserves momentum). For realistic galaxy models with N > 10^6, it is faster than the fastest comparable momentum-conserving tree code by a factor ranging from ~2 (using single timesteps) to ~10 (multiple timesteps in a concentrated galaxy).

Dynamical models for 17 Coma early-type galaxies are presented. The galaxy sample consists of flattened, rotating as well as non-rotating early-types including cD and S0 galaxies with luminosities between M=-18.79 and M=-22.56. Kinematical long-slit observations cover at least the major and minor axis and extend to 1-4 effective radii. Axisymmetric Schwarzschild models are used to derive stellar mass-to-light ratios and dark halo parameters. In every galaxy models with a dark matter halo match the data better than models without. The statistical significance is over 95 percent for 8 galaxies, around 90 percent for 5 galaxies and for four galaxies it is not significant. For the highly significant cases systematic deviations between observed and modelled kinematics are clearly seen; for the remaining galaxies differences are more statistical in nature. Best-fit models contain 10-50 percent dark matter inside the half-light radius. The central dark matter density is at least one order of magnitude lower than the luminous mass density. The central phase-space density of dark matter is often orders of magnitude lower than in the luminous component, especially when the halo core radius is large. The orbital system of the stars along the major-axis is slightly dominated by radial motions. Some galaxies show tangential anisotropy along the minor-axis, which is correlated with the minor-axis Gauss-Hermite coefficient H4. Changing the balance between data-fit and regularisation constraints does not change the reconstructed mass structure significantly. Model anisotropies tend to strengthen if the weight on regularisation is reduced, but the general property of a galaxy to be radially or tangentially anisotropic, respectively, does not change. (abridged)

Black hole masses predicted from the Mbh-sigma relationship conflict with those predicted from the Mbh-L relationship for the most luminous galaxies, such as brightest cluster galaxies (BCGs). This is because stellar velocity dispersion, sigma, increases only weakly with L for BCGs and other giant ellipticals. The Mbh-L relationship predicts that the most luminous BCGs may have Mbh approaching 10^{10}M_sol, while the M-sigma relationship always predicts Mbh<3X10^9M_sol. We argue that the Mbh-L relationship is a plausible or even preferred description for BCGs and other galaxies of similar luminosity. If cores in central stellar density are formed by binary BHs, then the inner-core cusp radius, r_gamma, may be an independent witness of Mbh. Using structural parameters derived from a large sample of early-type galaxies observed by HST, we argue that L is superior to sigma as an indicator of r_gamma in luminous galaxies. The observed r_gamma-Mbh relationship for 11 core galaxies with measured Mbh appears to be consistent with the Mbh-L relationship for BCGs. BCGs have large cores appropriate for their large luminosities that may be difficult to generate with the modest BH masses inferred from the Mbh-sigma relationship. Mbh~L may be expected to hold for BCGs, if they were formed in dissipationless mergers, which should preserve ratio of BH to stellar mass. This picture appears to be consistent with the slow increase in sigma with L and the more rapid increase in effective radii with L seen in BCGs. If BCGs have large BHs commensurate with their luminosities, then the local BH mass function for Mbh>3X10^9M_sol may be nearly an order of magnitude richer than that inferred from the Mbh-sigma relationship. The volume density of QSOs at earlier epochs may favor the predictions from the Mbh-L relationship.

When faced with the task of constraining a galaxy's potential given limited stellar kinematical information, what is the best way of treating the galaxy's unknown distribution function (DF)? Using the example of estimating black hole (BH) masses, I argue that the correct approach is to consider all possible DFs for each trial potential, marginalizing the DF using an infinitely divisible prior. Alternative approaches, such as the widely used maximum penalized likelihood method, neglect the huge degeneracies inherent in the problem and simply identify a single, special DF for each trial potential. Using simulated observations of toy galaxies with realistic amounts of noise, I find that this marginalization procedure yields significantly tighter constraints on BH masses than the conventional maximum-likelihood method, although it does pose a computational challenge which might be solved with the development of a suitable algorithm for massively parallel machines. I show that in practice the conventional maximum-likelihood method yields reliable BH masses with well-defined minima in their chi^2 distributions, contrary to claims made by Valluri, Merritt & Emsellem.

We present near infra-red (NIR) adaptive optics assisted spectroscopic observations of the CO ($\Delta\mu=2$) absorption bands towards the centre of the giant elliptical galaxy NGC 1399. The observations were made with NAOS-CONICA (ESO VLT) and have a FWHM resolution of 0\farcs15 (14pc). Kinematic analysis of the observations reveals a decoupled core and strongly non-Gaussian line-of-sight velocity profiles (VPs) in the central 0.2 arcsec (19pc). NIR imaging also indicates an asymmetric elongation of the central isophotes in the same region. We use spherical orbit-superposition models to interpret the kinematics, using a set of orthogonal ``eigenVPs'' that allow us to fit models directly to spectra. The models require a central black hole of mass $1.2^{+0.5}_{-0.6}\times10^9M_\odot$, with a strongly tangentially biased orbit distribution in the inner 40pc.

This is the second in a series of papers dedicated to unveil the mass structure and orbital content of a sample of flattened early-type galaxies in the Coma cluster. The ability of our orbit libraries to reconstruct internal stellar motions and the mass composition of a typical elliptical in the sample is investigated by means of Monte-Carlo simulations of isotropic rotator models. The simulations allow a determination of the optimal amount of regularization needed in the orbit superpositions. It is shown that under realistic observational conditions and with the appropriate regularization internal velocity moments can be reconstructed to an accuracy of about 15 per cent; the same accuracy can be achieved for the circular velocity and dark matter fraction. In contrast, the flattening of the halo remains unconstrained. Regularized orbit superpositions are applied to a first galaxy in our sample, NGC 4807, for which stellar kinematical observations extend to 3 Reff. The galaxy seems dark matter dominated outside 2 Reff. Logarithmic dark matter potentials are consistent with the data, as well as NFW-profiles, mimicking logarithmic potentials over the observationally sampled radial range. In both cases, the derived stellar mass-to-light ratio agrees well with independently obtained mass-to-light ratios from stellar population analysis. Kinematically, NGC 4807 is characterized by mild radial anisotropy outside r>0.5 Reff, becoming isotropic towards the center. Our orbit models hint at either a distinct stellar component or weak triaxiality in the outer parts of the galaxy.

We analyze HST+WFPC2 images of 77 early-type galaxies. Brightness profiles are classed into "core" or "power-law" forms. Cores are typically rounder than power-law galaxies. Nearly all power-laws with central ellipticity >=0.3 have stellar disks, implying that disks are present in power-laws with epsilon <0.3, but are not visible due to unfavorable geometry. A few low-luminosity core galaxies also have disks; these may be transition forms from power-laws. Cores and power-laws both have twisting isophotes at small radii. Core galaxies have somewhat weaker color gradients than power-laws. Nuclei are found in 29% of the cores and 60% of the power-laws. Nuclei are typically bluer than the surrounding galaxy. NGC 4073 and 4382 have central minima in their intrinsic starlight distributions; NGC 4382 resembles the double nucleus of M31. In general, the peak brightness location is coincident with the photocenter of the core to <1 pc. Five galaxies, however, have centers significantly displaced from their cores; these may be unresolved asymmetric double nuclei. Central dust is visible in half of the galaxies. The presence and strength of dust correlates with nuclear emission. The prevalence of dust and its morphology suggests that dust clouds form, settle to the center, and disappear repeatedly on ~10^8 yr timescales. We discuss the hypothesis that cores are created by the decay of a massive black hole binary. Apart from their brightness profiles, there are no strong differences between cores and power-laws that demand this scenario; however, the rounder shapes of cores, their lack of disks, and their reduced color gradients may be consistent with it.

This is the first of a series of papers dedicated to unveil the mass composition and dynamical structure of a sample of flattened early type galaxies in the Coma cluster. We describe our modifications to the Schwarzschild code of Richstone et al. (in preparation). Applying a Voronoi tessellation in the surface of section we are able to assign accurate phase-space volumes to individual orbits and to reconstruct the full three-integral phase-space distribution function (DF) of any axisymmetric orbit library. Two types of tests have been performed to check the accuracy with which DFs can be represented by appropriate orbit libraries. First, by mapping DFs of spherical gamma-models and flattened Plummer models on to the library we show that the resulting line-of-sight velocity distributions and internal velocity moments of the library match those directly derived from the DF to a precision better than that of present day observational errors. Second, by fitting libraries to the projected kinematics of the same DFs we show that the distribution function reconstructed from the fitted library matches the input DF to a rms of about 15 per cent over a region in phase-space covering 90 per cent of the mass of the library. The achieved accuracy allows us to implement effective entropy-based regularisation to fit real, noisy and spatially incomplete data.

We show that orbit-superposition dynamical models (Schwarzschild's method) provide reliable estimates of nuclear black hole masses and errors when constructed from adequate orbit libraries and kinematic data. We thus rebut two recent papers that argue that BH masses obtained from this method are unreliable. These papers claim to demonstrate that the range of allowable BH masses derived from a given dataset is artificially too narrow as a result of an inadequate number of orbits in the library used to construct dynamical models. This is an elementary error that is easily avoided. We describe a method to estimate the number and nature of orbits needed for the library. We provide an example that shows that this prescription is adequate, in the sense that the range of allowable BH masses is not artificially narrowed by use of too few orbits. A second point raised by critics is that kinematic data are generally obtained with insufficient spatial resolution to obtain a reliable mass. We make the distinction between unreliable determinations and imprecise ones. We show that there are several different properties of a kinematic dataset that can lead to imprecise BH determinations, but none of the attributes we have investigated leads to an unreliable determination. In short, the degree to which the BH radius of influence is resolved by spectroscopic observations is already reflected in the BH-mass error envelope, and is not a hidden source of error. The BH masses published by our group and the Leiden group are reliable.

We describe an algorithm for constructing N-body realizations of equilibrium spherical systems. A general form for the mass density rho(r) is used, making it possible to represent most of the popular density profiles found in the literature, including the cuspy density profiles found in high-resolution cosmological simulations. We demonstrate explicitly that our models are in equilibrium. In contrast, many existing N-body realizations of isolated systems have been constructed under the assumption that the local velocity distribution is Maxwellian. We show that a Maxwellian halo with an initial r^{-1} central density cusp immediately develops a constant-density core. Moreover, after just one crossing time the orbital anisotropy has changed over the entire system, and the initially isotropic model becomes radially anisotropic. These effects have important implications for many studies, including the survival of substructure in cold dark matter (CDM) models. Comparing the evolution and mass-loss rate of isotropic Maxwellian and self-consistent Navarro, Frenk, & White (NFW) satellites orbiting inside a static host CDM potential, we find that the former are unrealistically susceptible to tidal disruption. Thus, recent studies of the mass-loss rate and disruption timescales of substructure in CDM models may be compromized by using the Maxwellian approximation. We also demonstrate that a radially anisotropic, self-consistent NFW satellite loses mass at a rate several times higher than that of its isotropic counterpart on the same external tidal field and orbit.

We present preliminary results from spectroscopic observations of a sample of 11 low surface brightness galaxies (LSB). We measured the stellar and gaseous kinematics along their major and minor axes. Such information will allow us to accurately investigate the dark matter (DM) content within their optical regions, providing further constraints on the predictions of standard CDM models. Although dynamical modeling is still in progress, our observations already show that the derived stellar kinematics is more regular than the ionized gas one, which often shows evidence for non-circular and asymmetric motions.

We present stellar kinematics for a sample of 10 early-type galaxies observed using the STIS aboard the Hubble Space Telescope, and the Modular Spectrograph on the MDM Observatory 2.4-m telescope. The spectra are used to derive line-of-sight velocity distributions (LOSVDs) of the stars using a Maximum Penalized Likelihood method. We use Gauss-Hermite polynomials to parameterize the LOSVDs and find predominantly negative h4 values (boxy distributions) in the central regions of our galaxies. One galaxy, NGC 4697, has significantly positive central h4 (high tail weight). The majority of galaxies have a central velocity dispersion excess in the STIS kinematics over ground-based velocity dispersions. The galaxies with the strongest rotational support, as quantified with v_MAX/sigma_STIS, have the smallest dispersion excess at STIS resolution. The best-fitting, general, axisymmetric dynamical models (described in a companion paper) require black holes in all cases, with masses ranging from 10^6.5 to 10^9.3 Msun. We replot these updated masses on the BH/sigma relation, and show that the fit to only these 10 galaxies has a slope consistent with the fits to larger samples. The greatest outlier is NGC 2778, a dwarf elliptical with relatively poorly constrained black hole mass. The two best candidates for pseudobulges, NGC 3384 and 7457, do not deviate significantly from the established relation between black hole and sigma. Neither do the three galaxies which show the most evidence of a recent merger, NGC 3608, 4473, and 4697.

We have investigated the dynamics of the merging process in the minor merger hypothesis for active galactic nuclei. We find that for a satellite galaxy to be able to merge directly with the nucleus of the host galaxy (for example, to give rise to the compact dust discs which are seen in early type active galaxies) requires the initial orbit of the satellite to be well aimed. For the case of the host galaxy being a disc galaxy, if the initial orbits of the satellites are randomly oriented with respect to the host galaxy, then the orbits of those which reach the host nuclear regions in a reasonable time, are also fairly randomly oriented once they reach the nucleus. We note that this result might be able to provide an explanation of why the jet directions in the nuclei of Seyfert galaxies are apparently unrelated to the plane of the galaxy discs.

We present axisymmetric, orbit superposition models for 12 galaxies using data taken with the Hubble Space Telescope (HST) and ground-based observatories. In each galaxy, we detect a central black hole (BH) and measure its mass to accuracies ranging from 10% to 70%. We demonstrate that in most cases the BH detection requires both the HST and ground-based data. Using the ground-based data alone does provide an unbiased measure of the BH mass (provided they are fit with fully general models), but at a greatly reduced significance. The most significant correlation with host galaxy properties is the relation between the BH mass and the velocity dispersion of the host galaxy; we find no other equally strong correlation, and no second parameter that improves the quality of the mass-dispersion relation. We are also able to measure the stellar orbital properties from these general models. The most massive galaxies are strongly biased to tangential orbits near the BH, consistent with binary BH models, while lower-mass galaxies have a range of anisotropies, consistent with an adiabatic growth of the BH.

We used HST WFPC2 images to identify six early-type galaxies with surface- brightness profiles that decrease inward over a limited range of radii near their centers. The implied luminosity density profiles of these galaxies have local minima interior to their core break radii. NGC 3706 harbors a high surface brightness ring of starlight with radius ~20 pc. Its central structure may be related to that in the double-nucleus galaxies M31 and NGC 4486B. NGC 4406 and NGC 6876 have nearly flat cores that on close inspection are centrally depressed. Colors for both galaxies imply that this is not due to dust absorption. The surface brightness distributions of both galaxies are consistent with stellar tori that are more diffuse than the sharply defined system in NGC 3706. The remaining three galaxies are the brightest cluster galaxies in A260, A347, and A3574. Color information is not available for these objects, but they strongly resemble NGC 4406 and NGC 6876 in their cores. The thin ring in NGC 3706 may have formed dissipatively. The five other galaxies resemble the endpoints of some simulations of the merging of two gas-free stellar systems, each harboring a massive nuclear black hole. In one version of this scenario, diffuse stellar tori are produced when stars initially bound to one black hole are tidally stripped away by the second black hole. Alternatively, some inward-decreasing surface-brightness profiles may reflect the ejection of stars from a core during the hardening of the binary black hole created during the merger.

Observations of nearby galaxies reveal a strong correlation between the mass of the central dark object M and the velocity dispersion sigma of the host galaxy, of the form log(M/M_sun) = a + b*log(sigma/sigma_0); however, published estimates of the slope b span a wide range (3.75 to 5.3). Merritt & Ferrarese have argued that low slopes (<4) arise because of neglect of random measurement errors in the dispersions and an incorrect choice for the dispersion of the Milky Way Galaxy. We show that these explanations account for at most a small part of the slope range. Instead, the range of slopes arises mostly because of systematic differences in the velocity dispersions used by different groups for the same galaxies. The origin of these differences remains unclear, but we suggest that one significant component of the difference results from Ferrarese & Merritt's extrapolation of central velocity dispersions to r_e/8 (r_e is the effective radius) using an empirical formula. Another component may arise from dispersion-dependent systematic errors in the measurements. A new determination of the slope using 31 galaxies yields b=4.02 +/- 0.32, a=8.13 +/- 0.06, for sigma_0=200 km/s. The M-sigma relation has an intrinsic dispersion in log M that is no larger than 0.3 dex. In an Appendix, we present a simple model for the velocity-dispersion profile of the Galactic bulge.

We discuss the problem of using stellar kinematics of early-type galaxies to constrain the galaxies' orbital anisotropies and radial mass profiles. We demonstrate that compressing a galaxy's light distribution along the line of sight produces approximately the same signature in the line-of-sight velocity profiles as radial anisotropy. In particular, fitting spherically symmetric dynamical models to apparently round, isotropic face-on flattened galaxies leads to a spurious bias towards radial orbits in the models, especially if the galaxy has a weak face-on stellar disk. Such face-on stellar disks could plausibly be the cause of the radial anisotropy found in spherical models of intermediate luminosity ellipticals such as NGC 2434, NGC 3379 and NGC 6703. In the light of this result, we use simple dynamical models to constrain the outer mass profiles of a sample of 18 round, early-type galaxies. The galaxies follow a Tully-Fisher relation parallel to that for spiral galaxies, but fainter by at least 0.8 mag (I-band) for a given mass. The most luminous galaxies show clear evidence for the presence of a massive dark halo, but the case for dark haloes in fainter galaxies is more ambiguous. We discuss the observations that would be required to resolve this ambiguity.

We analyze the nuclear stellar dynamics of the SB0 galaxy NGC 1023, utilizing observational data both from the Space Telescope Imaging Spectrograph aboard the Hubble Space Telescope and from the ground. The stellar kinematics measured from these long-slit spectra show rapid rotation (V = 70 km/s at a distance of 0.1 arcsec = 4.9 pc from the nucleus) and increasing velocity dispersion toward the nucleus (where sigma = 295 +/- 30 km/s). We model the observed stellar kinematics assuming an axisymmetric mass distribution with both two and three integrals of motion. Both modeling techniques point to the presence of a central dark compact mass (which presumably is a supermassive black hole) with confidence > 99%. The isotropic two-integral models yield a best-fitting black hole mass of (6.0 +/- 1.4) x 10^7 M_sun and mass-to-light ratio (M/L_V) of 5.38 +/- 0.08, and the goodness-of-fit (chi^2) is insensitive to reasonable values for the galaxy's inclination. The three-integral models, which non-parametrically fit the observed line-of-sight velocity distribution as a function of position in the galaxy, suggest a black hole mass of (3.9 +/- 0.4) x 10^7 M_sun and M/L_V of 5.56 +/- 0.02 (internal errors), and the edge-on models are vastly superior fits over models at other inclinations. The internal dynamics in NGC 1023 as suggested by our best-fit three-integral model shows that the velocity distribution function at the nucleus is tangentially anisotropic, suggesting the presence of a nuclear stellar disk. The nuclear line of sight velocity distribution has enhanced wings at velocities >= 600 km/s from systemic, suggesting that perhaps we have detected a group of stars very close to the central dark mass.

Black hole (BH) masses that have been measured by reverberation mapping in active galaxies fall significantly below the correlation between bulge luminosity and BH mass determined from spatially resolved kinematics of nearby normal galaxies. This discrepancy has created concern that one or both techniques suffer from systematic errors. We show that BH masses from reverberation mapping are consistent with the recently discovered relationship between BH mass and galaxy velocity dispersion. Therefore the bulge luminosities are the probable source of the disagreement, not problems with either method of mass measurement. This result underscores the utility of the BH mass -- velocity dispersion relationship. Reverberation mapping can now be applied with increased confidence to galaxies whose active nuclei are too bright or whose distances are too large for BH searches based on spatially resolved kinematics.

We describe a correlation between the mass M_BH of a galaxy's central black hole and the luminosity-weighted line-of-sight velocity dispersion sigma_e within the half-light radius. The result is based on a sample of 26 galaxies, including 13 galaxies with new determinations of black hole masses from Hubble Space Telescope measurements of stellar kinematics. The best-fit correlation is M_BH = 1.2 (+-0.2) x 10^8 M_sun (sigma_e/200 km/s)^(3.75 (+-0.3))over almost three orders of magnitude in M_BH; the scatter in M_BH at fixed sigma_e is only 0.30 dex and most of this is due to observational errors. The M_BH-sigma_e relation is of interest not only for its strong predictive power but also because it implies that central black hole mass is constrained by and closely related to properties of the host galaxy's bulge.

In the context of the CNOC1 cluster survey, redshifts were obtained for galaxies in 16 clusters. The resulting sample is ideally suited for an analysis of the internal velocity and mass distribution of clusters. Previous analyses of this dataset used the Jeans equation to model the projected velocity dispersion profile. However, the results of such an analysis always yield a strong degeneracy between the mass density profile and the velocity dispersion anisotropy profile. Here we analyze the full (R,v) dataset of galaxy positions and velocities in an attempt to break this degeneracy. We build an `ensemble cluster' from the individual clusters under the assumption that they form a homologous sequence. To interpret the data we study a one-parameter family of spherical models with different constant velocity dispersion anisotropy. The best-fit model is sought using a variety of statistics, including the overall likelihood of the dataset. Although the results of our analysis depend slightly on which statistic is used to judge the models, all statistics agree that the best-fit model is close to isotropic. This result derives primarily from the fact that the observed grand-total velocity histogram is close to Gaussian, which is not expected to be the case for a strongly anisotropic model. The best-fitting models have a mass-to-number-density ratio that is approximately independent of radius over the range constrained by the data. They also have a mass-density profile that is consistent with the dark matter halo profile advocated by Navarro, Frenk & White, in terms of both the profile shape and the characteristic scale length. This adds important new weight to the evidence that clusters do indeed follow this proposed universal mass density profile. [Abridged]

We fit axisymmetric 3-integral dynamical models to NGC3379 using the line-of-sight velocity distribution obtained from HST/FOS spectra of the galaxy center and ground-based long-slit spectroscopy along four position angles, with the light distribution constrained by WFPC2 and ground-based images. We have fitted models with inclinations from 29 (intrinsic galaxy type E5) to 90 degrees (intrinsic E1) and black hole masses from 0 to 1e9 M_solar. The best-fit black hole masses range from 6e7 to 2e8 M_solar, depending on inclination. The velocity ellipsoid of the best model is not consistent with either isotropy or a two-integral distribution function. Along the major axis, the velocity ellipsoid becomes tangential at the innermost bin, radial in the mid-range radii, and tangential again at the outermost bins. For the acceptable models, the radial to tangential dispersion in the mid-range radii ranges from 1.1 < sigma_r / sigma_t < 1.7. Compared with these 3-integral models, 2-integral isotropic models overestimate the black hole mass since they cannot provide adequate radial motion. However, the models presented in this paper still contain restrictive assumptions-namely assumptions of constant M/L and spheroidal symmetry-requiring yet more models to study black hole properties in complete generality.

There is strong evidence for some kind of massive dark object in the centres of many galaxy bulges. The detection of flares from tidally disrupted stars could confirm that these objects are black holes (BHs). Here we present calculations of the stellar disruption rates in detailed dynamical models of real galaxies, taking into account the refilling of the loss cone of stars on disruptable orbits by two-body relaxation and tidal forces in non-spherical galaxies. The highest disruption rates (one star per 10^4 yr) occur in faint (L>10^10 L_sun) galaxies, which have steep central density cusps. More luminous galaxies are less dense and have much longer relaxation times and more massive BHs. Dwarf stars in such galaxies are swallowed whole by the BH and hence do not emit flares; giant stars could produce flares as often as every 10^5 yr, although the rate depends sensitively on the shape of the stellar distribution function. We discuss the possibility of detecting disruption flares in current supernova searches. The total mass of stars consumed over the lifetime of the galaxy is of order 10^6 M_sun, independent of galaxy luminosity; thus disrupted stars may contribute significantly to the present BH mass in galaxies fainter than about 10^9 L_sun.

The deprojection of the surface brightness distribution of an axisymmetric galaxy does not have a unique solution unless the galaxy is viewed precisely edge-on. I present an algorithm that finds the full range of smooth axisymmetric density distributions consistent with a given surface brightness distribution and inclination angle, and use it to investigate the effects of this non-uniqueness on the line-of-sight velocity profiles (VPs) of two-integral models of both real and toy disky galaxies viewed at a range of inclination angles. Photometrically invisible face-on disks leave very clear signatures in the minor-axis VPs of the models (Gauss--Hermite coefficients h_4>0.1), provided the disk-to-bulge ratio is greater than about 3%. I discuss the implications of these hitherto neglected disks for dynamical modelling.

Black holes, an extreme consequence of the mathematics of General Relativity, have long been suspected of being the prime movers of quasars, which emit more energy than any other objects in the Universe. Recent evidence indicates that supermassive black holes, which are probably quasar remnants, reside at the centers of most galaxies. As our knowledge of the demographics of these relics of a violent earlier Universe improve, we see tantalizing clues that they participated intimately in the formation of galaxies and have strongly influenced their present-day structure.

We construct dynamical models for a sample of 36 nearby galaxies with Hubble Space Telescope photometry and ground-based kinematics. The models assume that each galaxy is axisymmetric, with a two-integral distribution function, arbitrary inclination angle, a position-independent stellar mass-to-light ratio Upsilon, and a central massive dark object (MDO) of arbitrary mass M_bh. They provide acceptable fits to 32 of the galaxies for some value of M_bh and Upsilon; the four galaxies that cannot be fit have kinematically decoupled cores. The mass-to-light ratios inferred for the 32 well-fit galaxies are consistent with the fundamental plane correlation Upsilon \propto L^0.2, where L is galaxy luminosity. In all but six galaxies the models require at the 95% confidence level an MDO of mass M_bh ~ 0.006 M_bulge = 0.006 Upsilon L. Five of the six galaxies consistent with M_bh=0 are also consistent with this correlation. The other (NGC 7332) has a much stronger upper limit on M_bh. We consider various parameterizations for the probability distribution describing the correlation of the masses of these MDOs with other galaxy properties. One of the best models can be summarized thus: a fraction f ~0.97 of galaxies have MDOs, whose masses are well described by a Gaussian distribution in log (M_bh/M_bulge) of mean -2.27 and width ~0.07.

The stellar kinematics of the dwarf elliptical galaxy NGC 4486B have been measured in seeing sigma_* = .22 arcsec with the Canada-France-Hawaii Telescope. Lauer et al. 1996, ApJ, 471, L79 have shown that NGC 4486B is similar to M31 in having a double nucleus. We show that it also resembles M31 in its kinematics. The velocity dispersion gradient is very steep: sigma increases from 116 +- 6 km/s at r = 2" - 6" to 281 +- 11 km/s at the center. This is much higher than expected for an elliptical galaxy of absolute magnitude M_B = -16.8: NGC 4486B is far above the scatter in the Faber-Jackson correlation between sigma and bulge luminosity. Therefore the King core mass-to-light ratio, M/L_V = 20, is unusually high compared with normal values for old stellar populations. We construct dynamical models with isotropic velocity dispersions and show that they reproduce black hole (BH) masses derived by more detailed methods. We also fit axisymmetric, three-integral models. Isotropic models imply that NGC 4486B contains a central dark object, probably a BH, of mass M_BH = 6^{+3}_{-2} x 10^8 M_sun. However, anisotropic models fit the data without a BH if the ratio of radial to azimuthal dispersions is ~ 2 at 1". Therefore this is a less strong BH detection than the ones in M31, M32, and NGC 3115. A 6 x 10^8 M_sun BH is 9 % of the mass M_bulge in stars; even if M_BH is smaller than the isotropic value, M_BH/M_bulge is likely to be unusually large. Double nuclei are a puzzle because the dynamical friction timescales for self-gravitating star clusters in orbit around each other are short. Since both M31 and NGC 4486B contain central dark objects, our results support models in which the survival of double nuclei is connected with the presence of a BH (e. g., Tremaine 1995, AJ, 110, 628).