Search

The nuclear stellar disc (NSD) is a flattened high-density stellar structure that dominates the gravitational field of the Milky Way at Galactocentric radius 30pc≲R≲300pc⁠. We construct axisymmetric self-consistent equilibrium dynamical models of the NSD in which the distribution function is an analytic function of the action variables. We fit the models to the normalized kinematic distributions (line-of-sight velocities + VIRAC2 proper motions) of stars in the NSD survey of Fritz et al., taking the foreground contamination due to the Galactic Bar explicitly into account using an N-body model. The posterior marginalized probability distributions give a total mass of MNSD=10.5+1.1−1.0×108M⊙⁠, roughly exponential radial and vertical scale lengths of Rdisc=88.6+9.2−6.9pc and Hdisc=28.4+5.5−5.5pc⁠, respectively, and a velocity dispersion σ≃70kms−1 that decreases with radius. We find that the assumption that the NSD is axisymmetric provides a good representation of the data. We quantify contamination from the Galactic Bar in the sample, which is substantial in most observed fields. Our models provide the full 6D (position + velocity) distribution function of the NSD, which can be used to generate predictions for future surveys. We make the models publicly available as part of the software package AGAMA. © 2022 The Author(s) Published by Oxford University Press on behalf of Royal Astronomical Society., MCS acknowledges financial support from the European Research Council via the ERC Synergy Grant ‘ECOGAL – Understanding our Galactic ecosystem: from the disk of the Milky Way to the formation sites of stars and planets’ (grant 855130). JLS thanks the support of the Royal Society (URF\R1\191555). EV acknowledges support from STFC via the consolidated grant to the Institute of Astronomy. DM acknowledges support from project CATA FB210003. AMB acknowledges funding from the EU Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie grant agreement No 895174. JM acknowledges support from the UK Science and Technology Facilities Council under grant number ST/S000488/I. MCS, NN, FN-L, and RSK gratefully acknowledge support by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) - Project-ID 138713538 - SFB 881 (‘The Milky Way System’, subprojects B1, B2, B8). FN-L acknowledges the sponsorship provided by the Federal Ministry for Education and Research of Germany through the Alexander von Humboldt Foundation. RS and BS acknowledge financial support from the State Agency for Research of the Spanish MCIU through the ‘Center of Excellence Severo Ochoa’ award for the Instituto de Astrofísica de Andalucía (SEV-2017-0709) and financial support from national project PGC2018-095049-B-C21 (MCIU/AEI/FEDER, UE). RSK acknowledges support from the Heidelberg cluster of excellence EXC 2181 (Project-ID 390900948) ‘STRUCTURES: A unifying approach to emergent phenomena in the physical world, mathematics, and complex data’ funded by the German Excellence Strategy. This work used computing resources provided by the state of Baden-Württemberg through bwHPC and the German Research Foundation (DFG) through grant INST 35/1134-1 FUGG. Data are stored at SDS@hd supported by the Ministry of Science, Research and the Arts Baden-Württemberg (MWK) and DFG through grant INST 35/1314-1 FUGG.

