Showing 1–50 of 56 results for author: Magorrian, J

Mapping dust in the giant molecular cloud Orion A

Authors: Amery Gration, John Magorrian

Abstract: The Sun is located close to the Galactic mid-plane, meaning that we observe the Galaxy through significant quantities of dust. Moreover, the vast majority of the Galaxy's stars also lie in the disc, meaning that dust has an enormous impact on the massive astrometric, photometric and spectroscopic surveys of the Galaxy that are currently underway. To exploit the data from these surveys we require g… ▽ More The Sun is located close to the Galactic mid-plane, meaning that we observe the Galaxy through significant quantities of dust. Moreover, the vast majority of the Galaxy's stars also lie in the disc, meaning that dust has an enormous impact on the massive astrometric, photometric and spectroscopic surveys of the Galaxy that are currently underway. To exploit the data from these surveys we require good three-dimensional maps of the Galaxy's dust. We present a new method for making such maps in which we form the best linear unbiased predictor of the extinction at an arbitrary point based on the extinctions for a set of observed stars. This method allows us to avoid the artificial inhomogeneities (so-called 'fingers of God') and resolution limits that are characteristic of many published dust maps. Moreover, it requires minimal assumptions about the statistical properties of the interstellar medium. In fact, we require only a model of the first and second moments of the dust density field. The method is suitable for use with directly measured extinctions, such as those provided by the Rayleigh-Jeans colour excess method, and inferred extinctions, such as those provided by hierarchical Bayesian models like StarHorse. We test our method by mapping dust in the region of the giant molecular cloud Orion A. Our results indicate a foreground dust cloud at a distance of 350 pc, which has been identified in work by another author. △ Less

Submitted 23 January, 2024; originally announced January 2024.

Comments: 21 pages, 30 figures

The JWST Galactic Center Survey -- A White Paper

Authors: Rainer Schoedel, Steve Longmore, Jonny Henshaw, Adam Ginsburg, John Bally, Anja Feldmeier, Matt Hosek, Francisco Nogueras Lara, Anna Ciurlo, Mélanie Chevance, J. M. Diederik Kruijssen, Ralf Klessen, Gabriele Ponti, Pau Amaro-Seoane, Konstantina Anastasopoulou, Jay Anderson, Maria Arias, Ashley T. Barnes, Cara Battersby, Giuseppe Bono, Lucía Bravo Ferres, Aaron Bryant, Miguel Cano Gonzáalez, Santi Cassisi, Leonardo Chaves-Velasquez , et al. (85 additional authors not shown)

Abstract: The inner hundred parsecs of the Milky Way hosts the nearest supermassive black hole, largest reservoir of dense gas, greatest stellar density, hundreds of massive main and post main sequence stars, and the highest volume density of supernovae in the Galaxy. As the nearest environment in which it is possible to simultaneously observe many of the extreme processes shaping the Universe, it is one of… ▽ More The inner hundred parsecs of the Milky Way hosts the nearest supermassive black hole, largest reservoir of dense gas, greatest stellar density, hundreds of massive main and post main sequence stars, and the highest volume density of supernovae in the Galaxy. As the nearest environment in which it is possible to simultaneously observe many of the extreme processes shaping the Universe, it is one of the most well-studied regions in astrophysics. Due to its proximity, we can study the center of our Galaxy on scales down to a few hundred AU, a hundred times better than in similar Local Group galaxies and thousands of times better than in the nearest active galaxies. The Galactic Center (GC) is therefore of outstanding astrophysical interest. However, in spite of intense observational work over the past decades, there are still fundamental things unknown about the GC. JWST has the unique capability to provide us with the necessary, game-changing data. In this White Paper, we advocate for a JWST NIRCam survey that aims at solving central questions, that we have identified as a community: i) the 3D structure and kinematics of gas and stars; ii) ancient star formation and its relation with the overall history of the Milky Way, as well as recent star formation and its implications for the overall energetics of our galaxy's nucleus; and iii) the (non-)universality of star formation and the stellar initial mass function. We advocate for a large-area, multi-epoch, multi-wavelength NIRCam survey of the inner 100\,pc of the Galaxy in the form of a Treasury GO JWST Large Program that is open to the community. We describe how this survey will derive the physical and kinematic properties of ~10,000,000 stars, how this will solve the key unknowns and provide a valuable resource for the community with long-lasting legacy value. △ Less

