Source author record

Glenn van de Ven

Glenn van de Ven 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

89works
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

89 published item(s)

preprint2026arXiv

Distance measurements from the internal dynamics of globular clusters: Application to the Sombrero galaxy (M104)

Globular clusters (GCs) are dense star clusters found in all massive galaxies. Recent work has established that they follow a tight relation between their internal stellar velocity dispersion $σ$ and luminosity, enabling accurate distance measurements. In this work, we aim to apply this GC velocity dispersion (GCVD) distance method to measure the distance to M104 (NGC 4594, the Sombrero galaxy). We have measured internal stellar velocity dispersions for 85 globular clusters (GCs) and one ultra-compact dwarf galaxy around M104 using high-resolution multi-object integrated-light spectroscopy with FLAMES/GIRAFFE on the Very Large Telescope. The measured velocity dispersions range from $σ= 4 - 30$ km s$^{-1}$, with a mean uncertainty of $Δσ= 2.5$ km s$^{-1}$. For a subset of 77 GCs with $V$-band magnitudes and reliable velocity dispersion measurements above $σ> 4$ km s$^{-1}$, we constructed the $M_V$-$σ$ relation to measure the distance to M104, finding $D=9.00\pm0.29$ (stat.)~$\pm0.26$ (sys.) Mpc. The GCs follow the Milky Way and M31 $M_V-σ$ relation closely, with the exception of the luminous ultra-compact dwarf SUCD1, which is nearly one magnitude brighter than the mean relation. 29 GCs in the sample have sizes determined from Hubble Space Telescope imaging which allowed us to determine their masses and $V$-band dynamical mass-to-light ratios (M/L$_V$). We find a mean $<M/L_V> = 2.6 M_{\odot}/L_{\odot}$ for the luminous ($M_V < -8$ mag) M 104 GCs, which is higher than the Milky Way GCs, but is reminiscent of the brightest GCs in Centaurus A. With the exception of SUCD1, the GCs of M104 follow the GCVD relation irrespective of their mass-to-light ratio.

preprint2026arXiv

The MAGPI Survey: co-evolution of baryons and dark matter in star-forming disk-like galaxies at $0.1 \lesssim z \lesssim 0.85$

We present a comprehensive analysis of the dark matter (DM) content and its structural dependence in star-forming disk-like galaxies at intermediate redshifts ($0.1 \lesssim z \lesssim 0.85$), utilizing spatially resolved kinematic data from the MAGPI survey. We report the following: (1) Low stellar mass galaxies ($M_{\rm star} < 10^{9.5}\, M_\odot$) are strongly DM dominated across all radii, with average $\langle f_{_{\rm DM}} \rangle \sim 0.85$, while high-mass ($M_{\rm star} > 10^{10.5}\, M_\odot$) systems exhibit relatively low DM fractions in their inner regions ($\langle f_{_{\rm DM}} \rangle \sim 0.47$) which is equivalent to local massive disk galaxies (e.g., Milky Way and Andromeda). This suggests a mass-dependent structural dichotomy, most-likely governed by a combination of internal galactic processes and environmental influences. (2) A tight inverse correlation between $f_{_{\rm DM}}$ and baryon mass surface density ($Σ_{\rm bar}$), with intrinsic scatter of $\sim 0.11$ dex. This is consistent with an inside-out baryon assembly scenario and suggests that the fundamental structural correlations of galaxies were already established by $z\sim 0.85$. (3) No significant evolution in $f_{_{\rm DM}}$ with redshift across the MAGPI window, and when combined with higher-redshift ($0.6 \leq z \leq 1.5$) data from Sharma et al. 2025, we quantitatively show that the reported decline in $f_{_{\rm DM}}(z)$ is most-likely due to observational biases against low-mass systems at $z > 1$. These results offer empirical evidence for a scenario in which disk-like galaxies evolve through a co-regulated build-up of baryonic and DM components, preserving internal structural regularities (such as the total mass distribution and rotation-curve shape) throughout cosmic time.

preprint2022arXiv

Gas inflows in the polar ring of NGC 4111: the birth of an AGN

We have used Hubble Space Telescope (HST) images, SAURON Integral Field Spectroscopy (IFS) and adaptative optics assisted Gemini NIFS near-infrared K-band IFS to map the stellar and gas distribution, excitation and kinematics of the inner few kpc of the nearby edge-on S0 galaxy NGC 4111. The HST images map its $\approx$ 450 pc diameter dusty polar ring, with an estimated gas mass $\ge10^7$ M$_\odot$. The NIFS datacube maps the inner 110 pc radius at $\approx$ 7 pc spatial resolution revealing a $\approx$ 220 pc diameter polar ring in hot ($2267\pm166$ K) molecular H$_2$ 1-0 S(1) gas embedded in the polar ring. The stellar velocity field shows disk-dominated kinematics along the galaxy plane both in the SAURON large-scale and in the NIFS nuclear-scale data. The large-scale [O III] $\lambda5007$ Åvelocity field shows a superposition of two disk kinematics: one similar to that of the stars and another along the polar ring, showing non-circular motions that seem to connect with the velocity field of the nuclear H$_2$ ring, whose kinematics indicate accelerated inflow to the nucleus. The estimated mass inflow rate is enough not only to feed an Active Galactic Nucleus (AGN) but also to trigger circumnuclear star formation in the near future. We propose a scenario in which gas from the polar ring, which probably originated from the capture of a dwarf galaxy, is moving inwards and triggering an AGN, as supported by the local X-ray emission, which seems to be the source of the H$_2$ 1-0 S(1) excitation. The fact that we see neither near-UV nor Br$γ$ emission suggests that the nascent AGN is still deeply buried under the optically thick dust of the polar ring.

preprint2022arXiv

Mass of the dynamically hot inner stellar halo predicts the ancient accreted stellar mass

Galactic dynamical structures are fossil records of the assembly histories of galaxies. By analyzing the cosmological hydrodynamical simulation TNG50, we find that a dynamical structure that we call the "hot inner stellar halo," defined by stars on dynamically hot orbits with circularity $λ_z < 0.5$ at $3.5\,{\rm kpc}<r \lesssim 2\,R_e$, is a strong indicator of the mass of accreted satellite galaxies. We find a strong correlation between the mass of this hot inner stellar halo and the total ex situ stellar mass. There is a similarly strong correlation with the stellar mass of the most massive secondary galaxy ever merged. These TNG50 correlations are compatible with those predicted by other simulations, for example by TNG100 across the whole mass range under study (galaxy stellar masses, $M_*$, in the $10^{10.3-11.6}$\,\Msun\, range) and by EAGLE for $M_* \gtrsim 10^{10.6} $\,\Msun\, galaxies.\ This shows that our predictions are robust across different galaxy formation and feedback models and hold across a wide range of numerical resolution. The hot inner stellar halo is a product of massive and typically ancient mergers, with inner-halo stars exhibiting three main physical origins: accreted and stripped from massive satellites, dynamically heated by mergers from the bulge and/or disk in the main progenitor, and formed from star formation triggered during mergers. The mass of the hot inner stellar halo defined in this paper is a quantity that can be robustly obtained for real galaxies by applying a population-orbit superposition method to integral-field-unit spectroscopy data, out to a distance of $\sim2\,R_e$, which is possible with current observations. Hence, this paper shows that integral-field-unit observations and dynamical models of the inner regions of galaxies provide a way to quantitatively determine the mass of ancient accreted satellites.

preprint2022arXiv

Observational Evidence of Evolving Dark Matter Profiles at $z\leq 1$

We investigate the dark matter halos of 256 star-forming disc-like galaxies at $z\sim 1$ using the KMOS redshift one spectroscopic survey (KROSS). This sample covers the redshifts $0.6 \leq z \leq 1.04$, effective radii $0.69 \leq R_e [\mathrm{kpc}] \leq 7.76$, and total stellar masses $8.7 \leq log(M_{\mathrm{star}} \ [\mathrm{M_\odot}]) \leq 11.32$. We present a mass modelling approach to study the rotation curves of these galaxies, which allow us to dynamically calculate the physical properties associated with the baryons and the dark matter halo. For the former we assume a Freeman disc, while for the latter we employ the NFW and the Burkert halo profiles, separately. At the end, we compare the results of both cases with state-of-the-art cosmological galaxy simulations (EAGLE, TNG100 and TNG50). We find that the {\em cored} dark matter halo emerged as the dominant quantity from a radius 1-3 times the effective radius. Its fraction to the total mass is in good agreement with the outcome of hydrodynamical galaxy simulations. Remarkably, we found that the dark matter core of $z\sim 1$ star-forming galaxies are smaller and denser than their local counterparts. We conclude that dark matter halos have gradually expanded over the past 6.5 Gyrs. That is, observations are capable of capturing the dark matter response to the baryonic processes (e.g. feedbacks), and thus giving us the first empirical evidence of {\em gravitational potential fluctuations} in the inner region of galaxies, which can be verified with deep surveys and future missions.

preprint2022arXiv

On the accretion of a new group of galaxies onto Virgo: II. The effect of pre-processing on the stellar population content of dEs

Using MUSE spectra, we investigate how pre-processing and accretion onto a galaxy cluster affect the integrated stellar population properties of dwarf early-type galaxies (dEs). We analyze a sample of nine dEs with stellar masses of $\rm \sim 10^9 \, M_\odot$, which were accreted ($\sim$ 2-3 Gyr ago) onto the Virgo cluster as members of a massive galaxy group. We derive their stellar population properties, namely age, metallicity ([M/H]), and the abundance ratio of $α$ elements ([$α$/Fe]), by fitting observed spectral indices with a robust, iterative procedure, and infer their star formation history (SFH) by means of full spectral fitting. We find that these nine dEs are more metal-poor (at the 2-3$σ$ level) and significantly more $α$-enhanced than dEs in the Virgo and Coma clusters with similar stellar mass, cluster-centric distance, and infall time. Moreover, for six dEs, we find evidence for a recent episode of star formation during or right after the time of accretion onto Virgo. We interpret the high [$α$/Fe] of our sample of dEs as the result of the previous exposure of these galaxies to an environment hostile to star formation, and/or the putative short burst of star formation they underwent after infall into Virgo. Our results suggest that the stellar population properties of low-mass galaxies may be the result of the combined effect of pre-processing in galaxy groups and environmental processes (such as ram-pressure triggering star formation) acting during the early phases of accretion onto a cluster.

preprint2022arXiv

The Fornax3D project: discovery of ancient massive merger events in the Fornax cluster galaxies NGC 1380 and NGC 1427

We report the discovery of ancient massive merger events in the early-type galaxies NGC 1380 and NGC 1427, members of the Fornax galaxy cluster. Both galaxies have been observed by the MUSE IFU instrument on the VLT, as part of the Fornax3D project. By fitting recently-developed population-orbital superposition models to the observed surface brightness as well as stellar kinematic, age, and metallicity maps, we obtain the stellar orbits, age and metallicity distributions of each galaxy. We then decompose each galaxy into multiple orbital-based components, including a dynamically hot inner stellar halo component which is identified as the relic of past massive mergers. By comparing to analogues from cosmological galaxy simulations, chiefly TNG50, we find that the formation of such a hot inner stellar halo requires the merger with a now-destroyed massive satellite galaxy of $3.7_{-1.5}^{+2.7} \times 10^{10}$ Msun (about $1/5$ of its current stellar mass) in the case of NGC 1380 and of $1.5_{-0.7}^{+1.6} \times10^{10}$ Msun (about $1/4$ of its current stellar mass) in the case of NGC 1427. Moreover, we infer that the last massive merger in NGC 1380 happened $\sim10$ Gyr ago based on the stellar age distribution of the re-grown dynamically cold disk, whereas the merger in NGC 1427 ended $t\lesssim 8$ Gyr ago based on the stellar populations in its hot inner stellar halo. The major merger event in NGC 1380 is the first one with both merger mass and merger time quantitatively inferred in a galaxy beyond the Local Volume. Moreover, it is the oldest and most massive merger uncovered in nearby galaxies so far.

preprint2022arXiv

The SAMI Galaxy Survey: The internal orbital structure and mass distribution of passive galaxies from triaxial orbit-superposition Schwarzschild models

Dynamical models are crucial for uncovering the internal dynamics of galaxies, however, most of the results to date assume axisymmetry, which is not representative for a significant fraction of massive galaxies. Here, we build triaxial Schwarschild orbit-superposition models of galaxies taken from the SAMI Galaxy Survey, in order to reconstruct their inner orbital structure and mass distribution. The sample consists of 161 passive galaxies with total stellar masses in the range $10^{9.5}$ to $10^{12} M_{\odot}$. We find that the changes in internal structures within 1$R_{\rm e}$ are correlated with the total stellar mass of the individual galaxies. The majority of the galaxies in the sample ($73\% \pm 3\%$) are oblate, while $19\% \pm 3\%$ are mildly triaxial and $8\% \pm 2\%$ have triaxial/prolate shape. Galaxies with $\log M_{\star}/M_{\odot} > 10.50$ are more likely to be non-oblate. We find a mean dark matter fraction of $f_{\rm{DM}} = 0.28 \pm 0.20$, within 1$R_{\rm e}$. Galaxies with higher intrinsic ellipticity (flatter) are found to have more negative velocity anisotropy $β_r$ (tangential anisotropy). $β_r$ also shows an anti-correlation with the edge-on spin parameter \lam, so that $β_r$ decreases with increasing \lam. We see evidence of an increasing fraction of hot orbits with increasing stellar mass, while warm and cold orbits show a decreasing trend. We also find that galaxies with different ($V/σ$ - $h_3$) kinematic signatures have distinct combinations of orbits. These results are in agreement with a formation scenario in which slow- and fast-rotating galaxies form through two main channels.

preprint2021arXiv

Disentangling the formation mechanisms of nuclear star clusters

Nuclear star clusters (NSCs) are massive star clusters found ubiquitously in the centres of galaxies, from the dwarf regime to massive ellipticals and spirals. The fraction of nucleated galaxies is as high as $>$ 90 % at $M_{\text{gal}} \sim 10^9 M_\odot$. However, how NSC formation mechanisms work in different regimes and what determines galaxy nucleation is still unclear. The dissipationless accretion of infalling globular clusters (GCs) and the in situ formation of stars directly at the galactic centre likely operate to grow NSCs in most galaxies; however, their efficiency has been difficult to assess observationally. Here, we provide, for the first time, a quantitative determination of the relative strength of these processes in the build-up of individual NSCs. Using a semi-analytical model of NSC formation based on the orbital evolution of inspiraling GCs, together with observed NSC and GC system properties, we derived the mass fraction of in situ born stars $f_\text{in, NSC}$ for 119 galaxies with masses from $3 \times 10^{7}$ to $3 \times 10^{11} M_\odot$, in the Local Volume, the Fornax, and Virgo galaxy clusters. Our analysis reveals that the NSC mass, as well as the ratio of NSC to the total GC system mass, are strong indicators of the dominant NSC formation channel, and not the total galaxy stellar mass as previously suggested. More massive NSCs formed predominantly via the in situ formation of stars ($f_\text{in, NSC} \sim 0.9$), while the lower-mass NSCs are expected to have formed predominantly through the merger of GCs ($f_\text{in, NSC} \sim 0.2$). The results of this simple model are in agreement with recent independent estimates of the dominant NSC formation channel from recent stellar population analysis.

