Source author record

Ortwin Gerhard

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

53works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

53 published item(s)

preprint2026arXiv

Extragalactic planetary nebulae -- tracers of kinematics and stellar populations out to 100 Mpc

Extragalactic planetary nebulae (xPNe) in galaxies beyond the Local Universe serve as discrete tracers for studying the element abundances and kinematics of galaxies covering a wide range of morphologies and masses at a variety of angular distances, from the centre well out into their haloes. They are direct stellar probes to identify the galaxy progenitors of haloes and the intracluster light. Even with new facilities, reaching higher angular resolution and sensitivity, xPNe are the only stellar tracers that can be directly and singularly detected and characterised out to 100 Mpc distance, making them crucial for tracing halo and intracluster light assembly. New wide-field spectroscopic instruments at 10+meter-class telescopes will allow the unprecedented characterisation of xPN populations from galaxy centres to their diffuse outskirts. Panoramic integral-field spectroscopy will enable the simultaneous study of xPN and stellar population properties, establishing their use as age- and metallicity-tracers while also improving post-AGB stellar evolution models.

preprint2022arXiv

Large-scale Hydrodynamical Shocks as the Smoking Gun Evidence for a Bar in M31

The formation and evolutionary history of M31 are closely related to its dynamical structures, which remain unclear due to its high inclination. Gas kinematics could provide crucial evidence for the existence of a rotating bar in M31. Using the position-velocity diagram of [OIII] and HI, we are able to identify clear sharp velocity jump (shock) features with a typical amplitude over 100 km/s in the central region of M31 (4.6 kpc X 2.3 kpc, or 20 arcmin X 10 arcmin). We also simulate gas morphology and kinematics in barred M31 potentials and find that the bar-induced shocks can produce velocity jumps similar to those in [OIII]. The identified shock features in both [OIII] and HI are broadly consistent, and they are found mainly on the leading sides of the bar/bulge, following a hallmark pattern expected from the bar-driven gas inflow. Shock features on the far side of the disk are clearer than those on the near side, possibly due to limited data coverage on the near side, as well as obscuration by the warped gas and dust layers. Further hydrodynamical simulations with more sophisticated physics are desired to fully understand the observed gas features and to better constrain the parameters of the bar in M31.

preprint2022arXiv

The Pattern Speed of the Milky Way Bar/Bulge from VIRAC & Gaia

We compare distance resolved, absolute proper motions in the Milky Way bar/bulge region to a grid of made-to-measure dynamical models with well defined pattern speeds. The data are obtained by combining the relative VVV Infrared Astrometric Catalog v1 proper motions with the Gaia DR2 absolute reference frame. We undertake a comprehensive analysis of the various errors in our comparison, from both the data and the models, and allow for additional, unknown, contributions by using an outlier-tolerant likelihood function to evaluate the best fitting model. We quantify systematic effects such as the region of data included in the comparison, the possible overlap from spiral arms, and the choice of synthetic luminosity function and bar angle used to predict the data from the models. Resulting variations in the best-fit parameters are included in their final errors. We thus measure the bar pattern speed to be $Ω_b=33.29\pm1.81 \mathrm{km \, s^{-1} \, kpc^{-1}}$ and the azimuthal solar velocity to be $V_{ϕ,\,\odot} = 251.31 \pm 1.95 \mathrm{km \, s^{-1}}$. These values, when combined with recent measurements of the Galactic rotation curve, yield the distance of corotation, 6.5 < $R_{CR}$ [kpc] < 7.5, the outer Lindblad resonance (OLR), 10.7 < $R_{OLR}$ [kpc] < 12.4, and the higher order, m = 4, OLR, 8.7 < $R_{{OLR}_4}$ [kpc] < 10.0. The measured pattern speed provides strong evidence for the "long-slow" bar scenario.

preprint2022arXiv

The stellar mass distribution of the Milky Way's bar: an analytic model

We present an analytic model of the stellar mass distribution of the Milky Way bar. The model is obtained by fitting a multi-component parametric density distribution to a made-to-measure N-body model of Portail et al., constructed to match a variety of density and kinematics observational data. The analytic model reproduces in detail the 3D density distribution of the N-body bar including the X-shape. The model and the gravitational potential it generates are available as part of the software package AGAMA for galactic dynamics, and can be readily used for orbit integrations, hydrodynamical simulations or other applications.

preprint2021arXiv

Chemo-kinematics of the Milky Way spiral arms and bar resonances: connection to ridges and moving groups in the Solar vicinity

Using of a new N-body simulation, as well as Gaia DR2/EDR3 data complemented by data from the Galah, APOGEE, LAMOST surveys, we explore the possible link between the MW spiral arms, RVphi-ridges and moving groups in local UV space. We show that the tightly wound main spiral arms in the N-body model can be successfully identified using overdensities in angular momentum (AM) or guiding space, as well as in the distribution of dynamically cold stars close to their guiding centers. Stars in the AM overdensities that travel over many kpc in radius trace extended density ridges in RVphi space and overdensities in the UV plane of a SNd-like region, similar to those observed in the Gaia data. Similarly, the AM space of the MW contains several overdensities which correlate with a wave-like radial velocity pattern; this pattern is also reproduced by stars well beyond the SNd. We find that the fraction of Gaia stars located near their guiding centers shows three large-scale structures that approximately coincide with the MW spiral arms traced by distributions of maser sources in the Sagittarius, Local, and Perseus arms. Similar to the model, the stars in the AM overdensities follow the main RVphi ridges. When these ridges cross the SNd, they can be matched with the main UV features. Thus we suggest that the Hat is the inner tail of the Perseus arm, one of the Hercules components is the Sagittarius arm and the Arcturus stream is likely to be the outermost tail of the Scutum-Centaurus arm. Based on previous work, the bar corotation is suggested to coincide with the Hercules stream, and the OLR with the Sirius stream. The latter is supported by a sharp decrease of the mean metallicity beyond the Sirius stream, which is an expected behaviour of the OLR. In various phase-space coordinates the AM overdensities stars have systematically higher mean metallicity...(Abridged)

preprint2021arXiv

Gas Dynamics in the Galaxy: Total Mass Distribution and the Bar Pattern Speed

Gas morphology and kinematics in the Milky Way contain key information for understanding the formation and evolution of our Galaxy. We present a high resolution hydrodynamical simulation based on a realistic barred Milky Way potential constrained by recent observations. Our model can reproduce most features in the observed longitude-velocity diagram, including the Central Molecular Zone, the Near and Far 3-kpc arms, the Molecular Ring, and the spiral arm tangents. It can also explain the non-circular motions of masers obtained by the recent BeSSeL2 survey. The central gas kinematics are consistent with a mass of $6.9\times10^8\; {\rm M}_{\odot}$ in the Nuclear Stellar Disk. Our model predicts the formation of an elliptical gaseous ring surrounding the bar, which is composed of the 3-kpc arms, Norma arm, and the bar-spiral interfaces. This ring is similar to those "inner" rings in some Milky Way analogs with a boxy/peanut-shaped bulge. The kinematics of gas near the solar neighbourhood are governed by the Local arm, which is induced by the four major stellar spiral arms. The bar pattern speed constrained by our gas model is $37.5-40\; {\rm km}\;{\rm s}^{-1}\;{\rm kpc}^{-1}$, corresponding to a corotation radius of $R_{\rm CR}=6.0-6.4\;{\rm kpc}$. The rotation curve of our model rises gently within the central $\sim5\;{\rm kpc}$, which is significantly less steep than those predicted by modern zoom-in cosmological simulations such as Auriga.

preprint2020arXiv

The Sixth Data Release of the Radial Velocity Experiment (RAVE) -- I: Survey Description, Spectra and Radial Velocities

