Jens Thomas

HST STIS spectroscopy of the triple nucleus of M31: two nested disks in keplerian rotation around a supermassive black hole

We present Hubble Space Telescope (HST) spectroscopy of the nucleus of M31 obtained with the Space TelescopeImagingSpectrograph(STIS).SpectrathatincludetheCaiiinfraredtriplet(k ’ 85008)seeonlythered giant stars in the double brightness peaks P1 and P2. In contrast, spectra taken atk ’ 3600 51008 are sensitive to thetinybluenucleusembeddedinP2,thelowersurfacebrightnessnucleusofthegalaxy.P2 hasaK-typespectrum, but we find that the blue nucleus has an A-type spectrum: it shows strong Balmer absorption lines. Hence, the blue nucleus is blue not because of AGN light but rather because it is dominated by hot stars. We show that the spectrum is well described by A0 giant stars, A0 dwarf stars, or a 200 Myr old, single-burst stellar population. White dwarfs, in contrast, cannot fit the blue nucleus spectrum. Given the small likelihood for stellar collisions, recent star formation appears to be the most plausible origin of the blue nucleus. In stellar population, size, and velocity dispersion, the blue nucleus is so different from P1 and P2 that we call it P3 and refer to the nucleus of M31 as triple. Because P2 and P3 have very different spectra, we can make a clean decomposition of the red and blue stars and hence measure the light distribution and kinematics of eachuncontaminated by the other. The line-of-sight velocity distributions of the red stars near P2 strengthen the support for Tremaine’s eccentric disk model. Their wings indicate the presence of stars with velocities of up to 1000 km s � 1 on the anti-P1 side of P2. The kinematics of P3 are consistent with a circular stellar disk in Keplerian rotation around a supermassive black hole.If the P3 diskis perfectlythin,thentheinclination anglei ’ 55 � isidentical withinthe errorsto theinclination of the eccentric disk models for P1+P2 by Peiris & Tremaine and by Salow & Statler. Both disks rotate in the same sense and are almost coplanar. The observed velocity dispersion of P3 is largely caused by blurred rotation and has a maximum value of � ¼ 1183 � 201 km s � 1 . This is much larger than the dispersion � ’ 250 km s � 1 of the red stars along the same line of sight and is the largest integrated velocity dispersion observed in any galaxy. The rotation curve of P3 is symmetric around its center. It reaches an observed velocity of V ¼ 618 � 81 km s � 1 at radius 0B05 ¼ 0:19 pc, where the observed velocity dispersion is � ¼ 674 � 95 km s � 1 . The corresponding circular rotation velocity at this radius is � 1700 km s � 1 . We therefore confirm earlier suggestions that the central dark object

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

Dynamical models for 17 early-type galaxies in the Coma cluster are presented. The galaxy sample consists of flattened, rotating as well as non-rotating early-types including cD and S0 galaxies with luminosities between MB = 18.79 and MB = 22.56. Kinematical long-slit observations cover at least the major and minor axis and extend to 1 4 reff. Axisymmetric Schwarzschild models are used to derive stellar mass-tolight ratios and dark halo parameters. In every galaxy the best fit with dark matter matches the data better than the best fit 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, independent of the assumed dark matter density profile. 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 coefficientH4. 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. This paper is aimed to set the basis for a subsequent detailed analysis of luminous and dark matter scaling relations, orbital dynamics and stellar populations.

Mapping stationary axisymmetric phase-space distribution functions by orbit libraries

This is the first of a series of papers dedicated to unveiling the mass composition and dynamical structure of a sample of flattened early-type galaxies in the Coma cluster. We describe our modifications to the Schwarzschild code of Richstone et al. Applying a Voronoi tessellation in the surface of section, we are able to assign accurate phase-space volumes to individual orbits and to reconstruct the full three-integral phase-space distribution function (DF) of any axisymmetric orbit library. Two types of tests have been performed to check the accuracy with which DFs can be represented by appropriate orbit libraries. First, by mapping DFs of spherical γ-models and flattened Plummer models onto the library, we show that the resulting line-of-sight velocity distributions and internal velocity moments of the library match those derived directly from the DF to a precision better than that of present-day observational errors. Secondly, by fitting libraries to the projected kinematics of the same DFs, we show that the DF reconstructed from the fitted library matches the input DF to a rms of about 15 per cent over a region in phase space covering 90 per cent of the mass of the library. The accuracy achieved allows us to implement effective entropy-based regularization to fit real, noisy and spatially incomplete data.

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