preprint2021arXiv

The Fornax Cluster VLT Spectroscopic Survey III -- Kinematical characterisation of globular clusters across the Fornax galaxy cluster

The Fornax cluster provides an unparalleled opportunity to investigate in detail the formation and evolution of early-type galaxies in a dense environment. We aim at kinematically characterizing photometrically detected globular cluster (GC) candidates in the core of the cluster. We used the VLT/VIMOS spectroscopic data from the FVSS survey in the Fornax cluster, covering one square degree around the central massive galaxy NGC 1399. We confirmed a total of 777 GCs, almost doubling the previously detected GCs, using the same dataset by Pota et al. (2018). Combined with previous literature radial velocity measurements of GCs in Fornax, we compiled the most extensive spectroscopic GC sample of 2341 objects in this environment. We found that red GCs are mostly concentrated around major galaxies, while blue GCs are kinematically irregular and are widely spread throughout the core region of the cluster. The velocity dispersion profiles of blue and red GCs show a quite distinct behaviour. Blue GCs exhibit a sharp increase in the velocity dispersion profile from 250 to 400km/s within 5 arcminutes (29 kpc/1 reff of NGC 1399) from the central galaxy. The velocity dispersion profile of red GCs follows a constant value in between 200-300km/s until 8 arcminutes (46 kpc/1.6reff), and then rises to 350km/s at 10 arcminutes (58 kpc/2 reff). Beyond 10 arcminutes and out to 40 arcminutes (230 kpc/8 reff), blue and red GCs show a constant velocity dispersion of 300+/-50km/s, indicating that both GC populations are tracing the cluster potential. We have kinematically confirmed and characterized the previously photometrically discovered overdensities of intra-cluster GCs. We found that those substructured intra-cluster regions in Fornax are dominated mostly by blue GCs.

preprint2021arXiv

What to expect when using globular clusters as tracers of the total mass distribution in Milky Way-mass galaxies

Dynamical models allow us to connect the motion of a set of tracers to the underlying gravitational potential, and thus to the total (luminous and dark) matter distribution. They are particularly useful for understanding the mass and spatial distribution of dark matter (DM) in a galaxy. Globular clusters (GCs) are an ideal tracer population in dynamical models, since they are bright and can be found far out into the halo of galaxies. We aim to test how well Jeans-Anisotropic-MGE (JAM) models using GCs (positions and line-of-sight velocities) as tracers can constrain the mass and radial distribution of DM halos. For this, we use the E-MOSAICS suite of 25 zoom-in simulations of L* galaxies. We find that the DM halo properties are reasonably well recovered by the JAM models. There is, however, a strong correlation between how well we recover the mass and the radial distribution of the DM and the number of GCs in the galaxy: the constraints get exponentially worse with fewer GCs, and at least 150 GCs are needed in order to guarantee that the JAM model will perform well. We find that while the data quality (uncertainty on the radial velocities) can be important, the number of GCs is the dominant factor in terms of the accuracy and precision of the measurements. This work shows promising results for these models to be used in extragalactic systems with a sample of more than 150 GCs.

preprint2020arXiv

A deep view into the nucleus of the Sagittarius Dwarf Spheroidal Galaxy with MUSE. II. Kinematic characterization of the stellar populations

The Sagittarius dwarf spheroidal galaxy (Sgr dSph) is in an advanced stage of disruption but still hosts its nuclear star cluster (NSC), M54, at its center. In this paper, we present a detailed kinematic characterization of the three stellar populations present in M54: young metal-rich (YMR); intermediate-age metal-rich (IMR); and old metal-poor (OMP), based on the spectra of $\sim6500$ individual M54 member stars extracted from a large MUSE/VLT dataset. We find that the OMP population is slightly flattened with a low amount of rotation ($\sim0.8$ km s$^{-1}$) and with a velocity dispersion that follows a Plummer profile. The YMR population displays a high amount of rotation ($\sim5$ km s$^{-1}$) and a high degree of flattening, with a lower and flat velocity dispersion profile. The IMR population shows a high but flat velocity dispersion profile, with some degree of rotation ($\sim2$ km s$^{-1}$). We complement our MUSE data with information from \textit{Gaia DR2} and confirm that the stars from the OMP and YMR populations are comoving in 3D space, suggesting that they are dynamically bound. While dynamical evolutionary effects (e.g. energy equipartition) are able to explain the differences in velocity dispersion between the stellar populations, the strong differences in rotation indicate different formation paths for the populations, as supported by an $N$-body simulation tailored to emulate the YMR-OMP system. This study provides additional evidence for the M54 formation scenario proposed in our previous work, where this NSC formed via GC accretion (OMP) and in situ formation from gas accretion in a rotationally supported disc (YMR).

preprint2020arXiv

A discrete chemo-dynamical model of M87's globular clusters: Kinematics extending to ~ 400 kpc

We study the mass distribution and kinematics of the giant elliptical galaxy M87 (NGC 4486) using discrete chemo-dynamical, axisymmetric Jeans equation modelling. Our catalogue comprises 894 globular clusters (GCs) extending to a projected radius of $\sim 430$ kpc with line-of-sight velocities and colours, and Multi Unit Spectroscopic Explorer (MUSE) integral field unit data within the central $2.4$ kpc of the main galaxy. The gravitational potential for our models is a combination of a luminous matter potential with a varying mass-to-light ratio for the main galaxy, a supermassive black hole and a dark matter (DM) potential with a cusped or cored DM halo. The best-fitting models with either a cusped or a cored DM halo show no significant differences and both are acceptable. We obtain a total mass of $(2.16 \pm 0.38) \times 10^{13} M_{\odot}$ within $\sim$ 400 kpc. By including the stellar mass-to-light ratio gradient, the DM fraction increases from $\sim$ 26 percent (with no gradient) to $\sim$ 73 percent within $1\,R_e^{\rm maj}$ (major axis of half-light isophote, 14.2 kpc), and from $\sim$ 84 percent to $\sim$ 94 percent within $5\,R_e^{\rm maj}$ (71.2 kpc). Red GCs have moderate rotation with $V_{\rm max}/σ\sim$ 0.4, and blue GCs have weak rotation with $V_{\rm max}/σ\sim$ 0.1. Red GCs have tangential velocity dispersion anisotropy, while blue GCs are consistent with being nearly isotropic. Our results suggest that red GCs are more likely to be born in-situ, while blue GCs are more likely to be accreted.

preprint2020arXiv

A Universal fundamental plane and the $M_{dyn}-M_{\star}$ relation for galaxies with CALIFA and MaNGA

We use the stellar kinematics for $2458$ galaxies from the MaNGA survey to explore dynamical scaling relations between the stellar mass $M_{\star}$ and the total velocity parameter at the effective radius, $R_e$, defined as $S_{K}^{2}=KV_{R_e}^{2}+σ_{\star_e}^{2}$, which combines rotation velocity $V_{R_e}$, and velocity dispersion $σ_{\star_e}$. We confirm that spheroidal and spiral galaxies follow the same $M_{\star}-S_{0.5}$ scaling relation with lower scatter than the $M_{\star}-V_{R_e}$ and $M_{\star}-σ_{\star_e}$ ones. We also explore a more general Universal Fundamental Plane described by the equation $log(Υ_{e}) = log (S_{0.5}^{2}) - log (I_{e}) - log (R_{e}) + C$, which in addition to kinematics, $S_{0.5}$, and effective radius, $R_e$, it includes surface brightness, $I_e$, and dynamical mass-to-light ratio, $Υ_e$. We use sophisticated Schwarzschild dynamical models for a sub-sample of 300 galaxies from the CALIFA survey to calibrate the so called Universal Fundamental Plane. That calibration allows us to propose both: (i) a parametrization to estimate the difficult-to-measure dynamical mass-to-light ratio at the effective radius; and (ii) a new dynamical mass proxy consistent with dynamical models within $0.09\ dex$. We reproduce the relation between the dynamical mass and the stellar mass in the inner regions of galaxies. We use the estimated dynamical mass-to-light ratio from our analysis, $Υ_{e}^{fit}$, to explore the Universal Fundamental Plane with the MaNGA data set. We find that all classes of galaxies, from spheroids to disks, follow this Universal Fundamental Plane with a scatter significantly smaller $(0.05\ dex)$ than the one reported for the $M_{\star}-S_{0.5}$ relation $(0.1\ dex)$, the Fundamental Plane $(\sim 0.09\ dex)$ and comparable with Tully-Fisher studies $(\sim 0.05\ dex)$, but for a wider range of galaxy types.

preprint2020arXiv

Disentangling the formation history of galaxies via population-orbit superposition: method validation

We present population-orbit superposition models for external galaxies based on Schwarzschild's orbit-superposition method, by tagging the orbits with age and metallicity. The models fit the density distributions, as well as kinematic, age and metallicity maps from Integral Field Unit (IFU) spectroscopy observations. We validate the method and demonstrate its power by applying it to mock data, similar to those obtained by the Multi-Unit Spectroscopic Explorer (MUSE) IFU on the Very Large Telescope (VLT). These mock data are created from Auriga galaxy simulations, viewed at three different inclination angles ($\vartheta=40^o, 60^o, 80^o$). Constrained by MUSE-like mock data, our model can recover the galaxy's stellar orbit distribution projected in orbital circularity $λ_z$ vs. radius $r$, the intrinsic stellar population distribution in age $t$ vs. metallicity $Z$, and the correlation between orbits' circularity $λ_z$ and stellar age $t$. A physically motivated age-metallicity relation improves recovering the intrinsic stellar population distributions. We decompose galaxies into cold, warm and hot + counter-rotating components based on their orbit circularity distribution, and find that the surface density, mean velocity, velocity dispersion, age and metallicity maps of each component from our models well reproduce those from simulation, especially for projections close to edge-on. These galaxies exhibit strong global age vs. $σ_z$ relation, which is well recovered by our model. The method has the power to reveal the detailed build-up of stellar structures in galaxies, and offers a complement to local resolved, and high-redshift studies of galaxy evolution.

preprint2020arXiv

Distributed peer review enhanced with natural language processing and machine learning