Submitted 14 March, 2024; v1 submitted 18 October, 2023; originally announced October 2023.

Comments: This White Paper will be updated when required (e.g. new authors joining, editing of content). Most recent update: 24 Oct 2023

On the self-consistent time-dependent linearized response of stellar discs to external perturbations

Authors: Dominic Dootson, John Magorrian

Abstract: We study the explicitly time-dependent response of a razor-thin axisymmetric disc to externally imposed perturbations by recasting the linearized Collisionless Boltzmann equation as an integral equation and applying Kalnajs' matrix method. As an application we consider the idealized problem of calculating the dynamical friction torque on a steadily rotating, two-dimensional bar. We consider two ch… ▽ More We study the explicitly time-dependent response of a razor-thin axisymmetric disc to externally imposed perturbations by recasting the linearized Collisionless Boltzmann equation as an integral equation and applying Kalnajs' matrix method. As an application we consider the idealized problem of calculating the dynamical friction torque on a steadily rotating, two-dimensional bar. We consider two choices of basis functions in the matrix method, showing that both lead to comparable results. The torques from our linearised calculation are in excellent agreement with those measured from $N$-body simulation, as long as the bar perturbation does not resonate with a significant fraction of the disc's stars. △ Less

Submitted 31 May, 2022; originally announced May 2022.

Comments: 15 pages, 8 figures

Self-consistent modelling of the Milky Way's Nuclear Stellar Disc

Authors: Mattia C. Sormani, Jason L. Sanders, Tobias K. Fritz, Leigh C. Smith, Ortwin Gerhard, Rainer Schoedel, John Magorrian, Nadine Neumayer, Francisco Nogueras-Lara, Anja Feldmeier-Krause, Alessandra Mastrobuono-Battisti, Mathias Schultheis, Banafsheh Shahzamanian, Eugene Vasiliev, Ralf S. Klessen, Philip Lucas, Dante Minniti

Abstract: The Nuclear Stellar Disc (NSD) is a flattened high-density stellar structure that dominates the gravitational field of the Milky Way at Galactocentric radius $30\lesssim R\lesssim 300$ pc. 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 normalised kinematic… ▽ More The Nuclear Stellar Disc (NSD) is a flattened high-density stellar structure that dominates the gravitational field of the Milky Way at Galactocentric radius $30\lesssim R\lesssim 300$ pc. 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 normalised 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 marginalised probability distributions give a total mass of $M_{\rm NSD} = 10.5^{+1.1}_{-1.0} \times10^8 \,{\rm M_\odot}$, roughly exponential radial and vertical scale-lengths of $R_{\rm disc} = 88.6^{+9.2}_{-6.9}$ pc and $H_{\rm disc}=28.4^{+5.5}_{-5.5}$ pc respectively, and a velocity dispersion $σ\simeq 70$ km/s 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. △ Less

Submitted 4 March, 2022; v1 submitted 24 November, 2021; originally announced November 2021.

Comments: Accepted for publication in MNRAS

Stellar dynamics in the periodic cube

Authors: John Magorrian

Abstract: We use the problem of dynamical friction within the periodic cube to illustrate the application of perturbation theory in stellar dynamics, testing its predictions against measurements from $N$-body simulation. Our development is based on the explicitly time-dependent Volterra integral equation for the cube's linear response, which avoids the subtleties encountered in analyses based on complex fre… ▽ More We use the problem of dynamical friction within the periodic cube to illustrate the application of perturbation theory in stellar dynamics, testing its predictions against measurements from $N$-body simulation. Our development is based on the explicitly time-dependent Volterra integral equation for the cube's linear response, which avoids the subtleties encountered in analyses based on complex frequency. We obtain an expression for the self-consistent response of the cube to steady stirring by an external perturber. From this we show how to obtain the familiar Chandrasekhar dynamical friction formula and construct an elementary derivation of the Lenard--Balescu equation for the secular quasilinear evolution of an isolated cube composed of $N$ equal-mass stars. We present an alternative expression for the (real-frequency) van Kampen modes of the cube and show explicitly how to decompose any linear perturbation of the cube into a superposition of such modes. △ Less