The Radial Velocity Experiment (RAVE) is a magnitude-limited (9<I<12) spectroscopic survey of Galactic stars randomly selected in the southern hemisphere. The RAVE medium-resolution spectra (R~7500) cover the Ca-triplet region (8410-8795A). The 6th and final data release (DR6 or FDR) is based on 518387 observations of 451783 unique stars. RAVE observations were taken between 12 April 2003 and 4 April 2013. Here we present the genesis, setup and data reduction of RAVE as well as wavelength-calibrated and flux-normalized spectra and error spectra for all observations in RAVE DR6. Furthermore, we present derived spectral classification and radial velocities for the RAVE targets, complemented by cross matches with Gaia DR2 and other relevant catalogs. A comparison between internal error estimates, variances derived from stars with more than one observing epoch and a comparison with radial velocities of Gaia DR2 reveals consistently that 68% of the objects have a velocity accuracy better than 1.4 km/s, while 95% of the objects have radial velocities better than 4.0 km/s. Stellar atmospheric parameters, abundances and distances are presented in subsequent publication. The data can be accessed via the RAVE Web (http://rave-survey.org) or the Vizier database.

preprint2020arXiv

The Sixth Data Release of the Radial Velocity Experiment (RAVE) -- II: Stellar Atmospheric Parameters, Chemical Abundances and Distances

We present part 2 of the 6th and final Data Release (DR6 or FDR) of the Radial Velocity Experiment (RAVE), a magnitude-limited (9<I<12) spectroscopic survey of Galactic stars randomly selected in the southern hemisphere. The RAVE medium-resolution spectra (R~7500) cover the Ca-triplet region (8410-8795A) and span the complete time frame from the start of RAVE observations on 12 April 2003 to their completion on 4 April 2013. In the second of two publications, we present the data products derived from 518387 observations of 451783 unique stars using a suite of advanced reduction pipelines focussing on stellar atmospheric parameters, in particular purely spectroscopically derived stellar atmospheric parameters (Teff, log(g), and the overall metallicity), enhanced stellar atmospheric parameters inferred via a Bayesian pipeline using Gaia DR2 astrometric priors, and asteroseismically calibrated stellar atmospheric parameters for giant stars based on asteroseismic observations for 699 K2 stars. In addition, we provide abundances of the elements Fe, Al, and Ni, as well as an overall [alpha/Fe] ratio obtained using a new pipeline based on the GAUGUIN optimization method that is able to deal with variable signal-to-noise ratios. The RAVE DR6 catalogs are cross matched with relevant astrometric and photometric catalogs, and are complemented by orbital parameters and effective temperatures based on the infrared flux method. The data can be accessed via the RAVE Web site (http://rave-survey.org) or the Vizier database.

preprint2020arXiv

Three Mechanisms for Bar Thickening

We present simulations of bar-unstable stellar discs in which the bars thicken into box/peanut shapes. Detailed analysis of the evolution of each model revealed three different mechanisms for thickening the bars. The first mechanism is the well-known buckling instability, the second is the vertical excitation of bar orbits by passage through the 2:1 vertical resonance, and the third is a gradually increasing fraction of bar orbits trapped into this resonance. Since bars in many galaxies may have formed and thickened long ago, we have examined the models for fossil evidence in the velocity distribution of stars in the bar, finding a diagnostic to discriminate between a bar that had buckled from the other two mechanisms.

preprint2019arXiv

Elliptical galaxies/stellar halos connection

Cosmological simulations predict that early-type galaxies (ETGs) are the results of extended mass accretion histories. The latter are characterized by different numbers of mergers, mergers mass ratios and gas fractions, and timing. Depending on the sequence and nature of these mergers that follow the first phase of the in-situ star formation, these accretion histories may lead to ETGs that have low or high mass halos, and that rotate fast or slow. Since the stellar halos maintain the fossil records of the events that led to their formation, a discontinuity may be in place between the inner regions of ETGs and their outer halos, because the time required for the halos stars to exchange their energies and momenta is very long compared with the age of these systems. Exquisite deep photometry and extended spectroscopy for significant samples of ETGs are then used to quantify the occurrence and significance of such a transition in the galaxies structural and kinematical parameters. Once this transition radius is measured, its dependency with the effective radius of the galaxies light distribution and total stellar masses can be investigated. Such correlations can then be compared with the predictions of accreted, i.e. ex-situ vs. in-situ components from cosmological simulations to validate such models.

preprint2016arXiv

Dynamical modelling of the Galactic bulge and bar: the Milky Way's bar pattern speed, stellar, and dark matter mass distribution

We construct a large set of dynamical models of the galactic bulge, bar and inner disk using the Made-to-Measure method. Our models are constrained to match the red clump giant density from a combination of the VVV, UKIDSS and 2MASS infrared surveys together with stellar kinematics in the bulge from the BRAVA and OGLE surveys, and in the entire bar region from the ARGOS survey. We are able to recover the bar pattern speed and the stellar and dark matter mass distributions in the bar region, thus recovering the entire galactic effective potential. We find a bar pattern speed of $39.0 \pm 3.5 \,\rm{km\,s^{-1}\,kpc^{-1}}$, placing the bar corotation radius at $6.1 \pm 0.5 \rm{kpc}$ and making the Milky Way bar a typical fast rotator. We evaluate the stellar mass of the long bar and bulge structure to be $M_{\rm{bar/bulge}} = 1.88 \pm 0.12 \times 10^{10} \, \rm{M}_{\odot}$, larger than the mass of disk in the bar region, $M_{\rm{inner\ disk}} = 1.29\pm0.12 \times 10^{10} \, \rm{M}_{\odot}$. The total dynamical mass in the bulge volume is $1.85\pm0.05\times 10^{10} \, \rm{M}_{\odot}$. Thanks to more extended kinematic data sets and recent measurement of the bulge IMF our models have a low dark matter fraction in the bulge of $17\%\pm2\%$. We find a dark matter density profile which flattens to a shallow cusp or core in the bulge region. Finally, we find dynamical evidence for an extra central mass of $\sim0.2\times10^{10} \,\rm{M}_{\odot}$, probably in a nuclear disk or disky pseudobulge.

preprint2016arXiv

Gas dynamics in the Milky Way: a low pattern speed model

We present gas flow models for the Milky Way based on high-resolution grid-based hydrodynamical simulations. The basic galactic potential we use is from a N-body model constrained by the density of red clump giants in the Galactic bulge. We augment this potential with a nuclear bulge, two pairs of spiral arms and additional mass at the bar end to represent the long bar component. With this combined model we can reproduce many features in the observed ($l,v$) diagram with a bar pattern speed of $33\; {\rm km}\;{\rm s}^{-1}\;{\rm kpc}^{-1}$ and a spiral pattern speed of $23\; {\rm km}\;{\rm s}^{-1}\;{\rm kpc}^{-1}$. The shape and kinematics of the nuclear ring, Bania's Clump 2, the Connecting arm, the Near and Far 3-kpc arms, the Molecular Ring, and the spiral arm tangent points in our simulations are comparable to those in the observations. Our results imply that a low pattern speed model for the bar in our Milky Way reproduces the observations for a suitable Galactic potential. Our best model gives a better match to the ($l,v$) diagram than previous high pattern speed hydrodynamical simulations.

preprint2016arXiv

MOA-II Galactic Microlensing Constraints: The Inner Milky Way has a Low Dark Matter Fraction and a Near Maximal Disk

Microlensing provides a unique tool to break the stellar to dark matter degeneracy in the inner Milky Way. We combine N-body dynamical models fitted to the Milky Way's Boxy/Peanut bulge with exponential disk models outside this, and compute the microlensing properties. Considering the range of models consistent with the revised MOA-II data, we find low dark matter fractions in the inner Galaxy: at the peak of their stellar rotation curve a fraction $f_v=(0.88\pm0.07)$ of the circular velocity is baryonic (at $1σ$, $f_v > 0.72$ at $2σ$). These results are in agreement with constraints from the EROS-II microlensing survey of brighter resolved stars, where we find $f_v=(0.9\pm0.1)$ at $1σ$. Our fiducial model of a disk with scale length 2.6kpc, and a bulge with a low dark matter fraction of 12%, agrees with both the revised MOA-II and EROS-II microlensing data. The required baryonic fractions, and the resultant low contribution from dark matter, are consistent with the NFW profiles produced by dissipationless cosmological simulations in Milky Way mass galaxies. They are also consistent with recent prescriptions for the mild adiabatic contraction of Milky Way mass haloes without the need for strong feedback, but there is some tension with recent measurements of the local dark matter density. Microlensing optical depths from the larger OGLE-III sample could improve these constraints further when available.

preprint2016arXiv

Testing General Relativity with the Shadow Size of Sgr A*

In general relativity, the angular radius of the shadow of a black hole is primarily determined by its mass-to-distance ratio and depends only weakly on its spin and inclination. If general relativity is violated, however, the shadow size may also depend strongly on parametric deviations from the Kerr metric. Based on a reconstructed image of Sagittarius A* (Sgr A*) from a simulated one-day observing run of a seven-station Event Horizon Telescope (EHT) array, we employ a Markov chain Monte Carlo algorithm to demonstrate that such an observation can measure the angular radius of the shadow of Sgr A* with an uncertainty of ~1.5 uas (6%). We show that existing mass and distance measurements can be improved significantly when combined with upcoming EHT measurements of the shadow size and that tight constraints on potential deviations from the Kerr metric can be obtained.

preprint2016arXiv

The Stellar Halo in the Inner Milky Way: Predicted Shape and Kinematics

We have used N-body simulations for the Milky Way to investigate the kinematic and structural properties of the old metal-poor stellar halo in the barred inner region of the Galaxy. We find that the extrapolation of the density distribution for bulge RR Lyrae stars, $ρ\sim r^{-3}$, approximately matches the number density of RR Lyrae in the nearby stellar halo. We follow the evolution of such a tracer population through the formation and evolution of the bar and box/peanut bulge in the N-body model. We find that its density distribution changes from oblate to triaxial, and that it acquires slow rotation in agreement with recent measurements. The maximum radial velocity is $\sim15-25$ km/s at $| l|\!=10^\circ-30^\circ$, and the velocity dispersion is $\sim120$ km/s. Even though the simulated metal-poor halo in the bulge has a barred shape, just $12\%$ of the orbits follow the bar, and it does not trace the peanut/X structure. With these properties, the RR Lyrae population in the Galactic bulge is consistent with being the inward extension of the Galactic metal-poor stellar halo.

preprint2015arXiv

26Al kinematics: superbubbles following the spiral arms? Constraints from the statistics of star clusters and HI supershells

High energy resolution spectroscopy of the 1.8 MeV radioactive decay line of 26Al with the SPI instrument on board the INTEGRAL satellite has recently revealed that diffuse 26Al has large velocities in comparison to other components of the interstellar medium in the Milky Way. 26Al shows Galactic rotation in the same sense as the stars and other gas tracers, but reaches excess velocities up to 300 km/s. We investigate if this result can be understood in the context of superbubbles, taking into account the statistics of young star clusters and H I supershells, as well as the association of young star clusters with spiral arms. We derive energy output and 26Al mass of star clusters as a function of the cluster mass via population synthesis from stellar evolutionary tracks of massive stars. [...] We link this to the size distribution of HI supershells and assess the properties of likely 26Al-carrying superbubbles. 26Al is produced by star clusters of all masses above about 200 solar masses, roughly equally contributed over a logarithmic star cluster mass scale, and strongly linked to the injection of feedback energy. The observed superbubble size distribution cannot be related to the star cluster mass function in a straight forward manner. In order to avoid that the added volume of all superbubbles exceeds the volume of the Milky Way, individual superbubbles have to merge frequently. If any two superbubbles merge, or if 26Al is injected off-centre in a bigger HI supershell we expect the hot 26Al-carrying gas to obtain velocities of the order of the typical sound speed in superbubbles, about 300 km/s before decay. [...]

preprint2015arXiv

Dust within the old nuclear star cluster in the Milky Way

The mean absolute extinction towards the central parsec of the Milky Way is A_K~3 mag, including both foreground and Galactic center dust. Here we present a measurement of dust extinction within the Galactic old nuclear star cluster (NSC), based on combining differential extinctions of NSC stars with their u_l proper motions along Galactic longitude. Extinction within the NSC preferentially affects stars at its far side, and because the NSC rotates, this causes higher extinctions for NSC stars with negative u_l, as well as an asymmetry in the u_l-histograms. We model these effects using an axisymmetric dynamical model of the NSC in combination with simple models for the dust distribution. Comparing the predicted asymmetry to data for ~7100 stars in several NSC fields, we find that dust associated with the Galactic center mini-spiral with extinction A_K~=0.15-0.8 mag explains most of the data. The largest extinction A_K~=0.8 mag is found in the region of the Western arm of the mini-spiral. Comparing with total A_K determined from stellar colors, we determine the extinction in front of the NSC. Finally, we estimate that for a typical extinction of A_K~=0.4 the statistical parallax of the NSC changes by ~0.4%.

preprint2015arXiv

Planetary Nebulae and their parent stellar populations. Tracing the mass assembly of M87 and Intracluster light in the Virgo cluster core

The diffuse extended outer regions of galaxies are hard to study because they are faint, with typical surface brightness of 1% of the dark night sky. We can tackle this problem by using resolved star tracers which remain visible at large distances from the galaxy centres. This article describes the use of Planetary Nebulae as tracers and the calibration of their properties as indicators of the star formation history, mean age and metallicity of the parent stars in the Milky Way and Local Group galaxies . We then report on the results from a deep, extended, planetary nebulae survey in a 0.5 sqdeg region centred on the brightest cluster galaxy NGC 4486 (M87) in the Virgo cluster core, carried out with SuprimeCam@Subaru and FLAMES-GIRAFFE@VLT. Two PN populations are identified out to 150 kpc distance from the centre of M87. One population is associated with the M87 halo and the second one with the intracluster light in the Virgo cluster core. They have different line-of-sight velocity and spatial distributions, as well as different planetary nebulae specific frequencies and luminosity functions. The intracluster planetary nebulae in the surveyed region correspond to a luminosity of four times the luminosity of the Large Magellanic Cloud. The M87 halo planetary nebulae trace an older, more metal-rich, parent stellar population. A substructure detected in the projected phase-space of the line-of-sight velocity vs. major axis distance for the M87 halo planetary nebulae provides evidence for the recent accretion event of a satellite galaxy with luminosity twice that of M33. The satellite stars were tidally stripped about 1 Gyr ago, and reached apocenter at a major axis distance of 60-90 kpc from the centre of M87. The M87 halo is still growing significantly at the distances where the substructure is detected.

preprint2015arXiv

Spin-up of massive classical bulges during secular evolution

Classical bulges in spiral galaxies are known to rotate but the origin of this observed rotational motion is not well understood. It has been shown recently that a low-mass classical bulge (ClB) in a barred galaxy can acquire rotation from absorbing a significant fraction of the angular momentum emitted by the bar. Our aim here is to investigate whether bars can spin up also more massive ClBs during the secular evolution of the bar, and to study the kinematics and dynamics of these ClBs. We use a set of self-consistent N-body simulations to study the interaction of ClBs with a bar that forms self-consistently in the disk. We use orbital spectral analysis to investigate the angular momentum gain by the classical bulge stars. We show that the ClBs gain, on average, about 2 - 6% of the disk's initial angular momentum within the bar region. Most of this angular momentum gain occurs via low-order resonances, particularly 5:2 resonant orbits. A density wake forms in the ClB which corotates and aligns with the bar at the end of the evolution. The spin-up process creates a characteristic linear rotation profile and mild tangential anisotropy in the ClB. The induced rotation is small in the centre but significant beyond $\sim2$ bulge half mass radii, where it leads to mass-weighted $V/σ\sim 0.2$, and reaches a local $V_{max}/σ_{in} \sim 0.5$ at around the scale of the bar. The resulting $V/σ$ is tightly correlated with the ratio of the bulge size to the bar size. In all models, a box/peanut bulge forms suggesting that composite bulges may be common. Bar-bulge resonant interaction in barred galaxies can provide some spin up of massive ClBs, but the process appears to be less efficient than for low-mass ClBs. Further angular momentum transfer due to nuclear bars or gas inflow would be required to explain the observed rotation if it is not primordial.

preprint2015arXiv

The outer regions of the giant Virgo galaxy M87. Kinematic separation of stellar halo and intracluster light

We present a spectroscopic study of 287 Planetary Nebulas (PNs) in a total area of ~0.4 deg^2 around the BCG M87 in Virgo A. With these data we can distinguish the stellar halo from the co-spatial intracluster light (ICL). PNs were identified from their narrow and symmetric redshifted lambda 5007\4959 Angstrom [OIII] emission lines, and the absence of significant continuum. We implement a robust technique to measure the halo velocity dispersion from the projected phase-space to identify PNs associated with the M87 halo and ICL. The velocity distribution of the spectroscopically confirmed PNs is bimodal, containing a narrow component centred on the systemic velocity of the BCG and an off-centred broader component, that we identify as halo and ICL, respectively. Halo and ICPN have different spatial distributions: the halo PNs follow the galaxy's light, whereas the ICPNs are characterised by a shallower power-law profile. The composite PN number density profile shows the superposition of different PN populations associated with the M87 halo and the ICL, characterised by different PN alpha-parameters, the ICL contributing ~3 times more PNs per unit light. Down to m_5007=28.8, the M87 halo PN luminosity function (PNLF) has a steeper slope towards faint magnitudes than the IC PNLF, and both are steeper than the standard PNLF for the M31 bulge. Moreover, the IC PNLF has a dip at ~1-1.5 mag fainter than the bright cutoff, reminiscent of the PNLFs of systems with extended star formation history. The M87 halo and the Virgo ICL are dynamically distinct components with different density profiles and velocity distribution. The different alpha values and PNLF shapes of the halo and ICL indicate distinct parent stellar populations, consistent with the existence of a gradient towards bluer colours at large radii. These results reflect the hierarchical build-up of the Virgo cluster.

preprint2015arXiv

The Structure of the Milky Way's Bar Outside the Bulge

While it is incontrovertible that the inner Galaxy contains a bar, its structure near the Galactic plane has remained uncertain, where extinction from intervening dust is greatest. We investigate here the Galactic bar outside the bulge, the long bar, using red clump giant (RCG) stars from UKIDSS, 2MASS, VVV, and GLIMPSE. We match and combine these surveys to investigate a wide area in latitude and longitude, |b|<9deg and |l|<40deg. We find: (1) The bar extends to l~25deg at |b|~5deg from the Galactic plane, and to l~30deg at lower latitudes. (2) The long bar has an angle to the line-of-sight in the range (28-33)deg, consistent with studies of the bulge at |l|<10deg. (3) The scale-height of RCG stars smoothly transitions from the bulge to the thinner long bar. (4) There is evidence for two scale heights in the long bar. We find a ~180pc thin bar component reminiscent of the old thin disk near the sun, and a ~45pc super-thin bar component which exists predominantly towards the bar end. (5) Constructing parametric models for the RC magnitude distributions, we find a bar half length of 5.0+-0.2kpc for the 2-component bar, and 4.6+-0.3kpc for the thin bar component alone. We conclude that the Milky Way contains a central box/peanut bulge which is the vertical extension of a longer, flatter bar, similar as seen in both external galaxies and N-body models.

preprint2014arXiv

Formation Models of the Galactic Bulge

The Galactic bulge is now considered to be the inner three-dimensional part of the Milky Way's bar. It has a peanut shape and is characterized by cylindrical rotation. In N-body simulations, box/peanut bulges arise from disks through bar and buckling instabilities. Models of this kind explain much of the structure and kinematics of the Galactic bulge and, in principle, also its vertical metallicity gradient. Cosmological disk galaxy formation models with high resolution and improved feedback models are now able to generate late-type disk galaxies with disk-like or barred bulges. These bulges often contain an early collapse stellar population and a population driven by later disk instabilities. Due to the inside-out disk formation, these bulges can be predominantly old, similar to the Milky Way bulge.

preprint2014arXiv

The old nuclear star cluster in the Milky Way: dynamics, mass, statistical parallax, and black hole mass

We derive new constraints on the mass, rotation, orbit structure and statistical parallax of the Galactic old nuclear star cluster (NSC) and the mass of the supermassive black hole. We combine star counts and kinematic data from Fritz et al (2014), including 2'500 line-of-sight velocities and 10'000 proper motions. We show that the difference between the proper motion dispersions sigma_l and sigma_b cannot be explained by rotation, but is a consequence of the flattening of the NSC. We fit the surface density distribution of stars in the central 1000" by a spheroidal cluster with scale ~100" and a much larger nuclear disk component. We compute the two-integral distribution function f(E,Lz) for this density model, and add rotation self-consistently. We find that: (i) The orbit structure of the f(E,Lz) gives an excellent match to the observed velocity dispersion profiles as well as the proper motion and line-of-sight velocity histograms, including the double-peak in the v_l-histograms. (ii) This requires an axial ratio of q= 0.73+-0.04 for r<70" consistent with our determination from star counts. (iii) The NSC is approximately described by an isotropic rotator model. (iv) Using the corresponding Jeans equations to fit the proper motion and line-of-sight velocity dispersions, we obtain best estimates for the NSC mass, black hole mass, and distance M*(r<100")=(8.94+-0.31|stat+-0.9|syst)x10^6Msun, Mbh=(3.86+-0.14|stat+-0.4|syst)x10^6Msun, and R0=8.27+-0.09|stat+-0.1|syst kpc, where the systematic errors estimate additional uncertainties in the dynamical modeling. (v) The combination of the cluster dynamics with the S-star orbits around Sgr A* strongly reduces the degeneracy between black hole mass and Galactic centre distance present in previous S-star studies. A joint statistical analysis with the results of Gillessen et al (2009) gives Mbh=(4.23+-0.14)x10^6Msun and R0=8.33+-0.11kpc.

preprint2013arXiv

Elliptical galaxies with rapidly decreasing velocity dispersion profiles: NMAGIC models and dark halo parameter estimates for NGC 4494

NGC 4494 is one of several intermediate-luminosity elliptical galaxies inferred to have an unusually diffuse dark matter halo. We use the chi^2-made-to-measure particle code NMAGIC to construct axisymmetric models of NGC 4494 from photometric and various kinematic data. The extended kinematics include light spectra in multiple slitlets out to 3.5 R_e, and hundreds of planetary nebulae velocities out to ~7 R_e, thus allowing us to probe the dark matter content and orbital structure in the halo. We use Monte Carlo simulations to estimate confidence boundaries for the halo parameters, given our data and modelling set-up. We find that the true potential of the dark matter halo is recovered within Delta G (merit function)<26 (Delta chi^2<59) at 70% confidence level (C.L.), and within Delta G<32 (Delta chi^2<70) at 90% C.L.. These numbers are much larger than the usually assumed Delta chi^2=2.3 (4.6) for 70% (90%) C.L. for two free parameters, perhaps case-dependent, but calling into question the general validity of the standard assumptions used for halo and black hole mass determinations. The best-fitting models for NGC 4494 have a dark matter fraction of about 0.6\pm0.1 at 5R_e (70% C.L.), and are embedded in a dark matter halo with circular velocity ~200 km/s. The total circular velocity curve (CVC) is approximately flat at v_c=220 km/s outside ~0.5R_e. The orbital anisotropy of the stars is moderately radial. These results are independent of the assumed inclination of the galaxy, and edge-on models are preferred. Comparing with the halos of NGC 3379 and NGC 4697, whose velocity dispersion profiles also decrease rapidly from the center outwards, the outer CVCs and dark matter halos are quite similar. NGC 4494 shows a particularly high dark matter fraction inside ~3R_e, and a strong concentration of baryons in the center.

preprint2013arXiv

Kinematics of massive star ejecta in the Milky Way as traced by 26Al

Massive stars form in groups and their winds and supernova explosions create superbubbles up to kpc in size. Their ejecta are important for the dynamics of the interstellar medium and chemical evolution models. However, ejecta kinematics and the characteristic scales in space and time are rather unexplored beyond ~10 ka. Through measurement of radioactive 26Al with its decay time constant of ~10^6 years, we aim to trace the kinematics of massive-star and supernova ejecta over million-year time scales. Our goal is to identify the mixing time scale and the spatio-kinematics of such ejecta from the pc to kpc scale. We use the SPI spectrometer on INTEGRAL and its observations along the Galactic ridge to trace the line shape of the 1808.63 keV gamma-ray line from 26Al decay. We determine line centroids and compare these to Doppler shift expectations from large-scale systematic rotation around the Galaxy's center, as observed in other Galactic objects. We measure the radial velocities of gas traced by 26Al, averaged over the line of sight, as a function of Galactic longitude. We find substantially higher velocities than expected from Galactic rotation, the average bulk velocity being ~200 km s^-1 larger than the Galactic-rotation prediction. The observed radial velocity spread implies a Doppler-broadening of the gamma-ray line that is consistent with our measurements of the overall line width. We can reproduce the observed characteristics with 26Al sources located along the inner spiral arms, when we add a global blow-out into the forward direction away from arms, such as expected from massive stars offset towards the spiral-arm leading edge. The superbubbles-halo connection implies angular-momentum transfer in the disk-halo system and also radial gas flows.

preprint2013arXiv

Mapping the Three-Dimensional Density of the Galactic Bulge with VVV Red Clump Stars

The inner Milky Way is dominated by a boxy, triaxial bulge which is believed to have formed through disk instability processes. Despite its proximity, its large-scale properties are still not very well known, due to our position in the obscuring Galactic disk. Here we make a measurement of the three-dimensional density distribution of the Galactic bulge using red clump giants identified in DR1 of the VVV survey. Our density map covers the inner (2.2x1.4x1.1)kpc of the bulge/bar. Line-of-sight density distributions are estimated by deconvolving extinction and completeness corrected K-band magnitude distributions. In constructing our measurement, we assume that the three-dimensional bulge is 8-fold mirror triaxially symmetric. In doing so we measure the angle of the bar-bulge to the line-of-sight to be (27+- 2)deg, where the dominant error is systematic arising from the details of the deconvolution process. The resulting density distribution shows a highly elongated bar with projected axis ratios ~(1:2.1) for isophotes reaching ~2kpc along the major axis. Along the bar axes the density falls off roughly exponentially, with axis ratios (10:6.3:2.6) and exponential scale-lengths (0.70:0.44:0.18)kpc. From about 400pc above the Galactic plane, the bulge density distribution displays a prominent X-structure. Overall, the density distribution of the Galactic bulge is characteristic for a strongly boxy/peanut shaped bulge within a barred galaxy.

preprint2013arXiv

Metallicity gradients through disk Instability: A simple model for the Milky Way's boxy bulge

Observations show a clear vertical metallicity gradient in the Galactic bulge, which is often taken as a signature of dissipative processes in the formation of a classical bulge. Various evidence shows, however, that the Milky Way is a barred galaxy with a boxy bulge representing the inner three-dimensional part of the bar. Here we show with a secular evolution N-body model that a boxy bulge formed through bar and buckling instabilities can show vertical metallicity gradients similar to the observed gradient, if the initial axisymmetric disk had a comparable radial metallicity gradient. In this framework the range of metallicities in bulge fields constrains the chemical structure of the Galactic disk at early times, before bar formation. Our secular evolution model was previously shown to reproduce inner Galaxy star counts and we show here that it also has cylindrical rotation. We use it to predict a full mean metallicity map across the Galactic bulge from a simple metallicity model for the initial disk. This map shows a general outward gradient on the sky as well as longitudinal perspective asymmetries. We also briefly comment on interpreting metallicity gradient observations in external boxy bulges.

preprint2013arXiv

The mass and angular momentum distribution of simulated massive early-type galaxies to large radii

We study the dark and luminous mass distributions, circular velocity curves (CVC), line-of-sight kinematics, and angular momenta for a sample of 42 cosmological zoom simulations of massive galaxies. Using a temporal smoothing technique, we are able to reach large radii. We find that: (i)The dark matter halo density profiles outside a few kpc follow simple power-law models, with flat dark matter CVCs for lower-mass systems, and rising CVCs for high-mass haloes. The projected stellar density distributions at large radii can be fitted by Sersic functions with n>10, larger than for typical ETGs. (ii)The massive systems have nearly flat total CVCs at large radii, while the less massive systems have mildly decreasing CVCs. The slope of the CVC at large radii correlates with v_circ itself. (iii)The dark matter fractions within Re are in the range 15-30% and increase to 40-65% at 5Re. Larger and more massive galaxies have higher dark matter fractions. (iv)The short axes of simulated galaxies and their host dark matter haloes are well aligned and their short-to-long axis ratios are correlated. (v)The stellar vrms(R) profiles are slowly declining, in agreement with planetary nebulae observations in the outer haloes of most ETGs. (vi)The line-of-sight velocity fields v show that rotation properties at small and large radii are correlated. Most radial profiles for the cumulative specific angular momentum parameter lambda(R) are nearly flat or slightly rising, with values in [0.06,0.75] from 2Re to 5Re. (vii)Stellar mass, ellipticity at 5Re, and lambda(5Re) are correlated: the more massive systems have less angular momentum and are rounder, as for observed ETGs. (viii)More massive galaxies with a large fraction of accreted stars have radially anisotropic velocity distributions outside Re. Tangential anisotropy is seen only for galaxies with high fraction of in-situ stars. (Full abstract in PDF)

preprint2013arXiv

The Nuclear Cluster of the Milky Way: Total Mass and Luminosity

Here we present the fundamental properties of the nuclear cluster of the Milky Way. First, we derive its structural properties by constructing a density map of the central 1000'' using extinction-corrected star counts. We can describe the data with a two-component model built from Sersic profiles. The inner nearly spherical component is the nuclear cluster. The outer, strongly flattened component can be identified with the stellar component of the circumnuclear zone. Second, we enlarge the radius inside which detailed dynamics are available from $1\,$pc to $4\,$pc. We use more than 10000 individual proper motions and more than 2700 radial velocities. We determine the cluster mass by means of isotropic spherical Jeans modeling. We get a nuclear cluster mass within 100'' of $M_{100''}=(6.11 \pm 0.52|_{\mathrm{fix} R_0}\pm 0.97|_{R_0}) \times 10^6 $ M$_{\odot}$, which corresponds to a total cluster mass of M$_{\mathrm{NC}}=(13.08 \pm 2.51|_{\mathrm{fix} R_0}\pm 2.08|_{R_0}) \times 10^6 $ M$_{\odot}$. By combination of our mass with the flux we calculate $M/L=0.50 \pm 0.12 M_{\odot}/L_{\odot,\mathrm{Ks}}$ for the central 100''. That is broadly consistent with a Chabrier IMF. With its mass and a luminosity of M$_{\mathrm{Ks}}=-15.30\pm0.26$ the nuclear cluster is a bright and massive specimen with a typical size.

preprint2012arXiv

A tale of two tails and an off-centered envelope: diffuse light around the cD galaxy NGC 3311 in the Hydra I cluster

The formation of intracluster light and of the extended halos around brightest cluster galaxies is closely related to morphological transformation, tidal stripping, and disruption of galaxies in clusters. We analyze Ks- and V-band surface photometry as well as deep long-slit spectra, and establish a link between the structures in the light distribution, the absorption line kinematics, and the LOS velocity distributions of nearby galaxies and planetary nebulae (PNs). The central galaxy NGC 3311 is surrounded by an extended symmetric outer halo with n=10 and an additional, off-centered envelope ~ 50" to the North-East. Its luminosity L_V= 1.2x10^{10} +/- 6.0 x 10^8 L_sun corresponds to ~50 % of the luminosity of the symmetric halo in the same region. Based on measured PN velocities, at least part of the off-centered envelope consists of high-velocity accreted stars. We have also discovered two tidal streams in the cluster center, emerging from the dwarf galaxy HCC 026 and from the S0 galaxy HCC 007. The HCC 026 stream is redshifted by ~1200 km/s with respect to NGC 3311, similarly as HCC 026 itself, a fraction of PNs in the off-centered envelope, and several other dwarf galaxies nearby. The stars in one of the HCC 026 tails are known to be consistent with the low-metallicity population of HCC 026, and our photometry shows that this galaxy is already mostly dissolved in the tidal field. The tidal stream around HCC 007 extends over at least 110 kpc. It is fairly thick and is brighter on the side of the asymmetric outer halo of NGC 3311, which it may join. Its luminosity is several 10^9 L_sun, similar to the luminosity of the stripped-down galaxy HCC 007. The redshift of the stream is determined from a few PN velocities and is similar to that for HCC 007 and HCC 026.

preprint2012arXiv

NMAGIC made-to-measure particle models of galaxies

The parallel code NMAGIC is an implementation of a particle-based method to create made-to-measure models in agreement with observations of galaxies. It works by slowly correcting the particle weights of an evolving N-body system, until a satisfactory compromise is achieved between the goodness of the fit to a given set of observational data, and some degree of smoothness (regularization) of the underlying particle model. We briefly describe the method together with a new regularization scheme in phase-space, which improves recovering the correct orbit structure in the models. We also mention some practical applications showing the power of the technique in investigating the dynamics of galaxies.

preprint2012arXiv

Regularizing made-to-measure particle models of galaxies

Made-to-measure methods such as the parallel code NMAGIC are powerful tools to build galaxy models reproducing observational data. They work by adapting the particle weights in an N-body system until the target observables are well matched. Here we introduce a moving prior regularization (MPR) method for such particle models. It is based on determining from the particles a distribution of priors in phase-space, which are updated in parallel with the weight adaptation. This method allows one to construct smooth models from noisy data without erasing global phase-space gradients. We first apply MPR to a spherical system for which the distribution function can in theory be uniquely recovered from idealized data. We show that NMAGIC with MPR indeed converges to the true solution with very good accuracy, independent of the initial particle model. Compared to the standard weight entropy regularization, biases in the anisotropy structure are removed and local fluctuations in the intrinsic distribution function are reduced. We then investigate how the uncertainties in the inferred dynamical structure increase with less complete and noisier kinematic data, and how the dependence on the initial particle model also increases. Finally, we apply the MPR technique to the two intermediate-luminosity elliptical galaxies NGC 4697 and NGC 3379, obtaining smoother dynamical models in luminous and dark matter potentials.

preprint2012arXiv

Rotation of classical bulges during secular evolution of barred galaxies

Bar driven secular evolution plays a key role in changing the morphology and kinematics of disk galaxies, leading to the formation of rapidly rotating boxy/peanut bulges. If these disk galaxies also hosted a preexisting classical bulge, how would the secular evolution influence the classical bulge, and also the observational properties. We first study the co-evolution of a bar and a preexisting non-rotating low-mass classical bulge such as might be present in galaxies like the Milky Way. It is shown with N-body simulations that during the secular evolution, such a bulge can gain significant angular momentum emitted by the bar through resonant and stochastic orbits. Thereby it transforms into a cylindrically rotating, anisotropic and triaxial object, embedded in the fast rotating boxy bulge that forms via disk instability (Saha et al. 2012). The composite boxy/peanut bulge also rotates cylindrically. We then show that the growth of the bar depends only slightly on the rotation properties of the preexisting classical bulge. For the initially rotating small classical bulge, cylindrical rotation in the resulting composite boxy/peanut bulge extends to lower heights (Saha & Gerhard 2012). More massive classical bulges also gain angular momentum emitted by the bar, inducing surprisingly large rotational support within about 4 Gyrs (Saha et al. in prep).

preprint2012arXiv

The Apache Point Observatory Galactic Evolution Experiment: First Detection of High Velocity Milky Way Bar Stars

Commissioning observations with the Apache Point Observatory Galactic Evolution Experiment (APOGEE), part of the Sloan Digital Sky Survey III, have produced radial velocities (RVs) for ~4700 K/M-giant stars in the Milky Way bulge. These high-resolution (R \sim 22,500), high-S/N (>100 per resolution element), near-infrared (1.51-1.70 um; NIR) spectra provide accurate RVs (epsilon_v~0.2 km/s) for the sample of stars in 18 Galactic bulge fields spanning -1<l<20 deg, |b|<20 deg, and dec>-32 deg. This represents the largest NIR high-resolution spectroscopic sample of giant stars ever assembled in this region of the Galaxy. A cold (sigma_v~30 km/s), high-velocity peak (V_GSR \sim +200 km/s) is found to comprise a significant fraction (~10%) of stars in many of these fields. These high RVs have not been detected in previous MW surveys and are not expected for a simple, circularly rotating disk. Preliminary distance estimates rule out an origin from the background Sagittarius tidal stream or a new stream in the MW disk. Comparison to various Galactic models suggests that these high RVs are best explained by stars in orbits of the Galactic bar potential, although some observational features remain unexplained.

preprint2011arXiv

Dynamical evolution of a bulge in an N-body model of the Milky Way

The detailed dynamical structure of the bulge in the Milky Way is currently under debate. Although kinematics of the bulge stars can be well reproduced by a boxy-bulge, the possible existence of a small embedded classical bulge can not be ruled out. We study the dynamical evolution of a small classical bulge in a model of the Milky Way using a self-consistent high resolution N-body simulation. Detailed kinematics and dynamical properties of such a bulge are presented.

preprint2011arXiv

Spin-up of low mass classical bulges in barred galaxies

Secular evolution is one of the key routes through which galaxies evolve along the Hubble sequence. Not only the disk undergoes morphological and kinematic changes, but also a preexisting classical bulge may be dynamically changed by the secular processes driven primarily by the bar. We study the influence of a growing bar on the dynamical evolution of a low mass classical bulge such as might be present in galaxies like the Milky Way. Using self-consistent high resolution {\it N}-body simulations, we study how an initially isotropic non-rotating small classical bulge absorbs angular momentum emitted by the bar. The basic mechanism of this angular momentum exchange is through resonances and a considerable fraction of the angular momentum is channeled through Lagrange point (-1:1) and ILR (2:1) orbits. In the phase of rapid dynamical growth, also retrograde non-resonant orbits absorb significant angular momentum. As a result of this angular momentum gain, the initially non-rotating classical bulge transforms into a fast rotating, radially anisotropic and triaxial object, embedded in the similarly fast rotating boxy bulge formed from the disk. Towards the end of the evolution, the classical bulge develops cylindrical rotation. By that time, its inner regions host a "classical bulge-bar" whose distinct kinematics could serve as direct observational evidence for the secular evolution in the galaxy. Some implications of these results are discussed briefly.

preprint2011arXiv

The inner Galactic bulge: evidence for a nuclear bar?

Recent data from the VVV survey have strengthened evidence for a structural change in the Galactic bulge inwards of |l|<=4 deg. Here we show with an N-body barred galaxy simulation that a boxy bulge formed through the bar and buckling instabilities effortlessly matches measured bulge longitude profiles for red clump stars. The same simulation snapshot was earlier used to clarify the apparent boxy bulge - long bar dichotomy, for the same orientation and scaling. The change in the slope of the model longitude profiles in the inner few degrees is caused by a transition from highly elongated to more nearly axisymmetric isodensity contours in the inner boxy bulge. This transition is confined to a few degrees from the Galactic plane, thus the change of slope is predicted to disappear at higher Galactic latitudes. We also show that the nuclear star count map derived from this simulation snapshot displays a longitudinal asymmetry similar to that observed in the 2MASS data, but is less flattened to the Galactic plane than the 2MASS map. These results support the interpretation that the Galactic bulge originated from disk evolution, and question the evidence advanced from star count data for the existence of a secondary nuclear bar in the Milky Way.

preprint2011arXiv

Unifying boxy bulge and planar long bar in the Milky Way

It has been known for some time that the Milky Way is a barred disk galaxy. More recently several studies inferred from starcount observations that the Galaxy must contain a separate, new, flat long bar component, twisted relative to the barred bulge. Here we use a simulation with a boxy bulge and bar to suggest that these observations can be reproduced with a single structure. In this simulation a stellar bar evolved from the disk, and the boxy bulge originated from it through secular evolution and the buckling instability. We calculate starcount distributions for this model at different longitudes and latitudes, in a similar way as observers have done for resolved starcounts. Good agreement between the simulation and the observations can be achieved for a suitable snapshot, even though the simulation has a single boxy bulge and bar} structure. In this model, part of the long bar signature is due to a volume effect in the starcounts, and another part is due to chosing a snapshot in which the planar part of the boxy bulge and bar has developed leading ends through interaction with the adjacent spiral arm heads. We also provide predictions from this model for the line-of-sight velocity distributions at the longitudes with the long bar signature, for comparison with upcoming surveys.