While ancient scientists often had patrons to fund their work, peer review of proposals for the allocation of resources is a foundation of modern science. A very common method is that proposals are evaluated by a small panel of experts (due to logistics and funding limitations) nominated by the grant-giving institutions. The expert panel process introduces several issues - most notably: 1) biases introduced in the selection of the panel. 2) experts have to read a very large number of proposals. Distributed Peer Review promises to alleviate several of the described problems by distributing the task of reviewing among the proposers. Each proposer is given a limited number of proposals to review and rank. We present the result of an experiment running a machine-learning enhanced distributed peer review process for allocation of telescope time at the European Southern Observatory. In this work, we show that the distributed peer review is statistically the same as a `traditional' panel, that our machine learning algorithm can predict expertise of reviewers with a high success rate, and we find that seniority and reviewer expertise have an influence on review quality. The general experience has been overwhelmingly praised from the participating community (using an anonymous feedback mechanism).

preprint2020arXiv

Galaxies within galaxies in the TIMER survey: stellar populations of inner bars are scaled replicas of main bars

Inner bars are frequent structures in the local Universe and thought to substantially influence the nuclear regions of disc galaxies. In this study we explore the structure and dynamics of inner bars by deriving maps and radial profiles of their mean stellar population content and comparing them to previous findings in the context of main bars. To this end, we exploit observations obtained with the integral-field spectrograph MUSE of three double-barred galaxies in the TIMER sample. The results indicate that inner bars can be distinguished based on their stellar population properties alone. More precisely, inner bars show elevated metallicities and depleted [$α$/Fe] abundances. Although they exhibit slightly younger stellar ages compared to the nuclear disc, the typical age differences are small, except at their outer ends. These ends of the inner bars are clearly younger compared to their inner parts, an effect known from main bars as orbital age separation. In particular, the youngest stars (i.e. those with the lowest radial velocity dispersion) seem to occupy the most elongated orbits along the (inner) bar major axis. We speculate that these distinct ends of bars could be connected to the morphological feature of ansae. Radial profiles of metallicity and [$α$/Fe] enhancements are flat along the inner bar major axis, but show significantly steeper slopes along the minor axis. This radial mixing in the inner bar is also known from main bars and indicates that inner bars significantly affect the radial distribution of stars. In summary, based on maps and radial profiles of the mean stellar population content and in line with previous TIMER results, inner bars appear to be scaled down versions of the main bars seen in galaxies. This suggests the picture of a "galaxy within a galaxy", with inner bars in nuclear discs being dynamically equivalent to main bars in main galaxy discs.

preprint2020arXiv

Metal-poor nuclear star clusters in two dwarf galaxies near Centaurus A suggesting formation from the in-spiraling of globular clusters

Studies of nucleated dwarf galaxies can constrain the scenarios for the formation and evolution of nuclear star clusters (NSC) in low-mass galaxies and give us insights on the origin of ultra compact dwarf galaxies (UCDs). We report the discovery of a NSC in the dwarf galaxy KKs58 and investigate its properties together with those of another NSC in KK197. Both NSCs are hosted by dwarf elliptical galaxies of the Centaurus group. Combining ESO VLT MUSE data with photometry from VLT FORS2, CTIO Blanco DECam, and HST ACS, as well as high-resolution spectroscopy from VLT UVES, we analyse the photometric, kinematic and stellar population properties of the NSCs and their host galaxies. We confirm membership of the NSCs based on their radial velocities and location close to the galaxy centres. We also confirm the membership of two globular clusters (GCs) and detect oblate rotation in the main body of KK197. Based on high signal-to-noise spectra taken with MUSE of the NSCs of both KKs58 and KK197 we measure low metallicities, [Fe/H] = $-1.75 \pm 0.06$ dex and [Fe/H] = $-1.84 \pm 0.05$ dex, and stellar masses of $7.3 \times 10^5 M_\odot$ and $1.0 \times 10^6 M_\odot$, respectively. Both NSCs are more metal-poor than their hosts that have metallicities of $-1.35 \pm 0.23$ dex (KKs58) and $-0.84 \pm 0.12$ dex (KK197). This can be interpreted as NSC formation via the in-spiral of GCs. The masses, sizes and metallicities of the two NSCs place them among other NSCs, but also among the known UCDs of the Centaurus group. This indicates that NSCs might constitute the progenitors of a part of the low-mass UCDs, although their properties are almost indistinguishable from typical GCs.

preprint2020arXiv

On the accretion of a new group of galaxies onto Virgo: I. Internal kinematics of nine in-falling dEs

Galaxy environment has been shown to play an important role in transforming late-type, star-forming galaxies to quiescent spheroids. This transformation is expected to be more severe for low-mass galaxies (M < $10^{10}$ \(M_\odot\)) in dense galaxy groups and clusters, mostly due to the influence of their past host halos (also known as pre-processing) and their present-day environments. For the first time, in this study, we investigate a sample of nine early-type dwarf galaxies (dEs) that were accreted as a likely bound group onto the Virgo galaxy cluster about 2-3 Gyr ago. Considering this special condition, these nine dEs may provide a test bed for distinguishing between the influence of the Virgo galaxy cluster and the effects of the previous host halo on their current properties. Specifically, we use VLT/MUSE integral-field unit spectra to derive their kinematics and specific angular momentum ($λ_{R}$) profiles. We observe a spread in the $λ_{R}$ profiles of our sample dEs, finding that the $λ_{R}$ profiles of half of them are as high as those of low-mass field galaxies. The remaining dEs exhibit $λ_{R}$ profiles as low as those of Virgo dEs that were likely accreted longer ago. Moreover, we detect nebular emission in one dE with a gas velocity offset suggesting ongoing gas stripping in Virgo. We suggest that the low-$λ_{R}$ dEs in our sample were processed by their previous host halo, prior to their infall to Virgo, and that the high-$λ_{R}$ dEs may be experiencing ram pressure stripping in Virgo.

preprint2020arXiv

Stellar populations across galaxy bars in the MUSE TIMER project

Stellar populations in barred galaxies save an imprint of the influence of the bar on the host galaxy's evolution. We present a detailed analysis of star formation histories (SFHs) and chemical enrichment of stellar populations in nine nearby barred galaxies from the TIMER project. We use integral field observations with the MUSE instrument to derive unprecedented spatially resolved maps of stellar ages, metallicities, [Mg/Fe] abundances and SFHs, as well as H$α$ as a tracer of ongoing star formation. We find a characteristic V-shaped signature in the SFH perpendicular to the bar major axis which supports the scenario where intermediate age stars ($\sim 2$-$6\ \mathrm{Gyr}$) are trapped on more elongated orbits shaping a thinner part of the bar, while older stars ($> 8\ \mathrm{Gyr}$) are trapped on less elongated orbits shaping a rounder and thicker part of the bar. We compare our data to state-of-the-art cosmological magneto-hydrodynamical simulations of barred galaxies and show that such V-shaped SFHs arise naturally due to the dynamical influence of the bar on stellar populations with different ages and kinematic properties. Additionally, we find an excess of very young stars ($< 2\ \mathrm{Gyr}$) on the edges of the bars, predominantly on the leading side, confirming typical star formation patterns in bars. Furthermore, mass-weighted age and metallicity gradients are slightly shallower along the bar than in the disc likely due to orbital mixing in the bar. Finally, we find that bars are mostly more metal-rich and less [Mg/Fe]-enhanced than the surrounding discs. We interpret this as a signature that the bar quenches star formation in the inner region of discs, usually referred to as star formation deserts. We discuss these results and their implications on two different scenarios of bar formation and evolution.

preprint2020arXiv

The Fornax Deep Survey data release 1

We present the first data release of the Fornax Deep Survey (FDS), an imaging survey using using the wide-field imager OmegaCAM mounted on the VST in the SDSS u', g', r', and i'-bands covering the Fornax Galaxy Cluster and the infalling Fornax A Group. FDS is a joint project between NOVA (previously called FOCUS - PI: R. F. Peletier) and INAF (as part of VEGAS - PIs: M. Capaccioli and E. Iodice). With exposure times of about 9 hours over an area of ~28 square degrees, this survey is a legacy dataset for studies of members of the Fornax Galaxy Cluster and the infalling Fornax A Group down to a surface brightness limit of ~28 mag/arcsec^2 (1-sigma surface brightness over a 1 arcsecond^2 area) and opens a new parameter regime to investigate the role of the cluster environment in shaping the properties of its galaxy population. After the Virgo cluster,Fornax is the second nearest galaxy cluster to us, and with its different mass and evolutionary state, it provides a valuable comparison that makes it possible to understand the various evolutionary effects on galaxies and galaxy clusters. Details about the survey can be found in A. Venhola, R. F. Peletier, E. Laurikainen et al., 2018, A&A 620, 165. In this release, 181 Gb of (compressed) fits files reduced using the system are present. Catalogues with the complete sample of sources including dwarf galaxies part of the cluster, globular clusters, and background galaxies will be provided in forthcoming releases. The data products are available via the ESO Science Portal at https://archive.eso.org/scienceportal/home?publ_date=2020-08-26

preprint2020arXiv

The Fornax Deep Survey with VST. IX. The catalog of sources in the FDS area, with an example study for globular clusters and background galaxies

This paper continues the series of the Fornax Deep Survey (FDS). Following the previous studies dedicated to extended Fornax cluster members, we present the catalogs of compact stellar systems in the Fornax cluster as well as extended background sources and point-like sources. We derive ugri photometry of ~1.7 million sources over the $\sim$21 sq. degree area of FDS centered on NGC1399. For a wider area, of $\sim$27 sq. degs extending in the direction of NGC1316, we provide gri data for ~3.1 million sources. To improve the morphological characterization of sources we generate multi-band image stacks by coadding the best seeing gri-band single exposures with a cut at FWHM<=0.9 arcsec. We use the multi-band stacks as detection frames. The identification of compact sources is obtained from a combination of photometric and morphometric selection criteria taking as reference the properties of sources with well-defined classification from the literature. We present a preliminary analysis of globular cluster (GC) distributions in the Fornax area. The study confirms and extends further previous results. We observe the inter-galactic population of GCs, a population of mainly blue GCs centered on NGC1399, extends over $\sim$0.9Mpc, with an ellipticity $\sim$0.65. Several sub-structures extend over $\sim$0.5Mpc along various directions. Two of these structures do not cross any bright galaxy; one of them appears to be connected to NGC1404, a bright galaxy close to the cluster core and particularly poor of GCs. Using the gri catalogs we analyze the GC distribution over the extended FDS area, and do not find any obvious GC sub-structure bridging the two brightest cluster galaxies, NGC1316 and NGC1399. Although NGC1316 is twice brighter of NGC1399 in optical bands we estimate a factor of 3-4 richer GC population around NGC1399 compared to NGC1316, out to galactocentric distances of 40 arcmin

preprint2020arXiv

The SAMI -- Fornax Dwarfs Survey I: Sample, observations and the specific stellar angular momentum of dwarf elliptical galaxies

Dwarf ellipticals are the most common galaxy type in cluster environments, however the challenges associated with their observation mean their formation mechanisms are still poorly understood. To address this, we present deep integral field observations of a sample of 31 low-mass ($10^{7.5} <$ M$_\star < 10^{9.5}$ M$_\odot$) early-type galaxies in the Fornax cluster with the SAMI instrument. For 21 galaxies our observations are sufficiently deep to construct spatially resolved maps of the stellar velocity and velocity dispersion - for the remaining galaxies we extract global velocities and dispersions from aperture spectra only. From the kinematic maps we measure the specific stellar angular momentum $λ_R$ of the lowest mass dE galaxies to date. Combining our observations with early-type galaxy data from the literature spanning a large range in stellar mass, we find that $λ_R$ decreases towards lower stellar mass, with a corresponding increase in the proportion of slowly rotating galaxies in this regime. The decrease of $λ_R$ with mass in our sample dE galaxies is consistent with a similar trend seen in somewhat more massive spiral galaxies from the CALIFA survey. This suggests that the degree of dynamical heating required to produce dEs from low-mass starforming progenitors may be relatively modest, and consistent with a broad range of formation mechanisms.

preprint2019arXiv

A deep view into the nucleus of the Sagittarius Dwarf Spheroidal Galaxy with MUSE. I. Data and stellar population characterization

The center of the Sagittarius dwarf spheroidal galaxy (Sgr dSph) hosts a nuclear star cluster, M54, which is the only galaxy nucleus that can be resolved into individual stars at optical wavelengths. It is thus a key target for understanding the formation of nuclear star clusters and their relation to globular clusters. We present a large Multi-Unit Spectroscopic Explorer (MUSE) data set that covers M54 out to $\sim$2.5 half-light radius, from which we extracted the spectra of $\sim$6600 cluster member stars. We use these data in combination with HST photometry to derive age and metallicity for each star. The stellar populations show a well defined age-metallicity relation, implying an extended formation history for the central region of Sgr dSph. We classify these populations into three groups, all with the same systemic velocity: young metal-rich (YMR; 2.2\,Gyr, \mbox{[Fe/H]$=-0.04$}); intermediate-age metal-rich (IMR; 4.3\,Gyr, \mbox{[Fe/H]$=-0.29$}); and old metal-poor (OMP; 12.2\,Gyr, \mbox{[Fe/H]$=-1.41$}). The YMR and OMP populations are more centrally concentrated than the IMR population, which are likely stars of the Sgr dSph. We suggest the OMP population is the result of accretion and merging of two or more old and metal-poor globular clusters dragged to the center by dynamical friction. The YMR is consistent with being formed by in situ star formation in the nucleus. The ages of the YMR population suggest that it may have been triggered into forming when the Sgr dSph began losing its gas during the most recent interaction with the Milky Way, $\sim$3\,Gyr ago.

preprint2019arXiv

A dwarf-dwarf merger and dark matter core as a solution to the Globular Cluster problems in the Fornax dSph

The five globular clusters (GCs) of the Fornax dSph are puzzling for two reasons; the mass in GCs is high with respect to the galaxy's old stellar mass, and their survival and large distance (> 1 kpc) is at odds with naive expectations of dynamical friction. We present here a semi-analytic model, simultaneously addressing both problems in a comprehensive evolutionary framework for Fornax. Key to the model is inclusion of: 1) hydrodynamical constraints on the GC formation locations, 2) self-consistent velocity distribution functions in the dynamical friction calculations and 3) expansion of GC orbits due to a past dwarf-dwarf merger in the orbit integrations. The latter is crucial for reconciling the dynamical survival of the clusters, and their chemical properties with respect to the Fornax field stars. We find that in order for four of the GCs to survive at their observed projected location, a dark matter core of radius > 1.5 kpc and a dwarf merger with dynamical mass ratio of 1:5 to 1:2 with Fornax is required. We support the merger scenario by showing that aspects of the field star metallicity distribution function and anomalous chemical properties of GC5, are representative of a merging galaxy which is ~1/3 less massive than Fornax. Together the chemical and dynamical models suggest a scenario where three in-situ GCs in proto-Fornax were ejected to the outskirts during the merger, a GC4 formed during the merger at about 10 Gyrs ago, with GC5 being brought in by the merging galaxy to Fornax.

preprint2019arXiv

A galaxy's accretion history unveiled from its integrated spectrum

We present a new method of quantifying a galaxy's accretion history from its integrated spectrum alone. Using full spectral fitting and calibrated regularization techniques we show how we can accurately derive a galaxy's mass distribution in age-metallicity space and further separate this into stellar populations from different chemical enrichment histories. By exploiting the fact that accreted lower mass galaxies will exhibit an offset to lower metallicities at fixed age compared to the in-situ stellar population, we quantify the fraction of light that comes from past merger events, that are long since mixed in phase-space and otherwise indistinguishable. Empirical age-metallicity relations (AMRs) parameterized for different galaxy masses are used to identify the accreted stellar populations and link them back to the progenitor galaxy's stellar mass. This allows us to not only measure the host galaxy's total ex-situ mass fraction ($f_{acc}$), but also quantify the relative amount of accreted material deposited by satellite galaxies of different masses, i.e. the accreted satellite mass function in analogy to the subhalo mass function. Using mock spectra of simulated, present-day galaxies from the EAGLE suite we demonstrate that our method can recover the total accreted fraction to within $\approx 12 \%$, the stellar mass of the most massive accreted subhalo to within $\approx 26 \%$ and the slope of the accreted satellite mass function to within $\approx 16 \%$ of the true values from the EAGLE merger trees. Future application of this method to observations could potentially provide us accretion histories of hundreds of individual galaxies, for which deep integrated light spectroscopy is available.

preprint2019arXiv

Accreted Globular Clusters in External Galaxies: Why Adaptive Dynamics is not the Solution

Many astrophysical and galaxy-scale cosmological problems require a well determined gravitational potential which is often modeled by observers under strong assumptions. Globular clusters (GCs) surrounding galaxies can be used as dynamical tracers of the luminous and dark matter distribution at large (kpc) scales. A natural assumption for modeling the gravitational potential is that GCs accreted in the same dwarf galaxy merger event move at the present time on similar orbits in the host galaxy and should therefore have similar actions. We investigate this idea in one realistic Milky Way like galaxy of the cosmological N-body simulation suite Auriga. We show how the actions of accreted stellar particles in the simulation evolve and that minimizing the standard deviation of GCs in action space, however, cannot constrain the true potential. This approach known as adaptive dynamics does therefore not work for accreted GCs.

preprint2019arXiv

SDSS-IV MaNGA: Internal mass distributions and orbital structures of early-type galaxies and their dependence on environment

In our earlier 2019 paper, we evaluated the reliability of Schwarzschild's orbit-superposition dynamical modelling method in estimating the internal mass distribution, intrinsic stellar shapes and orbit distributions of early-type galaxies (ETGs) taken from the Illustris cosmological simulation. We now apply the same techniques to galaxies taken from the integral-field survey Mapping Nearby Galaxies with APO (MaNGA), using a sample of 149 ETGs in the mass range of $10^{9.90}\sim10^{11.80} M_{\odot}$ and made up of 105 central and 44 satellite galaxies. We find that low-mass ETGs with $\log(M_*/M_{\odot})<11.1$ have an average dark matter fraction of $\sim0.2$ within one effective radius $R_{\rm e}$, tend to be oblate-like, and are dominated by rotation about their minor axis. High-mass ETGs with $\log(M_*/M_{\odot})>11.1$ have an average dark matter fraction of $\sim0.4$ within one effective radius $R_{\rm e}$, tend to be prolate-like, and are dominated by rotation about their major axis and by centrophilic orbits. The changes of internal structures within one $R_{\rm e}$ are dominated by the total stellar mass of the individual galaxies. We find no differences of internal structures between central and satellite ETGs for the same stellar masses. However, for similar stellar mass and colour distributions, we find that ETGs more prolate-like, or with more hot orbits, tend to have higher close neighbour counts at $r_p\sim40$ kpc.

preprint2019arXiv

The peculiar kinematics of the multiple populations in the globular cluster Messier 80 (NGC 6093)

We combine MUSE spectroscopy and Hubble Space Telescope ultraviolet (UV) photometry to perform a study of the chemistry and dynamics of the Galactic globular cluster Messier 80 (M80, NGC 6093). Previous studies have revealed three stellar populations that not only vary in their light-element abundances, but also in their radial distributions, with concentration decreasing with increasing nitrogen enrichment. This remarkable trend, which sets M80 apart from the other Galactic globular clusters, points towards a complex formation and evolutionary history. To better understand how M80 formed and evolved, revealing its internal kinematics is key. We find that the most N-enriched population rotates faster than the other two populations at a 2 sigma confidence level. While our data further suggest that the intermediate population shows the least amount of rotation, this trend is rather marginal (1 - 2 sigma). Using axisymmetric Jeans models, we show that these findings can be explained from the radial distributions of the populations if they possess different angular momenta. Our findings suggest that the populations formed with primordial kinematical differences.

preprint2019arXiv

Using MUSE-AO observations to constrain the formation of the large nuclear star cluster in FCC47

Nuclear star clusters (NSCs) are found in at least 70% of all galaxies, but their formation path is still unclear. In the most common scenarios, NSCs form in-situ from the galaxy's central gas reservoir, through merging of globular clusters (GCs), or through a combination of the two. As the scenarios pose different expectations for angular momentum and stellar population properties of the NSC in comparison to the host galaxy and the GC system, it is necessary to characterise the stellar light, NSC, and GCs simultaneously. Wide-field observations with modern integral field units such as the Multi Unit Spectroscopic Explorer (MUSE) allow to perform such studies. However, at large distances, NSCs usually are not resolved in MUSE observations. The particularly large NSC (R_eff $\sim$ 66 pc) of the early-type galaxy FCC47 at distance of $\sim$ 20 Mpc is an exception and is therefore an ideal laboratory to constrain NSC formation of external galaxies.

preprint2016arXiv

A discrete chemo-dynamical model of the dwarf spheroidal galaxy Sculptor: mass profile, velocity anisotropy and internal rotation

We present a new discrete chemo-dynamical axisymmetric modeling technique, which we apply to the dwarf spheroidal galaxy Sculptor. The major improvement over previous Jeans models is that realistic chemical distributions are included directly in the dynamical modelling of the discrete data. This avoids loss of information due to spatial binning and eliminates the need for hard cuts to remove contaminants and to separate stars based on their chemical properties. Using a combined likelihood in position, metallicity and kinematics, we find that our models naturally separate Sculptor stars into a metal-rich and a metal-poor population. Allowing for non-spherical symmetry, our approach provides a central slope of the dark matter density of $γ= 0.5 \pm 0.3$. The metal-rich population is nearly isotropic (with $β_r^{red} = 0.0\pm0.1$) while the metal-poor population is tangentially anisotropic (with $β_r^{blue} = -0.2\pm0.1$) around the half light radius of $0.26$ kpc. A weak internal rotation of the metal-rich population is revealed with $v_{max}/σ_0 = 0.15 \pm 0.15$. We run tests using mock data to show that a discrete dataset with $\sim 6000$ stars is required to distinguish between a core ($γ= 0$) and cusp ($γ= 1$), and to constrain the possible internal rotation to better than $1\,σ$ confidence with our model. We conclude that our discrete chemo-dynamical modelling technique provides a flexible and powerful tool to robustly constrain the internal dynamics of multiple populations, and the total mass distribution in a stellar system.

preprint2016arXiv

A discrete chemo-dynamical model of the giant elliptical galaxy NGC 5846: dark matter fraction, internal rotation and velocity anisotropy out to six effective radii

We construct a suite of discrete chemo-dynamical models of the giant elliptical galaxy NGC 5846. These models are a powerful tool to constrain both the mass distribution and internal dynamics of multiple tracer populations. We use Jeans models to simultaneously fit stellar kinematics within the effective radius $R_{\rm e}$, planetary nebula (PN) radial velocities out to $3\, R_{\rm e}$, and globular cluster (GC) radial velocities and colours out to $6\,R_{\rm e}$. The best-fitting model is a cored DM halo which contributes $\sim 10\%$ of the total mass within $1\,R_{\rm e}$, and $67\% \pm 10\%$ within $6\,R_{\rm e}$, although a cusped DM halo is also acceptable. The red GCs exhibit mild rotation with $v_{\rm max}/σ_0 \sim 0.3$ in the region $R > \,R_{\rm e}$, aligned with but counter-rotating to the stars in the inner parts, while the blue GCs and PNe kinematics are consistent with no rotation. The red GCs are tangentially anisotropic, the blue GCs are mildly radially anisotropic, and the PNe vary from radially to tangentially anisotropic from the inner to the outer region. This is confirmed by general made-to-measure models. The tangential anisotropy of the red GCs in the inner regions could stem from the preferential destruction of red GCs on more radial orbits, while their outer tangential anisotropy -- similar to the PNe in this region -- has no good explanation. The mild radial anisotropy of the blue GCs is consistent with an accretion scenario.

preprint2016arXiv

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

We present our mass estimate of the central black hole in the isolated spiral galaxy NGC 4414. Using natural guide star adaptive optics assisted observations with the Gemini Near-Infrared Integral Field Spectrometer (NIFS) and the natural seeing Gemini Multi-Object Spectrographs-North (GMOS), we derived two-dimensional stellar kinematic maps of NGC 4414 covering the central 1.5 arcsec and 10 arcsec, respectively, at a NIFS spatial resolution of 0.13 arcsec. The kinematic maps reveal a regular rotation pattern and a central velocity dispersion dip down to around 105 km/s. We constructed dynamical methods using two different methods: Jeans anisotropic dynamical modeling and axisymmetric Schwarzschild modeling. Both modeling methods give consistent results, but we cannot constrain the lower mass limit and only measure an upper limit for the black hole mass of Mbh= 1.56 x 10^6 Msun(at 3 sigma level) which is at least 1 sigma below the recent Mbh-sigma_e relations. Further tests with dark matter, mass-to-light ratio variation and different light models confirm that our results are not dominated by uncertainties. The derived upper limit mass is not only below the Mbh-sigma_e relation, but is also five times lower than the lower limit black hole mass anticipated from the resolution limit of the sphere of influence. This proves that via high quality integral field data we are now able to push black hole measurements down to at least five times less than the resolution limit.

preprint2016arXiv

A spiral galaxy's mass distribution uncovered through lensing and dynamics

We investigate the matter distribution of a spiral galaxy with a counter-rotating stellar core, SDSS J1331+3628 (J1331), independently with gravitational lensing and stellar dynamical modelling. By fitting a gravitational potential model to a quadruplet of lensing images around J1331's bulge, we tightly constrain the mass inside the Einstein radius R_ein = (0.91 +/- 0.02)'' (~= 1.83 +/- 0.04 kpc) to within 4%: M_ein = (7.8 +/- 0.3) x 10^10 M_Sun. We model observed long-slit major axis stellar kinematics in J1331's central regions by finding Multi-Gaussian Expansion (MGE) models for the stellar and dark matter distribution that solve the axisymmetric Jeans equations. The lens and dynamical model are independently derived, but in very good agreement with each other around ~R_ein. We find that J1331's center requires a steep total mass-to-light ratio gradient. A dynamical model including an NFW halo (with virial velocity v_200 ~= 240 +/- 40 km/s and concentration c_200 ~= 8 +/- 2) and moderate tangential velocity anisotropy (beta_z ~= -0.4 +/- 0.1) can reproduce the signatures of J1331's counter-rotating core and predict the stellar and gas rotation curve at larger radii. However, our models do not agree with the observed velocity dispersion at large radii. We speculate that the reason could be a non-trivial change in structure and kinematics due to a possible merger event in J1331's recent past.

preprint2016arXiv

Being WISE II: Reducing the Influence of Star Formation History on the Mass-to-Light Ratio of Quiescent Galaxies

Stellar population synthesis models can now reproduce the photometry of old stellar systems (age $>$ 2 Gyr) in the near-infrared (NIR) bands at 3.4 and 4.6$μ$m (WISE W1 $\&$ W2 or IRAC 1 $\&$ 2). In this paper we derive stellar mass-to-light ratios for these and optical bands, and confirm that the NIR M/L shows dramatically reduced sensitivity to both age and metallicity compared to optical bands, and further, that this behavior leads to significantly more robust stellar masses for quiescent galaxies with [Fe/H] > -0.5 regardless of star formation history (SFH). We then use realistic early-type galaxy SFHs and metallicity distributions from the EAGLE simulations of galaxy formation to investigate two methods to determine the appropriate M/L for a galaxy: 1) We show that the uncertainties introduced by an unknown SFH can be largely removed using a spectroscopically inferred luminosity-weighted age and metallicity for the population to select the appropriate single stellar population (SSP) equivalent M/L. Using this method, the maximum systematic error due to SFH on the M/L of an early-type galaxy is $<$ 4$\%$ at 3.4 $μ$m and typical uncertainties due to errors in the age and metallicity create scatter of $\lesssim$13$\%$. The equivalent values for optical bands are more than 2-3 times greater, even before considering uncertainties associated with internal dust extinction. 2) We demonstrate that if the EAGLE SFHs and metallicities accurately reproduce the true properties of early-type galaxies, the use of an iterative approach to select a mass dependent M/L can provide even more accurate stellar masses for early-type galaxies, with typical uncertainties $<$ 9$\%$.

preprint2016arXiv

H$α$ Imaging spectroscopy of Balmer-dominated shocks in Tycho's supernova remnant

We present Fabry-Pérot interferometric observations of the narrow H$α$ component in the shock front of the historical supernova remnant Tycho (SN 1572). Using GH$α$FaS (Galaxy H$α$ Fabry-Pérot Spectrometer) on the William Herschel Telescope, we observed a great portion of the shock front in the northeastern (NE) region of the remnant. The angular resolution of $\sim$1$^{\prime\prime}$ and spectral resolving power of R$\sim$21 000 together with the large field-of-view (3.4$^{\prime}$ $\times$ 3.4$^{\prime}$) of the instrument allow us to measure the narrow H$α$-line width in 73 bins across individual parts of the shock simultaneously and thereby study the indicators of several shock precursors in a large variety of shock front conditions. Compared to previous studies, the detailed spatial resolution of the filament also allows us to mitigate possible artificial broadening of the line from unresolved differential motion and projection. Covering one quarter of the remnant's shell, we confirm the broadening of the narrow H$α$ line beyond its intrinsic width of $\sim$20 km/s and report it to extend over most of the filament, not only the previously investigated dense 'knot g'. Similarly, we confirm and find additional strong evidence for wide-spread intermediate-line ($\sim$150 km/s) emission. Our Bayesian analysis approach allows us to quantify the evidence for this intermediate component as well as a possible split in the narrow line. Suprathermal narrow line widths point toward an additional heating mechanism in the form of a cosmic-ray precursor, while the intermediate component, previously only qualitatively reported as a small non-Gaussian contribution to the narrow component, reveals a broad-neutral precursor.

preprint2016arXiv

The Black Hole - Bulge Mass Relation in Megamaser Host Galaxies

We present HST images for nine megamaser disk galaxies with the primary goal of studying photometric BH-galaxy scaling relations. The megamaser disks provide the highest-precision extragalactic BH mass measurements, while our high-resolution HST imaging affords us the opportunity to decompose the complex nuclei of their late-type hosts in detail. Based on the morphologies and shapes of the galaxy nuclei, we argue that most of these galaxies' central regions contain secularly evolving components (pseudo-bulges), and in many cases we photometrically identify co-existing "classical" bulge components as well. Using these decompositions, we draw the following conclusions: (1) The megamaser BH masses span two orders of magnitude ($10^6$ -- $10^8 M_\odot$) while the stellar mass of their spiral host galaxies are all $\sim 10^{11} M_\odot$ within a factor of three; (2) the BH masses at a given bulge mass or total stellar mass in the megamaser host spiral galaxies tend to be lower than expected, when compared to an extrapolation of the BH-bulge relation based on early-type galaxies; (3) the observed large intrinsic scatter of BH masses in the megamaser host galaxies raises the question of whether scaling relations exist in spiral galaxies.

preprint2016arXiv

The inner mass distribution of late-type spiral galaxies from SAURON stellar kinematic maps

We infer the central mass distributions within 0.4-1.2 disc scale lengths of 18 late-type spiral galaxies using two different dynamical modelling approaches - the Asymmetric Drift Correction (ADC) and axisymmetric Jeans Anisotropic Multi-gaussian expansion (JAM) model. ADC adopts a thin disc assumption, whereas JAM does a full line-of-sight velocity integration. We use stellar kinematics maps obtained with the integral-field spectrograph SAURON to derive the corresponding circular velocity curves from the two models. To find their best-fit values, we apply Markov Chain Monte Carlo (MCMC) method. ADC and JAM modelling approaches are consistent within 5% uncertainty when the ordered motions are significant comparable to the random motions, i.e, $\overline{v_ϕ}/σ_R$ is locally greater than 1.5. Below this value, the ratio $v_\mathrm{c,JAM}/v_\mathrm{c,ADC}$ gradually increases with decreasing $\overline{v_ϕ}/σ_R$, reaching $v_\mathrm{c,JAM}\approx 2 \times v_\mathrm{c,ADC}$. Such conditions indicate that the stellar masses of the galaxies in our sample are not confined to their disk planes and likely have a non-negligible contribution from their bulges and thick disks.

preprint2015arXiv

A New Channel for the Formation of Kinematically Decoupled Cores in Early-type galaxies

We present the formation of a Kinematically Decoupled Core (KDC) in an elliptical galaxy, resulting from a major merger simulation of two disk galaxies. We show that although the two progenitor galaxies are initially following a prograde orbit, strong reactive forces during the merger can cause a short-lived change of their orbital spin; the two progenitors follow a retrograde orbit right before their final coalescence. This results in a central kinematic decoupling and the formation of a large-scale (~2 kpc radius) counter-rotating core (CRC) at the center of the final elliptical-like merger remnant (M*=1.3x10^11 Msun), while its outer parts keep the rotation direction of the initial orbital spin. The stellar velocity dispersion distribution of the merger remnant galaxy exhibits two symmetrical off-centered peaks, comparable to the observed "2-sigma galaxies". The KDC/CRC consists mainly of old, pre-merger population stars (older than 5 Gyr), remaining prominent in the center of the galaxy for more than 2 Gyr after the coalescence of its progenitors. Its properties are consistent with KDCs observed in massive elliptical galaxies. This new channel for the formation of KDCs from prograde mergers is in addition to previously known formation scenarios from retrograde mergers and can help towards explaining the substantial fraction of KDCs observed in early-type galaxies.

preprint2015arXiv

An Integral View of Fast Shocks around Supernova 1006

Supernova remnants are among the most spectacular examples of astrophysical pistons in our cosmic neighborhood. The gas expelled by the supernova explosion is launched with velocities ~ 1000 km/s into the ambient, tenuous interstellar medium, producing shocks that excite hydrogen lines. We have used an optical integral-field spectrograph to obtain high-resolution spatial-spectral maps that allow us to study in detail the shocks in the northwestern rim of supernova 1006. The two-component Hα line is detected at 133 sky locations. Variations in the broad line widths and the broad-to-narrow line intensity ratios across tens of atomic mean free paths suggest the presence of suprathermal protons, the potential seed particles for generating high-energy cosmic-rays.

preprint2015arXiv

Creating S0s with major mergers: a 3D view

A number of simulators have argued that major mergers can sometimes preserve discs (e.g. Springel & Hernquist 2005), but the possibility that they could explain the emergence of lenticular galaxies (S0s) has been generally neglected. In fact, observations of S0s reveal a strong structural coupling between their bulges and discs, which seems difficult to reconcile with the idea that they come from major mergers. However, in Querejeta et al. (2015a) we have used N-body simulations of binary mergers to show that, under favourable conditions, discs are first destroyed but soon regrow out of the leftover debris, matching observational photometric scaling relations (e.g. Laurikainen et al. 2010). Additionally, in Querejeta et al. (2015b) we have shown how the merger scenario agrees with the recent discovery that S0s and most spirals are not compatible in an angular momentum--concentration plane. This important result from CALIFA constitutes a serious objection to the idea that spirals transform into S0s mainly by fading (e.g. via ram-pressure stripping, as that would not explain the observed simultaneous change in $λ_\mathrm{Re}$ and concentration), but our simulations of major mergers do explain that mismatch. From such a 3D comparison we conclude that mergers must be a relevant process in the build-up of the current population of S0s.

preprint2015arXiv

Formation of S0 galaxies through mergers: Explaining angular momentum and concentration change from spirals to S0s

The CALIFA team has recently found that the stellar angular momentum and concentration of late-type spiral galaxies are incompatible with those of lenticular galaxies (S0s), concluding that fading alone cannot satisfactorily explain the evolution from spirals into S0s. Here we explore whether major mergers can provide an alternative way to transform spirals into S0s by analysing the spiral-spiral major mergers from the GalMer database that lead to realistic, relaxed S0-like galaxies. We find that the change in stellar angular momentum and concentration can explain the differences in the $λ_\mathrm{Re}$--$R_{90}/R_{50}$ plane found by the CALIFA team. Major mergers thus offer a feasible explanation for the transformation of spirals into S0s.

preprint2015arXiv

IMF - metallicity: a tight local relation revealed by the CALIFA survey

Variations in the stellar initial mass function (IMF) have been invoked to explain the spectroscopic and dynamical properties of early-type galaxies. However, no observations have yet been able to disentangle the physical driver. We analyse here a sample of 24 early-type galaxies drawn from the CALIFA survey, deriving in a homogeneous way their stellar population and kinematic properties. We find that the local IMF is tightly related to the local metallicity, becoming more bottom-heavy towards metal-rich populations. Our result, combined with the galaxy mass-metallicity relation, naturally explains previous claims of a galaxy mass-IMF relation, derived from non-IFU spectra. If we assume that - within the star formation environment of early-type galaxies - metallicity is the main driver of IMF variations, a significant revision of the interpretation of galaxy evolution observables is necessary.

preprint2015arXiv

MRK 1216 & NGC 1277 - An orbit-based dynamical analysis of compact, high velocity dispersion galaxies

We present a dynamical analysis to infer the structural parameters and properties of the two nearby, compact, high velocity dispersion galaxies MRK1216 & NGC1277. Combining deep HST imaging, wide-field IFU stellar kinematics, and complementary long-slit spectroscopic data out to 3 R_e, we construct orbit-based models to constrain their black hole masses, dark matter content and stellar mass-to-light ratios. We obtain a black hole mass of log(Mbh/Msun) = 10.1(+0.1/-0.2) for NGC1277 and an upper limit of log(Mbh/Msun) = 10.0 for MRK1216, within 99.7 per cent confidence. The stellar mass-to-light ratios span a range of Upsilon_V = 6.5(+1.5/-1.5) in NGC1277 and Upsilon_H = 1.8(+0.5/-0.8) in MRK1216 and are in good agreement with SSP models of a single power-law Salpeter IMF. Even though our models do not place strong constraints on the dark halo parameters, they suggest that dark matter is a necessary ingredient in MRK1216, with a dark matter contribution of 22(+30/-20) per cent to the total mass budget within 1 R_e. NGC1277, on the other hand, can be reproduced without the need for a dark halo, and a maximal dark matter fraction of 13 per cent within the same radial extent. In addition, we investigate the orbital structures of both galaxies, which are rotationally supported and consistent with photometric multi-Sérsic decompositions, indicating that these compact objects do not host classical, non-rotating bulges formed during recent (z <= 2) dissipative events or through violent relaxation. Finally, both MRK 1216 and NGC 1277 are anisotropic, with a global anisotropy parameter delta of 0.33 and 0.58, respectively. While MRK 1216 follows the trend of fast-rotating, oblate galaxies with a flattened velocity dispersion tensor in the meridional plane of the order of beta_z = delta, NGC 1277 is highly tangentially anisotropic and seems to belong kinematically to a distinct class of objects.

preprint2015arXiv

The tilt of the velocity ellipsoid in the Milky Way disk

Accurate determination of the local dark matter density is important for understanding the nature and distribution of dark matter in the universe. This requires that the local velocity distribution is characterised correctly. Here, we present a kinematic study of 16,276 SEGUE G-type dwarf stars in the solar neighbourhood, with which we determine the shape of the velocity ellipsoid in the meridional plane. We separate our G-dwarf stars based on their [Fe/H] and [alpha/Fe] abundances and infer the local velocity distribution independently for each sub-sample using a maximum-likelihood method that accounts for possible contaminants. We show by constructing vertical Jeans models that the different sub-samples yield consistent results only when we allow the velocity ellipsoid in the disk to be tilted, demonstrating that the common assumption of decoupled radial and vertical motions in the disk is incorrect. Further, we obtain that the tilt of the velocity ellipsoid is consistent among the different sub-samples. We find that increase in the tilt with height is well described by the relation alpha_tilt = (-0.90 +- 0.04) arctan(|z|/R_sun) - (0.01 +- 0.005), which is close to alignment with the spherical coordinate system and hence a velocity ellipsoid pointing to the Galactic centre. We also confirm earlier findings that the sub-samples behave almost isothermally with both radial and vertical velocity dispersion approximately constant with height. We conclude that the coupling between radial and vertical motion captured in the velocity ellipsoid tilt cannot be ignored when considering dynamical models of the solar neighbourhood. In a subsequent paper, we will develop a new modelling scheme informed by these results and make an improved determination of the local dark matter density.

preprint2015arXiv

Understanding the central kinematics of globular clusters with simulated integrated-light IFU observations

The detection of intermediate mass black holes in the centres of globular clusters is highly controversial, as complementary observational methods often deliver significantly different results. In order to understand these discrepancies, we develop a procedure to simulate integral field unit (IFU) observations of globular clusters: Simulating IFU Star Cluster Observations (SISCO). The input of our software are realistic dynamical models of globular clusters that are then converted in a spectral data cube. We apply SISCO to Monte Carlo cluster simulations from Downing et al. (2010), with a realistic number of stars and concentrations. Using independent realisations of a given simulation we are able to quantify the stochasticity intrinsic to the problem of observing a partially resolved stellar population with integrated-light spectroscopy. We show that the luminosity-weighted IFU observations can be strongly biased by the presence of a few bright stars that introduce a scatter in the velocity dispersion measurements up to $\simeq$40% around the expected value, preventing any sound assessment of the central kinematic and a sensible interpretation of the presence/absence of an intermediate mass black hole. Moreover, we illustrate that, in our mock IFU observations, the average kinematic tracer has a mass of $\simeq$0.75 solar masses, only slightly lower than the mass of the typical stars examined in studies of resolved line-of-sight velocities of giant stars. Finally, in order to recover unbiased kinematic measurements we test different masking techniques that allow us to remove the spaxels dominated by bright stars, bringing the scatter down to a level of only a few percent. The application of SISCO will allow to investigate state-of-the-art simulations as realistic observations.

preprint2015arXiv

Virgo cluster and field dwarf ellipticals in 3D: III. Spatially and temporally resolved stellar populations

We present the stellar population analysis of a sample of 12 dwarf elliptical galaxies, observed with the SAURON integral field unit, using the full-spectrum fitting method. We show that star formation histories (SFHs) resolved into two populations can be recovered even within a limited wavelength range, provided that high S/N data is used. We confirm that dEs have had complex SFHs, with star formation extending to (more) recent epochs: for the majority of our galaxies star formation activity was either still strong a few ($\lesssim$ 5) Gyr ago or they experienced a secondary burst of star formation roughly at that time. This latter possibility is in agreement with the proposed dE formation scenario where tidal harassment drives the gas remaining in their progenitors inwards and induces a star formation episode. For one of our field galaxies, ID0918, we find a correlation between its stellar population and kinematic properties, pointing to a possible merger origin of its kinematically-decoupled core. One of our cluster objects, VCC1431, appears to be composed exclusively of an old population ($\gtrsim$ 10-12 Gyr). Combining this with our earlier dynamical results, we conclude that the galaxy was either ram-pressure stripped early on in its evolution in a group environment and subsequently tidally heated, or that it evolved in situ in the cluster's central parts, compact enough to avoid tidal disruption. These are only two of the examples illustrating the SFH richness of these objects confirmed with our data.

preprint2014arXiv

Angular Momentum across the Hubble sequence from the CALIFA survey

We investigate the stellar angular momentum of galaxies across the Hubble sequence from the CALIFA survey. The distribution of CALIFA elliptical and lenticular galaxies in the $λ_{\rm Re}-ε_{\rm e}$ diagram is consistent with that shown by the Atlas$^\mathrm{3D}$ survey. Our data, however, show that the location of spiral galaxies in this diagram is significantly different. We have found two families of spiral galaxies with particularly peculiar properties: (a) spiral galaxies with much higher $λ_{\rm Re}$ values than any elliptical and lenticular galaxy; (b) low-mass spiral galaxies with observed $λ_{\rm Re}$ values much lower than expected for their apparent flattening. We use these two families of objects to argue that (1) fading alone cannot explain the transformation of spiral to lenticular galaxies, and (2) that those low-mass spiral galaxies are in fact dark matter dominated, which explains the unusually low angular momentum.

preprint2014arXiv

Being WISE I: Validating Stellar Population Models and M/L ratios at 3.4 and 4.6 microns

Using data from the WISE mission, we have measured near infra-red (NIR) photometry of a diverse sample of dust-free stellar systems (globular clusters, dwarf and giant early-type galaxies) which have metallicities that span the range -2.2 < [Fe/H] (dex) < 0.3. This dramatically increases the sample size and broadens the metallicity regime over which the 3.4 (W1) and 4.6 micron (W2) photometry of stellar populations have been examined. We find that the W1 - W2 colors of intermediate and old (> 2 Gyr) stellar populations are insensitive to the age of the stellar population, but that the W1 - W2 colors become bluer with increasing metallicity, a trend not well reproduced by most stellar population synthesis (SPS) models. In common with previous studies, we attribute this behavior to the increasing strength of the CO absorption feature located in the 4.6 micron bandpass with metallicity. Having used our sample to validate the efficacy of some of the SPS models, we use these models to derive stellar mass-to-light ratios in the W1 and W2 bands. Utilizing observational data from the SAURON and ATLAS3D surveys, we demonstrate that these bands provide extremely simple, yet robust stellar mass tracers for dust free older stellar populations that are freed from many of the uncertainties common among optical estimators.

preprint2014arXiv

Circumnuclear Molecular Gas in Megamaser Disk Galaxies NGC 4388 and NGC 1194

We explore the warm molecular and ionized gas in the centers of two megamaser disk galaxies using K-band spectroscopy. Our ultimate goal is to determine how gas is funneled onto the accretion disk, here traced by megamaser spots on sub-pc scales. We present NIR IFU data with a resolution of ~50 pc for two galaxies: NGC 4388 with VLT/SINFONI and NGC 1194 with Keck/OSIRIS+AO. The high spatial resolution and rich spectral diagnostics allow us to study both the stellar and gas kinematics as well as gas excitation on scales only an order of magnitude larger than the maser disk. We find a drop in the stellar velocity dispersion in the inner ~100 pc of NGC 4388, a common signature of a dynamically cold central component seen in many active nuclei. We also see evidence for non-circular gas motions in the molecular hydrogen on similar scales, with the gas kinematics on 100-pc scales aligned with the megamaser disk. In contrast, the high ionization lines and Br-gamma trace outflow along the 100 pc-scale jet. In NGC 1194, the continuum from the accreting black hole is very strong, making it difficult to measure robust two-dimensional kinematics, but the spatial distribution and line ratios of the molecular hydrogen and Br-gamma have consistent properties between the two galaxies.

preprint2014arXiv

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 2/7 galaxies show some localized differences in their rotation curves between the two tracers, estimates of the central logarithmic slope of the DM density profile, gamma, are generally robust. The mean and standard deviation of the logarithmic slope for the population are gamma=0.67+/-0.10 when measured in the stars and gamma=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. 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. Determining the importance of these correlations will require further model developments and larger observational samples. (Abridged)

preprint2014arXiv

Reconstructing the stellar mass distributions of galaxies using S4G IRAC 3.6 and 4.5 micron images: II. The conversion from light to mass

We present a new approach for estimating the 3.6 micron stellar mass-to-light ratio in terms of the [3.6]-[4.5] colors of old stellar populations. Our approach avoids several of the largest sources of uncertainty in existing techniques. By focusing on mid-IR wavelengths, we gain a virtually dust extinction-free tracer of the old stars, avoiding the need to adopt a dust model to correctly interpret optical or optical/NIR colors normally leveraged to assign M/L. By calibrating a new relation between NIR and mid-IR colors of GLIMPSE giant stars we also avoid discrepancies in model predictions for the [3.6]-[4.5] colors of old stellar populations due to uncertainties in molecular line opacities. We find that the [3.6]-[4.5] color, which is driven primarily by metallicity, provides a tight constraint on M/L_3.6, which varies intrinsically less than at optical wavelengths. The uncertainty on M/L_3.6 of ~0.07 dex due to unconstrained age variations marks a significant improvement on existing techniques for estimating the stellar M/L with shorter wavelength data. A single M/L_3.6=0.6 (assuming a Chabrier IMF), independent of [3.6]-[4.5] color, is also feasible as it can be applied simultaneously to old, metal-rich and young, metal-poor populations, and still with comparable (or better) accuracy (~0.1 dex) as alternatives. We expect our M/L_3.6 to be optimal for mapping the stellar mass distributions in S4G galaxies, for which we have developed an Independent Component Analysis technique to first isolate the old stellar light at 3.6 micron from non-stellar emission (e.g. hot dust and the 3.3 PAH feature). Our estimate can also be used to determine the fractional contribution of non-stellar emission to global (rest-frame) 3.6 micron fluxes, e.g. in WISE imaging, and establishes a reliable basis for exploring variations in the stellar IMF.

preprint2014arXiv

The Mice at play in the CALIFA survey: A case study of a gas-rich major merger between first passage and coalescence

We present optical integral field spectroscopy (IFS) observations of the Mice, a major merger between two massive (>10^11Msol) gas-rich spirals NGC4676A and B, observed between first passage and final coalescence. The spectra provide stellar and gas kinematics, ionised gas properties and stellar population diagnostics, over the full optical extent of both galaxies. The Mice provide a perfect case study highlighting the importance of IFS data for improving our understanding of local galaxies. The impact of first passage on the kinematics of the stars and gas has been significant, with strong bars likely induced in both galaxies. The barred spiral NGC4676B exhibits a strong twist in both its stellar and ionised gas disk. On the other hand, the impact of the merger on the stellar populations has been minimal thus far: star formation induced by the recent close passage has not contributed significantly to the global star formation rate or stellar mass of the galaxies. Both galaxies show bicones of high ionisation gas extending along their minor axes. In NGC4676A the high gas velocity dispersion and Seyfert-like line ratios at large scaleheight indicate a powerful outflow. Fast shocks extend to ~6.6kpc above the disk plane. The measured ram pressure and mass outflow rate (~8-20Msol/yr) are similar to superwinds from local ULIRGs, although NGC4676A has only a moderate infrared luminosity of 3x10^10Lsol. Energy beyond that provided by the mechanical energy of the starburst appears to be required to drive the outflow. We compare the observations to mock kinematic and stellar population maps from a merger simulation. The models show little enhancement in star formation during and following first passage, in agreement with the observations. We highlight areas where IFS data could help further constrain the models.

preprint2013arXiv

ALMA follows streaming of dense gas down to 40 pc from the supermassive black hole in NGC1097

We present a kinematic analysis of the dense molecular gas in the central 200 parsecs of the nearby galaxy NGC1097, based on Cycle 0 observations with the Atacama Large Millimeter/sub-millimeter Array (ALMA). We use the HCN(4-3) line to trace the densest interstellar molecular gas, and quantify its kinematics, and estimate an inflow rate for the molecular gas. We find a striking similarity between the ALMA kinematic data and the analytic spiral inflow model that we have previously constructed based on ionized gas velocity fields on larger scales. We are able to follow dense gas streaming down to 40 pc distance from the supermassive black hole in this Seyfert 1 galaxy. In order to fulfill marginal stability, we deduce that the dense gas is confined to a very thin disc, and we derive a dense gas inflow rate of 0.09 Msun/yr at 40 pc radius. Combined with previous values from the Ha and CO gas, we calculate a combined molecular and ionized gas inflow rate of 0.2 Msun/yr at 40 pc distance from the central supermassive black hole of NGC1097.

preprint2013arXiv

An Integral View of Balmer-dominated Shocks in Supernova Remnants

We present integral-field spectroscopic observations with the VIMOS-IFU at the VLT of fast (2000-3000 km/s) Balmer-dominated shocks surrounding the northwestern rim of the remnant of supernova 1006. The high spatial and spectral resolution of the instrument enable us to show that the physical characteristics of the shocks exhibit a strong spatial variation over few atomic scale lengths across 133 sky locations. Our results point to the presence of a population of non-thermal protons (10-100 keV) which might well be the seed particles for generating high-energy cosmic rays. We also present observations of Tycho's supernova remnant taken with the narrow-band tunable filter imager OSIRIS at the GTC and the Fabry-Perot interferometer GHaFaS at the WHT to resolve respectively the broad and narrow Hα lines across a large part of the remnant.

preprint2013arXiv

Bottom-heavy initial mass function in a nearby compact L*-galaxy

We present orbit-based dynamical models and stellar population analysis of galaxy SDSS J151741.75-004217.6, a low-redshift (z=0.116) early-type galaxy (ETG) which, for its moderate luminosity, has an exceptionally high velocity dispersion. We aim to determine the central black hole mass (M_bh), the i-band stellar mass-to-light ratio, and the low-mass slope of the initial mass function (IMF). Combining constraints from HST imaging and longslit kinematic data with those from fitting the SDSS spectrum with stellar populations models of varying IMF, we show that this galaxy has a large fraction of low-mass stars, significantly higher than implied even by a Salpeter IMF. We exclude a Chabrier/Kroupa as well as a unimodal (i.e. single-segment) IMF, while a bimodal (low-mass tapered) shape is consistent with the dynamical constraints. Thereby, our study demonstrates that a very bottom-heavy IMF can exist even in an L* ETG. We place an upper limit of ~10^{10.5}M_sun on M_bh, which still leaves open the possibility of an extremely massive BH.

preprint2013arXiv

Constraining the Galactic potential via action-based distribution functions for mono-abundance stellar populations

We present a rigorous and practical way of constraining the Galactic potential based on the phase-space information for many individual stars. Such an approach is needed to dynamically model the data from ongoing spectroscopic surveys of the Galaxy and in the future Gaia. This approach describes the orbit distribution of stars by a family of parametrized distribution function (DF) proposed by McMillan and Binney, which are based on actions. We find that these parametrized DFs are flexible enough to capture well the observed phase-space distributions of individual abundance-selected Galactic subpopulations of stars (`mono-abundance populations') for a disc-like gravitational potential, which enables independent dynamical constraints from each of the Galactic mono-abundance populations. We lay out a statistically rigorous way to constrain the Galactic potential parameters by constructing the joint likelihood of potential and DF parameters, and subsequently marginalizing over the DF parameters. This approach explicitly incorporates the spatial selection function inherent to all Galactic surveys, and can account for the uncertainties of the individual position--velocity observations. On that basis, we study the precision of the parameters of the Galactic potential that can be reached with various sample sizes and realistic spatial selection functions. By creating mock samples from the DF, we show that, even under a restrictive and realistic spatial selection function, given a two-parameter gravitational potential, one can recover the true potential parameters to a few per cent with sample sizes of a few thousands. The assumptions of axisymmetry, of DFs that are smooth in the actions and of no time variation remain important limitations in our current study.