Submitted 16 September, 2021; v1 submitted 18 June, 2021; originally announced June 2021.

Comments: 13 pages. References added to match accepted version

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 507, Issue 4, November 2021, Pages 4840-4851

Jeans modelling of the Milky Way's nuclear stellar disc

Authors: Mattia C. Sormani, John Magorrian, Francisco Nogueras-Lara, Nadine Neumayer, Ralph Schoenrich, Ralf S. Klessen, Alessandra Mastrobuono-Battisti

Abstract: 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) th… ▽ More 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(r<100{\, \rm pc}) = 3.9 \pm 1 \times 10^8 {\, \rm M_\odot}$ and a total mass of $M_{\rm NSD} = 6.9 \pm 2 \times 10^8 {\, \rm M_\odot}$. (ii) The NSD might be the first example of a vertically biased disc, i.e. with ratio between the vertical and radial velocity dispersion $σ_z/σ_R>1$. 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 $σ_z/σ_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. △ Less

Submitted 9 September, 2020; v1 submitted 13 July, 2020; originally announced July 2020.

Comments: Accepted for publication in MNRAS

Orbit-superposition models of discrete, incomplete stellar kinematics: application to the Galactic centre

Authors: John Magorrian

Abstract: 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 kine… ▽ More 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 $ρ\propto r^{-γ}$ with $γ=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 $γ$ 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. △ Less

Submitted 21 January, 2019; v1 submitted 16 September, 2018; originally announced September 2018.

Comments: 17 pages, MNRAS accepted

Large-scale three-dimensional Gaussian process extinction mapping

Authors: S. E. Sale, J. Magorrian

Abstract: 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… ▽ More 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. △ Less

Submitted 31 August, 2017; originally announced August 2017.

Comments: 15 pages, 6 figures, submitted for publication in MNRAS

A theoretical explanation for the Central Molecular Zone asymmetry

Authors: Mattia C. Sormani, Robin G. Tress, Matthew Ridley, Simon C. O. Glover, Ralf S. Klessen, James Binney, John Magorrian, Rowan Smith

Abstract: 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 rema… ▽ More 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. △ Less

Submitted 15 December, 2017; v1 submitted 12 July, 2017; originally announced July 2017.

Comments: Accepted for publication in MNRAS. Videos of the simulations are available at http://www.ita.uni-heidelberg.de/~mattia/download.html

Nuclear Spirals in the inner Milky Way

Authors: Matthew Ridley, Mattia C. Sormani, Robin G. Treß, John Magorrian, Ralf S. Klessen

Abstract: 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… ▽ More 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. △ Less

Submitted 21 April, 2017; v1 submitted 12 April, 2017; originally announced April 2017.

Comments: Accepted for publication in MNRAS, 12 pages, 12 figures

Journal ref: Mon Not R Astron Soc (2017) 469 (2): 2251-2262

A low upper mass limit for the central black hole in the late-type galaxy NGC 4414

Authors: Sabine Thater, Davor Krajnović, Martin A. Bourne, Michele Cappellari, Tim de Zeeuw, Eric Emsellem, John Magorrian, Richard M. McDermid, Marc Sarzi, Glenn van de Ven

Abstract: 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 arcse… ▽ More 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. △ Less

Submitted 24 October, 2016; v1 submitted 13 October, 2016; originally announced October 2016.

Comments: 20 pages, 25 figures, accepted in A&A

Journal ref: A&A 597, A18 (2017)

The secular evolution of discrete quasi-Keplerian systems. I. Kinetic theory of stellar clusters near black holes

