-
The two phases of core formation -- orbital evolution in the centres of ellipticals with supermassive black hole binaries
Authors:
Matteo Frigo,
Thorsten Naab,
Antti Rantala,
Peter H. Johansson,
Bianca Neureiter,
Jens Thomas,
Francesco Rizzuto
Abstract:
The flat stellar density cores of massive elliptical galaxies form rapidly due to sinking supermassive black holes (SMBH) in gas-poor galaxy mergers. After the SMBHs form a bound binary, gravitational slingshot interactions with nearby stars drive the core regions towards a tangentially biased stellar velocity distribution. We use collisionless galaxy merger simulations with accurate collisional o…
▽ More
The flat stellar density cores of massive elliptical galaxies form rapidly due to sinking supermassive black holes (SMBH) in gas-poor galaxy mergers. After the SMBHs form a bound binary, gravitational slingshot interactions with nearby stars drive the core regions towards a tangentially biased stellar velocity distribution. We use collisionless galaxy merger simulations with accurate collisional orbit integration around the central SMBHs to demonstrate that the removal of stars from the centre by slingshot kicks accounts for the entire change in velocity anisotropy. The rate of strong (unbinding) kicks is constant over several hundred Myr at $\sim 3 \ M_\odot \rm yr^{-1}$ for our most massive SMBH binary ($M_{\rm BH} = 1.7 \times 10^{10} M_\odot$). Using a frequency-based orbit classification scheme (box, x-tube, z-tube, rosette) we demonstrate that slingshot kicks mostly affect box orbits with small pericentre distances, leading to a velocity anisotropy of $β\lesssim -0.6$ within several hundred Myr as observed in massive ellipticals with large cores. We show how different SMBH masses affect the orbital structure of the merger remnants and present a kinematic tomography connecting orbit families to integral field kinematic features. Our direct orbit classification agrees remarkably well with a modern triaxial Schwarzschild analysis applied to simulated mock kinematic maps.
△ Less
Submitted 21 September, 2021;
originally announced September 2021.
-
SMART: A new implementation of Schwarzschild's Orbit Superposition technique for triaxial galaxies and its application to an N-body merger simulation
Authors:
Bianca Neureiter,
Jens Thomas,
Roberto Saglia,
Ralf Bender,
Fabrizio Finozzi,
Aliaksei Krukau,
Thorsten Naab,
Antti Rantala,
Matteo Frigo
Abstract:
We present SMART, a new 3D implementation of the Schwarzschild Method and its application to a triaxial N-body merger simulation. SMART fits full line-of-sight velocity distributions (LOSVDs) to determine the viewing angles, black hole, stellar and dark matter (DM) masses and the stellar orbit distribution of galaxies. Our model uses a 5D orbital starting space to ensure a representative set of st…
▽ More
We present SMART, a new 3D implementation of the Schwarzschild Method and its application to a triaxial N-body merger simulation. SMART fits full line-of-sight velocity distributions (LOSVDs) to determine the viewing angles, black hole, stellar and dark matter (DM) masses and the stellar orbit distribution of galaxies. Our model uses a 5D orbital starting space to ensure a representative set of stellar trajectories adaptable to the integrals-of-motion space and it is designed to deal with non-parametric stellar and DM densities. SMART's efficiency is demonstrated by application to a realistic N-body merger simulation including supermassive black holes which we model from five different projections. When providing the true viewing angles, 3D stellar luminosity profile and normalized DM halo, we can (i) reproduce the intrinsic velocity moments and anisotropy profile with a precision of ~1% and (ii) recover the black hole mass, stellar mass-to-light ratio and DM normalization to better than a few percent accuracy. This precision is smaller than the currently discussed differences between initial-stellar-mass functions and scatter in black hole scaling relations. Further tests with toy models suggest that the recovery of the anisotropy in triaxial galaxies is almost unique when the potential is known and full LOSVDs are fitted. We show that orbit models even allow the reconstruction of full intrinsic velocity distributions, which contain more information than the classical anisotropy parameter. Surprisingly, the orbit library for the analysed N-body simulation's gravitational potential contains orbits with net rotation around the intermediate axis that is stable over some Gyrs.
△ Less
Submitted 18 September, 2020;
originally announced September 2020.
-
Ionized gas kinematics of massive elliptical galaxies in CALIFA and in cosmological zoom-in simulations
Authors:
Jan Florian,
Bodo Ziegler,
Michaela Hirschmann,
Polychronis Papaderos,
Ena Choi,
Matteo Frigo,
Jean-Michel Gomes,
Rachel S. Somerville
Abstract:
(Abridged) We present an investigation of kinematical imprints of AGN feedback on the Warm Ionized gas Medium (WIM) of massive early-type galaxies (ETGs). To this end, we take a two-fold approach that involves a comparative analysis of Halpha velocity fields in 123 local ETGs from the CALIFA integral field spectroscopy survey with 20 simulated galaxies from high-resolution hydrodynamic cosmologica…
▽ More
(Abridged) We present an investigation of kinematical imprints of AGN feedback on the Warm Ionized gas Medium (WIM) of massive early-type galaxies (ETGs). To this end, we take a two-fold approach that involves a comparative analysis of Halpha velocity fields in 123 local ETGs from the CALIFA integral field spectroscopy survey with 20 simulated galaxies from high-resolution hydrodynamic cosmological SPHgal simulations. The latter were re-simulated for two modeling setups, one with and another without AGN feedback. In order to quantify the effects of AGN feedback on gas kinematics we measure three parameters that probe deviations from simple regular rotation using the kinemetry package. These indicators trace the possible presence of distinct kinematic components in Fourier space (k3,5/k1), variations in the radial profile of the kinematic major axis (sigma_PA), and offsets between the stellar and gas velocity fields (Delta Phi). These quantities are monitored in the simulations from a redshift 3 to 0.2 to assess the connection between black hole accretion history, stellar mass growth and kinematical perturbation of the WIM. Observed local massive galaxies show a broad range of irregularities, indicating disturbed warm gas motions, irrespective of being classified via diagnostic lines as AGN or not. Simulations of massive galaxies with AGN feedback generally exhibit higher irregularity parameters than without AGN feedback, more consistent with observations. Besides AGN feedback, other processes like major merger events or infalling gas clouds can lead to elevated irregularity parameters, but they are typically of shorter duration. More specifically, k3,5/k1 is most sensitive to AGN feedback, whereas Delta Phi is most strongly affected by gas infall.
△ Less
Submitted 13 December, 2019;
originally announced December 2019.
-
The simultaneous formation of cored, tangentially biased, and kinematically decoupled centers in massive early-type galaxies
Authors:
Antti Rantala,
Peter H. Johansson,
Thorsten Naab,
Jens Thomas,
Matteo Frigo
Abstract:
We study the impact of merging supermassive black holes (SMBHs) on the central regions of massive early-type galaxies (ETGs) using a series of merger simulations with varying mass ratios. The ETG models include realistic stellar and dark matter components and are evolved with the GADGET based regularized tree code KETJU. We show that observed key properties of the nuclear stellar populations of ma…
▽ More
We study the impact of merging supermassive black holes (SMBHs) on the central regions of massive early-type galaxies (ETGs) using a series of merger simulations with varying mass ratios. The ETG models include realistic stellar and dark matter components and are evolved with the GADGET based regularized tree code KETJU. We show that observed key properties of the nuclear stellar populations of massive ETGs, namely flat stellar density distributions (cores), tangentially biased velocity distributions and kinematically decoupled (counter-)rotation can naturally result from a single process $-$ the scouring by SMBHs. Major mergers with mass ratios of $q>1/3$ produce flat, tangentially biased cores with kinematically distinct components. These features originate from spin reversals during the decay of the SMBH orbits caused by gravitational torques from tidally expelled nuclear material. Minor mergers ($q\lesssim1/3$) on the other hand, form non-rotating cores and the tidal reversal becomes less important. Low-density stellar cores scoured in (multiple) minor mergers are less tangentially biased. This implies that the nuclear stellar properties of massive ETGs can be solely explained by stellar dynamical processes during their final assembly without any need for `feedback' from accreting black holes. We predict a strong correlation between decoupled cores, central anisotropy and merger history: decoupled cores form in binary mergers and we predict them to occur in elliptical galaxies with the strongest central anisotropy. Measurements of the central orbital structure are the key to understanding the number of mergers a given galaxy has experienced.
△ Less
Submitted 5 February, 2019; v1 submitted 6 December, 2018;
originally announced December 2018.
-
The impact of AGN on stellar kinematics and orbits in simulated massive galaxies
Authors:
Matteo Frigo,
Thorsten Naab,
Michaela Hirschmann,
Ena Choi,
Rachel Somerville,
Davor Krajnovic,
Romeel Davé,
Michele Cappellari
Abstract:
We present a series of 20 cosmological zoom simulations of the formation of massive galaxies with and without a model for AGN feedback. Differences in stellar population and kinematic properties are evaluated by constructing mock integral field unit (IFU) maps. The impact of the AGN is weak at high redshift when all systems are mostly fast-rotating and disc-like. After $z \sim 1$ the AGN simulatio…
▽ More
We present a series of 20 cosmological zoom simulations of the formation of massive galaxies with and without a model for AGN feedback. Differences in stellar population and kinematic properties are evaluated by constructing mock integral field unit (IFU) maps. The impact of the AGN is weak at high redshift when all systems are mostly fast-rotating and disc-like. After $z \sim 1$ the AGN simulations result in lower mass, older, less metal rich and slower rotating systems with less disky isophotes - in general agreement with observations. Two-dimensional kinematic maps of in-situ and accreted stars show that these differences result from reduced in-situ star formation due to AGN feedback. A full analysis of stellar orbits indicates that galaxies simulated with AGN are typically more triaxial and have higher fractions of x-tubes and box orbits and lower fractions of z-tubes. This trend can also be explained by reduced late in-situ star formation. We introduce a global parameter, $ξ_3$ , to characterise the anti-correlation between the third-order kinematic moment $h_3$ and the line-of-sight velocity ($v_{los}/σ$), and compare to ATLAS$^{3D}$ observations. The kinematic asymmetry parameter $ξ_3$ might be a useful diagnostic for large integral field surveys as it is a kinematic indicator for intrinsic shape and orbital content.
△ Less
Submitted 27 November, 2018;
originally announced November 2018.
-
The formation of extremely diffuse galaxy cores by merging supermassive black holes
Authors:
Antti Rantala,
Peter H. Johansson,
Thorsten Naab,
Jens Thomas,
Matteo Frigo
Abstract:
Given its velocity dispersion, the early-type galaxy NGC 1600 has an unusually massive ($M_\bullet = 1.7 \times 10^{10} M_\odot$) central supermassive black hole (SMBH), surrounded by a large core ($r_\mathrm{b} = 0.7$ kpc) with a tangentially biased stellar distribution. We present high-resolution equal-mass merger simulations including SMBHs to study the formation of such systems. The structural…
▽ More
Given its velocity dispersion, the early-type galaxy NGC 1600 has an unusually massive ($M_\bullet = 1.7 \times 10^{10} M_\odot$) central supermassive black hole (SMBH), surrounded by a large core ($r_\mathrm{b} = 0.7$ kpc) with a tangentially biased stellar distribution. We present high-resolution equal-mass merger simulations including SMBHs to study the formation of such systems. The structural parameters of the progenitor ellipticals were chosen to produce merger remnants resembling NGC 1600. We test initial stellar density slopes of $ρ\propto r^{-1}$ and $ρ\propto r^{-3/2}$ and vary the initial SMBH masses from $8.5 \times 10^8$ to $8.5 \times 10^9$ $M_\odot$. With increasing SMBH mass the merger remnants show a systematic decrease in central surface brightness, an increasing core size, and an increasingly tangentially biased central velocity anisotropy. Two-dimensional kinematic maps reveal decoupled, rotating core regions for the most massive SMBHs. The stellar cores form rapidly as the SMBHs become bound, while the velocity anisotropy develops more slowly after the SMBH binaries become hard. The simulated merger remnants follow distinct relations between the core radius and the sphere-of-influence, and the SMBH mass, similar to observed systems. We find a systematic change in the relations as a function of the progenitor density slope, and present a simple scouring model reproducing this behavior. Finally, we find the best agreement with NGC 1600 using SMBH masses totaling the observed value of $M_\bullet = 1.7 \times 10^{10} M_\odot$. In general, density slopes of $ρ\propto r^{-3/2}$ for the progenitor galaxies are strongly favored for the equal-mass merger scenario.
△ Less
Submitted 6 September, 2018; v1 submitted 25 May, 2018;
originally announced May 2018.
-
Time Delay Lens Modeling Challenge: I. Experimental Design
Authors:
Xuheng Ding,
Tommaso Treu,
Anowar J. Shajib,
Dandan Xu,
Geoff C. -F. Chen,
Anupreeta More,
Giulia Despali,
Matteo Frigo,
Christopher D. Fassnacht,
Daniel Gilman,
Stefan Hilbert,
Philip J. Marshall,
Dominique Sluse,
Simona Vegetti
Abstract:
Strong gravitational lenses with measured time delay are a powerful tool to measure cosmological parameters, especially the Hubble constant ($H_0$). Recent studies show that by combining just three multiply-imaged AGN systems, one can determine $H_0$ to 2.4% precision. Furthermore, the number of time-delay lens systems is growing rapidly, enabling the determination of $H_0$ to 1% precision in the…
▽ More
Strong gravitational lenses with measured time delay are a powerful tool to measure cosmological parameters, especially the Hubble constant ($H_0$). Recent studies show that by combining just three multiply-imaged AGN systems, one can determine $H_0$ to 2.4% precision. Furthermore, the number of time-delay lens systems is growing rapidly, enabling the determination of $H_0$ to 1% precision in the near future. However, as the precision increases it is important to ensure that systematic errors and biases remain subdominant. For this purpose, challenges with simulated datasets are a key component in this process. Following the experience of the past challenge on time delay, where it was shown that time delays can indeed be measured precisely and accurately at the sub-percent level, we now present the "Time Delay Lens Modeling Challenge" (TDLMC). The goal of this challenge is to assess the present capabilities of lens modeling codes and assumptions and test the level of accuracy of inferred cosmological parameters given realistic mock datasets. We invite scientists to model a set of simulated HST observations of 50 mock lens systems. The systems are organized in rungs, with the complexity and realism increasing going up the ladder. The goal of the challenge is to infer $H_0$ for each rung, given the HST images, the time delay, and stellar velocity dispersion of the deflector for a fixed background cosmology. The TDLMC challenge starts with the mock data release on 2018 January 8th. The deadline for blind submission is different for each rung. The deadline for Rung 0-1 is 2018 September 8; the deadline for Rung 2 is 2019 April 8 and the one for Rung 3 is 2019 September 8. This first paper gives an overview of the challenge including the data design, and a set of metrics to quantify the modeling performance and challenge details.
△ Less
Submitted 15 June, 2020; v1 submitted 4 January, 2018;
originally announced January 2018.
-
The fate of the Antennae galaxies
Authors:
Natalia Lahén,
Peter H. Johansson,
Antti Rantala,
Thorsten Naab,
Matteo Frigo
Abstract:
We present a high-resolution smoothed particle hydrodynamics simulation of the Antennae galaxies (NGC 4038/4039) and follow the evolution $3$ Gyrs beyond the final coalescence. The simulation includes metallicity dependent cooling, star formation, and both stellar feedback and chemical enrichment. The simulated best-match Antennae reproduces well both the observed morphology and the off-nuclear st…
▽ More
We present a high-resolution smoothed particle hydrodynamics simulation of the Antennae galaxies (NGC 4038/4039) and follow the evolution $3$ Gyrs beyond the final coalescence. The simulation includes metallicity dependent cooling, star formation, and both stellar feedback and chemical enrichment. The simulated best-match Antennae reproduces well both the observed morphology and the off-nuclear starburst. We also produce for the first time a simulated two-dimensional metallicity map of the Antennae and find good agreement with the observed metallicity of off-nuclear stellar clusters, however the nuclear metallicities are overproduced by $\sim 0.5$ dex. Using the radiative transfer code SKIRT we produce multi-wavelength observations of both the Antennae and the merger remnant. The $1$ Gyr old remnant is well fitted with a Sérsic profile of $n=7.07$, and with an $r$-band effective radius of $r_{\mathrm{e}}= 1.6$ kpc and velocity dispersion of $σ_{\mathrm{e}}=180$ km$/$s the remnant is located on the fundamental plane of early-type galaxies (ETGs). The initially blue Antennae remnant evolves onto the red sequence after $\sim 2.5$ Gyr of secular evolution. The remnant would be classified as a fast rotator, as the specific angular momentum evolves from $λ_{\mathrm{Re}}\approx0.11$ to $λ_{\mathrm{Re}}\approx0.14$ during its evolution. The remnant shows ordered rotation and a double peaked maximum in the mean 2D line-of-sight velocity. These kinematical features are relatively common among local ETGs and we specifically identify three local ETGs (NGC 3226, NGC 3379 and NGC 4494) in the ATLAS$^\mathrm{3D}$ sample, whose photometric and kinematic properties most resemble the Antennae remnant.
△ Less
Submitted 7 January, 2018; v1 submitted 31 August, 2017;
originally announced September 2017.
-
Dynamical masses and non-homology of massive elliptical galaxies grown by dry mergers
Authors:
Matteo Frigo,
Marc Balcells
Abstract:
We study whether dry merger-driven size growth of massive elliptical galaxies depends on their initial structural concentration, and analyse the validity of the homology hypothesis for virial mass determination in massive ellipticals grown by dry mergers. High-resolution simulations of a few realistic merger trees, starting with compact progenitors of different structural concentrations (Sérsic in…
▽ More
We study whether dry merger-driven size growth of massive elliptical galaxies depends on their initial structural concentration, and analyse the validity of the homology hypothesis for virial mass determination in massive ellipticals grown by dry mergers. High-resolution simulations of a few realistic merger trees, starting with compact progenitors of different structural concentrations (Sérsic indices n), show that galaxy growth has little dependence on the initial Sérsic index (larger n leads to slightly larger size growth), and depends more on other particulars of the merger history. We show that the deposition of accreted matter in the outer parts leads to a systematic and predictable breaking of the homology between remnants and progenitors, which we characterize through the evolution, during the course of the merger history, of virial coefficients K = GM/Re σ^2 associated to the most commonly-used dynamical and stellar mass parameters. The virial coefficient for the luminous mass, K , is about 50 per cent larger at the z = 2 start of the merger evolution than in z = 0 remnants. Ignoring virial evolution leads to biased virial mass estimates. We provide K corresponding to a variety of dynamical and stellar mass parameters, and provide recipes for the dynamical determination of galaxy masses. For massive, non-compact ellipticals, the popular expression M = 5 Re σ^2 / G underestimates the dynamical mass within the luminous body by factors of up to 4; it instead provides an approximation to the total stellar mass with smaller uncertainty than current stellar population models.
△ Less
Submitted 19 April, 2017; v1 submitted 5 November, 2016;
originally announced November 2016.