preprint2013arXiv

Discrete dynamical modelling of omega Centauri

We present a new framework for modelling discrete kinematic data. Current techniques typically involve binning. Our approach works directly with the discrete data and uses maximum-likelihood methods to assess the probability of the dataset given model predictions. We avoid making hard cuts on the datasets by allowing for a contaminating population in our models. We apply our models to discrete proper-motion and line-of-sight-velocity data of Galactic globular cluster omega Centauri and find a mildly radial velocity anisotropy beta = 0.10 +/- 0.02, an inclination angle i = 50 +/- 1 deg, a V-band mass-to-light ratio Upsilon = 2.71 +/- 0.05 Msun/Lsun and a distance d = 4.59 +/- 0.08 kpc. All parameters are in agreement with previous studies, demonstrating the feasibility of our methods. We find that the models return lower distances and higher mass-to-light ratios than expected when we include proper motion stars with high errors or for which there is some blending. We believe this not a fault of our models but is instead due to underestimates or missing systematic uncertainties in the provided errors.

preprint2013arXiv

Supermassive Black Holes and Their Host Galaxies - I. Bulge luminosities from dedicated near-infrared data

In an effort to secure, refine and supplement the relation between central Supermassive Black Hole masses (Mbh), and the bulge luminosities of their host galaxies, (Lbul), we obtained deep, high spatial resolution K-band images of 35 nearby galaxies with securely measured Mbh, using the wide-field WIRCam imager at the Canada-France-Hawaii-Telescope (CFHT). A dedicated data reduction and sky subtraction strategy was adopted to estimate the brightness and structure of the sky, a critical step when tracing the light distribution of extended objects in the near-infrared. From the final image product, bulge and total magnitudes were extracted via two-dimensional profile fitting. As a first order approximation, all galaxies were modeled using a simple Sersic-bulge + exponential-disk decomposition. However, we found that such models did not adequately describe the structure that we observe in a large fraction of our sample galaxies which often include cores, bars, nuclei, inner disks, spiral arms, rings and envelopes. In such cases, we adopted profile modifications and/or more complex models with additional components. The derived bulge magnitudes are very sensitive to the details and number of components used in the models, although total magnitudes remain almost unaffected. Usually, but not always, the luminosities and sizes of the bulges are overestimated when a simple bulge+disk decomposition is adopted in lieu of a more complex model. Furthermore we found that some spheroids are not well fit when the ellipticity of the Sersic model is held fixed. This paper presents the details of the image processing and analysis, while in a companion paper we discuss how model-induced biases and systematics in bulge magnitudes impact the Mbh-Lbul relation.