Authors: Jean-Baptiste Fouvry, Christophe Pichon, John Magorrian

Abstract: 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… ▽ More 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. △ Less

Submitted 22 February, 2017; v1 submitted 17 June, 2016; originally announced June 2016.

Comments: 24 pages, 3 figures

Journal ref: A&A 598, A71 (2017)

Secular diffusion in discrete self-gravitating tepid discs II: accounting for swing amplification via the matrix method

Authors: Jean-Baptiste Fouvry, Christophe Pichon, John Magorrian, Pierre-Henri Chavanis

Abstract: 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 fluctuati… ▽ More 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. △ Less

Submitted 24 July, 2015; originally announced July 2015.

Comments: 21 pages, 19 figures

Journal ref: A&A 584, A129 (2015)

Gas flow in barred potentials - III. Effects of varying the Quadrupole

Authors: Mattia C. Sormani, James Binney, John Magorrian

Abstract: 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 W… ▽ More 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 $Ω=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. △ Less

Submitted 2 November, 2017; v1 submitted 11 July, 2015; originally announced July 2015.

Comments: accepted for publication in MNRAS. This version corrects a typo in the normalization factor that is present in the published version

Gas flow in barred potentials II. Bar Driven Spiral Arms

Authors: Mattia C. Sormani, James Binney, John Magorrian

Abstract: 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,… ▽ More 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. △ Less

Submitted 16 May, 2015; v1 submitted 19 April, 2015; originally announced April 2015.

Comments: Accepted for publication in MNRAS

Gas flow in barred potentials

Authors: Mattia C. Sormani, James Binney, John Magorrian

Abstract: 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 w… ▽ More 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. △ Less

Submitted 19 April, 2015; v1 submitted 9 February, 2015; originally announced February 2015.

Comments: Accepted for publication in MNRAS

Marginal likelihoods of distances and extinctions to stars: computation and compact representation

Authors: S. E. Sale, J. Magorrian

Abstract: 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 appl… ▽ More 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. △ Less

Submitted 12 January, 2015; v1 submitted 18 November, 2014; originally announced November 2014.

Comments: 14 pages, 12 figures. Accepted for publication in MNRAS. Source code is available at https://github.com/stuartsale

Recognizing the fingerprints of the Galactic bar: a quantitative approach to comparing model (l,v) distributions to observation

Authors: Mattia C. Sormani, John Magorrian

Abstract: 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-… ▽ More 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 χ^2 or envelope distances at identifying the correct underlying galaxy model. △ Less

Submitted 19 April, 2015; v1 submitted 29 October, 2014; originally announced October 2014.

Comments: Accepted for publication in MNRAS

Three-dimensional extinction mapping using Gaussian random fields

Authors: S. E. Sale, J. Magorrian

Abstract: 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… ▽ More 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. △ Less

Submitted 19 August, 2014; v1 submitted 4 June, 2014; originally announced June 2014.

Comments: 16 pages, 12 figures. Submitted for publication in MNRAS. Text revised

Bayes versus the virial theorem: inferring the potential of a galaxy from a kinematical snapshot

Authors: John Magorrian

Abstract: 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… ▽ More 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. △ Less

Submitted 21 October, 2013; v1 submitted 25 March, 2013; originally announced March 2013.

Comments: 19 pages, 7 figures. Updated to match accepted version

Three-dimensional Keplerian orbit-superposition models of the nucleus of M31

Authors: Calum K. Brown, John Magorrian

Abstract: 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… ▽ More 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. △ Less

Submitted 14 March, 2013; originally announced March 2013.

Comments: 13 pages, 14 figures. MNRAS, in press

Dynamical masses of early-type galaxies: a comparison to lensing results and implications for the stellar IMF and the distribution of dark matter

Authors: J. Thomas, R. P. Saglia, R. Bender, D. Thomas, K. Gebhardt, J. Magorrian, E. M. Corsini, G. Wegner, S. Seitz