This is the second in a series of papers dedicated to unveiling 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 ≈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 r eff. The galaxy seems dark-matter dominated outside r > 2 r eff. 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. The achieved accuracy is �ϒ ≈ 0.5. Kinematically, NGC 4807 is characterized by mild radial anisotropy outside r > 0.5 r eff, becoming isotropic towards the centre. Our orbit models hint at either a distinct stellar component or weak triaxiality in the outer parts of the galaxy. Ke yw ords: stellar dynamics ‐ galaxies: elliptical and lenticular, cD ‐ galaxies: kinematics and dynamics ‐ galaxies: structure.

Relaxation and Stripping: The Evolution of Sizes, Dispersions and Dark Matter Fractions in Major and Minor Mergers of Elliptical Galaxies

We revisit collisionless major and minor mergers of spheroidal galaxies in the context of recent observational insights on compact massive early-type galaxies at high redshift and their rapid evolution on cosmological timescales. The simulations are performed as a series of mergers with mass-ratios of 1:1 and 1:10 for models representing pure bulges as well as bulges embedded in dark matter halos. For major and minor mergers, respectively, we identify and analyze two different processes, violent relaxation and stripping, leading to size evolution and a change of the dark matter fraction within the observable effective radius re. Violent relaxation - which is the dominant mixing process for major mergers but less important for minor mergers - scatters relatively more dark matter particles than bulge particles to radii r < re. Stripping in minor mergers assembles stellar satellite particles at large radii in halo dominated regions of the massive host. This strongly increases the size of the bulge into regions with higher dark matter fractions leaving the inner host structure almost unchanged. A factor of two mass increase by minor mergers increases the dark matter fraction within the effective radius by 80 per cent whereas relaxation in one equal-mass merger only leads to an increase of 25 percent. We present analytic corrections to simple one-component virial estimates for the evolution of the gravitational radii. These estimates are shown to underpredict the evolution of the effective radii for parabolic minor mergers of bulges embedded in massive dark matter halos. If such a two-component system grows by minor mergers alone its size growth, re / M � , reaches values of α � 2.4, significantly exceeding the simple theoretical limit of α = 2. For major mergers the sizes grow with α . 1. In addition, we discuss the velocity dispersion evolution and velocity anisotropy profiles. Our results indicate that minor mergers of galaxies embedded in massive dark matter halos provide a potential mechanism for explaining the rapid size growth and the build-up of massive elliptical systems predicting significant dark matter fractions and radially biased velocity dispersions at large radii.

Dynamical Masses of Early-Type Galaxies: A Comparison to Lensing Results and Implications for the Stellar Initial Mass Function and the Distribution of Dark Matter

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.

DARK MATTER SCALING RELATIONS AND THE ASSEMBLY EPOCH OF COMA EARLY-TYPE GALAXIES

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-type galaxies?those with old stellar populations?have similar halo core radii but more than two times larger asymptotic halo velocities. The average dark matter density inside 2 r eff 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.

Do black hole masses scale with classical bulge luminosities only? The case of the two composite pseudo-bulge galaxies NGC 3368 and NGC 3489