preprint2013arXiv

Supermassive Black Holes and Their Host Galaxies - II. The correlation with near-infrared luminosity revisited

We present an investigation of the scaling relations between Supermassive Black Hole (SMBH) masses (Mbh), and their host galaxies' K-band bulge (Lbul) and total (Ltot) luminosities. The wide-field WIRCam imager at the Canada-France-Hawaii-Telescope (CFHT) was used to obtain the deepest and highest resolution near infrared images available for a sample of 35 galaxies with securely measured Mbh, selected irrespective of Hubble type. For each galaxy, we derive bulge and total magnitudes using a two-dimensional image decomposition code that allows us to account, if necessary, for large- and small-scale disks, cores, bars, nuclei, rings, envelopes and spiral arms. We find that the present-day Mbh-Lbul and Mbh-Ltot relations have consistent intrinsic scatter, suggesting that Mbh correlates equally well with bulge and total luminosity of the host. Our analysis provides only mild evidence of a decreased scatter if the fit is restricted to elliptical galaxies. The log-slopes of the Mbh-Lbul and Mbh-Ltot relations are 0.75+/-0.10 and 0.92+/-0.14, respectively. However, while the slope of the Mbh-Lbul relation depends on the detail of the image decomposition, the characterization of Mbh-Ltot does not. Given the difficulties and ambiguities of decomposing galaxy images into separate components, our results indicate that Ltot is more suitable as a tracer of SMBH mass than Lbul, and that the Mbh-Ltot relation should be used when studying the co-evolution of SMBHs and galaxies.