Abstract: 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 ste… ▽ More 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. △ Less

Submitted 20 December, 2011; v1 submitted 17 March, 2011; originally announced March 2011.

Comments: 21 pages, 22 figures (MNRAS, 2011, 415, 545), corrected version of fig. 5 in the appendix (cf. MNRAS, 2011, 418, 2815)

The M-sigma and M-L Relations in Galactic Bulges and Determinations of their Intrinsic Scatter

Authors: Kayhan Gultekin, Douglas O. Richstone, Karl Gebhardt, Tod R. Lauer, Scott Tremaine, M. C. Aller, Ralf Bender, Alan Dressler, S. M. Faber, Alexei V. Filippenko, Richard Green, Luis C. Ho, John Kormendy, John Magorrian, Jason Pinkney, Christos Siopis

Abstract: 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… ▽ More 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. △ Less

Submitted 27 March, 2009; originally announced March 2009.

Comments: 27 pages, 18 figures, 7 tables, ApJ accepted

Journal ref: Astrophys.J.698:198-221,2009

The flattening and the orbital structure of early-type galaxies and collisionless N-body binary disk mergers

Authors: J. Thomas, R. Jesseit, R. P. Saglia, R. Bender, A. Burkert, E. M. Corsini, K. Gebhardt, J. Magorrian, T. Naab, D. Thomas, G. Wegner

Abstract: 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 s… ▽ More 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. △ Less

Submitted 13 November, 2008; originally announced November 2008.

Comments: 14 pages, 8 figures; accepted for publication in MNRAS

Dark matter scaling relations and the assembly epoch of Coma early-type galaxies

Authors: J. Thomas, R. P. Saglia, R. Bender, D. Thomas, K. Gebhardt, J. Magorrian, E. M. Corsini, G. Wegner

Abstract: 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 asymptoti… ▽ More 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. △ Less

Submitted 23 September, 2008; originally announced September 2008.

Comments: 15 pages, 11 figures; accepted for publication in ApJ

Journal ref: Astrophys.J.691:770-782,2009

A Stellar Dynamical Measurement of the Black Hole Mass in the Maser Galaxy NGC 4258

Authors: Christos Siopis, Karl Gebhardt, Tod R. Lauer, John Kormendy, Jason Pinkney, Douglas Richstone, S. M. Faber, Scott Tremaine, M. C. Aller, Ralf Bender, Gary Bower, Alan Dressler, Alexei V. Filippenko, Richard Green, Luis C. Ho, John Magorrian

Abstract: 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… ▽ More 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. △ Less

Submitted 28 August, 2008; originally announced August 2008.

Comments: 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

Journal ref: Astrophys.J.693:946-969,2009

Structure and dynamics of galaxies with a low surface brightness disc. I. The stellar and ionised-gas kinematics

Authors: A. Pizzella, E. M. Corsini, M. Sarzi, J. Magorrian, J. Mendez-Abreu, L. Coccato, L. Morelli, F. Bertola

Abstract: (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 obje… ▽ More (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. △ Less

Submitted 18 March, 2008; originally announced March 2008.

Comments: Accepted for pubblication on MNRAS (11 figures and 15 tables)

Tidal Disruption of Stellar Objects by Hard Supermassive Black Hole Binaries

Authors: Xian Chen, F. K. Liu, John Magorrian

Abstract: Supermassive black hole binaries (SMBHBs) are expected by the hierarchical galaxy formation model in $Λ$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.… ▽ More Supermassive black hole binaries (SMBHBs) are expected by the hierarchical galaxy formation model in $Λ$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. △ Less

Submitted 3 December, 2007; originally announced December 2007.

Comments: 43 pages, 12 figures, 2 tables; accepted for publication in ApJ

Multi-mass schemes for collisionless N-body simulations

Authors: Mimi Zhang, John Magorrian

Abstract: 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 ``interestin… ▽ More 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. △ Less

Submitted 6 August, 2008; v1 submitted 17 November, 2007; originally announced November 2007.

Comments: 8 pages, 8 figures. Added discussion on Nbody codes updated to match published version

Journal ref: Mon. Not. R. Astron. Soc. 387, 1719-1726 (2008)

Dynamical modelling of luminous and dark matter in 17 Coma early-type galaxies

Authors: J. Thomas, R. P. Saglia, R. Bender, D. Thomas, K. Gebhardt, J. Magorrian, E. M. Corsini, G. Wegner

Abstract: 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… ▽ More 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) △ Less

Submitted 5 September, 2007; originally announced September 2007.

Comments: 31 pages, 34 figures; accepted for publication in MNRAS

Grommet: an N-body code for high-resolution simulations of individual galaxies

Authors: John Magorrian

Abstract: 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 be… ▽ More 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). △ Less