preprint2011arXiv

Using NMAGIC to probe the dark matter halo and orbital structure of the X-ray bright, massive elliptical galaxy, NGC 4649

We create dynamical models of the massive elliptical galaxy, NGC 4649, using the N-body made-to-measure code, NMAGIC, and kinematic constraints from long-slit and planetary nebula (PN) data. We explore a range of potentials based on previous determinations from X-ray observations and a dynamical model fitting globular cluster (GC) velocities and a stellar density profile. The X-ray mass distributions are similar in the central region but have varying outer slopes, while the GC mass profile is higher in the central region and on the upper end of the range further out. Our models cannot differentiate between the potentials in the central region, and therefore if non-thermal pressures or multi-phase components are present in the hot gas, they must be smaller than previously inferred. In the halo, we find that the PN velocities are sensitive tracers of the mass, preferring a less massive halo than that derived from the GC mass profile, but similar to one of the mass distributions derived from X-rays. Our results show that the GCs may form a dynamically distinct system, and that the properties of the hot gas derived from X-rays in the outer halo have considerable uncertainties that need to be better understood. Estimating the mass in stars using photometric information and a stellar population mass-to-light ratio, we infer a dark matter mass fraction in NGC 4649 of ~0.39 at 1Re (10.5 kpc) and ~0.78 at 4Re. We find that the stellar orbits are isotropic to mildly radial in the central ~6 kpc depending on the potential assumed. Further out, the orbital structure becomes slightly more radial along R and more isotropic along z, regardless of the potential assumed. In the equatorial plane, azimuthal velocity dispersions dominate over meridional velocity dispersions, implying that meridional velocity anisotropy is the mechanism for flattening the stellar system.