preprint2013arXiv

The central mass and mass-to-light profile of the Galactic globular cluster M15

We analyze line-of-sight velocity and proper motion data of stars in the Galactic globular cluster M15 using a new method to fit dynamical models to discrete kinematic data. Our fitting method maximizes the likelihood for individual stars and, as such, does not suffer the same loss of spatial and velocity information incurred when spatially binning data or measuring velocity moments. In this paper, we show that the radial variation in M15 of the mass-to-light ratio is consistent with previous estimates and theoretical predictions, which verifies our method. Our best-fitting axisymmetric Jeans models do include a central dark mass of $\sim2 \pm 1\cdot 10^3M_\odot$, which can be explained by a high concentration of stellar remnants at the cluster center. This paper shows that, from a technical point of view and with current computing power, spatial binning of data is no longer necessary. This not only leads to more accurate fits, but also avoids biased mass estimates due to the loss of resolution. Furthermore, we find that the mass concentration in M15 is significantly higher than previously measured, and is in close agreement with theoretical predictions for core-collapsed globular clusters without a central intermediate-mass black hole.

preprint2013arXiv

The complex nature of the nuclear star cluster in FCC 277

Recent observations have shown that compact nuclear star clusters (NSCs) are present in up to 80% of galaxies. However, detailed studies of their dynamical and chemical properties are confined mainly to spiral galaxy hosts, where they are more easily observed. In this paper we present our study of the NSC in FCC 277, a nucleated elliptical galaxy in the Fornax cluster. We use a combination of adaptive optics assisted near-infrared integral field spectroscopy, Hubble Space Telescope imaging, and literature long slit data. We show that while the NSC does not appear to rotate within our detection limit of ~6 km/s, rotation is detected at larger radii, where the isophotes appear to be disky, suggesting the presence of a nuclear disk. We also observe a distinct central velocity dispersion drop that is indicative of a dynamically cold rotating sub-system. Following the results of orbit-based dynamical modelling, co-rotating as well as counter-rotating stellar orbits are simultaneously needed to reproduce the observed kinematics. We find evidence for varying stellar populations, with the NSC and nuclear disk hosting younger and more metal rich stars than the main body of the galaxy. We argue that gas dissipation and some level of merging have likely played an important role in the formation of the nucleus of this intermediate-mass galaxy. This is in contrast to NSCs in low-mass early- type galaxies, which may have been formed primarily through the infall of star clusters.

preprint2013arXiv

The Gravitational Potential Near the Sun From SEGUE K-dwarf Kinematics

To constrain the Galactic gravitational potential near the Sun ($\sim$1.5 kpc), we derive and model the spatial and velocity distribution for a sample of 9000 K-dwarfs that have spectra from SDSS/SEGUE, which yield radial velocities and abundances ([Fe/H] & [$α$/Fe]). We first derive the spatial density distribution for stars of three abundance-selected sub-populations by accounting for the survey's selection function. The vertical profile of these sub-populations are simple exponentials and their vertical dispersion profile is nearly isothermal. To model these data, we apply the `vertical' Jeans Equation, which relates the observable tracer number density and vertical velocity dispersion to the gravitational potential or vertical force. We explore a number of functional forms for the vertical force law, and fit the dispersion and density profiles of all abundance selected sub-populations simultaneously in the same potential, and explore all parameter co-variances using MCMC. Our fits constrain a disk {\it mass} scale height $\lesssim$ 300 pc and the total surface mass density to be $67 \pm 6 M_{\odot} {\rm pc^{-2}}$ at $|z| = 1.0$ kpc of which the contribution from all stars is $42 \pm 5 M_{\odot} {\rm pc^{-2}}$ (presuming a contribution from cold gas of $13 M_{\odot} {\rm pc^{-2}}$). We find significant constraints on the local dark matter density of $0.0065\pm0.0023 M_{\odot} {\rm pc^{-3}}$ ($0.25\pm0.09 {\rm GeV cm^{-3}} $). Together with recent experiments this firms up the best estimate of $0.0075\pm0.0021 M_{\odot} {\rm pc^{-3}}$ ($0.28\pm0.08 {\rm GeV cm^{-3}} $), consistent with global fits of approximately round dark matter halos to kinematic data in the outskirts of the Galaxy.

preprint2013arXiv

Virgo Cluster and field dwarf ellipticals in 3D - II. Internal dynamics points to tidal harassment?

We present the dynamical analysis of a sample of 12 dwarf elliptical (dE) galaxies for which we have obtained SAURON large-scale two-dimensional spectroscopic data. We construct Jeans axisymmetric models and obtain total dynamical masses enclosed within one effective radius. We use the obtained values to show that the validity of the dynamical scaling relations of massive early-type galaxies can be extended to these low-mass systems, except that dEs seem to contain relatively larger fraction of dark matter in their inner parts. We then demonstrate that dE galaxies have lower angular momenta than the present-day analogues of their presumed late-type progenitors and we show that dE circular velocity curves are steeper than the rotation curves of galaxies with equal and up to an order of magnitude higher luminosity. This requires a transformation mechanism that is not only able to lower the angular momentum but also one that needs to account for increased stellar concentration. Additionally, we match the dark matter fraction of our galaxies to their location in the Virgo Cluster and see that galaxies in the cluster outskirts tend to have a higher dark-to-stellar matter ratio. Transformation due to tidal harassment is able to expain all of the above, unless the dE progenitors were already compact and had lower angular momenta at higher redshifts.

preprint2012arXiv

3D view on Virgo and field dwarf elliptical galaxies: late-type origin and environmental transformations

We show the effects of environmental evolution on Virgo cluster and field dwarf elliptical galaxies (dEs), presenting the first large-scale integral-field spectroscopic data for this galaxy class. The great variety of morphological, kinematic, and stellar population parameters seen in our data supports the claim that dEs are defunct dwarf spiral/irregular galaxies and points to a formation scenario that allows for a stochastic shaping of galaxy properties. We further investigate the properties of our sample by analyzing its kinematic and dynamical properties. We compare the level of rotational support of dEs and giant early-type galaxies and show that the properties of the former largely resemble those of giant fast-rotators. Based on our data, no trend exists between the level of rotational support in dEs and their location in the cluster. However, a tentative trend is seen in dark matter fraction: it increases for larger Virgocentric distances.

preprint2012arXiv

A census of orbital properties of the M31 satellites

We present an analysis of the dynamics of the M31 satellite system. Proper motion data exist for only two of the M31 satellites. We account for this incompleteness in velocity data by a statistical analysis using a combination of the timing argument and phase-space distribution functions. The bulk of the M31 satellites are well fit by these models and we offer a table of orbital properties, including period, eccentricity and semi-major axis. This enables us to search for evidence of group infall based on orbital similarity rather than propinquity on the sky. Our results favour an association between Cass II and NGC 185, as the orbital parameters are in close agreement, but not for NGC 185 and NGC 147, which have often been associated in the past. Other possible satellite groupings include the pair And I and And XVII; the pair And IX and And X; and the triple And V, And XXV and NGC 147. And XXII has been claimed as a satellite of M33; we find that they are not moving independently along the same orbit, but cannot determine whether they are orbiting each other or are unrelated. Two satellites, And XII and And XIV, have high line-of-sight velocities, consistent with very recent infall from the edge of the Local Group. They are not well described by our underlying smooth phase space distribution function, and are reanalysed without priors on their orbital parameters. For And XIV, multiple pericentric passages are possible and improved distance information is needed to draw further conclusions. For And XII, orbits which assume at least one pericentric passage can be ruled out and it must be on its first infall into the M31 system.

preprint2012arXiv