Submitted 7 November, 2007; v1 submitted 15 March, 2007; originally announced March 2007.

Comments: 10 pages, 7 figures. Simplified integrator updated to match published version

Journal ref: Mon.Not.Roy.Astron.Soc. 381 (2007) 1663-1671

Constraining black hole masses from stellar kinematics by summing over all possible distribution functions

Authors: John Magorrian

Abstract: 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. A… ▽ More 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. △ Less

Submitted 8 September, 2006; originally announced September 2006.

Comments: 11 pages, 7 figures, 1 table

Journal ref: Mon.Not.Roy.Astron.Soc.373:425-434,2006

The Masses of Nuclear Black Holes in Luminous Elliptical Galaxies and Implications for the Space Density of the Most Massive Black Holes

Authors: Tod R. Lauer, S. M. Faber, Douglas Richstone, Karl Gebhardt, Scott Tremaine, Marc Postman, Alan Dressler, M. C. Aller, Alexei V. Filippenko, Richard Green, Luis C. Ho, John Kormendy, John Magorrian, Jason Pinkney

Abstract: 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 1… ▽ More 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. △ Less

Submitted 5 January, 2007; v1 submitted 29 June, 2006; originally announced June 2006.

Comments: 62 pages, 18 figures, submitted to the Astrophysical Journal; revised after initial review

Journal ref: Astrophys.J.662:808-834,2007

The central kinematics of NGC 1399 measured with 14pc resolution

Authors: R. C. W. Houghton, J. Magorrian, M. Sarzi, N. Thatte, R. L. Davies, D. Krajnovic