It is now well established that all galaxies with a massive bulge component harbour a central supermassive black hole (SMBH). The mass of the SMBH correlates with bulge properties such as the bulge mass and the velocity dispersion, which implies that the bulge and the central black hole of a galaxy have grown together during the formation process. As part of an investigation of the dependence of the SMBH mass on bulge types and formation mechanisms, we present measurements of SMBH masses in two pseudobulge galaxies. The spiral galaxy NGC 3368 is double-barred and hosts a large pseudobulge with a tiny classical bulge component at the very centre. The S0 galaxy NGC 3489 has only a weak large-scale bar, a small pseudobulge and a small classical bulge. Both galaxies show weak nuclear activity in the optical, indicative of the presence of a supermassive black hole. We present high resolution, adaptive-optics-assisted, near-infrared integral field data of these two galaxies, taken with SINFONI at the Very Large Telescope, and use axisymmetric orbit models to determine the masses of the SMBHs. The SMBH mass of NGC 3368, averaged over the four quadrants, ishM i = 7:5 10 6 M with an error of 1:5 10 6 M , which mostly comes from the non-axisymmetry in the data. For NGC 3489, a solution without black hole cannot be excluded when modelling the SINFONI data alone, but can be clearly ruled out when modelling a combination of SINFONI, OASIS and SAURON data, for which we obtain M = (6:00 +0:56 0:54 jstat 0:64jsys) 10 6 M . Although both galaxies seem to be consistent with the M - relation, at face value they do not agree with the relation between bulge magnitude and black hole mass when the total bulge magnitude (i.e., including both classical bulge and pseudobulge) is considered; the agreement is better when only the small classical bulge components are considered. However, taking into account the ageing of the stellar population could change this conclusion.

DWARF GALAXY DARK MATTER DENSITY PROFILES INFERRED FROM STELLAR AND GAS KINEMATICS

We present new constraints on the density profiles of dark matter (DM) halos in seven nearby dwarf galaxies from measurements of their integrated stellar light and gas kinematics. The gas kinematics of low-mass galaxies frequently suggest that they contain constant density DM cores, while N-body simulations instead predict a cuspy profile. We present a data set of high-resolution integral-field spectroscopy on seven galaxies and measure the stellar and gas kinematics simultaneously. Using Jeans modeling on our full sample, we examine whether gas kinematics in general produce shallower density profiles than are derived from the stars. Although two of the seven galaxies show some localized differences in their rotation curves between the two tracers, estimates of the central logarithmic slope of the DM density profile, γ, are generally robust. The mean and standard deviation of the logarithmic slope for the population are γ = 0.67 ± 0.10 when measured in the stars and γ = 0.58 ± 0.24 when measured in the gas. We also find that the halos are not under-concentrated at the radii of half their maximum velocities. Finally, we search for correlations of the DM density profile with stellar velocity anisotropy and other baryonic properties. Two popular mechanisms to explain cored DM halos are an exotic DM component or feedback models that strongly couple the energy of supernovae into repeatedly driving out gas and dynamically heating the DM halos. While such models do not yet have falsifiable predictions that we can measure, we investigate correlations that may eventually be used to test models. We do not find a secondary parameter that strongly correlates with the central DM density slope, but we do find some weak correlations. The central DM density slope weakly correlates with the abundance of α elements in the stellar population, anti-correlates with H I fraction, and anti-correlates with vertical orbital anisotropy. We expect, if anything, the opposite of these three trends for feedback models. Determining the importance of these correlations will require further model developments and larger observational samples.

The MASSIVE survey. I. A volume-limited integral-field spectroscopic study of the most massive early-type galaxies within 108 Mpc

Massive early-type galaxies represent the modern-day remnants of the earliest major star formation episodes in the history of the universe. These galaxies are central to our understanding of the evolution of cosmic structure, stellar populations, and supermassive black holes, but the details of their complex formation histories remain uncertain. To address this situation, we have initiated the MASSIVE Survey, a volume-limited, multi-wavelength, integral-field spectroscopic (IFS) and photometric survey of the structure and dynamics of the ~100 most massive early-type galaxies within a distance of 108 Mpc. This survey probes a stellar mass range M* > 10^{11.5} Msun and diverse galaxy environments that have not been systematically studied to date. Our wide-field IFS data cover about two effective radii of individual galaxies, and for a subset of them, we are acquiring additional IFS observations on sub-arcsecond scales with adaptive optics. We are also acquiring deep K-band imaging to trace the extended halos of the galaxies and measure accurate total magnitudes. Dynamical orbit modeling of the combined data will allow us to simultaneously determine the stellar, black hole, and dark matter halo masses. The primary goals of the project are to constrain the black hole scaling relations at high masses, investigate systematically the stellar initial mass function and dark matter distribution in massive galaxies, and probe the late-time assembly of ellipticals through stellar population and kinematical gradients. In this paper, we describe the MASSIVE sample selection, discuss the distinct demographics and structural and environmental properties of the selected galaxies, and provide an overview of our basic observational program, science goals and early survey results.