A resonant feature near the Perseus arm revealed by red clump stars

We investigate the extinction together with the radial velocity dispersion and distribution of red clump stars in the anti-center direction using spectra obtained with Hectospec on the MMT. We find that extinction peaks at Galactocentric radii of about 9.5 and 12.5 kpc, right in front of the locations of the Perseus and Outer arms and in line with the relative position of dust and stars in external spiral galaxies. The radial velocity dispersion peaks around 10kpc, which coincides with the location of the Perseus arm, yields an estimated arm-interarm density contrast of 1.3-1.5 and is in agreement with previous studies. Finally, we discover that the radial velocity distribution bifurcates around 10-11 kpc into two peaks at +27 km/s and -4 km/s. This seems to be naturally explained by the presence of the outer Lindblad resonance of the Galactic bar, but further observations will be needed to understand if the corotation resonance of the spirals arms also plays a role.

preprint2012arXiv

An Over-Massive Black Hole in the Compact Lenticular Galaxy NGC1277

All massive galaxies likely have supermassive black holes at their centers, and the masses of the black holes are known to correlate with properties of the host galaxy bulge component. Several explanations have been proposed for the existence of these locally-established empirical relationships; they include the non-causal, statistical process of galaxy-galaxy merging, direct feedback between the black hole and its host galaxy, or galaxy-galaxy merging and the subsequent violent relaxation and dissipation. The empirical scaling relations are thus important for distinguishing between various theoretical models of galaxy evolution, and they further form the basis for all black hole mass measurements at large distances. In particular, observations have shown that the mass of the black hole is typically 0.1% of the stellar bulge mass of the galaxy. The small galaxy NGC4486B currently has the largest published fraction of its mass in a black hole at 11%. Here we report observations of the stellar kinematics of NGC 1277, which is a compact, disky galaxy with a mass of 1.2 x 10^11 Msun. From the data, we determine that the mass of the central black hole is 1.7 x 10^10 Msun, or 59% its bulge mass. Five other compact galaxies have properties similar to NGC 1277 and therefore may also contain over-sized black holes. It is not yet known if these galaxies represent a tail of a distribution, or if disk-dominated galaxies fail to follow the normal black hole mass scaling relations.

preprint2012arXiv

Chemo-orbital evidence from SDSS/SEGUE G-type dwarf stars for a mixed origin of the Milky Way's thick disk

We combine the estimated metallicities [Fe/H], abundances [α/Fe], positions and motions of a sample of 27,500 local (7<R/kpc<9, 0.5<|z|/kpc<2.5) SDSS/SEGUE G-type dwarf stars to investigate the chemo-orbital properties of the Milky Way's disk around the Sun. When we derive the orbital properties reflecting angular momentum, circularity, and thickness as function of [α/Fe] vs. [Fe/H], we find that there is a smooth variation with [α/Fe], a proxy for age. At the same time, the orbital properties of the old stars with [α/Fe]$\gtrsim$0.25 do show a transition with [Fe/H]: below [Fe/H]$\simeq$-0.6 the orbital angular momentum decreases, and the orbits become significantly non-circular and thicker. Radial migration of stars into the Solar neighborhood would naturally result in a smooth variation in the orbital properties, but the latter old metal-poor stars form a clear challenge, in particular because a basic feature of radial migration is that stars remain on near-circular orbits. When we next select stars on near-circular orbits, we indeed find besides the α-young 'thin-disk' stars a significant contribution to the α-old 'thick-disk' metal-rich stars. However, the remaining α-old 'thick-disk' stars on eccentric orbits, including nearly all old metal-poor stars, are difficult to explain with radial migration alone, but might have formed through early-on gas-rich mergers. We thus find chemo-orbital evidence that the thicker component of the Milky Way disk is not distinct from the thin component as expected from smooth internal evolution through radial migration, except for the old metal-poor stars with different orbital properties which could be part of a distinct thick-disk component formed through an external mechanism.

preprint2012arXiv

The SAURON project-XXI. The spatially-resolved UV-line strength relations of early-type galaxies

The unexpected rising flux of early-type galaxies at decreasing ultraviolet (UV) wavelengths is a long-standing mystery. One important observational constraint is the correlation between UV-optical colours and Mg2 line strengths found by Burstein et al. (1988). The simplest interpretation of this phenomenon is that the UV strength is related to the Mg line strength. Under this assumption, we expect galaxies with larger Mg gradients to have larger UV colour gradients. By combining UV imaging from GALEX, optical imaging from MDM and SAURON integral-field spectroscopy, we investigate the spatially-resolved relationships between UV colours and stellar population properties of 34 early-type galaxies from the SAURON survey sample. We find that galaxies with old stellar populations show tight correlations between the FUV colours (FUV-V and FUV-NUV) and the Mgb index, Hβ index and metallicity [Z/H]. The equivalent correlations for the Fe5015 index, α-enhancement [α/Fe] and age are present but weaker. We have also derived logarithmic internal radial colour, measured line strength and derived stellar population gradients for each galaxy and again found a strong dependence of the FUV-V and FUV-NUV colour gradients on both the Mg b line strength and the metallicity gradients for galaxies with old stellar populations. In particular, global gradients of Mg b and [Z/H] with respect to the UV colour across galaxies are consistent with their local gradients within galaxies, suggesting that the global correlations also hold locally. From a simple model based on multi-band colour fits of UV upturn and UV-weak galaxies, we have identified a plausible range of parameters that reproduces the observed radial colour profiles. In these models, the centers of elliptical galaxies, where the UV flux is strong, are enhanced in metals by roughly 60% compared to UV-weak regions.

preprint2012arXiv

Virgo Cluster and field dwarf ellipticals in 3D: I. On the variety of stellar kinematic and line-strength properties

We present the first large-scale stellar kinematic and line-strength maps for dwarf elliptical galaxies (9 in the Virgo Cluster and 3 in the field environment) obtained with the SAURON integral-field unit. No two galaxies in our sample are alike: we see that the level of rotation is not tied to flattening (we have, e.g. round rotators and flattened nonrotators); we observe kinematic twists in 1 Virgo and 1 field object; we discover large-scale kinematically-decoupled components in 2 field galaxies; we see varying gradients in line-strength maps, from nearly flat to strongly peaked in the center. The great variety of morphological, kinematic, and stellar population parameters seen in our data points to a formation scenario in which properties are shaped stochastically. A combined effect of ram-pressure stripping and galaxy harassment is the most probable explanation. We show the need for a comprehensive analysis of kinematic, dynamical, and stellar population properties which will enable us to place dwarf ellipticals and processes that govern their evolution in the wider context of galaxy formation.

preprint2011arXiv

Physical Condition of Molecular Gas at the Centre of NGC 1097

We have used the Xco conversion factor, Local Thermal Equilibrium and Large Velocity Gradient approximation to parametrize the cold and warm phase of the interstellar medium from five different low transitions of CO in the central 21"(kpc) of NGC 1097. We have applied a one-component model and derived a typical kinetic temperature of ~33K, and its molecular Hydrogen density and a CO column density. A two-component model results in 85% cold-to-total gas fraction with a 90K warm counterpart. Furthermore, we "resolve" the spatially unresolved single dish observations by selecting velocity channels that in an interferometric velocity map correspond to specific regions. We have selected five such regions and found that the physical properties in these regions are comparable to those derived from the full line profile. This implies that the central kpc of NGC 1097 is rather homogeneous in nature, and, although the regions are not uniquely located within the ring, the star formation along the ring is homogeneously distributed (in agreement with recent Herschel observations). We have further revised the mass inflow rate onto the Supermassive Black Hole in this prototype LINER/Sy1 galaxy and found that, accounting for the total interstellar medium and applying a careful contribution of the disc thickness and corresponding stability criterion, increases the previous estimations by a factor 10. Finally we have calculated the Xco for the centre of NGC 1097 using an independent estimation of the surface density to the CO emission, and obtained Xco=(2.8+-0.5)x10^20 cm^-2(K km s^-1)^-1 at radius 10.5" and Xco=(5.0+-0.5)x10^20 cm^-2(K km s^-1)^-1 at 7.5". With the approach and analysis described, we have demonstrated that important physical properties can be derived to a resolution beyond the single dish resolution element, however, caution is necessary for interpreting the results.(Abriged)

preprint2011arXiv

Resolving the dynamical mass of a z~1.3 QSO host galaxy using SINFONI and Laser Guide Star assisted Adaptive Optics

Recent studies of the tight scaling relations between the masses of supermassive black holes and their host galaxies have suggested that in the past black holes constituted a larger fraction of their host galaxies' mass. However, these arguments are limited by selection effects and difficulties in determining robust host galaxy masses at high redshifts. Here we report the first results of a new, complementary diagnostic route: we directly determine a dynamical host galaxy mass for the z=1.3 luminous quasar J090543.56+043347.3 through high-spatial-resolution (0.47", 4kpc FWHM) observations of the host galaxy gas kinematics over 30x40 kpc using ESO/VLT/SINFONI with LGS/AO. Combining our result of M_dyn = 2.05+1.68_0.74 x 10^11 M_sun (within a radius 5.25 +- 1.05 kpc) with M_BH,MgII = 9.02 \pm 1.43 x 10^8 M_sun, M_BH,Halpha = 2.83 +1.93-1.13 x 10^8 M_sun, we find that the ratio of black hole mass to host galaxy dynamical mass for J090543.56+043347.3 matches the present-day relation for M_BH vs. M_Bulge,Dyn, well within the IR scatter, deviating at most a factor of two from the mean. J090543.56+043347.3 displays clear signs of an ongoing tidal interaction and of spatially extended star formation at a rate of 50-100 M_sun/yr, above the cosmic average for a galaxy of this mass and redshift. We argue that its subsequent evolution may move J090543.56+043347.3 even closer to the z=0 relation for M_BH vs. M_Bulge,Dyn. Our results support the picture where any substantive evolution in these relations must occur prior to z~1.3. Having demonstrated the power of this modelling approach we are currently analyzing similar data on seven further objects to better constrain such evolution.

preprint2011arXiv

The ATLAS3D project - III. A census of the stellar angular momentum within the effective radius of early-type galaxies: unveiling the distribution of Fast and Slow Rotators

We provide a census of the apparent stellar angular momentum within 1 Re of a volume-limited sample of 260 early-type galaxies (ETGs) in the nearby Universe, using integral-field spectroscopy obtained in the course of the ATLAS3D project. We exploit the LambdaR parameter to characterise the existence of two families of ETGs: Slow Rotators which exhibit complex stellar velocity fields and often include stellar kinematically Distinct Cores (KDCs), and Fast Rotators which have regular velocity fields. Our complete sample of 260 ETGs leads to a new criterion to disentangle Fast and Slow Rotators which now includes a dependency on the apparent ellipticity (Epsilon). It separates the two classes significantly better than the previous prescription, and than a criterion based on V/Sigma: Slow Rotators and Fast Rotators have LambdaR lower and larger than kFSxSQRT(Epsilon), respectively, where kFS=0.31 for measurements made within 1 Re. We show that the vast majority of early-type galaxies are Fast Rotators: these have regular stellar rotation, with aligned photometric and kinematic axes (Paper II, Krajnovic et al. 2011}, include discs and often bars and represent 86% (224/260) of all early-type galaxies in the volume-limited ATLAS3D sample. Fast Rotators span the full range of apparent ellipticities from 0 to 0.85, and we suggest that they cover intrinsic ellipticities from about 0.35 to 0.85, the most flattened having morphologies consistent with spiral galaxies. Only a small fraction of ETGs are Slow Rotators representing 14% (36/260) of the ATLAS3D sample of ETGs. Of all Slow Rotators, 11% (4/36) exhibit two counter-rotating stellar disc-like components and are rather low mass objects (Mdyn<10^10.5 M_Sun). All other Slow Rotators (32/36) appear relatively round on the sky (Epsilon_e<0.4), tend to be massive (Mdyn>10^10.5 M_Sun), and often (17/32) exhibit KDCs.

preprint2011arXiv

The central dark matter distribution of NGC 2976

We study the mass distribution in the late-type dwarf galaxy NGC 2976 through stellar kinematics obtained with the VIRUS-P integral-field spectrograph and anisotropic Jeans models as a test of cosmological simulations and baryonic processes that putatively alter small-scale structure. Previous measurements of the H-alpha emission-line kinematics have determined that the dark matter halo of NGC 2976 is most consistent with a cored density profile. We find that the stellar kinematics are best fit with a cuspy halo. Cored dark matter halo fits are only consistent with the stellar kinematics if the stellar mass-to-light ratio is significantly larger than that derived from stellar population synthesis, while the best-fitting cuspy model has no such conflict. The inferred mass distribution from a harmonic decomposition of the gaseous kinematics is inconsistent with that of the stellar kinematics. This difference is likely due to the gas disk not meeting the assumptions that underlie the analysis such as no pressure support, a constant kinematic axis, and planar orbits. By relaxing some of these assumptions, in particular the form of the kinematic axis with radius, the gas-derived solution can be made consistent with the stellar kinematic models. A strong kinematic twist in the gas of NGC 2976's center suggests caution, and we advance the mass model based on the stellar kinematics as more reliable. The analysis of this first galaxy shows promising evidence that dark matter halos in late-type dwarfs may in fact be more consistent with cuspy dark matter distributions than earlier work has claimed.

preprint2011arXiv

The SAURON project - XVIII. The integrated UV-linestrength relations of early-type galaxies

Using far (FUV) and near (NUV) ultraviolet photometry from guest investigator programmes on the Galaxy Evolution Explorer (GALEX) satellite, optical photometry from the MDM Observatory and optical integral-field spectroscopy from SAURON, we explore the UV-linestrength relations of the 48 nearby early-type galaxies in the SAURON sample. Identical apertures are used for all quantities, avoiding aperture mismatch. We show that galaxies with purely old stellar populations show well-defined correlations of the integrated FUV-V and FUV-NUV colours with the integrated Mgb and Hbeta absorption linestrength indices, strongest for FUV-NUV. Correlations with the NUV-V colour, Fe5015 index and stellar velocity dispersion are much weaker. These correlations put stringent constraints on the origin of the UV-upturn phenomenon in early-type galaxies, and highlight its dependence on age and metallicity. In particular, despite recent debate, we recover the negative correlation between FUV-V colour and Mg linestrength originally publicised by Burstein et al. (1988), which we refer to as the "Burstein relation", suggesting a positive dependence of the UV-upturn on metallicity. We argue that the scatter in the correlations is real, and present mild evidence that a strong UV excess is preferentially present in slow-rotating galaxies. We also demonstrate that most outliers in the correlations are galaxies with current or recent star formation, some at very low levels. We believe that this sensitivity to weak star formation, afforded by the deep and varied data available for the SAURON sample, explains why our results are occasionally at odds with other recent but shallower surveys. This is supported by the analysis of a large, carefully-crafted sample of more distant early-type galaxies from the Sloan Digital Sky Survey (SDSS), more easily comparable with current and future large surveys.

preprint2011arXiv

The SAURON Project - XX. The Spitzer [3.6] - [4.5] colour in early-type galaxies: colours, colour gradients and inverted scaling relations