Abstract: We present near infra-red (NIR) adaptive optics assisted spectroscopic observations of the CO ($Δμ=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 pro… ▽ More We present near infra-red (NIR) adaptive optics assisted spectroscopic observations of the CO ($Δμ=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. △ Less

Submitted 10 October, 2005; originally announced October 2005.

Comments: 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

Journal ref: Mon.Not.Roy.Astron.Soc.367:2-18,2006

Regularized orbit models unveiling the stellar structure and dark matter halo of the Coma elliptical NGC 4807

Authors: J. Thomas, R. P. Saglia, R. Bender, D. Thomas, K. Gebhardt, J. Magorrian, E. M. Corsini, G. Wegner

Abstract: 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 simulatio… ▽ More 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. △ Less

Submitted 21 April, 2005; originally announced April 2005.

Comments: 20 pages, 25 figures, accepted for publication in MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc. 360 (2005) 1355-1372

The Centers of Early-Type Galaxies with HST. V. New WFPC2 Photometry

Authors: Tod R. Lauer, S. M. Faber, Karl Gebhardt, Douglas Richstone, Scott Tremaine, Edward A. Ajhar, M. C. Aller, Ralf Bender, Alan Dressler, Alexei V. Filippenko, Richard Green, Carl J. Grillmair, Luis C. Ho, John Kormendy, John Magorrian, Jason Pinkney, Christos Siopis

Abstract: 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 al… ▽ More 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. △ Less

Submitted 1 December, 2004; originally announced December 2004.

Comments: 76 pages, 40 figures, submitted to the Astronomical Journal. The complete postscript version of the paper is available at http://www.noao.edu/noao/staff/lauer/nuker.html (the Nuker Team page)

Mapping stationary axisymmetric phase-space distribution functions by orbit libraries

Authors: J. Thomas, R. P. Saglia, R. Bender, D. Thomas, K. Gebhardt, J. Magorrian, D. Richstone

Abstract: 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… ▽ More 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. △ Less

Submitted 1 June, 2004; originally announced June 2004.

Comments: 15 pages, 13 figures, accepted for publication in MNRAS

Black Hole Mass Determinations From Orbit Superposition Models are Reliable

Authors: Douglas Richstone, Karl Gebhardt, Monique Aller, Ralf Bender, Gary Bower, Alan Dressler, S. M. Faber, Alexei V. Filippenko, Richard Green, Luis C. Ho, John Kormendy, Tod R. Lauer, John Magorrian, Jason Pinkney, Christos Siopis, Scott Tremaine

Abstract: 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 gi… ▽ More 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. △ Less

Submitted 10 March, 2004; originally announced March 2004.

Comments: submitted to the Astrophysical Journal Letters

Central DM density cuspiness in LSB's: a stellar kinematics approach

Authors: A. Pizzella, E. M. Corsini, F. Bertola, L. Coccato, J. Magorrian, M. Sarzi, J. G. Funes S. J.

Abstract: 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… ▽ More 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. △ Less

Submitted 3 November, 2003; originally announced November 2003.

Comments: 2 pages, 2 figures, to appear in Proc. IAU Symp. 220, "Dark Matter in Galaxies" eds. S. Ryder, D.J. Pisano, M. Walker, & K. Freeman (San Francisco: ASP)

Generating Equilibrium Dark Matter Halos: Inadequacies of the Local Maxwellian Approximation

Authors: Stelios Kazantzidis, John Magorrian, Ben Moore

Abstract: 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 con… ▽ More 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. △ Less

Submitted 28 January, 2004; v1 submitted 18 September, 2003; originally announced September 2003.

Comments: Accepted for publication in ApJ, 10 pages, 6 figures, LaTeX (uses emulateapj5.sty)

Journal ref: Astrophys.J.601:37-46,2004

Kinematics of Ten Early-Type Galaxies from HST and Ground-Based Spectroscopy

Authors: Jason Pinkney, Karl Gebhardt, Ralf Bender, Gary Bower, Alan Dressler, S. M. Faber, Alexei V. Filippenko, Richard Green, Luis C. Ho, John Kormendy, Tod Lauer, John Magorrian, Douglas Richstone, Scott Tremaine

Abstract: 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 fi… ▽ More 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. △ Less

Submitted 23 June, 2003; originally announced June 2003.

Comments: 43 pages, accepted for publication in the Astrophysical Journal, high resolution version found at http://hoku.as.utexas.edu/~gebhardt/pinkney.ps

Journal ref: Astrophys.J.596:903-929,2003

Active galactic nuclei and the minor merger hypothesis

Authors: Philip Kendall, John Magorrian, J. E. Pringle

Abstract: 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… ▽ More 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. △ Less

Submitted 6 May, 2003; originally announced May 2003.

Comments: Accepted for publication by MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc. 346 (2003) 1078

Axisymmetric Dynamical Models of the Central Regions of Galaxies

Authors: Karl Gebhardt, Douglas Richstone, Scott Tremaine, Tod R. Lauer, Ralf Bender, Gary Bower, Alan Dressler, S. M. Faber, Alexei V. Filippenko, Richard Green, Carl Grillmair, Luis C. Ho, John Kormendy, John Magorrian, Jason Pinkney

Abstract: 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 d… ▽ More 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. △ Less

Submitted 8 October, 2002; v1 submitted 23 September, 2002; originally announced September 2002.

Comments: 25 pages, accepted for publication in The Astrophysical Journal

Journal ref: Astrophys.J. 583 (2003) 92-115

Galaxies with a Central Minimum in Stellar Luminosity Density

Authors: Tod R. Lauer, Karl Gebhardt, Douglas Richstone, Scott Tremaine, Ralf Bender, Gary Bower, Alan Dressler, S. M. Faber, Alexei V. Filippenko, Richard Green, Carl J. Grillmair, Luis C. Ho, John Kormendy, John Magorrian, Jason Pinkney, S. Laine, Marc Postman, Roeland P. van der Marel

Abstract: 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… ▽ More 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. △ Less

Submitted 7 June, 2002; originally announced June 2002.

Comments: 35 pages, 17 postscript figures, submitted to the Astronomical Journal

The slope of the black-hole mass versus velocity dispersion correlation

Authors: Scott Tremaine, Karl Gebhardt, Ralf Bender, Gary Bower, Alan Dressler, S. M. Faber, Alexei V. Filippenko, Richard Green, Carl Grillmair, Luis C. Ho, John Kormendy, Tod R. Lauer, John Magorrian, Jason Pinkney, Douglas Richstone

Abstract: 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 i… ▽ More 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. △ Less

Submitted 26 March, 2002; originally announced March 2002.

Comments: 37 pages, 9 figures

Journal ref: Astrophys.J.574:740-753,2002

Mass profiles and anisotropies of early-type galaxies

Authors: John Magorrian, David Ballantyne

Abstract: 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… ▽ More 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. △ Less

Submitted 4 December, 2000; originally announced December 2000.

Comments: 14 pages, 5 figures, accepted for publication in MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc. 322 (2001) 702

Evidence of a Supermassive Black Hole in the Galaxy NGC 1023 from the Nuclear Stellar Dynamics

Authors: G. A. Bower, R. F. Green, R. Bender, K. Gebhardt, T. R. Lauer, J. Magorrian, D. O. Richstone, A. Danks, T. Gull, J. Hutchings, C. Joseph, M. E. Kaiser, D. Weistrop, B. Woodgate, C. Nelson, E. M. Malumuth

Abstract: 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… ▽ More 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. △ Less

Submitted 9 November, 2000; originally announced November 2000.

Comments: 21 pages, 12 figures, accepted in the Astrophysical Journal

Black Hole Mass Estimates from Reverberation Mapping and from Spatially Resolved Kinematics

Authors: Karl Gebhardt, John Kormendy, Luis Ho, Ralf Bender, Gary Bower, Alan Dressler, S. M. Faber, Alexei Filippenko, Richard Green, Carl Grillmair, Tod Lauer, John Magorrian, Jason Pinkney, Douglas Richstone, Scott Tremaine

Abstract: 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 consis… ▽ More 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. △ Less

Submitted 5 September, 2000; v1 submitted 10 July, 2000; originally announced July 2000.

Comments: 4 pages, 1 figure, ApJ Letters accepted, minor revisions

Journal ref: Astrophys.J. 543 (2000) L5

A Relationship Between Nuclear Black Hole Mass and Galaxy Velocity Dispersion

Authors: Karl Gebhardt, Ralf Bender, Gary Bower, Alan Dressler, S. M. Faber, Alexei V. Filippenko, Richard Green, Carl Grillmair, Luis C. Ho, John Kormendy, Tod R. Lauer, John Magorrian, Jason Pinkney, Douglas Richstone, Scott Tremaine

Abstract: 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… ▽ More 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. △ Less

Submitted 29 June, 2000; v1 submitted 20 June, 2000; originally announced June 2000.

Comments: 4 pages, accepted ApJ Letters, figure added and minor corrections

Journal ref: Astrophys.J.539:L13,2000

Axisymmetric, 3-Integral Models of Galaxies: A Massive Black Hole in NGC3379

Authors: Karl Gebhardt, Douglas Richstone, John Kormendy, Tod Lauer, Edward Ajhar, Ralf Bender, Alan Dressler, Sandra Faber, Carl Grillmair, John Magorrian, Scott Tremaine

Abstract: 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… ▽ More 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. △ Less

Submitted 1 December, 1999; originally announced December 1999.

Comments: 15 pages, emulateapj sty, accepted for publication in The Astronomical Journal