preprint2010arXiv

Counter-dispersed slitless-spectroscopy technique: planetary nebula velocities in the halo of NGC 1399

Using a counter-dispersed slitless spectroscopy technique, we detect and measure the line-of-sight velocities of 187 planetary nebulae (PNe) around one of the nearest cD galaxies, NGC 1399, with FORS1 on the VLT. We describe the method for identifying and classifying the emission-line sources and the procedure for computing their J2000 coordinates and velocities. The number of PN detections and the errors in the velocity measurements (37 km/s indicate that this technique is comparable to other methods, such as that described by Teodorescu et al. (2005). We present the spatial distribution of the PNe and a basic analysis of their velocities. The PN two-dimensional velocity field shows marginal rotation consistent with other studies. We also find a low-velocity substructure in the halo and a flatter velocity-dispersion profile compared to previous observations that extends to ~400 arcsec. The detection of a low-velocity subcomponent underscores the importance of discrete velocity tracers for the detection of un-mixed components. The new velocity-dispersion profile is in good agreement with revised velocity dispersions for the red globular clusters in NGC 1399, using the data of Schuberth et al. (2009). The outer parts of this profile are consistent with one of the dynamical models of Kronawitter et al. (2000), which corresponds to a circular velocity of ~340 km/s and a rescaled B-band mass-to-light ratio of ~20 at 7' radius. These measurements trace the kinematics of the outer halo and disentangle the heterogenous populations in the Fornax Cluster core. The new data set the stage for a revised dynamical model of the outer halo of NGC 1399.

preprint2010arXiv

Diffuse Light in Galaxy Clusters

Diffuse intracluster light (ICL) has now been observed in nearby and in intermediate redshift clusters. Individual intracluster stars have been detected in the Virgo and Coma clusters and the first color-magnitude diagram and velocity measurements have been obtained. Recent studies show that the ICL contains of the order of 10% and perhaps up to 30% of the stellar mass in the cluster, but in the cores of some dense and rich clusters like Coma, the local ICL fraction can be high as 40%-50%. What can we learn from the ICL about the formation of galaxy clusters and the evolution of cluster galaxies? How and when did the ICL form? What is the connection to the central brightest cluster galaxy? Cosmological N-body and hydrodynamical simulations are beginning to make predictions for the kinematics and origin of the ICL. The ICL traces the evolution of baryonic substructures in dense environments and can thus be used to constrain some aspects of cosmological simulations that are most uncertain, such as the modeling of star formation and the mass distribution of the baryonic component in galaxies.

preprint2010arXiv

Distinct core and halo stellar populations and the formation history of the bright Coma cluster early-type galaxy NGC 4889

We study the stellar population far into the halo of one of the two brightest galaxies in the Coma cluster, NGC 4889, based on deep medium resolution spectroscopy with FOCAS at the Subaru 8.2m telescope. We fit single stellar population models to the measured line-strength (Lick) indices (Hbeta, Mgb, [MgFe]' and <Fe>). Combining with literature data, we construct radial profiles of metallicity, [alpha/Fe] element abundance ratio and age for NGC 4889, from the center out to ~60 kpc (~4Re). We find evidence for different chemical and star formation histories for stars inside and outside 1.2Re = 18 kpc radius. The inner regions are characterized by a steep [Z/H] gradient and high [alpha/Fe] at ~2.5 times solar value. In the halo, between 18 and 60 kpc, the [Z/H] is near-solar with a shallow gradient, while [alpha/Fe] shows a strong negative gradient, reaching solar values at ~60 kpc. We interpret these data in terms of different formation histories for both components. The data for the inner galaxy are consistent with a rapid, quasi-monolithic, dissipative merger origin at early redshifts, followed by one or at most a few dry mergers. Those for the halo argue for later accretion of stars from old systems with more extended star formation histories. The half-light radius of the inner component alone is estimated as ~6 kpc, suggesting a significantly smaller size of this galaxy in the past. This may be the local stellar population signature of the size evolution found for early-type galaxies from high-redshift observations.

preprint2010arXiv

Pattern speeds in the Milky Way

A brief review is given of different methods used to determine the pattern speeds of the Galactic bar and spiral arms. The Galactic bar rotates rapidly, with corotation about halfway between the Galactic center and the Sun, and outer Lindblad resonance not far from the solar orbit, R0. The Galactic spiral arms currently rotate with a distinctly slower pattern speed, such that corotation is just outside R0. Both structures therefore seem dynamically decoupled.

preprint2010arXiv

Steepening mass profiles, dark matter and environment of X-ray bright elliptical galaxies

We use a new non-parametric Bayesian approach to obtain the most probable mass distributions and circular velocity curves along with their confidence ranges, given deprojected density and temperature profiles of the hot gas surrounding X-ray bright elliptical galaxies. For a sample of six X-ray bright ellipticals, we find that all circular velocity curves are rising in the outer parts due to a combination of a rising temperature profile and a logarithmic pressure gradient that increases in magnitude. Comparing the circular velocity curves we obtain from X-rays to those obtained from dynamical models, we find that the former are often lower in the central ~10 kpc. This is probably due to a combination of: i) Non-thermal contributions of up to ~35% in the pressure (with stronger effects in NGC 4486), ii) multiple-temperature components in the hot gas, iii) incomplete kinematic spatial coverage in the dynamical models, and iv) mass profiles that are insufficiently general in the dynamical modelling. Complementing the total mass information from the X-rays with photometry and stellar population models to infer the dark matter content, we find evidence for massive dark matter haloes with dark matter mass fractions of ~35-80% at 2Re, rising to a maximum of 80-90% at the outermost radii. We also find that the six galaxies follow a Tully-Fisher relation with slope ~4 and that their circular velocities at 1Re correlate strongly with the velocity dispersion of the local environment. As a result, the galaxy luminosity at 1Re also correlates with the velocity dispersion of the environment. These relations suggest a close link between the properties of central X-ray bright elliptical galaxies and their environments (abridged).