We investigate the [3.6] - [4.5] Spitzer-IRAC colour behaviour of the early-type galaxies of the SAURON survey, a representative sample of 48 nearby ellipticals and lenticulars. We investigate how this colour, which is unaffected by dust extinction, can be used to constrain the stellar populations in these galaxies. We find a tight relation between the [3.6]-[4.5] colour and effective velocity dispersion, a good mass-indicator in early-type galaxies. Contrary to other colours in the optical and near-infrared, we find that the colours become bluer for larger galaxies. The relations are tighter when using the colour inside r_e, rather than the much smaller r_e/8 aperture, due to the presence of young populations in the central regions. We also obtain strong correlations between the [3.6]-[4.5] colour and 3 strong absorption lines (H beta, Mg b and Fe 5015). Comparing our data with the models of Marigo et al., which show that more metal rich galaxies are bluer, we can explain our results in a way consistent with results from the optical, by stating that larger galaxies are more metal rich. The blueing is caused by a strong CO absorption band, whose line strength increases strongly with decreasing temperature and which covers a considerable fraction of the 4.5 micron filter. In galaxies that contain a compact radio source, the [3.6]-[4.5] colour is generally slightly redder than in the other galaxies, indicating small amounts of either hot dust, non-thermal emission, or young stars near the center. We find that the large majority of the galaxies show redder colours with increasing radius. abbreviated...

preprint2010arXiv

Kinematic analysis of nuclear spirals: feeding the black hole in NGC1097

We present a harmonic expansion of the observed line-of-sight velocity field as a method to recover and investigate spiral structures in the nuclear regions of galaxies. We apply it to the emission-line velocity field within the circumnuclear starforming ring of NGC1097, obtained with the GMOS-IFU spectrograph. The radial variation of the third harmonic terms are well described by a logarithmic spiral, from which we interpret that the gravitational potential is weakly perturbed by a two-arm spiral density wave with inferred pitch angle of of 52+/-4 degrees. This interpretation predicts a two-arm spiral distortion in the surface brightness, as hinted by the dust structures in central images of NGC1097, and predicts a combined one-arm and three-arm spiral structure in the velocity field, as revealed in the non-circular motions of the ionised gas within the circumnuclear region of this galaxy. Next, we use a simple spiral perturbation model to constrain the fraction of the measured non-circular motions that is due to radial inflow. We combine the resulting inflow velocity with the gas density in the spiral arms, inferred from emission line ratios, to estimate the mass inflow rate as a function of radius, which reaches about 0.011 Msun/yr at a distance of 70 pc from the center. This value corresponds to a fraction of about 4.2 x 10^{-3} of the Eddington mass accretion rate onto the central black hole in this LINER/Seyfert1 galaxy. We conclude that the line-of-sight velocity not only can provide a cleaner view of nuclear spirals than the associated dust, but that the presented method also allows the quantitative study of these possibly important links in fueling the centers of galaxies, including providing a handle on the mass inflow rate as a function of radius.

preprint2010arXiv

Short gamma-ray bursts from dynamically-assembled compact binaries in globular clusters: pathways, rates, hydrodynamics and cosmological setting

We present a detailed assessment of the dynamical pathways leading to the coalescence of compact objects in Globular Clusters (GCs) and Short Gamma-Ray Burst (SGRB) production. We consider primordial binaries, dynamically formed binaries (through tidal two-body and three-body exchange interactions) and direct impacts of compact objects (WD/NS/BH). We show that if the primordial binary fraction is small, close encounters dominate the production rate of coalescing compact systems. We find that the two dominant channels are the interaction of field NSs with dynamically formed binaries, and two-body encounters. We then estimate the redshift distribution and host galaxy demographics of SGRB progenitors, and find that GCs can provide a significant contribution to the overall observed rate. We have carried out hydrodynamical modeling of evolution of close stellar encounters with WD/NS/BH, and show that there is no problem in accounting for the energy budget of a typical SGRB. The particulars of each encounter are variable and lead to interesting diversity: the encounter characteristics are dependent on the impact parameter, in contrast to the merger scenario; the nature of the compact star itself can produce very different outcomes; the presence of tidal tails in which material falls back onto the central object at later times is a robust feature of these calculations, with the mass involved being larger than for binary mergers. It is thus possible to account generically in this scenario for a prompt episode of energy release, as well as for activity many dynamical time scales later (abridged).

preprint2010arXiv

The Einstein Cross: constraint on dark matter from stellar dynamics and gravitational lensing

[Abridged] We present two-dimensional line-of-sight stellar kinematics of the lens galaxy in the Einstein Cross, obtained with the GEMINI 8m telescope, using the GMOS integral-field spectrograph. The velocity map shows regular rotation up to ~100 km/s around the minor axis of the bulge, consistent with axisymmetry. The velocity dispersion map shows a weak gradient increasing towards a central (R<1") value of sigma_0=170+/-9 km/s. We deproject the observed surface brightness from HST imaging to obtain a realistic luminosity density of the lens galaxy, which in turn is used to build axisymmetric dynamical models that fit the observed kinematic maps. We also construct a gravitational lens model that accurately fits the positions and relative fluxes of the four quasar images. We find that the resulting luminous and total mass distribution are nearly identical around the Einstein radius R_E = 0.89", with a slope that is close to isothermal, but which becomes shallower towards the center if indeed mass follows light. The dynamical model fits to the observed kinematic maps result in a total mass-to-light ratio (M/L)_dyn=3.7+/-0.5 M_sun/L_sun,I (in the I-band). This is consistent with the Einstein mass M_E = 1.54 x 10^10 M_sun divided by the (projected) luminosity within R_E, which yields a total mass-to-light ratio of (M/L)_E=3.4 M_sun/L_sun,I, with an error of at most a few per cent. We estimate from stellar populations model fits to colors of the lens galaxy a stellar mass-to-light ratio (M/L)_* from 2.8 to 4.1 M_sun/L_sun,I. Although a constant dark matter fraction of 20 per cent is not excluded, dark matter may play no significant role in the bulge of this ~L* early-type spiral galaxy.

preprint2010arXiv

The SAURON project - XVII. Stellar population analysis of the absorption line strength maps of 48 early-type galaxies

(Abridged) We present a stellar population analysis of the absorption line strength maps for 48 early-type galaxies from the SAURON sample. Using the line strength index maps of Hbeta, Fe5015, and Mgb, measured in the Lick/IDS system and spatially binned to a constant signal-to-noise, together with predictions from up-to-date stellar population models, we estimate the simple stellar population-equivalent (SSP-equivalent) age, metallicity and abundance ratio [alpha/Fe] over a two-dimensional field extending up to approximately one effective radius. We find a large range of SSP-equivalent ages in our sample, of which ~40% of the galaxies show signs of a contribution from a young stellar population. The most extreme cases of post-starburst galaxies, with SSP-equivalent ages of <=3 Gyr observed over the full field-of-view, and sometimes even showing signs of residual star-formation, are restricted to low mass systems(sigma_e <= 100 k/ms or ~2x10^10 M_sol). Spatially restricted cases of young stellar populations in circumnuclear regions can almost exclusively be linked to the presence of star-formation in a thin, dusty disk/ring, also seen in the near-UV or mid-IR. The flattened components with disk-like kinematics previously identified in all fast rotators (Krajnović et al.) are shown to be connected to regions of distinct stellar populations. These range from the young, still star-forming circumnuclear disks and rings with increased metallicity preferentially found in intermediate-mass fast rotators, to apparently old structures with extended disk-like kinematics, which are observed to have an increased metallicity and mildly depressed [alpha/Fe] ratio compared to the main body of the galaxy. The slow rotators generally show no stellar population signatures over and above the well known metallicity gradients and are largely consistent with old (>=10 Gyr) stellar populations.

preprint2009arXiv

Galaxy density profiles and shapes -- I. simulation pipeline for lensing by realistic galaxy models

Studies of strong gravitational lensing in current and upcoming wide and deep photometric surveys, and of stellar kinematics from (integral-field) spectroscopy at increasing redshifts, promise to provide valuable constraints on galaxy density profiles and shapes. However, both methods are affected by various selection and modelling biases, whch we aim to investigate in a consistent way. In this first paper in a series we develop a flexible but efficient pipeline to simulate lensing by realistic galaxy models. These galaxy models have separate stellar and dark matter components, each with a range of density profiles and shapes representative of early-type, central galaxies without significant contributions from other nearby galaxies. We use Fourier methods to calculate the lensing properties of galaxies with arbitrary surface density distributions, and Monte Carlo methods to compute lensing statistics such as point-source lensing cross-sections. Incorporating a variety of magnification bias modes lets us examine different survey limitations in image resolution and flux. We rigorously test the numerical methods for systematic errors and sensitivity to basic assumptions. We also determine the minimum number of viewing angles that must be sampled in order to recover accurate orientation-averaged lensing quantities. We find that for a range of non-isothermal stellar and dark matter density profiles typical of elliptical galaxies, the combined density profile and corresponding lensing properties are surprisingly close to isothermal around the Einstein radius. The converse implication is that constraints from strong lensing and/or stellar kinematics, which are indeed consistent with isothermal models near the Einstein radius, cannot trivially be extrapolated to smaller and larger radii.

preprint2009arXiv

Galaxy density profiles and shapes -- II. selection biases in strong lensing surveys

[Abridged] Many current and future astronomical surveys will rely on samples of strong gravitational lens systems to draw conclusions about galaxy mass distributions. We use a new strong lensing pipeline (presented in Paper I of this series) to explore selection biases that may cause the population of strong lensing systems to differ from the general galaxy population. Our focus is on point-source lensing by early-type galaxies with two mass components (stellar and dark matter) that have a variety of density profiles and shapes motivated by observational and theoretical studies of galaxy properties. We seek not only to quantify but also to understand the physics behind selection biases related to: galaxy mass, orientation and shape; dark matter profile parameters such as inner slope and concentration; and adiabatic contraction. We study how all of these properties affect the lensing Einstein radius, total cross-section, quad/double ratio, and image separation distribution. We find significant (factors of several) selection biases with mass; orientation, for a given galaxy shape at fixed mass; cusped dark matter profile inner slope and concentration; concentration of the stellar and dark matter deprojected Sersic models. Interestingly, the intrinsic shape of a galaxy does not strongly influence its lensing cross-section when we average over viewing angles. Our results are an important first step towards understanding how strong lens systems relate to the general galaxy population.

preprint2009arXiv

Migration of Star Clusters and Nuclear Rings

Star clusters that form in nuclear rings appear to be at slightly larger radii than the gas. We argue that the star clusters move out from the gas in which they are formed because of satellite-disk tidal interactions. In calculating the dynamics of this star cluster and gas ring system, we include the effects of dynamical friction of the background stars in the host galaxy on the star cluster, and inflowing gas along the bar onto the nuclear ring at the two contact points. We show that the final separation is of the order of the Hill radius of the nuclear ring, which is typically 20-30% of its radius. Massive star clusters can reach half of this separation very quickly and produce a factor of a few enhancement in the gas surface density. If this leads to star formation in addition to the (ongoing) formation of star clusters near the contact points, a possible (initial) azimuthal age gradient may become diluted or even disappear. Finally, if the star cluster are massive and/or numerous enough, we expect the nuclear ring to migrate inward, away from the (possibly) associated (inner) Lindblad resonance. We discuss how these predictions may be tested observationally.

preprint2009arXiv

The SAURON project - XIII. SAURON-GALEX study of early-type galaxies: the ultraviolet colour-magnitude relations and Fundamental Planes

We present GALEX far (FUV) and near (NUV) ultraviolet imaging of 34 nearby early-type galaxies from the SAURON representative sample of 48 E/S0 galaxies, all of which have ground-based optical imaging from the MDM Observatory. The surface brightness profiles of nine galaxies show regions with blue UV-optical colours suggesting recent star formation. Five of these show blue integrated UV-optical colours that set them aside in the NUV integrated colour-magnitude relation. They also have other properties confirming they have had recent star formation, in particular H_beta absorption higher than expected for a quiescent population and a higher CO detection rate. NUV-blue galaxies are generally drawn from the lower stellar velocity dispersion and thus lower dynamical mass part of the sample. We have also constructed the first UV Fundamental Planes and show that NUV blue galaxies bias the slopes and increase the scatters. If they are eliminated the fits get closer to expectations from the virial theorem. Although our analysis is based on a limited sample, it seems that a dominant fraction of the tilt and scatter of the UV Fundamental Planes is due to the presence of young stars in preferentially low-mass early-type galaxies. Interestingly, the UV-optical radial colour profiles reveal a variety of behaviours, with many galaxies showing signs of recent star formation, a central UV-upturn phenomenon, smooth but large-scale age and metallicity gradients, and in many cases a combination of these. In addition, FUV-NUV and FUV-V colours even bluer than those normally associated with UV-upturn galaxies are observed at the centre of some quiescent galaxies. Four out of the five UV-upturn galaxies are slow rotators. These objects should thus pose interesting challenges to stellar evolutionary models of the UV-upturn.

preprint2005arXiv

A Bar Signature and Central Disk in the Gaseous and Stellar Velocity Fields of NGC5448

We analyse SAURON kinematic maps of the inner kpc of the early-type (Sa) barred spiral galaxy NGC5448. The observed morphology and kinematics of the emission-line gas is patchy and perturbed, indicating clear departures from circular motion. The kinematics of the stars is more regular, and display a small inner disk-like system embedded in a large-scale rotating structure. We focus on the [OIII] gas, and use a harmonic decomposition formalism to analyse the gas velocity field. The higher-order harmonic terms and the main kinematic features of the observed data are consistent with an analytically constructed simple bar model. The bar model is derived using linear theory, considering an m=2 perturbation mode, and with bar parameters which are consistent with the large-scale bar detected via imaging. We also study optical and near infra-red images to reveal the asymmetric extinction in NGC5448, and we recognise that some of the deviations between the data and the analytical bar model may be due to these complex dust features. Our study illustrates how the harmonic decomposition formalism can be used as a powerful tool to quantify non-circular motions in observed gas velocity fields.

preprint2005arXiv

The dynamical M/L-profile and distance of the globular cluster M15

We construct orbit-based axisymmetric dynamical models for the globular cluster M15 which fit groundbased line-of-sight velocities and Hubble Space Telescope line-of-sight velocities and proper motions. This allows us to constrain the variation of the mass-to-light ratio M/L as a function of radius in the cluster, and to measure the distance and inclination of the cluster. We obtain a best-fitting inclination of 60+/-15 degrees, a dynamical distance of 10.3+/-0.4 kpc and an M/L profile with a central peak. The inferred mass in the central 0.05 parsec is 3400 Msun, implying a central density of at least 7.4x10^6 Msun pc^-3. We cannot distinguish the nature of the central mass concentration. It could be an IMBH or it could be large number of compact objects, or it could be a combination. The central 4 arcsec of M15 appears to contain a rapidly spinning core, and we speculate on its origin.

preprint2003arXiv

Jeans solutions for triaxial galaxies

The Jeans equations relate the second-order velocity moments to the density and potential of a stellar system. For general three-dimensional stellar systems, there are three equations, but these are not very helpful, as they contain six independent moments. By assuming that the potential is triaxial and of separable Staeckel form, the mixed moments vanish in confocal ellipsoidal coordinates. The three Jeans equations and three remaining non-vanishing moments form a closed system of three highly-symmetric coupled first-order partial differential equations in three variables. They were first derived by Lynden-Bell in 1960, but have resisted solution by standard methods. Here we present the general solution by superposition of singular solutions.