-
Strong atomic reconstruction in twisted bilayers of highly flexible InSe: Machine-Learned Interatomic Potential and continuum model approaches
Authors:
Samuel J. Magorrian,
Anas Siddiqui,
Nicholas D. M. Hine
Abstract:
The relaxation of atomic positions to their optimal structural arrangement is crucial for understanding the emergence of new physical behavior in long scale superstructures in twisted bilayers of two-dimensional materials. The amount of deviation from a rigid moiré structure will depend on the elastic properties of the constituent monolayers which for the twisted bilayer - the more flexible the mo…
▽ More
The relaxation of atomic positions to their optimal structural arrangement is crucial for understanding the emergence of new physical behavior in long scale superstructures in twisted bilayers of two-dimensional materials. The amount of deviation from a rigid moiré structure will depend on the elastic properties of the constituent monolayers which for the twisted bilayer - the more flexible the monolayers are, the lower the energy required to deform the layers to maximize the areas with an energetically optimal interlayer arrangement of atoms. We investigate this atomic reconstruction for twisted bilayers of highly flexible InSe. Results using two methods are demonstrated - first we train a machine-learned interatomic potential (MLIP) to enable fully atomistic relaxations of small-twist-angle large-length-scale moiré supercells while retaining density functional theory (DFT) level accuracy. We find substantial out-of-plane corrugation and in-plane domain formation for a wide range of twist angles and moiré length scales. We then adapt an existing continuum approach and show that it can reproduce some, but not all, features of the fully atomistic calculations.
△ Less
Submitted 27 June, 2024;
originally announced June 2024.
-
Strain-dependent one-dimensional confinement channels in twisted bilayer 1T$'$-WTe$_2$
Authors:
Samuel J. Magorrian,
Nicholas D. M. Hine
Abstract:
The low symmetry and anistropic lattice of 1T$'$ WTe$_2$ is responsible for the existence of parallel one-dimensional channels in the moiré patterns of twisted bilayers. This gives the opportunity to explore moiré physics of a different nature to that widely observed in twisted bilayers of materials with hexagonal symmetries. Here, we combine plane-wave and linear-scaling density functional theory…
▽ More
The low symmetry and anistropic lattice of 1T$'$ WTe$_2$ is responsible for the existence of parallel one-dimensional channels in the moiré patterns of twisted bilayers. This gives the opportunity to explore moiré physics of a different nature to that widely observed in twisted bilayers of materials with hexagonal symmetries. Here, we combine plane-wave and linear-scaling density functional theory calculations to describe the electronic properties of twisted bilayer 1T$'$ WTe$_2$. For a small change in the lattice parameters of the constituent 1T$'$ WTe$_2$ monolayers, we find a substantial moiré-induced striped electrostatic potential landscape in the twisted bilayer, with a peak-to-trough magnitude $>$200~meV.
△ Less
Submitted 16 April, 2024;
originally announced April 2024.
-
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.
-
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.
-
Band alignment and interlayer hybridisation in transition metal dichalcogenide/hexagonal boron nitride heterostructures
Authors:
S. J. Magorrian,
A. J. Graham,
N. Yeung,
F. Ferreira,
P. V. Nguyen,
A. Barinov,
V. I. Fal'ko,
N. R. Wilson,
N. D. M. Hine
Abstract:
In van der Waals heterostructures, the relative alignment of bands between layers, and the resulting band hybridisation, are key factors in determining a range of electronic properties. This work examines these effects for heterostructures of transition metal dichalcogenides (TMDs) and hexagonal boron nitride (hBN), an ubiquitous combination given the role of hBN as an encapsulating material. By c…
▽ More
In van der Waals heterostructures, the relative alignment of bands between layers, and the resulting band hybridisation, are key factors in determining a range of electronic properties. This work examines these effects for heterostructures of transition metal dichalcogenides (TMDs) and hexagonal boron nitride (hBN), an ubiquitous combination given the role of hBN as an encapsulating material. By comparing results of density functional calculations with experimental angle-resolved photoemission spectroscopy (ARPES) results, we explore the hybridisation between the valence states of the TMD and hBN layers, and show that it introduces avoided crossings between the TMD and hBN bands, with umklapp processes opening `ghost' avoided crossings in individual bands. Comparison between DFT and ARPES spectra for the MoSe$_2$/hBN heterostructure shows that the valence bands of MoSe$_2$ and hBN are significantly further separated in energy in experiment as compared to DFT. We then show that a novel scissor operator can be applied to the hBN valence states in the DFT calculations, to correct the band alignment and enable quantitative comparison to ARPES, explaining avoided crossings and other features of band visibility in the ARPES spectra.
△ Less
Submitted 5 October, 2022; v1 submitted 19 July, 2022;
originally announced July 2022.
-
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.
-
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.
-
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.
-
Multifaceted moiré superlattice physics in twisted WSe$_2$ bilayers
Authors:
S. J. Magorrian,
V. V. Enaldiev,
V. Zólyomi,
Fábio Ferreira,
Vladimir I. Fal'ko,
David A. Ruiz-Tijerina
Abstract:
Lattice reconstruction in twisted transition-metal dichalcogenide (TMD) bilayers gives rise to piezo- and ferroelectric moiré potentials for electrons and holes, as well as a modulation of the hybridisation across the bilayer. Here, we develop hybrid $\mathbf{k}\cdot \mathbf{p}$ tight-binding models to describe electrons and holes in the relevant valleys of twisted TMD homobilayers with parallel (…
▽ More
Lattice reconstruction in twisted transition-metal dichalcogenide (TMD) bilayers gives rise to piezo- and ferroelectric moiré potentials for electrons and holes, as well as a modulation of the hybridisation across the bilayer. Here, we develop hybrid $\mathbf{k}\cdot \mathbf{p}$ tight-binding models to describe electrons and holes in the relevant valleys of twisted TMD homobilayers with parallel (P) and anti-parallel (AP) orientations of the monolayer unit cells. We apply these models to describe moiré superlattice effects in twisted WSe${}_2$ bilayers, in conjunction with microscopic \emph{ab initio} calculations, and considering the influence of encapsulation, pressure and an electric displacement field. Our analysis takes into account mesoscale lattice relaxation, interlayer hybridisation, piezopotentials, and a weak ferroelectric charge transfer between the layers, and describes a multitude of possibilities offered by this system, depending on the choices of P or AP orientation, twist angle magnitude, and electron/hole valley.
△ Less
Submitted 28 September, 2021; v1 submitted 10 June, 2021;
originally announced June 2021.
-
Tunable spin-orbit coupling in two-dimensional InSe
Authors:
A. Ceferino,
S. J. Magorrian,
V. Zólyomi,
D. A. Bandurin,
A. K. Geim,
A. Patanè,
Z. D. Kovalyuk,
Z. R. Kudrynskyi,
I. V. Grigorieva,
V. I. Fal'ko
Abstract:
We demonstrate that spin-orbit coupling (SOC) strength for electrons near the conduction band edge in few-layer $γ$-InSe films can be tuned over a wide range. This tunability is the result of a competition between film-thickness-dependent intrinsic and electric-field-induced SOC, potentially, allowing for electrically switchable spintronic devices. Using a hybrid $\mathbf{k\cdot p}$ tight-binding…
▽ More
We demonstrate that spin-orbit coupling (SOC) strength for electrons near the conduction band edge in few-layer $γ$-InSe films can be tuned over a wide range. This tunability is the result of a competition between film-thickness-dependent intrinsic and electric-field-induced SOC, potentially, allowing for electrically switchable spintronic devices. Using a hybrid $\mathbf{k\cdot p}$ tight-binding model, fully parameterized with the help of density functional theory computations, we quantify SOC strength for various geometries of InSe-based field-effect transistors. The theoretically computed SOC strengths are compared with the results of weak antilocalization measurements on dual-gated multilayer InSe films, interpreted in terms of Dyakonov-Perel spin relaxation due to SOC, showing a good agreement between theory and experiment.
△ Less
Submitted 8 June, 2021;
originally announced June 2021.
-
Weak ferroelectric charge transfer in layer-asymmetric bilayers of 2D semiconductors
Authors:
F. Ferreira,
V. V. Enaldiev,
V. I. Fal'ko,
S. J. Magorrian
Abstract:
In bilayers of two-dimensional (2D) semiconductors with stacking arrangements which lack inversion symmetry charge transfer between the layers due to layer-asymmetric interband hybridisation can generate a potential difference between the layers. We analyse bilayers of transition metal dichalcogenides (TMDs) - in particular, WSe$_2$ - for which we find a substantial stacking-dependent charge trans…
▽ More
In bilayers of two-dimensional (2D) semiconductors with stacking arrangements which lack inversion symmetry charge transfer between the layers due to layer-asymmetric interband hybridisation can generate a potential difference between the layers. We analyse bilayers of transition metal dichalcogenides (TMDs) - in particular, WSe$_2$ - for which we find a substantial stacking-dependent charge transfer, and InSe, for which the charge transfer is found to be negligibly small. The information obtained about TMDs is then used to map potentials generated by the interlayer charge transfer across the moiré superlattice in twistronic bilayers.
△ Less
Submitted 7 June, 2021; v1 submitted 10 March, 2021;
originally announced March 2021.
-
Structures of bulk hexagonal post-transition-metal chalcogenides from dispersion-corrected density-functional theory
Authors:
S. J. Magorrian,
V. Zolyomi,
N. D. Drummond
Abstract:
We use dispersion-corrected density-functional theory to determine the relative energies of competing polytypes of bulk layered hexagonal post-transition-metal chalcogenides, to search for the most stable structures of these potentially technologically important semiconductors. We show that there is some degree of consensus among dispersion-corrected exchange-correlation functionals regarding the…
▽ More
We use dispersion-corrected density-functional theory to determine the relative energies of competing polytypes of bulk layered hexagonal post-transition-metal chalcogenides, to search for the most stable structures of these potentially technologically important semiconductors. We show that there is some degree of consensus among dispersion-corrected exchange-correlation functionals regarding the energetic orderings of polytypes, but we find that for each material there are multiple stacking orders with relative energies of less than 1 meV per monolayer unit cell, implying that stacking faults are expected to be abundant in all post-transition-metal chalcogenides. By fitting a simple model to all our energy data, we predict that the most stable hexagonal structure has P$6_3$/mmc space group in each case, but that the stacking order differs between GaS, GaSe, GaTe, and InS on the one hand and InSe and InTe on the other. At zero pressure, the relative energies obtained with different functionals disagree by around 1-5 meV per monolayer unit cell, which is not sufficient to identify the most stable structure unambiguously; however, multi-GPa pressures reduce the number of competing phases significantly. At higher pressures, an AB$'$-stacked structure of the most stable monolayer polytype is found to be the most stable bulk structure; this structure has not been reported in experiments thus far.
△ Less
Submitted 18 March, 2021; v1 submitted 27 January, 2021;
originally announced January 2021.
-
Piezoelectric networks and ferroelectric moiré superlattice domains in twistronic WS$_2$/MoS$_2$ and WSe$_2$/MoSe$_2$ bilayers
Authors:
V. V. Enaldiev,
F. Ferreira,
S. J. Magorrian,
V. I. Fal'ko
Abstract:
Twistronic van der Waals heterostrutures offer exciting opportunities for engineering optoelectronic properties of nanomaterials. Here, we use multiscale modeling to study trapping of charge carriers and excitons by ferroelectric polarisation and piezoelectric charges by domain structures in twistronic WX$_2$/MoX$_2$ bilayers (X=S,Se). For almost aligned 2H-type bilayers, we find that holes and el…
▽ More
Twistronic van der Waals heterostrutures offer exciting opportunities for engineering optoelectronic properties of nanomaterials. Here, we use multiscale modeling to study trapping of charge carriers and excitons by ferroelectric polarisation and piezoelectric charges by domain structures in twistronic WX$_2$/MoX$_2$ bilayers (X=S,Se). For almost aligned 2H-type bilayers, we find that holes and electrons are trapped in the opposite -- WMo and XX (tungsten over molybdenum {\it versus} overlaying chalcogens) -- corners of the honeycomb domain wall network, swapping their position at a twist angle $0.2^{\circ}$, with XX corners providing $30$\,meV deep traps for the interlayer excitons for all angles. In 3R-type bilayers, both electrons and holes are trapped in triangular "3R stacking" domains, where WX$_2$ chalcogens set over MoX$_2$ molybdenums, which act as $130$\,meV deep quantum boxes for interlayer excitons for twist angles $\lesssim 1^{\circ}$, for larger angles shifting towards domain wall network XX stacking sites.
△ Less
Submitted 4 February, 2021; v1 submitted 9 November, 2020;
originally announced November 2020.
-
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.
-
Crossover from weakly indirect to direct excitons in atomically thin films of InSe
Authors:
Adrián Ceferino,
Kok Wee Song,
Samuel J. Magorrian,
Viktor Zólyomi,
Vladimir I. Fal'ko
Abstract:
We perform a $\mathbf{k \cdot p}$ theory analysis of the spectra of the lowest energy and excited states of the excitons in few-layer atomically thin films of InSe taking into account in-plane electric polarizability of the film and the influence of the encapsulation environment. For the thinner films, the lowest-energy state of the exciton is weakly indirect in momentum space, with its dispersion…
▽ More
We perform a $\mathbf{k \cdot p}$ theory analysis of the spectra of the lowest energy and excited states of the excitons in few-layer atomically thin films of InSe taking into account in-plane electric polarizability of the film and the influence of the encapsulation environment. For the thinner films, the lowest-energy state of the exciton is weakly indirect in momentum space, with its dispersion showing minima at a layer-number-dependent wave number, due to an inverted edge of a relatively flat topmost valence band branch of the InSe film spectrum and we compute the activation energy from the momentum dark exciton ground state into the bright state. For the films with more than seven In$_2$Se$_2$ layers, the exciton dispersion minimum shifts to $Γ$-point.
△ Less
Submitted 22 June, 2020; v1 submitted 14 January, 2020;
originally announced January 2020.
-
Stacking domains and dislocation networks in marginally twisted bilayers of transition metal dichalcogenides
Authors:
V. V. Enaldiev,
V. Zólyomi,
C. Yelgel,
S. J. Magorrian,
V. I. Fal'ko
Abstract:
We apply a multiscale modeling approach to study lattice reconstruction in marginally twisted bilayers of transition metal dichalcogenides (TMD). For this, we develop DFT-parametrized interpolation formulae for interlayer adhesion energies of MoSe$_2$, WSe$_2$, MoS$_2$, and WS$_2$, combine those with elasticity theory, and analyze the bilayer lattice relaxation into mesoscale domain structures. Pa…
▽ More
We apply a multiscale modeling approach to study lattice reconstruction in marginally twisted bilayers of transition metal dichalcogenides (TMD). For this, we develop DFT-parametrized interpolation formulae for interlayer adhesion energies of MoSe$_2$, WSe$_2$, MoS$_2$, and WS$_2$, combine those with elasticity theory, and analyze the bilayer lattice relaxation into mesoscale domain structures. Paying particular attention to the inversion asymmetry of TMD monolayers, we show that 3R and 2H stacking domains, separated by a network of dislocations develop for twist angles $θ^{\circ}<θ^{\circ}_P\sim 2.5^{\circ}$ and $θ^{\circ}<θ^{\circ}_{AP}\sim 1^{\circ}$ for, respectively, bilayers with parallel (P) and antiparallel (AP) orientation of the monolayer unit cells and suggest how the domain structures would manifest itself in local probe scanning of marginally twisted P- and AP-bilayers.
△ Less
Submitted 5 April, 2020; v1 submitted 28 November, 2019;
originally announced November 2019.
-
Infrared-to-violet tunable optical activity in atomic films of GaSe, InSe, and their heterostructures
Authors:
Daniel J. Terry,
Viktor Zólyomi,
Matthew Hamer,
Anastasia V. Tyurnina,
David G. Hopkinson,
Alexander M. Rakowski,
Samuel J. Magorrian,
Nick Clark,
Yuri M. Andreev,
Olga Kazakova,
Konstantin Novoselov,
Sarah J. Haigh,
Vladimir I. Fal'ko,
Roman Gorbachev
Abstract:
Two-dimensional semiconductors - atomic layers of materials with covalent intra-layer bonding and weak (van der Waals or quadrupole) coupling between the layers - are a new class of materials with great potential for optoelectronic applications. Among those, a special position is now being taken by post-transition metal chalcogenides (PTMC), InSe and GaSe. It has recently been found that the band…
▽ More
Two-dimensional semiconductors - atomic layers of materials with covalent intra-layer bonding and weak (van der Waals or quadrupole) coupling between the layers - are a new class of materials with great potential for optoelectronic applications. Among those, a special position is now being taken by post-transition metal chalcogenides (PTMC), InSe and GaSe. It has recently been found that the band gap in 2D crystals of InSe more than doubles in the monolayer compared to thick multilayer crystals, while the high mobility of conduction band electrons is promoted by their light in-plane mass. Here, we use Raman and PL measurements of encapsulated few layer samples, coupled with accurate atomic force and transmission electron microscope structural characterisation to reveal new optical properties of atomically thin GaSe preserved by hBN encapsulation. The band gaps we observe complement the spectral range provided by InSe films, so that optical activity of these two almost lattice-matched PTMC films and their heterostructures densely cover the spectrum of photons from violet to infrared. We demonstrate the realisation of the latter by the first observation of interlayer excitonic photoluminescence in few-layer InSe-GaSe heterostructures. The spatially indirect transition is direct in k-space and therefore is bright, while its energy can be tuned in a broad range by the number of layers.
△ Less
Submitted 3 October, 2018;
originally announced October 2018.
-
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.
-
Hybrid $\mathbf{k\cdot p}$-tight-binding model for intersubband optics in atomically thin InSe films
Authors:
Samuel J. Magorrian,
Adrian Ceferino,
Viktor Zólyomi,
Vladimir I. Fal'ko
Abstract:
We propose atomic films of n-doped $γ$-InSe as a platform for intersubband optics in the infrared (IR) and far infrared (FIR) range, coupled to out-of-plane polarized light. Depending on the film thickness (number of layers) of the InSe film these transitions span from $\sim 0.7$ eV for bilayer to $\sim 0.05$ eV for 15-layer InSe. We use a hybrid $\mathbf{k} \cdot \mathbf{p}$ theory and tight-bind…
▽ More
We propose atomic films of n-doped $γ$-InSe as a platform for intersubband optics in the infrared (IR) and far infrared (FIR) range, coupled to out-of-plane polarized light. Depending on the film thickness (number of layers) of the InSe film these transitions span from $\sim 0.7$ eV for bilayer to $\sim 0.05$ eV for 15-layer InSe. We use a hybrid $\mathbf{k} \cdot \mathbf{p}$ theory and tight-binding model, fully parametrized using density functional theory, to predict their oscillator strengths and thermal linewidths at room temperature.
△ Less
Submitted 12 April, 2018; v1 submitted 22 January, 2018;
originally announced January 2018.
-
Spin-orbit coupling, optical transitions, and spin pumping in mono- and few-layer InSe
Authors:
S. J. Magorrian,
V. Zólyomi,
V. I. Fal'ko
Abstract:
We show that spin-orbit coupling (SOC) in InSe enables the optical transition across the principal band gap to couple with in-plane polarized light. This transition, enabled by $p_{x,y}\leftrightarrow p_z$ hybridization due to intra-atomic SOC in both In and Se, can be viewed as a transition between two dominantly $s$- and $p_z$-orbital based bands, accompanied by an electron spin-flip. Having par…
▽ More
We show that spin-orbit coupling (SOC) in InSe enables the optical transition across the principal band gap to couple with in-plane polarized light. This transition, enabled by $p_{x,y}\leftrightarrow p_z$ hybridization due to intra-atomic SOC in both In and Se, can be viewed as a transition between two dominantly $s$- and $p_z$-orbital based bands, accompanied by an electron spin-flip. Having parametrized $\mathbf{k\cdot p}$ theory using first principles density functional theory we estimate the absorption for $σ^{\pm}$ circularly polarized photons in the monolayer as $\sim 1.5\%$, which saturates to $\sim 0.3\%$ in thicker films ($3-5$ layers). Circularly polarized light can be used to selectively excite electrons into spin-polarized states in the conduction band, which permits optical pumping of the spin polarization of In nuclei through the hyperfine interaction.
△ Less
Submitted 9 November, 2017;
originally announced November 2017.
-
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.
-
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.
-
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.
-
Electronic and optical properties of two-dimensional InSe from a DFT-parameterized tight-binding model
Authors:
S. J. Magorrian,
V. Zólyomi,
V. I. Fal'ko
Abstract:
We present a tight-binding (TB) model and $\mathbf{k\cdot p}$ theory for electrons in monolayer and few-layer InSe. The model is constructed from a basis of all $s$ and $p$ valence orbitals on both indium and selenium atoms, with tight-binding parameters obtained from fitting to independently computed density functional theory (DFT) band structures for mono- and bilayer InSe. For the valence and c…
▽ More
We present a tight-binding (TB) model and $\mathbf{k\cdot p}$ theory for electrons in monolayer and few-layer InSe. The model is constructed from a basis of all $s$ and $p$ valence orbitals on both indium and selenium atoms, with tight-binding parameters obtained from fitting to independently computed density functional theory (DFT) band structures for mono- and bilayer InSe. For the valence and conduction band edges of few-layer InSe, which appear to be in the vicinity of the $Γ$ point, we calculate the absorption coefficient for the principal optical transitions as a function of the number of layers, $N$. We find a strong dependence on $N$ of the principal optical transition energies, selection rules, and optical oscillation strengths, in agreement with recent observations \cite{Bandurin2016}. Also, we find that the conduction band electrons are relatively light ($m \propto 0.14-0.18 m_e$), in contrast to an almost flat, and slightly inverted, dispersion of valence band holes near the $Γ$-point, which is found for up to $N \propto 6$.
△ Less
Submitted 2 December, 2016; v1 submitted 1 November, 2016;
originally announced November 2016.
-
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.
-
Polarization memory in the nonpolar magnetic ground state of multiferroic CuFeO2
Authors:
J. Beilsten-Edmands,
S. J. Magorrian,
F. R. Foronda,
D. Prabhakaran,
P. G. Radaelli,
R. D. Johnson
Abstract:
We investigate polarization memory effects in single-crystal CuFeO2, which has a magnetically-induced ferroelectric phase at low temperatures and applied B fields between 7.5 and 13 T. Following electrical poling of the ferroelectric phase, we find that the nonpolar collinear antiferromagnetic ground state at B = 0 T retains a strong memory of the polarization magnitude and direction, such that up…
▽ More
We investigate polarization memory effects in single-crystal CuFeO2, which has a magnetically-induced ferroelectric phase at low temperatures and applied B fields between 7.5 and 13 T. Following electrical poling of the ferroelectric phase, we find that the nonpolar collinear antiferromagnetic ground state at B = 0 T retains a strong memory of the polarization magnitude and direction, such that upon re-entering the ferroelectric phase a net polarization of comparable magnitude to the initial polarization is recovered in the absence of external bias. This memory effect is very robust: in pulsed-magnetic-field measurements, several pulses into the ferroelectric phase with reverse bias are required to switch the polarization direction, with significant switching only seen after the system is driven out of the ferroelectric phase and ground state either magnetically (by application of B > 13 T) or thermally. The memory effect is also largely insensitive to the magnetoelastic domain composition, since no change in the memory effect is observed for a sample driven into a single-domain state by application of stress in the [1-10] direction. On the basis of Monte Carlo simulations of the ground state spin configurations, we propose that the memory effect is due to the existence of helical domain walls within the nonpolar collinear antiferromagnetic ground state, which would retain the helicity of the polar phase for certain magnetothermal histories.
△ Less
Submitted 7 October, 2016;
originally announced October 2016.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.