preprint2010arXiv

The outer halos of elliptical galaxies

Recent progress is summarized on the determination of the density distributions of stars and dark matter, stellar kinematics, and stellar population properties, in the extended, low surface brightness halo regions of elliptical galaxies. With integral field absorption spectroscopy and with planetary nebulae as tracers, velocity dispersion and rotation profiles have been followed to ~4 and ~5-8 effective radii, respectively, and in M87 to the outer edge at ~150 kpc. The results are generally consistent with the known dichotomy of elliptical galaxy types, but some galaxies show more complex rotation profiles in their halos and there is a higher incidence of misalignments, indicating triaxiality. Dynamical models have shown a range of slopes for the total mass profiles, and that the inner dark matter densities in ellipticals are higher than in spiral galaxies, indicating earlier assembly redshifts. Analysis of the hot X-ray emitting gas in X-ray bright ellipticals and comparison with dynamical mass determinations indicates that non-thermal components to the pressure may be important in the inner ~10 kpc, and that the properties of these systems are closely related to their group environments. First results on the outer halo stellar population properties do not yet give a clear picture. In the halo of one bright galaxy, lower [alpha/Fe] abundances indicate longer star formation histories pointing towards late accretion of the halo. This is consistent with independent evidence for on-going accretion, and suggests a connection to the observed size evolution of elliptical galaxies with redshift.

