Source author record

John Magorrian

John Magorrian appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

13works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

13 published item(s)

preprint2022arXiv

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

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.

preprint2016arXiv

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

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.

preprint2015arXiv

Gas flow in barred potentials

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.

preprint2015arXiv

Gas flow in barred potentials II. Bar Driven Spiral Arms

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.

preprint2015arXiv

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

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.

preprint2015arXiv

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

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.

preprint2013arXiv

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

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.

preprint2013arXiv

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

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.

preprint2009arXiv

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

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.

preprint2008arXiv

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

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.

preprint2007arXiv

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

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.

preprint2000arXiv

A Relationship Between Nuclear Black Hole Mass and Galaxy Velocity Dispersion

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.

preprint1999arXiv

Rates of tidal disruption of stars by massive central black holes

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