The nuclear stellar disc (NSD) is a flattened stellar structure that dominates the gravitational potential of the Milky Way at Galactocentric radii $30 \lesssim R \lesssim 300{\, \rm pc}$. In this paper, we construct axisymmetric Jeans dynamical models of the NSD based on previous photometric studies and we fit them to line-of-sight kinematic data of APOGEE and SiO maser stars. We find that (i) the NSD mass is lower but consistent with the mass independently determined from photometry by Launhardt et al. (2002). Our fiducial model has a mass contained within spherical radius $r=100{\, \rm pc}$ of $M(r1$. Observations and theoretical models of the star-forming molecular gas in the central molecular zone suggest that large vertical oscillations may be already imprinted at stellar birth. However, the finding $\sigma_z/\sigma_R > 1$ depends on a drop in the velocity dispersion in the innermost few tens of parsecs, on our assumption that the NSD is axisymmetric, and that the available (extinction corrected) stellar samples broadly trace the underlying light and mass distributions, all of which need to be established by future observations and/or modelling. (iii) We provide the most accurate rotation curve to date for the innermost $500 {\, \rm pc}$ of our Galaxy., Comment: Accepted for publication in MNRAS

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., Comment: 17 pages, MNRAS accepted

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., Comment: 15 pages, 6 figures, submitted for publication in MNRAS

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., Comment: Accepted for publication in MNRAS. Videos of the simulations are available at http://www.ita.uni-heidelberg.de/~mattia/download.html

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., Accepted for publication in MNRAS

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, Accepted for publication in MNRAS

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., Comment: 16 pages, 12 figures. Submitted for publication in MNRAS. Text revised

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., Comment: Accepted for publication in MNRAS, 12 pages, 12 figures

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., Comment: 19 pages, 7 figures. Updated to match accepted version

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., 14 pages, 5 figures, 1 landscape table. Submitted to MNRAS

(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., Accepted for pubblication on MNRAS (11 figures and 15 tables)

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., Comment: 24 pages, 3 figures

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., Accepted for publication in MNRAS

We present observations of 77 early-type galaxies imaged with the PC1 CCD of the Hubble Space Telescope (HST) WFPC2. "Nuker-law" parametric fits to the surface brightness profiles are used to classify the central structure into "core" or "power-law" forms. Core galaxies are typically rounder than power-law galaxies. Nearly all power-law galaxies with central ellipticities ε ≥ 0.3 have stellar disks, implying that disks are present in power-law galaxies with ε < 0.3 but are not visible because of unfavorable geometry. A few low-luminosity flattened core galaxies also have disks; these may be transition forms from power-law galaxies to more luminous core galaxies, which lack disks. Several core galaxies have strong isophote twists interior to their break radii, although power-law galaxies have interior twists of similar physical significance when the photometric perturbations implied by the twists are evaluated. Central color gradients are typically consistent with the envelope gradients; core galaxies have somewhat weaker color gradients than power-law galaxies. Nuclei are found in 29% of the core galaxies and 60% of the power-law galaxies. Nuclei are typically bluer than the surrounding galaxy. While some nuclei are associated with active galactic nuclei (AGNs), just as many are not; conversely, not all galaxies known to have a low-level AGN exhibit detectable nuclei in the broadband filters. NGC 4073 and 4382 are found to 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 a typical physical scale of < 1 pc. Five galaxies, however, have centers significantly displaced from their surrounding cores; these may be unresolved asymmetric double nuclei. Finally, as noted by previous authors, central dust is visible in about half of the galaxies. The presence and strength of dust correlates with nuclear emission; thus, dust may outline gas that is falling into the central black hole. The prevalence of dust and its morphology suggest that dust clouds form, settle to the center, and disappear repeatedly on ∼ 108 yr timescales. We discuss the hypothesis that cores are created by the decay of a massive black hole binary formed in a merger. Apart from their brightness profiles, there are no strong differences between core galaxies and power-law galaxies that demand this scenario; however, the rounder shapes of core, their lack of disks, and their reduced color gradients may be consistent with it. © 2005. The American Astronomical Society. All rights reserved.

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., Comment: accepted for publication in MNRAS. This version corrects a typo in the normalization factor that is present in the published version

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 (, 37 pages, 9 figures

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., Comment: 4 pages, 1 figure, ApJ Letters accepted, minor revisions

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., Comment: 4 pages, accepted ApJ Letters, figure added and minor corrections

We fit axisymmetric three-integral dynamical models to NGC 3379 using the line-of-sight velocity distribution obtained from Hubble Space Telescope 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° (intrinsic E1) and black hole masses from 0 to 109 M⊙. The best-fit black hole masses range from 6 × 107 to 2 × 108 M⊙, depending on inclination. The preferred inclination is 90° (edge-on); however, the constraints on allowed inclination are not very strong, owing to our assumption of constant M/LV. 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 midrange radii, and tangential again at the outermost bins. The rotation rises quickly at small radii owing to the presence of the black hole. For the acceptable models, the radial-to-tangential [(σ2θ + σ2φ)/2] dispersion in the midrange radii ranges over 1.1 < σr/σt < 1.7, with the smaller black holes requiring larger radial anisotropy. Compared with these three-integral models, two-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/LV and spheroidal symmetry - requiring yet more models to study black hole properties in complete generality.

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., Comment: 14 pages, 12 figures. Accepted for publication in MNRAS. Source code is available at https://github.com/stuartsale

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., 13 pages, 14 figures. MNRAS, in press

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., 15 pages, 11 figures; accepted for publication in ApJ

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., Comment: 27 pages, 18 figures, 7 tables, ApJ accepted

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., 58 pages, submitted to ApJ. Some figures excluded due to size. The entire paper is at http://www.noao.edu/noao/staff/lauer/nuker_papers.html

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., Comment: 43 pages, 12 figures, 2 tables; accepted for publication in ApJ

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), 31 pages, 34 figures; accepted for publication in MNRAS

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., Comment: 8 pages, 8 figures. Added discussion on Nbody codes updated to match published version

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)., 10 pages, 7 figures. Simplified integrator updated to match published version

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., Comment: 21 pages, 19 figures

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., 11 pages, 7 figures, 1 table

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 Mbh3X10^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., 62 pages, 18 figures, submitted to the Astrophysical Journal; revised after initial review

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., 20 pages, 25 figures, accepted for publication in MNRAS

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., Comment: 19 pages, 9 figures. Accepted for publication in MNRAS. High resolution version available at http://www-astro.physics.ox.ac.uk/~rcwh/ngc1399.ps.gz

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., 15 pages, 13 figures, accepted for publication in MNRAS

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., Accepted for publication in ApJ, 10 pages, 6 figures, LaTeX (uses emulateapj5.sty)

We present stellar kinematics for a sample of 10 early-type galaxies observed using the Space Telescope Imaging Spectrograph (STIS) aboard the Hubble Space Telescope and the Modular Spectrograph on the MDM Observatory 2.4 m telescope. These observations are a part of an ongoing program to understand the coevolution of supermassive black holes and their host galaxies. Our spectral ranges include either the calcium triplet absorption lines at 8498, 8542, and 8662 Åor the Mg b absorption at 5175 Å. The lines 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 vmax/σ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 106.5 to 10 9.3 M⊙. We replot these updated masses on the M BH-σ 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 NGC 7457, do not deviate significantly from the established relation between MBH and σ. Neither do the three galaxies that show the most evidence of a recent merger, NGC 3608, NGC 4473, and NGC 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., Accepted for publication by MNRAS

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., 35 pages, 17 postscript figures, submitted to the Astronomical Journal

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., 25 pages, accepted for publication in The Astrophysical Journal

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., Comment: 14 pages, 5 figures, accepted for publication in MNRAS

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., 21 pages, 12 figures, accepted in the Astrophysical Journal

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., 7 pages, 7 figures, uses mn.tex

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., Comment: 20 pages, 25 figures, accepted in A&A