preprint2009arXiv

Intracluster Light in the Virgo Cluster: Large Scale Distribution

The intracluster light (ICL) is a faint diffuse stellar component in clusters made of stars not bound to individual galaxies. We have carried out a large scale study of this component in the nearby Virgo cluster. The diffuse light is traced using planetary nebulae (PNe). The PNe are detected in the on-band image due to their strong emission in the [OIII] 5007 line, but disappear in the off-band image. The contribution of Ly-alpha emitters at z=3.14 are corrected statistically using blank field surveys. We have surveyed a total area of 3.3 square degrees in the Virgo cluster with eleven fields located at different radial distances. Those fields located at smaller radii than 80 arcmin from the cluster center contain most of the detected diffuse light. In this central region of the cluster, the ICL has a surface brightness in the range 28.8 - 30 mag per sqarsec in the B band, it is not uniformly distributed, and represents about 7% of the total galaxy light in this area. At distances larger than 80 arcmin the ICL is confined to single fields and individual sub-structures, e.g. in the Virgo sub-clump B, the M60/M59 group. For several fields at 2 and 3 degrees from the Virgo cluster center we set only upper limits. These results indicate that the ICL is not homogeneously distributed in the Virgo core, and it is concentrated in the high density regions of the Virgo cluster, e.g. the cluster core and other sub-structures. Outside these regions, the ICL is confined within areas of 100 kpc in size, where tidal effects may be at work. These observational results link the formation of the ICL with the formation history of the most luminous cluster galaxies.

preprint2009arXiv

The edge of the M87 halo and the kinematics of the diffuse light in the Virgo cluster core

We present high resolution FLAMES/VLT spectroscopy of intracluster planetary nebula (ICPN) candidates, targeting three new fields in the Virgo cluster core with surface brightness down to mu_B = 28.5. Based on the projected phase space information we separate the old and 12 newly-confirmed PNs into galaxy and intracluster components. The M87 PNs are confined to the extended stellar envelope of M87, within a projected radius of ~ 160 kpc, while the ICPNs are scattered across the whole surveyed region between M87 and M86. The velocity dispersions determined from the M87 PNs at projected radii of 60 kpc and 144 kpc show that the galaxy's velocity dispersion profile decreases in the outer halo, down to 78 +/- 25 km/s. A Jeans model for the M87 halo stars in the gravitational potential traced by the X-ray emission fits the observed velocity dispersion profile only if the stellar orbits are strongly radially anisotropic (beta ~= 0.4 at r ~= 10 kpc increasing to 0.8 at the outer edge), and if additionally the stellar halo is truncated at ~= 150 kpc average elliptical radius. From the spatial and velocity distribution of the ICPNs we infer that M87 and M86 are falling towards each other and that we may be observing them just before the first close pass. The inferred luminosity-specific PN numbers for the M87 halo and the ICL are in the range of values observed for old (> 10 Gyr) stellar populations (abridged).

preprint1997arXiv

Breaking the degeneracy between anisotropy and mass: The dark halo of the E0 galaxy NGC 6703

(abridged) We have measured line-of-sight velocity profiles (VPs) in the E0 galaxy NGC 6703 out to 2.6 R_e. From these data we constrain the mass distribution and the anisotropy of the stellar orbits in this galaxy. We have developed a non-parametric technique to determine the DF f(E,L^2) directly from the kinematic data. From Monte Carlo tests using the spatial extent, sampling, and error bars of the NGC 6703 data we find that smooth underlying DFs can be recovered to an rms accuracy of 12%, and the anisotropy parameter beta(r) to an accuracy of 0.1, in a given potential. An asymptotically constant halo circular velocity v_0 can be determined with an accuracy of +- \lta 50km/s. For NGC 6703 we determine the true circular velocity at 2.6 R_e to be 250 +- 40km/s at 95% c.l., corresponding to a total mass in NGC 6703 inside 78'' (13.5 h_50^-1 kpc), of 1.6-2.6 x 10^11 h_50^-1 Msun. No model without dark matter will fit the data; however, a maximum stellar mass model in which the luminous component provides nearly all the mass in the centre does. In such a model, the total luminous mass inside 78'' is 9 x 10^10 Msun and the integrated M/L_B=5.3-10, corresponding to a rise from the center by at least a factor of 1.6. The anisotropy of the stellar distribution function in NGC 6703 changes from near-isotropic at the centre to slightly radially anisotropic (beta=0.3-0.4 at 30'', beta=0.2-0.4 at 60'') and is not well-constrained at the outer edge of the data. Our results suggest that also elliptical galaxies begin to be dominated by dark matter at radii of \sim 10kpc.

preprint1996arXiv

The photometric structure of the inner Galaxy

The light distribution in the inner few kiloparsecs of the Milky Way is recovered non-parametrically from a dust-corrected near-infrared COBE/DIRBE surface brightness map of the inner Galaxy. The best fits to the photometry are obtained when the Sun is assumed to lie $\sim14\pm4\pc$ below the plane. The recovered density distributions clearly show an elongated three-dimensional bulge set in a highly non-axisymmetric disk. In the favoured models, the bulge has axis ratios $1{:}0.6{:}0.4$ and semi-major axis length $\sim2\kpc$. Its nearer long axis lies in the first quadrant. The bulge is surrounded by an elliptical disk that extends to $\sim2\kpc$ on the minor axis and $\sim3.5\kpc$ on the major axis. In all models there is a local density minimum $\sim2.2\kpc$ down the minor axis. The subsequent maximum $\sim3\kpc$ down the minor axis (corresponding to $l\simeq-22°$ and $l\simeq 17°$) may be associated with the Lagrange point L$_4$. From this identification and the length of the bulge-bar, we infer a pattern speed $Ω_b\simeq 60-70\kms\kpc^{-1}$ for the bar. Experiments in which pseudo-data derived from models with spiral structure were deprojected under the assumption that the Galaxy is either eight-fold or four-fold symmetric, indicate that the highly non-axisymmetric disks recovered from the COBE data could reflect spiral structure within the Milky Way if that structure involves density contrasts greater than $\gta 3$ at NIR wavelengths. These experiments indicate that the angle $ϕ_0$ between the Sun--centre line and a major axis of the bulge lies near $20°$.

preprint1996arXiv

Two modes of gas flow in a single barred galaxy

We investigate stationary gas flows in a fixed, rotating barred potential. The gas is assumed to be isothermal with an effective sound speed c_s, and the equations of motion are solved with smoothed particle hydrodynamics (SPH). Since the thermal energy in cloud random motions is negligible compared to the orbital kinetic energy, no dependence of the flow on c_s is expected. However, this is not the case when shocks are involved. For low values of c_s an open, off-axis shock flow forms that is characteristic for potentials with an inner Lindblad resonance (ILR). Through this shock the gas streams inwards from x_1 to x_2-orbits. At high sound speeds the gas arranges itself in a different, on-axis shock flow pattern. In this case, there is no gas on x_2-orbits, demonstrating that the gas can behave as if there were no ILR. The critical effective sound speed dividing the two regimes is in the range of values observed in the Milky Way. We give a heuristic explanation for this effect. A possible consequence is that star formation may change the structure of the flow by which it was initiated. Low-mass galaxies should predominantly be in the on-axis regime. A brief comparison of our SPH results with those from a grid-based hydrodynamic code is also given.

preprint1995arXiv

On the deprojection of axisymmetric bodies

Axisymmetric density distributions are constructed which are invisible when viewed from a range of inclination angles $i$. By adding such distributions to a model galaxy, it can be made either disky or boxy without in any way affecting its projected image. As the inclination of a galaxy decreases from edge-on to face on, the range of `invisible' densities, the uncertainty in the deprojection, and the sensitivity of the deprojection to noise all increase. The relation between these phenomena is clarified by an analysis of Palmer's deprojection algorithm. These results imply that disk-to-bulge ratios are in principle ill-determined from photometry unless the disk is strong or the system is seen precisely edge-on. The uncertain role of third integrals in galaxies makes it unclear to what degree this indeterminacy can be resolved by kinematic studies.

preprint1995arXiv

On the deprojection of the Galactic bulge

An algorithm is developed and tested for the problem posed by photometric observations of the bulge of the Milky Way. The latter subtends a non-trivial solid angle on the sky, and we show that this permits inversion of the projected brightness distribution under the assumption that the bulge has three orthogonal mirror planes of specified orientation. A serious error in the assumed orientation of the mirror planes should be detectable.